Multi-scale semi-supervised clustering of brain images: Deriving disease subtypes

Analyzing heterogeneity in Alzheimer Disease using multimodal normative modeling on imaging-based ATN biomarkers

INTRODUCTION

Previous studies have applied normative modeling on a single neuroimaging modality to investigate Alzheimer Disease (AD) heterogeneity. We employed a deep learning-based multimodal normative framework to analyze individual-level variation across ATN (amyloid-tau-neurodegeneration) imaging biomarkers.

METHODS

We selected cross-sectional discovery (n = 665) and replication cohorts (n = 430) with available T1-weighted MRI, amyloid and tau PET. Normative modeling estimated individual-level abnormal deviations in amyloid-positive individuals compared to amyloid-negative controls. Regional abnormality patterns were mapped at different clinical group levels to assess intra-group heterogeneity. An individual-level disease severity index (DSI) was calculated using both the spatial extent and magnitude of abnormal deviations across ATN.

RESULTS

Greater intra-group heterogeneity in ATN abnormality patterns was observed in more severe clinical stages of AD. Higher DSI was associated with worse cognitive function and increased risk of disease progression.

DISCUSSION

Subject-specific abnormality maps across ATN reveal the heterogeneous impact of AD on the brain.

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.

Identifying Progression-Specific Alzheimer's Subtypes Using Multimodal Transformer

Alzheimer's disease (AD) is the most prevalent neurodegenerative disease, yet its current treatments are limited to stopping disease progression. Moreover, the effectiveness of these treatments remains uncertain due to the heterogeneity of the disease. Therefore, it is essential to identify disease subtypes at a very early stage. Current data-driven approaches can be used to classify subtypes during later stages of AD or related disorders, but making predictions in the asymptomatic or prodromal stage is challenging. Furthermore, the classifications of most existing models lack explainability, and these models rely solely on a single modality for assessment, limiting the scope of their analysis. Thus, we propose a multimodal framework that utilizes early-stage indicators, including imaging, genetics, and clinical assessments, to classify AD patients into progression-specific subtypes at an early stage. In our framework, we introduce a tri-modal co-attention mechanism (Tri-COAT) to explicitly capture cross-modal feature associations. Data from the Alzheimer's Disease Neuroimaging Initiative (ADNI) (slow progressing = 177, intermediate = 302, and fast = 15) were used to train and evaluate Tri-COAT using a 10-fold stratified cross-testing approach. Our proposed model outperforms baseline models and sheds light on essential associations across multimodal features supported by known biological mechanisms. The multimodal design behind Tri-COAT allows it to achieve the highest classification area under the receiver operating characteristic curve while simultaneously providing interpretability to the model predictions through the co-attention mechanism.

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

BACKGROUND

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

METHODS

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

RESULTS

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

LIMITATIONS

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

CONCLUSIONS

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

A review of neuroimaging-based data-driven approach for Alzheimer's disease heterogeneity analysis

