Effects of gene dosage and development on subcortical nuclei volumes in individuals with 22q11.2 copy number variations

The 22q11.2 locus contains genes critical for brain development. Reciprocal Copy Number Variations (CNVs) at this locus impact risk for neurodevelopmental and psychiatric disorders. Both 22q11.2 deletions (22qDel) and duplications (22qDup) are associated with autism, but 22qDel uniquely elevates schizophrenia risk. Understanding brain phenotypes associated with these highly penetrant CNVs can provide insights into genetic pathways underlying neuropsychiatric disorders. Human neuroimaging and animal models indicate subcortical brain alterations in 22qDel, yet little is known about developmental differences across specific nuclei between reciprocal 22q11.2 CNV carriers and typically developing (TD) controls. We conducted a longitudinal MRI study in a total of 385 scans from 22qDel (n = 96, scans = 191, 53.1% female), 22qDup (n = 37, scans = 64, 45.9% female), and TD controls (n = 80, scans = 130, 51.2% female), across a wide age range (5.5-49.5 years). Volumes of the thalamus, hippocampus, amygdala, and anatomical subregions were estimated using FreeSurfer, and the linear effects of 22q11.2 gene dosage and non-linear effects of age were characterized with generalized additive mixed models (GAMMs). Positive gene dosage effects (volume increasing with copy number) were observed for total intracranial and whole hippocampus volumes, but not whole thalamus or amygdala volumes. Several amygdala subregions exhibited similar positive effects, with bi-directional effects found across thalamic nuclei. Distinct age-related trajectories were observed across the three groups. Notably, both 22qDel and 22qDup carriers exhibited flattened development of hippocampal CA2/3 subfields relative to TD controls. This study provides novel insights into the impact of 22q11.2 CNVs on subcortical brain structures and their developmental trajectories.

Large-Scale Neuroimaging of Mental Illness

Neuroimaging has provided important insights into the brain variations related to mental illness. Inconsistencies in prior studies, however, call for methods that lead to more replicable and generalizable brain markers that can reliably predict illness severity, treatment course, and prognosis. A paradigm shift is underway with large-scale international research teams actively pooling data and resources to drive consensus findings and test emerging methods aimed at achieving the goals of precision psychiatry. In parallel with large-scale psychiatric genomics studies, international consortia combining neuroimaging data are mapping the transdiagnostic brain signatures of mental illness on an unprecedented scale. This chapter discusses the major challenges, recent findings, and a roadmap for developing better neuroimaging-based tools and markers for mental illness.

Examining the association between posttraumatic stress disorder and disruptions in cortical networks identified using data-driven methods

Posttraumatic stress disorder (PTSD) is associated with lower cortical thickness (CT) in prefrontal, cingulate, and insular cortices in diverse trauma-affected samples. However, some studies have failed to detect differences between PTSD patients and healthy controls or reported that PTSD is associated with greater CT. Using data-driven dimensionality reduction, we sought to conduct a well-powered study to identify vulnerable networks without regard to neuroanatomic boundaries. Moreover, this approach enabled us to avoid the excessive burden of multiple comparison correction that plagues vertex-wise methods. We derived structural covariance networks (SCNs) by applying non-negative matrix factorization (NMF) to CT data from 961 PTSD patients and 1124 trauma-exposed controls without PTSD. We used regression analyses to investigate associations between CT within SCNs and PTSD diagnosis (with and without accounting for the potential confounding effect of trauma type) and symptom severity in the full sample. We performed additional regression analyses in subsets of the data to examine associations between SCNs and comorbid depression, childhood trauma severity, and alcohol abuse. NMF identified 20 unbiased SCNs, which aligned closely with functionally defined brain networks. PTSD diagnosis was most strongly associated with diminished CT in SCNs that encompassed the bilateral superior frontal cortex, motor cortex, insular cortex, orbitofrontal cortex, medial occipital cortex, anterior cingulate cortex, and posterior cingulate cortex. CT in these networks was significantly negatively correlated with PTSD symptom severity. Collectively, these findings suggest that PTSD diagnosis is associated with widespread reductions in CT, particularly within prefrontal regulatory regions and broader emotion and sensory processing cortical regions.

Beyond the Global Brain Differences: Intraindividual Variability Differences in 1q21.1 Distal and 15q11.2 BP1-BP2 Deletion Carriers

BACKGROUND

Carriers of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants exhibit regional and global brain differences compared with noncarriers. However, interpreting regional differences is challenging if a global difference drives the regional brain differences. Intraindividual variability measures can be used to test for regional differences beyond global differences in brain structure.

METHODS

Magnetic resonance imaging data were used to obtain regional brain values for 1q21.1 distal deletion (n = 30) and duplication (n = 27) and 15q11.2 BP1-BP2 deletion (n = 170) and duplication (n = 243) carriers and matched noncarriers (n = 2350). Regional intra-deviation scores, i.e., the standardized difference between an individual's regional difference and global difference, were used to test for regional differences that diverge from the global difference.

RESULTS

For the 1q21.1 distal deletion carriers, cortical surface area for regions in the medial visual cortex, posterior cingulate, and temporal pole differed less and regions in the prefrontal and superior temporal cortex differed more than the global difference in cortical surface area. For the 15q11.2 BP1-BP2 deletion carriers, cortical thickness in regions in the medial visual cortex, auditory cortex, and temporal pole differed less and the prefrontal and somatosensory cortex differed more than the global difference in cortical thickness.

CONCLUSIONS

We find evidence for regional effects beyond differences in global brain measures in 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants. The results provide new insight into brain profiling of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants, with the potential to increase understanding of the mechanisms involved in altered neurodevelopment.

Multi-site benchmark classification of major depressive disorder using machine learning on cortical and subcortical measures

Machine learning (ML) techniques have gained popularity in the neuroimaging field due to their potential for classifying neuropsychiatric disorders. However, the diagnostic predictive power of the existing algorithms has been limited by small sample sizes, lack of representativeness, data leakage, and/or overfitting. Here, we overcome these limitations with the largest multi-site sample size to date (N = 5365) to provide a generalizable ML classification benchmark of major depressive disorder (MDD) using shallow linear and non-linear models. Leveraging brain measures from standardized ENIGMA analysis pipelines in FreeSurfer, we were able to classify MDD versus healthy controls (HC) with a balanced accuracy of around 62%. But after harmonizing the data, e.g., using ComBat, the balanced accuracy dropped to approximately 52%. Accuracy results close to random chance levels were also observed in stratified groups according to age of onset, antidepressant use, number of episodes and sex. Future studies incorporating higher dimensional brain imaging/phenotype features, and/or using more advanced machine and deep learning methods may yield more encouraging prospects.

Source-based morphometry reveals structural brain pattern abnormalities in 22q11.2 deletion syndrome

22q11.2 deletion syndrome (22q11DS) is the most frequently occurring microdeletion in humans. It is associated with a significant impact on brain structure, including prominent reductions in gray matter volume (GMV), and neuropsychiatric manifestations, including cognitive impairment and psychosis. It is unclear whether GMV alterations in 22q11DS occur according to distinct structural patterns. Then, 783 participants (470 with 22q11DS: 51% females, mean age [SD] 18.2 [9.2]; and 313 typically developing [TD] controls: 46% females, mean age 18.0 [8.6]) from 13 datasets were included in the present study. We segmented structural T1-weighted brain MRI scans and extracted GMV images, which were then utilized in a novel source-based morphometry (SBM) pipeline (SS-Detect) to generate structural brain patterns (SBPs) that capture co-varying GMV. We investigated the impact of the 22q11.2 deletion, deletion size, intelligence quotient, and psychosis on the SBPs. Seventeen GMV-SBPs were derived, which provided spatial patterns of GMV covariance associated with a quantitative metric (i.e., loading score) for analysis. Patterns of topographically widespread differences in GMV covariance, including the cerebellum, discriminated individuals with 22q11DS from healthy controls. The spatial extents of the SBPs that revealed disparities between individuals with 22q11DS and controls were consistent with the findings of the univariate voxel-based morphometry analysis. Larger deletion size was associated with significantly lower GMV in frontal and occipital SBPs; however, history of psychosis did not show a strong relationship with these covariance patterns. 22q11DS is associated with distinct structural abnormalities captured by topographical GMV covariance patterns that include the cerebellum. Findings indicate that structural anomalies in 22q11DS manifest in a nonrandom manner and in distinct covarying anatomical patterns, rather than a diffuse global process. These SBP abnormalities converge with previously reported cortical surface area abnormalities, suggesting disturbances of early neurodevelopment as the most likely underlying mechanism.

Effects of Gene Dosage and Development on Subcortical Nuclei Volumes in Individuals with 22q11.2 Copy Number Variations

The 22q11.2 locus contains genes critical for brain development. Reciprocal Copy Number Variations (CNVs) at this locus impact risk for neurodevelopmental and psychiatric disorders. Both 22q11.2 deletions (22qDel) and duplications (22qDup) are associated with autism, but 22qDel uniquely elevates schizophrenia risk. Understanding brain phenotypes associated with these highly penetrant CNVs can provide insights into genetic pathways underlying neuropsychiatric disorders. Human neuroimaging and animal models indicate subcortical brain alterations in 22qDel, yet little is known about developmental differences across specific nuclei between reciprocal 22q11.2 CNV carriers and typically developing (TD) controls. We conducted a longitudinal MRI study in 22qDel (n=96, 53.1% female), 22qDup (n=37, 45.9% female), and TD controls (n=80, 51.2% female), across a wide age range (5.5-49.5 years). Volumes of the thalamus, hippocampus, amygdala, and anatomical subregions were estimated using FreeSurfer, and the effect of 22q11.2 gene dosage was examined using linear mixed models. Age-related changes were characterized with general additive mixed models (GAMMs). Positive gene dosage effects (22qDel < TD < 22qDup) were observed for total intracranial and whole hippocampus volumes, but not whole thalamus or amygdala volumes. Several amygdala subregions exhibited similar positive effects, with bi-directional effects found across thalamic nuclei. Distinct age-related trajectories were observed across the three groups. Notably, both 22qDel and 22qDup carriers exhibited flattened development of hippocampal CA2/3 subfields relative to TD controls. This study provides novel insights into the impact of 22q11.2 CNVs on subcortical brain structures and their developmental trajectories.

