Replicating extensive brain structural heterogeneity in individuals with schizophrenia and bipolar disorder

Polyconnectomic scoring of functional connectivity patterns across eight neuropsychiatric and three neurodegenerative disorders

BACKGROUND

Neuropsychiatric and neurodegenerative disorders involve diverse changes in brain functional connectivity. As an alternative to approaches searching for specific mosaic patterns of affected connections and networks, we used polyconnectomic scoring to quantify disorder-related whole-brain connectivity signatures into interpretable, personalized scores.

METHODS

The polyconnectomic score (PCS) measures the extent to which an individual's functional connectivity (FC) mirrors the whole-brain circuitry characteristics of a trait. We computed PCS for eight neuropsychiatric conditions (attention-deficit/hyperactivity disorder, anxiety-related disorders, autism spectrum disorder, obsessive-compulsive disorder, bipolar disorder, major depressive disorder, schizoaffective disorder, and schizophrenia) and three neurodegenerative conditions (Alzheimer's disease, frontotemporal dementia, and Parkinson's disease) across 22 datasets with resting-state functional MRI of 10,667 individuals (5,325 patients, 5,342 controls). We further examined PCS in 26,673 individuals from the population-based UK Biobank cohort.

RESULTS

PCS was consistently higher in out-of-sample patients across six of the eight neuropsychiatric and across all three investigated neurodegenerative disorders ([min, max]: AUC = [0.55, 0.73], pFDR = [1.8 x 10-16, 4.5 x 10-2]). Individuals with elevated PCS levels for neuropsychiatric conditions exhibited higher neuroticism (pFDR < 9.7 x 10-5), lower cognitive performance (pFDR < 5.3 x 10-5), and lower general wellbeing (pFDR < 9.7 x 10-4).

CONCLUSIONS

Our findings reveal generalizable whole-brain connectivity alterations in brain disorders. PCS effectively aggregates disorder-related signatures across the entire brain into an interpretable, subject-specific metric. A toolbox is provided for PCS computation.

A generalizable normative deep autoencoder for brain morphological anomaly detection: application to the multi-site StratiBip dataset on bipolar disorder in an external validation framework

The heterogeneity of psychiatric disorders makes researching disorder-specific neurobiological markers an ill-posed problem. Here, we face the need for disease stratification models by presenting a generalizable multivariate normative modelling framework for characterizing brain morphology, applied to bipolar disorder (BD). We employed deep autoencoders in an anomaly detection framework, combined with a confounder removal step integrating training and external validation. The model was trained with healthy control (HC) data from the human connectome project and applied to multi-site external data of HC and BD individuals. We found that brain deviating scores were greater, more heterogeneous, and with increased extreme values in the BD group, with volumes prominently from the basal ganglia, hippocampus and adjacent regions emerging as significantly deviating. Similarly, individual brain deviating maps based on modified z scores expressed higher abnormalities occurrences, but their overall spatial overlap was lower compared to HCs. Our generalizable framework enabled the identification of subject- and group-level brain normative-deviating patterns, a step forward towards the development of more effective and personalized clinical decision support systems and patient stratification in psychiatry.

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.

Assessing regional intracortical myelination in schizophrenia spectrum and bipolar disorders using the optimized T1w/T2w-ratio

BACKGROUND

Dysmyelination could be part of the pathophysiology of schizophrenia spectrum (SCZ) and bipolar disorders (BPD), yet few studies have examined myelination of the cerebral cortex. The ratio of T1- and T2-weighted magnetic resonance images (MRI) correlates with intracortical myelin. We investigated the T1w/T2w-ratio and its age trajectories in patients and healthy controls (CTR) and explored associations with antipsychotic medication use and psychotic symptoms.

METHODS

Patients with SCZ (n = 64; mean age = 30.4 years, s.d. = 9.8), BPD (n = 91; mean age 31.0 years, s.d. = 10.2), and CTR (n = 155; mean age = 31.9 years, s.d. = 9.1) who participated in the TOP study (NORMENT, University of Oslo, Norway) were clinically assessed and scanned using a General Electric 3 T MRI system. T1w/T2w-ratio images were computed using an optimized pipeline with intensity normalization and field inhomogeneity correction. Vertex-wise regression models were used to compare groups and examine group × age interactions. In regions showing significant differences, we explored associations with antipsychotic medication use and psychotic symptoms.

