A genetic analysis of cortical thickness in 372 twins

Imaging genetics is a new field of neuroscience that blends methods from computational anatomy and quantitative genetics to identify genetic influences on brain structure and function. Here we analyzed brain MRI data from 372 young adult twins to identify cortical regions in which gray matter volume is influenced by genetic differences across subjects. Thickness maps, reconstructed from surface models of the cortical gray/white and gray/CSF interfaces, were smoothed with a 25 mm FWHM kernel and automatically parcellated into 34 regions of interest per hemisphere. In structural equation models fitted to volume values at each surface vertex, we computed components of variance due to additive genetic (A), shared (C) and unique (E) environmental factors, and tested their significance. Cortical regions in the vicinity of the perisylvian language cortex, and at the frontal and temporal poles, showed significant additive genetic variance, suggesting that volume measures from these regions may provide quantitative phenotypes to narrow the search for quantitative trait loci that influence brain structure.

The contribution of genes to cortical thickness and volume

We analyzed brain MRI data from 372 young adult twins toidentify cortical regions in which gray matter thickness and volume are influenced by genetics. This was achieved using an A/C/E structural equation model that divides the variance of these traits, at each point on the cortex, into additive genetic (A), shared (C), and unique environmental (E) components. A strong genetic influencewas found in frontal and parietal regions. Inaddition, we correlated cortical thickness with full-scale intelligence quotient for comparison with the A/C/E maps, and several regions where cortical structure was correlated with intelligence quotient are under genetic control. These cortical measures may be useful phenotypes to narrow the searchfor quantitative trait lociinfluencing brain structure.

Heritability analysis of surface-based cortical thickness estimation on a large twin cohort

The aim of this paper is to assess the heritability of cerebral cortex, based on measurements of grey matter (GM) thickness derived from structural MR images (sMRI). With data acquired from a large twin cohort (328 subjects), an automated method was used to estimate the cortical thickness, and EM-ICP surface registration algorithm was used to establish the correspondence of cortex across the population. An ACE model was then employed to compute the heritability of cortical thickness. Heritable cortical thickness measures various cortical regions, especially in frontal and parietal lobes, such as bilateral postcentral gyri, superior occipital gyri, superior parietal gyri, precuneus, the orbital part of the right frontal gyrus, right medial superior frontal gyrus, right middle occipital gyrus, right paracentral lobule, left precentral gyrus, and left dorsolateral superior frontal gyrus.

Mapping cortical change in Alzheimer's disease, brain development, and schizophrenia

This paper describes algorithms that can identify patterns of brain structure and function associated with Alzheimer's disease, schizophrenia, normal aging, and abnormal brain development based on imaging data collected in large human populations. Extraordinary information can be discovered with these techniques: dynamic brain maps reveal how the brain grows in childhood, how it changes in disease, and how it responds to medication. Genetic brain maps can reveal genetic influences on brain structure, shedding light on the nature-nurture debate, and the mechanisms underlying inherited neurobehavioral disorders. Recently, we created time-lapse movies of brain structure for a variety of diseases. These identify complex, shifting patterns of brain structural deficits, revealing where, and at what rate, the path of brain deterioration in illness deviates from normal. Statistical criteria can then identify situations in which these changes are abnormally accelerated, or when medication or other interventions slow them. In this paper, we focus on describing our approaches to map structural changes in the cortex. These methods have already been used to reveal the profile of brain anomalies in studies of dementia, epilepsy, depression, childhood- and adult-onset schizophrenia, bipolar disorder, attention-deficit/hyperactivity disorder, fetal alcohol syndrome, Tourette syndrome, Williams syndrome, and in methamphetamine abusers. Specifically, we describe an image analysis pipeline known as cortical pattern matching that helps compare and pool cortical data over time and across subjects. Statistics are then defined to identify brain structural differences between groups, including localized alterations in cortical thickness, gray matter density (GMD), and asymmetries in cortical organization. Subtle features, not seen in individual brain scans, often emerge when population-based brain data are averaged in this way. Illustrative examples are presented to show the profound effects of development and various diseases on the human cortex. Dynamically spreading waves of gray matter loss are tracked in dementia and schizophrenia, and these sequences are related to normally occurring changes in healthy subjects of various ages.

The relationship between stress, HPA axis functioning and brain structure in first episode psychosis over the first 12 weeks of treatment