Multi-Site Statistical Mapping of Along-Tract Microstructural Abnormalities in Bipolar Disorder with Diffusion MRI Tractometry

Investigating alterations in brain circuitry associated with bipolar disorder (BD) may offer a valuable approach to discover brain biomarkers for genetic and interventional studies of the disorder and related mental illnesses. Some diffusion MRI studies report evidence of microstructural abnormalities in white matter regions of interest, but we lack a fine-scale spatial mapping of brain microstructural differences along tracts in BD. We also lack large-scale studies that integrate tractometry data from multiple sites, as larger datasets can greatly enhance power to detect subtle effects and assess whether effects replicate across larger international datasets. In this multisite diffusion MRI study, we used BUndle ANalytics (BUAN, Chandio 2020), a recently developed analytic approach for tractography, to extract, map, and visualize profiles of microstructural abnormalities on 3D models of fiber tracts in 148 participants with BD and 259 healthy controls from 6 independent scan sites. Modeling site differences as random effects, we investigated along-tract white matter (WM) microstructural differences between diagnostic groups. QQ plots showed that group differences were gradually enhanced as more sites were added. Using the BUAN pipeline, BD was associated with lower mean fractional anisotropy (FA) in fronto-limbic, interhemispheric, and posterior pathways; higher FA was also noted in posterior bundles, relative to controls. By integrating tractography and anatomical information, BUAN effectively captures unique effects along white matter (WM) tracts, providing valuable insights into anatomical variations that may assist in the classification of diseases.

Subcortical Brain Alterations in Carriers of Genomic Copy Number Variants

OBJECTIVE

Copy number variants (CNVs) are well-known genetic pleiotropic risk factors for multiple neurodevelopmental and psychiatric disorders (NPDs), including autism (ASD) and schizophrenia. Little is known about how different CNVs conferring risk for the same condition may affect subcortical brain structures and how these alterations relate to the level of disease risk conferred by CNVs. To fill this gap, the authors investigated gross volume, vertex-level thickness, and surface maps of subcortical structures in 11 CNVs and six NPDs.

METHODS

Subcortical structures were characterized using harmonized ENIGMA protocols in 675 CNV carriers (CNVs at 1q21.1, TAR, 13q12.12, 15q11.2, 16p11.2, 16p13.11, and 22q11.2; age range, 6-80 years; 340 males) and 782 control subjects (age range, 6-80 years; 387 males) as well as ENIGMA summary statistics for ASD, schizophrenia, attention deficit hyperactivity disorder, obsessive-compulsive disorder, bipolar disorder, and major depression.

RESULTS

All CNVs showed alterations in at least one subcortical measure. Each structure was affected by at least two CNVs, and the hippocampus and amygdala were affected by five. Shape analyses detected subregional alterations that were averaged out in volume analyses. A common latent dimension was identified, characterized by opposing effects on the hippocampus/amygdala and putamen/pallidum, across CNVs and across NPDs. Effect sizes of CNVs on subcortical volume, thickness, and local surface area were correlated with their previously reported effect sizes on cognition and risk for ASD and schizophrenia.

CONCLUSIONS

The findings demonstrate that subcortical alterations associated with CNVs show varying levels of similarities with those associated with neuropsychiatric conditions, as well distinct effects, with some CNVs clustering with adult-onset conditions and others with ASD. These findings provide insight into the long-standing questions of why CNVs at different genomic loci increase the risk for the same NPD and why a single CNV increases the risk for a diverse set of NPDs.

Preliminary validation of a structural magnetic resonance imaging metric for tracking dementia-related neurodegeneration and future decline

Alzheimer's disease (AD) is a progressive neurodegenerative disease characterized by cognitive decline and atrophy in the medial temporal lobe (MTL) and subsequent brain regions. Structural magnetic resonance imaging (sMRI) has been widely used in research and clinical care for diagnosis and monitoring AD progression. However, atrophy patterns are complex and vary by patient. To address this issue, researchers have made efforts to develop more concise metrics that can summarize AD-specific atrophy. Many of these methods can be difficult to interpret clinically, hampering adoption. In this study, we introduce a novel index which we call an "AD-NeuroScore," that uses a modified Euclidean-inspired distance function to calculate differences between regional brain volumes associated with cognitive decline. The index is adjusted for intracranial volume (ICV), age, sex, and scanner model. We validated AD-NeuroScore using 929 older adults from the Alzheimer's Disease Neuroimaging Initiative (ADNI) study, with a mean age of 72.7 years (SD = 6.3; 55.1-91.5) and cognitively normal (CN), mild cognitive impairment (MCI), or AD diagnoses. Our validation results showed that AD-NeuroScore was significantly associated with diagnosis and disease severity scores (measured by MMSE, CDR-SB, and ADAS-11) at baseline. Furthermore, baseline AD-NeuroScore was associated with both changes in diagnosis and disease severity scores at all time points with available data. The performance of AD-NeuroScore was equivalent or superior to adjusted hippocampal volume (AHV), a widely used metric in AD research. Further, AD-NeuroScore typically performed as well as or sometimes better when compared to other existing sMRI-based metrics. In conclusion, we have introduced a new metric, AD-NeuroScore, which shows promising results in detecting AD, benchmarking disease severity, and predicting disease progression. AD-NeuroScore differentiates itself from other metrics by being clinically practical and interpretable.

Volume of subcortical brain regions in social anxiety disorder: mega-analytic results from 37 samples in the ENIGMA-Anxiety Working Group

There is limited convergence in neuroimaging investigations into volumes of subcortical brain regions in social anxiety disorder (SAD). The inconsistent findings may arise from variations in methodological approaches across studies, including sample selection based on age and clinical characteristics. The ENIGMA-Anxiety Working Group initiated a global mega-analysis to determine whether differences in subcortical volumes can be detected in adults and adolescents with SAD relative to healthy controls. Volumetric data from 37 international samples with 1115 SAD patients and 2775 controls were obtained from ENIGMA-standardized protocols for image segmentation and quality assurance. Linear mixed-effects analyses were adjusted for comparisons across seven subcortical regions in each hemisphere using family-wise error (FWE)-correction. Mixed-effects d effect sizes were calculated. In the full sample, SAD patients showed smaller bilateral putamen volume than controls (left: d = -0.077, pFWE = 0.037; right: d = -0.104, pFWE = 0.001), and a significant interaction between SAD and age was found for the left putamen (r = -0.034, pFWE = 0.045). Smaller bilateral putamen volumes (left: d = -0.141, pFWE < 0.001; right: d = -0.158, pFWE < 0.001) and larger bilateral pallidum volumes (left: d = 0.129, pFWE = 0.006; right: d = 0.099, pFWE = 0.046) were detected in adult SAD patients relative to controls, but no volumetric differences were apparent in adolescent SAD patients relative to controls. Comorbid anxiety disorders and age of SAD onset were additional determinants of SAD-related volumetric differences in subcortical regions. To conclude, subtle volumetric alterations in subcortical regions in SAD were detected. Heterogeneity in age and clinical characteristics may partly explain inconsistencies in previous findings. The association between alterations in subcortical volumes and SAD illness progression deserves further investigation, especially from adolescence into adulthood.

Mega-analysis of association between obesity and cortical morphology in bipolar disorders: ENIGMA study in 2832 participants

BACKGROUND

Obesity is highly prevalent and disabling, especially in individuals with severe mental illness including bipolar disorders (BD). The brain is a target organ for both obesity and BD. Yet, we do not understand how cortical brain alterations in BD and obesity interact.

METHODS

We obtained body mass index (BMI) and MRI-derived regional cortical thickness, surface area from 1231 BD and 1601 control individuals from 13 countries within the ENIGMA-BD Working Group. We jointly modeled the statistical effects of BD and BMI on brain structure using mixed effects and tested for interaction and mediation. We also investigated the impact of medications on the BMI-related associations.

RESULTS

BMI and BD additively impacted the structure of many of the same brain regions. Both BMI and BD were negatively associated with cortical thickness, but not surface area. In most regions the number of jointly used psychiatric medication classes remained associated with lower cortical thickness when controlling for BMI. In a single region, fusiform gyrus, about a third of the negative association between number of jointly used psychiatric medications and cortical thickness was mediated by association between the number of medications and higher BMI.

CONCLUSIONS

We confirmed consistent associations between higher BMI and lower cortical thickness, but not surface area, across the cerebral mantle, in regions which were also associated with BD. Higher BMI in people with BD indicated more pronounced brain alterations. BMI is important for understanding the neuroanatomical changes in BD and the effects of psychiatric medications on the brain.

Subcortical brain alterations in carriers of genomic copy number variants

Objectives

Copy number variants (CNVs) are well-known genetic pleiotropic risk factors for multiple neurodevelopmental and psychiatric disorders (NPDs) including autism (ASD) and schizophrenia (SZ). Overall, little is known about how different CNVs conferring risk for the same condition may affect subcortical brain structures and how these alterations relate to the level of disease risk conferred by CNVs. To fill this gap, we investigated gross volume, and vertex level thickness and surface maps of subcortical structures in 11 different CNVs and 6 different NPDs.

Methods

Subcortical structures were characterized using harmonized ENIGMA protocols in 675 CNV carriers (at the following loci: 1q21.1, TAR, 13q12.12, 15q11.2, 16p11.2, 16p13.11, and 22q11.2) and 782 controls (Male/Female: 727/730; age-range: 6-80 years) as well as ENIGMA summary-statistics for ASD, SZ, ADHD, Obsessive-Compulsive-Disorder, Bipolar-Disorder, and Major-Depression.

