Genomic loci influence patterns of structural covariance in the human brain.

Normal and pathologic neurobiological processes influence brain morphology in coordinated ways that give rise to patterns of structural covariance (PSC) across brain regions and individuals during brain aging and diseases. The genetic underpinnings of these patterns remain largely unknown. We apply a stochastic multivariate factorization method to a diverse population of 50,699 individuals (12 studies and 130 sites) and derive data-driven, multi-scale PSCs of regional brain size. PSCs were significantly correlated with 915 genomic loci in the discovery set, 617 of which are newly identified, and 72% were independently replicated. Key pathways influencing PSCs involve reelin signaling, apoptosis, neurogenesis, and appendage development, while pathways of breast cancer indicate potential interplays between brain metastasis and PSCs associated with neurodegeneration and dementia. Using support vector machines, multi-scale PSCs effectively derive imaging signatures of several brain diseases. Our results elucidate genetic and biological underpinnings that influence structural covariance patterns in the human brain.

AmyloidPETNet: Classification of Amyloid Positivity in Brain PET Imaging Using End-to-End Deep Learning.

Background Visual assessment of amyloid PET scans relies on the availability of radiologist expertise, whereas quantification of amyloid burden typically involves MRI for processing and analysis, which can be computationally expensive. Purpose To develop a deep learning model to classify minimally processed brain PET scans as amyloid positive or negative, evaluate its performance on independent data sets and different tracers, and compare it with human visual reads. Materials and Methods This retrospective study used 8476 PET scans (6722 patients) obtained from late 2004 to early 2023 that were analyzed across five different data sets. A deep learning model, AmyloidPETNet, was trained on 1538 scans from 766 patients, validated on 205 scans from 95 patients, and internally tested on 184 scans from 95 patients in the Alzheimer's Disease Neuroimaging Initiative (ADNI) fluorine 18 (18F) florbetapir (FBP) data set. It was tested on ADNI scans using different tracers and scans from independent data sets. Scan amyloid positivity was based on mean cortical standardized uptake value ratio cutoffs. To compare with model performance, each scan from both the Centiloid Project and a subset of the Anti-Amyloid Treatment in Asymptomatic Alzheimer's Disease (A4) study were visually interpreted with a confidence level (low, intermediate, high) of amyloid positivity/negativity. The area under the receiver operating characteristic curve (AUC) and other performance metrics were calculated, and Cohen κ was used to measure physician-model agreement. Results The model achieved an AUC of 0.97 (95% CI: 0.95, 0.99) on test ADNI 18F-FBP scans, which generalized well to 18F-FBP scans from the Open Access Series of Imaging Studies (AUC, 0.95; 95% CI: 0.93, 0.97) and the A4 study (AUC, 0.98; 95% CI: 0.98, 0.98). Model performance was high when applied to data sets with different tracers (AUC ≥ 0.97). Other performance metrics provided converging evidence. Physician-model agreement ranged from fair (Cohen κ = 0.39; 95% CI: 0.16, 0.60) on a sample of mostly equivocal cases from the A4 study to almost perfect (Cohen κ = 0.93; 95% CI: 0.86, 1.0) on the Centiloid Project. Conclusion The developed model was capable of automatically and accurately classifying brain PET scans as amyloid positive or negative without relying on experienced readers or requiring structural MRI. Clinical trial registration no. NCT00106899 © RSNA, 2024 Supplemental material is available for this article. See also the editorial by Bryan and Forghani in this issue.

Psychosis brain subtypes validated in first-episode cohorts and related to illness remission: results from the PHENOM consortium.