Stress and abnormal hypothalamic-pituitary-adrenal axis functioning have been implicated in the early phase of psychosis and may partly explain reported changes in brain structure. This study used magnetic resonance imaging to investigate whether biological measures of stress were related to brain structure at baseline and to structural changes over the first 12 weeks of treatment in first episode patients (n=22) compared with matched healthy controls (n=22). At baseline, no significant group differences in biological measures of stress, cortical thickness or hippocampal volume were observed, but a significantly stronger relationship between baseline levels of cortisol and smaller white matter volumes of the cuneus and anterior cingulate was found in patients compared with controls. Over the first 12 weeks of treatment, patients showed a significant reduction in thickness of the posterior cingulate compared with controls. Patients also showed a significant positive relationship between baseline cortisol and increases in hippocampal volume over time, suggestive of brain swelling in association with psychotic exacerbation, while no such relationship was observed in controls. The current findings provide some support for the involvement of stress mechanisms in the pathophysiology of early psychosis, but the changes are subtle and warrant further investigation.

Genetic and environmental influences on neuroimaging phenotypes: A meta-analytical perspective on twin imaging studies

Because brain structure and function are affected in neurological and psychiatric disorders, it is important to disentangle the sources of variation in these phenotypes. Over the past 15 years, twin studies have found evidence for both genetic and environmental influences on neuroimaging phenotypes, but considerable variation across studies makes it difficult to draw clear conclusions about the relative magnitude of these influences. Here we performed the first meta-analysis of structural MRI data from 48 studies on >1,250 twin pairs, and diffusion tensor imaging data from 10 studies on 444 twin pairs. The proportion of total variance accounted for by genes (A), shared environment (C), and unshared environment (E), was calculated by averaging A, C, and E estimates across studies from independent twin cohorts and weighting by sample size. The results indicated that additive genetic estimates were significantly different from zero for all metaanalyzed phenotypes, with the exception of fractional anisotropy (FA) of the callosal splenium, and cortical thickness (CT) of the uncus, left parahippocampal gyrus, and insula. For many phenotypes there was also a significant influence of C. We now have good estimates of heritability for many regional and lobar CT measures, in addition to the global volumes. Confidence intervals are wide and number of individuals small for many of the other phenotypes. In conclusion, while our meta-analysis shows that imaging measures are strongly influenced by genes, and that novel phenotypes such as CT measures, FA measures, and brain activation measures look especially promising, replication across independent samples and demographic groups is necessary.

Tracking Alzheimer's disease

Population-based brain mapping provides great insight into the trajectory of aging and dementia, as well as brain changes that normally occur over the human life span.We describe three novel brain mapping techniques, cortical thickness mapping, tensor-based morphometry (TBM), and hippocampal surface modeling, which offer enormous power for measuring disease progression in drug trials, and shed light on the neuroscience of brain degeneration in Alzheimer's disease (AD) and mild cognitive impairment (MCI).We report the first time-lapse maps of cortical atrophy spreading dynamically in the living brain, based on averaging data from populations of subjects with Alzheimer's disease and normal subjects imaged longitudinally with MRI. These dynamic sequences show a rapidly advancing wave of cortical atrophy sweeping from limbic and temporal cortices into higher-order association and ultimately primary sensorimotor areas, in a pattern that correlates with cognitive decline. A complementary technique, TBM, reveals the 3D profile of atrophic rates, at each point in the brain. A third technique, hippocampal surface modeling, plots the profile of shape alterations across the hippocampal surface. The three techniques provide moderate to highly automated analyses of images, have been validated on hundreds of scans, and are sensitive to clinically relevant changes in individual patients and groups undergoing different drug treatments. We compare time-lapse maps of AD, MCI, and other dementias, correlate these changes with cognition, and relate them to similar time-lapse maps of childhood development, schizophrenia, and HIV-associated brain degeneration. Strengths and weaknesses of these different imaging measures for basic neuroscience and drug trials are discussed.

Sensitivity of proximal femoral stiffness and areal bone mineral density to changes in bone geometry and density