Results

Nine of the 11 CNVs affected volume of at least one subcortical structure. The hippocampus and amygdala were affected by five CNVs. Effect sizes of CNVs on subcortical volume, thickness and local surface area were correlated with their previously reported effect sizes on cognition and risk for ASD and SZ. Shape analyses were able to identify subregional alterations that were averaged out in volume analyses. We identified a common latent dimension - characterized by opposing effects on basal ganglia and limbic structures - across CNVs and across NPDs.

Conclusion

Our findings demonstrate that subcortical alterations associated with CNVs show varying levels of similarities with those associated with neuropsychiatric conditions. We also observed distinct effects with some CNVs clustering with adult conditions while others clustered with ASD. This large cross-CNV and NPDs analysis provide insight into the long-standing questions of why CNVs at different genomic loci increase the risk for the same NPD, as well as why a single CNV increases the risk for a diverse set of NPDs.

Harmonization of multi-scanner in vivo magnetic resonance spectroscopy: ENIGMA consortium task group considerations

Magnetic resonance spectroscopy is a powerful, non-invasive, quantitative imaging technique that allows for the measurement of brain metabolites that has demonstrated utility in diagnosing and characterizing a broad range of neurological diseases. Its impact, however, has been limited due to small sample sizes and methodological variability in addition to intrinsic limitations of the method itself such as its sensitivity to motion. The lack of standardization from a data acquisition and data processing perspective makes it difficult to pool multiple studies and/or conduct multisite studies that are necessary for supporting clinically relevant findings. Based on the experience of the ENIGMA MRS work group and a review of the literature, this manuscript provides an overview of the current state of MRS data harmonization. Key factors that need to be taken into consideration when conducting both retrospective and prospective studies are described. These include (1) MRS acquisition issues such as pulse sequence, RF and B0 calibrations, echo time, and SNR; (2) data processing issues such as pre-processing steps, modeling, and quantitation; and (3) biological factors such as voxel location, age, sex, and pathology. Various approaches to MRS data harmonization are then described including meta-analysis, mega-analysis, linear modeling, ComBat and artificial intelligence approaches. The goal is to provide both novice and experienced readers with the necessary knowledge for conducting MRS data harmonization studies.

A comparison of methods to harmonize cortical thickness measurements across scanners and sites

Results of neuroimaging datasets aggregated from multiple sites may be biased by site-specific profiles in participants' demographic and clinical characteristics, as well as MRI acquisition protocols and scanning platforms. We compared the impact of four different harmonization methods on results obtained from analyses of cortical thickness data: (1) linear mixed-effects model (LME) that models site-specific random intercepts (LMEINT), (2) LME that models both site-specific random intercepts and age-related random slopes (LMEINT+SLP), (3) ComBat, and (4) ComBat with a generalized additive model (ComBat-GAM). Our test case for comparing harmonization methods was cortical thickness data aggregated from 29 sites, which included 1,340 cases with posttraumatic stress disorder (PTSD) (6.2-81.8 years old) and 2,057 trauma-exposed controls without PTSD (6.3-85.2 years old). We found that, compared to the other data harmonization methods, data processed with ComBat-GAM was more sensitive to the detection of significant case-control differences (Χ2(3) = 63.704, p < 0.001) as well as case-control differences in age-related cortical thinning (Χ2(3) = 12.082, p = 0.007). Both ComBat and ComBat-GAM outperformed LME methods in detecting sex differences (Χ2(3) = 9.114, p = 0.028) in regional cortical thickness. ComBat-GAM also led to stronger estimates of age-related declines in cortical thickness (corrected p-values < 0.001), stronger estimates of case-related cortical thickness reduction (corrected p-values < 0.001), weaker estimates of age-related declines in cortical thickness in cases than controls (corrected p-values < 0.001), stronger estimates of cortical thickness reduction in females than males (corrected p-values < 0.001), and stronger estimates of cortical thickness reduction in females relative to males in cases than controls (corrected p-values < 0.001). Our results support the use of ComBat-GAM to minimize confounds and increase statistical power when harmonizing data with non-linear effects, and the use of either ComBat or ComBat-GAM for harmonizing data with linear effects.

Site effects how-to and when: An overview of retrospective techniques to accommodate site effects in multi-site neuroimaging analyses

Site differences, or systematic differences in feature distributions across multiple data-acquisition sites, are a known source of heterogeneity that may adversely affect large-scale meta- and mega-analyses of independently collected neuroimaging data. They influence nearly all multi-site imaging modalities and biomarkers, and methods to compensate for them can improve reliability and generalizability in the analysis of genetics, omics, and clinical data. The origins of statistical site effects are complex and involve both technical differences (scanner vendor, head coil, acquisition parameters, imaging processing) and differences in sample characteristics (inclusion/exclusion criteria, sample size, ancestry) between sites. In an age of expanding international consortium research, there is a growing need to disentangle technical site effects from sample characteristics of interest. Numerous statistical and machine learning methods have been developed to control for, model, or attenuate site effects - yet to date, no comprehensive review has discussed the benefits and drawbacks of each for different use cases. Here, we provide an overview of the different existing statistical and machine learning methods developed to remove unwanted site effects from independently collected neuroimaging samples. We focus on linear mixed effect models, the ComBat technique and its variants, adjustments based on image quality metrics, normative modeling, and deep learning approaches such as generative adversarial networks. For each method, we outline the statistical foundation and summarize strengths and weaknesses, including their assumptions and conditions of use. We provide information on software availability and comment on the ease of use and the applicability of these methods to different types of data. We discuss validation and comparative reports, mention caveats and provide guidance on when to use each method, depending on context and specific research questions.

RORA rs8042149 polymorphism moderates the association between PTSD symptom severity and transverse temporal gyrus thickness in Han Chinese adults who lost their only child

BACKGROUND

The G allele in retinoid-related orphan receptor alpha (RORA, rs8042149) gene is associated with post-traumatic stress disorder (PTSD) diagnosis and more severe symptoms, reported in the first genome-wide association study of PTSD and subsequent replication studies. Although recent MRI studies identified brain structural deficits in RORA rs8042149 risk G allele carriers, the neural mechanism underlying RORA-related brain structural changes in PTSD remains poorly understood.

METHODS

This study included 227 Han Chinese adults who lost their only child. Cortical thickness and subcortical volume were extracted using FreeSurfer, and PTSD severity was assessed using the Clinician-Administered PTSD Scale. Hierarchical linear regression was used to assess the interaction effect between RORA genotypes (T/T, G/T, and G/G) and PTSD severity on cortical and subcortical structures.

RESULTS

Significant genotype × PTSD symptom severity interaction effects were found for bilateral transverse temporal gyrus thickness. For individuals with the homozygous T/T genotype, current PTSD symptom severity was positively associated with bilateral transverse temporal gyrus thickness. For individuals with heterozygous G/T genotype, current PTSD symptom severity was negatively associated with the left transverse temporal gyrus thickness. No significant main or interaction effects were found in any subcortical regions.

LIMITATION

Cross-sectional design of this study.

CONCLUSION

These findings suggest that the non-risk T/T genotype - but not the risk G allele carriers - has a potentially protective or compensatory role on temporal gyrus thickness in adults who lost their only child. These results highlight the moderation effect of RORA polymorphism on the relationship between PTSD symptom severity and cortical structural changes.

Shape analysis of subcortical structures in obsessive-compulsive disorder and the relationship with comorbid anxiety, depression, and medication use: A meta-analysis by the OCD Brain Imaging Consortium

OBJECTIVE

Neuroimaging studies of obsessive-compulsive disorder (OCD) patients have highlighted the important role of deep gray matter structures. Less work has however focused on subcortical shape in OCD patients.

METHODS

Here we pooled brain MRI scans from 412 OCD patients and 368 controls to perform a meta-analysis utilizing the ENIGMA-Shape protocol. In addition, we investigated modulating effects of medication status, comorbid anxiety or depression, and disease duration on subcortical shape.

RESULTS

There was no significant difference in shape thickness or surface area between OCD patients and healthy controls. For the subgroup analyses, OCD patients with comorbid depression or anxiety had lower thickness of the hippocampus and caudate nucleus and higher thickness of the putamen and pallidum compared to controls. OCD patients with comorbid depression had lower shape surface area in the thalamus, caudate nucleus, putamen, hippocampus, and nucleus accumbens and higher shape surface area in the pallidum. OCD patients with comorbid anxiety had lower shape surface area in the putamen and the left caudate nucleus and higher shape surface area in the pallidum and the right caudate nucleus. Further, OCD patients on medication had lower shape thickness of the putamen, thalamus, and hippocampus and higher thickness of the pallidum and caudate nucleus, as well as lower shape surface area in the hippocampus and amygdala and higher surface area in the putamen, pallidum, and caudate nucleus compared to controls. There were no significant differences between OCD patients without co-morbid anxiety and/or depression and healthy controls on shape measures. In addition, there were also no significant differences between OCD patients not using medication and healthy controls.

CONCLUSIONS

The findings here are partly consistent with prior work on brain volumes in OCD, insofar as they emphasize that alterations in subcortical brain morphology are associated with comorbidity and medication status. Further work is needed to understand the biological processes contributing to subcortical shape.

Longitudinal trajectories of cortical development in 22q11.2 copy number variants and typically developing controls

