Understanding Heterogeneity in Clinical Cohorts Using Normative Models: Beyond Case-Control Studies

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 used deep autoencoders in an anomaly detection framework, combined for the first time with a confounder removal step that integrates 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 brain deviating patterns differing between the subject and the group levels, a step forward towards the development of more effective and personalized clinical decision support systems and patient stratification in psychiatry.

Smooth Normative Brain Mapping of Three-Dimensional Morphometry Imaging Data Using Skew-Normal Regression

Tensor-based morphometry (TBM) aims at showing local differences in brain volumes with respect to a common template. TBM images are smooth, but they exhibit (especially in diseased groups) higher values in some brain regions called lateral ventricles. More specifically, our voxelwise analysis shows both a mean-variance relationship in these areas and evidence of spatially dependent skewness. We propose a model for three-dimensional imaging data where mean, variance and skewness functions vary smoothly across brain locations. We model the voxelwise distributions as skew-normal. We illustrate an interpolation-based approach to obtain smooth parameter functions based on a subset of voxels. The effects of age and sex are estimated on a reference population of cognitively normal subjects from the Alzheimer's Disease Neuroimaging Initiative (ADNI) data set and mapped across the whole brain. The three parameter functions allow transforming each TBM image (in the reference population as well as in a test set) into a normative map based on Gaussian distributions. These subject-specific normative maps are used to derive indices of deviation from a healthy condition to assess the individual risk of pathological degeneration.

The Future of Disentangling the Heterogeneity of Autism With Neuroimaging Studies

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

Brain-Charting Autism and Attention-Deficit/Hyperactivity Disorder Reveals Distinct and Overlapping Neurobiology

BACKGROUND

Autism and attention-deficit/hyperactivity disorder (ADHD) are heterogeneous neurodevelopmental conditions with complex underlying neurobiology that is still poorly understood. Despite overlapping presentation and sex-biased prevalence, autism and ADHD are rarely studied together and sex differences are often overlooked. Population modeling, often referred to as normative modeling, provides a unified framework for studying age-specific and sex-specific divergences in brain development.

METHODS

Here, we used population modeling and a large, multisite neuroimaging dataset (N = 4255 after quality control) to characterize cortical anatomy associated with autism and ADHD, benchmarked against models of average brain development based on a sample of more than 75,000 individuals. We also examined sex and age differences and relationship with autistic traits and explored the co-occurrence of autism and ADHD.

RESULTS

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

CONCLUSIONS

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

Investigating dopaminergic abnormalities in schizophrenia and first-episode psychosis with normative modelling and multisite molecular neuroimaging

Molecular neuroimaging techniques, like PET and SPECT, offer invaluable insights into the brain's in-vivo biology and its dysfunction in neuropsychiatric patients. However, the transition of molecular neuroimaging into diagnostics and precision medicine has been limited to a few clinical applications, hindered by issues like practical feasibility, high costs, and high between-subject heterogeneity of neuroimaging measures. In this study, we explore the use of normative modelling (NM) to identify individual patient alterations by describing the physiological variability of molecular functions. NM potentially addresses challenges such as small sample sizes and diverse acquisition protocols typical of molecular neuroimaging studies. We applied NM to two PET radiotracers targeting the dopaminergic system ([11C]-(+)-PHNO and [18F]FDOPA) to create a reference-cohort model of healthy controls. The models were subsequently utilized on different independent cohorts of patients with psychosis in different disease stages and treatment outcomes. Our results showed that patients with psychosis exhibited a higher degree of extreme deviations (~3-fold increase) than controls, although this pattern was heterogeneous, with minimal overlap of extreme deviations topology (max 20%). We also confirmed that striatal [18F]FDOPA signal, when referenced to a normative distribution, can predict treatment response (striatal AUC ROC: 0.77-0.83). In conclusion, our results indicate that normative modelling can be effectively applied to molecular neuroimaging after proper harmonization, enabling insights into disease mechanisms and advancing precision medicine. In addition, the method is valuable in understanding the heterogeneity of patient populations and can contribute to maximising cost efficiency in studies aimed at comparing cases and controls.

Neurofind: using deep learning to make individualised inferences in brain-based disorders

Within precision psychiatry, there is a growing interest in normative models given their ability to parse heterogeneity. While they are intuitive and informative, the technical expertise and resources required to develop normative models may not be accessible to most researchers. Here we present Neurofind, a new freely available tool that bridges this gap by wrapping sound and previously tested methods on data harmonisation and advanced normative models into a web-based platform that requires minimal input from the user. We explain how Neurofind was developed, how to use the Neurofind website in four simple steps ( www.neurofind.ai ), and provide exemplar applications. Neurofind takes as input structural MRI images and outputs two main metrics derived from independent normative models: (1) Outlier Index Score, a deviation score from the normative brain morphology, and (2) Brain Age, the predicted age based on an individual's brain morphometry. The tool was trained on 3362 images of healthy controls aged 20-80 from publicly available datasets. The volume of 101 cortical and subcortical regions was extracted and modelled with an adversarial autoencoder for the Outlier index model and a support vector regression for the Brain age model. To illustrate potential applications, we applied Neurofind to 364 images from three independent datasets of patients diagnosed with Alzheimer's disease and schizophrenia. In Alzheimer's disease, 55.2% of patients had very extreme Outlier Index Scores, mostly driven by larger deviations in temporal-limbic structures and ventricles. Patients were also homogeneous in how they deviated from the norm. Conversely, only 30.1% of schizophrenia patients were extreme outliers, due to deviations in the hippocampus and pallidum, and patients tended to be more heterogeneous than controls. Both groups showed signs of accelerated brain ageing.

Beyond case-control study in neuroimaging for psychiatric disorders: Harmonizing and utilizing the brain images from multiple sites

Recent magnetic resonance imaging (MRI) research has advanced our understanding of brain pathophysiology in psychiatric disorders. This progress necessitates re-evaluation of the diagnostic system for psychiatric disorders based on MRI-based biomarkers, with implications for precise clinical diagnosis and optimal therapeutics. To achieve this goal, large-scale multi-site studies are essential to develop a standardized MRI database, with the analysis of several thousands of images and the incorporation of new data. A critical challenge in these studies is to minimize sampling and measurement biases in MRI studies to accurately capture the diversity of disease-derived biomarkers. Various techniques have been employed to consolidate datasets from multiple sites in case-control studies. Traveling subject harmonization stands out as a powerful tool that can differentiate measurement bias from sample variety and sampling bias. A non-linear statistical model for a normative trajectory across the lifespan also strengthens the database to mitigate sampling bias from known factors such as age and sex. These approaches can enhance the alterations between psychiatric disorders and integrate new data and follow-up scans into existing life-course trajectory, enhancing the reliability of machine learning classification and subtyping. Although this approach has been developed using T1-weighted structural image features, future research may extend this framework to other modalities and measures. The required sample size and methodological establishment are needed for future investigations, leading to novel insights into the brain pathophysiology of psychiatric disorders and the development of optimal therapeutics for bedside clinical applications. Sharing big data and their findings also need to be considered.