Using machine learning, we recently decomposed the neuroanatomical heterogeneity of established schizophrenia to discover two volumetric subgroups-a 'lower brain volume' subgroup (SG1) and an 'higher striatal volume' subgroup (SG2) with otherwise normal brain structure. In this study, we investigated whether the MRI signatures of these subgroups were also already present at the time of the first-episode of psychosis (FEP) and whether they were related to clinical presentation and clinical remission over 1-, 3-, and 5-years. We included 572 FEP and 424 healthy controls (HC) from 4 sites (Sao Paulo, Santander, London, Melbourne) of the PHENOM consortium. Our prior MRI subgrouping models (671 participants; USA, Germany, and China) were applied to both FEP and HC. Participants were assigned into 1 of 4 categories: subgroup 1 (SG1), subgroup 2 (SG2), no subgroup membership ('None'), and mixed SG1 + SG2 subgroups ('Mixed'). Voxel-wise analyses characterized SG1 and SG2 subgroups. Supervised machine learning analyses characterized baseline and remission signatures related to SG1 and SG2 membership. The two dominant patterns of 'lower brain volume' in SG1 and 'higher striatal volume' (with otherwise normal neuromorphology) in SG2 were identified already at the first episode of psychosis. SG1 had a significantly higher proportion of FEP (32%) vs. HC (19%) than SG2 (FEP, 21%; HC, 23%). Clinical multivariate signatures separated the SG1 and SG2 subgroups (balanced accuracy = 64%; p < 0.0001), with SG2 showing higher education but also greater positive psychosis symptoms at first presentation, and an association with symptom remission at 1-year, 5-year, and when timepoints were combined. Neuromorphological subtypes of schizophrenia are already evident at illness onset, separated by distinct clinical presentations, and differentially associated with subsequent remission. These results suggest that the subgroups may be underlying risk phenotypes that could be targeted in future treatment trials and are critical to consider when interpreting neuroimaging literature.

Data-driven decomposition and staging of flortaucipir uptake in Alzheimer's disease.

INTRODUCTION: Previous approaches pursuing in vivo staging of tau pathology in Alzheimer's disease (AD) have typically relied on neuropathologically defined criteria. In using predefined systems, these studies may miss spatial deposition patterns which are informative of disease progression. METHODS: We selected discovery (n = 418) and replication (n = 132) cohorts with flortaucipir imaging. Non-negative matrix factorization (NMF) was applied to learn tau covariance patterns and develop a tau staging system. Flortaucipir components were also validated by comparison with amyloid burden, gray matter loss, and the expression of AD-related genes. RESULTS: We found eight flortaucipir covariance patterns which were reproducible and overlapped with relevant gene expression maps. Tau stages were associated with AD severity as indexed by dementia status and neuropsychological performance. Comparisons of flortaucipir uptake with amyloid and atrophy also supported our model of tau progression. DISCUSSION: Data-driven decomposition of flortaucipir uptake provides a novel framework for tau staging which complements existing systems. HIGHLIGHTS: NMF reveals patterns of tau deposition in AD. Data-driven staging of flortaucipir tracks AD severity. Learned flortaucipir patterns overlap with AD-related gene expression.

Benchmarking the generalizability of brain age models: Challenges posed by scanner variance and prediction bias.

Machine learning has been increasingly applied to neuroimaging data to predict age, deriving a personalized biomarker with potential clinical applications. The scientific and clinical value of these models depends on their applicability to independently acquired scans from diverse sources. Accordingly, we evaluated the generalizability of two brain age models that were trained across the lifespan by applying them to three distinct early-life samples with participants aged 8-22 years. These models were chosen based on the size and diversity of their training data, but they also differed greatly in their processing methods and predictive algorithms. Specifically, one brain age model was built by applying gradient tree boosting (GTB) to extracted features of cortical thickness, surface area, and brain volume. The other model applied a 2D convolutional neural network (DBN) to minimally preprocessed slices of T1-weighted scans. Additional model variants were created to understand how generalizability changed when each model was trained with data that became more similar to the test samples in terms of age and acquisition protocols. Our results illustrated numerous trade-offs. The GTB predictions were relatively more accurate overall and yielded more reliable predictions when applied to lower quality scans. In contrast, the DBN displayed the most utility in detecting associations between brain age gaps and cognitive functioning. Broadly speaking, the largest limitations affecting generalizability were acquisition protocol differences and biased brain age estimates. If such confounds could eventually be removed without post-hoc corrections, brain age predictions may have greater utility as personalized biomarkers of healthy aging.

Brain deficit patterns of metabolic illnesses overlap with those for major depressive disorder: A new metric of brain metabolic disease.

