Evidence for embracing normative modeling

COVID-19 lockdown effects on adolescent brain structure suggest accelerated maturation that is more pronounced in females than in males

Adolescence is a period of substantial social-emotional development, accompanied by dramatic changes to brain structure and function. Social isolation due to lockdowns that were imposed because of the COVID-19 pandemic had a detrimental impact on adolescent mental health, with the mental health of females more affected than males. We assessed the impact of the COVID-19 pandemic lockdowns on adolescent brain structure with a focus on sex differences. We collected MRI structural data longitudinally from adolescents prior to and after the pandemic lockdowns. The pre-COVID data were used to create a normative model of cortical thickness change with age during typical adolescent development. Cortical thickness values in the post-COVID data were compared to this normative model. The analysis revealed accelerated cortical thinning in the post-COVID brain, which was more widespread throughout the brain and greater in magnitude in females than in males. When measured in terms of equivalent years of development, the mean acceleration was found to be 4.2 y in females and 1.4 y in males. Accelerated brain maturation as a result of chronic stress or adversity during development has been well documented. These findings suggest that the lifestyle disruptions associated with the COVID-19 pandemic lockdowns caused changes in brain biology and had a more severe impact on the female than the male brain.

Personalizing progressive changes to brain structure in Alzheimer's disease using normative modeling

INTRODUCTION

Neuroanatomical normative modeling captures individual variability in Alzheimer's disease (AD). Here we used normative modeling to track individuals' disease progression in people with mild cognitive impairment (MCI) and patients with AD.

METHODS

Cortical and subcortical normative models were generated using healthy controls (n ≈ 58k). These models were used to calculate regional z scores in 3233 T1-weighted magnetic resonance imaging time-series scans from 1181 participants. Regions with z scores < -1.96 were classified as outliers mapped on the brain and summarized by total outlier count (tOC).

RESULTS

tOC increased in AD and in people with MCI who converted to AD and also correlated with multiple non-imaging markers. Moreover, a higher annual rate of change in tOC increased the risk of progression from MCI to AD. Brain outlier maps identified the hippocampus as having the highest rate of change.

DISCUSSION

Individual patients' atrophy rates can be tracked by using regional outlier maps and tOC.

HIGHLIGHTS

Neuroanatomical normative modeling was applied to serial Alzheimer's disease (AD) magnetic resonance imaging (MRI) data for the first time. Deviation from the norm (outliers) of cortical thickness or brain volume was computed in 3233 scans. The number of brain-structure outliers increased over time in people with AD. Patterns of change in outliers varied markedly between individual patients with AD. People with mild cognitive impairment whose outliers increased over time had a higher risk of progression from AD.

Brain topology underlying executive functions across the lifespan: focus on the default mode network

Introduction

While traditional neuroimaging approaches to the study of executive functions (EFs) have typically employed task-evoked paradigms, resting state studies are gaining popularity as a tool for investigating inter-individual variability in the functional connectome and its relationship to cognitive performance outside of the scanner.

Method

Using resting state functional magnetic resonance imaging data from the Human Connectome Project Lifespan database, the present study capitalized on graph theory to chart cross-sectional variations in the intrinsic functional organization of the frontoparietal (FPN) and the default mode (DMN) networks in 500 healthy individuals (from 10 to 100 years of age), to investigate the neural underpinnings of EFs across the lifespan.

Results

Topological properties of both the FPN and DMN were associated with EF performance but not with a control task of picture naming, providing specificity in support for a tight link between neuro-functional and cognitive-behavioral efficiency within the EF domain. The topological organization of the DMN, however, appeared more sensitive to age-related changes relative to that of the FPN.

Discussion

The DMN matures earlier in life than the FPN and it ıs more susceptible to neurodegenerative changes. Because its activity is stronger in conditions of resting state, the DMN might be easier to measure in noncompliant populations and in those at the extremes of the life-span curve, namely very young or elder participants. Here, we argue that the study of its functional architecture in relation to higher order cognition across the lifespan might, thus, be of greater interest compared with what has been traditionally thought.

