Neurostructural subgroup in 4291 individuals with schizophrenia identified using the subtype and stage inference algorithm.

Machine learning can be used to define subtypes of psychiatric conditions based on shared biological foundations of mental disorders. Here we analyzed cross-sectional brain images from 4,222 individuals with schizophrenia and 7038 healthy subjects pooled across 41 international cohorts from the ENIGMA, non-ENIGMA cohorts and public datasets. Using the Subtype and Stage Inference (SuStaIn) algorithm, we identify two distinct neurostructural subgroups by mapping the spatial and temporal 'trajectory' of gray matter change in schizophrenia. Subgroup 1 was characterized by an early cortical-predominant loss with enlarged striatum, whereas subgroup 2 displayed an early subcortical-predominant loss in the hippocampus, striatum and other subcortical regions. We confirmed the reproducibility of the two neurostructural subtypes across various sample sites, including Europe, North America and East Asia. This imaging-based taxonomy holds the potential to identify individuals with shared neurobiological attributes, thereby suggesting the viability of redefining existing disorder constructs based on biological factors.

Admixture mapping of cognitive function in diverse Hispanic and Latino adults: Results from the Hispanic Community Health Study/Study of Latinos.

INTRODUCTION: We conducted admixture mapping and fine-mapping analyses to identify ancestry-of-origin loci influencing cognitive abilities. METHODS: We estimated the association of local ancestry intervals across the genome with five neurocognitive measures in 7140 diverse Hispanic and Latino adults (mean age 55 years). We prioritized genetic variants in associated loci and tested them for replication in four independent cohorts. RESULTS: We identified nine local ancestry-associated regions for the five neurocognitive measures. There was strong biological support for the observed associations to cognitive function at all loci and there was statistical evidence of independent replication at 4q12, 9p22.1, and 13q12.13. DISCUSSION: Our study identified multiple novel loci harboring genes implicated in cognitive functioning and dementia, and uncovered ancestry-relevant genetic variants. It adds to our understanding of the genetic architecture of cognitive function in Hispanic and Latino adults and demonstrates the power of admixture mapping to discover unique haplotypes influencing cognitive function, complementing genome-wide association studies. HIGHLIGHTS: We identified nine ancestry-of-origin chromosomal regions associated with five neurocognitive traits. In each associated region, we identified single nucleotide polymorphisms (SNPs) that explained, at least in part, the admixture signal and were tested for replication in independent samples of Black, non-Hispanic White, and Hispanic/Latino adults with the same or similar neurocognitive tests. Statistical evidence of independent replication of the prioritized SNPs was observed for three of the nine associations, at chr4q12, chr9p22.1, and chr13q12.13. At all loci, there was strong biological support for the observed associations to cognitive function and dementia, prioritizing genes such as KIT, implicated in autophagic clearance of neurotoxic proteins and on mast cell and microglial-mediated inflammation; SLC24A2, implicated in synaptic plasticity associated with learning and memory; and MTMR6, implicated in phosphoinositide lipids metabolism.

Brain-age prediction: Systematic evaluation of site effects, and sample age range and size.

