The two hemispheres of the human brain differ functionally and structurally. Despite over a century of research, the extent to which brain asymmetry is influenced by sex, handedness, age, and genetic factors is still controversial. Here we present the largest ever analysis of subcortical brain asymmetries, in a harmonized multi-site study using meta-analysis methods. Volumetric asymmetry of seven subcortical structures was assessed in 15,847 MRI scans from 52 datasets worldwide. There were sex differences in the asymmetry of the globus pallidus and putamen. Heritability estimates, derived from 1170 subjects belonging to 71 extended pedigrees, revealed that additive genetic factors influenced the asymmetry of these two structures and that of the hippocampus and thalamus. Handedness had no detectable effect on subcortical asymmetries, even in this unprecedented sample size, but the asymmetry of the putamen varied with age. Genetic drivers of asymmetry in the hippocampus, thalamus and basal ganglia may affect variability in human cognition, including susceptibility to psychiatric disorders.

Alcohol use disorder (AUD) and cannabis use disorder (CUD) are associated with brain alterations particularly involving fronto-cerebellar and meso-cortico-limbic circuitry. However, such abnormalities have additionally been reported in other psychiatric conditions, and until recently there has been few large-scale investigations to compare such findings. The current study uses the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) consortium method of standardising structural brain measures to quantify case-control differences and to compare brain-correlates of substance use disorders with those published in relation to other psychiatric disorders. Using the ENIGMA protocols, we report effect sizes derived from a meta-analysis of alcohol (seven studies, N = 798, 54% are cases) and cannabis (seven studies, N = 447, 45% are cases) dependent cases and age- and sex-matched controls. We conduct linear analyses using harmonised methods to process and parcellate brain data identical to those reported in the literature for ENIGMA case-control studies of major depression disorder (MDD), schizophrenia (SCZ) and bipolar disorder so that effect sizes are optimally comparable across disorders. R elationships between substance use disorder diagnosis and subcortical grey matter volumes and cortical thickness were assessed with intracranial volume, age and sex as co-variates . After correcting for multiple comparisons, AUD case-control meta-analysis of subcortical regions indicated significant differences in the thalamus, hippocampus, amygdala and accumbens, with effect sizes (0.23) generally equivalent to, or larger than |0.23| those previously reported for other psychiatric disorders (except for the pallidum and putamen). On measures of cortical thickness, AUD was associated with significant differences bilaterally in the fusiform gyrus, inferior temporal gyrus, temporal pole, superior frontal gyrus, and rostral and caudal anterior cingulate gyri. Meta-analysis of CUD case-control studies indicated reliable reductions in amygdala, accumbens and hippocampus volumes, with the former effect size comparable to, and the latter effect size around half of that reported for alcohol and SCZ. CUD was associated with lower cortical thickness in the frontal regions, particularly the medial orbitofrontal region, but this effect was not significant after correcting for multiple testing. This study allowed for an unbiased cross-disorder comparison of brain correlates of substance use disorders and showed alcohol-related brain anomalies equivalent in effect size to that found in SCZ in several subcortical and cortical regions and significantly greater alterations than those found in MDD in several subcortical and cortical regions. Although modest, CUD results overlapped with findings reported for AUD and other psychiatric conditions, but appear to be most robustly related to reduce thickness of the medial orbitofrontal cortex.

Understanding the modularity of functional magnetic resonance imaging (fMRI)-derived brain networks or "connectomes" can inform the study of brain function organization. However, fMRI connectomes additionally involve negative edges, which may not be optimally accounted for by existing approaches to modularity that variably threshold, binarize, or arbitrarily weight these connections. Consequently, many existing Q maximization-based modularity algorithms yield variable modular structures. Here, we present an alternative complementary approach that exploits how frequent the blood-oxygen-level-dependent (BOLD) signal correlation between two nodes is negative. We validated this novel probability-based modularity approach on two independent publicly-available resting-state connectome data sets (the Human Connectome Project [HCP] and the 1,000 functional connectomes) and demonstrated that negative correlations alone are sufficient in understanding resting-state modularity. In fact, this approach (a) permits a dual formulation, leading to equivalent solutions regardless of whether one considers positive or negative edges; (b) is theoretically linked to the Ising model defined on the connectome, thus yielding modularity result that maximizes data likelihood. Additionally, we were able to detect novel and consistent sex differences in modularity in both data sets. As data sets like HCP become widely available for analysis by the neuroscience community at large, alternative and perhaps more advantageous computational tools to understand the neurobiological information of negative edges in fMRI connectomes are increasingly important.