Dissecting Heterogeneity in Functional Network Connectivity Aberrations in Antipsychotic Medication-Naïve First Episode Psychosis Patients - A Normative Modeling Study

Importance

While there is a general consensus that functional connectome pathology is a key mechanism underlying psychosis spectrum disorders, the literature is plagued with inconsistencies and translation into clinical practice is non-existent. This is perhaps because group-level findings may not be accurate reflections of pathology at the individual patient level.

Objective

To characterize inter-individual heterogeneity in functional networks and investigate if normative values can be leveraged to identify biologically less heterogeneous subgroups of patients.

Design Setting and Participants

We used data collected in a case-control study conducted at the University of Alabama at Birmingham (UAB). We recruited antipsychotic medication-naïve first-episode psychosis patients from UAB outpatient, inpatient, and emergency room settings.

Main Outcomes and Measures

Individual-level patterns of deviations from a normative reference range in resting-state functional networks using the Yeo-17 atlas for parcellations.

Results

Statistical analyses included 108 medication-naïve first-episode psychosis patients. We found that there is a high level of inter-individual heterogeneity in resting-state network connectivity deviations from the normative reference range. Interestingly 48% of patients did not have any functional connectivity deviations, and no more than 11.1% of patients shared functional deviations between the same regions of interest. In a post hoc analysis, we grouped patients based on deviations into four theoretically possible groups. We discovered that all four groups do exist in our experimental data and showed that subgroups based on deviation profiles were significantly less heterogeneous compared to the overall group (positive deviation group: z= -2.88, p = 0.002; negative deviation group: z= -3.36, p<0.001).

Conclusions and Relevance

Our findings experimentally demonstrate that there is a high level of inter-individual heterogeneity in resting-state network pathology in first-episode psychosis patients which support the idea that group-level findings are not accurate reflections of pathology at the individual level. We also demonstrated that normative functional connectivity deviations may have utility for identifying biologically less heterogeneous subgroups of patients, even though they are not distinguishable clinically. Our findings constitute a significant step towards making precision psychiatry a reality, where patients are selected for treatments based on their individual biological characteristics.

KEY POINTS

Question: How heterogeneous is individual-level resting-state functional network pathology in patients suffering from a first psychotic episode? Can normative reference values in functional network connectivity be leveraged to identify biologically more homogenous subgroups of patients?Findings: We report that functional network pathology is highly heterogeneous, with no more than 11% of patients sharing functional deviations between the same regions of interest.Meaning: Normative modeling is a tool that can map individual neurobiological differences and enables the classification of a clinically heterogenous patient group into subgroups that are neurobiologically less heterogenous.

Normative Modeling of Thalamic Nuclear Volumes and Characterization of Lateralized Volume Alterations in Alzheimer's Disease Versus Schizophrenia

BACKGROUND

Thalamic nuclei facilitate a wide range of complex behaviors, emotions, and cognition and have been implicated in neuropsychiatric disorders including Alzheimer's disease (AD) and schizophrenia. The aim of this work was to establish novel normative models of thalamic nuclear volumes and their laterality indices and investigate their changes in schizophrenia and AD.

METHODS

Volumes of bilateral whole thalami and 10 thalamic nuclei were generated from T1 MRI data using a state-of-the-art novel segmentation method in healthy control subjects (n=2374) and early mild cognitive impairment (MCI, n=211), late MCI (n=113), AD (n=88), and schizophrenia (n=168). Normative models for each nucleus were generated from healthy control subjects while controlling for sex, intracranial volume, and site. Extreme z-score deviations (|z|>1.96) and z-score distributions were compared across phenotypes. Z-scores were associated with clinical descriptors.

RESULTS

