Reply to: New Meta- and Mega-analyses of Magnetic Resonance Imaging Findings in Schizophrenia: Do They Really Increase Our Knowledge About the Nature of the Disease Process?

This work was supported by National Institute of Biomedical Imaging and Bioengineering Grant No. U54EB020403 (to the ENIGMA consortium)

Using brain structural neuroimaging measures to predict psychosis onset for individuals at clinical high-risk

Machine learning approaches using structural magnetic resonance imaging (sMRI) can be informative for disease classification, although their ability to predict psychosis is largely unknown. We created a model with individuals at CHR who developed psychosis later (CHR-PS+) from healthy controls (HCs) that can differentiate each other. We also evaluated whether we could distinguish CHR-PS+ individuals from those who did not develop psychosis later (CHR-PS-) and those with uncertain follow-up status (CHR-UNK). T1-weighted structural brain MRI scans from 1165 individuals at CHR (CHR-PS+, n = 144; CHR-PS-, n = 793; and CHR-UNK, n = 228), and 1029 HCs, were obtained from 21 sites. We used ComBat to harmonize measures of subcortical volume, cortical thickness and surface area data and corrected for non-linear effects of age and sex using a general additive model. CHR-PS+ (n = 120) and HC (n = 799) data from 20 sites served as a training dataset, which we used to build a classifier. The remaining samples were used external validation datasets to evaluate classifier performance (test, independent confirmatory, and independent group [CHR-PS- and CHR-UNK] datasets). The accuracy of the classifier on the training and independent confirmatory datasets was 85% and 73% respectively. Regional cortical surface area measures-including those from the right superior frontal, right superior temporal, and bilateral insular cortices strongly contributed to classifying CHR-PS+ from HC. CHR-PS- and CHR-UNK individuals were more likely to be classified as HC compared to CHR-PS+ (classification rate to HC: CHR-PS+, 30%; CHR-PS-, 73%; CHR-UNK, 80%). We used multisite sMRI to train a classifier to predict psychosis onset in CHR individuals, and it showed promise predicting CHR-PS+ in an independent sample. The results suggest that when considering adolescent brain development, baseline MRI scans for CHR individuals may be helpful to identify their prognosis. Future prospective studies are required about whether the classifier could be actually helpful in the clinical settings.</p

Increased power by harmonizing structural MRI site differences with the ComBat batch adjustment method in ENIGMA

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)

Large-scale analysis of structural brain asymmetries in schizophrenia via the ENIGMA consortium

BACKGROUND Left-right asymmetry is an important organizing feature of the healthy brain that may be altered in schizophrenia, but most studies have used relatively small samples and heterogeneous approaches, resulting in equivocal findings. We carried out the largest case-control study of structural brain asymmetries in schizophrenia (N = 11,095), using a single image analysis protocol. METHODS We included T1-weighted data from 46 datasets (5,080 affected individuals and 6,015 controls) from the ENIGMA Consortium. Asymmetry indexes were calculated for global and regional cortical thickness, surface area, and subcortical volume measures. Differences of asymmetry were calculated between affected individuals and controls per dataset, and effect sizes were meta-analyzed across datasets. Analyses were also performed with respect to the use of antipsychotic medication and other clinical variables, as well as age and sex. Case-control differences in a multivariate context were assessed in a subset of the data (N = 2,029). RESULTS Small average differences between cases and controls were observed for asymmetries in cortical thickness, specifically of the rostral anterior cingulate (d = −0.08, pFDR = 0.047) and the middle temporal gyrus (d = −0.07, pFDR = 0.048), both driven primarily by thinner cortices in the left hemisphere in schizophrenia. These asymmetries were not significantly associated with the use of antipsychotic medication or other clinical variables. Older individuals with schizophrenia showed a stronger average leftward asymmetry of pallidum volume than older controls (d = 0.08, pFDR = 9.0 × 10−3). The multivariate analysis revealed that 7% of the variance across all structural asymmetries was explained by case-control status (F = 1.87, p = 1.25 × 10−5). CONCLUSIONS Altered trajectories of asymmetrical brain development and/or lifespan asymmetry may contribute to schizophrenia pathophysiology. Small case-control differences of brain macro-structural asymmetry may manifest due to more substantial differences at the molecular, cytoarchitectonic or circuit levels, with functional relevance for lateralized cognitive processes

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