Probing naturally-occurring, reciprocal genomic copy number variations (CNVs) may help us understand mechanisms that underlie deviations from typical brain development. Cross-sectional studies have identified prominent reductions in cortical surface area (SA) and increased cortical thickness (CT) in 22q11.2 deletion carriers (22qDel), with the opposite pattern in duplication carriers (22qDup), but the longitudinal trajectories of these anomalies-and their relationship to clinical symptomatology-are unknown. Here, we examined neuroanatomic changes within a longitudinal cohort of 261 22q11.2 CNV carriers and demographically-matched typically developing (TD) controls (84 22qDel, 34 22qDup, and 143 TD; mean age 18.35, ±10.67 years; 50.47% female). A total of 431 magnetic resonance imaging scans (164 22qDel, 59 22qDup, and 208 TD control scans; mean interscan interval = 20.27 months) were examined. Longitudinal FreeSurfer analysis pipelines were used to parcellate the cortex and calculate average CT and SA for each region. First, general additive mixed models (GAMMs) were used to identify regions with between-group differences in developmental trajectories. Secondly, we investigated whether these trajectories were associated with clinical outcomes. Developmental trajectories of CT were more protracted in 22qDel relative to TD and 22qDup. 22qDup failed to show normative age-related SA decreases. 22qDel individuals with psychosis spectrum symptoms showed two distinct periods of altered CT trajectories relative to 22qDel without psychotic symptoms. In contrast, 22q11.2 CNV carriers with autism spectrum diagnoses showed early alterations in SA trajectories. Collectively, these results provide new insights into altered neurodevelopment in 22q11.2 CNV carriers, which may shed light on neural mechanisms underlying distinct clinical outcomes.

Remodeling of the Cortical Structural Connectome in Posttraumatic Stress Disorder: Results From the ENIGMA-PGC Posttraumatic Stress Disorder Consortium

BACKGROUND

Posttraumatic stress disorder (PTSD) is accompanied by disrupted cortical neuroanatomy. We investigated alteration in covariance of structural networks associated with PTSD in regions that demonstrate the case-control differences in cortical thickness (CT) and surface area (SA).

METHODS

Neuroimaging and clinical data were aggregated from 29 research sites in >1300 PTSD cases and >2000 trauma-exposed control subjects (ages 6.2-85.2 years) by the ENIGMA-PGC (Enhancing Neuro Imaging Genetics through Meta Analysis-Psychiatric Genomics Consortium) PTSD working group. Cortical regions in the network were rank ordered by the effect size of PTSD-related cortical differences in CT and SA. The top-n (n = 2-148) regions with the largest effect size for PTSD > non-PTSD formed hypertrophic networks, the largest effect size for PTSD < non-PTSD formed atrophic networks, and the smallest effect size of between-group differences formed stable networks. The mean structural covariance (SC) of a given n-region network was the average of all positive pairwise correlations and was compared with the mean SC of 5000 randomly generated n-region networks.

RESULTS

Patients with PTSD, relative to non-PTSD control subjects, exhibited lower mean SC in CT-based and SA-based atrophic networks. Comorbid depression, sex, and age modulated covariance differences of PTSD-related structural networks.

CONCLUSIONS

Covariance of structural networks based on CT and cortical SA are affected by PTSD and further modulated by comorbid depression, sex, and age. The SC networks that are perturbed in PTSD comport with converging evidence from resting-state functional connectivity networks and networks affected by inflammatory processes and stress hormones in PTSD.

Virtual Ontogeny of Cortical Growth Preceding Mental Illness

BACKGROUND

Morphology of the human cerebral cortex differs across psychiatric disorders, with neurobiology and developmental origins mostly undetermined. Deviations in the tangential growth of the cerebral cortex during pre/perinatal periods may be reflected in individual variations in cortical surface area later in life.

METHODS

Interregional profiles of group differences in surface area between cases and controls were generated using T1-weighted magnetic resonance imaging from 27,359 individuals including those with attention-deficit/hyperactivity disorder, autism spectrum disorder, bipolar disorder, major depressive disorder, schizophrenia, and high general psychopathology (through the Child Behavior Checklist). Similarity of interregional profiles of group differences in surface area and prenatal cell-specific gene expression was assessed.

RESULTS

Across the 11 cortical regions, group differences in cortical area for attention-deficit/hyperactivity disorder, schizophrenia, and Child Behavior Checklist were dominant in multimodal association cortices. The same interregional profiles were also associated with interregional profiles of (prenatal) gene expression specific to proliferative cells, namely radial glia and intermediate progenitor cells (greater expression, larger difference), as well as differentiated cells, namely excitatory neurons and endothelial and mural cells (greater expression, smaller difference). Finally, these cell types were implicated in known pre/perinatal risk factors for psychosis. Genes coexpressed with radial glia were enriched with genes implicated in congenital abnormalities, birth weight, hypoxia, and starvation. Genes coexpressed with endothelial and mural genes were enriched with genes associated with maternal hypertension and preterm birth.

CONCLUSIONS

Our findings support a neurodevelopmental model of vulnerability to mental illness whereby prenatal risk factors acting through cell-specific processes lead to deviations from typical brain development during pregnancy.

Local molecular and global connectomic contributions to cross-disorder cortical abnormalities

Numerous brain disorders demonstrate structural brain abnormalities, which are thought to arise from molecular perturbations or connectome miswiring. The unique and shared contributions of these molecular and connectomic vulnerabilities to brain disorders remain unknown, and has yet to be studied in a single multi-disorder framework. Using MRI morphometry from the ENIGMA consortium, we construct maps of cortical abnormalities for thirteen neurodevelopmental, neurological, and psychiatric disorders from N = 21,000 participants and N = 26,000 controls, collected using a harmonised processing protocol. We systematically compare cortical maps to multiple micro-architectural measures, including gene expression, neurotransmitter density, metabolism, and myelination (molecular vulnerability), as well as global connectomic measures including number of connections, centrality, and connection diversity (connectomic vulnerability). We find a relationship between molecular vulnerability and white-matter architecture that drives cortical disorder profiles. Local attributes, particularly neurotransmitter receptor profiles, constitute the best predictors of both disorder-specific cortical morphology and cross-disorder similarity. Finally, we find that cross-disorder abnormalities are consistently subtended by a small subset of network epicentres in bilateral sensory-motor, inferior temporal lobe, precuneus, and superior parietal cortex. Collectively, our results highlight how local molecular attributes and global connectivity jointly shape cross-disorder cortical abnormalities.

Diagnosis of bipolar disorders and body mass index predict clustering based on similarities in cortical thickness-ENIGMA study in 2436 individuals

AIMS

Rates of obesity have reached epidemic proportions, especially among people with psychiatric disorders. While the effects of obesity on the brain are of major interest in medicine, they remain markedly under-researched in psychiatry.

METHODS

We obtained body mass index (BMI) and magnetic resonance imaging-derived regional cortical thickness, surface area from 836 bipolar disorders (BD) and 1600 control individuals from 14 sites within the ENIGMA-BD Working Group. We identified regionally specific profiles of cortical thickness using K-means clustering and studied clinical characteristics associated with individual cortical profiles.

RESULTS

We detected two clusters based on similarities among participants in cortical thickness. The lower thickness cluster (46.8% of the sample) showed thinner cortex, especially in the frontal and temporal lobes and was associated with diagnosis of BD, higher BMI, and older age. BD individuals in the low thickness cluster were more likely to have the diagnosis of bipolar disorder I and less likely to be treated with lithium. In contrast, clustering based on similarities in the cortical surface area was unrelated to BD or BMI and only tracked age and sex.

CONCLUSIONS

We provide evidence that both BD and obesity are associated with similar alterations in cortical thickness, but not surface area. The fact that obesity increased the chance of having low cortical thickness could explain differences in cortical measures among people with BD. The thinner cortex in individuals with higher BMI, which was additive and similar to the BD-associated alterations, may suggest that treating obesity could lower the extent of cortical thinning in BD.

Cross disorder comparisons of brain structure in schizophrenia, bipolar disorder, major depressive disorder, and 22q11.2 deletion syndrome: A review of ENIGMA findings

This review compares the main brain abnormalities in schizophrenia (SZ), bipolar disorder (BD), major depressive disorder (MDD), and 22q11.2 Deletion Syndrome (22q11DS) determined by ENIGMA (Enhancing Neuro Imaging Genetics through Meta Analysis) consortium investigations. We obtained ranked effect sizes for subcortical volumes, regional cortical thickness, cortical surface area, and diffusion tensor imaging abnormalities, comparing each of these disorders relative to healthy controls. In addition, the studies report on significant associations between brain imaging metrics and disorder-related factors such as symptom severity and treatments. Visual comparison of effect size profiles shows that effect sizes are generally in the same direction and scale in severity with the disorders (in the order SZ > BD > MDD). The effect sizes for 22q11DS, a rare genetic syndrome that increases the risk for psychiatric disorders, appear to be much larger than for either of the complex psychiatric disorders. This is consistent with the idea of generally larger effects on the brain of rare compared to common genetic variants. Cortical thickness and surface area effect sizes for 22q11DS with psychosis compared to 22q11DS without psychosis are more similar to those of SZ and BD than those of MDD; a pattern not observed for subcortical brain structures and fractional anisotropy effect sizes. The observed similarities in effect size profiles for cortical measures across the psychiatric disorders mimic those observed for shared genetic variance between these disorders reported based on family and genetic studies and are consistent with shared genetic risk for SZ and BD and structural brain phenotypes.

Longitudinal Structural Brain Changes in Bipolar Disorder: A Multicenter Neuroimaging Study of 1232 Individuals by the ENIGMA Bipolar Disorder Working Group

BACKGROUND

Bipolar disorder (BD) is associated with cortical and subcortical structural brain abnormalities. It is unclear whether such alterations progressively change over time, and how this is related to the number of mood episodes. To address this question, we analyzed a large and diverse international sample with longitudinal magnetic resonance imaging (MRI) and clinical data to examine structural brain changes over time in BD.

METHODS

Longitudinal structural MRI and clinical data from the ENIGMA (Enhancing Neuro Imaging Genetics through Meta Analysis) BD Working Group, including 307 patients with BD and 925 healthy control subjects, were collected from 14 sites worldwide. Male and female participants, aged 40 ± 17 years, underwent MRI at 2 time points. Cortical thickness, surface area, and subcortical volumes were estimated using FreeSurfer. Annualized change rates for each imaging phenotype were compared between patients with BD and healthy control subjects. Within patients, we related brain change rates to the number of mood episodes between time points and tested for effects of demographic and clinical variables.