Estimating Multimodal Structural Brain Variability in Schizophrenia Spectrum Disorders: A Worldwide ENIGMA Study

OBJECTIVE

The clinical diversity of schizophrenia is reflected by structural brain variability. It remains unclear how this variability manifests across different gray and white matter features. In this meta- and mega-analysis, the authors investigated how brain heterogeneity in schizophrenia is distributed across multimodal structural indicators.

METHODS

The authors used the ENIGMA dataset of MRI-based brain measures from 22 international sites with up to 6,037 individuals for a given brain measure. Variability and mean values of cortical thickness, cortical surface area, cortical folding index, subcortical volume, and fractional anisotropy were examined in individuals with schizophrenia and healthy control subjects.

RESULTS

Individuals with schizophrenia showed greater variability in cortical thickness, cortical surface area, subcortical volume, and fractional anisotropy within the frontotemporal and subcortical network. This increased structural variability was mainly associated with psychopathological symptom domains, and the schizophrenia group frequently displayed lower mean values in the respective structural measures. Unexpectedly, folding patterns were more uniform in individuals with schizophrenia, particularly in the right caudal anterior cingulate region. The mean folding values of the right caudal anterior cingulate region did not differ between the schizophrenia and healthy control groups, and folding patterns in this region were not associated with disease-related parameters.

CONCLUSIONS

In patients with schizophrenia, uniform folding patterns in the right caudal anterior cingulate region contrasted with the multimodal variability in the frontotemporal and subcortical network. While variability in the frontotemporal and subcortical network was associated with disease-related diversity, uniform folding may indicate a less flexible interplay between genetic and environmental factors during neurodevelopment.

Judged by your neighbors: Brain structural normativity profiles for large and heterogeneous samples

The detection of norm deviations is fundamental to clinical decision making and impacts our ability to diagnose and treat diseases effectively. Current normative modeling approaches rely on generic comparisons and quantify deviations in relation to the population average. However, generic models interpolate subtle nuances and risk the loss of critical information, thereby compromising effective personalization of health care strategies. To acknowledge the substantial heterogeneity among patients and support the paradigm shift of precision medicine, we introduce Nearest Neighbor Normativity (N 3 ), which is a strategy to refine normativity evaluations in diverse and heterogeneous clinical study populations. We address current methodological shortcomings by accommodating several equally normative population prototypes, comparing individuals from multiple perspectives and designing specifically tailored control groups. Applied to brain structure in 36,896 individuals, the N 3 framework provides empirical evidence for its utility and significantly outperforms traditional methods in the detection of pathological alterations. Our results underscore N 3 's potential for individual assessments in medical practice, where norm deviations are not merely a benchmark, but an important metric supporting the realization of personalized patient care.

Common neuroanatomical differential factors underlying heterogeneous gray matter volume variations in five common psychiatric disorders

Multifaceted evidence has shown that psychiatric disorders share common neurobiological mechanisms. However, the tremendous inter-individual heterogeneity among patients with psychiatric disorders limits trans-diagnostic studies with case-control designs, aimed at identifying clinically promising neuroimaging biomarkers. This study aims to identify neuroanatomical differential factors (ND factors) underlying gray matter volume variations in five psychiatric disorders. We leverage 4 independent datasets of 878 patients diagnosed with psychiatric disorders and 585 healthy controls (HCs) to identify shared ND factors underlying individualized gray matter volume variations. Individualized gray matter volume variations are represented with the linear weighted sum of ND factors, and each case is assigned a unique factor composition, thus preserving interindividual variation. We identify four robust ND factors that can be generalized to unseen disorders. ND factors show significant association with group-level morphological abnormalities, reconciling individual- and group-level morphological abnormalities, and are characterized by dissociable cognitive processes, molecular signatures, and connectome-informed epicenters. Moreover, using factor compositions as features, we discover two robust transdiagnostic subtypes with opposite gray matter volume variations relative to HCs. In conclusion, we identify four reproducible and shared neuroanatomical factors that underlie the highly heterogeneous morphological abnormalities in psychiatric disorders.

Individualized resting-state functional connectivity abnormalities unveil two major depressive disorder subtypes with contrasting abnormal patterns of abnormality

Modern neuroimaging research has recognized that major depressive disorder (MDD) is a connectome disorder, characterized by altered functional connectivity across large-scale brain networks. However, the clinical heterogeneity, likely stemming from diverse neurobiological disturbances, complicates findings from standard group comparison methods. This variability has driven the search for MDD subtypes using objective neuroimaging markers. In this study, we sought to identify potential MDD subtypes from subject-level abnormalities in functional connectivity, leveraging a large multi-site dataset of resting-state MRI from 1276 MDD patients and 1104 matched healthy controls. Subject-level extreme functional connections, determined by comparing against normative ranges derived from healthy controls using tolerance intervals, were used to identify biological subtypes of MDD. We identified a set of extreme functional connections that were predominantly between the visual network and the frontoparietal network, the default mode network and the ventral attention network, with the key regions in the anterior cingulate cortex, bilateral orbitofrontal cortex, and supramarginal gyrus. In MDD patients, these extreme functional connections were linked to age of onset and reward-related processes. Using these features, we identified two subtypes with distinct patterns of functional connectivity abnormalities compared to healthy controls (p < 0.05, Bonferroni correction). When considering all patients together, no significant differences were found. These subtypes significantly enhanced case-control discriminability and showed strong internal discriminability between subtypes. Furthermore, the subtypes were reproducible across varying parameters, study sites, and in untreated patients. Our findings provide new insights into the taxonomy and have potential implications for both diagnosis and treatment of MDD.

Machine learning reveals sex differences in distinguishing between conduct-disordered and neurotypical youth based on emotion processing dysfunction

BACKGROUND

Theoretical models of conduct disorder (CD) highlight that deficits in emotion recognition, learning, and regulation play a pivotal role in CD etiology. With CD being more prevalent in boys than girls, various theories aim to explain this sex difference. The "differential threshold" hypothesis suggests greater emotion dysfunction in conduct-disordered girls than boys, but previous research using conventional statistical analyses has failed to support this hypothesis. Here, we used novel analytic techniques such as machine learning (ML) to uncover potentially sex-specific differences in emotion dysfunction among girls and boys with CD compared to their neurotypical peers.

METHODS