Increased infranormal and decreased supranormal z-scores were observed in schizophrenia and AD. Z-score shifts representing reduced volumes were observed in most nuclei in schizophrenia and AD with strong overlap in the bilateral pulvinar, medial dorsal, and centromedian nuclei. Shifts were larger in AD with evidence of a left-sided preference in early MCI while a predilection for right thalamic nuclei was observed in schizophrenia. The right medial dorsal nucleus was associated with disorganized thought and daily auditory verbal hallucinations.

CONCLUSION

In AD, thalamic nuclei are more severely and symmetrically affected while in schizophrenia, the right thalamic nuclei are more affected. We highlight the right medial dorsal nucleus, which may mediate multiple symptoms of schizophrenia and is affected early in the disease course.

Multivariate brain-cognition associations in euthymic bipolar disorder

BACKGROUND

People with bipolar disorder (BD) tend to show widespread cognitive impairment compared to healthy controls. Impairments in processing speed (PS), attention and executive function (EF) may represent 'core' impairments that have a role in wider cognitive dysfunction. Cognitive impairments appear to relate to structural brain abnormalities in BD, but whether core deficits are related to particular brain regions is unclear and much of the research on brain-cognition associations is limited by univariate analysis and small samples.

METHODS

Euthymic BD patients (n = 56) and matched healthy controls (n = 26) underwent T1-weighted MRI scans and completed neuropsychological tests of PS, attention and EF. We utilised public datasets to develop normative models of cortical thickness (n = 5977) to generate robust estimations of cortical abnormalities in patients. Canonical correlation analysis was used to assess multivariate brain-cognition associations in BD, controlling for age, sex and premorbid IQ.

RESULTS

BD showed impairments on tests of PS, attention and EF, and abnormal cortical thickness in several brain regions compared to healthy controls. Impairments in tests of PS and EF were most strongly associated with cortical thickness in the left inferior temporal, right entorhinal and right temporal pole areas.

CONCLUSION

Impairments in PS, attention and EF can be observed in euthymic BD and may be related to abnormal cortical thickness in temporal regions. Future research should continue to leverage normative modelling and multivariate methods to examine complex brain-cognition associations in BD. Future research may benefit from exploring covariance between traditional brain structural morphological metrics such as cortical thickness, cortical volume and surface area.

Identifying epileptogenic abnormality by decomposing intracranial EEG and MEG power spectra

BACKGROUND

Accurate identification of abnormal electroencephalographic (EEG) activity is pivotal for diagnosing and treating epilepsy. Recent studies indicate that decomposing brain activity into periodic (oscillatory) and aperiodic (trend across all frequencies) components can illuminate the drivers of spectral activity changes.

NEW METHODS

We analysed intracranial EEG (iEEG) data from 234 subjects, creating a normative map. This map was compared to a cohort of 63 patients with refractory focal epilepsy under consideration for neurosurgery. The normative map was computed using three approaches: (i) relative complete band power, (ii) relative band power with the aperiodic component removed, and (iii) the aperiodic exponent. Abnormalities were calculated for each approach in the patient cohort. We evaluated the spatial profiles, assessed their ability to localize abnormalities, and replicated the findings using magnetoencephalography (MEG).

RESULTS

Normative maps of relative complete band power and relative periodic band power exhibited similar spatial profiles, while the aperiodic normative map revealed higher exponent values in the temporal lobe. Abnormalities estimated through complete band power effectively distinguished between good and bad outcome patients. Combining periodic and aperiodic abnormalities enhanced performance, like the complete band power approach.

COMPARISON WITH EXISTING METHODS AND CONCLUSIONS

Sparing cerebral tissue with abnormalities in both periodic and aperiodic activity may result in poor surgical outcomes. Both periodic and aperiodic components do not carry sufficient information in isolation. The relative complete band power solution proved to be the most reliable method for this purpose. Future studies could investigate how cerebral location or pathology influences periodic or aperiodic abnormalities.

ComBatLS: A location- and scale-preserving method for multi-site image harmonization