Metabolic illnesses (MET) are detrimental to brain integrity and are common comorbidities in patients with mental illnesses, including major depressive disorder (MDD). We quantified effects of MET on standard regional brain morphometric measures from 3D brain MRI as well as diffusion MRI in a large sample of UK BioBank participants. The pattern of regional effect sizes of MET in non-psychiatric UKBB subjects was significantly correlated with the spatial profile of regional effects reported by the largest meta-analyses in MDD but not in bipolar disorder, schizophrenia or Alzheimer's disease. We used a regional vulnerability index (RVI) for MET (RVI-MET) to measure individual's brain similarity to the expected patterns in MET in the UK Biobank sample. Subjects with MET showed a higher effect size for RVI-MET than for any of the individual brain measures. We replicated elevation of RVI-MET in a sample of MDD participants with MET versus non-MET. RVI-MET scores were significantly correlated with the volume of white matter hyperintensities, a neurological consequence of MET and age, in both groups. Higher RVI-MET in both samples was associated with obesity, tobacco smoking and frequent alcohol use but was unrelated to antidepressant use. In summary, MET effects on the brain were regionally specific and individual similarity to the pattern was more strongly associated with MET than any regional brain structural metric. Effects of MET overlapped with the reported brain differences in MDD, likely due to higher incidence of MET, smoking and alcohol use in subjects with MDD.

Novel Gyrification Networks Reveal Links with Psychiatric Risk Factors in Early Illness.

Adult gyrification provides a window into coordinated early neurodevelopment when disruptions predispose individuals to psychiatric illness. We hypothesized that the echoes of such disruptions should be observed within structural gyrification networks in early psychiatric illness that would demonstrate associations with developmentally relevant variables rather than specific psychiatric symptoms. We employed a new data-driven method (Orthogonal Projective Non-Negative Matrix Factorization) to delineate novel gyrification-based networks of structural covariance in 308 healthy controls. Gyrification within the networks was then compared to 713 patients with recent onset psychosis or depression, and at clinical high-risk. Associations with diagnosis, symptoms, cognition, and functioning were investigated using linear models. Results demonstrated 18 novel gyrification networks in controls as verified by internal and external validation. Gyrification was reduced in patients in temporal-insular, lateral occipital, and lateral fronto-parietal networks (pFDR < 0.01) and was not moderated by illness group. Higher gyrification was associated with better cognitive performance and lifetime role functioning, but not with symptoms. The findings demonstrated that gyrification can be parsed into novel brain networks that highlight generalized illness effects linked to developmental vulnerability. When combined, our study widens the window into the etiology of psychiatric risk and its expression in adulthood.

White matter hyperintensity longitudinal morphometric analysis in association with Alzheimer disease.

INTRODUCTION: Vascular damage in Alzheimer's disease (AD) has shown conflicting findings particularly when analyzing longitudinal data. We introduce white matter hyperintensity (WMH) longitudinal morphometric analysis (WLMA) that quantifies WMH expansion as the distance from lesion voxels to a region of interest boundary. METHODS: WMH segmentation maps were derived from 270 longitudinal fluid-attenuated inversion recovery (FLAIR) ADNI images. WLMA was performed on five data-driven WMH patterns with distinct spatial distributions. Amyloid accumulation was evaluated with WMH expansion across the five WMH patterns. RESULTS: The preclinical group had significantly greater expansion in the posterior ventricular WM compared to controls. Amyloid significantly associated with frontal WMH expansion primarily within AD individuals. WLMA outperformed WMH volume changes for classifying AD from controls primarily in periventricular and posterior WMH. DISCUSSION: These data support the concept that localized WMH expansion continues to proliferate with amyloid accumulation throughout the entirety of the disease in distinct spatial locations.

Synthesizing pseudo-T2w images to recapture missing data in neonatal neuroimaging with applications in rs-fMRI.