Multi-site data from 542 youth with CD and 710 neurotypical controls (64% girls, 9-18 years) who completed emotion recognition, learning, and regulation tasks were analyzed using a multivariate ML classifier to distinguish between youth with CD and controls separately by sex.

RESULTS

Both female and male ML classifiers accurately predicted (above chance level) individual CD status based solely on the neurocognitive features of emotion dysfunction. Notably, the female classifier outperformed the male classifier in identifying individuals with CD. However, the classification and identification performance of both classifiers was below the clinically relevant 80% accuracy threshold (although they still provided relatively fair and realistic estimates of ~ 60% classification performance), probably due to the substantial neurocognitive heterogeneity within such a large and diverse, multi-site sample of youth with CD (and neurotypical controls).

CONCLUSIONS

These findings confirm the close association between emotion dysfunction and CD in both sexes, with a stronger association observed in affected girls, which aligns with the "differential threshold" hypothesis. However, the data also underscore the heterogeneity of CD, namely that only a subset of those affected are likely to have emotion dysfunction and that other neurocognitive domains (not tested here) probably also contribute to CD etiology.

CLINICAL TRIAL NUMBER

Not applicable.

Predicting individual traits from models of brain dynamics accurately and reliably using the Fisher kernel

Predicting an individual's cognitive traits or clinical condition using brain signals is a central goal in modern neuroscience. This is commonly done using either structural aspects, such as structural connectivity or cortical thickness, or aggregated measures of brain activity that average over time. But these approaches are missing a central aspect of brain function: the unique ways in which an individual's brain activity unfolds over time. One reason why these dynamic patterns are not usually considered is that they have to be described by complex, high-dimensional models; and it is unclear how best to use these models for prediction. We here propose an approach that describes dynamic functional connectivity and amplitude patterns using a Hidden Markov model (HMM) and combines it with the Fisher kernel, which can be used to predict individual traits. The Fisher kernel is constructed from the HMM in a mathematically principled manner, thereby preserving the structure of the underlying model. We show here, in fMRI data, that the HMM-Fisher kernel approach is accurate and reliable. We compare the Fisher kernel to other prediction methods, both time-varying and time-averaged functional connectivity-based models. Our approach leverages information about an individual's time-varying amplitude and functional connectivity for prediction and has broad applications in cognitive neuroscience and personalised medicine.

Individualized gray matter morphological abnormalities unveil two neuroanatomical obsessive-compulsive disorder subtypes

Obsessive-compulsive disorder (OCD) is a highly heterogeneous disorder, with notable variations among cases in structural brain abnormalities. To address this heterogeneity, our study aimed to delineate OCD subtypes based on individualized gray matter morphological differences. We recruited 100 untreated, first-episode OCD patients and 106 healthy controls for structural imaging scans. Utilizing normative models of gray matter volume, we identified subtypes based on individual morphological abnormalities. Sensitivity analyses were conducted to validate the reproducibility of clustering outcomes. To gain deeper insights into the connectomic and molecular underpinnings of structural brain abnormalities in the identified subtypes, we investigated their associations with normal brain network architecture and the distribution of neurotransmitter receptors/transporters. Our findings revealed two distinct OCD subtypes exhibiting divergent patterns of structural brain abnormalities. Sensitivity analysis results confirmed the robustness of the identified subtypes. Subtype 1 displayed significantly increased gray matter volume in regions including the frontal gyrus, precuneus, insula, hippocampus, parahippocampal gyrus, amygdala, and temporal gyrus, while subtype 2 exhibited decreased gray matter volume in the frontal gyrus, precuneus, insula, superior parietal gyrus, temporal gyrus, and fusiform gyrus. When considering all patients collectively, structural brain abnormalities nullified. The identified subtypes were characterized by divergent disease epicenters. Specifically, subtype 1 showed disease epicenters in the middle frontal gyrus, while subtype 2 displayed disease epicenters in the striatum, thalamus and hippocampus. Furthermore, structural brain abnormalities in these subtypes displayed distinct associations with neurotransmitter receptors/transporters. The identified subtypes offer novel insights into nosology and the heterogeneous nature of OCD.

Spectral normative modeling of brain structure

Normative modeling in neuroscience aims to characterize interindividual variation in brain phenotypes and thus establish reference ranges, or brain charts, against which individual brains can be compared. Normative models are typically limited to coarse spatial scales due to computational constraints, limiting their spatial specificity. They additionally depend on fixed regions from fixed parcellation atlases, restricting their adaptability to alternative parcellation schemes. To overcome these key limitations, we propose spectral normative modeling (SNM), which leverages brain eigenmodes for efficient spatial reconstruction to generate normative ranges for arbitrary new regions of interest. Benchmarking against conventional counterparts, SNM achieves a 98.3% speedup in computing accurate normative ranges across spatial scales, from millimeters to the whole brain. We demonstrate its utility by elucidating high-resolution individual cortical atrophy patterns and characterizing the heterogeneous nature of neurodegeneration in Alzheimer's disease. SNM lays the groundwork for a new generation of spatially precise brain charts, offering substantial potential to drive advances in individualized precision medicine.

Rightward brain structural asymmetry in young children with autism

To understand the neural mechanism of autism spectrum disorder (ASD) and developmental delay/intellectual disability (DD/ID) that can be associated with ASD, it is important to investigate individuals at an early stage with brain, behavioural and also genetic measures, but such research is still lacking. Here, using the cross-sectional sMRI data of 1030 children under 8 years old, we employed developmental normative models to investigate the atypical development of gray matter volume (GMV) asymmetry in individuals with ASD without DD/ID, ASD with DD/ID and individuals with only DD/ID, and their associations with behavioral and clinical measures and transcription profiles. By extracting the individual deviations of patients from the typical controls with normative models, we found a commonly abnormal pattern of GMV asymmetry across all ASD children: more rightward laterality in the inferior parietal lobe and precentral gyrus, and higher individual variability in the temporal pole. Specifically, ASD with DD/ID children showed a severer and more extensive abnormal pattern in GMV asymmetry deviation values, which was linked with both ASD symptoms and verbal IQ. The abnormal pattern of ASD without DD/ID children showed higher and more extensive individual variability, which was linked with ASD symptoms only. DD/ID children showed no significant differences from healthy population in asymmetry. Lastly, the GMV laterality patterns of all patient groups were significantly associated with both shared and unique gene expression profiles. Our findings provide evidence for rightward GMV asymmetry of some cortical regions in young ASD children (1-7 years) in a large sample (1030 cases), show that these asymmetries are related to ASD symptoms, and identify genes that are significantly associated with these differences.

A multimodal neural signature of face processing in autism within the fusiform gyrus