Alzheimer's disease (AD) is a complex form of dementia and due to its high phenotypic variability, its diagnosis and monitoring can be quite challenging. Biomarkers play a crucial role in AD diagnosis and monitoring, but interpreting these biomarkers can be problematic due to their spatial and temporal heterogeneity. Therefore, researchers are increasingly turning to imaging-based biomarkers that employ data-driven computational approaches to examine the heterogeneity of AD. In this comprehensive review article, we aim to provide health professionals with a comprehensive view of past applications of data-driven computational approaches in studying AD heterogeneity and planning future research directions. We first define and offer basic insights into different categories of heterogeneity analysis, including spatial heterogeneity, temporal heterogeneity, and spatial-temporal heterogeneity. Then, we scrutinize 22 articles relating to spatial heterogeneity, 14 articles relating to temporal heterogeneity, and five articles relating to spatial-temporal heterogeneity, highlighting the strengths and limitations of these strategies. Furthermore, we discuss the importance of understanding spatial heterogeneity in AD subtypes and their clinical manifestations, biomarkers for abnormal orderings and AD stages, the recent advancements in spatial-temporal heterogeneity analysis for AD, and the emerging role of omics data integration in advancing personalized diagnosis and treatment for AD patients. By emphasizing the significance of understanding AD heterogeneity, we hope to stimulate further research in this field to facilitate the development of personalized interventions for AD patients.

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 that present significant differences in various brain phenotypic measures. In this review, we first present a systematic literature overview of studies using machine learning and multimodal MRI to unravel disease heterogeneity in various neuropsychiatric and neurodegenerative disorders, including Alzheimer's disease, schizophrenia, major depressive disorder, autism spectrum disorder, multiple sclerosis, as well as their potential in transdiagnostic settings. Subsequently, we summarize relevant machine learning methodologies and discuss an emerging paradigm which we call dimensional neuroimaging endophenotype (DNE). DNE dissects the neurobiological heterogeneity of neuropsychiatric and neurodegenerative disorders into a low-dimensional yet informative, quantitative brain phenotypic representation, serving as a robust intermediate phenotype (i.e., endophenotype) largely reflecting underlying genetics and etiology. Finally, we discuss the potential clinical implications of the current findings and envision future research avenues.

Five dominant dimensions of brain aging are identified via deep learning: associations with clinical, lifestyle, and genetic measures

Brain aging is a complex process influenced by various lifestyle, environmental, and genetic factors, as well as by age-related and often co-existing pathologies. MRI and, more recently, AI methods have been instrumental in understanding the neuroanatomical changes that occur during aging in large and diverse populations. However, the multiplicity and mutual overlap of both pathologic processes and affected brain regions make it difficult to precisely characterize the underlying neurodegenerative profile of an individual from an MRI scan. Herein, we leverage a state-of-the art deep representation learning method, Surreal-GAN, and present both methodological advances and extensive experimental results that allow us to elucidate the heterogeneity of brain aging in a large and diverse cohort of 49,482 individuals from 11 studies. Five dominant patterns of neurodegeneration were identified and quantified for each individual by their respective (herein referred to as) R-indices. Significant associations between R-indices and distinct biomedical, lifestyle, and genetic factors provide insights into the etiology of observed variances. Furthermore, baseline R-indices showed predictive value for disease progression and mortality. These five R-indices contribute to MRI-based precision diagnostics, prognostication, and may inform stratification into clinical trials.

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

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

Functional Imaging Derived ADHD Biotypes Based on Deep Clustering May Guide Personalized Medication Therapy

Attention deficit hyperactivity disorder (ADHD) is one prevalent neurodevelopmental disorder with childhood onset, however, there is no clear correspondence established between clinical ADHD subtypes and primary medications. Identifying objective and reliable neuroimaging markers for categorizing ADHD biotypes may lead to more individualized, biotype-guided treatment. Here we proposed graph convolutional network plus deep clustering for ADHD biotype detection using functional network connectivity (FNC), resulting in two biotypes based on 1069 ADHD patients selected from Adolescent Brain and Cognitive Development (ABCD) study, which were well replicated on independent ADHD adolescents undergoing longitudinal medication treatment (n=130). Interestingly, in addition to differences in cognitive performance and hyperactivity/impulsivity symptoms, biotype 1 treated with methylphenidate demonstrated significantly better recovery than biotype 2 treated with atomoxetine (p<0.05, FDR corrected). This imaging-driven, biotype-guided approach holds promise for facilitating personalized treatment of ADHD, exploring possible boundaries through innovative deep learning algorithms aimed at improving medication treatment effectiveness.

Dissociation between neuroanatomical and symptomatic subtypes in schizophrenia

BACKGROUND

Schizophrenia is a complex and heterogeneous syndrome with high clinical and biological stratification. Identifying distinctive subtypes can improve diagnostic accuracy and help precise therapy. A key challenge for schizophrenia subtyping is understanding the subtype-specific biological underpinnings of clinical heterogeneity. This study aimed to investigate if the machine learning (ML)-based neuroanatomical and symptomatic subtypes of schizophrenia are associated.

