In this study we introduce the regional flux analysis, a novel approach to deformation based morphometry based on the Helmholtz decomposition of deformations parameterized by stationary velocity fields. We use the scalar pressure map associated to the irrotational component of the deformation to discover the critical regions of volume change. These regions are used to consistently quantify the associated measure of volume change by the probabilistic integration of the flux of the longitudinal deformations across the boundaries. The presented framework unifies voxel-based and regional approaches, and robustly describes the volume changes at both group-wise and subject-specific level as a spatial process governed by consistently defined regions. Our experiments on the large cohorts of the ADNI dataset show that the regional flux analysis is a powerful and flexible instrument for the study of Alzheimer's disease in a wide range of scenarios: cross-sectional deformation based morphometry, longitudinal discovery and quantification of group-wise volume changes, and statistically powered and robust quantification of hippocampal and ventricular atrophy.

Background

Digital cognitive assessments, particularly those that can be done at home, present as low-burden biomarkers for participants and patients alike, but their effectiveness in the diagnosis of Alzheimer's disease (AD) or predicting its trajectory is still unclear. Here, we assessed what utility or added value these digital cognitive assessments provide for identifying those at high risk of cognitive decline.

Methods

We analyzed >500 Alzheimer's Disease Neuroimaging Initiative participants who underwent a brief digital cognitive assessment and amyloid beta (Aβ)/tau positron emission tomography scans, examining their ability to distinguish cognitive status and predict cognitive decline.

Results

Performance on the digital cognitive assessment was superior to both cortical Aβ and entorhinal tau in detecting mild cognitive impairment and future cognitive decline, with mnemonic discrimination deficits emerging as the most critical measure for predicting decline and future tau accumulation.

Discussion

Digital assessments are effective at identifying at-risk individuals, supporting their utility as low-burden tools for early AD detection and monitoring.

Highlights

Performance on digital cognitive assessments predicts progression to mild cognitive impairment at a higher proficiency compared to amyloid beta and tau. Deficits in mnemonic discrimination are indicative of future cognitive decline. Impaired mnemonic discrimination predicts future entorhinal and inferior temporal tau.

Resting-state functional magnetic resonance imaging (rs-fMRI) and other technologies have been offering evidence and insights showing that altered brain functional networks are associated with neurological illnesses such as Alzheimer's disease. Exploring brain networks of clinical populations compared to those of controls would be a key inquiry to reveal underlying neurological processes related to such illnesses. For such a purpose, group-level inference is a necessary first step in order to establish whether there are any genuinely disrupted brain subnetworks. Such an analysis is also challenging due to the high dimensionality of the parameters in a network model and high noise levels in neuroimaging data. We are still in the early stage of method development as highlighted by Varoquaux and Craddock (2013) that "there is currently no unique solution, but a spectrum of related methods and analytical strategies" to learn and compare brain connectivity. In practice the important issue of how to choose several critical parameters in estimating a network, such as what association measure to use and what is the sparsity of the estimated network, has not been carefully addressed, largely because the answers are unknown yet. For example, even though the choice of tuning parameters in model estimation has been extensively discussed in the literature, as to be shown here, an optimal choice of a parameter for network estimation may not be optimal in the current context of hypothesis testing. Arbitrarily choosing or mis-specifying such parameters may lead to extremely low-powered tests. Here we develop highly adaptive tests to detect group differences in brain connectivity while accounting for unknown optimal choices of some tuning parameters. The proposed tests combine statistical evidence against a null hypothesis from multiple sources across a range of plausible tuning parameter values reflecting uncertainty with the unknown truth. These highly adaptive tests are not only easy to use, but also high-powered robustly across various scenarios. The usage and advantages of these novel tests are demonstrated on an Alzheimer's disease dataset and simulated data.

Background

APOE is the largest genetic risk factor for Alzheimer's disease (AD), but there is a substantial polygenic component. Polygenic risk scores (PRS) can summarize small effects across the genome but may obscure differential risk across molecular processes and pathways that contribute to heterogeneity of disease presentation.

Objective

We examined polygenic risk impacting specific AD-associated pathways and its relationship with clinical status and biomarkers of amyloid, tau, and neurodegeneration (A/T/N).

Methods

We analyzed data from 1,411 participants from the Alzheimer's Disease Neuroimaging Initiative (ADNI). We applied pathway analysis and clustering to identify AD-associated "pathway clusters" and construct pathway-specific PRSs (excluding the APOE region). We tested associations with diagnostic status, abnormal levels of amyloid and ptau, and hippocampal volume.