Atypical face processing is commonly reported in autism. Its neural correlates have been explored extensively across single neuroimaging modalities within key regions of the face processing network, such as the fusiform gyrus (FFG). Nonetheless, it is poorly understood how variation in brain anatomy and function jointly impacts face processing and social functioning. Here we leveraged a large multimodal sample to study the cross-modal signature of face processing within the FFG across four imaging modalities (structural magnetic resonance imaging (MRI), resting-state functional magnetic resonance imaging, task-functional magnetic resonance imaging and electroencephalography) in 204 autistic and nonautistic individuals aged 7-30 years (case-control design). We combined two methodological innovations-normative modeling and linked independent component analysis-to integrate individual-level deviations across modalities and assessed how multimodal components differentiated groups and informed social functioning in autism. Groups differed significantly in a multimodal component driven by bilateral resting-state functional MRI, bilateral structure, right task-functional MRI and left electroencephalography loadings in face-selective and retinotopic FFG. Multimodal components outperformed unimodal ones in differentiating groups. In autistic individuals, multimodal components were associated with cognitive and clinical features linked to social, but not nonsocial, functioning. These findings underscore the importance of elucidating multimodal neural associations of social functioning in autism, offering potential for the identification of mechanistic and prognostic biomarkers.

The Transition From Homogeneous to Heterogeneous Machine Learning in Neuropsychiatric Research

Despite the advantage of neuroimaging-based machine learning (ML) models as pivotal tools for investigating brain-behavior relationships in neuropsychiatric studies, these data-driven predictive approaches have yet to yield substantial, clinically actionable insights for mental health care. A notable impediment lies in the inadequate accommodation of most ML research to the natural heterogeneity within large samples. Although commonly thought of as individual-level analyses, many ML algorithms are unimodal and homogeneous and thus incapable of capturing the potentially heterogeneous relationships between biology and psychopathology. We review the current landscape of computational research targeting population heterogeneity and argue that there is a need to expand from brain subtyping and behavioral phenotyping to analyses that focus on heterogeneity at the relational level. To this end, we review and suggest several existing ML models with the capacity to discern how external environmental and sociodemographic factors moderate the brain-behavior mapping function in a data-driven fashion. These heterogeneous ML models hold promise for enhancing the discovery of individualized brain-behavior associations and advancing precision psychiatry.

A Multiform Heterogeneity Framework for Alzheimer's Disease Based on Multimodal Neuroimaging

Understanding the heterogeneity of Alzheimer's disease (AD) is crucial for advancing precision medicine specifically tailored to this disorder. Recent research has deepened our understanding of AD heterogeneity; however, translating these insights from bench to bedside via neuroimaging heterogeneity frameworks presents significant challenges. In this review, we systematically revisit prior studies and summarize the existing methodology of data-driven neuroimaging studies for AD heterogeneity. We organized the current methodology into 1) a subtyping clustering strategy for patients with AD, and we also subdivided it into subtyping analysis based on cross-sectional multimodal neuroimaging profiles and the identification of long-term disease progression from short-term datasets; 2) a stratified strategy that integrates neuroimaging measures with biomarkers; and 3) individual-specific abnormal patterns based on the normative model. Then, we evaluated the characteristics of these studies along 2 dimensions: 1) the understanding of pathology and 2) clinical application. We systematically address the limitations, challenges, and future directions of research into AD heterogeneity. Our goal is to enhance the neuroimaging heterogeneity framework for AD, thereby facilitating its transition from bench to bedside.

Individual suicide risk factors with resting-state brain functional connectivity patterns in bipolar disorder patients based on latent Dirichlet allocation model

BACKGROUND

The widespread problem of suicide and its severe burden in bipolar disorder (BD) necessitate the development of objective risk markers, aiming to enhance individual suicide risk prediction in BD.

METHODS

This study recruited 123 BD patients (61 patients with prior suicide attempted history (PSAs), 62 without (NSAs)) and 68 healthy controls (HEs). The Latent Dirichlet Allocation (LDA) model was used to decompose the resting state functional connectivity (RSFC) into multiple hyper/hypo-RSFC patterns. Thereafter, according to the quantitative results of individual heterogeneity over latent factor dimensions, the correlations were analyzed to test prediction ability.

RESULTS

Model constructed without introducing suicide-related labels yielded three latent factors with dissociable hyper/hypo-RSFC patterns. In the subsequent analysis, significant differences in the factor distributions of PSAs and NSAs showed biases on the default-mode network (DMN) hyper-RSFC factor (factor 3) and the salience network (SN) and central executive network (CEN) hyper-RSFC factor (factor 1), indicating predictive value. Correlation analysis of the individuals' expressions with their Nurses' Global Assessment of Suicide Risk (NGASR) revealed factor 3 positively correlated (r = 0.4180, p < 0.0001) and factor 1 negatively correlated (r = - 0.2492, p = 0.0055) with suicide risk. Therefore, it could be speculated that patterns more associated with suicide reflected hyper-connectivity in DMN and hypo-connectivity in SN, CEN.

CONCLUSIONS

This study provided individual suicide-associated risk factors that could reflect the abnormal RSFC patterns, and explored the suicide related brain mechanisms, which is expected to provide supports for clinical decision-making and timely screening and intervention for individuals at high risks of suicide.

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

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

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.

Lifespan Normative Models of White Matter Fractional Anisotropy: Applications to Early Psychosis

This study presents large-scale normative models of white matter (WM) organization across the lifespan, using diffusion MRI data from over 25,000 healthy individuals aged 0-100 years. These models capture lifespan trajectories and inter-individual variation in fractional anisotropy (FA), a marker of white matter integrity. By addressing non-Gaussian data distributions, race, and site effects, the models offer reference baselines across diverse ages, ethnicities, and scanning conditions. We applied these FA models to the HCP Early Psychosis cohort and performed a multivariate analysis to map symptoms onto deviations from multimodal normative models using multi-view sparse canonical correlation analysis (msCCA). Our results reveal extensive white matter heterogeneity in psychosis, which is not captured by group-level analyses, with key regions identified, including the right uncinate fasciculus and thalami. These normative models offer valuable tools for individualized WM deviation identification, improving precision in psychiatric assessments. All models are publicly available for community use.

Neurobiology of attention-deficit hyperactivity disorder: historical challenges and emerging frontiers