METHODS

A total of 314 schizophrenia patients and 257 healthy controls from four sites were recruited. Gray matter volume (GMV) and Positive and Negative Syndrome Scale (PANSS) scores were employed to recognize schizophrenia neuroanatomical and symptomatic subtypes using K-means and hierarchical methods, respectively.

RESULTS

Patients with ML-based neuroanatomical subtype-1 had focally increased GMV, and subtype-2 had widespread reduced GMV than the healthy controls based on either K-means or Hierarchical methods. In contrast, patients with symptomatic subtype-1 had severe PANSS scores than subtype-2. No differences in PANSS scores were shown between the two neuroanatomical subtypes; similarly, no GMV differences were found between the two symptomatic subtypes. Cohen's Kappa test further demonstrated an apparent dissociation between the ML-based neuroanatomical and symptomatic subtypes (P > 0.05). The dissociation patterns were validated in four independent sites with diverse disease progressions (chronic vs. first episodes) and ancestors (Chinese vs. Western).

CONCLUSIONS

These findings revealed a replicable dissociation between ML-based neuroanatomical and symptomatic subtypes of schizophrenia, which provides a new viewpoint toward understanding the heterogeneity of schizophrenia.

Neuroimaging-AI Endophenotypes of Brain Diseases in the General Population: Towards a Dimensional System of Vulnerability

Disease heterogeneity poses a significant challenge for precision diagnostics in both clinical and sub-clinical stages. Recent work leveraging artificial intelligence (AI) has offered promise to dissect this heterogeneity by identifying complex intermediate phenotypes - herein called dimensional neuroimaging endophenotypes (DNEs) - which subtype various neurologic and neuropsychiatric diseases. We investigate the presence of nine such DNEs derived from independent yet harmonized studies on Alzheimer's disease (AD1-2)1, autism spectrum disorder (ASD1-3)2, late-life depression (LLD1-2)3, and schizophrenia (SCZ1-2)4, in the general population of 39,178 participants in the UK Biobank study. Phenome-wide associations revealed prominent associations between the nine DNEs and phenotypes related to the brain and other human organ systems. This phenotypic landscape aligns with the SNP-phenotype genome-wide associations, revealing 31 genomic loci associated with the nine DNEs (Bonferroni corrected P-value < 5×10-8/9). The DNEs exhibited significant genetic correlations, colocalization, and causal relationships with multiple human organ systems and chronic diseases. A causal effect (odds ratio=1.25 [1.11, 1.40], P-value=8.72×1-4) was established from AD2, characterized by focal medial temporal lobe atrophy, to AD. The nine DNEs and their polygenic risk scores significantly improved the prediction accuracy for 14 systemic disease categories and mortality. These findings underscore the potential of the nine DNEs to identify individuals at a high risk of developing the four brain diseases during preclinical stages for precision diagnostics. All results are publicly available at: http://labs.loni.usc.edu/medicine/.

An Alzheimer's disease category progression sub-grouping analysis using manifold learning on ADNI

Many current statistical and machine learning methods have been used to explore Alzheimer's disease (AD) and its associated patterns that contribute to the disease. However, there has been limited success in understanding the relationship between cognitive tests, biomarker data, and patient AD category progressions. In this work, we perform exploratory data analysis of AD health record data by analyzing various learned lower dimensional manifolds to separate early-stage AD categories further. Specifically, we used Spectral embedding, Multidimensional scaling, Isomap, t-Distributed Stochastic Neighbour Embedding, Uniform Manifold Approximation and Projection, and sparse denoising autoencoder based manifolds on the Alzheimer's Disease Neuroimaging Initiative (ADNI) dataset. We then determine the clustering potential of the learned embeddings and then determine if category sub-groupings or sub-categories can be found. We then used a Kruskal-sWallis H test to determine the statistical significance of the discovered AD subcategories. Our results show that the existing AD categories do exhibit sub-groupings, especially in mild cognitive impairment transitions in many of the tested manifolds, showing there may be a need for further subcategories to describe AD progression.