Results

Thirteen pathway clusters were identified, and eight pathway-specific PRSs were significantly associated with AD diagnosis. Amyloid-positivity was associated with endocytosis and fibril formation, response misfolded protein, and regulation protein tyrosine PRSs. Ptau positivity and hippocampal volume were both related to protein localization and mitophagy PRS, and ptau-positivity was also associated with an immune signaling PRS. A global AD PRS showed stronger associations with diagnosis and all biomarkers compared to pathway PRSs.

Conclusions

Pathway PRS may contribute to understanding separable disease processes, but do not add significant power for predictive purposes. These findings demonstrate that AD-phenotypes may be preferentially associated with risk in specific pathways, and defining genetic risk along multiple dimensions may clarify etiological heterogeneity in AD. This approach to delineate pathway-specific PRS can be used to study other complex diseases.

Alzheimer's disease pathology (AD) originates in the hippocampus and subsequently spreads to temporal, parietal, and prefrontal association cortices in a relatively stereotyped progression. Current evidence attributes this orderly progression to transneuronal transmission of misfolded proteins along the projection pathways of affected neurons. A network diffusion model was recently proposed to mathematically predict disease topography resulting from transneuronal transmission on the brain's connectivity network. Here, we use this model to predict future patterns of regional atrophy and metabolism from baseline regional patterns of 418 subjects. The model accurately predicts end-of-study regional atrophy and metabolism starting from baseline data, with significantly higher correlation strength than given by the baseline statistics directly. The model's rate parameter encapsulates overall atrophy progression rate; group analysis revealed this rate to depend on diagnosis as well as baseline cerebrospinal fluid (CSF) biomarker levels. This work helps validate the model as a prognostic tool for Alzheimer's disease assessment.

This article assesses the feasibility of using shape information to detect and quantify the subcortical and ventricular structural changes in mild cognitive impairment (MCI) and Alzheimer's disease (AD) patients. We first demonstrate structural shape abnormalities in MCI and AD as compared with healthy controls (HC). Exploring the development to AD, we then divide the MCI participants into two subgroups based on longitudinal clinical information: (1) MCI patients who remained stable; (2) MCI patients who converted to AD over time. We focus on seven structures (amygdala, hippocampus, thalamus, caudate, putamen, globus pallidus, and lateral ventricles) in 754 MR scans (210 HC, 369 MCI of which 151 converted to AD over time, and 175 AD). The hippocampus and amygdala were further subsegmented based on high field 0.8 mm isotropic 7.0T scans for finer exploration. For MCI and AD, prominent ventricular expansions were detected and we found that these patients had strongest hippocampal atrophy occurring at CA1 and strongest amygdala atrophy at the basolateral complex. Mild atrophy in basal ganglia structures was also detected in MCI and AD. Stronger atrophy in the amygdala and hippocampus, and greater expansion in ventricles was observed in MCI converters, relative to those MCI who remained stable. Furthermore, we performed principal component analysis on a linear shape space of each structure. A subsequent linear discriminant analysis on the principal component values of hippocampus, amygdala, and ventricle leads to correct classification of 88% HC subjects and 86% AD subjects.

We proposed a diffeomorphometry-based statistical pipeline to study the regional shape change rates of the bilateral hippocampus, amygdala, and ventricle in mild cognitive impairment (MCI) and Alzheimer's disease (AD) compared with healthy controls (HC), using sequential magnetic resonance imaging (MRI) scans of 713 subjects (3,123 scans in total). The subgroup shape atrophy rates of the bilateral hippocampus and amygdala, as well as the expansion rates of the bilateral ventricles, for a majority of vertices were found to follow the order of AD>MCI>HC. The bilateral hippocampus and the left amygdala were subsegmented into multiple functionally meaningful subregions with the help of high-field MRI scans. The largest group differences in localized shape atrophy rates on the hippocampus were found to occur in CA1, followed by subiculum, CA2, and finally CA3/dentate gyrus, which is consistent with the neurofibrillary tangle accumulation trajectory. Highly nonuniform group differences were detected on the amygdala; vertices on the core amygdala (basolateral and lateral nucleus) revealed much larger atrophy rates, whereas those on the noncore amygdala (mainly centromedial) displayed similar or even smaller atrophy rates in AD relative to HC. The temporal horns of the ventricles were observed to have the largest localized ventricular expansion rate differences; with the AD group showing larger localized expansion rates on the anterior horn and the body part of the ventricles as well. Significant correlations were observed between the localized shape change rates of each of these six structures and the cognitive deterioration rates as quantified by the Alzheimer's Disease Assessment Scale-Cognitive Behavior Section increase rate and the Mini Mental State Examination decrease rate.

