The hippocampus consists of anatomically and functionally distinct subfields that may be differentially involved in the pathophysiology of bipolar disorder (BD). Here we, the Enhancing NeuroImaging Genetics through Meta-Analysis Bipolar Disorder workinggroup, study hippocampal subfield volumetry in BD. T1-weighted magnetic resonance imaging scans from 4,698 individuals (BD = 1,472, healthy controls [HC] = 3,226) from 23 sites worldwide were processed with FreeSurfer. We used linear mixed-effects models and mega-analysis to investigate differences in hippocampal subfield volumes between BD and HC, followed by analyses of clinical characteristics and medication use. BD showed significantly smaller volumes of the whole hippocampus (Cohen's d = -0.20), cornu ammonis (CA)1 (d = -0.18), CA2/3 (d = -0.11), CA4 (d = -0.19), molecular layer (d = -0.21), granule cell layer of dentate gyrus (d = -0.21), hippocampal tail (d = -0.10), subiculum (d = -0.15), presubiculum (d = -0.18), and hippocampal amygdala transition area (d = -0.17) compared to HC. Lithium users did not show volume differences compared to HC, while non-users did. Antipsychotics or antiepileptic use was associated with smaller volumes. In this largest study of hippocampal subfields in BD to date, we show widespread reductions in nine of 12 subfields studied. The associations were modulated by medication use and specifically the lack of differences between lithium users and HC supports a possible protective role of lithium in BD.

Genomewide association studies have found significant genetic correlations among many neuropsychiatric disorders. In contrast, we know much less about the degree to which structural brain alterations are similar among disorders and, if so, the degree to which such similarities have a genetic etiology. From the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) consortium, we acquired standardized mean differences (SMDs) in regional brain volume and cortical thickness between cases and controls. We had data on 41 brain regions for: attention-deficit/hyperactivity disorder (ADHD), autism spectrum disorder (ASD), bipolar disorder (BD), epilepsy, major depressive disorder (MDD), obsessive compulsive disorder (OCD), and schizophrenia (SCZ). These data had been derived from 24,360 patients and 37,425 controls. The SMDs were significantly correlated between SCZ and BD, OCD, MDD, and ASD. MDD was positively correlated with BD and OCD. BD was positively correlated with OCD and negatively correlated with ADHD. These pairwise correlations among disorders were correlated with the corresponding pairwise correlations among disorders derived from genomewide association studies (r = 0.494). Our results show substantial similarities in sMRI phenotypes among neuropsychiatric disorders and suggest that these similarities are accounted for, in part, by corresponding similarities in common genetic variant architectures.

A common limitation of neuroimaging studies is their small sample sizes. To overcome this hurdle, the Enhancing Neuro Imaging Genetics through Meta-Analysis (ENIGMA) Consortium combines neuroimaging data from many institutions worldwide. However, this introduces heterogeneity due to different scanning devices and sequences. ENIGMA projects commonly address this heterogeneity with random-effects meta-analysis or mixed-effects mega-analysis. Here we tested whether the batch adjustment method, ComBat, can further reduce site-related heterogeneity and thus increase statistical power. We conducted random-effects meta-analyses, mixed-effects mega-analyses and ComBat mega-analyses to compare cortical thickness, surface area and subcortical volumes between 2897 individuals with a diagnosis of schizophrenia and 3141 healthy controls from 33 sites. Specifically, we compared the imaging data between individuals with schizophrenia and healthy controls, covarying for age and sex. The use of ComBat substantially increased the statistical significance of the findings as compared to random-effects meta-analyses. The findings were more similar when comparing ComBat with mixed-effects mega-analysis, although ComBat still slightly increased the statistical significance. ComBat also showed increased statistical power when we repeated the analyses with fewer sites. Results were nearly identical when we applied the ComBat harmonization separately for cortical thickness, cortical surface area and subcortical volumes. Therefore, we recommend applying the ComBat function to attenuate potential effects of site in ENIGMA projects and other multi-site structural imaging work. We provide easy-to-use functions in R that work even if imaging data are partially missing in some brain regions, and they can be trained with one data set and then applied to another (a requirement for some analyses such as machine learning).