Revealing Individual Neuroanatomical Heterogeneity in Alzheimer Disease Using Neuroanatomical Normative Modeling

BACKGROUND AND OBJECTIVES

Alzheimer disease (AD) is highly heterogeneous, with marked individual differences in clinical presentation and neurobiology. To explore this, we used neuroanatomical normative modeling to index regional patterns of variability in cortical thickness. We aimed to characterize individual differences and outliers in cortical thickness in patients with AD, people with mild cognitive impairment (MCI), and controls. Furthermore, we assessed the relationships between cortical thickness heterogeneity and cognitive function, β-amyloid, phosphorylated-tau, and ApoE genotype. Finally, we examined whether cortical thickness heterogeneity was predictive of conversion from MCI to AD.

METHODS

Cortical thickness measurements across 148 brain regions were obtained from T1-weighted MRI scans from 62 sites of the Alzheimer's Disease Neuroimaging Initiative. AD was determined by clinical and neuropsychological examination with no comorbidities present. Participants with MCI had reported memory complaints, and controls were cognitively normal. A neuroanatomical normative model indexed cortical thickness distributions using a separate healthy reference data set (n = 33,072), which used hierarchical Bayesian regression to predict cortical thickness per region using age and sex, while adjusting for site noise. Z-scores per region were calculated, resulting in a Z-score brain map per participant. Regions with Z-scores <-1.96 were classified as outliers.

RESULTS

Patients with AD (n = 206) had a median of 12 outlier regions (out of a possible 148), with the highest proportion of outliers (47%) in the parahippocampal gyrus. For 62 regions, over 90% of these patients had cortical thicknesses within the normal range. Patients with AD had more outlier regions than people with MCI (n = 662) or controls (n = 159) (F(2, 1,022) = 95.39, p = 2.0 × 10-16). They were also more dissimilar to each other than people with MCI or controls (F(2, 1,024) = 209.42, p = 2.2 × 10-16). A greater number of outlier regions were associated with worse cognitive function, CSF protein concentrations, and an increased risk of converting from MCI to AD within 3 years (hazard ratio 1.028, 95% CI 1.016-1.039, p = 1.8 × 10-16).

DISCUSSION

Individualized normative maps of cortical thickness highlight the heterogeneous effect of AD on the brain. Regional outlier estimates have the potential to be a marker of disease and could be used to track an individual's disease progression or treatment response in clinical trials.

Current Approaches in Computational Psychiatry for the Data-Driven Identification of Brain-Based Subtypes

The ability of our current psychiatric nosology to accurately delineate clinical populations and inform effective treatment plans has reached a critical point with only moderately successful interventions and high relapse rates. These challenges continue to motivate the search for approaches to better stratify clinical populations into more homogeneous delineations, to better inform diagnosis and disease evaluation, and prescribe and develop more precise treatment plans. The promise of brain-based subtyping based on neuroimaging data is that finding subgroups of individuals with a common biological signature will facilitate the development of biologically grounded, targeted treatments. This review provides a snapshot of the current state of the field in empirical brain-based subtyping studies in child, adolescent, and adult psychiatric populations published between 2019 and March 2022. We found that there is vast methodological exploration and a surprising number of new methods being created for the specific purpose of brain-based subtyping. However, this methodological exploration and advancement is not being met with rigorous validation approaches that assess both reproducibility and clinical utility of the discovered brain-based subtypes. We also found evidence for a collaboration crisis, in which methodological exploration and advancements are not clearly grounded in clinical goals. We propose several steps that we believe are crucial to address these shortcomings in the field. We conclude, and agree with the authors of the reviewed studies, that the discovery of biologically grounded subtypes would be a significant advancement for treatment development in psychiatry.