Recent work has leveraged massive datasets and advanced harmonization methods to construct normative models of neuroanatomical features and benchmark individuals' morphology. However, current harmonization tools do not preserve the effects of biological covariates including sex and age on features' variances; this failure may induce error in normative scores, particularly when such factors are distributed unequally across sites. Here, we introduce a new extension of the popular ComBat harmonization method, ComBatLS, that preserves biological variance in features' locations and scales. We use UK Biobank data to show that ComBatLS robustly replicates individuals' normative scores better than other ComBat methods when subjects are assigned to sex-imbalanced synthetic "sites". Additionally, we demonstrate that ComBatLS significantly reduces sex biases in normative scores compared to traditional methods. Finally, we show that ComBatLS successfully harmonizes consortium data collected across over 50 studies. R implementation of ComBatLS is available at https://github.com/andy1764/ComBatFamily.

Multimodal Neuroimaging Summary Scores as Neurobiological Markers of Psychosis

BACKGROUND AND HYPOTHESIS

Structural brain alterations are well-established features of schizophrenia but they do not effectively predict disease/disease risk. Similar to polygenic risk scores in genetics, we integrated multifactorial aspects of brain structure into a summary "Neuroscore" and examined its potential as a marker of disease.

STUDY DESIGN

We extracted measures from T1-weighted scans and diffusion tensor imaging (DTI) models from three studies with schizophrenia and healthy individuals. We calculated individual-level summary scores (Neuroscores) for T1-weighted and DTI measures and a combined score (Multimodal Neuroscore-MM). We assessed each score's ability to differentiate schizophrenia cases from controls and its relationship to clinical symptomatology, intelligence quotient (IQ), and medication dosage. We assessed Neuroscore specificity by performing all analyses in a more inclusive psychosis sample and by using scores generated from MDD effect sizes.

STUDY RESULTS

All Neuroscores significantly differentiated schizophrenia cases from controls (T1 d = 0.56, DTI d = 0.29, MM d = 0.64) to a greater degree than individual brain regions. Higher Neuroscores (ie, increased liability) were associated with lower IQ (T1 β = -0.26, DTI β = -0.15, MM β = -0.30). Higher T1-weighted Neuroscores were associated with higher positive and negative symptom severity (Positive β = 0.21, Negative β = 0.16); Higher Multimodal Neuroscores were associated with higher positive symptom severity (β = 0.30). SZ Neuroscores outperformed MDD Neuroscores in predicting IQ (T1: z = 3.5, q = 0.0007; MM: z = 1.8, q = 0.05).

CONCLUSIONS

Neuroscores are a step toward leveraging widespread structural brain alterations in psychosis to identify robust neurobiological markers of disease. Future studies will assess ways to improve neuroscore calculation, including developing the optimal methods to calculate neuroscores and considering disorder overlap.

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.

Quantifying impairment and disease severity using AI models trained on healthy subjects

Automatic assessment of impairment and disease severity is a key challenge in data-driven medicine. We propose a framework to address this challenge, which leverages AI models trained exclusively on healthy individuals. The COnfidence-Based chaRacterization of Anomalies (COBRA) score exploits the decrease in confidence of these models when presented with impaired or diseased patients to quantify their deviation from the healthy population. We applied the COBRA score to address a key limitation of current clinical evaluation of upper-body impairment in stroke patients. The gold-standard Fugl-Meyer Assessment (FMA) requires in-person administration by a trained assessor for 30-45 minutes, which restricts monitoring frequency and precludes physicians from adapting rehabilitation protocols to the progress of each patient. The COBRA score, computed automatically in under one minute, is shown to be strongly correlated with the FMA on an independent test cohort for two different data modalities: wearable sensors (ρ = 0.814, 95% CI [0.700,0.888]) and video (ρ = 0.736, 95% C.I [0.584, 0.838]). To demonstrate the generalizability of the approach to other conditions, the COBRA score was also applied to quantify severity of knee osteoarthritis from magnetic-resonance imaging scans, again achieving significant correlation with an independent clinical assessment (ρ = 0.644, 95% C.I [0.585,0.696]).

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.