Extensive investigations spanning multiple levels of inquiry, from genetic to behavioural studies, have sought to unravel the mechanistic foundations of attention-deficit hyperactivity disorder (ADHD), with the aspiration of developing efficacious treatments for this condition. Despite these efforts, the pathogenesis of ADHD remains elusive. In this Review, we reflect on what has been learned about ADHD while also providing a framework that may serve as a roadmap for future investigations. We emphasize that ADHD is a highly heterogeneous disorder with multiple aetiologies that necessitates a multifactorial dimensional phenotype, rather than a fixed dichotomous conceptualization. We highlight new findings that suggest a more brain-wide, 'global' view of the disorder, rather than the traditional localizationist framework, which asserts that a limited set of brain regions or networks underlie ADHD. Last, we underscore how underpowered studies that have aimed to associate neurobiology with ADHD phenotypes have long precluded the field from making progress. However, a new age of ADHD research with refined phenotypes, advanced methods, creative study designs and adequately powered investigations is beginning to put the field on a good footing. Indeed, the field is at a promising juncture to advance the neurobiological understanding of ADHD and fulfil the promise of clinical utility.

Normative Neuroimaging Library: Designing a Comprehensive and Demographically Diverse Dataset of Healthy Controls to Support Traumatic Brain Injury Diagnostic and Therapeutic Development

The past decade has seen impressive advances in neuroimaging, moving from qualitative to quantitative outputs. Available techniques now allow for the inference of microscopic changes occurring in white and gray matter, along with alterations in physiology and function. These existing and emerging techniques hold the potential of providing unprecedented capabilities in achieving a diagnosis and predicting outcomes for traumatic brain injury (TBI) and a variety of other neurological diseases. To see this promise move from the research lab into clinical care, an understanding is needed of what normal data look like for all age ranges, sex, and other demographic and socioeconomic categories. Clinicians can only use the results of imaging scans to support their decision-making if they know how the results for their patient compare with a normative standard. This potential for utilizing magnetic resonance imaging (MRI) in TBI diagnosis motivated the American College of Radiology and Cohen Veterans Bioscience to create a reference database of healthy individuals with neuroimaging, demographic data, and characterization of psychological functioning and neurocognitive data that will serve as a normative resource for clinicians and researchers for development of diagnostics and therapeutics for TBI and other brain disorders. The goal of this article is to introduce the large, well-curated Normative Neuroimaging Library (NNL) to the research community. NNL consists of data collected from ∼1900 healthy participants. The highlights of NNL are (1) data are collected across a diverse population, including civilians, veterans, and active-duty service members with an age range (18-64 years) not well represented in existing datasets; (2) comprehensive structural and functional neuroimaging acquisition with state-of-the-art sequences (including structural, diffusion, and functional MRI; raw scanner data are preserved, allowing higher quality data to be derived in the future; standardized imaging acquisition protocols across sites reflect sequences and parameters often recommended for use with various neurological and psychiatric conditions, including TBI, post-traumatic stress disorder, stroke, neurodegenerative disorders, and neoplastic disease); and (3) the collection of comprehensive demographic details, medical history, and a broad structured clinical assessment, including cognition and psychological scales, relevant to multiple neurological conditions with functional sequelae. Thus, NNL provides a demographically diverse population of healthy individuals who can serve as a comparison group for brain injury study and clinical samples, providing a strong foundation for precision medicine. Use cases include the creation of imaging-derived phenotypes (IDPs), derivation of reference ranges of imaging measures, and use of IDPs as training samples for artificial intelligence-based biomarker development and for normative modeling to help identify injury-induced changes as outliers for precision diagnosis and targeted therapeutic development. On its release, NNL is poised to support the use of advanced imaging in clinician decision support tools, the validation of imaging biomarkers, and the investigation of brain-behavior anomalies, moving the field toward precision medicine.

Shared differential factors underlying individual spontaneous neural activity abnormalities in major depressive disorder

BACKGROUND

In contemporary neuroimaging studies, it has been observed that patients with major depressive disorder (MDD) exhibit aberrant spontaneous neural activity, commonly quantified through the amplitude of low-frequency fluctuations (ALFF). However, the substantial individual heterogeneity among patients poses a challenge to reaching a unified conclusion.

METHODS

To address this variability, our study adopts a novel framework to parse individualized ALFF abnormalities. We hypothesize that individualized ALFF abnormalities can be portrayed as a unique linear combination of shared differential factors. Our study involved two large multi-center datasets, comprising 2424 patients with MDD and 2183 healthy controls. In patients, individualized ALFF abnormalities were derived through normative modeling and further deconstructed into differential factors using non-negative matrix factorization.

RESULTS

Two positive and two negative factors were identified. These factors were closely linked to clinical characteristics and explained group-level ALFF abnormalities in the two datasets. Moreover, these factors exhibited distinct associations with the distribution of neurotransmitter receptors/transporters, transcriptional profiles of inflammation-related genes, and connectome-informed epicenters, underscoring their neurobiological relevance. Additionally, factor compositions facilitated the identification of four distinct depressive subtypes, each characterized by unique abnormal ALFF patterns and clinical features. Importantly, these findings were successfully replicated in another dataset with different acquisition equipment, protocols, preprocessing strategies, and medication statuses, validating their robustness and generalizability.

CONCLUSIONS

This research identifies shared differential factors underlying individual spontaneous neural activity abnormalities in MDD and contributes novel insights into the heterogeneity of spontaneous neural activity abnormalities in MDD.

Individualized gray matter morphological abnormalities uncover two robust transdiagnostic biotypes

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

Dynamic Real-Time Biosensing Enabled Biorhythm Tracking for Psychiatric Disorders

This review article explores the transformative potential of dynamic, real-time biosensing in biorhythm tracking for psychiatric disorders. Psychiatric diseases, characterized by a complex, heterogeneous, and multifactorial pathophysiology, pose challenges in both diagnosis and treatment. Common denominators in the pathophysiology of psychiatric diseases include disruptions in the stress response, sleep-wake cycle, energy metabolism, and immune response: all of these are characterized by a strong biorhythmic regulation (e.g., circadian), leading to dynamic changes in the levels of biomarkers involved. Technological and practical limitations have hindered the analysis of such dynamic processes to date. The integration of biosensors marks a paradigm shift in psychiatric research. These advanced technologies enable multiplex, non-invasive, and near-continuous analysis of biorhythmic biomarkers in real time, overcoming the constraints of conventional approaches. Focusing on the regulation of the stress response, sleep/wake cycle, energy metabolism, and immune response, biosensing allows for a deeper understanding of the heterogeneous and multifactorial pathophysiology of psychiatric diseases. The potential applications of nanobiosensing in biorhythm tracking, however, extend beyond observation. Continuous monitoring of biomarkers can provide a foundation for personalized medicine in Psychiatry, and allow for the transition from syndromal diagnostic entities to pathophysiology-based psychiatric diagnoses. This evolution promises enhanced disease tracking, early relapse prediction, and tailored disease management and treatment strategies. As non-invasive biosensing continues to advance, its integration into biorhythm tracking holds promise not only to unravel the intricate etiology of psychiatric disorders but also for ushering in a new era of precision medicine, ultimately improving the outcomes and quality of life for individuals grappling with these challenging conditions.