RESULTS

Compared with healthy control subjects, patients with BD showed faster enlargement of ventricular volumes and slower thinning of the fusiform and parahippocampal cortex (0.18 <d < 0.22). More (hypo)manic episodes were associated with faster cortical thinning, primarily in the prefrontal cortex.

CONCLUSIONS

In the hitherto largest longitudinal MRI study on BD, we did not detect accelerated cortical thinning but noted faster ventricular enlargements in BD. However, abnormal frontocortical thinning was observed in association with frequent manic episodes. Our study yields insights into disease progression in BD and highlights the importance of mania prevention in BD treatment.

Genomic and neuroimaging approaches to bipolar disorder

BACKGROUND

To date, besides genome-wide association studies, a variety of other genetic analyses (e.g. polygenic risk scores, whole-exome sequencing and whole-genome sequencing) have been conducted, and a large amount of data has been gathered for investigating the involvement of common, rare and very rare types of DNA sequence variants in bipolar disorder. Also, non-invasive neuroimaging methods can be used to quantify changes in brain structure and function in patients with bipolar disorder.

AIMS

To provide a comprehensive assessment of genetic findings associated with bipolar disorder, based on the evaluation of different genomic approaches and neuroimaging studies.

METHOD

We conducted a PubMed search of all relevant literatures from the beginning to the present, by querying related search strings.

RESULTS

ANK3, CACNA1C, SYNE1, ODZ4 and TRANK1 are five genes that have been replicated as key gene candidates in bipolar disorder pathophysiology, through the investigated studies. The percentage of phenotypic variance explained by the identified variants is small (approximately 4.7%). Bipolar disorder polygenic risk scores are associated with other psychiatric phenotypes. The ENIGMA-BD studies show a replicable pattern of lower cortical thickness, altered white matter integrity and smaller subcortical volumes in bipolar disorder.

CONCLUSIONS

The low amount of explained phenotypic variance highlights the need for further large-scale investigations, especially among non-European populations, to achieve a more complete understanding of the genetic architecture of bipolar disorder and the missing heritability. Combining neuroimaging data with genetic data in large-scale studies might help researchers acquire a better knowledge of the engaged brain regions in bipolar disorder.

What we learn about bipolar disorder from large-scale neuroimaging: Findings and future directions from the ENIGMA Bipolar Disorder Working Group

MRI-derived brain measures offer a link between genes, the environment and behavior and have been widely studied in bipolar disorder (BD). However, many neuroimaging studies of BD have been underpowered, leading to varied results and uncertainty regarding effects. The Enhancing Neuro Imaging Genetics through Meta-Analysis (ENIGMA) Bipolar Disorder Working Group was formed in 2012 to empower discoveries, generate consensus findings and inform future hypothesis-driven studies of BD. Through this effort, over 150 researchers from 20 countries and 55 institutions pool data and resources to produce the largest neuroimaging studies of BD ever conducted. The ENIGMA Bipolar Disorder Working Group applies standardized processing and analysis techniques to empower large-scale meta- and mega-analyses of multimodal brain MRI and improve the replicability of studies relating brain variation to clinical and genetic data. Initial BD Working Group studies reveal widespread patterns of lower cortical thickness, subcortical volume and disrupted white matter integrity associated with BD. Findings also include mapping brain alterations of common medications like lithium, symptom patterns and clinical risk profiles and have provided further insights into the pathophysiological mechanisms of BD. Here we discuss key findings from the BD working group, its ongoing projects and future directions for large-scale, collaborative studies of mental illness.

Effects of copy number variations on brain structure and risk for psychiatric illness: Large-scale studies from the ENIGMA working groups on CNVs

The Enhancing NeuroImaging Genetics through Meta-Analysis copy number variant (ENIGMA-CNV) and 22q11.2 Deletion Syndrome Working Groups (22q-ENIGMA WGs) were created to gain insight into the involvement of genetic factors in human brain development and related cognitive, psychiatric and behavioral manifestations. To that end, the ENIGMA-CNV WG has collated CNV and magnetic resonance imaging (MRI) data from ~49,000 individuals across 38 global research sites, yielding one of the largest studies to date on the effects of CNVs on brain structures in the general population. The 22q-ENIGMA WG includes 12 international research centers that assessed over 533 individuals with a confirmed 22q11.2 deletion syndrome, 40 with 22q11.2 duplications, and 333 typically developing controls, creating the largest-ever 22q11.2 CNV neuroimaging data set. In this review, we outline the ENIGMA infrastructure and procedures for multi-site analysis of CNVs and MRI data. So far, ENIGMA has identified effects of the 22q11.2, 16p11.2 distal, 15q11.2, and 1q21.1 distal CNVs on subcortical and cortical brain structures. Each CNV is associated with differences in cognitive, neurodevelopmental and neuropsychiatric traits, with characteristic patterns of brain structural abnormalities. Evidence of gene-dosage effects on distinct brain regions also emerged, providing further insight into genotype-phenotype relationships. Taken together, these results offer a more comprehensive picture of molecular mechanisms involved in typical and atypical brain development. This "genotype-first" approach also contributes to our understanding of the etiopathogenesis of brain disorders. Finally, we outline future directions to better understand effects of CNVs on brain structure and behavior.

Subcortical volumes across the lifespan: Data from 18,605 healthy individuals aged 3-90 years

Age has a major effect on brain volume. However, the normative studies available are constrained by small sample sizes, restricted age coverage and significant methodological variability. These limitations introduce inconsistencies and may obscure or distort the lifespan trajectories of brain morphometry. In response, we capitalized on the resources of the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) Consortium to examine age-related trajectories inferred from cross-sectional measures of the ventricles, the basal ganglia (caudate, putamen, pallidum, and nucleus accumbens), the thalamus, hippocampus and amygdala using magnetic resonance imaging data obtained from 18,605 individuals aged 3-90 years. All subcortical structure volumes were at their maximum value early in life. The volume of the basal ganglia showed a monotonic negative association with age thereafter; there was no significant association between age and the volumes of the thalamus, amygdala and the hippocampus (with some degree of decline in thalamus) until the sixth decade of life after which they also showed a steep negative association with age. The lateral ventricles showed continuous enlargement throughout the lifespan. Age was positively associated with inter-individual variability in the hippocampus and amygdala and the lateral ventricles. These results were robust to potential confounders and could be used to examine the functional significance of deviations from typical age-related morphometric patterns.

ENIGMA-DTI: Translating reproducible white matter deficits into personalized vulnerability metrics in cross-diagnostic psychiatric research

The ENIGMA-DTI (diffusion tensor imaging) workgroup supports analyses that examine the effects of psychiatric, neurological, and developmental disorders on the white matter pathways of the human brain, as well as the effects of normal variation and its genetic associations. The seven ENIGMA disorder-oriented working groups used the ENIGMA-DTI workflow to derive patterns of deficits using coherent and coordinated analyses that model the disease effects across cohorts worldwide. This yielded the largest studies detailing patterns of white matter deficits in schizophrenia spectrum disorder (SSD), bipolar disorder (BD), major depressive disorder (MDD), obsessive-compulsive disorder (OCD), posttraumatic stress disorder (PTSD), traumatic brain injury (TBI), and 22q11 deletion syndrome. These deficit patterns are informative of the underlying neurobiology and reproducible in independent cohorts. We reviewed these findings, demonstrated their reproducibility in independent cohorts, and compared the deficit patterns across illnesses. We discussed translating ENIGMA-defined deficit patterns on the level of individual subjects using a metric called the regional vulnerability index (RVI), a correlation of an individual's brain metrics with the expected pattern for a disorder. We discussed the similarity in white matter deficit patterns among SSD, BD, MDD, and OCD and provided a rationale for using this index in cross-diagnostic neuropsychiatric research. We also discussed the difference in deficit patterns between idiopathic schizophrenia and 22q11 deletion syndrome, which is used as a developmental and genetic model of schizophrenia. Together, these findings highlight the importance of collaborative large-scale research to provide robust and reproducible effects that offer insights into individual vulnerability and cross-diagnosis features.

A meta-analysis of deep brain structural shape and asymmetry abnormalities in 2,833 individuals with schizophrenia compared with 3,929 healthy volunteers via the ENIGMA Consortium

Schizophrenia is associated with widespread alterations in subcortical brain structure. While analytic methods have enabled more detailed morphometric characterization, findings are often equivocal. In this meta-analysis, we employed the harmonized ENIGMA shape analysis protocols to collaboratively investigate subcortical brain structure shape differences between individuals with schizophrenia and healthy control participants. The study analyzed data from 2,833 individuals with schizophrenia and 3,929 healthy control participants contributed by 21 worldwide research groups participating in the ENIGMA Schizophrenia Working Group. Harmonized shape analysis protocols were applied to each site's data independently for bilateral hippocampus, amygdala, caudate, accumbens, putamen, pallidum, and thalamus obtained from T1-weighted structural MRI scans. Mass univariate meta-analyses revealed more-concave-than-convex shape differences in the hippocampus, amygdala, accumbens, and thalamus in individuals with schizophrenia compared with control participants, more-convex-than-concave shape differences in the putamen and pallidum, and both concave and convex shape differences in the caudate. Patterns of exaggerated asymmetry were observed across the hippocampus, amygdala, and thalamus in individuals with schizophrenia compared to control participants, while diminished asymmetry encompassed ventral striatum and ventral and dorsal thalamus. Our analyses also revealed that higher chlorpromazine dose equivalents and increased positive symptom levels were associated with patterns of contiguous convex shape differences across multiple subcortical structures. Findings from our shape meta-analysis suggest that common neurobiological mechanisms may contribute to gray matter reduction across multiple subcortical regions, thus enhancing our understanding of the nature of network disorganization in schizophrenia.