In a recent machine learning study classifying "brain age" based on cross-sectional neuroanatomical data, clinical high-risk (CHR) individuals were observed to show deviation from the normal neuromaturational pattern, which in turn was predictive of greater risk of conversion to psychosis and a pattern of stably poor functional outcome. These effects were unique to cases who were between 12 and 17 years of age when their prodromal and psychotic symptoms began, suggesting that neuroanatomical deviance observable at the point of ascertainment of a CHR syndrome marks risk for an early onset form of psychosis. In the present study, we sought to clarify the pattern of neuroanatomical deviance linked to this "early onset" form of psychosis and whether this deviance is associated with poorer premorbid functioning. T1 MRI scans from 378 CHR individuals and 190 healthy controls (HC) from the North American Prodrome Longitudinal Study (NAPLS2) were analyzed. Widespread smaller cortical volume was observed among CHR individuals compared with HC at baseline evaluation, particularly among the younger group (i.e., those who were 12 to 17 years of age). Moreover, the younger CHR individuals who converted or presented worsened clinical symptoms at follow-up (within 2 years) exhibited smaller surface area in rostral anterior cingulate, lateral and medial prefrontal regions, and parahippocampal gyrus relative to the younger CHR individuals who remitted or presented a stable pattern of prodromal symptoms at follow-up. In turn, poorer premorbid functioning in childhood was associated with smaller surface area in medial orbitofrontal, lateral frontal, rostral anterior cingulate, precuneus, and temporal regions. Together with our prior report, these results are consistent with the view that neuroanatomical deviance manifesting in early adolescence marks vulnerability to a form of psychosis presenting with poor premorbid adjustment, an earlier age of onset (generally prior to the age of 18 years), and poor long-term outcome.

Alzheimer's disease (AD) is characterized by cortical atrophy and disrupted anatomic connectivity, and leads to abnormal interactions between neural systems. Diffusion-weighted imaging (DWI) and graph theory can be used to evaluate major brain networks and detect signs of a breakdown in network connectivity. In a longitudinal study using both DWI and standard magnetic resonance imaging (MRI), we assessed baseline white-matter connectivity patterns in 30 subjects with mild cognitive impairment (MCI, mean age 71.8 ± 7.5 years, 18 males and 12 females) from the Alzheimer's Disease Neuroimaging Initiative. Using both standard MRI-based cortical parcellations and whole-brain tractography, we computed baseline connectivity maps from which we calculated global "small-world" architecture measures, including mean clustering coefficient and characteristic path length. We evaluated whether these baseline network measures predicted future volumetric brain atrophy in MCI subjects, who are at risk for developing AD, as determined by 3-dimensional Jacobian "expansion factor maps" between baseline and 6-month follow-up anatomic scans. This study suggests that DWI-based network measures may be a novel predictor of AD progression.