T1- and T2-weighted (T1w and T2w) images are essential for tissue classification and anatomical localization in Magnetic Resonance Imaging (MRI) analyses. However, these anatomical data can be challenging to acquire in non-sedated neonatal cohorts, which are prone to high amplitude movement and display lower tissue contrast than adults. As a result, one of these modalities may be missing or of such poor quality that they cannot be used for accurate image processing, resulting in subject loss. While recent literature attempts to overcome these issues in adult populations using synthetic imaging approaches, evaluation of the efficacy of these methods in pediatric populations and the impact of these techniques in conventional MR analyses has not been performed. In this work, we present two novel methods to generate pseudo-T2w images: the first is based in deep learning and expands upon previous models to 3D imaging without the requirement of paired data, the second is based in nonlinear multi-atlas registration providing a computationally lightweight alternative. We demonstrate the anatomical accuracy of pseudo-T2w images and their efficacy in existing MR processing pipelines in two independent neonatal cohorts. Critically, we show that implementing these pseudo-T2w methods in resting-state functional MRI analyses produces virtually identical functional connectivity results when compared to those resulting from T2w images, confirming their utility in infant MRI studies for salvaging otherwise lost subject data.

Age-dependent white matter disruptions after military traumatic brain injury: Multivariate analysis results from ENIGMA brain injury.

Mild Traumatic brain injury (mTBI) is a signature wound in military personnel, and repetitive mTBI has been linked to age-related neurogenerative disorders that affect white matter (WM) in the brain. However, findings of injury to specific WM tracts have been variable and inconsistent. This may be due to the heterogeneity of mechanisms, etiology, and comorbid disorders related to mTBI. Non-negative matrix factorization (NMF) is a data-driven approach that detects covarying patterns (components) within high-dimensional data. We applied NMF to diffusion imaging data from military Veterans with and without a self-reported TBI history. NMF identified 12 independent components derived from fractional anisotropy (FA) in a large dataset (n = 1,475) gathered through the ENIGMA (Enhancing Neuroimaging Genetics through Meta-Analysis) Military Brain Injury working group. Regressions were used to examine TBI- and mTBI-related associations in NMF-derived components while adjusting for age, sex, post-traumatic stress disorder, depression, and data acquisition site/scanner. We found significantly stronger age-dependent effects of lower FA in Veterans with TBI than Veterans without in four components (q < 0.05), which are spatially unconstrained by traditionally defined WM tracts. One component, occupying the most peripheral location, exhibited significantly stronger age-dependent differences in Veterans with mTBI. We found NMF to be powerful and effective in detecting covarying patterns of FA associated with mTBI by applying standard parametric regression modeling. Our results highlight patterns of WM alteration that are differentially affected by TBI and mTBI in younger compared to older military Veterans.

Differential cortical microstructural maturation in the preterm human brain with diffusion kurtosis and tensor imaging.

During the third trimester, the human brain undergoes rapid cellular and molecular processes that reshape the structural architecture of the cerebral cortex. Knowledge of cortical differentiation obtained predominantly from histological studies is limited in localized and small cortical regions. How cortical microstructure is differentiated across cortical regions in this critical period is unknown. In this study, the cortical microstructural architecture across the entire cortex was delineated with non-Gaussian diffusion kurtosis imaging as well as conventional diffusion tensor imaging of 89 preterm neonates aged 31-42 postmenstrual weeks. The temporal changes of cortical mean kurtosis (MK) or fractional anisotropy (FA) were heterogeneous across the cortical regions. Cortical MK decreases were observed throughout the studied age period, while cortical FA decrease reached its plateau around 37 weeks. More rapid decreases in MK were found in the primary visual region, while faster FA declines were observed in the prefrontal cortex. We found that distinctive cortical microstructural changes were coupled with microstructural maturation of associated white matter tracts. Both cortical MK and FA measurements predicted the postmenstrual age of preterm infants accurately. This study revealed a differential 4D spatiotemporal cytoarchitectural signature inferred by non-Gaussian diffusion barriers inside the cortical plate during the third trimester. The cytoarchitectural processes, including dendritic arborization and neuronal density decreases, were inferred by regional cortical FA and MK measurements. The presented findings suggest that cortical MK and FA measurements could be used as effective imaging markers for cortical microstructural changes in typical and potentially atypical brain development.

Two distinct neuroanatomical subtypes of schizophrenia revealed using machine learning.

