Mapping the Heterogeneous Phenotype of Schizophrenia and Bipolar Disorder Using Normative Models

Lifespan development of thalamic nuclei and characterizing thalamic nuclei abnormalities in schizophrenia using normative modeling

Thalamic abnormalities have been repeatedly implicated in the pathophysiology of schizophrenia and other neurodevelopmental disorders. Uncovering the etiology of thalamic abnormalities and how they may contribute to illness phenotypes faces at least two obstacles. First, the typical developmental trajectories of thalamic nuclei and their association with cognition across the lifespan are largely unknown. Second, modest effect sizes indicate marked individual differences and pose a significant challenge to personalized medicine. To address these knowledge gaps, we characterized the development of thalamic nuclei volumes using normative models generated from the Human Connectome Project Lifespan datasets (5-100+ years), then applied them to an independent clinical cohort to determine the frequency of thalamic volume deviations in people with schizophrenia (17-61 years). Normative models revealed diverse non-linear age effects across the lifespan. Association nuclei exhibited negative age effects during youth but stabilized in adulthood until turning negative again with older age. Sensorimotor nuclei volumes remained relatively stable through youth and adulthood until also turning negative with older age. Up to 18% of individuals with schizophrenia exhibited abnormally small (i.e., below the 5th centile) mediodorsal and pulvinar volumes, and the degree of deviation, but not raw volumes, correlated with the severity of cognitive impairment. While case-control differences are robust, only a minority of patients demonstrate unusually small thalamic nuclei volumes. Normative modeling enables the identification of these individuals, which is a necessary step toward precision medicine.

The future of disentangling the heterogeneity of autism with neuroimaging studies

Autism spectrum disorder (ASD) is a lifelong neurodevelopmental condition. Over the past decade, a considerable number of approaches have been developed to identify potential neuroimaging-based biomarkers of ASD and have uncovered specific neural mechanisms underlying behaviors associated with ASD. However, the substantial heterogeneity among those diagnosed with ASD hinders the development of biomarkers. Disentangling the heterogeneity of ASD is pivotal to improve quality of life for individuals with ASD by facilitating early diagnosis and individualized interventions for those who need support. In this Review, we discuss recent advances in neuroimaging that have facilitated the characterization of the heterogeneity of this condition from three frameworks: neurosubtyping, dimensional models, and normative models. In addition, we discuss the challenges, possible solutions, and clinical utility of these three frameworks. We argue that several factors need to be considered when parsing heterogeneity using neuroimaging, including co-occurring conditions, neurodevelopment, heredity and environment, and multi-site and multi-modality data. We close with a discussion of future directions for achieving a better understanding of the neural mechanisms underlying neurodevelopmental heterogeneity and the future of precision medicine in ASD.

Individualized functional connectome abnormalities obtained using two normative model unveil neurophysiological subtypes of obsessive compulsive disorder

The high heterogeneity observed among patients with obsessive-compulsive disorder (OCD) underscores the need to identify neurophysiological OCD subtypes to facilitate personalized diagnosis and treatment. In this study, our aim was to identify potential OCD subtypes based on individualized functional connectome abnormalities. We recruited a total of 99 patients with OCD and 104 healthy controls (HCs) matched for demographic characteristics. Individualized functional connectome abnormalities were obtained using normative models of functional connectivity strength (FCS) and used as features to unveil OCD subtypes. Sensitivity analyses were conducted to assess the reproducibility and robustness of the clustering outcomes. Patients exhibited significant intersubject heterogeneity in individualized functional connectome abnormalities. Two subtypes with distinct patterns of FCS abnormalities relative to HCs were identified. Subtype 1 patients primarily exhibited significantly decreased FCS in regions including the frontal gyrus, insula, hippocampus, and precentral/postcentral gyrus, whereas subtype 2 patients demonstrated increased FCS in widespread brain regions. When all patients were combined, no significant differences were observed. Additionally, the identified subtypes showed significant differences in age of onset. Furthermore, sensitivity analyses confirmed the reproducibility of our subtyping results. In conclusion, the identified OCD subtypes shed new light on the taxonomy and neurophysiological heterogeneity of OCD.

Individualized gray matter morphological abnormalities uncover two robust transdiagnostic biotypes

Psychiatric disorders exhibit a shared neuropathology, yet the diverse presentations among patients necessitate the identification of transdiagnostic subtypes to enhance diagnostic and treatment strategies. This study aims to unveil potential transdiagnostic subtypes based on personalized gray matter morphological abnormalities. A total of 496 patients with psychiatric disorders and 255 healthy controls (HCs) from three distinct datasets (one for discovery and two for validation) were enrolled. Individualized gray matter morphological abnormalities were determined using normative modeling to identify transdiagnostic subtypes. In the discovery dataset, two transdiagnostic subtypes with contrasting patterns of structural abnormalities compared to HCs were identified. Reproducibility and generalizability analyses demonstrated that these subtypes could be generalized to new patients and even to new disorders in the validation datasets. These subtypes were characterized by distinct disease epicenters. The gray matter abnormal pattern in subtype 1 was mainly linked to excitatory receptors, whereas subtype 2 showed a predominant association with inhibitory receptors. Furthermore, we observed that the gray matter abnormal pattern in subtype 2 was correlated with transcriptional profiles of inflammation-related genes, while subtype 1 did not show this association. Our findings reveal two robust transdiagnostic biotypes, offering novel insights into psychiatric nosology.

Multiscale heterogeneity of white matter morphometry in psychiatric disorders

Background

Inter-individual variability in neurobiological and clinical characteristics in mental illness is often overlooked by classical group-mean case-control studies. Studies using normative modelling to infer person-specific deviations of grey matter volume have indicated that group means are not representative of most individuals. The extent to which this variability is present in white matter morphometry, which is integral to brain function, remains unclear.

Methods

We applied Warped Bayesian Linear Regression normative models to T1-weighted magnetic resonance imaging data and mapped inter-individual variability in person-specific white matter volume deviations in 1,294 cases (58% male) diagnosed with one of six disorders (attention-deficit/hyperactivity, autism, bipolar, major depressive, obsessive-compulsive and schizophrenia) and 1,465 matched controls (54% male) recruited across 25 scan sites. We developed a framework to characterize deviation heterogeneity at multiple spatial scales, from individual voxels, through inter-regional connections, specific brain regions, and spatially extended brain networks.

Results

The specific locations of white matter volume deviations were highly heterogeneous across participants, affecting the same voxel in fewer than 8% of individuals with the same diagnosis. For autism and schizophrenia, negative deviations (i.e., areas where volume is lower than normative expectations) aggregated into common tracts, regions and large-scale networks in up to 35% of individuals.

Conclusions

The prevalence of white matter volume deviations was lower than previously observed in grey matter, and the specific location of these deviations was highly heterogeneous when considering voxel-wise spatial resolution. Evidence of aggregation within common pathways and networks was apparent in schizophrenia and autism but not other disorders.

Individualized Texture Similarity Network in Schizophrenia

BACKGROUND

Structural covariance network disruption has been considered an important pathophysiological indicator for schizophrenia. Here, we introduced a novel individualized structural covariance network measure, referred to as a texture similarity network (TSN), and hypothesized that the TSN could reliably reveal unique intersubject heterogeneity and complex dysconnectivity patterns in schizophrenia.

METHODS

The TSN was constructed by measuring the covariance of 180 three-dimensional voxelwise gray-level co-occurrence matrix feature maps between brain areas in each participant. We first tested the validity and reproducibility of the TSN in characterizing the intersubject variability in 2 longitudinal test-retest healthy cohorts. The TSN was further applied to elucidate intersubject variability and dysconnectivity patterns in 10 schizophrenia case-control datasets (609 schizophrenia cases vs. 579 controls) as well as in a first-episode depression dataset (69 patients with depression vs. 69 control participants).

RESULTS

The test-retest analysis demonstrated higher TSN intersubject than intrasubject variability. Moreover, the TSN reliably revealed higher intersubject variability in both chronic and first-episode schizophrenia, but not in depression. The TSN also reproducibly detected coexistent increased and decreased TSN strength in widespread brain areas, increased global small-worldness, and the coexistence of both structural hyposynchronization in the central networks and hypersynchronization in peripheral networks in patients with schizophrenia but not in patients with depression. Finally, aberrant intersubject variability and covariance strength patterns revealed by the TSN showed a missing or weak correlation with other individualized structural covariance network measures, functional connectivity, and regional volume changes.