Applications of MRI in Schizophrenia: Current Progress in Establishing Clinical Utility

Schizophrenia is a severe mental illness that significantly impacts the lives of affected individuals and with increasing mortality rates. Early detection and intervention are crucial for improving outcomes but the lack of validated biomarkers poses great challenges in such efforts. The use of magnetic resonance imaging (MRI) in schizophrenia enables the investigation of the disorder's etiological and neuropathological substrates in vivo. After decades of research, promising findings of MRI have been shown to aid in screening high-risk individuals and predicting illness onset, and predicting symptoms and treatment outcomes of schizophrenia. The integration of machine learning and deep learning techniques makes it possible to develop intelligent diagnostic and prognostic tools with extracted or selected imaging features. In this review, we aimed to provide an overview of current progress and prospects in establishing clinical utility of MRI in schizophrenia. We first provided an overview of MRI findings of brain abnormalities that might underpin the symptoms or treatment response process in schizophrenia patients. Then, we summarized the ongoing efforts in the computer-aided utility of MRI in schizophrenia and discussed the gap between MRI research findings and real-world applications. Finally, promising pathways to promote clinical translation were provided. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY: Stage 3.

Spatiotemporal patterns in cortical development: Age, puberty, and individual variability from 9 to 13 years of age

Humans and nonhuman primate studies suggest that timing and tempo of cortical development varies neuroanatomically along a sensorimotor-to-association (S-A) axis. Prior human studies have reported a principal S-A axis across various modalities, but largely rely on cross-sectional samples with wide age-ranges. Here, we investigate developmental changes and individual variability in cortical organization along the S-A axis between the ages of 9-13 years using a large, longitudinal sample (N = 2487-3747, 46-50% female) from the Adolescent Brain Cognitive Development Study (ABCD Study®). This work assesses multiple aspects of neurodevelopment indexed by changes in cortical thickness, cortical microarchitecture, and resting-state functional fluctuations. First, we evaluated S-A organization in age-related changes and, then, computed individual-level S-A alignment in brain changes and assessing differences therein due to age, sex, and puberty. Varying degrees of linear and quadratic age-related brain changes were identified along the S-A axis. Yet, these patterns of cortical development were overshadowed by considerable individual variability in S-A alignment. Even within individuals, there was little correspondence between S-A patterning across the different aspects of neurodevelopment investigated (i.e., cortical morphology, microarchitecture, function). Some of the individual variation in developmental patterning of cortical morphology and microarchitecture was explained by age, sex, and pubertal development. Altogether, this work contextualizes prior findings that regional age differences do progress along an S-A axis at a group level, while highlighting broad variation in developmental change between individuals and between aspects of cortical development, in part due to sex and puberty.

Significance Statement

Understanding normative patterns of adolescent brain change, and individual variability therein, is crucial for disentangling healthy and abnormal development. We used longitudinal human neuroimaging data to study several aspects of neurodevelopment during early adolescence and assessed their organization along a sensorimotor-to-association (S-A) axis across the cerebral cortex. Age differences in brain changes were linear and curvilinear along this S-A axis. However, individual-level sensorimotor-association alignment varied considerably, driven in part by differences in age, sex, and pubertal development.

Amygdala Volume is Associated with ADHD Risk and Severity Beyond Comorbidities in Adolescents: Clinical Testing of Brain Chart Reference Standards

Understanding atypicalities in ADHD brain correlates is a step towards better understanding ADHD etiology. Efforts to map atypicalities at the level of brain structure have been hindered by the absence of normative reference standards. Recent publication of brain charts allows for assessment of individual variation relative to age- and sex-adjusted reference standards and thus estimation not only of case-control differences but also of intraindividual prediction.