Brain structural covariance networks reflect covariation in morphology of different brain areas and are thought to reflect common trajectories in brain development and maturation. Large-scale investigation of structural covariance networks in obsessive-compulsive disorder (OCD) may provide clues to the pathophysiology of this neurodevelopmental disorder. Using T1-weighted MRI scans acquired from 1616 individuals with OCD and 1463 healthy controls across 37 datasets participating in the ENIGMA-OCD Working Group, we calculated intra-individual brain structural covariance networks (using the bilaterally-averaged values of 33 cortical surface areas, 33 cortical thickness values, and six subcortical volumes), in which edge weights were proportional to the similarity between two brain morphological features in terms of deviation from healthy controls (i.e. z-score transformed). Global networks were characterized using measures of network segregation (clustering and modularity), network integration (global efficiency), and their balance (small-worldness), and their community membership was assessed. Hub profiling of regional networks was undertaken using measures of betweenness, closeness, and eigenvector centrality. Individually calculated network measures were integrated across the 37 datasets using a meta-analytical approach. These network measures were summated across the network density range of K = 0.10-0.25 per participant, and were integrated across the 37 datasets using a meta-analytical approach. Compared with healthy controls, at a global level, the structural covariance networks of OCD showed lower clustering (P < 0.0001), lower modularity (P < 0.0001), and lower small-worldness (P = 0.017). Detection of community membership emphasized lower network segregation in OCD compared to healthy controls. At the regional level, there were lower (rank-transformed) centrality values in OCD for volume of caudate nucleus and thalamus, and surface area of paracentral cortex, indicative of altered distribution of brain hubs. Centrality of cingulate and orbito-frontal as well as other brain areas was associated with OCD illness duration, suggesting greater involvement of these brain areas with illness chronicity. In summary, the findings of this study, the largest brain structural covariance study of OCD to date, point to a less segregated organization of structural covariance networks in OCD, and reorganization of brain hubs. The segregation findings suggest a possible signature of altered brain morphometry in OCD, while the hub findings point to OCD-related alterations in trajectories of brain development and maturation, particularly in cingulate and orbitofrontal regions.

Carriers of large recurrent copy number variants (CNVs) have a higher risk of developing neurodevelopmental disorders. The 16p11.2 distal CNV predisposes carriers to e.g., autism spectrum disorder and schizophrenia. We compared subcortical brain volumes of 12 16p11.2 distal deletion and 12 duplication carriers to 6882 non-carriers from the large-scale brain Magnetic Resonance Imaging collaboration, ENIGMA-CNV. After stringent CNV calling procedures, and standardized FreeSurfer image analysis, we found negative dose-response associations with copy number on intracranial volume and on regional caudate, pallidum and putamen volumes (β = -0.71 to -1.37; P < 0.0005). In an independent sample, consistent results were obtained, with significant effects in the pallidum (β = -0.95, P = 0.0042). The two data sets combined showed significant negative dose-response for the accumbens, caudate, pallidum, putamen and ICV (P = 0.0032, 8.9 × 10-6, 1.7 × 10-9, 3.5 × 10-12 and 1.0 × 10-4, respectively). Full scale IQ was lower in both deletion and duplication carriers compared to non-carriers. This is the first brain MRI study of the impact of the 16p11.2 distal CNV, and we demonstrate a specific effect on subcortical brain structures, suggesting a neuropathological pattern underlying the neurodevelopmental syndromes.

Altered structural brain asymmetry in autism spectrum disorder (ASD) has been reported. However, findings have been inconsistent, likely due to limited sample sizes. Here we investigated 1,774 individuals with ASD and 1,809 controls, from 54 independent data sets of the ENIGMA consortium. ASD was significantly associated with alterations of cortical thickness asymmetry in mostly medial frontal, orbitofrontal, cingulate and inferior temporal areas, and also with asymmetry of orbitofrontal surface area. These differences generally involved reduced asymmetry in individuals with ASD compared to controls. Furthermore, putamen volume asymmetry was significantly increased in ASD. The largest case-control effect size was Cohen's d = -0.13, for asymmetry of superior frontal cortical thickness. Most effects did not depend on age, sex, IQ, severity or medication use. Altered lateralized neurodevelopment may therefore be a feature of ASD, affecting widespread brain regions with diverse functions. Large-scale analysis was necessary to quantify subtle alterations of brain structural asymmetry in ASD.