Neurobiological heterogeneity in schizophrenia is poorly understood and confounds current analyses. We investigated neuroanatomical subtypes in a multi-institutional multi-ethnic cohort, using novel semi-supervised machine learning methods designed to discover patterns associated with disease rather than normal anatomical variation. Structural MRI and clinical measures in established schizophrenia (n = 307) and healthy controls (n = 364) were analysed across three sites of PHENOM (Psychosis Heterogeneity Evaluated via Dimensional Neuroimaging) consortium. Regional volumetric measures of grey matter, white matter, and CSF were used to identify distinct and reproducible neuroanatomical subtypes of schizophrenia. Two distinct neuroanatomical subtypes were found. Subtype 1 showed widespread lower grey matter volumes, most prominent in thalamus, nucleus accumbens, medial temporal, medial prefrontal/frontal and insular cortices. Subtype 2 showed increased volume in the basal ganglia and internal capsule, and otherwise normal brain volumes. Grey matter volume correlated negatively with illness duration in Subtype 1 (r = -0.201, P = 0.016) but not in Subtype 2 (r = -0.045, P = 0.652), potentially indicating different underlying neuropathological processes. The subtypes did not differ in age (t = -1.603, df = 305, P = 0.109), sex (chi-square = 0.013, df = 1, P = 0.910), illness duration (t = -0.167, df = 277, P = 0.868), antipsychotic dose (t = -0.439, df = 210, P = 0.521), age of illness onset (t = -1.355, df = 277, P = 0.177), positive symptoms (t = 0.249, df = 289, P = 0.803), negative symptoms (t = 0.151, df = 289, P = 0.879), or antipsychotic type (chi-square = 6.670, df = 3, P = 0.083). Subtype 1 had lower educational attainment than Subtype 2 (chi-square = 6.389, df = 2, P = 0.041). In conclusion, we discovered two distinct and highly reproducible neuroanatomical subtypes. Subtype 1 displayed widespread volume reduction correlating with illness duration, and worse premorbid functioning. Subtype 2 had normal and stable anatomy, except for larger basal ganglia and internal capsule, not explained by antipsychotic dose. These subtypes challenge the notion that brain volume loss is a general feature of schizophrenia and suggest differential aetiologies. They can facilitate strategies for clinical trial enrichment and stratification, and precision diagnostics.

The Brain Chart of Aging: Machine-learning analytics reveals links between brain aging, white matter disease, amyloid burden, and cognition in the iSTAGING consortium of 10,216 harmonized MR scans.

INTRODUCTION: Relationships between brain atrophy patterns of typical aging and Alzheimer's disease (AD), white matter disease, cognition, and AD neuropathology were investigated via machine learning in a large harmonized magnetic resonance imaging database (11 studies; 10,216 subjects). METHODS: Three brain signatures were calculated: Brain-age, AD-like neurodegeneration, and white matter hyperintensities (WMHs). Brain Charts measured and displayed the relationships of these signatures to cognition and molecular biomarkers of AD. RESULTS: WMHs were associated with advanced brain aging, AD-like atrophy, poorer cognition, and AD neuropathology in mild cognitive impairment (MCI)/AD and cognitively normal (CN) subjects. High WMH volume was associated with brain aging and cognitive decline occurring in an ≈10-year period in CN subjects. WMHs were associated with doubling the likelihood of amyloid beta (Aß) positivity after age 65. Brain aging, AD-like atrophy, and WMHs were better predictors of cognition than chronological age in MCI/AD. DISCUSSION: A Brain Chart quantifying brain-aging trajectories was established, enabling the systematic evaluation of individuals' brain-aging patterns relative to this large consortium.

Gestational Age is Dimensionally Associated with Structural Brain Network Abnormalities Across Development.

Prematurity is associated with diverse developmental abnormalities, yet few studies relate cognitive and neurostructural deficits to a dimensional measure of prematurity. Leveraging a large sample of children, adolescents, and young adults (age 8-22 years) studied as part of the Philadelphia Neurodevelopmental Cohort, we examined how variation in gestational age impacted cognition and brain structure later in development. Participants included 72 preterm youth born before 37 weeks' gestation and 206 youth who were born at term (37 weeks or later). Using a previously-validated factor analysis, cognitive performance was assessed in three domains: (1) executive function and complex reasoning, (2) social cognition, and (3) episodic memory. All participants completed T1-weighted neuroimaging at 3 T to measure brain volume. Structural covariance networks were delineated using non-negative matrix factorization, an advanced multivariate analysis technique. Lower gestational age was associated with both deficits in executive function and reduced volume within 11 of 26 structural covariance networks, which included orbitofrontal, temporal, and parietal cortices as well as subcortical regions including the hippocampus. Notably, the relationship between lower gestational age and executive dysfunction was accounted for in part by structural network deficits. Together, these findings emphasize the durable impact of prematurity on cognition and brain structure, which persists across development.

