INTRODUCTION: The cerebral cortex underlies our complex cognitive capabilities. Variations in human cortical surface area and thickness are associated with neurological, psychological, and behavioral traits and can be measured in vivo by magnetic resonance imaging (MRI). Studies in model organisms have identified genes that influence cortical structure, but little is known about common genetic variants that affect human cortical structure. RATIONALE: To identify genetic variants associated with human cortical structure at both global and regional levels, we conducted a genome-wide association meta-analysis of brain MRI data from 51,665 individuals across 60 cohorts. We analyzed the surface area and average thickness of the whole cortex and 34 cortical regions with known functional specializations. RESULTS: We identified 369 nominally genome-wide significant loci ( P < 5 × 10 ⁻⁸ ) associated with cortical structure in a discovery sample of 33,992 participants of European ancestry. Of the 360 loci for which replication data were available, 241 loci influencing surface area and 66 influencing thickness remained significant after replication, with 237 loci passing multiple testing correction ( P < 8.3 × 10 ⁻¹⁰ ; 187 influencing surface area and 50 influencing thickness). Common genetic variants explained 34% (SE = 3%) of the variation in total surface area and 26% (SE = 2%) in average thickness; surface area and thickness showed a negative genetic correlation ( r _G = −0.32, SE = 0.05, P = 6.5 × 10 ⁻¹² ), which suggests that genetic influences have opposing effects on surface area and thickness. Bioinformatic analyses showed that total surface area is influenced by genetic variants that alter gene regulatory activity in neural progenitor cells during fetal development. By contrast, average thickness is influenced by active regulatory elements in adult brain samples, which may reflect processes that occur after mid-fetal development, such as myelination, branching, or pruning. When considered together, these results support the radial unit hypothesis that different developmental mechanisms promote surface area expansion and increases in thickness. To identify specific genetic influences on individual cortical regions, we controlled for global measures (total surface area or average thickness) in the regional analyses. After multiple testing correction, we identified 175 loci that influence regional surface area and 46 that influence regional thickness. Loci that affect regional surface area cluster near genes involved in the Wnt signaling pathway, which is known to influence areal identity. We observed significant positive genetic correlations and evidence of bidirectional causation of total surface area with both general cognitive functioning and educational attainment. We found additional positive genetic correlations between total surface area and Parkinson’s disease but did not find evidence of causation. Negative genetic correlations were evident between total surface area and insomnia, attention deficit hyperactivity disorder, depressive symptoms, major depressive disorder, and neuroticism. CONCLUSION: This large-scale collaborative work enhances our understanding of the genetic architecture of the human cerebral cortex and its regional patterning. The highly polygenic architecture of the cortex suggests that distinct genes are involved in the development of specific cortical areas. Moreover, we find evidence that brain structure is a key phenotype along the causal pathway that leads from genetic variation to differences in general cognitive function. Identifying genetic influences on human cortical structure.: ( A ) Measurement of cortical surface area and thickness from MRI. ( B ) Genomic locations of common genetic variants that influence global and regional cortical structure. ( C ) Our results support the radial unit hypothesis that the expansion of cortical surface area is driven by proliferating neural progenitor cells. ( D ) Cortical surface area shows genetic correlation with psychiatric and cognitive traits. Error bars indicate SE. IMAGE CREDITS: (A) K. COURTNEY; (C) M. R. GLASS

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).

Schizophrenia (SZ) is associated with an increased risk of life-long cognitive impairments, age-related chronic disease, and premature mortality. We investigated evidence for advanced brain ageing in adult SZ patients, and whether this was associated with clinical characteristics in a prospective meta-analytic study conducted by the ENIGMA Schizophrenia Working Group. The study included data from 26 cohorts worldwide, with a total of 2803 SZ patients (mean age 34.2 years; range 18–72 years; 67% male) and 2598 healthy controls (mean age 33.8 years, range 18–73 years, 55% male). Brain-predicted age was individually estimated using a model trained on independent data based on 68 measures of cortical thickness and surface area, 7 subcortical volumes, lateral ventricular volumes and total intracranial volume, all derived from T1-weighted brain magnetic resonance imaging (MRI) scans. Deviations from a healthy brain ageing trajectory were assessed by the difference between brain-predicted age and chronological age (brain-predicted age difference [brain-PAD]). On average, SZ patients showed a higher brain-PAD of +3.55 years (95% CI: 2.91, 4.19; I² = 57.53%) compared to controls, after adjusting for age, sex and site (Cohen’s d = 0.48). Among SZ patients, brain-PAD was not associated with specific clinical characteristics (age of onset, duration of illness, symptom severity, or antipsychotic use and dose). This large-scale collaborative study suggests advanced structural brain ageing in SZ. Longitudinal studies of SZ and a range of mental and somatic health outcomes will help to further evaluate the clinical implications of increased brain-PAD and its ability to be influenced by interventions.