RESULTS

No main effect of diagnosis was found. However, age slopes of the T1w/T2w-ratio differed significantly between SCZ and CTR, predominantly in frontal and temporal lobe regions: Lower T1w/T2w-ratio values with higher age were found in CTR, but not in SCZ. Follow-up analyses revealed a more positive age slope in patients who were using antipsychotics and patients using higher chlorpromazine-equivalent doses.

CONCLUSIONS

While we found no evidence of reduced intracortical myelin in SCZ or BPD relative to CTR, different regional age trajectories in SCZ may suggest a promyelinating effect of antipsychotic medication.

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.

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.

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.

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.

Regional, circuit and network heterogeneity of brain abnormalities in psychiatric disorders

The substantial individual heterogeneity that characterizes people with mental illness is often ignored by classical case-control research, which relies on group mean comparisons. Here we present a comprehensive, multiscale characterization of the heterogeneity of gray matter volume (GMV) differences in 1,294 cases diagnosed with one of six conditions (attention-deficit/hyperactivity disorder, autism spectrum disorder, bipolar disorder, depression, obsessive-compulsive disorder and schizophrenia) and 1,465 matched controls. Normative models indicated that person-specific deviations from population expectations for regional GMV were highly heterogeneous, affecting the same area in <7% of people with the same diagnosis. However, these deviations were embedded within common functional circuits and networks in up to 56% of cases. The salience-ventral attention system was implicated transdiagnostically, with other systems selectively involved in depression, bipolar disorder, schizophrenia and attention-deficit/hyperactivity disorder. Phenotypic differences between cases assigned the same diagnosis may thus arise from the heterogeneous localization of specific regional deviations, whereas phenotypic similarities may be attributable to the dysfunction of common functional circuits and networks.

Abnormal temporal variability of rich-club organization in three major psychiatric conditions

Introduction

Convergent evidence has demonstrated a shared rich-club reorganization across multiple major psychiatric conditions. However, previous studies assessing altered functional couplings between rich-club regions have typically focused on the mean time series from entire functional magnetic resonance imaging (fMRI) scanning session, neglecting their time-varying properties.

Methods

In this study, we aim to explore the common and/or unique alterations in the temporal variability of rich-club organization among schizophrenia (SZ), bipolar disorder (BD), and attention deficit/hyperactivity disorder (ADHD). We employed a temporal rich-club (TRC) approach to quantitatively assess the propensity of well-connected nodes to form simultaneous and stable structures in a temporal network derived from resting-state fMRI data of 156 patients with major psychiatric disorders (SZ/BD/ADHD = 71/45/40) and 172 healthy controls. We executed the TRC workflow at both whole-brain and subnetwork scales across varying network sparsity, sliding window strategies, lengths and steps of sliding windows, and durations of TRC coefficients.

Results

The SZ and BD groups displayed significantly decreased TRC coefficients compared to corresponding HC groups at the whole-brain scale and in most subnetworks. In contrast, the ADHD group exhibited reduced TRC coefficients in longer durations, as opposed to shorter durations, which markedly differs from the SZ and BD groups. These findings reveal both transdiagnostic and illness-specific patterns in temporal variability of rich-club organization across SZ, BD, and ADHD.

Discussion

TRC may serve as an effective metric for detecting brain network disruptions in particular states, offering novel insights and potential biomarkers into the neurobiological basis underpinning the behavioral and cognitive deficits observed in these disorders.

Genetic architecture of brain age and its causal relations with brain and mental disorders