This Special Issue of Human Brain Mapping is dedicated to a 10-year anniversary of the Enhancing NeuroImaging Genetics through Meta-Analysis (ENIGMA) Consortium. It reports updates from a broad range of international neuroimaging projects that pool data from around the world to answer fundamental questions in neuroscience. Since ENIGMA was formed in December 2009, the initiative grew into a worldwide effort with over 2,000 participating scientists from 45 countries, and over 50 working groups leading large-scale studies of human brain disorders. Over the last decade, many lessons were learned on how best to pool brain data from diverse sources. Working groups were created to develop methods to analyze worldwide data from anatomical and diffusion magnetic resonance imaging (MRI), resting state and task-based functional MRI, electroencephalography (EEG), magnetoencephalography (MEG), and magnetic resonance spectroscopy (MRS). The quest to understand genetic effects on human brain development and disease also led to analyses of brain scans on an unprecedented scale. Genetic roadmaps of the human cortex were created by researchers worldwide who collaborated to perform statistically well-powered analyses of common and rare genetic variants on brain measures and rates of brain development and aging. Here, we summarize the 31 papers in this Special Issue, covering: (a) technical approaches to harmonize analysis of different types of brain imaging data, (b) reviews of the last decade of work by several of ENIGMA's clinical and technical working groups, and (c) new empirical papers reporting large-scale international brain mapping analyses in patients with substance use disorders, schizophrenia, bipolar disorders, major depression, posttraumatic stress disorder, obsessive compulsive disorder, epilepsy, and stroke.

The lack of robust neuroanatomical markers of psychosis risk has been traditionally attributed to heterogeneity. A complementary hypothesis is that variation in neuroanatomical measures in individuals at psychosis risk may be nested within the range observed in healthy individuals.