CONCLUSIONS

These findings support the reliability of a TSN in revealing unique structural heterogeneity and complex dysconnectivity in patients with schizophrenia.

Exploring the potential of representation and transfer learning for anatomical neuroimaging: Application to psychiatry

The perspective of personalized medicine for brain disorders requires efficient learning models for anatomical neuroimaging-based prediction of clinical conditions. There is now a consensus on the benefit of deep learning (DL) in addressing many medical imaging tasks, such as image segmentation. However, for single-subject prediction problems, recent studies yielded contradictory results when comparing DL with Standard Machine Learning (SML) on top of classical feature extraction. Most existing comparative studies were limited in predicting phenotypes of little clinical interest, such as sex and age, and using a single dataset. Moreover, they conducted a limited analysis of the employed image pre-processing and feature selection strategies. This paper extensively compares DL and SML prediction capacity on five multi-site problems, including three increasingly complex clinical applications in psychiatry namely schizophrenia, bipolar disorder, and Autism Spectrum Disorder (ASD) diagnosis. To compensate for the relative scarcity of neuroimaging data on these clinical datasets, we also evaluate three pre-training strategies for transfer learning from brain imaging of the general healthy population: self-supervised learning, generative modeling and supervised learning with age. Overall, we find similar performance between randomly initialized DL and SML for the three clinical tasks and a similar scaling trend for sex prediction. This was replicated on an external dataset. We also show highly correlated discriminative brain regions between DL and linear ML models in all problems. Nonetheless, we demonstrate that self-supervised pre-training on large-scale healthy population imaging datasets (N≈10k), along with Deep Ensemble, allows DL to learn robust and transferable representations to smaller-scale clinical datasets (N≤1k). It largely outperforms SML on 2 out of 3 clinical tasks both in internal and external test sets. These findings suggest that the improvement of DL over SML in anatomical neuroimaging mainly comes from its capacity to learn meaningful and useful abstract representations of the brain anatomy, and it sheds light on the potential of transfer learning for personalized medicine in psychiatry.

Dissecting task-based fMRI activity using normative modelling: an application to the Emotional Face Matching Task

Functional neuroimaging has contributed substantially to understanding brain function but is dominated by group analyses that index only a fraction of the variation in these data. It is increasingly clear that parsing the underlying heterogeneity is crucial to understand individual differences and the impact of different task manipulations. We estimate large-scale (N = 7728) normative models of task-evoked activation during the Emotional Face Matching Task, which enables us to bind heterogeneous datasets to a common reference and dissect heterogeneity underlying group-level analyses. We apply this model to a heterogenous patient cohort, to map individual differences between patients with one or more mental health diagnoses relative to the reference cohort and determine multivariate associations with transdiagnostic symptom domains. For the face>shapes contrast, patients have a higher frequency of extreme deviations which are spatially heterogeneous. In contrast, normative models for faces>baseline have greater predictive value for individuals' transdiagnostic functioning. Taken together, we demonstrate that normative modelling of fMRI task-activation can be used to illustrate the influence of different task choices and map replicable individual differences, and we encourage its application to other neuroimaging tasks in future studies.

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.

A 10-Year Longitudinal Study of Brain Cortical Thickness in People with First-Episode Psychosis Using Normative Models

BACKGROUND

Clinical forecasting models have potential to optimize treatment and improve outcomes in psychosis, but predicting long-term outcomes is challenging and long-term follow-up data are scarce. In this 10-year longitudinal study, we aimed to characterize the temporal evolution of cortical correlates of psychosis and their associations with symptoms.

DESIGN

Structural magnetic resonance imaging (MRI) from people with first-episode psychosis and controls (n = 79 and 218) were obtained at enrollment, after 12 months (n = 67 and 197), and 10 years (n = 23 and 77), within the Thematically Organized Psychosis (TOP) study. Normative models for cortical thickness estimated on public MRI datasets (n = 42 983) were applied to TOP data to obtain deviation scores for each region and timepoint. Positive and Negative Syndrome Scale (PANSS) scores were acquired at each timepoint along with registry data. Linear mixed effects models assessed effects of diagnosis, time, and their interactions on cortical deviations plus associations with symptoms.

RESULTS

LMEs revealed conditional main effects of diagnosis and time × diagnosis interactions in a distributed cortical network, where negative deviations in patients attenuate over time. In patients, symptoms also attenuate over time. LMEs revealed effects of anterior cingulate on PANSS total, and insular and orbitofrontal regions on PANSS negative scores.

CONCLUSIONS

This long-term longitudinal study revealed a distributed pattern of cortical differences which attenuated over time together with a reduction in symptoms. These findings are not in line with a simple neurodegenerative account of schizophrenia, and deviations from normative models offer a promising avenue to develop biomarkers to track clinical trajectories over time.

Heterogeneity of morphometric similarity networks in health and schizophrenia

Introduction

Morphometric similarity is a recently developed neuroimaging phenotype of inter-regional connectivity by quantifying the similarity of a region to other regions based on multiple MRI parameters. Altered average morphometric similarity has been reported in psychotic disorders at the group level, with considerable heterogeneity across individuals. We used normative modeling to address cross-sectional and longitudinal inter-individual heterogeneity of morphometric similarity in health and schizophrenia.

Methods

Morphometric similarity for 62 cortical regions was obtained from baseline and follow-up T1-weighted scans of healthy individuals and patients with chronic schizophrenia. Cortical regions were classified into seven predefined brain functional networks. Using Bayesian Linear Regression and taking into account age, sex, image quality and scanner, we trained and validated normative models in healthy controls from eleven datasets (n = 4310). Individual deviations from the norm (z-scores) in morphometric similarity were computed for each participant for each network and region at both timepoints. A z-score ≧ than 1.96 was considered supra-normal and a z-score ≦ -1.96 infra-normal. As a longitudinal metric, we calculated the change over time of the total number of infra- or supra-normal regions per participant.

Results

At baseline, patients with schizophrenia had decreased morphometric similarity of the default mode network and increased morphometric similarity of the somatomotor network when compared with healthy controls. The percentage of patients with infra- or supra-normal values for any region at baseline and follow-up was low (<6%) and did not differ from healthy controls. Mean intra-group changes over time in the total number of infra- or supra-normal regions were small in schizophrenia and healthy control groups (<1) and there were no significant between-group differences.

Conclusions

In a case-control setting, a decrease of morphometric similarity within the default mode network may be a robust finding implicated in schizophrenia. However, normative modeling suggests that significant reductions and changes over time of regional morphometric similarity are evident only in a minority of patients.

Topological abnormalities of the morphometric similarity network of the cerebral cortex in schizophrenia

A morphometric similarity (MS) network can be constructed using multiple magnetic resonance imaging parameters of each cortical region. An MS network can be used to assess the similarity between cortical regions. Although MS networks can detect microstructural alterations and capture connections between histologically similar cortical areas, the influence of schizophrenia on the topological characteristics of MS networks remains unclear. We obtained T1- and diffusion-weighted images of 239 healthy controls and 190 individuals with schizophrenia to construct the MS network. Group comparisons of the mean MS of the cortical regions and subnetworks were performed. The strengths of the connections between the cortical regions and the global and nodal network indices were compared between the groups. Clinical associations with the network indices were tested using Spearman's rho. Compared with healthy controls, individuals with schizophrenia had significant group differences in the mean MS of several cortical regions and subnetworks. Individuals with schizophrenia had both superior and inferior strengths of connections between cortical regions compared with those of healthy controls. We observed regional abnormalities of the MS network in individuals with schizophrenia regarding lower centrality values of the pars opercularis, superior frontal, and superior temporal areas. Specific nodal network measures of the right pars opercularis and left superior temporal areas were associated with illness duration in individuals with schizophrenia. We identified regional abnormalities of the MS network in schizophrenia with the left superior temporal area possibly being a key region in topological organization and cortical connections.

Functional connectome through the human life span