Patterns of coordinated cortical remodeling during adolescence and their associations with functional specialization and evolutionary expansion.

During adolescence, the human cortex undergoes substantial remodeling to support a rapid expansion of behavioral repertoire. Accurately quantifying these changes is a prerequisite for understanding normal brain development, as well as the neuropsychiatric disorders that emerge in this vulnerable period. Past accounts have demonstrated substantial regional heterogeneity in patterns of brain development, but frequently have been limited by small samples and analytics that do not evaluate complex multivariate imaging patterns. Capitalizing on recent advances in multivariate analysis methods, we used nonnegative matrix factorization (NMF) to uncover coordinated patterns of cortical development in a sample of 934 youths ages 8-20, who completed structural neuroimaging as part of the Philadelphia Neurodevelopmental Cohort. Patterns of structural covariance (PSCs) derived by NMF were highly reproducible over a range of resolutions, and differed markedly from common gyral-based structural atlases. Moreover, PSCs were largely symmetric and showed correspondence to specific large-scale functional networks. The level of correspondence was ordered according to their functional role and position in the evolutionary hierarchy, being high in lower-order visual and somatomotor networks and diminishing in higher-order association cortex. Furthermore, PSCs showed divergent developmental associations, with PSCs in higher-order association cortex networks showing greater changes with age than primary somatomotor and visual networks. Critically, such developmental changes within PSCs were significantly associated with the degree of evolutionary cortical expansion. Together, our findings delineate a set of structural brain networks that undergo coordinated cortical thinning during adolescence, which is in part governed by evolutionary novelty and functional specialization.

Diminished Cortical Thickness Is Associated with Impulsive Choice in Adolescence.

Adolescence is characterized by both maturation of brain structure and increased risk of negative outcomes from behaviors associated with impulsive decision-making. One important index of impulsive choice is delay discounting (DD), which measures the tendency to prefer smaller rewards available soon over larger rewards delivered after a delay. However, it remains largely unknown how individual differences in structural brain development may be associated with impulsive choice during adolescence. Leveraging a unique large sample of 427 human youths (208 males and 219 females) imaged as part of the Philadelphia Neurodevelopmental Cohort, we examined associations between delay discounting and cortical thickness within structural covariance networks. These structural networks were derived using non-negative matrix factorization, an advanced multivariate technique for dimensionality reduction, and analyzed using generalized additive models with penalized splines to capture both linear and nonlinear developmental effects. We found that impulsive choice, as measured by greater discounting, was most strongly associated with diminished cortical thickness in structural brain networks that encompassed the ventromedial prefrontal cortex, orbitofrontal cortex, temporal pole, and temporoparietal junction. Furthermore, structural brain networks predicted DD above and beyond cognitive performance. Together, these results suggest that reduced cortical thickness in regions known to be involved in value-based decision-making is a marker of impulsive choice during the critical period of adolescence.SIGNIFICANCE STATEMENT Risky behaviors during adolescence, such as initiation of substance use or reckless driving, are a major source of morbidity and mortality. In this study, we present evidence from a large sample of youths that diminished cortical thickness in specific structural brain networks is associated with impulsive choice. Notably, the strongest association between impulsive choice and brain structure was seen in regions implicated in value-based decision-making; namely, the ventromedial prefrontal and orbitofrontal cortices. Moving forward, such neuroanatomical markers of impulsivity may aid in the development of personalized interventions targeted to reduce risk of negative outcomes resulting from impulsivity during adolescence.

MIDAS: Regionally linear multivariate discriminative statistical mapping.