Small average differences in the left-right asymmetry of cerebral cortical thickness have been reported in individuals with autism spectrum disorder (ASD) compared to typically developing controls, affecting widespread cortical regions. The possible impacts of these regional alterations in terms of structural network effects have not previously been characterized. Inter-regional morphological covariance analysis can capture network connectivity between different cortical areas at the macroscale level. Here, we used cortical thickness data from 1455 individuals with ASD and 1560 controls, across 43 independent datasets of the ENIGMA consortium's ASD Working Group, to assess hemispheric asymmetries of intra-individual structural covariance networks, using graph theory-based topological metrics. Compared with typical features of small-world architecture in controls, the ASD sample showed significantly altered average asymmetry of networks involving the fusiform, rostral middle frontal, and medial orbitofrontal cortex, involving higher randomization of the corresponding right-hemispheric networks in ASD. A network involving the superior frontal cortex showed decreased right-hemisphere randomization. Based on comparisons with meta-analyzed functional neuroimaging data, the altered connectivity asymmetry particularly affected networks that subserve executive functions, language-related and sensorimotor processes. These findings provide a network-level characterization of altered left-right brain asymmetry in ASD, based on a large combined sample. Altered asymmetrical brain development in ASD may be partly propagated among spatially distant regions through structural connectivity.

Understanding the fundamental alterations in brain functioning that lead to psychotic disorders remains a major challenge in clinical neuroscience. In particular, it is unknown whether any state-independent biomarkers can potentially predict the onset of psychosis and distinguish patients from healthy controls, regardless of paradigm. Here, using multi-paradigm fMRI data from the North American Prodrome Longitudinal Study consortium, we show that individuals at clinical high risk for psychosis display an intrinsic "trait-like" abnormality in brain architecture characterized as increased connectivity in the cerebello-thalamo-cortical circuitry, a pattern that is significantly more pronounced among converters compared with non-converters. This alteration is significantly correlated with disorganization symptoms and predictive of time to conversion to psychosis. Moreover, using an independent clinical sample, we demonstrate that this hyperconnectivity pattern is reliably detected and specifically present in patients with schizophrenia. These findings implicate cerebello-thalamo-cortical hyperconnectivity as a robust state-independent neural signature for psychosis prediction and characterization.

Alterations in subcortical brain structures have been reported in adults with HIV and, to a lesser extent, pediatric cohorts. The extent of longitudinal structural abnormalities in children with perinatal HIV infection (PaHIV) remains unclear. We modeled subcortical morphometry from whole brain structural magnetic resonance imaging (1.5 T) scans of 43 Thai children with PaHIV (baseline age = 11.09±2.36 years) and 50 HIV- children (11.26±2.80 years) using volumetric and surface-based shape analyses. The PaHIV sample were randomized to initiate combination antiretroviral treatment (cART) when CD4 counts were 15-24% (immediate: n = 22) or when CD4 < 15% (deferred: n = 21). Follow-up scans were acquired approximately 52 weeks after baseline. Volumetric and shape descriptors capturing local thickness and surface area dilation were defined for the bilateral accumbens, amygdala, putamen, pallidum, thalamus, caudate, and hippocampus. Regression models adjusting for clinical and demographic variables examined between and within group differences in morphometry associated with HIV. We assessed whether baseline CD4 count and cART status or timing associated with brain maturation within the PaHIV group. All models were adjusted for multiple comparisons using the false discovery rate. A pallidal subregion was significantly thinner in children with PaHIV. Regional thickness, surface area, and volume of the pallidum was associated with CD4 count in children with PaHIV. Longitudinal morphometry was not associated with HIV or cART status or timing, however, the trajectory of the left pallidum volume was positively associated with baseline CD4 count. Our findings corroborate reports in adult cohorts demonstrating a high predilection for HIV-mediated abnormalities in the basal ganglia, but suggest the effect of stable PaHIV infection on morphological aspects of brain development may be subtle.