"Pscore": A Novel Percentile-Based Metric to Accurately Assess Individual Deviations in Non-Gaussian Distributions of Quantitative MRI Metrics

BACKGROUND

Quantitative magnetic resonance imaging (MRI) metrics could be used in personalized medicine to assess individuals against normative distributions. Conventional Zscore analysis is inadequate in the presence of non-Gaussian distributions. Therefore, if quantitative MRI metrics deviate from normality, an alternative is needed.

PURPOSE

To confirm non-Gaussianity of diffusion MRI (dMRI) metrics on a publicly available dataset, and to propose a novel percentile-based method, "Pscore" to address this issue.

STUDY TYPE

Retrospective cohort.

POPULATION

Nine hundred and sixty-one healthy young adults (age: 22-35 years, females: 53%) from the Human Connectome Project.

FIELD STRENGTH/SEQUENCE

3-T, spin-echo diffusion echo-planar imaging, T1-weighted: MPRAGE.

ASSESSMENT

The dMRI data were preprocessed using the TORTOISE pipeline. Forty-eight regions of interest (ROIs) from the JHU atlas were redrawn on a study-specific diffusion tensor (DT) template and average values were computed from various DT and mean apparent propagator (MAP) metrics. For each ROI, percentile ranks across participants were computed to generate "Pscores"-which normalized the difference between the median and a participant's value with the corresponding difference between the median and the 5th/95th percentile values.

STATISTICAL TESTS

ROI-wise distributions were assessed using log transformations, Zscore, and the "Pscore" methods. The percentages of extreme values above-95th and below-5th percentile boundaries (PEV>95(%), PEV<5(%)) were also assessed in the overall white matter. Bootstrapping was performed to test the reliability of Pscores in small samples (N = 100) using 100 iterations.

RESULTS

The dMRI metric distributions were systematically non-Gaussian, including positively skewed (eg, mean and radial diffusivity) and negatively skewed (eg, fractional and propagator anisotropy) metrics. This resulted in unbalanced tails in Zscore distributions (PEV>95 ≠ 5%, PEV<5 ≠ 5%) whereas "Pscore" distributions were symmetric and balanced (PEV>95 = PEV<5 = 5%); even for small bootstrapped samples (averagePEV > 95 ¯ = PEV < 5 ¯ = 5 ± 0 % [SD]).

DATA CONCLUSION

The inherent skewness observed for dMRI metrics may preclude the use of conventional Zscore analysis. The proposed "Pscore" method may help estimating individual deviations more accurately in skewed normative data, even from small datasets.

LEVEL OF EVIDENCE

1 TECHNICAL EFFICACY: Stage 1.

Running in the FAMILY: understanding and predicting the intergenerational transmission of mental illness

Over 50% of children with a parent with severe mental illness will develop mental illness by early adulthood. However, intergenerational transmission of risk for mental illness in one's children is insufficiently considered in clinical practice, nor is it sufficiently utilised into diagnostics and care for children of ill parents. This leads to delays in diagnosing young offspring and missed opportunities for protective actions and resilience strengthening. Prior twin, family, and adoption studies suggest that the aetiology of mental illness is governed by a complex interplay of genetic and environmental factors, potentially mediated by changes in epigenetic programming and brain development. However, how these factors ultimately materialise into mental disorders remains unclear. Here, we present the FAMILY consortium, an interdisciplinary, multimodal (e.g., (epi)genetics, neuroimaging, environment, behaviour), multilevel (e.g., individual-level, family-level), and multisite study funded by a European Union Horizon-Staying-Healthy-2021 grant. FAMILY focuses on understanding and prediction of intergenerational transmission of mental illness, using genetically informed causal inference, multimodal normative prediction, and animal modelling. Moreover, FAMILY applies methods from social sciences to map social and ethical consequences of risk prediction to prepare clinical practice for future implementation. FAMILY aims to deliver: (i) new discoveries clarifying the aetiology of mental illness and the process of resilience, thereby providing new targets for prevention and intervention studies; (ii) a risk prediction model within a normative modelling framework to predict who is at risk for developing mental illness; and (iii) insight into social and ethical issues related to risk prediction to inform clinical guidelines.

Comparative analysis of cortical anatomy in male participants with internet gaming disorder or tobacco use disorder: Insights from normative modeling

Background

Research on individual differences in brain structural features of internet gaming disorder (IGD) and established addictions such as tobacco use disorder (TUD) is currently limited. This study utilized normative modeling to analyze the cortical thickness (CT) development patterns of male patients with IGD and TUD, aiming to provide further insights into whether IGD qualifies as an addiction.

Methods

Surface-based brain morphometry (SBM) was used to calculate CT from T1-weighted magnetic resonance imaging data of 804 male participants (665 healthy individuals, 68 IGD and 71 TUD). Gaussian process regression was employed to generate normative models of CT development. Deviation maps were produced to depict deviations of IGD and TUD participants from the typical developmental patterns.

Results

Both addiction groups exhibited widespread cortical thinning, particularly in regions such as the bilateral temporal pole and medial orbitofrontal cortex. The TUD group demonstrated a higher degree of individualization and limited spatial overlap compared to the IGD group. Opposite trends in CT changes were observed between the two groups in the bilateral pericalcarine cortex and pars triangularis.

Conclusions

These findings regarding the similarities and differences between IGD and TUD provide support for the idea that IGD shares common features with substance-related addictions and contribute to a deeper understanding of the neural mechanisms underlying IGD.

Three autism subtypes based on single-subject gray matter network revealed by semi-supervised machine learning

Autism spectrum disorder (ASD) is a heterogeneous, early-onset neurodevelopmental condition characterized by persistent impairments in social interaction and communication. This study aims to delineate ASD subtypes based on individual gray matter brain networks and provide new insights from a graph theory perspective. In this study, we extracted and normalized single-subject gray matter networks and calculated each network's topological properties. The heterogeneity through discriminative analysis (HYDRA) method was utilized to subtype all patients based on network properties. Next, we explored the differences among ASD subtypes in terms of network properties and clinical measures. Our investigation identified three distinct ASD subtypes. In the case-control study, these subtypes exhibited significant differences, particularly in the precentral gyrus, lingual gyrus, and middle frontal gyrus. In the case analysis, significant differences in global and nodal properties were observed between any two subtypes. Clinically, subtype 1 showed lower VIQ and PIQ compared to subtype 3, but exhibited higher scores in ADOS-Communication and ADOS-Total compared to subtype 2. The results highlight the distinct brain network properties and behaviors among different subtypes of male patients with ASD, providing valuable insights into the neural mechanisms underlying ASD heterogeneity.

A perspective on brain-age estimation and its clinical promise