Statistical parametric maps formed via voxel-wise mass-univariate tests, such as the general linear model, are commonly used to test hypotheses about regionally specific effects in neuroimaging cross-sectional studies where each subject is represented by a single image. Despite being informative, these techniques remain limited as they ignore multivariate relationships in the data. Most importantly, the commonly employed local Gaussian smoothing, which is important for accounting for registration errors and making the data follow Gaussian distributions, is usually chosen in an ad hoc fashion. Thus, it is often suboptimal for the task of detecting group differences and correlations with non-imaging variables. Information mapping techniques, such as searchlight, which use pattern classifiers to exploit multivariate information and obtain more powerful statistical maps, have become increasingly popular in recent years. However, existing methods may lead to important interpretation errors in practice (i.e., misidentifying a cluster as informative, or failing to detect truly informative voxels), while often being computationally expensive. To address these issues, we introduce a novel efficient multivariate statistical framework for cross-sectional studies, termed MIDAS, seeking highly sensitive and specific voxel-wise brain maps, while leveraging the power of regional discriminant analysis. In MIDAS, locally linear discriminative learning is applied to estimate the pattern that best discriminates between two groups, or predicts a variable of interest. This pattern is equivalent to local filtering by an optimal kernel whose coefficients are the weights of the linear discriminant. By composing information from all neighborhoods that contain a given voxel, MIDAS produces a statistic that collectively reflects the contribution of the voxel to the regional classifiers as well as the discriminative power of the classifiers. Critically, MIDAS efficiently assesses the statistical significance of the derived statistic by analytically approximating its null distribution without the need for computationally expensive permutation tests. The proposed framework was extensively validated using simulated atrophy in structural magnetic resonance imaging (MRI) and further tested using data from a task-based functional MRI study as well as a structural MRI study of cognitive performance. The performance of the proposed framework was evaluated against standard voxel-wise general linear models and other information mapping methods. The experimental results showed that MIDAS achieves relatively higher sensitivity and specificity in detecting group differences. Together, our results demonstrate the potential of the proposed approach to efficiently map effects of interest in both structural and functional data.

Evaluation of non-negative matrix factorization of grey matter in age prediction.

The relationship between grey matter volume (GMV) patterns and age can be captured by multivariate pattern analysis, allowing prediction of individuals' age based on structural imaging. Raw data, voxel-wise GMV and non-sparse factorization (with Principal Component Analysis, PCA) show good performance but do not promote relatively localized brain components for post-hoc examinations. Here we evaluated a non-negative matrix factorization (NNMF) approach to provide a reduced, but also interpretable representation of GMV data in age prediction frameworks in healthy and clinical populations. This examination was performed using three datasets: a multi-site cohort of life-span healthy adults, a single site cohort of older adults and clinical samples from the ADNI dataset with healthy subjects, participants with Mild Cognitive Impairment and patients with Alzheimer's disease (AD) subsamples. T1-weighted images were preprocessed with VBM8 standard settings to compute GMV values after normalization, segmentation and modulation for non-linear transformations only. Non-negative matrix factorization was computed on the GM voxel-wise values for a range of granularities (50-690 components) and LASSO (Least Absolute Shrinkage and Selection Operator) regression were used for age prediction. First, we compared the performance of our data compression procedure (i.e., NNMF) to various other approaches (i.e., uncompressed VBM data, PCA-based factorization and parcellation-based compression). We then investigated the impact of the granularity on the accuracy of age prediction, as well as the transferability of the factorization and model generalization across datasets. We finally validated our framework by examining age prediction in ADNI samples. Our results showed that our framework favorably compares with other approaches. They also demonstrated that the NNMF based factorization derived from one dataset could be efficiently applied to compress VBM data of another dataset and that granularities between 300 and 500 components give an optimal representation for age prediction. In addition to the good performance in healthy subjects our framework provided relatively localized brain regions as the features contributing to the prediction, thereby offering further insights into structural changes due to brain aging. Finally, our validation in clinical populations showed that our framework is sensitive to deviance from normal structural variations in pathological aging.

Heterogeneity of neuroanatomical patterns in prodromal Alzheimer's disease: links to cognition, progression and biomarkers.