METHODS

Aim was to examine, whether brain charts can be applied in a sample of adolescents (N = 140, 38% female) to determine whether atypical brain subcortical and total volumes are associated with ADHD at-risk status and severity of parent-rated symptoms, accounting for self-rated anxiety and depression, and parent-rated oppositional defiant disorder (ODD) as well as motion.

RESULTS

Smaller bilateral amygdala volume was associated with ADHD at-risk status, beyond effects of comorbidities and motion, and smaller bilateral amygdala volume was associated with inattention and hyperactivity/impulsivity, beyond effects of comorbidities except for ODD symptoms, and motion.

CONCLUSIONS

Individual differences in amygdala volume meaningfully add to estimating ADHD risk and severity. Conceptually, amygdalar involvement is consistent with behavioral and functional imaging data on atypical reinforcement sensitivity as a marker of ADHD-related risk. Methodologically, results show that brain chart reference standards can be applied to address clinically informative, focused and specific questions.

Age-atypical brain functional networks in autism spectrum disorder: a normative modeling approach

BACKGROUND

Despite extensive research into the neural basis of autism spectrum disorder (ASD), the presence of substantial biological and clinical heterogeneity among diagnosed individuals remains a major barrier. Commonly used case‒control designs assume homogeneity among subjects, which limits their ability to identify biological heterogeneity, while normative modeling pinpoints deviations from typical functional network development at individual level.

METHODS

Using a world-wide multi-site database known as Autism Brain Imaging Data Exchange, we analyzed individuals with ASD and typically developed (TD) controls (total n = 1218) aged 5-40 years, generating individualized whole-brain network functional connectivity (FC) maps of age-related atypicality in ASD. We then used local polynomial regression to estimate a networkwise normative model of development and explored correlations between ASD symptoms and brain networks.

RESULTS

We identified a subset exhibiting highly atypical individual-level FC, exceeding 2 standard deviation from the normative value. We also identified clinically relevant networks (mainly default mode network) at cohort level, since the outlier rates decreased with age in TD participants, but increased in those with autism. Moreover, deviations were linked to severity of repetitive behaviors and social communication symptoms.

CONCLUSIONS

Individuals with ASD exhibit distinct, highly individualized trajectories of brain functional network development. In addition, distinct developmental trajectories were observed among ASD and TD individuals, suggesting that it may be challenging to identify true differences in network characteristics by comparing young children with ASD to their TD peers. This study enhances understanding of the biological heterogeneity of the disorder and can inform precision medicine.

Population-level normative models reveal race- and socioeconomic-related variability in cortical thickness of threat neurocircuitry

The inequitable distribution of economic resources and exposure to adversity between racial groups contributes to mental health disparities within the United States. Consideration of the potential neurodevelopmental consequences, however, has been limited particularly for neurocircuitry known to regulate the emotional response to threat. Characterizing the consequences of inequity on threat neurocircuitry is critical for robust and generalizable neurobiological models of psychiatric illness. Here we use data from the Adolescent Brain and Cognitive Development Study 4.0 release to investigate the contributions of individual and neighborhood-level economic resources and exposure to discrimination. We investigate the potential appearance of race-related differences using both standard methods and through population-level normative modeling. We show that, in a sample of white and Black adolescents, racial inequities in socioeconomic factors largely contribute to the appearance of race-related differences in cortical thickness of threat neurocircuitry. The race-related differences are preserved through the use of population-level models and such models also preserve associations between cortical thickness and specific socioeconomic factors. The present findings highlight that such socioeconomic inequities largely underlie race-related differences in brain morphology. The present findings provide important new insight for the generation of generalizable neurobiological models of psychiatric illness.

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.

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.

Initiating PeriCBD to probe perinatal influences on neurodevelopment during 3-10 years in China