The lifespan growth of the functional connectome remains unknown. Here, we assemble task-free functional and structural magnetic resonance imaging data from 33,250 individuals aged 32 postmenstrual weeks to 80 years from 132 global sites. We report critical inflection points in the nonlinear growth curves of the global mean and variance of the connectome, peaking in the late fourth and late third decades of life, respectively. After constructing a fine-grained, lifespan-wide suite of system-level brain atlases, we show distinct maturation timelines for functional segregation within different systems. Lifespan growth of regional connectivity is organized along a primary-to-association cortical axis. These connectome-based normative models reveal substantial individual heterogeneities in functional brain networks in patients with autism spectrum disorder, major depressive disorder, and Alzheimer's disease. These findings elucidate the lifespan evolution of the functional connectome and can serve as a normative reference for quantifying individual variation in development, aging, and neuropsychiatric disorders.

Delineating the Heterogeneity of Alzheimer's Disease and Mild Cognitive Impairment Using Normative Models of Dynamic Brain Functional Networks

BACKGROUND

Alzheimer's disease (AD), which has been identified as the most common type of dementia, presents considerable heterogeneity in its clinical manifestations. Early intervention at the stage of mild cognitive impairment (MCI) holds potential in AD prevention. However, characterizing the heterogeneity of neurobiological abnormalities and identifying MCI subtypes pose significant challenges.

METHODS

We constructed sex-specific normative age models of dynamic brain functional networks and mapped the deviations of the brain characteristics for individuals from multiple datasets, including 295 patients with AD, 441 patients with MCI, and 1160 normal control participants. Then, based on these individual deviation patterns, subtypes for both AD and MCI were identified using the clustering method, and their similarities and differences were comprehensively assessed.

RESULTS

Individuals with AD and MCI were clustered into 2 subtypes, and these subtypes exhibited significant differences in their intrinsic brain functional phenotypes and spatial atrophy patterns, as well as in disease progression and cognitive decline trajectories. The subtypes with positive deviations in AD and MCI shared similar deviation patterns, as did those with negative deviations. There was a potential transformation of MCI with negative deviation patterns into AD, and participants with MCI had a more severe cognitive decline rate.

CONCLUSIONS

In this study, we quantified neurophysiological heterogeneity by analyzing deviation patterns from the dynamic functional connectome normative model and identified disease subtypes of AD and MCI using a comprehensive resting-state functional magnetic resonance imaging multicenter dataset. The findings provide new insights for developing early prevention and personalized treatment strategies for AD.

Normative modelling of molecular-based functional circuits captures clinical heterogeneity transdiagnostically in psychiatric patients

Advanced methods such as REACT have allowed the integration of fMRI with the brain's receptor landscape, providing novel insights transcending the multiscale organisation of the brain. Similarly, normative modelling has allowed translational neuroscience to move beyond group-average differences and characterise deviations from health at an individual level. Here, we bring these methods together for the first time. We used REACT to create functional networks enriched with the main modulatory, inhibitory, and excitatory neurotransmitter systems and generated normative models of these networks to capture functional connectivity deviations in patients with schizophrenia, bipolar disorder (BPD), and ADHD. Substantial overlap was seen in symptomatology and deviations from normality across groups, but these could be mapped into a common space linking constellations of symptoms through to underlying neurobiology transdiagnostically. This work provides impetus for developing novel biomarkers that characterise molecular- and systems-level dysfunction at the individual level, facilitating the transition towards mechanistically targeted treatments.

Dynamic structure-function coupling across three major psychiatric disorders

BACKGROUND

Convergent evidence has suggested atypical relationships between brain structure and function in major psychiatric disorders, yet how the abnormal patterns coincide and/or differ across different disorders remains largely unknown. Here, we aim to investigate the common and/or unique dynamic structure-function coupling patterns across major depressive disorder (MDD), bipolar disorder (BD), and schizophrenia (SZ).

METHODS

We quantified the dynamic structure-function coupling in 452 patients with psychiatric disorders (MDD/BD/SZ = 166/168/118) and 205 unaffected controls at three distinct brain network levels, such as global, meso-, and local levels. We also correlated dynamic structure-function coupling with the topological features of functional networks to examine how the structure-function relationship facilitates brain information communication over time.

RESULTS

The dynamic structure-function coupling is preserved for the three disorders at the global network level. Similar abnormalities in the rich-club organization are found in two distinct functional configuration states at the meso-level and are associated with the disease severity of MDD, BD, and SZ. At the local level, shared and unique alterations are observed in the brain regions involving the visual, cognitive control, and default mode networks. In addition, the relationships between structure-function coupling and the topological features of functional networks are altered in a manner indicative of state specificity.

CONCLUSIONS

These findings suggest both transdiagnostic and illness-specific alterations in the dynamic structure-function relationship of large-scale brain networks across MDD, BD, and SZ, providing new insights and potential biomarkers into the neurodevelopmental basis underlying the behavioral and cognitive deficits observed in these disorders.

Spinal cord and brain atrophy patterns in neuromyelitis optica spectrum disorder and multiple sclerosis

BACKGROUND

Spinal cord and brain atrophy are common in neuromyelitis optica spectrum disorder (NMOSD) and relapsing-remitting multiple sclerosis (RRMS) but harbor distinct patterns accounting for disability and cognitive impairment.

METHODS

This study included 209 NMOSD and 304 RRMS patients and 436 healthy controls. Non-negative matrix factorization was used to parse differences in spinal cord and brain atrophy at subject level into distinct patterns based on structural MRI. The weights of patterns were obtained using a linear regression model and associated with Expanded Disability Status Scale (EDSS) and cognitive scores. Additionally, patients were divided into cognitive impairment (CI) and cognitive preservation (CP) groups.

RESULTS

Three patterns were observed in NMOSD: (1) Spinal Cord-Deep Grey Matter (SC-DGM) pattern was associated with high EDSS scores and decline of visuospatial memory function; (2) Frontal-Temporal pattern was associated with decline of language learning function; and (3) Cerebellum-Brainstem pattern had no observed association. Patients with CI had higher weights of SC-DGM pattern than CP group. Three patterns were observed in RRMS: (1) DGM pattern was associated with high EDSS scores, decreased information processing speed, and decreased language learning and visuospatial memory functions; (2) Frontal-Temporal pattern was associated with overall cognitive decline; and (3) Occipital pattern had no observed association. Patients with CI trended to have higher weights of DGM and Frontal-Temporal patterns than CP group.

CONCLUSION

This study estimated the heterogeneity of spinal cord and brain atrophy patterns in NMOSD and RRMS patients at individual level, and evaluated the clinical relevance of these patterns, which may contribute to stratifying participants for targeted therapy.

Aberrant individual structure covariance network in patients with mesial temporal lobe epilepsy

Objective

Mesial temporal lobe epilepsy (mTLE) is a complex neurological disorder that has been recognized as a widespread global network disorder. The group-level structural covariance network (SCN) could reveal the structural connectivity disruption of the mTLE but could not reflect the heterogeneity at the individual level.

Methods

This study adopted a recently proposed individual structural covariance network (IDSCN) method to clarify the alternated structural covariance connection mode in mTLE and to associate IDSCN features with the clinical manifestations and regional brain atrophy.

Results

We found significant IDSCN abnormalities in the ipsilesional hippocampus, ipsilesional precentral gyrus, bilateral caudate, and putamen in mTLE patients than in healthy controls. Moreover, the IDSCNs of these areas were positively correlated with the gray matter atrophy rate. Finally, we identified several connectivities with weak associations with disease duration, frequency, and surgery outcome.

Significance

Our research highlights the role of hippo-thalamic-basal-cortical circuits in the pathophysiologic process of disrupted whole-brain morphological covariance networks in mTLE, and builds a bridge between brain-wide covariance network changes and regional brain atrophy.

Schizophrenia, the very idea: On self-disorder, hyperreflexivity, and the diagnostic concept