The difference between chronological age and the apparent age of the brain estimated from brain imaging data-the brain age gap (BAG)-is widely considered a general indicator of brain health. Converging evidence supports that BAG is sensitive to an array of genetic and nongenetic traits and diseases, yet few studies have examined the genetic architecture and its corresponding causal relationships with common brain disorders. Here, we estimate BAG using state-of-the-art neural networks trained on brain scans from 53,542 individuals (age range 3-95 years). A genome-wide association analysis across 28,104 individuals (40-84 years) from the UK Biobank revealed eight independent genomic regions significantly associated with BAG (p < 5 × 10-8) implicating neurological, metabolic, and immunological pathways - among which seven are novel. No significant genetic correlations or causal relationships with BAG were found for Parkinson's disease, major depressive disorder, or schizophrenia, but two-sample Mendelian randomization indicated a causal influence of AD (p = 7.9 × 10-4) and bipolar disorder (p = 1.35 × 10-2) on BAG. These results emphasize the polygenic architecture of brain age and provide insights into the causal relationship between selected neurological and neuropsychiatric disorders and BAG.

Evidence for embracing normative modeling

In this work, we expand the normative model repository introduced in Rutherford et al., 2022a to include normative models charting lifespan trajectories of structural surface area and brain functional connectivity, measured using two unique resting-state network atlases (Yeo-17 and Smith-10), and an updated online platform for transferring these models to new data sources. We showcase the value of these models with a head-to-head comparison between the features output by normative modeling and raw data features in several benchmarking tasks: mass univariate group difference testing (schizophrenia versus control), classification (schizophrenia versus control), and regression (predicting general cognitive ability). Across all benchmarks, we show the advantage of using normative modeling features, with the strongest statistically significant results demonstrated in the group difference testing and classification tasks. We intend for these accessible resources to facilitate the wider adoption of normative modeling across the neuroimaging community.

Visual processing deficits in patients with schizophrenia spectrum and bipolar disorders and associations with psychotic symptoms, and intellectual abilities

Objective

Low-level sensory disruption is hypothesized as a precursor to clinical and cognitive symptoms in severe mental disorders. We compared visual discrimination performance in patients with schizophrenia spectrum disorder or bipolar disorder with healthy controls, and investigated associations with clinical symptoms and IQ.

Methods

Patients with schizophrenia spectrum disorder (n = 32), bipolar disorder (n = 55) and healthy controls (n = 152) completed a computerized visual discrimination task. Participants responded whether the latter of two consecutive grids had higher or lower spatial frequency, and discrimination thresholds were estimated using an adaptive maximum likelihood procedure. Case-control differences in threshold were assessed using linear regression, F-test and post-hoc pair-wise comparisons. Linear models were used to test for associations between visual discrimination threshold and psychotic symptoms derived from the PANSS and IQ assessed using the Matrix Reasoning and Vocabulary subtests from the Wechsler Abbreviated Scale of Intelligence (WASI).

Results

Robust regression revealed a significant main effect of diagnosis on discrimination threshold (robust F = 6.76, p = .001). Post-hoc comparisons revealed that patients with a schizophrenia spectrum disorder (mean = 14%, SD = 0.08) had higher thresholds compared to healthy controls (mean = 10.8%, SD = 0.07, β = 0.35, t = 3.4, p = .002), as did patients with bipolar disorder (12.23%, SD = 0.07, β = 0.21, t = 2.42, p = .04). There was no significant difference between bipolar disorder and schizophrenia (β = -0.14, t = -1.2, p = .45). Linear models revealed negative associations between IQ and threshold across all participants when controlling for diagnostic group (β = -0.3, t = -3.43, p = .0007). This association was found within healthy controls (t = -3.72, p = .0003) and patients with bipolar disorder (t = -2.53, p = .015), and no significant group by IQ interaction on threshold (F = 0.044, p = .97). There were no significant associations between PANSS domain scores and discrimination threshold.

Conclusion

Patients with schizophrenia spectrum or bipolar disorders exhibited higher visual discrimination thresholds than healthy controls, supporting early visual deficits among patients with severe mental illness. Discrimination threshold was negatively associated with IQ among healthy controls and bipolar disorder patients. These findings elucidate perception-related disease mechanisms in severe mental illness, which warrants replication in independent samples.

Normative modeling of brain morphometry in Clinical High-Risk for Psychosis

Importance

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

Objective

To quantify deviations from the normative range of neuroanatomical variation in individuals at clinical high-risk for psychosis (CHR-P) and evaluate their overlap with healthy variation and their association with positive symptoms, cognition, and conversion to a psychotic disorder.

Design Setting and Participants