Neuroanatomical heterogeneity and homogeneity in individuals at clinical high risk for psychosis.

Individuals at Clinical High Risk for Psychosis (CHR-P) demonstrate heterogeneity in clinical profiles and outcome features. However, the extent of neuroanatomical heterogeneity in the CHR-P state is largely undetermined. We aimed to quantify the neuroanatomical heterogeneity in structural magnetic resonance imaging measures of cortical surface area (SA), cortical thickness (CT), subcortical volume (SV), and intracranial volume (ICV) in CHR-P individuals compared with healthy controls (HC), and in relation to subsequent transition to a first episode of psychosis. The ENIGMA CHR-P consortium applied a harmonised analysis to neuroimaging data across 29 international sites, including 1579 CHR-P individuals and 1243 HC, offering the largest pooled CHR-P neuroimaging dataset to date. Regional heterogeneity was indexed with the Variability Ratio (VR) and Coefficient of Variation (CV) ratio applied at the group level. Personalised estimates of heterogeneity of SA, CT and SV brain profiles were indexed with the novel Person-Based Similarity Index (PBSI), with two complementary applications. First, to assess the extent of within-diagnosis similarity or divergence of neuroanatomical profiles between individuals. Second, using a normative modelling approach, to assess the 'normativeness' of neuroanatomical profiles in individuals at CHR-P. CHR-P individuals demonstrated no greater regional heterogeneity after applying FDR corrections. However, PBSI scores indicated significantly greater neuroanatomical divergence in global SA, CT and SV profiles in CHR-P individuals compared with HC. Normative PBSI analysis identified 11 CHR-P individuals (0.70%) with marked deviation (>1.5 SD) in SA, 118 (7.47%) in CT and 161 (10.20%) in SV. Psychosis transition was not significantly associated with any measure of heterogeneity. Overall, our examination of neuroanatomical heterogeneity within the CHR-P state indicated greater divergence in neuroanatomical profiles at an individual level, irrespective of psychosis conversion. Further large-scale investigations are required of those who demonstrate marked deviation

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

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

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

White matter disturbances in major depressive disorder : a coordinated analysis across 20 international cohorts in the ENIGMA MDD working group