The purpose of the present article is to consider schizophrenia-the very idea-from the perspective of phenomenological psychopathology, with special attention to the problematic nature of the diagnostic concept as well as to the prospect and challenges inherent in focusing on subjective experience. First, we address historical and philosophical topics relevant to the legitimacy of diagnostic categorization-in general and regarding "schizophrenia" in particular. William James's pragmatist approach to categorization is discussed. Then we offer a version of the well-known basic-self or ipseity-disturbance model (IDM) of schizophrenia, but in a significantly revised form (IDMrevised). The revised model better acknowledges the diverse and even seemingly contradictory nature of schizophrenic symptoms while, at the same time, interpreting these in a more unitary fashion via the key concept of hyperreflexivity-a form of exaggerated self-awareness that tends to undermine normal world-directedness and the stability of self-experience. Particular attention is paid to forms of exaggerated "self-presence" that are sometimes neglected yet imbue classically schizophrenic experiences involving subjectivism or quasi-solipsism and/or all-inclusive or ontological forms of paranoia. We focus on the distinctively paradoxical nature of schizophrenic symptomatology. In concluding we consider precursors in the work of Klaus Conrad, Kimura Bin and Henri Grivois. Finally we defend the concept of schizophrenia by considering its distinctive way of altering certain core aspects of the human condition itself.

Gray matter atrophy is constrained by normal structural brain network architecture in depression

BACKGROUND

There is growing evidence that gray matter atrophy is constrained by normal brain network (or connectome) architecture in neuropsychiatric disorders. However, whether this finding holds true in individuals with depression remains unknown. In this study, we aimed to investigate the association between gray matter atrophy and normal connectome architecture at individual level in depression.

METHODS

In this study, 297 patients with depression and 256 healthy controls (HCs) from two independent Chinese dataset were included: a discovery dataset (105 never-treated first-episode patients and matched 130 HCs) and a replication dataset (106 patients and matched 126 HCs). For each patient, individualized regional atrophy was assessed using normative model and brain regions whose structural connectome profiles in HCs most resembled the atrophy patterns were identified as putative epicenters using a backfoward stepwise regression analysis.

RESULTS

In general, the structural connectome architecture of the identified disease epicenters significantly explained 44% (±16%) variance of gray matter atrophy. While patients with depression demonstrated tremendous interindividual variations in the number and distribution of disease epicenters, several disease epicenters with higher participation coefficient than randomly selected regions, including the hippocampus, thalamus, and medial frontal gyrus were significantly shared by depression. Other brain regions with strong structural connections to the disease epicenters exhibited greater vulnerability. In addition, the association between connectome and gray matter atrophy uncovered two distinct subgroups with different ages of onset.

CONCLUSIONS

These results suggest that gray matter atrophy is constrained by structural brain connectome and elucidate the possible pathological progression in depression.

Signature of Altered Retinal Microstructures and Electrophysiology in Schizophrenia Spectrum Disorders Is Associated With Disease Severity and Polygenic Risk

BACKGROUND

Optical coherence tomography and electroretinography studies have revealed structural and functional retinal alterations in individuals with schizophrenia spectrum disorders (SSDs). However, it remains unclear which specific retinal layers are affected; how the retina, brain, and clinical symptomatology are connected; and how alterations of the visual system are related to genetic disease risk.

METHODS

Optical coherence tomography, electroretinography, and brain magnetic resonance imaging were applied to comprehensively investigate the visual system in a cohort of 103 patients with SSDs and 130 healthy control individuals. The sparse partial least squares algorithm was used to identify multivariate associations between clinical disease phenotype and biological alterations of the visual system. The association of the revealed patterns with individual polygenic disease risk for schizophrenia was explored in a post hoc analysis. In addition, covariate-adjusted case-control comparisons were performed for each individual optical coherence tomography and electroretinography parameter.

RESULTS

The sparse partial least squares analysis yielded a phenotype-eye-brain signature of SSDs in which greater disease severity, longer duration of illness, and impaired cognition were associated with electrophysiological alterations and microstructural thinning of most retinal layers. Higher individual loading onto this disease-relevant signature of the visual system was significantly associated with elevated polygenic risk for schizophrenia. In case-control comparisons, patients with SSDs had lower macular thickness, thinner retinal nerve fiber and inner plexiform layers, less negative a-wave amplitude, and lower b-wave amplitude.

CONCLUSIONS

This study demonstrates multimodal microstructural and electrophysiological retinal alterations in individuals with SSDs that are associated with disease severity and individual polygenic burden.

Large-Scale Neuroimaging of Mental Illness

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

Population-wide cerebellar growth models of children and adolescents

In the past, the cerebellum has been best known for its crucial role in motor function. However, increasingly more findings highlight the importance of cerebellar contributions in cognitive functions and neurodevelopment. Using a total of 7240 neuroimaging scans from 4862 individuals, we describe and provide detailed, openly available models of cerebellar development in childhood and adolescence (age range: 6-17 years), an important time period for brain development and onset of neuropsychiatric disorders. Next to a traditionally used anatomical parcellation of the cerebellum, we generated growth models based on a recently proposed functional parcellation. In both, we find an anterior-posterior growth gradient mirroring the age-related improvements of underlying behavior and function, which is analogous to cerebral maturation patterns and offers evidence for directly related cerebello-cortical developmental trajectories. Finally, we illustrate how the current approach can be used to detect cerebellar abnormalities in clinical samples.

Identifying two distinct neuroanatomical subtypes of first-episode depression using heterogeneity through discriminative analysis

BACKGROUND

Neurobiological heterogeneity in depression remains largely unknown, leading to inconsistent neuroimaging findings.

METHODS

Here, we adopted a novel proposed machine learning method ground on gray matter volumes (GMVs) to investigate neuroanatomical subtypes of first-episode treatment-naïve depression. GMVs were obtained from high-resolution T1-weighted images of 195 patients with first-episode, treatment-naïve depression and 78 matched healthy controls (HCs). Then we explored distinct subtypes of depression by employing heterogeneity through discriminative analysis (HYDRA) with regional GMVs as features.

RESULTS

Two prominently divergent subtypes of first-episode depression were identified, exhibiting opposite structural alterations compared with HCs but no different demographic features. Subtype 1 presented widespread increased GMVs mainly located in frontal, parietal, temporal cortex and partially located in limbic system. Subtype 2 presented widespread decreased GMVs mainly located in thalamus, cerebellum, limbic system and partially located in frontal, parietal, temporal cortex. Subtype 2 had smaller TIV and longer illness duration than Subtype 1. And TIV in Subtype 1 was positively correlated with age of onset while not in Subtype 2, probably implying the different potential neuropathological mechanisms.

LIMITATIONS

Despite results obtained in this study were validated by employing another brain atlas, the conclusions were acquired from a single dataset.

CONCLUSIONS

This study revealed two distinguishing neuroanatomical subtypes of first-episode depression, which provides new insights into underlying biological mechanisms of the heterogeneity in depression and might be helpful for accurate clinical diagnosis and future treatment.

Capturing the Individual Deviations From Normative Models of Brain Structure for Depression Diagnosis and Treatment

BACKGROUND

The high heterogeneity of depression prevents us from obtaining reproducible and definite anatomical maps of brain structural changes associated with the disorder, which limits the individualized diagnosis and treatment of patients. In this study, we investigated the clinical issues related to depression according to individual deviations from normative ranges of gray matter volume.

METHODS

We enrolled 1092 participants, including 187 patients with depression and 905 healthy control participants. Structural magnetic resonance imaging data of healthy control participants from the Human Connectome Project (n = 510) and REST-meta-MDD Project (n = 229) were used to establish a normative model across the life span in adults 18 to 65 years old for each brain region. Deviations from the normative range for 187 patients and 166 healthy control participants recruited from two local hospitals were captured as normative probability maps, which were used to identify the disease risk and treatment-related latent factors.

RESULTS

In contrast to case-control results, our normative modeling approach revealed highly individualized patterns of anatomic abnormalities in depressed patients (less than 11% extreme deviation overlapping for any regions). Based on our classification framework, models trained with individual normative probability maps (area under the receiver operating characteristic curve range, 0.7146-0.7836) showed better performance than models trained with original gray matter volume values (area under the receiver operating characteristic curve range, 0.6800-0.7036), which was verified in an independent external test set. Furthermore, different latent brain structural factors in relation to antidepressant treatment were revealed by a Bayesian model based on normative probability maps, suggesting distinct treatment response and inclination.

CONCLUSIONS

Capturing personalized deviations from a normative range could help in understanding the heterogeneous neurobiology of depression and thus guide clinical diagnosis and treatment of depression.

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.

MVComp toolbox: MultiVariate Comparisons of brain MRI features accounting for common information across metrics