Clinical, IQ and FreeSurfer-derived regional measures of cortical thickness (CT), cortical surface area (SA), and subcortical volume (SV) from 1,340 CHR-P individuals [47.09% female; mean age: 20.75 (4.74) years] and 1,237 healthy individuals [44.70% female; mean age: 22.32 (4.95) years] from 29 international sites participating in the ENIGMA Clinical High Risk for Psychosis Working Group.

Main Outcomes and Measures

For each regional morphometric measure, z-scores were computed that index the degree of deviation from the normative means of that measure in a healthy reference population (N=37,407). Average deviation scores (ADS) for CT, SA, SV, and globally across all measures (G) were generated by averaging the respective regional z-scores. Regression analyses were used to quantify the association of deviation scores with clinical severity and cognition and two-proportion z-tests to identify case-control differences in the proportion of individuals with infranormal (z<-1.96) or supranormal (z>1.96) scores.

Results

CHR-P and healthy individuals overlapped in the distributions of the observed values, regional z-scores, and all ADS vales. The proportion of CHR-P individuals with infranormal or supranormal values in any metric was low (<12%) and similar to that of healthy individuals. CHR-P individuals who converted to psychosis compared to those who did not convert had a higher percentage of infranormal values in temporal regions (5-7% vs 0.9-1.4%). In the CHR-P group, only the ADSSA showed significant but weak associations (|β|<0.09; PFDR<0.05) with positive symptoms and IQ.

Conclusions and Relevance

The study findings challenge the usefulness of macroscale neuromorphometric measures as diagnostic biomarkers of psychosis risk and suggest that such measures do not provide an adequate explanation for psychosis risk.

Cumulative Impact of Morphometric Features in Schizophrenia in Two Independent Samples