Altres ajuts: The ENIGMA-Major Depressive Disorder working group gratefully acknowledges support from the NIH Big Data to Knowledge (BD2K) award (U54 EB020403 to PMT) and NIH grant R01 MH116147 (PMT). LS is supported by an NHMRC MRFF Career Development Fellowship (APP1140764). We wish to acknowledge the patients and control subjects that have particiaped int the study. We thank Rosa Schirmer, Elke Schreiter, Reinhold Borschke and Ines Eidner for image acquisition and data preparation, and Anna Oliynyk for quality checks. We thank Dorothee P. Auer and F. Holsboer for initiation of the RUD study. We wish to acknowledge the patients and control subjects that have particiaped int the study. We thank Rosa Schirmer, Elke Schreiter, Reinhold Borschke and Ines Eidner for image acquisition and data preparation, and Anna Oliynyk for quality checks. We thank Dorothee P. Auer and F. Holsboer for initiation of the RUD study. NESDA: The infrastructure for the NESDA study (www.nesda.nl) is funded through the Geestkracht program of the Netherlands Organisation for Health Research and Development (Zon-Mw, grant number 10-000-1002) and is supported by participating universities (VU University Medical Center, GGZ inGeest, Arkin, Leiden University Medical Center, GGZ Rivierduinen, University Medical Center Groningen) and mental health care organizations, see www.nesda.nl. M-JvT was supported by a VENI grant (NWO grant number 016.156.077). UCSF: This work was supported by the Brain and Behavior Research Foundation (formerly NARSAD) to TTY; the National Institute of Mental Health (R01MH085734 to TTY; K01MH117442 to TCH) and by the American Foundation for Suicide Prevention (PDF-1-064-13) to TCH. Stanford: This work was supported by NIMH Grants R01MH59259 and R37101495 to IHG. MS is partially supported by an award funded by the Phyllis and Jerome Lyle Rappaport Foundation. Muenster: This work was funded by the German Research Foundation (SFB-TRR58, Projects C09 and Z02 to UD) and the Interdisciplinary Center for Clinical Research (IZKF) of the medical faculty of Münster (grant Dan3/012/17 to UD). Marburg: This work was funded by the German Research Foundation (DFG, grant FOR2107 DA1151/5-1 and DA1151/5-2 to UD; KI 588/ 14-1, KI 588/14-2 to TK; KR 3822/7-1, KR 3822/7-2 to AK; JA 1890/ 7-1, JA 1890/7-2 to AJ). IMH-MDD: This work was supported by the National Healthcare Group Research Grant (SIG/15012) awarded to KS. Barcelona: This study was funded by two grants of the Fondo de Investigación Sanitaria from the Instituto de Salud Carlos III, by the Centro de Investigación Biomédica en Red de Salud Mental (CIBERSAM). The author is funded through 'Miguel Servet' research contract (CP16-0020), co-financed by the European Regional Development Fund (ERDF) (2016-2019). QTIM: We thank the twins and singleton siblings who gave generously of their time to participate in the QTIM study. We also thank the many research assistants, radiographers, and IT support staff for data acquisition and DNA sample preparation. This study was funded by White matter disturbances in major depressive disorder: a coordinated analysis across 20 international. . . 1521 the National Institute of Child Health & Human Development (RO1 HD050735); National Institute of Biomedical Imaging and Bioengineering (Award 1U54EB020403-01, Subaward 56929223); National Health and Medical Research Council, Australia (Project Grants 496682, 1009064). NIH ENIGMA-BD2K U54 EB020403 (Thompson); R01 MH117601 (Jahanshad/Schmaal). Magdeburg: M.L. and M.W. are funded by SFB 779. Bipolar Family Study: This study has received funding from the European Community's Seventh Framework Programme (FP7/2007-2013). This paper reflects only the author's views and the European Union is not liable for any use that may be made of the information contained therein. This work was also supported by a Wellcome Trust Strategic Award (104036/Z/14/Z). Minnesota Adolescent Depression Study: The study was funded by the National Institute of Mental Health (K23MH090421), the National Alliance for Research on Schizophrenia and Depression, the University of Minnesota Graduate School, the Minnesota Medical Foundation, and the Biotechnology Research Center (P41 RR008079 to the Center for Magnetic Resonance Research), University of Minnesota, and the Deborah E. Powell Center for Women's Health Seed Grant, University of Minnesota. Dublin: This study was supported by Science Foundation Ireland through a Stokes Professorhip grant to TF. MPIP: The MPIP Sample comprises patients included in the Recurrent Unipolar Depression (RUD) Case-Control study at the clinic of the Max Planck Institute of Psychiatry, Munich, German. The RUD study was supported by GlaxoSmithKline.Alterations in white matter (WM) microstructure have been implicated in the pathophysiology of major depressive disorder (MDD). However, previous findings have been inconsistent, partially due to low statistical power and the heterogeneity of depression. In the largest multi-site study to date, we examined WM anisotropy and diffusivity in 1305 MDD patients and 1602 healthy controls (age range 12-88 years) from 20 samples worldwide, which included both adults and adolescents, within the MDD Working Group of the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) consortium. Processing of diffusion tensor imaging (DTI) data and statistical analyses were harmonized across sites and effects were meta-analyzed across studies. We observed subtle, but widespread, lower fractional anisotropy (FA) in adult MDD patients compared with controls in 16 out of 25 WM tracts of interest (Cohen's d between 0.12 and 0.26). The largest differences were observed in the corpus callosum and corona radiata. Widespread higher radial diffusivity (RD) was also observed (all Cohen's d between 0.12 and 0.18). Findings appeared to be driven by patients with recurrent MDD and an adult age of onset of depression. White matter microstructural differences in a smaller sample of adolescent MDD patients and controls did not survive correction for multiple testing. In this coordinated and harmonized multisite DTI study, we showed subtle, but widespread differences in WM microstructure in adult MDD, which may suggest structural disconnectivity in MDD