Multivariate approaches have recently gained in popularity to address the physiological unspecificity of neuroimaging metrics and to better characterize the complexity of biological processes underlying behavior. However, commonly used approaches are biased by the intrinsic associations between variables, or they are computationally expensive and may be more complicated to implement than standard univariate approaches. Here, we propose using the Mahalanobis distance (D2), an individual-level measure of deviation relative to a reference distribution that accounts for covariance between metrics. To facilitate its use, we introduce an open-source python-based tool for computing D2 relative to a reference group or within a single individual: the MultiVariate Comparison (MVComp) toolbox. The toolbox allows different levels of analysis (i.e., group- or subject-level), resolutions (e.g., voxel-wise, ROI-wise) and dimensions considered (e.g., combining MRI metrics or WM tracts). Several example cases are presented to showcase the wide range of possible applications of MVComp and to demonstrate the functionality of the toolbox. The D2 framework was applied to the assessment of white matter (WM) microstructure at 1) the group-level, where D2 can be computed between a subject and a reference group to yield an individualized measure of deviation. We observed that clustering applied to D2 in the corpus callosum yields parcellations that highly resemble known topography based on neuroanatomy, suggesting that D2 provides an integrative index that meaningfully reflects the underlying microstructure. 2) At the subject level, D2 was computed between voxels to obtain a measure of (dis)similarity. The loadings of each MRI metric (i.e., its relative contribution to D2) were then extracted in voxels of interest to showcase a useful option of the MVComp toolbox. These relative contributions can provide important insights into the physiological underpinnings of differences observed. Integrative multivariate models are crucial to expand our understanding of the complex brain-behavior relationships and the multiple factors underlying disease development and progression. Our toolbox facilitates the implementation of a useful multivariate method, making it more widely accessible.

Functional magnetic resonance imaging in schizophrenia: current evidence, methodological advances, limitations and future directions

Functional neuroimaging emerged with great promise and has provided fundamental insights into the neurobiology of schizophrenia. However, it has faced challenges and criticisms, most notably a lack of clinical translation. This paper provides a comprehensive review and critical summary of the literature on functional neuroimaging, in particular functional magnetic resonance imaging (fMRI), in schizophrenia. We begin by reviewing research on fMRI biomarkers in schizophrenia and the clinical high risk phase through a historical lens, moving from case-control regional brain activation to global connectivity and advanced analytical approaches, and more recent machine learning algorithms to identify predictive neuroimaging features. Findings from fMRI studies of negative symptoms as well as of neurocognitive and social cognitive deficits are then reviewed. Functional neural markers of these symptoms and deficits may represent promising treatment targets in schizophrenia. Next, we summarize fMRI research related to antipsychotic medication, psychotherapy and psychosocial interventions, and neurostimulation, including treatment response and resistance, therapeutic mechanisms, and treatment targeting. We also review the utility of fMRI and data-driven approaches to dissect the heterogeneity of schizophrenia, moving beyond case-control comparisons, as well as methodological considerations and advances, including consortia and precision fMRI. Lastly, limitations and future directions of research in the field are discussed. Our comprehensive review suggests that, in order for fMRI to be clinically useful in the care of patients with schizophrenia, research should address potentially actionable clinical decisions that are routine in schizophrenia treatment, such as which antipsychotic should be prescribed or whether a given patient is likely to have persistent functional impairment. The potential clinical utility of fMRI is influenced by and must be weighed against cost and accessibility factors. Future evaluations of the utility of fMRI in prognostic and treatment response studies may consider including a health economics analysis.

Understanding divergence: Placing developmental neuroscience in its dynamic context

Neurodevelopment is not merely a process of brain maturation, but an adaptation to constraints unique to each individual and to the environments we co-create. However, our theoretical and methodological toolkits often ignore this reality. There is growing awareness that a shift is needed that allows us to study divergence of brain and behaviour across conventional categorical boundaries. However, we argue that in future our study of divergence must also incorporate the developmental dynamics that capture the emergence of those neurodevelopmental differences. This crucial step will require adjustments in study design and methodology. If our ultimate aim is to incorporate the developmental dynamics that capture how, and ultimately when, divergence takes place then we will need an analytic toolkit equal to these ambitions. We argue that the over reliance on group averages has been a conceptual dead-end with regard to the neurodevelopmental differences. This is in part because any individual differences and developmental dynamics are inevitably lost within the group average. Instead, analytic approaches which are themselves new, or simply newly applied within this context, may allow us to shift our theoretical and methodological frameworks from groups to individuals. Likewise, methods capable of modelling complex dynamic systems may allow us to understand the emergent dynamics only possible at the level of an interacting neural system.

The schizophrenia syndrome, circa 2024: What we know and how that informs its nature

With new data about different aspects of schizophrenia being continually generated, it becomes necessary to periodically revisit exactly what we know. Along with a need to review what we currently know about schizophrenia, there is an equal imperative to evaluate the construct itself. With these objectives, we undertook an iterative, multi-phase process involving fifty international experts in the field, with each step building on learnings from the prior one. This review assembles currently established findings about schizophrenia (construct, etiology, pathophysiology, clinical expression, treatment) and posits what they reveal about its nature. Schizophrenia is a heritable, complex, multi-dimensional syndrome with varying degrees of psychotic, negative, cognitive, mood, and motor manifestations. The illness exhibits a remitting and relapsing course, with varying degrees of recovery among affected individuals with most experiencing significant social and functional impairment. Genetic risk factors likely include thousands of common genetic variants that each have a small impact on an individual's risk and a plethora of rare gene variants that have a larger individual impact on risk. Their biological effects are concentrated in the brain and many of the same variants also increase the risk of other psychiatric disorders such as bipolar disorder, autism, and other neurodevelopmental conditions. Environmental risk factors include but are not limited to urban residence in childhood, migration, older paternal age at birth, cannabis use, childhood trauma, antenatal maternal infection, and perinatal hypoxia. Structural, functional, and neurochemical brain alterations implicate multiple regions and functional circuits. Dopamine D-2 receptor antagonists and partial agonists improve psychotic symptoms and reduce risk of relapse. Certain psychological and psychosocial interventions are beneficial. Early intervention can reduce treatment delay and improve outcomes. Schizophrenia is increasingly considered to be a heterogeneous syndrome and not a singular disease entity. There is no necessary or sufficient etiology, pathology, set of clinical features, or treatment that fully circumscribes this syndrome. A single, common pathophysiological pathway appears unlikely. The boundaries of schizophrenia remain fuzzy, suggesting the absence of a categorical fit and need to reconceptualize it as a broader, multi-dimensional and/or spectrum construct.

From Complexity to Precision-Charting Decision-Making Through Normative Modeling

This Viewpoint synthesizes data-driven and theory-driven approaches to normative modeling.

Using normative modeling to assess pharmacological treatment effect on brain state in patients with schizophrenia

Quantifying individual differences in neuroimaging metrics is attracting interest in clinical studies with mental disorders. Schizophrenia is diagnosed exclusively based on symptoms, and the biological heterogeneity makes it difficult to accurately assess pharmacological treatment effects on the brain state. Using the Cambridge Centre for Ageing and Neuroscience data set, we built normative models of brain states and mapped the deviations of the brain characteristics of each patient, to test whether deviations were related to symptoms, and further investigated the pharmacological treatment effect on deviation distributions. Specifically, we found that the patients can be divided into 2 groups: the normalized group had a normalization trend and milder symptoms at baseline, and the other group showed a more severe deviation trend. The baseline severity of the depression as well as the overall symptoms could predict the deviation of the static characteristics for the dorsal and ventral attention networks after treatment. In contrast, the positive symptoms could predict the deviations of the dynamic fluctuations for the default mode and dorsal attention networks after treatment. This work evaluates the effect of pharmacological treatment on static and dynamic brain states using an individualized approach, which may assist in understanding the heterogeneity of the illness pathology as well as the treatment response.

A multiscale characterization of cortical shape asymmetries in early psychosis