Microstructural alterations in cortico-subcortical connections are thought to be present in obsessive-compulsive disorder (OCD). However, prior studies have yielded inconsistent findings, perhaps because small sample sizes provided insufficient power to detect subtle abnormalities. Here we investigated microstructural white matter alterations and their relation to clinical features in the largest dataset of adult and pediatric OCD to date. We analyzed diffusion tensor imaging metrics from 700 adult patients and 645 adult controls, as well as 174 pediatric patients and 144 pediatric controls across 19 sites participating in the ENIGMA OCD Working Group, in a cross-sectional case-control magnetic resonance study. We extracted measures of fractional anisotropy (FA) as main outcome, and mean diffusivity, radial diffusivity, and axial diffusivity as secondary outcomes for 25 white matter regions. We meta-analyzed patient-control group differences (Cohen's d) across sites, after adjusting for age and sex, and investigated associations with clinical characteristics. Adult OCD patients showed significant FA reduction in the sagittal stratum (d = -0.21, z = -3.21, p = 0.001) and posterior thalamic radiation (d = -0.26, z = -4.57, p < 0.0001). In the sagittal stratum, lower FA was associated with a younger age of onset (z = 2.71, p = 0.006), longer duration of illness (z = -2.086, p = 0.036), and a higher percentage of medicated patients in the cohorts studied (z = -1.98, p = 0.047). No significant association with symptom severity was found. Pediatric OCD patients did not show any detectable microstructural abnormalities compared to controls. Our findings of microstructural alterations in projection and association fibers to posterior brain regions in OCD are consistent with models emphasizing deficits in connectivity as an important feature of this disorder.

Schizophrenia and bipolar disorder share genetic liability, and some structural brain abnormalities are common to both conditions. First-degree relatives of patients with schizophrenia (FDRs-SZ) show similar brain abnormalities to patients, albeit with smaller effect sizes. Imaging findings in first-degree relatives of patients with bipolar disorder (FDRs-BD) have been inconsistent in the past, but recent studies report regionally greater volumes compared with control subjects.

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.

We present an empirically benchmarked framework for sex-specific normative modeling of brain morphometry that can inform about the biological and behavioral significance of deviations from typical age-related neuroanatomical changes and support future study designs. This framework was developed using regional morphometric data from 37,407 healthy individuals (53% female; aged 3-90 years) following a comparative evaluation of eight algorithms and multiple covariate combinations pertaining to image acquisition and quality, parcellation software versions, global neuroimaging measures, and longitudinal stability. The Multivariate Factorial Polynomial Regression (MFPR) emerged as the preferred algorithm optimized using nonlinear polynomials for age and linear effects of global measures as covariates. The MFPR models showed excellent accuracy across the lifespan and within distinct age-bins, and longitudinal stability over a 2-year period. The performance of all MFPR models plateaued at sample sizes exceeding 3,000 study participants. The model and scripts described here are freely available through CentileBrain (https://centilebrain.org/).

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

Identifying brain alterations associated with suicidal thoughts and behaviors (STBs) in young people is critical to understanding their development and improving early intervention and prevention. The ENIGMA Suicidal Thoughts and Behaviours (ENIGMA-STB) consortium analyzed neuroimaging data harmonized across sites to examine brain morphology associated with STBs in youth. We performed analyses in three separate stages, in samples ranging from most to least homogeneous in terms of suicide assessment instrument and mental disorder. First, in a sample of 577 young people with mood disorders, in which STBs were assessed with the Columbia Suicide Severity Rating Scale (C-SSRS). Second, in a sample of young people with mood disorders, in which STB were assessed using different instruments, MRI metrics were compared among healthy controls without STBs (HC; N = 519), clinical controls with a mood disorder but without STBs (CC; N = 246) and young people with current suicidal ideation (N = 223). In separate analyses, MRI metrics were compared among HCs (N = 253), CCs (N = 217), and suicide attempters (N = 64). Third, in a larger transdiagnostic sample with various assessment instruments (HC = 606; CC = 419; Ideation = 289; HC = 253; CC = 432; Attempt=91). In the homogeneous C-SSRS sample, surface area of the frontal pole was lower in young people with mood disorders and a history of actual suicide attempts (N = 163) than those without a lifetime suicide attempt (N = 323; FDR-p = 0.035, Cohen's d = 0.34). No associations with suicidal ideation were found. When examining more heterogeneous samples, we did not observe significant associations. Lower frontal pole surface area may represent a vulnerability for a (non-interrupted and non-aborted) suicide attempt; however, more research is needed to understand the nature of its relationship to suicide risk.

For many traits, males show greater variability than females, with possible implications for understanding sex differences in health and disease. Here, the ENIGMA (Enhancing Neuro Imaging Genetics through Meta-Analysis) Consortium presents the largest-ever mega-analysis of sex differences in variability of brain structure, based on international data spanning nine decades of life. Subcortical volumes, cortical surface area and cortical thickness were assessed in MRI data of 16,683 healthy individuals 1-90 years old (47% females). We observed significant patterns of greater male than female between-subject variance for all subcortical volumetric measures, all cortical surface area measures, and 60% of cortical thickness measures. This pattern was stable across the lifespan for 50% of the subcortical structures, 70% of the regional area measures, and nearly all regions for thickness. Our findings that these sex differences are present in childhood implicate early life genetic or gene-environment interaction mechanisms. The findings highlight the importance of individual differences within the sexes, that may underpin sex-specific vulnerability to disorders.

The NTRK3 gene (also known as TRKC) encodes a high affinity receptor for the neurotrophin 3'-nucleotidase (NT3), which is implicated in oligodendrocyte and myelin development. We previously found that white matter integrity in young adults is related to common variants in genes encoding neurotrophins and their receptors. This underscores the importance of neurotrophins for white matter development. NTRK3 variants are putative risk factors for schizophrenia, bipolar disorder, and obsessive-compulsive disorder hoarding, suggesting that some NTRK3 variants may affect the brain. To test this, we scanned 392 healthy adult twins and their siblings (mean age, 23.6 ± 2.2 years; range: 20-29 years) with 105-gradient 4-Tesla diffusion tensor imaging (DTI). We identified 18 single nucleotide polymorphisms (SNPs) in the NTRK3 gene that have been associated with neuropsychiatric disorders. We used a multi-SNP model, adjusting for family relatedness, age, and sex, to relate these variants to voxelwise fractional anisotropy (FA) - a DTI measure of white matter integrity. FA was optimally predicted (based on the highest false discovery rate critical p), by five SNPs (rs1017412, rs2114252, rs16941261, rs3784406, and rs7176429; overall FDR critical p=0.028). Gene effects were widespread and included the corpus callosum genu and inferior longitudinal fasciculus - regions implicated in several neuropsychiatric disorders and previously associated with other neurotrophin-related genetic variants in an overlapping sample of subjects. NTRK3 genetic variants, and neurotrophins more generally, may influence white matter integrity in brain regions implicated in neuropsychiatric disorders.

Reciprocal chromosomal rearrangements at the 22q11.2 locus are associated with elevated risk of neurodevelopmental disorders. The 22q11.2 deletion confers the highest known genetic risk for schizophrenia, but a duplication in the same region is strongly associated with autism and is less common in schizophrenia cases than in the general population. Here we conducted the first study of 22q11.2 gene dosage effects on brain structure in a sample of 143 human subjects: 66 with 22q11.2 deletions (22q-del; 32 males), 21 with 22q11.2 duplications (22q-dup; 14 males), and 56 age- and sex-matched controls (31 males). 22q11.2 gene dosage varied positively with intracranial volume, gray and white matter volume, and cortical surface area (deletion < control < duplication). In contrast, gene dosage varied negatively with mean cortical thickness (deletion > control > duplication). Widespread differences were observed for cortical surface area with more localized effects on cortical thickness. These diametric patterns extended into subcortical regions: 22q-dup carriers had a significantly larger right hippocampus, on average, but lower right caudate and corpus callosum volume, relative to 22q-del carriers. Novel subcortical shape analysis revealed greater radial distance (thickness) of the right amygdala and left thalamus, and localized increases and decreases in subregions of the caudate, putamen, and hippocampus in 22q-dup relative to 22q-del carriers. This study provides the first evidence that 22q11.2 is a genomic region associated with gene-dose-dependent brain phenotypes. Pervasive effects on cortical surface area imply that this copy number variant affects brain structure early in the course of development.SIGNIFICANCE STATEMENT Probing naturally occurring reciprocal copy number variation in the genome may help us understand mechanisms underlying deviations from typical brain and cognitive development. The 22q11.2 genomic region is particularly susceptible to chromosomal rearrangements and contains many genes crucial for neuronal development and migration. Not surprisingly, reciprocal genomic imbalances at this locus confer some of the highest known genetic risks for developmental neuropsychiatric disorders. Here we provide the first evidence that brain morphology differs meaningfully as a function of reciprocal genomic variation at the 22q11.2 locus. Cortical thickness and surface area were affected in opposite directions with more widespread effects of gene dosage on cortical surface area.

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.

Larger thalamic volume has been found in children with obsessive-compulsive disorder (OCD) and children with clinical-level symptoms within the general population. Particular thalamic subregions may drive these differences. The ENIGMA-OCD working group conducted mega- and meta-analyses to study thalamic subregional volume in OCD across the lifespan. Structural T1-weighted brain magnetic resonance imaging (MRI) scans from 2649 OCD patients and 2774 healthy controls across 29 sites (50 datasets) were processed using the FreeSurfer built-in ThalamicNuclei pipeline to extract five thalamic subregions. Volume measures were harmonized for site effects using ComBat before running separate multiple linear regression models for children, adolescents, and adults to estimate volumetric group differences. All analyses were pre-registered ( https://osf.io/73dvy ) and adjusted for age, sex and intracranial volume. Unmedicated pediatric OCD patients (<12 years) had larger lateral (d = 0.46), pulvinar (d = 0.33), ventral (d = 0.35) and whole thalamus (d = 0.40) volumes at unadjusted p-values <0.05. Adolescent patients showed no volumetric differences. Adult OCD patients compared with controls had smaller volumes across all subregions (anterior, lateral, pulvinar, medial, and ventral) and smaller whole thalamic volume (d = -0.15 to -0.07) after multiple comparisons correction, mostly driven by medicated patients and associated with symptom severity. The anterior thalamus was also significantly smaller in patients after adjusting for thalamus size. Our results suggest that OCD-related thalamic volume differences are global and not driven by particular subregions and that the direction of effects are driven by both age and medication status.

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.

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.