Areal bone mineral density (aBMD) is the most common surrogate measurement for assessing the bone strength of the proximal femur associated with osteoporosis. Additional factors, however, contribute to the overall strength of the proximal femur, primarily the anatomical geometry. Finite element analysis (FEA) is an effective and widely used computerbased simulation technique for modeling mechanical loading of various engineering structures, providing predictions of displacement and induced stress distribution due to the applied load. FEA is therefore inherently dependent upon both density and anatomical geometry. FEA may be performed on both three-dimensional and two-dimensional models of the proximal femur derived from radiographic images, from which the mechanical stiffness may be redicted. It is examined whether the outcome measures of two-dimensional FEA, two-dimensional, finite element analysis of X-ray images (FEXI), and three-dimensional FEA computed stiffness of the proximal femur were more sensitive than aBMD to changes in trabecular bone density and femur geometry. It is assumed that if an outcome measure follows known trends with changes in density and geometric parameters, then an increased sensitivity will be indicative of an improved prediction of bone strength. All three outcome measures increased non-linearly with trabecular bone density, increased linearly with cortical shell thickness and neck width, decreased linearly with neck length, and were relatively insensitive to neck-shaft angle. For femoral head radius, aBMD was relatively insensitive, with two-dimensional FEXI and threedimensional FEA demonstrating a non-linear increase and decrease in sensitivity, respectively. For neck anteversion, aBMD decreased non-linearly, whereas both two-dimensional FEXI and three dimensional FEA demonstrated a parabolic-type relationship, with maximum stiffness achieved at an angle of approximately 15o. Multi-parameter analysis showed that all three outcome measures demonstrated their highest sensitivity to a change in cortical thickness. When changes in all input parameters were considered simultaneously, three and twodimensional FEA had statistically equal sensitivities (0.41±0.20 and 0.42±0.16 respectively, p = ns) that were significantly higher than the sensitivity of aBMD (0.24±0.07; p = 0.014 and 0.002 for three-dimensional and two-dimensional FEA respectively). This simulation study suggests that since mechanical integrity and FEA are inherently dependent upon anatomical geometry, FEXI stiffness, being derived from conventional two-dimensional radiographic images, may provide an improvement in the prediction of bone strength of the proximal femur than currently provided by aBMD.

Mice lacking β-Adrenergic receptors have increased bone mass but are not protected from deleterious skeletal effects of ovariectomy

Activation of β2-adrenergic receptors inhibits osteoblastic bone formation and enhances osteoclastic bone resorption. Whether β-blockers inhibit ovariectomy-induced bone loss and decrease fracture risk remains controversial. To further explore the role of β-adrenergic signaling in skeletal acquisition and response to estrogen deficiency, we evaluated mice lacking the three known β-adrenergic receptors (β-less). Body weight, percent fat, and bone mineral density were significantly higher in male β-less than wild-type (WT) mice, more so with increasing age. Consistent with their greater fat mass, serum leptin was significantly higher in β-less than WT mice. Mid-femoral cross-sectional area and cortical thickness were significantly higher in adult β-less than WT mice, as were femoral biomechanical properties (+28 to +49%, P < 0.01). Young male β-less had higher vertebral (1.3-fold) and distal femoral (3.5-fold) trabecular bone volume than WT (P < 0.001 for both) and lower osteoclast surface. With aging, these differences lessened, with histological evidence of increased osteoclast surface and decreased bone formation rate at the distal femur in β-less vs. WT mice. Serum tartrate-resistance alkaline phosphatase-5B was elevated in β-less compared with WT mice from 8–16 wk of age (P < 0.01). Ovariectomy inhibited bone mass gain and decreased trabecular bone volume/total volume similarly in β-less and WT mice. Altogether, these data indicate that absence of β-adrenergic signaling results in obesity and increased cortical bone mass in males but does not prevent deleterious effects of estrogen deficiency on trabecular bone microarchitecture. Our findings also suggest direct positive effects of weight and/or leptin on bone turnover and cortical bone structure, independent of adrenergic signaling.

Intermittent Fugu parathyroid hormone 1 (1–34) is an anabolic bone agent in young male rats and osteopenic ovariectomized rats