Structural neuroimaging data have been used to compute an estimate of the biological age of the brain (brain-age) which has been associated with other biologically and behaviorally meaningful measures of brain development and aging. The ongoing research interest in brain-age has highlighted the need for robust and publicly available brain-age models pre-trained on data from large samples of healthy individuals. To address this need we have previously released a developmental brain-age model. Here we expand this work to develop, empirically validate, and disseminate a pre-trained brain-age model to cover most of the human lifespan. To achieve this, we selected the best-performing model after systematically examining the impact of seven site harmonization strategies, age range, and sample size on brain-age prediction in a discovery sample of brain morphometric measures from 35,683 healthy individuals (age range: 5-90 years; 53.59% female). The pre-trained models were tested for cross-dataset generalizability in an independent sample comprising 2101 healthy individuals (age range: 8-80 years; 55.35% female) and for longitudinal consistency in a further sample comprising 377 healthy individuals (age range: 9-25 years; 49.87% female). This empirical examination yielded the following findings: (1) the accuracy of age prediction from morphometry data was higher when no site harmonization was applied; (2) dividing the discovery sample into two age-bins (5-40 and 40-90 years) provided a better balance between model accuracy and explained age variance than other alternatives; (3) model accuracy for brain-age prediction plateaued at a sample size exceeding 1600 participants. These findings have been incorporated into CentileBrain (https://centilebrain.org/#/brainAGE2), an open-science, web-based platform for individualized neuroimaging metrics.

Role of Neurocellular Endoplasmic Reticulum Stress Response in Alzheimer's Disease and Related Dementias Risk.

Currently, more than 55 million people around the world suffer from dementia, and Alzheimer's Disease and Related Dementias (ADRD) accounts for nearly 60-70% of all those cases. The spread of Alzheimer's Disease (AD) pathology and progressive neurodegeneration in the hippocampus and cerebral cortex is strongly correlated with cognitive decline in AD patients; however, the molecular underpinning of ADRD's causality is still unclear. Studies of postmortem AD brains and animal models of AD suggest that elevated endoplasmic reticulum (ER) stress may have a role in ADRD pathology through altered neurocellular homeostasis in brain regions associated with learning and memory. To study the ER stress-associated neurocellular response and its effects on neurocellular homeostasis and neurogenesis, we modeled an ER stress challenge using thapsigargin (TG), a specific inhibitor of sarco/endoplasmic reticulum Ca2+ ATPase (SERCA), in the induced pluripotent stem cell (iPSC)-derived neural stem cells (NSCs) of two individuals from our Mexican American Family Study (MAFS). High-content screening and transcriptomic analysis of the control and ER stress-challenged NSCs showed that the NSCs' ER stress response resulted in a significant decline in NSC self-renewal and an increase in apoptosis and cellular oxidative stress. A total of 2300 genes were significantly (moderated t statistics FDR-corrected p-value ≤ 0.05 and fold change absolute ≥ 2.0) differentially expressed (DE). The pathway enrichment and gene network analysis of DE genes suggests that all three unfolded protein response (UPR) pathways, protein kinase RNA-like ER kinase (PERK), activating transcription factor-6 (ATF-6), and inositol-requiring enzyme-1 (IRE1), were significantly activated and cooperatively regulated the NSCs' transcriptional response to ER stress. Our results show that IRE1/X-box binding protein 1 (XBP1) mediated transcriptional regulation of the E2F transcription factor 1 (E2F1) gene, and its downstream targets have a dominant role in inducing G1/S-phase cell cycle arrest in ER stress-challenged NSCs. The ER stress-challenged NSCs also showed the activation of C/EBP homologous protein (CHOP)-mediated apoptosis and the dysregulation of synaptic plasticity and neurotransmitter homeostasis-associated genes. Overall, our results suggest that the ER stress-associated attenuation of NSC self-renewal, increased apoptosis, and dysregulated synaptic plasticity and neurotransmitter homeostasis plausibly play a role in the causation of ADRD.

When virtual reality becomes psychoneuroendocrine reality: A stress(or) review.

This review article was awarded the Dirk Hellhammer award from ISPNE in 2023. It explores the dynamic relationship between stressors and stress from a historical view as well as a vision towards the future of stress research via virtual reality (VR). We introduce the concept of a "syncytium," a permeable boundary that blurs the distinction between stress and stressor, in order to understand why the field of stress biology continues to inadequately measure stress alone as a proxy for the force of external stressors. Using Virtual Reality (VR) as an illustrative example to explicate the black box of stressors, we examine the distinction between 'immersion' and 'presence' as analogous terms for stressor and stress, respectively. We argue that the conventional psychological approaches to stress measurement and appraisal theory unfortunately fall short in quantifying the force of the stressor, leading to reverse causality fallacies. Further, the concept of affordances is introduced as an ecological or holistic tool to measure and design a stressor's force, bridging the gap between the external environment and an individual's physiological response to stress. Affordances also serve to ameliorate shortcomings in stress appraisal by integrating ecological interdependencies. By combining VR and psychobiological measures, this paper aims to unravel the complexity of the stressor-stress syncytium, highlighting the necessity of assessing both the internal and external facets to gain a holistic understanding of stress physiology and shift away from reverse causality reasoning. We find that the utility of VR extends beyond presence to include affordance-based measures of immersion, which can effectively model stressor force. Future research should prioritize the development of tools that can measure both immersion and presence, thereby providing a more comprehensive understanding of how external stressors interact with individual psychological states.

Bridging the gap: improving correspondence between low-field and high-field magnetic resonance images in young people.

Background: Portable low-field-strength magnetic resonance imaging (MRI) systems represent a promising alternative to traditional high-field-strength systems with the potential to make MR technology available at scale in low-resource settings. However, lower image quality and resolution may limit the research and clinical potential of these devices. We tested two super-resolution methods to enhance image quality in a low-field MR system and compared their correspondence with images acquired from a high-field system in a sample of young people. Methods: T1- and T2-weighted structural MR images were obtained from a low-field (64mT) Hyperfine and high-field (3T) Siemens system in N = 70 individuals (mean age = 20.39 years, range 9-26 years). We tested two super-resolution approaches to improve image correspondence between images acquired at high- and low-field: (1) processing via a convolutional neural network ('SynthSR'), and (2) multi-orientation image averaging. We extracted brain region volumes, cortical thickness, and cortical surface area estimates. We used Pearson correlations to test the correspondence between these measures, and Steiger Z tests to compare the difference in correspondence between standard imaging and super-resolution approaches. Results: Single pairs of T1- and T2-weighted images acquired at low field showed high correspondence to high-field-strength images for estimates of total intracranial volume, surface area cortical volume, subcortical volume, and total brain volume (r range = 0.60-0.88). Correspondence was lower for cerebral white matter volume (r = 0.32, p = 0.007, q = 0.009) and non-significant for mean cortical thickness (r = -0.05, p = 0.664, q = 0.664). Processing images with SynthSR yielded significant improvements in correspondence for total brain volume, white matter volume, total surface area, subcortical volume, cortical volume, and total intracranial volume (r range = 0.85-0.97), with the exception of global mean cortical thickness (r = 0.14). An alternative multi-orientation image averaging approach improved correspondence for cerebral white matter and total brain volume. Processing with SynthSR also significantly improved correspondence across widespread regions for estimates of cortical volume, surface area and subcortical volume, as well as within isolated prefrontal and temporal regions for estimates of cortical thickness. Conclusion: Applying super-resolution approaches to low-field imaging improves regional brain volume and surface area accuracy in young people. Finer-scale brain measurements, such as cortical thickness, remain challenging with the limited resolution of low-field systems.

Using rare genetic mutations to revisit structural brain asymmetry.

Asymmetry between the left and right hemisphere is a key feature of brain organization. Hemispheric functional specialization underlies some of the most advanced human-defining cognitive operations, such as articulated language, perspective taking, or rapid detection of facial cues. Yet, genetic investigations into brain asymmetry have mostly relied on common variants, which typically exert small effects on brain-related phenotypes. Here, we leverage rare genomic deletions and duplications to study how genetic alterations reverberate in human brain and behavior. We designed a pattern-learning approach to dissect the impact of eight high-effect-size copy number variations (CNVs) on brain asymmetry in a multi-site cohort of 552 CNV carriers and 290 non-carriers. Isolated multivariate brain asymmetry patterns spotlighted regions typically thought to subserve lateralized functions, including language, hearing, as well as visual, face and word recognition. Planum temporale asymmetry emerged as especially susceptible to deletions and duplications of specific gene sets. Targeted analysis of common variants through genome-wide association study (GWAS) consolidated partly diverging genetic influences on the right versus left planum temporale structure. In conclusion, our gene-brain-behavior data fusion highlights the consequences of genetically controlled brain lateralization on uniquely human cognitive capacities.

Lifespan reference curves for harmonizing multi-site regional brain white matter metrics from diffusion MRI.

Age-related white matter (WM) microstructure maturation and decline occur throughout the human lifespan, complementing the process of gray matter development and degeneration. Here, we create normative lifespan reference curves for global and regional WM microstructure by harmonizing diffusion MRI (dMRI)-derived data from ten public datasets (N = 40,898 subjects; age: 3-95 years; 47.6% male). We tested three harmonization methods on regional diffusion tensor imaging (DTI) based fractional anisotropy (FA), a metric of WM microstructure, extracted using the ENIGMA-DTI pipeline. ComBat-GAM harmonization provided multi-study trajectories most consistent with known WM maturation peaks. Lifespan FA reference curves were validated with test-retest data and used to assess the effect of the ApoE4 risk factor for dementia in WM across the lifespan. We found significant associations between ApoE4 and FA in WM regions associated with neurodegenerative disease even in healthy individuals across the lifespan, with regional age-by-genotype interactions. Our lifespan reference curves and tools to harmonize new dMRI data to the curves are publicly available as eHarmonize (https://github.com/ahzhu/eharmonize).

Normative modelling of brain morphometry across the lifespan with CentileBrain: algorithm benchmarking and model optimisation.

The value of normative models in research and clinical practice relies on their robustness and a systematic comparison of different modelling algorithms and parameters; however, this has not been done to date. We aimed to identify the optimal approach for normative modelling of brain morphometric data through systematic empirical benchmarking, by quantifying the accuracy of different algorithms and identifying parameters that optimised model performance. We developed this framework with regional morphometric data from 37 407 healthy individuals (53% female and 47% male; aged 3-90 years) from 87 datasets from Europe, Australia, the USA, South Africa, and east Asia 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 fractional polynomial regression (MFPR) emerged as the preferred algorithm, optimised with non-linear 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 3000 study participants. This model can inform about the biological and behavioural implications of deviations from typical age-related neuroanatomical changes and support future study designs. The model and scripts described here are freely available through CentileBrain.

Connectome architecture shapes large-scale cortical alterations in schizophrenia: a worldwide ENIGMA study.

Schizophrenia is a prototypical network disorder with widespread brain-morphological alterations, yet it remains unclear whether these distributed alterations robustly reflect the underlying network layout. We tested whether large-scale structural alterations in schizophrenia relate to normative structural and functional connectome architecture, and systematically evaluated robustness and generalizability of these network-level alterations. Leveraging anatomical MRI scans from 2439 adults with schizophrenia and 2867 healthy controls from 26 ENIGMA sites and normative data from the Human Connectome Project (n = 207), we evaluated structural alterations of schizophrenia against two network susceptibility models: (i) hub vulnerability, which examines associations between regional network centrality and magnitude of disease-related alterations; (ii) epicenter mapping, which identifies regions whose typical connectivity profile most closely resembles the disease-related morphological alterations. To assess generalizability and specificity, we contextualized the influence of site, disease stages, and individual clinical factors and compared network associations of schizophrenia with that found in affective disorders. Our findings show schizophrenia-related cortical thinning is spatially associated with functional and structural hubs, suggesting that highly interconnected regions are more vulnerable to morphological alterations. Predominantly temporo-paralimbic and frontal regions emerged as epicenters with connectivity profiles linked to schizophrenia's alteration patterns. Findings were robust across sites, disease stages, and related to individual symptoms. Moreover, transdiagnostic comparisons revealed overlapping epicenters in schizophrenia and bipolar, but not major depressive disorder, suggestive of a pathophysiological continuity within the schizophrenia-bipolar-spectrum. In sum, cortical alterations over the course of schizophrenia robustly follow brain network architecture, emphasizing marked hub susceptibility and temporo-frontal epicenters at both the level of the group and the individual. Subtle variations of epicenters across disease stages suggest interacting pathological processes, while associations with patient-specific symptoms support additional inter-individual variability of hub vulnerability and epicenters in schizophrenia. Our work outlines potential pathways to better understand macroscale structural alterations, and inter- individual variability in schizophrenia.

Structural brain abnormalities and aggressive behaviour in schizophrenia: Mega-analysis of data from 2095 patients and 2861 healthy controls via the ENIGMA consortium.

Background: Schizophrenia is associated with an increased risk of aggressive behaviour, which may partly be explained by illness-related changes in brain structure. However, previous studies have been limited by group-level analyses, small and selective samples of inpatients and long time lags between exposure and outcome. Methods: This cross-sectional study pooled data from 20 sites participating in the international ENIGMA-Schizophrenia Working Group. Sites acquired T1-weighted and diffusion-weighted magnetic resonance imaging scans in a total of 2095 patients with schizophrenia and 2861 healthy controls. Measures of grey matter volume and white matter microstructural integrity were extracted from the scans using harmonised protocols. For each measure, normative modelling was used to calculate how much patients deviated (in z-scores) from healthy controls at the individual level. Ordinal regression models were used to estimate the associations of these deviations with concurrent aggressive behaviour (as odds ratios [ORs] with 99% confidence intervals [CIs]). Mediation analyses were performed for positive symptoms (i.e., delusions, hallucinations and disorganised thinking), impulse control and illness insight. Aggression and potential mediators were assessed with the Positive and Negative Syndrome Scale, Scale for the Assessment of Positive Symptoms or Brief Psychiatric Rating Scale. Results: Aggressive behaviour was significantly associated with reductions in total cortical volume (OR [99% CI] = 0.88 [0.78, 0.98], p = .003) and global white matter integrity (OR [99% CI] = 0.72 [0.59, 0.88], p = 3.50 × 10-5) and additional reductions in dorsolateral prefrontal cortex volume (OR [99% CI] = 0.85 [0.74, 0.97], p =.002), inferior parietal lobule volume (OR [99% CI] = 0.76 [0.66, 0.87], p = 2.20 × 10-7) and internal capsule integrity (OR [99% CI] = 0.76 [0.63, 0.92], p = 2.90 × 10-4). Except for inferior parietal lobule volume, these associations were largely mediated by increased severity of positive symptoms and reduced impulse control. Conclusions: This study provides evidence that the co-occurrence of positive symptoms, poor impulse control and aggressive behaviour in schizophrenia has a neurobiological basis, which may inform the development of therapeutic interventions.

Returning Individual Research Results from Digital Phenotyping in Psychiatry.

Psychiatry is rapidly adopting digital phenotyping and artificial intelligence/machine learning tools to study mental illness based on tracking participants' locations, online activity, phone and text message usage, heart rate, sleep, physical activity, and more. Existing ethical frameworks for return of individual research results (IRRs) are inadequate to guide researchers for when, if, and how to return this unprecedented number of potentially sensitive results about each participant's real-world behavior. To address this gap, we convened an interdisciplinary expert working group, supported by a National Institute of Mental Health grant. Building on established guidelines and the emerging norm of returning results in participant-centered research, we present a novel framework specific to the ethical, legal, and social implications of returning IRRs in digital phenotyping research. Our framework offers researchers, clinicians, and Institutional Review Boards (IRBs) urgently needed guidance, and the principles developed here in the context of psychiatry will be readily adaptable to other therapeutic areas.

Riemannian frameworks for the harmonization of resting-state functional MRI scans.

Magnetic Resonance Imaging provides unprecedented images of the brain. Unfortunately, scanners and acquisition protocols can significantly impact MRI scans. The development of statistical methods able to reduce this variability without altering the relevant information in the scans, often coined harmonization methods, has been the topic of an increasing research effort supported by the recent growth of publicly available neuroimaging data sets and new possibilities for combining them to achieve greater statistical power. In this work, we focus on the challenges specifically raised by the harmonization of resting-state functional MRI scans. We propose to harmonize resting-state fMRI scans by reducing the impact of covariates such as scanner differences and scanning protocols on their associated functional connectomes and then propagating the changes back to the rs-fMRI time series. We use Riemannian geometric frameworks to preserve the mathematical properties of functional connectomes during their harmonization, and we demonstrate how state-of-the-art harmonization methods can be embedded within these frameworks to reduce covariates effects while preserving the relevant clinical information associated with aging or brain disorders. During our experiments, a large set of synthetic data was generated and processed to compare eighty variants of the proposed approach. The framework achieving the best harmonization was then applied to three low-dimensional data sets made of 712 sets of fMRI time series provided by the ABIDE consortium and two high-dimensional data sets obtained by processing 1527 rs-fMRI scans provided by the Human Connectome Project, the Framingham Heart Study and the Genetics of Brain Structure and Function study. These experiments established that our new framework could successfully harmonize low-dimensional connectomes and voxelwise functional time series and confirmed the need for preserving connectomes properties during their harmonization.

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.

Normative Modeling of Brain Morphometry Across the Lifespan Using CentileBrain: Algorithm Benchmarking and Model Optimization.

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

Estimating multimodal brain variability in schizophrenia spectrum disorders: A worldwide ENIGMA study.

Objective: Schizophrenia is a multifaceted disorder associated with structural brain heterogeneity. Despite its relevance for identifying illness subtypes and informative biomarkers, structural brain heterogeneity in schizophrenia remains incompletely understood. Therefore, the objective of this study was to provide a comprehensive insight into the structural brain heterogeneity associated with schizophrenia. Methods: This meta- and mega-analysis investigated the variability of multimodal structural brain measures of white and gray matter in individuals with schizophrenia versus healthy controls. Using the ENIGMA dataset of MRI-based brain measures from 22 international sites with up to 6139 individuals for a given brain measure, we examined variability in cortical thickness, surface area, folding index, subcortical volume and fractional anisotropy. Results: We found that individuals with schizophrenia are distinguished by higher heterogeneity in the frontotemporal network with regard to multimodal structural measures. Moreover, individuals with schizophrenia showed higher homogeneity of the folding index, especially in the left parahippocampal region. Conclusions: Higher multimodal heterogeneity in frontotemporal regions potentially implies different subtypes of schizophrenia that converge on impaired frontotemporal interaction as a core feature of the disorder. Conversely, more homogeneous folding patterns in the left parahippocampal region might signify a consistent characteristic of schizophrenia shared across subtypes. These findings underscore the importance of structural brain variability in advancing our neurobiological understanding of schizophrenia, and aid in identifying illness subtypes as well as informative biomarkers.

Two neurostructural subtypes: results of machine learning on brain images from 4,291 individuals with schizophrenia.

Machine learning can be used to define subtypes of psychiatric conditions based on shared clinical and biological foundations, presenting a crucial step toward establishing biologically based subtypes of mental disorders. With the goal of identifying subtypes of disease progression in schizophrenia, here we analyzed cross-sectional brain structural magnetic resonance imaging (MRI) data from 4,291 individuals with schizophrenia (1,709 females, age=32.5 years±11.9) and 7,078 healthy controls (3,461 females, age=33.0 years±12.7) pooled across 41 international cohorts from the ENIGMA Schizophrenia Working Group, non-ENIGMA cohorts and public datasets. Using a machine learning approach known as Subtype and Stage Inference (SuStaIn), we implemented a brain imaging-driven classification that identifies two distinct neurostructural subgroups by mapping the spatial and temporal trajectory of gray matter (GM) loss in schizophrenia. Subgroup 1 (n=2,622) was characterized by an early cortical-predominant loss (ECL) with enlarged striatum, whereas subgroup 2 (n=1,600) displayed an early subcortical-predominant loss (ESL) in the hippocampus, amygdala, thalamus, brain stem and striatum. These reconstructed trajectories suggest that the GM volume reduction originates in the Broca's area/adjacent fronto-insular cortex for ECL and in the hippocampus/adjacent medial temporal structures for ESL. With longer disease duration, the ECL subtype exhibited a gradual worsening of negative symptoms and depression/anxiety, and less of a decline in positive symptoms. We confirmed the reproducibility of these imaging-based subtypes across various sample sites, independent of macroeconomic and ethnic factors that differed across these geographic locations, which include Europe, North America and East Asia. These findings underscore the presence of distinct pathobiological foundations underlying schizophrenia. This new imaging-based taxonomy holds the potential to identify a more homogeneous sub-population of individuals with shared neurobiological attributes, thereby suggesting the viability of redefining existing disorder constructs based on biological factors.

Multimodal Neuroimaging Summary Scores as Neurobiological Markers of Psychosis.

BACKGROUND AND HYPOTHESIS: Structural brain alterations are well-established features of schizophrenia but they do not effectively predict disease/disease risk. Similar to polygenic risk scores in genetics, we integrated multifactorial aspects of brain structure into a summary "Neuroscore" and examined its potential as a marker of disease. STUDY DESIGN: We extracted measures from T1-weighted scans and diffusion tensor imaging (DTI) models from three studies with schizophrenia and healthy individuals. We calculated individual-level summary scores (Neuroscores) for T1-weighted and DTI measures and a combined score (Multimodal Neuroscore-MM). We assessed each score's ability to differentiate schizophrenia cases from controls and its relationship to clinical symptomatology, intelligence quotient (IQ), and medication dosage. We assessed Neuroscore specificity by performing all analyses in a more inclusive psychosis sample and by using scores generated from MDD effect sizes. STUDY RESULTS: All Neuroscores significantly differentiated schizophrenia cases from controls (T1 d = 0.56, DTI d = 0.29, MM d = 0.64) to a greater degree than individual brain regions. Higher Neuroscores (ie, increased liability) were associated with lower IQ (T1 ß = -0.26, DTI ß = -0.15, MM ß = -0.30). Higher T1-weighted Neuroscores were associated with higher positive and negative symptom severity (Positive ß = 0.21, Negative ß = 0.16); Higher Multimodal Neuroscores were associated with higher positive symptom severity (ß = 0.30). SZ Neuroscores outperformed MDD Neuroscores in predicting IQ (T1: z = 3.5, q = 0.0007; MM: z = 1.8, q = 0.05). CONCLUSIONS: Neuroscores are a step toward leveraging widespread structural brain alterations in psychosis to identify robust neurobiological markers of disease. Future studies will assess ways to improve neuroscore calculation, including developing the optimal methods to calculate neuroscores and considering disorder overlap.

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.