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

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 (SCZ). DESIGN, SETTING, AND PARTICIPANTSProfiles 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. 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 number of cases and controls in each of the 6 disorders were as follows:

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

Complex systems are commonly characterized by the properties of their graph representation. Dynamical complex systems are then typically represented by a graph of temporal dependencies between time series of state variables of their subunits. It has been shown recently that graphs constructed in this way tend to have relatively clustered structure, potentially leading to spurious detection of small-world properties even in the case of systems with no or randomly distributed true interactions. However, the strength of this bias depends heavily on a range of parameters and its relevance for real-world data has not yet been established. In this work, we assess the relevance of the bias using two examples of multivariate time series recorded in natural complex systems. The first is the time series of local brain activity as measured by functional magnetic resonance imaging in resting healthy human subjects, and the second is the time series of average monthly surface air temperature coming from a large reanalysis of climatological data over the period 1948-2012. In both cases, the clustering in the thresholded correlation graph is substantially higher compared with a realization of a density-matched random graph, while the shortest paths are relatively short, showing thus distinguishing features of small-world structure. However, comparable or even stronger small-world properties were reproduced in correlation graphs of model processes with randomly scrambled interconnections. This suggests that the small-world properties of the correlation matrices of these real-world systems indeed do not reflect genuinely the properties of the underlying interaction structure, but rather result from the inherent properties of correlation matrix.