Brain-age estimation has gained increased attention in the neuroscientific community owing to its potential use as a biomarker of brain health. The difference between estimated and chronological age based on neuroimaging data enables a unique perspective on brain development and aging, with multiple open questions still remaining in the brain-age research field. This Perspective presents an overview of current advancements in the field and envisions the future evolution of the brain-age framework before its potential deployment in hospital settings.

Functional network characteristics in anxiety- and mania-based subgroups of bipolar I disorder

BACKGROUND

Bipolar disorder I (BD-I) is a heterogeneous disorder with a high prevalence of comorbid anxiety. The aim of this study was to investigate whether anxiety and mania symptoms define distinct subgroups within BD-I and to explore potential differences in functional network characteristics between these subgroups.

METHODS

Subgroups were identified using scores from clinical anxiety and mania scales. After dimension reduction of these scores, data-driven clustering analysis with cross-validation was employed to reveal the existence of subgroups. Resting-state functional magnetic resonance imaging (rs-fMRI) scans were pre-processed using fMRIPrep. After parcellation and network construction, global and regional graph theoretical measures were calculated per subgroup.

RESULTS

Clustering results revealed that, based on anxiety symptomatology, subjects fell into two distinct subgroups, whereas mania symptoms divided subjects into four unique subgroups. These subgroups varied notably on several symptom scales. Network assortativity was significantly associated with anxiety subgroups. Post-hoc pairwise comparisons did not reveal significant global functional network differences between the anxiety subgroups or between mania subgroups. Regional network differences between clinical subgroups were especially apparent for strength and degree in the temporal and frontal lobes.

LIMITATIONS

Small sample size of some subgroups is a limitation of this study as is the categorical rather than continuous representation of anxiety and mania symptoms.

CONCLUSIONS

BD-I populations may be stratified into robust subgroups based on anxiety and mania symptoms, showing differences in functional network connectivity. Our findings highlight new avenues of research for investigating heterogeneity in psychiatric populations.

Dimensional Neuroimaging Endophenotypes: Neurobiological Representations of Disease Heterogeneity Through Machine Learning

Machine learning has been increasingly used to obtain individualized neuroimaging signatures for disease diagnosis, prognosis, and response to treatment in neuropsychiatric and neurodegenerative disorders. Therefore, it has contributed to a better understanding of disease heterogeneity by identifying disease subtypes with different brain phenotypic measures. In this review, we first present a systematic literature overview of studies using machine learning and multimodal magnetic resonance imaging to unravel disease heterogeneity in various neuropsychiatric and neurodegenerative disorders, including Alzheimer's disease, schizophrenia, major depressive disorder, autism spectrum disorder, and multiple sclerosis, as well as their potential in a transdiagnostic framework, where neuroanatomical and neurobiological commonalities were assessed across diagnostic boundaries. Subsequently, we summarize relevant machine learning methodologies and their clinical interpretability. We discuss the potential clinical implications of the current findings and envision future research avenues. Finally, we discuss an emerging paradigm called dimensional neuroimaging endophenotypes. Dimensional neuroimaging endophenotypes dissects the neurobiological heterogeneity of neuropsychiatric and neurodegenerative disorders into low-dimensional yet informative, quantitative brain phenotypic representations, serving as robust intermediate phenotypes (i.e., endophenotypes), presumably reflecting the interplay of underlying genetic, lifestyle, and environmental processes associated with disease etiology.

The impact of psychosocial adversity on brain and behaviour: an overview of existing knowledge and directions for future research

Environmental experiences play a critical role in shaping the structure and function of the brain. Its plasticity in response to different external stimuli has been the focus of research efforts for decades. In this review, we explore the effects of adversity on brain's structure and function and its implications for brain development, adaptation, and the emergence of mental health disorders. We are focusing on adverse events that emerge from the immediate surroundings of an individual, i.e., microenvironment. They include childhood maltreatment, peer victimisation, social isolation, affective loss, domestic conflict, and poverty. We also take into consideration exposure to environmental toxins. Converging evidence suggests that different types of adversity may share common underlying mechanisms while also exhibiting unique pathways. However, they are often studied in isolation, limiting our understanding of their combined effects and the interconnected nature of their impact. The integration of large, deep-phenotyping datasets and collaborative efforts can provide sufficient power to analyse high dimensional environmental profiles and advance the systematic mapping of neuronal mechanisms. This review provides a background for future research, highlighting the importance of understanding the cumulative impact of various adversities, through data-driven approaches and integrative multimodal analysis techniques.

Anhedonia Across and Beyond the Schizophrenia Spectrum

Anhedonia refers to the diminished ability to experience pleasure, and is a core feature of schizophrenia (SCZ). The neurocognitive and neural correlates of anhedonia remain elusive. Based on several influential theoretical models for negative symptoms, this selective review proposed four important neurocognitive domains, which may unveil the neurobiological mechanisms of anhedonia. The authors critically reviewed the current evidence regarding value representation of reward, prospection, emotion-behavior decoupling, and belief updating in the Chinese setting, covering both behavioral and neuroimaging research. We observed a limited application of the transdiagnostic approach in previous studies on the four domains, and the lack of adequate measures to tap into the expressivity deficit in SCZ. Despite many behavioral paradigms for these four domains utilized both social and non-social stimuli, previous studies seldom focused on the social-versus-non-social differentiation. We further advocated several important directions for future research.

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.

Regional patterns of human cortex development correlate with underlying neurobiology

Human brain morphology undergoes complex changes over the lifespan. Despite recent progress in tracking brain development via normative models, current knowledge of underlying biological mechanisms is highly limited. We demonstrate that human cortical thickness development and aging trajectories unfold along patterns of molecular and cellular brain organization, traceable from population-level to individual developmental trajectories. During childhood and adolescence, cortex-wide spatial distributions of dopaminergic receptors, inhibitory neurons, glial cell populations, and brain-metabolic features explain up to 50% of the variance associated with a lifespan model of regional cortical thickness trajectories. In contrast, modeled cortical thickness change patterns during adulthood are best explained by cholinergic and glutamatergic neurotransmitter receptor and transporter distributions. These relationships are supported by developmental gene expression trajectories and translate to individual longitudinal data from over 8000 adolescents, explaining up to 59% of developmental change at cohort- and 18% at single-subject level. Integrating neurobiological brain atlases with normative modeling and population neuroimaging provides a biologically meaningful path to understand brain development and aging in living humans.

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.

Dataset factors associated with age-related changes in brain structure and function in neurodevelopmental conditions