AI-based dimensional neuroimaging system for characterizing heterogeneity in brain structure and function in major depressive disorder: COORDINATE-MDD consortium design and rationale

BACKGROUND

Efforts to develop neuroimaging-based biomarkers in major depressive disorder (MDD), at the individual level, have been limited to date. As diagnostic criteria are currently symptom-based, MDD is conceptualized as a disorder rather than a disease with a known etiology; further, neural measures are often confounded by medication status and heterogeneous symptom states.

METHODS

We describe a consortium to quantify neuroanatomical and neurofunctional heterogeneity via the dimensions of novel multivariate coordinate system (COORDINATE-MDD). Utilizing imaging harmonization and machine learning methods in a large cohort of medication-free, deeply phenotyped MDD participants, patterns of brain alteration are defined in replicable and neurobiologically-based dimensions and offer the potential to predict treatment response at the individual level. International datasets are being shared from multi-ethnic community populations, first episode and recurrent MDD, which are medication-free, in a current depressive episode with prospective longitudinal treatment outcomes and in remission. Neuroimaging data consist of de-identified, individual, structural MRI and resting-state functional MRI with additional positron emission tomography (PET) data at specific sites. State-of-the-art analytic methods include automated image processing for extraction of anatomical and functional imaging variables, statistical harmonization of imaging variables to account for site and scanner variations, and semi-supervised machine learning methods that identify dominant patterns associated with MDD from neural structure and function in healthy participants.

RESULTS

We are applying an iterative process by defining the neural dimensions that characterise deeply phenotyped samples and then testing the dimensions in novel samples to assess specificity and reliability. Crucially, we aim to use machine learning methods to identify novel predictors of treatment response based on prospective longitudinal treatment outcome data, and we can externally validate the dimensions in fully independent sites.

CONCLUSION

We describe the consortium, imaging protocols and analytics using preliminary results. Our findings thus far demonstrate how datasets across many sites can be harmonized and constructively pooled to enable execution of this large-scale project.

Transdiagnostic inflammatory subgroups among psychiatric disorders and their relevance to role functioning: a nested case-control study of the ALSPAC cohort

Individuals with psychotic disorders and depressive disorder exhibit altered concentrations of peripheral inflammatory markers. It has been suggested that clinical trials of anti-inflammatory therapies for psychiatric disorders should stratify patients by their inflammatory profile. Hence, we investigated whether different subgroups of individuals exist across psychiatric disorders, based on their inflammatory biomarker signatures. We measured the plasma concentrations of 17 inflammatory markers and receptors in 380 participants with psychotic disorder, depressive disorder or generalised anxiety disorder and 399 controls without psychiatric symptoms from the ALSPAC cohort at age 24. We employed a semi-supervised clustering algorithm, which discriminates multiple clusters of psychiatric disorder cases from controls. The best fit was for a two-cluster model of participants with psychiatric disorders (Adjusted Rand Index (ARI) = 0.52 ± 0.01) based on the inflammatory markers. Permutation analysis indicated the stability of the clustering solution performed better than chance (ARI = 0.43 ± 0.11; p < 0.001), and the clusters explained the inflammatory marker data better than a Gaussian distribution (p = 0.021). Cluster 2 exhibited marked increases in sTNFR1/2, suPAR, sCD93 and sIL-2RA, compared to cluster 1. Participants in the cluster exhibiting higher inflammation were less likely to be in employment, education or training, indicating poorer role functioning. This study found evidence for a novel pattern of inflammatory markers specific to psychiatric disorders and strongly associated with a transdiagnostic measure of illness severity. sTNFR1/2, suPAR, sCD93 and sIL-2RA could be used to stratify clinical trials of anti-inflammatory therapies for psychiatric disorders.

Pseudo-labeling generative adversarial networks for medical image classification