Psychosis has often been linked to abnormal cortical asymmetry, but prior results have been inconsistent. Here, we applied a novel spectral shape analysis to characterize cortical shape asymmetries in patients with early psychosis across different spatial scales. We used the Human Connectome Project for Early Psychosis dataset (aged 16-35), comprising 56 healthy controls (37 males, 19 females) and 112 patients with early psychosis (68 males, 44 females). We quantified shape variations of each hemisphere over different spatial frequencies and applied a general linear model to compare differences between healthy controls and patients with early psychosis. We further used canonical correlation analysis to examine associations between shape asymmetries and clinical symptoms. Cortical shape asymmetries, spanning wavelengths from about 22 to 75 mm, were significantly different between healthy controls and patients with early psychosis (Cohen's d = 0.28-0.51), with patients showing greater asymmetry in cortical shape than controls. A single canonical mode linked the asymmetry measures to symptoms (canonical correlation analysis r = 0.45), such that higher cortical asymmetry was correlated with more severe excitement symptoms and less severe emotional distress. Significant group differences in the asymmetries of traditional morphological measures of cortical thickness, surface area, and gyrification, at either global or regional levels, were not identified. Cortical shape asymmetries are more sensitive than other morphological asymmetries in capturing abnormalities in patients with early psychosis. These abnormalities are expressed at coarse spatial scales and are correlated with specific symptom domains.

Resolving heterogeneity in schizophrenia, bipolar I disorder, and attention-deficit/hyperactivity disorder through individualized structural covariance network analysis

Disruptions in large-scale brain connectivity are hypothesized to contribute to psychiatric disorders, including schizophrenia, bipolar I disorder, and attention-deficit/hyperactivity disorder. However, high inter-individual variation among patients with psychiatric disorders hinders achievement of unified findings. To this end, we adopted a newly proposed method to resolve heterogeneity of differential structural covariance network in schizophrenia, bipolar I disorder, and attention-deficit/hyperactivity disorder. This method could infer individualized structural covariance aberrance by assessing the deviation from healthy controls. T1-weighted anatomical images of 114 patients with psychiatric disorders (schizophrenia: n = 37; bipolar I disorder: n = 37; attention-deficit/hyperactivity disorder: n = 37) and 110 healthy controls were analyzed to obtain individualized differential structural covariance network. Patients exhibited tremendous heterogeneity in profiles of individualized differential structural covariance network. Despite notable heterogeneity, patients with the same disorder shared altered edges at network level. Moreover, individualized differential structural covariance network uncovered two distinct psychiatric subtypes with opposite differences in structural covariance edges, that were otherwise obscured when patients were merged, compared with healthy controls. These results provide new insights into heterogeneity and have implications for the nosology in psychiatric disorders.

Linking Polygenic Risk of Schizophrenia to Variation in Magnetic Resonance Imaging Brain Measures: A Comprehensive Systematic Review

BACKGROUND AND HYPOTHESIS

Schizophrenia is highly heritable, with a polygenic effect of many genes conferring risk. Evidence on whether cumulative risk also predicts alterations in brain morphology and function is inconsistent. This systematic review examined evidence for schizophrenia polygenic risk score (sczPRS) associations with commonly used magnetic resonance imaging (MRI) measures. We expected consistent evidence to emerge for significant sczPRS associations with variation in structure and function, specifically in frontal, temporal, and insula cortices that are commonly implicated in schizophrenia pathophysiology.

STUDY DESIGN

In accordance with Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines, we searched MEDLINE, Embase, and PsycINFO for peer-reviewed studies published between January 2013 and March 2022. Studies were screened against predetermined criteria and National Institutes of Health (NIH) quality assessment tools.

STUDY RESULTS

In total, 57 studies of T1-weighted structural, diffusion, and functional MRI were included (age range = 9-80 years, Nrange = 64-76 644). We observed moderate, albeit preliminary, evidence for higher sczPRS predicting global reductions in cortical thickness and widespread variation in functional connectivity, and to a lesser extent, region-specific reductions in frontal and temporal volume and thickness. Conversely, sczPRS does not predict whole-brain surface area or gray/white matter volume. Limited evidence emerged for sczPRS associations with diffusion tensor measures of white matter microstructure in a large community sample and smaller cohorts of children and young adults. These findings were broadly consistent across community and clinical populations.

CONCLUSIONS

Our review supports the hypothesis that schizophrenia is a disorder of disrupted within and between-region brain connectivity, and points to specific whole-brain and regional MRI metrics that may provide useful intermediate phenotypes.

Multi-polygenic scores in psychiatry: From disorder specific to transdiagnostic perspectives

The dense co-occurrence of psychiatric disorders questions the categorical classification tradition and motivates efforts to establish dimensional constructs with neurobiological foundations that transcend diagnostic boundaries. In this study, we examined the genetic liability for eight major psychiatric disorder phenotypes under both a disorder-specific and a transdiagnostic framework. The study sample (n = 513) was deeply phenotyped, consisting of 452 patients from tertiary care with mood disorders, anxiety disorders (ANX), attention-deficit/hyperactivity disorder (ADHD), autism spectrum disorders, and/or substance use disorders (SUD) and 61 unaffected comparison individuals. We computed subject-specific polygenic risk score (PRS) profiles and assessed their associations with psychiatric diagnoses, comorbidity status, as well as cross-disorder behavioral dimensions derived from a rich battery of psychopathology assessments. High PRSs for depression were unselectively associated with the diagnosis of SUD, ADHD, ANX, and mood disorders (p < 1e-4). In the dimensional approach, four distinct functional domains were uncovered, namely the negative valence, social, cognitive, and regulatory systems, closely matching the major functional domains proposed by the Research Domain Criteria (RDoC) framework. Critically, the genetic predisposition for depression was selectively reflected in the functional aspect of negative valence systems (R2  = 0.041, p = 5e-4) but not others. This study adds evidence to the ongoing discussion about the misalignment between current psychiatric nosology and the underlying psychiatric genetic etiology and underscores the effectiveness of the dimensional approach in both the functional characterization of psychiatric patients and the delineation of the genetic liability for psychiatric disorders.

Alzheimer's disease heterogeneity revealed by neuroanatomical normative modeling

INTRODUCTION

Overlooking the heterogeneity in Alzheimer's disease (AD) may lead to diagnostic delays and failures. Neuroanatomical normative modeling captures individual brain variation and may inform our understanding of individual differences in AD-related atrophy.

METHODS

We applied neuroanatomical normative modeling to magnetic resonance imaging from a real-world clinical cohort with confirmed AD (n = 86). Regional cortical thickness was compared to a healthy reference cohort (n = 33,072) and the number of outlying regions was summed (total outlier count) and mapped at individual- and group-levels.

RESULTS

The superior temporal sulcus contained the highest proportion of outliers (60%). Elsewhere, overlap between patient atrophy patterns was low. Mean total outlier count was higher in patients who were non-amnestic, at more advanced disease stages, and without depressive symptoms. Amyloid burden was negatively associated with outlier count.

DISCUSSION

Brain atrophy in AD is highly heterogeneous and neuroanatomical normative modeling can be used to explore anatomo-clinical correlations in individual patients.

Structural variability of the cerebral cortex in schizophrenia and its association with clinical symptoms

BACKGROUND

Substantial evidence indicates structural abnormalities in the cerebral cortex of patients with schizophrenia (SCZ), although their clinical implications remain unclear. Previous case-control studies have investigated group-level differences in structural abnormalities, although the study design cannot account for interindividual differences. Recent research has focused on the association between the heterogeneity of the cerebral cortex morphometric features and clinical heterogeneity.

METHODS

We used neuroimaging data from 420 healthy controls and 695 patients with SCZ from seven studies. Four cerebral cortex measures were obtained: surface area, gray matter volume, thickness, and local gyrification index. We calculated the coefficient of variation (CV) and person-based similarity index (PBSI) scores and performed group comparisons. Associations between the PBSI scores and cognitive functions were evaluated using Spearman's rho test and normative modeling.

RESULTS

Patients with SCZ had a greater CV of surface area and cortical thickness than those of healthy controls. All PBSI scores across cortical measures were lower in patients with SCZ than in HCs. In the patient group, the PBSI scores for gray matter volume and all cortical measures taken together positively correlated with the full-scale IQ scores. Patients with deviant PBSI scores for gray matter volume and all cortical measures taken together had lower full-scale IQ scores than those of other patients.

CONCLUSIONS

The cerebral cortex in patients with SCZ showed greater regional and global structural variability than that in healthy controls. Patients with deviant similarity of cortical structural profiles exhibited a lower general intelligence than those exhibited by the other patients.

What Came First, Mania or Depression? Polarity at Onset in Bipolar I and II: Temperament and Clinical Course

(1) Background: Bipolar disorder (BD) is divided into type I (BD-I) and type II (BD-II). Polarity at onset (PO) is a proposal to specify the clinical course of BD, based on the type of the first episode at disorder onset-depressive (D-PO) or manic (M-PO). At the same time, affective temperaments represent preexisting variants of the spectrum of affective disorders. Our objectives were to investigate the hypothesis that temperament may exert an influence on PO, and that this factor can serve as an indicator of the forthcoming course of the disorder, carrying significant therapeutic implications. (2) Methods: We included 191 patients with BD and examined clinical variables and temperament; the latter was assessed using the short version of the Temperament Evaluation of Memphis, Pisa, Paris, and San Diego-Auto-questionnaire (TEMPS-A-39-SV). We tested the associations between these variables and PO using standard univariate/bivariate methods followed by multivariate logistic regression models. (3) Results: 52.9% of the sample had D-PO and 47.1% had M-PO. D-PO and M-PO patients scored higher for dysthymic and hyperthymic temperaments, respectively (p < 0.001). Also, they differed in BD subtypes, age at first affective episode, illness duration, number of depressive episodes, seasonality, suicide risk, substance use, lithium, and benzodiazepine use (p < 0.05). Only BD-II and age at first depressive episode were predictors of D-PO, whereas BD-I, age at first manic/hypomanic episode, and hyperthymic temperament were predictors of M-PO (p < 0.01). (4) Conclusions: Our findings point to the importance of carefully assessing temperament and PO in patients with BD, to better predict the clinical course and tailor therapeutic interventions to individual patients' needs.

Intra- and inter-individual cognitive variability in schizophrenia and bipolar spectrum disorder: an investigation across multiple cognitive domains

There is substantial cognitive heterogeneity among patients with schizophrenia (SZ) and bipolar disorders (BD). More knowledge about the magnitude and clinical correlates of performance variability could improve our understanding of cognitive impairments. Using double generalized linear models (DGLMs) we investigated cognitive mean and variability differences between patients with SZ (n = 905) and BD spectrum disorders (n = 522), and healthy controls (HC, n = 1170) on twenty-two variables. The analysis revealed significant case-control differences on 90% of the variables. Compared to HC, patients showed larger intra-individual (within subject) variability across tests and larger inter-individual (between subject) variability in measures of fine-motor speed, mental processing speed, and inhibitory control (SZ and BD), and in verbal learning and memory and intellectual functioning (SZ). In SZ, we found that lager intra -and inter (on inhibitory control and speed functions) individual variability, was associated with lower functioning and more negative symptoms. Inter-individual variability on single measures of memory and intellectual function was additionally associated with disorganized and positive symptoms, and use of antidepressants. In BD, there were no within-subject associations with symptom severity. However, greater inter-individual variability (primarily on inhibitory control and speeded functions) was associated with lower functioning, more negative -and disorganized symptoms, earlier age at onset, longer duration of illness, and increased medication use. These results highlight larger individual differences in patients compared to controls on various cognitive domains. Further investigations of the causes and correlates of individual differences in cognitive function are warranted.

Mapping Neurophysiological Subtypes of Major Depressive Disorder Using Normative Models of the Functional Connectome

BACKGROUND

Major depressive disorder (MDD) is a highly heterogeneous disorder that typically emerges in adolescence and can occur throughout adulthood. Studies aimed at quantitatively uncovering the heterogeneity of individual functional connectome abnormalities in MDD and identifying reproducibly distinct neurophysiological MDD subtypes across the lifespan, which could provide promising insights for precise diagnosis and treatment prediction, are still lacking.

METHODS

Leveraging resting-state functional magnetic resonance imaging data from 1148 patients with MDD and 1079 healthy control participants (ages 11-93), we conducted the largest multisite analysis to date for neurophysiological MDD subtyping. First, we characterized typical lifespan trajectories of functional connectivity strength based on the normative model and quantitatively mapped the heterogeneous individual deviations among patients with MDD. Then, we identified neurobiological MDD subtypes using an unsupervised clustering algorithm and evaluated intersite reproducibility. Finally, we validated the subtype differences in baseline clinical variables and longitudinal treatment predictive capacity.

RESULTS

Our findings indicated great intersubject heterogeneity in the spatial distribution and severity of functional connectome deviations among patients with MDD, which inspired the identification of 2 reproducible neurophysiological subtypes. Subtype 1 showed severe deviations, with positive deviations in the default mode, limbic, and subcortical areas and negative deviations in the sensorimotor and attention areas. Subtype 2 showed a moderate but converse deviation pattern. More importantly, subtype differences were observed in depressive item scores and the predictive ability of baseline deviations for antidepressant treatment outcomes.

CONCLUSIONS

These findings shed light on our understanding of different neurobiological mechanisms underlying the clinical heterogeneity of MDD and are essential for developing personalized treatments for this disorder.

Brain-charting autism and attention deficit hyperactivity disorder reveals distinct and overlapping neurobiology

Background

Autism and attention deficit hyperactivity disorder (ADHD) are heterogeneous neurodevelopmental conditions with complex underlying neurobiology. Despite overlapping presentation and sex-biased prevalence, autism and ADHD are rarely studied together, and sex differences are often overlooked. Normative modelling provides a unified framework for studying age-specific and sex-specific divergences in neurodivergent brain development.

Methods

Here we use normative modelling and a large, multi-site neuroimaging dataset to characterise cortical anatomy associated with autism and ADHD, benchmarked against models of typical brain development based on a sample of over 75,000 individuals. We also examined sex and age differences, relationship with autistic traits, and explored the co-occurrence of autism and ADHD (autism+ADHD).

Results

We observed robust neuroanatomical signatures of both autism and ADHD. Overall, autistic individuals showed greater cortical thickness and volume localised to the superior temporal cortex, whereas individuals with ADHD showed more global effects of cortical thickness increases but lower cortical volume and surface area across much of the cortex. The autism+ADHD group displayed a unique pattern of widespread increases in cortical thickness, and certain decreases in surface area. We also found evidence that sex modulates the neuroanatomy of autism but not ADHD, and an age-by-diagnosis interaction for ADHD only.

Conclusions

These results indicate distinct cortical differences in autism and ADHD that are differentially impacted by age, sex, and potentially unique patterns related to their co-occurrence.

Extreme deviations from the normative model reveal cortical heterogeneity and associations with negative symptom severity in first-episode psychosis from the OPTiMiSE and GAP studies

There is currently no quantifiable method to predict long-term clinical outcomes in patients presenting with a first episode of psychosis. A major barrier to developing useful markers for this is biological heterogeneity, where many different pathological mechanisms may underly the same set of symptoms in different individuals. Normative modelling has been used to quantify this heterogeneity in established psychotic disorders by identifying regions of the cortex which are thinner than expected based on a normative healthy population range. These brain atypicalities are measured at the individual level and therefore potentially useful in a clinical setting. However, it is still unclear whether alterations in individual brain structure can be detected at the time of the first psychotic episode, and whether they are associated with subsequent clinical outcomes. We applied normative modelling of cortical thickness to a sample of first-episode psychosis patients, with the aim of quantifying heterogeneity and to use any pattern of cortical atypicality to predict symptoms and response to antipsychotic medication at timepoints from baseline up to 95 weeks (median follow-ups = 4). T1-weighted brain magnetic resonance images from the GAP and OPTiMiSE samples were processed with Freesurfer V6.0.0 yielding 148 cortical thickness features. An existing normative model of cortical thickness (n = 37,126) was adapted to integrate data from each clinical site and account for effects of gender and site. Our test sample consisted of control participants (n = 149, mean age = 26, SD = 6.7) and patient data (n = 295, mean age = 26, SD = 6.7), this sample was used for estimating deviations from the normative model and subsequent statistical analysis. For each individual, the 148 cortical thickness features were mapped to centiles of the normative distribution and converted to z-scores reflecting the distance from the population mean. Individual cortical thickness metrics of +/- 2.6 standard deviations from the mean were considered extreme deviations from the norm. We found that no more than 6.4% of psychosis patients had extreme deviations in a single brain region (regional overlap) demonstrating a high degree of heterogeneity. Mann-Whitney U tests were run on z-scores for each region and significantly lower z-scores were observed in FEP patients in the frontal, temporal, parietal and occipital lobes. Finally, linear mixed-effects modelling showed that negative deviations in cortical thickness in parietal and temporal regions at baseline are related to more severe negative symptoms over the medium-term. This study shows that even at the early stage of symptom onset normative modelling provides a framework to identify individualised cortical markers which can be used for early personalised intervention and stratification.

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

Network-Based Spreading of Gray Matter Changes Across Different Stages of Psychosis

Importance

Psychotic illness is associated with anatomically distributed gray matter reductions that can worsen with illness progression, but the mechanisms underlying the specific spatial patterning of these changes is unknown.

Objective

To test the hypothesis that brain network architecture constrains cross-sectional and longitudinal gray matter alterations across different stages of psychotic illness and to identify whether certain brain regions act as putative epicenters from which volume loss spreads.

Design, Settings, and Participants

This case-control study included 534 individuals from 4 cohorts, spanning early and late stages of psychotic illness. Early-stage cohorts included patients with antipsychotic-naive first-episode psychosis (n = 59) and a group of patients receiving medications within 3 years of psychosis onset (n = 121). Late-stage cohorts comprised 2 independent samples of people with established schizophrenia (n = 136). Each patient group had a corresponding matched control group (n = 218). A sample of healthy adults (n = 356) was used to derive representative structural and functional brain networks for modeling of network-based spreading processes. Longitudinal illness-related and antipsychotic-related gray matter changes over 3 and 12 months were examined using a triple-blind randomized placebo-control magnetic resonance imaging study of the antipsychotic-naive patients. All data were collected between April 29, 2008, and January 15, 2020, and analyses were performed between March 1, 2021, and January 14, 2023.

Main Outcomes and Measures

Coordinated deformation models were used to estimate the extent of gray matter volume (GMV) change in each of 332 parcellated areas by the volume changes observed in areas to which they were structurally or functionally coupled. To identify putative epicenters of volume loss, a network diffusion model was used to simulate the spread of pathology from different seed regions. Correlations between estimated and empirical spatial patterns of GMV alterations were used to quantify model performance.

Results

Of 534 included individuals, 354 (66.3%) were men, and the mean (SD) age was 28.4 (7.4) years. In both early and late stages of illness, spatial patterns of cross-sectional volume differences between patients and controls were more accurately estimated by coordinated deformation models constrained by structural, rather than functional, network architecture (r range, >0.46 to <0.57; P < .01). The same model also robustly estimated longitudinal volume changes related to illness (r ≥ 0.52; P < .001) and antipsychotic exposure (r ≥ 0.50; P < .004). Network diffusion modeling consistently identified, across all 4 data sets, the anterior hippocampus as a putative epicenter of pathological spread in psychosis. Epicenters of longitudinal GMV loss were apparent in posterior cortex early in the illness and shifted to the prefrontal cortex with illness progression.

Conclusion and Relevance

These findings highlight a central role for white matter fibers as conduits for the spread of pathology across different stages of psychotic illness, mirroring findings reported in neurodegenerative conditions. The structural connectome thus represents a fundamental constraint on brain changes in psychosis, regardless of whether these changes are caused by illness or medication. Moreover, the anterior hippocampus represents a putative epicenter of early brain pathology from which dysfunction may spread to affect connected areas.

Understanding the complexity of individual developmental pathways: A primer on metaphors, models, and methods to study resilience in development

The modern study of resilience in development is conceptually based on a complex adaptive system ontology in which many (intersystem) factors are involved in the emergence of resilient developmental pathways. However, the methods and models developed to study complex dynamical systems have not been widely adopted, and it has recently been noted this may constitute a problem moving the field forward. In the present paper, I argue that an ontological commitment to complex adaptive systems is not only possible, but highly recommended for the study of resilience in development. Such a commitment, however, also comes with a commitment to a different causal ontology and different research methods. In the first part of the paper, I discuss the extent to which current research on resilience in development conceptually adheres to the complex systems perspective. In the second part, I introduce conceptual tools that may help researchers conceptualize causality in complex systems. The third part discusses idiographic methods that could be used in a research program that embraces the interaction dominant causal ontology and idiosyncratic nature of the dynamics of complex systems. The conclusion is that a strong ontological commitment is warranted, but will require a radical departure from nomothetic science.

Predictive Clinical Neuroscience Portal (PCNportal): instant online access to research-grade normative models for clinical neuroscientists

Background

The neurobiology of mental disorders remains poorly understood despite substantial scientific efforts, due to large clinical heterogeneity and to a lack of tools suitable to map individual variability. Normative modeling is one recently successful framework that can address these problems by comparing individuals to a reference population. The methodological underpinnings of normative modelling are, however, relatively complex and computationally expensive. Our research group has developed the python-based normative modelling package Predictive Clinical Neuroscience toolkit (PCNtoolkit) which provides access to many validated algorithms for normative modelling. PCNtoolkit has since proven to be a strong foundation for large scale normative modelling, but still requires significant computation power, time and technical expertise to develop.

Methods

To address these problems, we introduce PCNportal. PCNportal is an online platform integrated with PCNtoolkit that offers access to pre-trained research-grade normative models estimated on tens of thousands of participants, without the need for computation power or programming abilities. PCNportal is an easy-to-use web interface that is highly scalable to large user bases as necessary. Finally, we demonstrate how the resulting normalized deviation scores can be used in a clinical application through a schizophrenia classification task applied to cortical thickness and volumetric data from the longitudinal Northwestern University Schizophrenia Data and Software Tool (NUSDAST) dataset.

Results

At each longitudinal timepoint, the transferred normative models achieved a mean[std. dev.] explained variance of 9.4[8.8]%, 9.2[9.2]%, 5.6[7.4]% respectively in the control group and 4.7[5.5]%, 6.0[6.2]%, 4.2[6.9]% in the schizophrenia group. Diagnostic classifiers achieved AUC of 0.78, 0.76 and 0.71 respectively.

Conclusions

This replicates the utility of normative models for diagnostic classification of schizophrenia and showcases the use of PCNportal for clinical neuroimaging. By facilitating and speeding up research with high-quality normative models, this work contributes to research in inter-individual variability, clinical heterogeneity and precision medicine.

A normative model of brain responses to social scenarios reflects the maturity of children and adolescents' social-emotional abilities

The rapid brain maturation in childhood and adolescence accompanies the development of socio-emotional functioning. However, it is unclear how the maturation of the neural activity drives the development of socio-emotional functioning and individual differences. This study aimed to reflect the age dependence of inter-individual differences in brain responses to socio-emotional scenarios and to develop naturalistic imaging indicators to assess the maturity of socio-emotional ability at the individual level. Using three independent naturalistic imaging datasets containing healthy participants (n = 111, 21 and 122), we found and validated that age-modulated inter-individual concordance of brain responses to socio-emotional movies in specific brain regions. The similarity of an individual's brain response to the average response of older participants was defined as response typicality, which predicted an individual's emotion regulation strategies in adolescence and theory of mind (ToM) in childhood. Its predictive power was not superseded by age, sex, cognitive performance or executive function. We further showed that the movie's valence and arousal ratings grounded the response typicality. The findings highlight that forming typical brain response patterns may be a neural phenotype underlying the maturation of socio-emotional ability. The proposed response typicality represents a neuroimaging approach to measure individuals' maturity of cognitive reappraisal and ToM.

Diurnal oscillations of MRI metrics in the brains of male participants

Regulation of biological processes according to a 24-hr rhythm is essential for the normal functioning of an organism. Temporal variation in brain MRI data has often been attributed to circadian or diurnal oscillations; however, it is not clear if such oscillations exist. Here, we provide evidence that diurnal oscillations indeed govern multiple MRI metrics. We recorded cerebral blood flow, diffusion-tensor metrics, T1 relaxation, and cortical structural features every three hours over a 24-hr period in each of 16 adult male controls and eight adult male participants with bipolar disorder. Diurnal oscillations are detected in numerous MRI metrics at the whole-brain level, and regionally. Rhythmicity parameters in the participants with bipolar disorder are similar to the controls for most metrics, except for a larger phase variation in cerebral blood flow. The ubiquitous nature of diurnal oscillations has broad implications for neuroimaging studies and furthers our understanding of the dynamic nature of the human brain.

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.