Human parathyroid hormone (hPTH) is currently the only treatment for osteoporosis that forms new bone. Previously we described a fish equivalent, Fugu parathyroid hormone 1 (fPth1) which has hPTH-like biological activity in vitro despite fPth1(1–34) sharing only 53% identity with hPTH(1–34). Here we demonstrate the in vivo actions of fPth1(1–34) on bone. In study 1, young male rats were injected intermittently for 30 days with fPth1 [30 μg–1000 μg/kg body weight (b.w.), (30fPth1–1000fPth1)] or hPTH [30 μg–100 μg/kg b.w. (30hPTH–100hPTH)]. In proximal tibiae at low doses, the fPth1 was positively correlated with trabecular bone volume/total volume (TbBV/TV) while hPTH increased TbBV/TV, trabecular thickness (TbTh) and trabecular number (TbN). 500fPth1 and 1000fPth1 increased TbBV/TV, TbTh, TbN, mineral apposition rate (MAR) and bone formation rate/bone surface (BFR/BS) with a concomitant decrease in osteoclast surface and number. In study 2 ovariectomized (OVX), osteopenic rats and sham operated (SHAM) rats were injected intermittently with 500 μg/kg b.w. of fPth1 (500fPth1) for 11 weeks. 500fPth1 treatment resulted in increased TbBV/TV (151%) and TbTh (96%) in the proximal tibiae due to increased bone formation as assessed by BFR/BS (490%) and MAR (131%). The effect was restoration of TbBV/TV to SHAM levels without any effect on bone resorption. 500fPth1 also increased TbBV/TV and TbTh in the vertebrae (L6) and cortical thickness in the mid-femora increasing bone strength at these sites. fPth1 was similarly effective in SHAM rats. Notwithstanding the low amino acid sequence homology with hPTH (1–34), we have clearly established the efficacy of fPth1 (1–34) as an anabolic bone agent.

Proximal femoral morphology in regions of hard and soft water

This study compared proximal femoral morphology in patients living in soft and hard water regions. The proximal femoral morphology of two groups of 70 patients living in hard and soft water regions with a mean age of 72.3 (range 50 to 87 years) were measured using an antero-posterior radiograph of the non-operated hip with magnification adjusted. The medullary canal diameter at the level of the lesser trochanter (LT) was significantly wider in patients living in the hard water region (mean width 1.9 mm wider; p= 0.003). No statistical significant difference was found in the medullary canal width at 10 cm below the level of LT, Dorr index, or Canal Bone Ratio (CBR). In conclusion, the proximal femoral morphology does differ in patients living in soft and hard water areas. These results may have an important clinical bearing in patients undergoing total hip replacement surgery. Further research is needed to determine whether implant survivorship is affected in patients living in hard and soft water regions.

Mutations in known monogenic high bone mass loci only explain a small proportion of high bone mass cases

High bone mass (HBM) can be an incidental clinical finding; however, monogenic HBM disorders (eg, LRP5 or SOST mutations) are rare. We aimed to determine to what extent HBM is explained by mutations in known HBM genes. A total of 258 unrelated HBM cases were identified from a review of 335,115 DXA scans from 13 UK centers. Cases were assessed clinically and underwent sequencing of known anabolic HBM loci: LRP5 (exons 2, 3, 4), LRP4 (exons 25, 26), SOST (exons 1, 2, and the van Buchem's disease [VBD] 52-kb intronic deletion 3'). Family members were assessed for HBM segregation with identified variants. Three-dimensional protein models were constructed for identified variants. Two novel missense LRP5 HBM mutations ([c.518C>T; p.Thr173Met], [c.796C>T; p.Arg266Cys]) were identified, plus three previously reported missense LRP5 mutations ([c.593A>G; p.Asn198Ser], [c.724G>A; p.Ala242Thr], [c.266A>G; p.Gln89Arg]), associated with HBM in 11 adults from seven families. Individuals with LRP5 HBM ( approximately prevalence 5/100,000) displayed a variable phenotype of skeletal dysplasia with increased trabecular BMD and cortical thickness on HRpQCT, and gynoid fat mass accumulation on DXA, compared with both non-LRP5 HBM and controls. One mostly asymptomatic woman carried a novel heterozygous nonsense SOST mutation (c.530C>A; p.Ser177X) predicted to prematurely truncate sclerostin. Protein modeling suggests the severity of the LRP5-HBM phenotype corresponds to the degree of protein disruption and the consequent effect on SOST-LRP5 binding. We predict p.Asn198Ser and p.Ala242Thr directly disrupt SOST binding; both correspond to severe HBM phenotypes (BMD Z-scores +3.1 to +12.2, inability to float). Less disruptive structural alterations predicted from p.Arg266Cys, p.Thr173Met, and p.Gln89Arg were associated with less severe phenotypes (Z-scores +2.4 to +6.2, ability to float). In conclusion, although mutations in known HBM loci may be asymptomatic, they only account for a very small proportion ( approximately 3%) of HBM individuals, suggesting the great majority are explained by either unknown monogenic causes or polygenic inheritance.