Machine learning can be used to define subtypes of psychiatric conditions based on shared clinical and biological foundations, presenting a crucial step toward establishing biologically based subtypes of mental disorders. With the goal of identifying subtypes of disease progression in schizophrenia, here we analyzed cross-sectional brain structural magnetic resonance imaging (MRI) data from 4,291 individuals with schizophrenia (1,709 females, age=32.5 years±11.9) and 7,078 healthy controls (3,461 females, age=33.0 years±12.7) pooled across 41 international cohorts from the ENIGMA Schizophrenia Working Group, non-ENIGMA cohorts and public datasets. Using a machine learning approach known as Subtype and Stage Inference (SuStaIn), we implemented a brain imaging-driven classification that identifies two distinct neurostructural subgroups by mapping the spatial and temporal trajectory of gray matter (GM) loss in schizophrenia. Subgroup 1 (n=2,622) was characterized by an early cortical-predominant loss (ECL) with enlarged striatum, whereas subgroup 2 (n=1,600) displayed an early subcortical-predominant loss (ESL) in the hippocampus, amygdala, thalamus, brain stem and striatum. These reconstructed trajectories suggest that the GM volume reduction originates in the Broca's area/adjacent fronto-insular cortex for ECL and in the hippocampus/adjacent medial temporal structures for ESL. With longer disease duration, the ECL subtype exhibited a gradual worsening of negative symptoms and depression/anxiety, and less of a decline in positive symptoms. We confirmed the reproducibility of these imaging-based subtypes across various sample sites, independent of macroeconomic and ethnic factors that differed across these geographic locations, which include Europe, North America and East Asia. These findings underscore the presence of distinct pathobiological foundations underlying schizophrenia. This new imaging-based taxonomy holds the potential to identify a more homogeneous sub-population of individuals with shared neurobiological attributes, thereby suggesting the viability of redefining existing disorder constructs based on biological factors.

Alzheimer's disease (AD) involves a gradual breakdown of brain connectivity, and network analyses offer a promising new approach to track and understand disease progression. Even so, our ability to detect degenerative changes in brain networks depends on the methods used. Here we compared several tractography and feature extraction methods to see which ones gave best diagnostic classification for 202 people with AD, mild cognitive impairment or normal cognition, scanned with 41-gradient diffusion-weighted magnetic resonance imaging as part of the Alzheimer's Disease Neuroimaging Initiative (ADNI) project. We computed brain networks based on whole brain tractography with nine different methods - four of them tensor-based deterministic (FACT, RK2, SL, and TL), two orientation distribution function (ODF)-based deterministic (FACT, RK2), two ODF-based probabilistic approaches (Hough and PICo), and one "ball-and-stick" approach (Probtrackx). Brain networks derived from different tractography algorithms did not differ in terms of classification performance on ADNI, but performing principal components analysis on networks helped classification in some cases. Small differences may still be detectable in a truly vast cohort, but these experiments help assess the relative advantages of different tractography algorithms, and different post-processing choices, when used for classification.

The goal of this study was to compare brain structure between individuals with generalized anxiety disorder (GAD) and healthy controls. Previous studies have generated inconsistent findings, possibly due to small sample sizes, or clinical/analytic heterogeneity. To address these concerns, we combined data from 28 research sites worldwide through the ENIGMA-Anxiety Working Group, using a single, pre-registered mega-analysis. Structural magnetic resonance imaging data from children and adults (5-90 years) were processed using FreeSurfer. The main analysis included the regional and vertex-wise cortical thickness, cortical surface area, and subcortical volume as dependent variables, and GAD, age, age-squared, sex, and their interactions as independent variables. Nuisance variables included IQ, years of education, medication use, comorbidities, and global brain measures. The main analysis (1020 individuals with GAD and 2999 healthy controls) included random slopes per site and random intercepts per scanner. A secondary analysis (1112 individuals with GAD and 3282 healthy controls) included fixed slopes and random intercepts per scanner with the same variables. The main analysis showed no effect of GAD on brain structure, nor interactions involving GAD, age, or sex. The secondary analysis showed increased volume in the right ventral diencephalon in male individuals with GAD compared to male healthy controls, whereas female individuals with GAD did not differ from female healthy controls. This mega-analysis combining worldwide data showed that differences in brain structure related to GAD are small, possibly reflecting heterogeneity or those structural alterations are not a major component of its pathophysiology.