See Coulthard and Knight (doi:10.1093/aww335) for a scientific commentary on this article.Individuals with mild cognitive impairment and Alzheimer's disease clinical diagnoses can display significant phenotypic heterogeneity. This variability likely reflects underlying genetic, environmental and neuropathological differences. Characterizing this heterogeneity is important for precision diagnostics, personalized predictions, and recruitment of relatively homogeneous sets of patients into clinical trials. In this study, we apply state-of-the-art semi-supervised machine learning methods to the Alzheimer's disease Neuroimaging cohort (ADNI) to elucidate the heterogeneity of neuroanatomical differences between subjects with mild cognitive impairment (n = 530) and Alzheimer's disease (n = 314) and cognitively normal individuals (n = 399), thereby adding to an increasing literature aiming to establish neuroanatomical and neuropathological (e.g. amyloid and tau deposition) dimensions in Alzheimer's disease and its prodromal stages. These dimensional approaches aim to provide surrogate measures of heterogeneous underlying pathologic processes leading to cognitive impairment. We relate these neuroimaging patterns to cerebrospinal fluid biomarkers, white matter hyperintensities, cognitive and clinical measures, and longitudinal trajectories. We identified four such atrophy patterns: (i) individuals with largely normal neuroanatomical profiles, who also turned out to have the least abnormal cognitive and cerebrospinal fluid biomarker profiles and the slowest clinical progression during follow-up; (ii) individuals with classical Alzheimer's disease neuroanatomical, cognitive, cerebrospinal fluid biomarkers and clinical profile, who presented the fastest clinical progression; (iii) individuals with a diffuse pattern of atrophy with relatively less pronounced involvement of the medial temporal lobe, abnormal cerebrospinal fluid amyloid-ß1-42 values, and proportionally greater executive impairment; and (iv) individuals with notably focal involvement of the medial temporal lobe and a slow steady progression, likely representing in early Alzheimer's disease stages. These four atrophy patterns effectively define a 4-dimensional categorization of neuroanatomical alterations in mild cognitive impairment and Alzheimer's disease that can complement existing dimensional approaches for staging Alzheimer's disease using a variety of biomarkers, which offer the potential for enabling precision diagnostics and prognostics, as well as targeted patient recruitment of relatively homogeneous subgroups of subjects for clinical trials.

HYDRA: Revealing heterogeneity of imaging and genetic patterns through a multiple max-margin discriminative analysis framework.

Multivariate pattern analysis techniques have been increasingly used over the past decade to derive highly sensitive and specific biomarkers of diseases on an individual basis. The driving assumption behind the vast majority of the existing methodologies is that a single imaging pattern can distinguish between healthy and diseased populations, or between two subgroups of patients (e.g., progressors vs. non-progressors). This assumption effectively ignores the ample evidence for the heterogeneous nature of brain diseases. Neurodegenerative, neuropsychiatric and neurodevelopmental disorders are largely characterized by high clinical heterogeneity, which likely stems in part from underlying neuroanatomical heterogeneity of various pathologies. Detecting and characterizing heterogeneity may deepen our understanding of disease mechanisms and lead to patient-specific treatments. However, few approaches tackle disease subtype discovery in a principled machine learning framework. To address this challenge, we present a novel non-linear learning algorithm for simultaneous binary classification and subtype identification, termed HYDRA (Heterogeneity through Discriminative Analysis). Neuroanatomical subtypes are effectively captured by multiple linear hyperplanes, which form a convex polytope that separates two groups (e.g., healthy controls from pathologic samples); each face of this polytope effectively defines a disease subtype. We validated HYDRA on simulated and clinical data. In the latter case, we applied the proposed method independently to the imaging and genetic datasets of the Alzheimer's Disease Neuroimaging Initiative (ADNI 1) study. The imaging dataset consisted of T1-weighted volumetric magnetic resonance images of 123 AD patients and 177 controls. The genetic dataset consisted of single nucleotide polymorphism information of 103 AD patients and 139 controls. We identified 3 reproducible subtypes of atrophy in AD relative to controls: (1) diffuse and extensive atrophy, (2) precuneus and extensive temporal lobe atrophy, as well some prefrontal atrophy, (3) atrophy pattern very much confined to the hippocampus and the medial temporal lobe. The genetics dataset yielded two subtypes of AD characterized mainly by the presence/absence of the apolipoprotein E (APOE) ε4 genotype, but also involving differential presence of risk alleles of CD2AP, SPON1 and LOC39095 SNPs that were associated with differences in the respective patterns of brain atrophy, especially in the precuneus. The results demonstrate the potential of the proposed approach to map disease heterogeneity in neuroimaging and genetic studies.