ENIGMA + COINSTAC: Improving Findability, Accessibility, Interoperability, and Re-usability

The FAIR principles, as applied to clinical and neuroimaging data, reflect the goal of making research products Findable, Accessible, Interoperable, and Reusable. The use of the Collaborative Informatics and Neuroimaging Suite Toolkit for Anonymized Computation (COINSTAC) platform in the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) consortium combines the technological approach of decentralized analyses with the sociological approach of sharing data. In addition, ENIGMA + COINSTAC provides a platform to facilitate the use of machine-actionable data objects. We first present how ENIGMA and COINSTAC support the FAIR principles, and then showcase their integration with a decentralized meta-analysis of sex differences in negative symptom severity in schizophrenia, and finally present ongoing activities and plans to advance FAIR principles in ENIGMA + COINSTAC. ENIGMA and COINSTAC currently represent efforts toward improved Access, Interoperability, and Reusability. We highlight additional improvements needed in these areas, as well as future connections to other resources for expanded Findability.

Cortical thickness across the lifespan: Data from 17,075 healthy individuals aged 3-90 years

Delineating the association of age and cortical thickness in healthy individuals is critical given the association of cortical thickness with cognition and behavior. Previous research has shown that robust estimates of the association between age and brain morphometry require large-scale studies. In response, we used cross-sectional data from 17,075 individuals aged 3-90 years from the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) Consortium to infer age-related changes in cortical thickness. We used fractional polynomial (FP) regression to quantify the association between age and cortical thickness, and we computed normalized growth centiles using the parametric Lambda, Mu, and Sigma method. Interindividual variability was estimated using meta-analysis and one-way analysis of variance. For most regions, their highest cortical thickness value was observed in childhood. Age and cortical thickness showed a negative association; the slope was steeper up to the third decade of life and more gradual thereafter; notable exceptions to this general pattern were entorhinal, temporopolar, and anterior cingulate cortices. Interindividual variability was largest in temporal and frontal regions across the lifespan. Age and its FP combinations explained up to 59% variance in cortical thickness. These results may form the basis of further investigation on normative deviation in cortical thickness and its significance for behavioral and cognitive outcomes.

Association between body mass index and subcortical brain volumes in bipolar disorders-ENIGMA study in 2735 individuals

Individuals with bipolar disorders (BD) frequently suffer from obesity, which is often associated with neurostructural alterations. Yet, the effects of obesity on brain structure in BD are under-researched. We obtained MRI-derived brain subcortical volumes and body mass index (BMI) from 1134 BD and 1601 control individuals from 17 independent research sites within the ENIGMA-BD Working Group. We jointly modeled the effects of BD and BMI on subcortical volumes using mixed-effects modeling and tested for mediation of group differences by obesity using nonparametric bootstrapping. All models controlled for age, sex, hemisphere, total intracranial volume, and data collection site. Relative to controls, individuals with BD had significantly higher BMI, larger lateral ventricular volume, and smaller volumes of amygdala, hippocampus, pallidum, caudate, and thalamus. BMI was positively associated with ventricular and amygdala and negatively with pallidal volumes. When analyzed jointly, both BD and BMI remained associated with volumes of lateral ventricles  and amygdala. Adjusting for BMI decreased the BD vs control differences in ventricular volume. Specifically, 18.41% of the association between BD and ventricular volume was mediated by BMI (Z = 2.73, p = 0.006). BMI was associated with similar regional brain volumes as BD, including lateral ventricles, amygdala, and pallidum. Higher BMI may in part account for larger ventricles, one of the most replicated findings in BD. Comorbidity with obesity could explain why neurostructural alterations are more pronounced in some individuals with BD. Future prospective brain imaging studies should investigate whether obesity could be a modifiable risk factor for neuroprogression.

Brain aging in major depressive disorder: results from the ENIGMA major depressive disorder working group