Diffusion MRI (dMRI) can be used to probe microstructural properties of brain tissue and holds great promise as a means to non-invasively map Alzheimer's disease (AD) pathology. Few studies have evaluated multi-shell dMRI models, such as neurite orientation dispersion and density imaging (NODDI) and mean apparent propagator (MAP)-MRI, in cortical gray matter where many of the earliest histopathological changes occur in AD. Here, we investigated the relationship between CSF pTau 181 and Aβ 1-42 burden and regional cortical NODDI and MAP-MRI indices in 46 cognitively unimpaired individuals, 18 with mild cognitive impairment, and two with dementia (mean age: 71.8±6.2 years) from the Alzheimer's Disease Neuroimaging Initiative. We compared findings to more conventional cortical thickness measures. Lower CSF Aβ 1-42 and higher pTau 181 were associated with cortical dMRI measures reflecting less hindered or restricted diffusion and greater diffusivity. Cortical dMRI measures were more widely associated with Aβ 1-42 than pTau 181 and better distinguished Aβ+ from Aβ-participants than pTau+/- participants. Conversely, cortical thickness was more tightly linked with pTau 181 . dMRI associations mediated the relationship between CSF markers and delayed logical memory performance, commonly impaired in early AD. dMRI measures sensitive to early AD pathogenesis and microstructural damage may elucidate mechanisms underlying cognitive decline.

Modern non-invasive brain imaging technologies, such as diffusion weighted magnetic resonance imaging (DWI), enable the mapping of neural fiber tracts in the white matter, providing a basis to reconstruct a detailed map of brain structural connectivity networks. Brain connectivity networks differ from random networks in their topology, which can be measured using small worldness, modularity, and high-degree nodes (hubs). Still, little is known about how individual differences in structural brain network properties relate to age, sex, or genetic differences. Recently, some groups have reported brain network biomarkers that enable differentiation among individuals, pairs of individuals, and groups of individuals. In addition to studying new topological features, here we provide a unifying general method to investigate topological brain networks and connectivity differences between individuals, pairs of individuals, and groups of individuals at several levels of the data hierarchy, while appropriately controlling false discovery rate (FDR) errors. We apply our new method to a large dataset of high quality brain connectivity networks obtained from High Angular Resolution Diffusion Imaging (HARDI) tractography in 303 young adult twins, siblings, and unrelated people. Our proposed approach can accurately classify brain connectivity networks based on sex (93% accuracy) and kinship (88.5% accuracy). We find statistically significant differences associated with sex and kinship both in the brain connectivity networks and in derived topological metrics, such as the clustering coefficient and the communicability matrix.

It is increasingly recognized that multiple psychiatric conditions are underpinned by shared neural pathways, affecting similar brain systems. Here, we carried out a multiscale neural contextualization of shared alterations of cortical morphology across six major psychiatric conditions (autism spectrum disorder, attention deficit/hyperactivity disorder, major depression disorder, obsessive-compulsive disorder, bipolar disorder, and schizophrenia). Our framework cross-referenced shared morphological anomalies with respect to cortical myeloarchitecture and cytoarchitecture, as well as connectome and neurotransmitter organization. Pooling disease-related effects on MRI-based cortical thickness measures across six ENIGMA working groups, including a total of 28,546 participants (12,876 patients and 15,670 controls), we identified a cortex-wide dimension of morphological changes that described a sensory-fugal pattern, with paralimbic regions showing the most consistent alterations across conditions. The shared disease dimension was closely related to cortical gradients of microstructure as well as neurotransmitter axes, specifically cortex-wide variations in serotonin and dopamine. Multiple sensitivity analyses confirmed robustness with respect to slight variations in analytical choices. Our findings embed shared effects of common psychiatric conditions on brain structure in multiple scales of brain organization, and may provide insights into neural mechanisms of transdiagnostic vulnerability.