Schizophrenia and bipolar disorder share a common structural brain alteration profile. However, there is considerable between- and within-diagnosis variability in these features, which may underestimate informative individual differences. Using a recently established morphometric risk score (MRS) approach, we aim to provide confirmation that individual MRS scores are higher in individuals with a psychosis diagnosis, helping to parse individual heterogeneity. Using the Human Connectome Project Early Psychosis (N = 124), we estimate MRS for psychosis and specifically for bipolar/schizophrenia using T1-weighted MRI data and prior meta-analysis effect sizes. We confirm associations in an independent replication sample (N = 69). We assess (1) the impact of diagnosis on these MRS, (2) compare effect sizes of MRS to all individual, cytoarchitecturally defined brain regions, and (3) perform negative control analyses to assess MRS specificity. The MRS specifically for SCZ was higher in the whole psychosis group (Cohen's d = 0.56; P = 0.003) and outperformed any single region of interest in standardized mean difference (ZMRS>75 ROIS = 2.597; P = 0.009) and correlated with previously reported effect sizes (PSPIN/SHUFFLE < 0.005). MRS without Enhancing Neuroimaging Genomics through Meta-Analysis weights did not delineate groups with empirically null associations (t = 2.29; P = 0.02). We replicate MRS specifically for SCZ associations in the independent sample. Akin to polygenic risk scoring and individual allele effect sizes, these observations suggest that assessing the combined impact of regional structural alterations may be more informative than any single cytoarchitecturally constrained cortical region, where well-powered, meta-analytical samples are informative in the delineation of diagnosis and within psychosis case differences, in smaller independent samples.

Closing the life-cycle of normative modeling using federated hierarchical Bayesian regression

Clinical neuroimaging data availability has grown substantially in the last decade, providing the potential for studying heterogeneity in clinical cohorts on a previously unprecedented scale. Normative modeling is an emerging statistical tool for dissecting heterogeneity in complex brain disorders. However, its application remains technically challenging due to medical data privacy issues and difficulties in dealing with nuisance variation, such as the variability in the image acquisition process. Here, we approach the problem of estimating a reference normative model across a massive population using a massive multi-center neuroimaging dataset. To this end, we introduce a federated probabilistic framework using hierarchical Bayesian regression (HBR) to complete the life-cycle of normative modeling. The proposed model provides the possibilities to learn, update, and adapt the model parameters on decentralized neuroimaging data. Our experimental results confirm the superiority of HBR in deriving more accurate normative ranges on large multi-site neuroimaging datasets compared to the current standard methods. In addition, our approach provides the possibility to recalibrate and reuse the learned model on local datasets and even on datasets with very small sample sizes. The proposed method will facilitate applications of normative modeling as a medical tool for screening the biological deviations in individuals affected by complex illnesses such as mental disorders.

Clinical and cortical similarities identified between bipolar disorder I and schizophrenia: A multivariate approach

BACKGROUND

Structural neuroimaging studies have identified similarities in the brains of individuals diagnosed with schizophrenia (SZ) and bipolar I disorder (BP), with overlap in regions of gray matter (GM) deficits between the two disorders. Recent studies have also shown that the symptom phenotypes associated with SZ and BP may allow for a more precise categorization than the current diagnostic criteria. In this study, we sought to identify GM alterations that were unique to each disorder and whether those alterations were also related to unique symptom profiles.

MATERIALS AND METHODS

We analyzed the GM patterns and clinical symptom presentations using independent component analysis (ICA), hierarchical clustering, and n-way biclustering in a large (N ∼ 3,000), merged dataset of neuroimaging data from healthy volunteers (HV), and individuals with either SZ or BP.

RESULTS

Component A showed a SZ and BP < HV GM pattern in the bilateral insula and cingulate gyrus. Component B showed a SZ and BP < HV GM pattern in the cerebellum and vermis. There were no significant differences between diagnostic groups in these components. Component C showed a SZ < HV and BP GM pattern bilaterally in the temporal poles. Hierarchical clustering of the PANSS scores and the ICA components did not yield new subgroups. N-way biclustering identified three unique subgroups of individuals within the sample that mapped onto different combinations of ICA components and symptom profiles categorized by the PANSS but no distinct diagnostic group differences.

CONCLUSION

These multivariate results show that diagnostic boundaries are not clearly related to structural differences or distinct symptom profiles. Our findings add support that (1) BP tend to have less severe symptom profiles when compared to SZ on the PANSS without a clear distinction, and (2) all the gray matter alterations follow the pattern of SZ < BP < HV without a clear distinction between SZ and BP.

Brain anomalies in early psychosis: From secondary to primary psychosis

Brain anomalies are frequently found in early psychoses. Although they may remain undetected for many years, their interpretation is critical for differential diagnosis. In secondary psychoses, their identification may allow specific management. They may also shed light on various pathophysiological aspects of primary psychoses. Here we reviewed cases of secondary psychoses associated with brain anomalies, reported over a 20-year period in adolescents and young adults aged 13-30 years old. We considered age at first psychotic symptoms, relevant medical history, the nature of psychiatric symptoms, clinical red flags, the nature of the brain anomaly reported, and the underlying disease. We discuss the relevance of each brain area in light of normal brain function, recent case-control studies, and postulated pathophysiology. We show that anomalies in all regions, whether diffuse, multifocal, or highly localized, may lead to psychosis, without necessarily being associated with non-psychiatric symptoms. This underlines the interest of neuroimaging in the initial workup, and supports the hypothesis of psychosis as a global network dysfunction that involves many different regions.

The normative modeling framework for computational psychiatry

Normative modeling is an emerging and innovative framework for mapping individual differences at the level of a single subject or observation in relation to a reference model. It involves charting centiles of variation across a population in terms of mappings between biology and behavior, which can then be used to make statistical inferences at the level of the individual. The fields of computational psychiatry and clinical neuroscience have been slow to transition away from patient versus 'healthy' control analytic approaches, probably owing to a lack of tools designed to properly model biological heterogeneity of mental disorders. Normative modeling provides a solution to address this issue and moves analysis away from case-control comparisons that rely on potentially noisy clinical labels. Here we define a standardized protocol to guide users through, from start to finish, normative modeling analysis using the Predictive Clinical Neuroscience toolkit (PCNtoolkit). We describe the input data selection process, provide intuition behind the various modeling choices and conclude by demonstrating several examples of downstream analyses that the normative model may facilitate, such as stratification of high-risk individuals, subtyping and behavioral predictive modeling. The protocol takes ~1-3 h to complete.

Dopaminergic dysfunction and excitatory/inhibitory imbalance in treatment-resistant schizophrenia and novel neuromodulatory treatment

Antipsychotic drugs are the mainstay in the treatment of schizophrenia. However, one-third of patients do not show adequate improvement in positive symptoms with non-clozapine antipsychotics. Additionally, approximately half of them show poor response to clozapine, electroconvulsive therapy, or other augmentation strategies. However, the development of novel treatment for these conditions is difficult due to the complex and heterogenous pathophysiology of treatment-resistant schizophrenia (TRS). Therefore, this review provides key findings, potential treatments, and a roadmap for future research in this area. First, we review the neurobiological pathophysiology of TRS, particularly the dopaminergic, glutamatergic, and GABAergic pathways. Next, the limitations of existing and promising treatments are presented. Specifically, this article focuses on the therapeutic potential of neuromodulation, including electroconvulsive therapy, repetitive transcranial magnetic stimulation, transcranial direct current stimulation, and deep brain stimulation. Finally, we propose multivariate analyses that integrate various perspectives of the pathogenesis, such as dopaminergic dysfunction and excitatory/inhibitory imbalance, thereby elucidating the heterogeneity of TRS that could not be obtained by conventional statistics. These analyses can in turn lead to a precision medicine approach with closed-loop neuromodulation targeting the detected pathophysiology of TRS.

Longitudinal Structural Brain Changes in Bipolar Disorder: A Multicenter Neuroimaging Study of 1232 Individuals by the ENIGMA Bipolar Disorder Working Group

BACKGROUND

Bipolar disorder (BD) is associated with cortical and subcortical structural brain abnormalities. It is unclear whether such alterations progressively change over time, and how this is related to the number of mood episodes. To address this question, we analyzed a large and diverse international sample with longitudinal magnetic resonance imaging (MRI) and clinical data to examine structural brain changes over time in BD.

METHODS

Longitudinal structural MRI and clinical data from the ENIGMA (Enhancing Neuro Imaging Genetics through Meta Analysis) BD Working Group, including 307 patients with BD and 925 healthy control subjects, were collected from 14 sites worldwide. Male and female participants, aged 40 ± 17 years, underwent MRI at 2 time points. Cortical thickness, surface area, and subcortical volumes were estimated using FreeSurfer. Annualized change rates for each imaging phenotype were compared between patients with BD and healthy control subjects. Within patients, we related brain change rates to the number of mood episodes between time points and tested for effects of demographic and clinical variables.

RESULTS

Compared with healthy control subjects, patients with BD showed faster enlargement of ventricular volumes and slower thinning of the fusiform and parahippocampal cortex (0.18 <d < 0.22). More (hypo)manic episodes were associated with faster cortical thinning, primarily in the prefrontal cortex.

CONCLUSIONS

In the hitherto largest longitudinal MRI study on BD, we did not detect accelerated cortical thinning but noted faster ventricular enlargements in BD. However, abnormal frontocortical thinning was observed in association with frequent manic episodes. Our study yields insights into disease progression in BD and highlights the importance of mania prevention in BD treatment.

Charting brain growth and aging at high spatial precision

Defining reference models for population variation, and the ability to study individual deviations is essential for understanding inter-individual variability and its relation to the onset and progression of medical conditions. In this work, we assembled a reference cohort of neuroimaging data from 82 sites (N=58,836; ages 2-100) and used normative modeling to characterize lifespan trajectories of cortical thickness and subcortical volume. Models are validated against a manually quality checked subset (N=24,354) and we provide an interface for transferring to new data sources. We showcase the clinical value by applying the models to a transdiagnostic psychiatric sample (N=1985), showing they can be used to quantify variability underlying multiple disorders whilst also refining case-control inferences. These models will be augmented with additional samples and imaging modalities as they become available. This provides a common reference platform to bind results from different studies and ultimately paves the way for personalized clinical decision-making.

Shared and distinct patterns of dynamical degree centrality in bipolar disorder across different mood states

BACKGROUND

Previous studies have probed the brain static activity pattern in bipolar disorder across different states. However, human intrinsic brain activity is time-varying and dynamic. There is a lack of knowledge about the brain dynamical pattern in bipolar disorder across different mood states.

METHODS

This study used the dynamical degree centrality (dDC) to investigate the resting-state whole-brain dynamical pattern voxel-wise in a total of 62 bipolar disorder [28 bipolar depression (BD), 13 bipolar mania (BM), 21 bipolar euthymia (BE)], and 30 healthy controls (HCs). One-way analysis of variance (ANOVA) was applied to explore the omnibus differences of the dDC pattern across all groups, and Pearson's correlation analysis was used to evaluate the relationship between the dDC variability in detected regions with clinical symptom severity.

RESULTS

One-way ANOVA analysis showed the omnibus differences in the left inferior parietal lobule/middle occipital gyrus (IPL/MOG) and right precuneus/posterior cingulate cortex (PCUN/PCC) across all groups. The post hoc analysis revealed that BD showed decreased dDC in the IPL/MOG compared with all other groups, and both BD and BM exhibited decreased dDC in the PCUN/PCC compared with BE and HCs. Furthermore, correlation analysis showed that the dDC variability of the IPL/MOG and PCUN/PCC negatively correlated with the depression symptom levels in all patients with bipolar disorder.

CONCLUSION

This study demonstrated the distinct and shared brain dynamical pattern of the depressive, manic, and euthymia states. Our findings provide new insights into the pathophysiology of bipolar disorder across different mood states from the dynamical brain network pattern perspective.

Thalamic Shape Abnormalities Differentially Relate to Cognitive Performance in Early-Onset and Adult-Onset Schizophrenia

Early-onset schizophrenia (EOS) shares many biological and clinical features with adult-onset schizophrenia (AOS), but may represent a unique subgroup with greater susceptibility for disease onset and worsened symptomatology and progression, which could potentially derive from exaggerated neurodevelopmental abnormalities. Neurobiological explanations of schizophrenia have emphasized the involvement of deep-brain structures, particularly alterations of the thalamus, which have been linked to core features of the disorder. The aim of this study was to compare thalamic shape abnormalities between EOS and AOS subjects and determine whether unique behavioral profiles related to these differences. It was hypothesized abnormal thalamic shape would be observed in anterior, mediodorsal and pulvinar regions in both schizophrenia groups relative to control subjects, but exacerbated in EOS. Magnetic resonance T1-weighted images were collected from adult individuals with EOS (n = 28), AOS (n = 33), and healthy control subjects (n = 60), as well as collection of clinical and cognitive measures. Large deformation high-dimensional brain mapping was used to obtain three-dimensional surfaces of the thalamus. General linear models were used to compare groups on surface shape features, and Pearson correlations were used to examine relationships between thalamic shape and behavioral measures. Results revealed both EOS and AOS groups demonstrated significant abnormal shape of anterior, lateral and pulvinar thalamic regions relative to CON (all p < 0.007). Relative to AOS, EOS exhibited exacerbated abnormalities in posterior lateral, mediodorsal and lateral geniculate thalamic regions (p = 0.003). Thalamic abnormalities related to worse episodic memory in EOS (p = 0.03) and worse working memory (p = 0.047) and executive functioning (p = 0003) in AOS. Overall, findings suggest thalamic abnormalities are a prominent feature in both early- and late-onset schizophrenia, but exaggerated in EOS and have different brain-behavior profiles for each. The persistence of these abnormalities in adult EOS patients suggests they may represent markers of disrupted neurodevelopment that uniquely relate to the clinical and cognitive aspects of the illness.

Basal ganglia shape features differentiate schizoaffective disorder from schizophrenia

There is growing evidence that schizophrenia and schizoaffective disorder represent closely related syndromes that vary in severity along a neurobiological continuum. In the present study, volume and shape of the basal ganglia was examined in people with schizophrenia and schizoaffective disorder relative to healthy controls and hypothesized that unique neuroanatomical differences would be observed in each patient group. Magnetic resonance 1.5T images were obtained from schizophrenia (n = 47), schizoaffective disorder (n = 15), and from healthy control (n = 42) participants, matched for age, gender, parental socioeconomic status, and race. The caudate, putamen, and globus pallidus were characterized using high-dimensional brain mapping procedures (Csernansky et al., 2004b). Results revealed significant shape deformations between schizophrenia and schizoaffective disorder that also differed from control subjects. Relative to schizophrenia, schizoaffective subjects showed exaggerated inward deformations indicative of localized volume loss in subregions of the caudate, putamen, and globus pallidus (all p < 0.001). These shape features correlated with mental flexibility and negative symptoms in schizophrenia (all p < 0.05), but not schizoaffective disorder. To the extent that differences in important basal ganglia substructures reflect biological heterogeneity among these two psychotic illnesses, this data could prove useful in improving diagnostic precision, as well as informing the affective component of mental illness.

Improving the predictive potential of diffusion MRI in schizophrenia using normative models-Towards subject-level classification

Diffusion MRI studies consistently report group differences in white matter between individuals diagnosed with schizophrenia and healthy controls. Nevertheless, the abnormalities found at the group-level are often not observed at the individual level. Among the different approaches aiming to study white matter abnormalities at the subject level, normative modeling analysis takes a step towards subject-level predictions by identifying affected brain locations in individual subjects based on extreme deviations from a normative range. Here, we leveraged a large harmonized diffusion MRI dataset from 512 healthy controls and 601 individuals diagnosed with schizophrenia, to study whether normative modeling can improve subject-level predictions from a binary classifier. To this aim, individual deviations from a normative model of standard (fractional anisotropy) and advanced (free-water) dMRI measures, were calculated by means of age and sex-adjusted z-scores relative to control data, in 18 white matter regions. Even though larger effect sizes are found when testing for group differences in z-scores than are found with raw values (p < .001), predictions based on summary z-score measures achieved low predictive power (AUC < 0.63). Instead, we find that combining information from the different white matter tracts, while using multiple imaging measures simultaneously, improves prediction performance (the best predictor achieved AUC = 0.726). Our findings suggest that extreme deviations from a normative model are not optimal features for prediction. However, including the complete distribution of deviations across multiple imaging measures improves prediction, and could aid in subject-level classification.

The Normative Modeling Framework for Computational Psychiatry.. [PMID: 35650452 PMCID: PMC7613648 DOI: 10.1101/2021.08.08.455583]

Normative modeling is an emerging and innovative framework for mapping individual differences at the level of a single subject or observation in relation to a reference model. It involves charting centiles of variation across a population in terms of mappings between biology and behavior which can then be used to make statistical inferences at the level of the individual. The fields of computational psychiatry and clinical neuroscience have been slow to transition away from patient versus “healthy” control analytic approaches, likely due to a lack of tools designed to properly model biological heterogeneity of mental disorders. Normative modeling provides a solution to address this issue and moves analysis away from case-control comparisons that rely on potentially noisy clinical labels. In this article, we define a standardized protocol to guide users through, from start to finish, normative modeling analysis using the Predictive Clinical Neuroscience toolkit (PCNtoolkit). We describe the input data selection process, provide intuition behind the various modeling choices, and conclude by demonstrating several examples of down-stream analyses the normative model results may facilitate, such as stratification of high-risk individuals, subtyping, and behavioral predictive modeling. The protocol takes approximately 1-3 hours to complete.

Replicating extensive brain structural heterogeneity in individuals with schizophrenia and bipolar disorder

Identifying brain processes involved in the risk and development of mental disorders is a major aim. We recently reported substantial interindividual heterogeneity in brain structural aberrations among patients with schizophrenia and bipolar disorder. Estimating the normative range of voxel-based morphometry (VBM) data among healthy individuals using a Gaussian process regression (GPR) enables us to map individual deviations from the healthy range in unseen datasets. Here, we aim to replicate our previous results in two independent samples of patients with schizophrenia (n1 = 94; n2 = 105), bipolar disorder (n1 = 116; n2 = 61), and healthy individuals (n1 = 400; n2 = 312). In line with previous findings with exception of the cerebellum our results revealed robust group level differences between patients and healthy individuals, yet only a small proportion of patients with schizophrenia or bipolar disorder exhibited extreme negative deviations from normality in the same brain regions. These direct replications support that group level-differences in brain structure disguise considerable individual differences in brain aberrations, with important implications for the interpretation and generalization of group-level brain imaging findings to the individual with a mental disorder.