Adverse perinatal factors can interfere with the normal development of the brain, potentially resulting in long-term effects on the comprehensive development of children. Presently, the understanding of cognitive and neurodevelopmental processes under conditions of adverse perinatal factors is substantially limited. There is a critical need for an open resource that integrates various perinatal factors with the development of the brain and mental health to facilitate a deeper understanding of these developmental trajectories. In this Data Descriptor, we introduce a multicenter database containing information on perinatal factors that can potentially influence children's brain-mind development, namely, periCBD, that combines neuroimaging and behavioural phenotypes with perinatal factors at county/region/central district hospitals. PeriCBD was designed to establish a platform for the investigation of individual differences in brain-mind development associated with perinatal factors among children aged 3-10 years. Ultimately, our goal is to help understand how different adverse perinatal factors specifically impact cognitive development and neurodevelopment. Herein, we provide a systematic overview of the data acquisition/cleaning/quality control/sharing, processes of periCBD.

Normative modeling of thalamic nuclear volumes

Thalamic nuclei have been implicated in neurodegenerative and neuropsychiatric disorders. Normative models for thalamic nuclear volumes have not been proposed thus far. The aim of this work was to establish normative models of thalamic nuclear volumes and subsequently investigate changes in thalamic nuclei in cognitive and psychiatric disorders. Volumes of the bilateral thalami and 12 nuclear regions were generated from T1 MRI data using a novel segmentation method (HIPS-THOMAS) in healthy control subjects (n=2374) and non-control subjects (n=695) with early and late mild cognitive impairment (EMCI, LMCI), Alzheimer's disease (AD), Early psychosis and Schizophrenia, Bipolar disorder, and Attention deficit hyperactivity disorder. Three different normative modelling methods were evaluated while controlling for sex, intracranial volume, and site. Z-scores and extreme z-score deviations were calculated and compared across phenotypes. GAMLSS models performed the best. Statistically significant shifts in z-score distributions consistent with atrophy were observed for most phenotypes. Shifts of progressively increasing magnitude were observed bilaterally from EMCI to AD with larger shifts in the left thalamic regions. Heterogeneous shifts were observed in psychiatric diagnoses with a predilection for the right thalamic regions. Here we present the first normative models of thalamic nuclear volumes and highlight their utility in evaluating heterogenous disorders such as Schizophrenia.

Treatment Response Prediction in Major Depressive Disorder Using Multimodal MRI and Clinical Data: Secondary Analysis of a Randomized Clinical Trial

OBJECTIVE

Response to antidepressant treatment in major depressive disorder varies substantially between individuals, which lengthens the process of finding effective treatment. The authors sought to determine whether a multimodal machine learning approach could predict early sertraline response in patients with major depressive disorder. They assessed the predictive contribution of MR neuroimaging and clinical assessments at baseline and after 1 week of treatment.

METHODS

This was a preregistered secondary analysis of data from the Establishing Moderators and Biosignatures of Antidepressant Response in Clinical Care (EMBARC) study, a multisite double-blind, placebo-controlled randomized clinical trial that included 296 adult outpatients with unmedicated recurrent or chronic major depressive disorder. MR neuroimaging and clinical data were collected before and after 1 week of treatment. Performance in predicting response and remission, collected after 8 weeks, was quantified using balanced accuracy (bAcc) and area under the receiver operating characteristic curve (AUROC) scores.

RESULTS

A total of 229 patients were included in the analyses (mean age, 38 years [SD=13]; 66% female). Internal cross-validation performance in predicting response to sertraline (bAcc=68% [SD=10], AUROC=0.73 [SD=0.03]) was significantly better than chance. External cross-validation on data from placebo nonresponders (bAcc=62%, AUROC=0.66) and placebo nonresponders who were switched to sertraline (bAcc=65%, AUROC=0.68) resulted in differences that suggest specificity for sertraline treatment compared with placebo treatment. Finally, multimodal models outperformed unimodal models.

CONCLUSIONS