Major depressive disorder (MDD) is associated with an increased risk of brain atrophy, aging-related diseases, and mortality. We examined potential advanced brain aging in adult MDD patients, and whether this process is associated with clinical characteristics in a large multicenter international dataset. We performed a mega-analysis by pooling brain measures derived from T1-weighted MRI scans from 19 samples worldwide. Healthy brain aging was estimated by predicting chronological age (18-75 years) from 7 subcortical volumes, 34 cortical thickness and 34 surface area, lateral ventricles and total intracranial volume measures separately in 952 male and 1236 female controls from the ENIGMA MDD working group. The learned model coefficients were applied to 927 male controls and 986 depressed males, and 1199 female controls and 1689 depressed females to obtain independent unbiased brain-based age predictions. The difference between predicted "brain age" and chronological age was calculated to indicate brain-predicted age difference (brain-PAD). On average, MDD patients showed a higher brain-PAD of +1.08 (SE 0.22) years (Cohen's d = 0.14, 95% CI: 0.08-0.20) compared with controls. However, this difference did not seem to be driven by specific clinical characteristics (recurrent status, remission status, antidepressant medication use, age of onset, or symptom severity). This highly powered collaborative effort showed subtle patterns of age-related structural brain abnormalities in MDD. Substantial within-group variance and overlap between groups were observed. Longitudinal studies of MDD and somatic health outcomes are needed to further assess the clinical value of these brain-PAD estimates.

Prioritizing Genetic Contributors to Cortical Alterations in 22q11.2 Deletion Syndrome Using Imaging Transcriptomics

22q11.2 deletion syndrome (22q11DS) results from a hemizygous deletion that typically spans 46 protein-coding genes and is associated with widespread alterations in brain morphology. The specific genetic mechanisms underlying these alterations remain unclear. In the 22q11.2 ENIGMA Working Group, we characterized cortical alterations in individuals with 22q11DS (n = 232) versus healthy individuals (n = 290) and conducted spatial convergence analyses using gene expression data from the Allen Human Brain Atlas to prioritize individual genes that may contribute to altered surface area (SA) and cortical thickness (CT) in 22q11DS. Total SA was reduced in 22q11DS (Z-score deviance = -1.04), with prominent reductions in midline posterior and lateral association regions. Mean CT was thicker in 22q11DS (Z-score deviance = +0.64), with focal thinning in a subset of regions. Regional expression of DGCR8 was robustly associated with regional severity of SA deviance in 22q11DS; AIFM3 was also associated with SA deviance. Conversely, P2RX6 was associated with CT deviance. Exploratory analysis of gene targets of microRNAs previously identified as down-regulated due to DGCR8 deficiency suggested that DGCR8 haploinsufficiency may contribute to altered corticogenesis in 22q11DS by disrupting cell cycle modulation. These findings demonstrate the utility of combining neuroanatomic and transcriptomic datasets to derive molecular insights into complex, multigene copy number variants.

Structural brain imaging studies offer clues about the effects of the shared genetic etiology among neuropsychiatric disorders

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.

1q21.1 distal copy number variants are associated with cerebral and cognitive alterations in humans

Low-frequency 1q21.1 distal deletion and duplication copy number variant (CNV) carriers are predisposed to multiple neurodevelopmental disorders, including schizophrenia, autism and intellectual disability. Human carriers display a high prevalence of micro- and macrocephaly in deletion and duplication carriers, respectively. The underlying brain structural diversity remains largely unknown. We systematically called CNVs in 38 cohorts from the large-scale ENIGMA-CNV collaboration and the UK Biobank and identified 28 1q21.1 distal deletion and 22 duplication carriers and 37,088 non-carriers (48% male) derived from 15 distinct magnetic resonance imaging scanner sites. With standardized methods, we compared subcortical and cortical brain measures (all) and cognitive performance (UK Biobank only) between carrier groups also testing for mediation of brain structure on cognition. We identified positive dosage effects of copy number on intracranial volume (ICV) and total cortical surface area, with the largest effects in frontal and cingulate cortices, and negative dosage effects on caudate and hippocampal volumes. The carriers displayed distinct cognitive deficit profiles in cognitive tasks from the UK Biobank with intermediate decreases in duplication carriers and somewhat larger in deletion carriers-the latter potentially mediated by ICV or cortical surface area. These results shed light on pathobiological mechanisms of neurodevelopmental disorders, by demonstrating gene dose effect on specific brain structures and effect on cognitive function.

Association of Immunosuppression and Viral Load With Subcortical Brain Volume in an International Sample of People Living With HIV

IMPORTANCE

Despite more widely accessible combination antiretroviral therapy (cART), HIV-1 infection remains a global public health challenge. Even in treated patients with chronic HIV infection, neurocognitive impairment often persists, affecting quality of life. Identifying the neuroanatomical pathways associated with infection in vivo may delineate the neuropathologic processes underlying these deficits. However, published neuroimaging findings from relatively small, heterogeneous cohorts are inconsistent, limiting the generalizability of the conclusions drawn to date.

OBJECTIVE

To examine structural brain associations with the most commonly collected clinical assessments of HIV burden (CD4+ T-cell count and viral load), which are generalizable across demographically and clinically diverse HIV-infected individuals worldwide.

DESIGN, SETTING, AND PARTICIPANTS

This cross-sectional study established the HIV Working Group within the Enhancing Neuro Imaging Genetics Through Meta Analysis (ENIGMA) consortium to pool and harmonize data from existing HIV neuroimaging studies. In total, data from 1295 HIV-positive adults were contributed from 13 studies across Africa, Asia, Australia, Europe, and North America. Regional and whole brain segmentations were extracted from data sets as contributing studies joined the consortium on a rolling basis from November 1, 2014, to December 31, 2019.

MAIN OUTCOMES AND MEASURES

Volume estimates for 8 subcortical brain regions were extracted from T1-weighted magnetic resonance images to identify associations with blood plasma markers of current immunosuppression (CD4+ T-cell counts) or detectable plasma viral load (dVL) in HIV-positive participants. Post hoc sensitivity analyses stratified data by cART status.

RESULTS

After quality assurance, data from 1203 HIV-positive individuals (mean [SD] age, 45.7 [11.5] years; 880 [73.2%] male; 897 [74.6%] taking cART) remained. Lower current CD4+ cell counts were associated with smaller hippocampal (mean [SE] β = 16.66 [4.72] mm3 per 100 cells/mm3; P < .001) and thalamic (mean [SE] β = 32.24 [8.96] mm3 per 100 cells/mm3; P < .001) volumes and larger ventricles (mean [SE] β = -391.50 [122.58] mm3 per 100 cells/mm3; P = .001); in participants not taking cART, however, lower current CD4+ cell counts were associated with smaller putamen volumes (mean [SE] β = 57.34 [18.78] mm3 per 100 cells/mm3; P = .003). A dVL was associated with smaller hippocampal volumes (d = -0.17; P = .005); in participants taking cART, dVL was also associated with smaller amygdala volumes (d = -0.23; P = .004).

CONCLUSIONS AND RELEVANCE

In a large-scale international population of HIV-positive individuals, volumes of structures in the limbic system were consistently associated with current plasma markers. Our findings extend beyond the classically implicated regions of the basal ganglia and may represent a generalizable brain signature of HIV infection in the cART era.

Virtual Histology of Cortical Thickness and Shared Neurobiology in 6 Psychiatric Disorders

IMPORTANCE

Large-scale neuroimaging studies have revealed group differences in cortical thickness across many psychiatric disorders. The underlying neurobiology behind these differences is not well understood.

OBJECTIVE

To determine neurobiologic correlates of group differences in cortical thickness between cases and controls in 6 disorders: attention-deficit/hyperactivity disorder (ADHD), autism spectrum disorder (ASD), bipolar disorder (BD), major depressive disorder (MDD), obsessive-compulsive disorder (OCD), and schizophrenia.

DESIGN, SETTING, AND PARTICIPANTS

Profiles of group differences in cortical thickness between cases and controls were generated using T1-weighted magnetic resonance images. Similarity between interregional profiles of cell-specific gene expression and those in the group differences in cortical thickness were investigated in each disorder. Next, principal component analysis was used to reveal a shared profile of group difference in thickness across the disorders. Analysis for gene coexpression, clustering, and enrichment for genes associated with these disorders were conducted. Data analysis was conducted between June and December 2019. The analysis included 145 cohorts across 6 psychiatric disorders drawn from the ENIGMA consortium. The numbers of cases and controls in each of the 6 disorders were as follows: ADHD: 1814 and 1602; ASD: 1748 and 1770; BD: 1547 and 3405; MDD: 2658 and 3572; OCD: 2266 and 2007; and schizophrenia: 2688 and 3244.

MAIN OUTCOMES AND MEASURES

Interregional profiles of group difference in cortical thickness between cases and controls.

RESULTS

A total of 12 721 cases and 15 600 controls, ranging from ages 2 to 89 years, were included in this study. Interregional profiles of group differences in cortical thickness for each of the 6 psychiatric disorders were associated with profiles of gene expression specific to pyramidal (CA1) cells, astrocytes (except for BD), and microglia (except for OCD); collectively, gene-expression profiles of the 3 cell types explain between 25% and 54% of variance in interregional profiles of group differences in cortical thickness. Principal component analysis revealed a shared profile of difference in cortical thickness across the 6 disorders (48% variance explained); interregional profile of this principal component 1 was associated with that of the pyramidal-cell gene expression (explaining 56% of interregional variation). Coexpression analyses of these genes revealed 2 clusters: (1) a prenatal cluster enriched with genes involved in neurodevelopmental (axon guidance) processes and (2) a postnatal cluster enriched with genes involved in synaptic activity and plasticity-related processes. These clusters were enriched with genes associated with all 6 psychiatric disorders.

CONCLUSIONS AND RELEVANCE

In this study, shared neurobiologic processes were associated with differences in cortical thickness across multiple psychiatric disorders. These processes implicate a common role of prenatal development and postnatal functioning of the cerebral cortex in these disorders.

Altered white matter microstructure in 22q11.2 deletion syndrome: a multisite diffusion tensor imaging study

22q11.2 deletion syndrome (22q11DS)-a neurodevelopmental condition caused by a hemizygous deletion on chromosome 22-is associated with an elevated risk of psychosis and other developmental brain disorders. Prior single-site diffusion magnetic resonance imaging (dMRI) studies have reported altered white matter (WM) microstructure in 22q11DS, but small samples and variable methods have led to contradictory results. Here we present the largest study ever conducted of dMRI-derived measures of WM microstructure in 22q11DS (334 22q11.2 deletion carriers and 260 healthy age- and sex-matched controls; age range 6-52 years). Using harmonization protocols developed by the ENIGMA-DTI working group, we identified widespread reductions in mean, axial and radial diffusivities in 22q11DS, most pronounced in regions with major cortico-cortical and cortico-thalamic fibers: the corona radiata, corpus callosum, superior longitudinal fasciculus, posterior thalamic radiations, and sagittal stratum (Cohen's d's ranging from -0.9 to -1.3). Only the posterior limb of the internal capsule (IC), comprised primarily of corticofugal fibers, showed higher axial diffusivity in 22q11DS. 22q11DS patients showed higher mean fractional anisotropy (FA) in callosal and projection fibers (IC and corona radiata) relative to controls, but lower FA than controls in regions with predominantly association fibers. Psychotic illness in 22q11DS was associated with more substantial diffusivity reductions in multiple regions. Overall, these findings indicate large effects of the 22q11.2 deletion on WM microstructure, especially in major cortico-cortical connections. Taken together with findings from animal models, this pattern of abnormalities may reflect disrupted neurogenesis of projection neurons in outer cortical layers.

In vivo hippocampal subfield volumes in bipolar disorder—A mega‐analysis from The Enhancing Neuro Imaging Genetics through Meta‐Analysis Bipolar Disorder Working Group

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.

Intelligence, educational attainment, and brain structure in those at familial high-risk for schizophrenia or bipolar disorder

First-degree relatives of patients diagnosed with schizophrenia (SZ-FDRs) show similar patterns of brain abnormalities and cognitive alterations to patients, albeit with smaller effect sizes. First-degree relatives of patients diagnosed with bipolar disorder (BD-FDRs) show divergent patterns; on average, intracranial volume is larger compared to controls, and findings on cognitive alterations in BD-FDRs are inconsistent. Here, we performed a meta-analysis of global and regional brain measures (cortical and subcortical), current IQ, and educational attainment in 5,795 individuals (1,103 SZ-FDRs, 867 BD-FDRs, 2,190 controls, 942 schizophrenia patients, 693 bipolar patients) from 36 schizophrenia and/or bipolar disorder family cohorts, with standardized methods. Compared to controls, SZ-FDRs showed a pattern of widespread thinner cortex, while BD-FDRs had widespread larger cortical surface area. IQ was lower in SZ-FDRs (d = -0.42, p = 3 × 10^-5 ), with weak evidence of IQ reductions among BD-FDRs (d = -0.23, p = .045). Both relative groups had similar educational attainment compared to controls. When adjusting for IQ or educational attainment, the group-effects on brain measures changed, albeit modestly. Changes were in the expected direction, with less pronounced brain abnormalities in SZ-FDRs and more pronounced effects in BD-FDRs. To conclude, SZ-FDRs and BD-FDRs show a differential pattern of structural brain abnormalities. In contrast, both had lower IQ scores and similar school achievements compared to controls. Given that brain differences between SZ-FDRs and BD-FDRs remain after adjusting for IQ or educational attainment, we suggest that differential brain developmental processes underlying predisposition for schizophrenia or bipolar disorder are likely independent of general cognitive impairment.

Large-scale mapping of cortical alterations in 22q11.2 deletion syndrome: Convergence with idiopathic psychosis and effects of deletion size

The 22q11.2 deletion (22q11DS) is a common chromosomal microdeletion and a potent risk factor for psychotic illness. Prior studies reported widespread cortical changes in 22q11DS, but were generally underpowered to characterize neuroanatomic abnormalities associated with psychosis in 22q11DS, and/or neuroanatomic effects of variability in deletion size. To address these issues, we developed the ENIGMA (Enhancing Neuro Imaging Genetics Through Meta-Analysis) 22q11.2 Working Group, representing the largest analysis of brain structural alterations in 22q11DS to date. The imaging data were collected from 10 centers worldwide, including 474 subjects with 22q11DS (age = 18.2 ± 8.6; 46.9% female) and 315 typically developing, matched controls (age = 18.0 ± 9.2; 45.9% female). Compared to controls, 22q11DS individuals showed thicker cortical gray matter overall (left/right hemispheres: Cohen's d = 0.61/0.65), but focal thickness reduction in temporal and cingulate cortex. Cortical surface area (SA), however, showed pervasive reductions in 22q11DS (left/right hemispheres: d = -1.01/-1.02). 22q11DS cases vs. controls were classified with 93.8% accuracy based on these neuroanatomic patterns. Comparison of 22q11DS-psychosis to idiopathic schizophrenia (ENIGMA-Schizophrenia Working Group) revealed significant convergence of affected brain regions, particularly in fronto-temporal cortex. Finally, cortical SA was significantly greater in 22q11DS cases with smaller 1.5 Mb deletions, relative to those with typical 3 Mb deletions. We found a robust neuroanatomic signature of 22q11DS, and the first evidence that deletion size impacts brain structure. Psychotic illness in this highly penetrant deletion was associated with similar neuroanatomic abnormalities to idiopathic schizophrenia. These consistent cross-site findings highlight the homogeneity of this single genetic etiology, and support the suitability of 22q11DS as a biological model of schizophrenia.

Mapping Subcortical Brain Alterations in 22q11.2 Deletion Syndrome: Effects of Deletion Size and Convergence With Idiopathic Neuropsychiatric Illness

OBJECTIVE

22q11.2 deletion syndrome (22q11DS) is among the strongest known genetic risk factors for schizophrenia. Previous studies have reported variable alterations in subcortical brain structures in 22q11DS. To better characterize subcortical alterations in 22q11DS, including modulating effects of clinical and genetic heterogeneity, the authors studied a large multicenter neuroimaging cohort from the ENIGMA 22q11.2 Deletion Syndrome Working Group.

METHODS

Subcortical structures were measured using harmonized protocols for gross volume and subcortical shape morphometry in 533 individuals with 22q11DS and 330 matched healthy control subjects (age range, 6-56 years; 49% female).

RESULTS

Compared with the control group, the 22q11DS group showed lower intracranial volume (ICV) and thalamus, putamen, hippocampus, and amygdala volumes and greater lateral ventricle, caudate, and accumbens volumes (Cohen's d values, -0.90 to 0.93). Shape analysis revealed complex differences in the 22q11DS group across all structures. The larger A-D deletion was associated with more extensive shape alterations compared with the smaller A-B deletion. Participants with 22q11DS with psychosis showed lower ICV and hippocampus, amygdala, and thalamus volumes (Cohen's d values, -0.91 to 0.53) compared with participants with 22q11DS without psychosis. Shape analysis revealed lower thickness and surface area across subregions of these structures. Compared with subcortical findings from other neuropsychiatric disorders studied by the ENIGMA consortium, significant convergence was observed between participants with 22q11DS with psychosis and participants with schizophrenia, bipolar disorder, major depressive disorder, and obsessive-compulsive disorder.

CONCLUSIONS

In the largest neuroimaging study of 22q11DS to date, the authors found widespread alterations to subcortical brain structures, which were affected by deletion size and psychotic illness. Findings indicate significant overlap between 22q11DS-associated psychosis, idiopathic schizophrenia, and other severe neuropsychiatric illnesses.

Association of Copy Number Variation of the 15q11.2 BP1-BP2 Region With Cortical and Subcortical Morphology and Cognition

Importance

Recurrent microdeletions and duplications in the genomic region 15q11.2 between breakpoints 1 (BP1) and 2 (BP2) are associated with neurodevelopmental disorders. These structural variants are present in 0.5% to 1.0% of the population, making 15q11.2 BP1-BP2 the site of the most prevalent known pathogenic copy number variation (CNV). It is unknown to what extent this CNV influences brain structure and affects cognitive abilities.

Objective

To determine the association of the 15q11.2 BP1-BP2 deletion and duplication CNVs with cortical and subcortical brain morphology and cognitive task performance.

Design, Setting, and Participants

In this genetic association study, T1-weighted brain magnetic resonance imaging were combined with genetic data from the ENIGMA-CNV consortium and the UK Biobank, with a replication cohort from Iceland. In total, 203 deletion carriers, 45 247 noncarriers, and 306 duplication carriers were included. Data were collected from August 2015 to April 2019, and data were analyzed from September 2018 to September 2019.

Main Outcomes and Measures

The associations of the CNV with global and regional measures of surface area and cortical thickness as well as subcortical volumes were investigated, correcting for age, age2, sex, scanner, and intracranial volume. Additionally, measures of cognitive ability were analyzed in the full UK Biobank cohort.

Results

Of 45 756 included individuals, the mean (SD) age was 55.8 (18.3) years, and 23 754 (51.9%) were female. Compared with noncarriers, deletion carriers had a lower surface area (Cohen d = -0.41; SE, 0.08; P = 4.9 × 10-8), thicker cortex (Cohen d = 0.36; SE, 0.07; P = 1.3 × 10-7), and a smaller nucleus accumbens (Cohen d = -0.27; SE, 0.07; P = 7.3 × 10-5). There was also a significant negative dose response on cortical thickness (β = -0.24; SE, 0.05; P = 6.8 × 10-7). Regional cortical analyses showed a localization of the effects to the frontal, cingulate, and parietal lobes. Further, cognitive ability was lower for deletion carriers compared with noncarriers on 5 of 7 tasks.

Conclusions and Relevance

These findings, from the largest CNV neuroimaging study to date, provide evidence that 15q11.2 BP1-BP2 structural variation is associated with brain morphology and cognition, with deletion carriers being particularly affected. The pattern of results fits with known molecular functions of genes in the 15q11.2 BP1-BP2 region and suggests involvement of these genes in neuronal plasticity. These neurobiological effects likely contribute to the association of this CNV with neurodevelopmental disorders.

Subcortical shape alterations in major depressive disorder: Findings from the ENIGMA major depressive disorder working group

Alterations in regional subcortical brain volumes have been investigated as part of the efforts of an international consortium, ENIGMA, to identify reliable neural correlates of major depressive disorder (MDD). Given that subcortical structures are comprised of distinct subfields, we sought to build significantly from prior work by precisely mapping localized MDD‐related differences in subcortical regions using shape analysis. In this meta‐analysis of subcortical shape from the ENIGMA‐MDD working group, we compared 1,781 patients with MDD and 2,953 healthy controls (CTL) on individual measures of shape metrics (thickness and surface area) on the surface of seven bilateral subcortical structures: nucleus accumbens, amygdala, caudate, hippocampus, pallidum, putamen, and thalamus. Harmonized data processing and statistical analyses were conducted locally at each site, and findings were aggregated by meta‐analysis. Relative to CTL, patients with adolescent‐onset MDD (≤ 21 years) had lower thickness and surface area of the subiculum, cornu ammonis (CA) 1 of the hippocampus and basolateral amygdala (Cohen's d = −0.164 to −0.180). Relative to first‐episode MDD, recurrent MDD patients had lower thickness and surface area in the CA1 of the hippocampus and the basolateral amygdala (Cohen's d = −0.173 to −0.184). Our results suggest that previously reported MDD‐associated volumetric differences may be localized to specific subfields of these structures that have been shown to be sensitive to the effects of stress, with important implications for mapping treatments to patients based on specific neural targets and key clinical features.

The genetic architecture of the human cerebral cortex

The cerebral cortex underlies our complex cognitive capabilities, yet little is known about the specific genetic loci that influence human cortical structure. To identify genetic variants that affect cortical structure, we conducted a genome-wide association meta-analysis of brain magnetic resonance imaging data from 51,665 individuals. We analyzed the surface area and average thickness of the whole cortex and 34 regions with known functional specializations. We identified 199 significant loci and found significant enrichment for loci influencing total surface area within regulatory elements that are active during prenatal cortical development, supporting the radial unit hypothesis. Loci that affect regional surface area cluster near genes in Wnt signaling pathways, which influence progenitor expansion and areal identity. Variation in cortical structure is genetically correlated with cognitive function, Parkinson's disease, insomnia, depression, neuroticism, and attention deficit hyperactivity disorder.

ENIGMA and global neuroscience: A decade of large-scale studies of the brain in health and disease across more than 40 countries

This review summarizes the last decade of work by the ENIGMA (Enhancing NeuroImaging Genetics through Meta Analysis) Consortium, a global alliance of over 1400 scientists across 43 countries, studying the human brain in health and disease. Building on large-scale genetic studies that discovered the first robustly replicated genetic loci associated with brain metrics, ENIGMA has diversified into over 50 working groups (WGs), pooling worldwide data and expertise to answer fundamental questions in neuroscience, psychiatry, neurology, and genetics. Most ENIGMA WGs focus on specific psychiatric and neurological conditions, other WGs study normal variation due to sex and gender differences, or development and aging; still other WGs develop methodological pipelines and tools to facilitate harmonized analyses of "big data" (i.e., genetic and epigenetic data, multimodal MRI, and electroencephalography data). These international efforts have yielded the largest neuroimaging studies to date in schizophrenia, bipolar disorder, major depressive disorder, post-traumatic stress disorder, substance use disorders, obsessive-compulsive disorder, attention-deficit/hyperactivity disorder, autism spectrum disorders, epilepsy, and 22q11.2 deletion syndrome. More recent ENIGMA WGs have formed to study anxiety disorders, suicidal thoughts and behavior, sleep and insomnia, eating disorders, irritability, brain injury, antisocial personality and conduct disorder, and dissociative identity disorder. Here, we summarize the first decade of ENIGMA's activities and ongoing projects, and describe the successes and challenges encountered along the way. We highlight the advantages of collaborative large-scale coordinated data analyses for testing reproducibility and robustness of findings, offering the opportunity to identify brain systems involved in clinical syndromes across diverse samples and associated genetic, environmental, demographic, cognitive, and psychosocial factors.

Altered Cortical Brain Structure and Increased Risk for Disease Seen Decades After Perinatal Exposure to Maternal Smoking: A Study of 9000 Adults in the UK Biobank

Secondhand smoke exposure is a major public health risk that is especially harmful to the developing brain, but it is unclear if early exposure affects brain structure during middle age and older adulthood. Here we analyzed brain MRI data from the UK Biobank in a population-based sample of individuals (ages 44-80) who were exposed (n = 2510) or unexposed (n = 6079) to smoking around birth. We used robust statistical models, including quantile regressions, to test the effect of perinatal smoke exposure (PSE) on cortical surface area (SA), thickness, and subcortical volumes. We hypothesized that PSE would be associated with cortical disruption in primary sensory areas compared to unexposed (PSE-) adults. After adjusting for multiple comparisons, SA was significantly lower in the pericalcarine (PCAL), inferior parietal (IPL), and regions of the temporal and frontal cortex of PSE+ adults; these abnormalities were associated with increased risk for several diseases, including circulatory and endocrine conditions. Sensitivity analyses conducted in a hold-out group of healthy participants (exposed, n = 109, unexposed, n = 315) replicated the effect of PSE on SA in the PCAL and IPL. Collectively our results show a negative, long term effect of PSE on sensory cortices that may increase risk for disease later in life.