With brain structure and function undergoing complex changes throughout childhood and adolescence, age is a critical consideration in neuroimaging studies, particularly for those of individuals with neurodevelopmental conditions. However, despite the increasing use of large, consortium-based datasets to examine brain structure and function in neurotypical and neurodivergent populations, it is unclear whether age-related changes are consistent between datasets and whether inconsistencies related to differences in sample characteristics, such as demographics and phenotypic features, exist. To address this, we built models of age-related changes of brain structure (regional cortical thickness and regional surface area; N = 1218) and function (resting-state functional connectivity strength; N = 1254) in two neurodiverse datasets: the Province of Ontario Neurodevelopmental Network and the Healthy Brain Network. We examined whether deviations from these models differed between the datasets, and explored whether these deviations were associated with demographic and clinical variables. We found significant differences between the two datasets for measures of cortical surface area and functional connectivity strength throughout the brain. For regional measures of cortical surface area, the patterns of differences were associated with race/ethnicity, while for functional connectivity strength, positive associations were observed with head motion. Our findings highlight that patterns of age-related changes in the brain may be influenced by demographic and phenotypic characteristics, and thus future studies should consider these when examining or controlling for age effects in analyses.

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.

Consistently lower volumes across thalamus nuclei in very premature-born adults

Lasting thalamus volume reduction after preterm birth is a prominent finding. However, whether thalamic nuclei volumes are affected differentially by preterm birth and whether nuclei aberrations are relevant for cognitive functioning remains unknown. Using T1-weighted MR-images of 83 adults born very preterm (≤ 32 weeks' gestation; VP) and/or with very low body weight (≤ 1,500 g; VLBW) as well as of 92 full-term born (≥ 37 weeks' gestation) controls, we compared thalamic nuclei volumes of six subregions (anterior, lateral, ventral, intralaminar, medial, and pulvinar) across groups at the age of 26 years. To characterize the functional relevance of volume aberrations, cognitive performance was assessed by full-scale intelligence quotient using the Wechsler Adult Intelligence Scale and linked to volume reductions using multiple linear regression analyses. Thalamic volumes were significantly lower across all examined nuclei in VP/VLBW adults compared to controls, suggesting an overall rather than focal impairment. Lower nuclei volumes were linked to higher intensity of neonatal treatment, indicating vulnerability to stress exposure after birth. Furthermore, we found that single results for lateral, medial, and pulvinar nuclei volumes were associated with full-scale intelligence quotient in preterm adults, albeit not surviving correction for multiple hypotheses testing. These findings provide evidence that lower thalamic volume in preterm adults is observable across all subregions rather than focused on single nuclei. Data suggest the same mechanisms of aberrant thalamus development across all nuclei after premature birth.

Regional patterns of human cortex development correlate with underlying neurobiology

Human brain morphology undergoes complex changes over the lifespan. Despite recent progress in tracking brain development via normative models, current knowledge of underlying biological mechanisms is highly limited. We demonstrate that human cortical thickness development and aging trajectories unfold along patterns of molecular and cellular brain organization, traceable from population-level to individual developmental trajectories. During childhood and adolescence, cortex-wide spatial distributions of dopaminergic receptors, inhibitory neurons, glial cell populations, and brain-metabolic features explain up to 50% of variance associated with a lifespan model of regional cortical thickness trajectories. In contrast, modeled cortical thickness change patterns during adulthood are best explained by cholinergic and glutamatergic neurotransmitter receptor and transporter distributions. These relationships are supported by developmental gene expression trajectories and translate to individual longitudinal data from over 8,000 adolescents, explaining up to 59% of developmental change at cohort- and 18% at single-subject level. Integrating neurobiological brain atlases with normative modeling and population neuroimaging provides a biologically meaningful path to understand brain development and aging in living humans.

Nonlinear latent representations of high-dimensional task-fMRI data: Unveiling cognitive and behavioral insights in heterogeneous spatial maps

Finding an interpretable and compact representation of complex neuroimaging data is extremely useful for understanding brain behavioral mapping and hence for explaining the biological underpinnings of mental disorders. However, hand-crafted representations, as well as linear transformations, may inadequately capture the considerable variability across individuals. Here, we implemented a data-driven approach using a three-dimensional autoencoder on two large-scale datasets. This approach provides a latent representation of high-dimensional task-fMRI data which can account for demographic characteristics whilst also being readily interpretable both in the latent space learned by the autoencoder and in the original voxel space. This was achieved by addressing a joint optimization problem that simultaneously reconstructs the data and predicts clinical or demographic variables. We then applied normative modeling to the latent variables to define summary statistics ('latent indices') and establish a multivariate mapping to non-imaging measures. Our model, trained with multi-task fMRI data from the Human Connectome Project (HCP) and UK biobank task-fMRI data, demonstrated high performance in age and sex predictions and successfully captured complex behavioral characteristics while preserving individual variability through a latent representation. Our model also performed competitively with respect to various baseline models including several variants of principal components analysis, independent components analysis and classical regions of interest, both in terms of reconstruction accuracy and strength of association with behavioral variables.

Multiscale heterogeneity of white matter morphometry in psychiatric disorders

Background

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

Methods

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

Results

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

Conclusions

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

Individualized Structural Perturbations on Normative Brain Connectome Restrict Deep Brain Stimulation Outcomes in Parkinson's Disease

BACKGROUND

Patients with Parkinson's disease (PD) respond to deep brain stimulation (DBS) variably. However, how brain substrates restrict DBS outcomes remains unclear.

OBJECTIVE

In this article, we aim to identify prognostic brain signatures for explaining the response variability.

METHODS

We retrospectively investigated a cohort of patients with PD (n = 141) between 2017 and 2022, and defined DBS outcomes as the improvement ratio of clinical motor scores. We used a deviation index to quantify individual perturbations on a reference structural covariance network acquired with preoperative T1-weighted magnetic resonance imaging. The neurobiological perturbations of patients were represented as z scored indices based on the chronological perturbations measured on a group of normal aging adults.

RESULTS

After applying stringent statistical tests (z > 2.5) and correcting for false discoveries (P < 0.01), we found that accelerated deviations mainly affected the prefrontal cortex, motor strip, limbic system, and cerebellum in PD. Particularly, a negative network within the accelerated deviations, expressed as "more preoperative deviations, less postoperative improvements," could predict DBS outcomes (mean absolute error = 0.09, R2 = 0.15). Moreover, a fusion of personal brain predictors and medical responses significantly improved traditional evaluations of DBS outcomes. Notably, the most important brain predictor, a pathway connecting the cognitive unit (prefrontal cortex) and motor control unit (cerebellum and motor strip), partially mediates DBS outcomes with the age at surgery.

CONCLUSIONS

Our findings suggest that individual structural perturbations on the cognitive motor control circuit are critical for modulating DBS outcomes. Interventions toward the circuit have the potential for additional clinical improvements. © 2024 The Author(s). Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

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.