The study results confirm that early sertraline treatment response can be predicted; that the models are sertraline specific compared with placebo; that prediction benefits from integrating multimodal MRI data with clinical data; and that perfusion imaging contributes most to these predictions. Using this approach, a lean and effective protocol could individualize sertraline treatment planning to improve psychiatric care.

Deep Normative Tractometry for Identifying Joint White Matter Macro- and Micro-structural Abnormalities in Alzheimer's Disease

This study introduces the Deep Normative Tractometry (DNT) framework, that encodes the joint distribution of both macrostructural and microstructural profiles of the brain white matter tracts through a variational autoencoder (VAE). By training on data from healthy controls, DNT learns the normative distribution of tract data, and can delineate along-tract micro-and macro-structural abnormalities. Leveraging a large sample size via generative pre-training, we assess DNT's generalizability using transfer learning on data from an independent cohort acquired in India. Our findings demonstrate DNT's capacity to detect widespread diffusivity abnormalities along tracts in mild cognitive impairment and Alzheimer's disease, aligning closely with results from the Bundle Analytics (BUAN) tractometry pipeline. By incorporating tract geometry information, DNT may be able to distinguish disease-related abnormalities in anisotropy from tract macrostructure, and shows promise in enhancing fine-scale mapping and detection of white matter alterations in neurodegenerative conditions.

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.

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.

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.

Unpacking the functional heterogeneity of the Emotional Face Matching Task: a normative modelling approach

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.

Effects of anterior temporal lobe resection on cortical morphology

Neuroimaging can capture brain restructuring after anterior temporal lobe resection (ATLR), a surgical procedure to treat drug-resistant temporal lobe epilepsy (TLE). Here, we examine the effects of this surgery on brain morphology measured in recently-proposed independent variables. We studied 101 individuals with TLE (55 left, 46 right onset) who underwent ATLR. For each individual we considered one pre-surgical MRI and one follow-up MRI 2-13 months after surgery. We used a surface-based method to locally compute traditional morphological variables, and the independent measures K, I, and S, where K measures white matter tension, I captures isometric scaling, and S contains the remaining information about cortical shape. A normative model trained on data from 924 healthy controls was used to debias the data and account for healthy ageing effects occurring during scans. A SurfStat random field theory clustering approach assessed changes across the cortex caused by ATLR. Compared to preoperative data, surgery had marked effects on all morphological measures. Ipsilateral effects were located in the orbitofrontal and inferior frontal gyri, the pre- and postcentral gyri and supramarginal gyrus, and the lateral occipital gyrus and lingual cortex. Contralateral effects were in the lateral occipital gyrus, and inferior frontal gyrus and frontal pole. The restructuring following ATLR is reflected in widespread morphological changes, mainly in regions near the resection, but also remotely in regions that are structurally connected to the anterior temporal lobe. The causes could include mechanical effects, Wallerian degeneration, or compensatory plasticity. The study of independent measures revealed additional effects compared to traditional measures.

Personalising Alzheimer's Disease progression using brain atrophy markers

INTRODUCTION

Neuroanatomical normative modelling can capture individual variability in Alzheimer's Disease (AD). We used neuroanatomical normative modelling to track individuals' disease progression in people with mild cognitive impairment (MCI) and patients with AD.

METHODS

Cortical thickness and subcortical volume neuroanatomical normative models were generated using healthy controls (n~58k). These models were used to calculate regional Z-scores in 4361 T1-weighted MRI time-series scans. Regions with Z-scores <-1.96 were classified as outliers and mapped on the brain, and also summarised by total outlier count (tOC).

RESULTS

Rate of change in tOC increased in AD and in people with MCI who converted to AD and correlated with multiple non-imaging markers. Moreover, a higher annual rate of change in tOC increased the risk of MCI progression to AD. Brain Z-score maps showed that the hippocampus had the highest rate of atrophy change.

CONCLUSIONS

Individual-level atrophy rates can be tracked by using regional outlier maps and tOC.