Objective

To identify brain atrophy from structural-MRI and cerebral blood flow(CBF) patterns from arterial spin labeling perfusion-MRI that are best predictors of the Aβ-burden, measured as composite ¹⁸F-AV45-PET uptake, in individuals with early mild cognitive impairment(MCI). Furthermore, to assess the relative importance of imaging modalities in classification of Aβ+/Aβ- early mild cognitive impairment.

Methods

Sixty-seven ADNI-GO/2 participants with early-MCI were included. Voxel-wise anatomical shape variation measures were computed by estimating the initial diffeomorphic mapping momenta from an unbiased control template. CBF measures normalized to average motor cortex CBF were mapped onto the template space. Using partial least squares regression, we identified the structural and CBF signatures of Aβ after accounting for normal cofounding effects of age, sex, and education.

Results

¹⁸F-AV45-positive early-MCIs could be identified with 83% classification accuracy, 87% positive predictive value, and 84% negative predictive value by multidisciplinary classifiers combining demographics data, ApoE ε4-genotype, and a multimodal MRI-based Aβ score.

Interpretation

Multimodal-MRI can be used to predict the amyloid status of early-MCI individuals. MRI is a very attractive candidate for the identification of inexpensive and non-invasive surrogate biomarkers of Aβ deposition. Our approach is expected to have value for the identification of individuals likely to be Aβ+ in circumstances where cost or logistical problems prevent Aβ detection using cerebrospinal fluid analysis or Aβ-PET. This can also be used in clinical settings and clinical trials, aiding subject recruitment and evaluation of treatment efficacy. Imputation of the Aβ-positivity status could also complement Aβ-PET by identifying individuals who would benefit the most from this assessment.

There is an increasing need to develop and apply powerful statistical tests to detect multiple traits-single locus associations, as arising from neuroimaging genetics and other studies. For example, in the Alzheimer's Disease Neuroimaging Initiative (ADNI), in addition to genome-wide single nucleotide polymorphisms (SNPs), thousands of neuroimaging and neuropsychological phenotypes as intermediate phenotypes for Alzheimer's disease, have been collected. Although some classic methods like MANOVA and newly proposed methods may be applied, they have their own limitations. For example, MANOVA cannot be applied to binary and other discrete traits. In addition, the relationships among these methods are not well understood. Importantly, since these tests are not data adaptive, depending on the unknown association patterns among multiple traits and between multiple traits and a locus, these tests may or may not be powerful. In this paper we propose a class of data-adaptive weights and the corresponding weighted tests in the general framework of generalized estimation equations (GEE). A highly adaptive test is proposed to select the most powerful one from this class of the weighted tests so that it can maintain high power across a wide range of situations. Our proposed tests are applicable to various types of traits with or without covariates. Importantly, we also analytically show relationships among some existing and our proposed tests, indicating that many existing tests are special cases of our proposed tests. Extensive simulation studies were conducted to compare and contrast the power properties of various existing and our new methods. Finally, we applied the methods to an ADNI dataset to illustrate the performance of the methods. We conclude with the recommendation for the use of the GEE-based Score test and our proposed adaptive test for their high and complementary performance.

With the rapid growth of modern technology, many biomedical studies are being conducted to collect massive datasets with volumes of multi-modality imaging, genetic, neurocognitive, and clinical information from increasingly large cohorts. Simultaneously extracting and integrating rich and diverse heterogeneous information in neuroimaging and/or genomics from these big datasets could transform our understanding of how genetic variants impact brain structure and function, cognitive function, and brain-related disease risk across the lifespan. Such understanding is critical for diagnosis, prevention, and treatment of numerous complex brain-related disorders (e.g., schizophrenia and Alzheimer's disease). However, the development of analytical methods for the joint analysis of both high-dimensional imaging phenotypes and high-dimensional genetic data, a big data squared (BD²) problem, presents major computational and theoretical challenges for existing analytical methods. Besides the high-dimensional nature of BD², various neuroimaging measures often exhibit strong spatial smoothness and dependence and genetic markers may have a natural dependence structure arising from linkage disequilibrium. We review some recent developments of various statistical techniques for imaging genetics, including massive univariate and voxel-wise approaches, reduced rank regression, mixture models, and group sparse multi-task regression. By doing so, we hope that this review may encourage others in the statistical community to enter into this new and exciting field of research.