Semi-supervised learning has become a popular technology in recent years. In this paper, we propose a novel semi-supervised medical image classification algorithm, called Pseudo-Labeling Generative Adversarial Networks (PLGAN), which only uses a small number of real images with few labels to generate fake images or mask images to enlarge the sample size of the labeled training set. First, we combine MixMatch to generate pseudo labels for the fake and unlabeled images to do the classification. Second, contrastive learning and self-attention mechanisms are introduced into PLGAN to exclude the influence of unimportant details. Third, the problem of mode collapse in contrastive learning is well addressed by cyclic consistency loss. Finally, we design global and local classifiers to complement each other with the key information needed for classification. The experimental results on four medical image datasets show that PLGAN can obtain relatively high learning performance by using few labeled and unlabeled data. For example, the classification accuracy of PLGAN is 11% higher than that of MixMatch with 100 labeled images and 1000 unlabeled images on the OCT dataset. In addition, we also conduct other experiments to verify the effectiveness of our algorithm.

Characterizing Heterogeneity in Neuroimaging, Cognition, Clinical Symptoms, and Genetics Among Patients With Late-Life Depression

Importance

Late-life depression (LLD) is characterized by considerable heterogeneity in clinical manifestation. Unraveling such heterogeneity might aid in elucidating etiological mechanisms and support precision and individualized medicine.

Objective

To cross-sectionally and longitudinally delineate disease-related heterogeneity in LLD associated with neuroanatomy, cognitive functioning, clinical symptoms, and genetic profiles.

Design, Setting, and Participants

The Imaging-Based Coordinate System for Aging and Neurodegenerative Diseases (iSTAGING) study is an international multicenter consortium investigating brain aging in pooled and harmonized data from 13 studies with more than 35 000 participants, including a subset of individuals with major depressive disorder. Multimodal data from a multicenter sample (N = 996), including neuroimaging, neurocognitive assessments, and genetics, were analyzed in this study. A semisupervised clustering method (heterogeneity through discriminative analysis) was applied to regional gray matter (GM) brain volumes to derive dimensional representations. Data were collected from July 2017 to July 2020 and analyzed from July 2020 to December 2021.

Main Outcomes and Measures

Two dimensions were identified to delineate LLD-associated heterogeneity in voxelwise GM maps, white matter (WM) fractional anisotropy, neurocognitive functioning, clinical phenotype, and genetics.

Results

A total of 501 participants with LLD (mean [SD] age, 67.39 [5.56] years; 332 women) and 495 healthy control individuals (mean [SD] age, 66.53 [5.16] years; 333 women) were included. Patients in dimension 1 demonstrated relatively preserved brain anatomy without WM disruptions relative to healthy control individuals. In contrast, patients in dimension 2 showed widespread brain atrophy and WM integrity disruptions, along with cognitive impairment and higher depression severity. Moreover, 1 de novo independent genetic variant (rs13120336; chromosome: 4, 186387714; minor allele, G) was significantly associated with dimension 1 (odds ratio, 2.35; SE, 0.15; P = 3.14 ×108) but not with dimension 2. The 2 dimensions demonstrated significant single-nucleotide variant-based heritability of 18% to 27% within the general population (N = 12 518 in UK Biobank). In a subset of individuals having longitudinal measurements, those in dimension 2 experienced a more rapid longitudinal change in GM and brain age (Cohen f2 = 0.03; P = .02) and were more likely to progress to Alzheimer disease (Cohen f2 = 0.03; P = .03) compared with those in dimension 1 (N = 1431 participants and 7224 scans from the Alzheimer's Disease Neuroimaging Initiative [ADNI], Baltimore Longitudinal Study of Aging [BLSA], and Biomarkers for Older Controls at Risk for Dementia [BIOCARD] data sets).

Conclusions and Relevance

This study characterized heterogeneity in LLD into 2 dimensions with distinct neuroanatomical, cognitive, clinical, and genetic profiles. This dimensional approach provides a potential mechanism for investigating the heterogeneity of LLD and the relevance of the latent dimensions to possible disease mechanisms, clinical outcomes, and responses to interventions.