Left-right asymmetry is an important organizing feature of the healthy brain that may be altered in schizophrenia, but most studies have used relatively small samples and heterogeneous approaches, resulting in equivocal findings. We carried out the largest case-control study of structural brain asymmetries in schizophrenia, with MRI data from 5,080 affected individuals and 6,015 controls across 46 datasets, using a single image analysis protocol. Asymmetry indexes were calculated for global and regional cortical thickness, surface area, and subcortical volume measures. Differences of asymmetry were calculated between affected individuals and controls per dataset, and effect sizes were meta-analyzed across datasets. Small average case-control differences were observed for thickness asymmetries of the rostral anterior cingulate and the middle temporal gyrus, both driven by thinner left-hemispheric cortices in schizophrenia. Analyses of these asymmetries with respect to the use of antipsychotic medication and other clinical variables did not show any significant associations. Assessment of age- and sex-specific effects revealed a stronger average leftward asymmetry of pallidum volume between older cases and controls. Case-control differences in a multivariate context were assessed in a subset of the data (N = 2,029), which revealed that 7% of the variance across all structural asymmetries was explained by case-control status. Subtle case-control differences of brain macrostructural asymmetry may reflect differences at the molecular, cytoarchitectonic, or circuit levels that have functional relevance for the disorder. Reduced left middle temporal cortical thickness is consistent with altered left-hemisphere language network organization in schizophrenia.

Brain lobar volumes are heritable but genetic studies are limited. We performed genome-wide association studies of frontal, occipital, parietal and temporal lobe volumes in 16,016 individuals, and replicated our findings in 8,789 individuals. We identified six genetic loci associated with specific lobar volumes independent of intracranial volume. Two loci, associated with occipital (6q22.32) and temporal lobe volume (12q14.3), were previously reported to associate with intracranial and hippocampal volume, respectively. We identified four loci previously unknown to affect brain volumes: 3q24 for parietal lobe volume, and 1q22, 4p16.3 and 14q23.1 for occipital lobe volume. The associated variants were located in regions enriched for histone modifications (DAAM1 and THBS3), or close to genes causing Mendelian brain-related diseases (ZIC4 and FGFRL1). No genetic overlap between lobar volumes and neurological or psychiatric diseases was observed. Our findings reveal part of the complex genetics underlying brain development and suggest a role for regulatory regions in determining brain volumes.

A commonly carried C677T polymorphism in a folate-related gene, MTHFR, is associated with higher plasma homocysteine, a well-known mediator of neuronal damage and brain atrophy. As homocysteine promotes brain atrophy, we set out to discover whether people carrying the C677T MTHFR polymorphism which increases homocysteine, might also show systematic differences in brain structure. Using tensor-based morphometry, we tested this association in 359 elderly Caucasian subjects with mild cognitive impairment (MCI) (mean age: 75 ± 7.1 years) scanned with brain MRI and genotyped as part of Alzheimer's Disease Neuroimaging Initiative. We carried out a replication study in an independent, non-overlapping sample of 51 elderly Caucasian subjects with MCI (mean age: 76 ± 5.5 years), scanned with brain MRI and genotyped for MTHFR, as part of the Cardiovascular Health Study. At each voxel in the brain, we tested to see where regional volume differences were associated with carrying one or more MTHFR 'T' alleles. In ADNI subjects, carriers of the MTHFR risk allele had detectable brain volume deficits, in the white matter, of up to 2-8% per risk T allele locally at baseline and showed accelerated brain atrophy of 0.5-1.5% per T allele at 1 year follow-up, after adjusting for age and sex. We replicated these brain volume deficits of up to 5-12% per MTHFR T allele in the independent cohort of CHS subjects. As expected, the associations weakened after controlling for homocysteine levels, which the risk gene affects. The MTHFR risk variant may thus promote brain atrophy by elevating homocysteine levels. This study aims to investigate the spatially detailed effects of this MTHFR polymorphism on brain structure in 3D, pointing to a causal pathway that may promote homocysteine-mediated brain atrophy in elderly people with MCI.