A data-driven study of Alzheimer's disease related amyloid and tau pathology progression

Amyloid-β is thought to facilitate the spread of tau throughout the neocortex in Alzheimer's disease, though how this occurs is not well understood. This is because of the spatial discordance between amyloid-β, which accumulates in the neocortex, and tau, which accumulates in the medial temporal lobe during ageing. There is evidence that in some cases amyloid-β-independent tau spreads beyond the medial temporal lobe where it may interact with neocortical amyloid-β. This suggests that there may be multiple distinct spatiotemporal subtypes of Alzheimer's-related protein aggregation, with potentially different demographic and genetic risk profiles. We investigated this hypothesis, applying data-driven disease progression subtyping models to post-mortem neuropathology and in vivo PET-based measures from two large observational studies: the Alzheimer's Disease Neuroimaging Initiative (ADNI) and the Religious Orders Study and Rush Memory and Aging Project (ROSMAP). We consistently identified 'amyloid-first' and 'tau-first' subtypes using cross-sectional information from both studies. In the amyloid-first subtype, extensive neocortical amyloid-β precedes the spread of tau beyond the medial temporal lobe, while in the tau-first subtype, mild tau accumulates in medial temporal and neocortical areas prior to interacting with amyloid-β. As expected, we found a higher prevalence of the amyloid-first subtype among apolipoprotein E (APOE) ε4 allele carriers while the tau-first subtype was more common among APOE ε4 non-carriers. Within tau-first APOE ε4 carriers, we found an increased rate of amyloid-β accumulation (via longitudinal amyloid PET), suggesting that this rare group may belong within the Alzheimer's disease continuum. We also found that tau-first APOE ε4 carriers had several fewer years of education than other groups, suggesting a role for modifiable risk factors in facilitating amyloid-β-independent tau. Tau-first APOE ε4 non-carriers, in contrast, recapitulated many of the features of primary age-related tauopathy. The rate of longitudinal amyloid-β and tau accumulation (both measured via PET) within this group did not differ from normal ageing, supporting the distinction of primary age-related tauopathy from Alzheimer's disease. We also found reduced longitudinal subtype consistency within tau-first APOE ε4 non-carriers, suggesting additional heterogeneity within this group. Our findings support the idea that amyloid-β and tau may begin as independent processes in spatially disconnected regions, with widespread neocortical tau resulting from the local interaction of amyloid-β and tau. The site of this interaction may be subtype-dependent: medial temporal lobe in amyloid-first, neocortex in tau-first. These insights into the dynamics of amyloid-β and tau may inform research and clinical trials that target these pathologies.

Comparison of genetic and health risk factors for mild cognitive impairment and Alzheimer's disease between Hispanic and non-Hispanic white participants

INTRODUCTION

The influence of apolipoprotein E (APOE) genotype on mild cognitive impairment (MCI) and Alzheimer's disease (AD) is well studied in the non-Hispanic white (NHW) population but not in the Hispanic population. Additionally, health risk factors such as hypertension, stroke, and depression may also differ between the two populations.

METHODS

We combined three data sets (National Alzheimer's Coordinating Center [NACC], Alzheimer's Disease Neuroimaging Initiative [ADNI], Health and Aging Brain Study: Health Disparities [HABS-HD]) and compared risk factors for MCI and AD between Hispanic and NHW participants, with a total of 24,268 participants (11.1% Hispanic).

RESULTS

APOEε4 was associated with fewer all-cause MCI cases in Hispanic participants (Hispanic odds ratio [OR]: 1.114; NHW OR: 1.453), and APOEε2 (Hispanic OR: 1.224; NHW OR: 0.592) and depression (Hispanic OR: 2.817; NHW OR: 1.847) were associated with more AD cases in Hispanic participants.

DISCUSSION

APOEε2 may not be protective for AD in Hispanic participants and Hispanic participants with depression may face a higher risk for AD.

HIGHLIGHTS

GAAIN allows for discovery of data sets to use in secondary analyses. APOEε2 was not protective for AD in Hispanic participants. APOEε4 was associated with fewer MCI cases in Hispanic participants. Depression was associated with more AD cases in Hispanic participants.

Event-based modeling in temporal lobe epilepsy demonstrates progressive atrophy from cross-sectional data

OBJECTIVE

Recent work has shown that people with common epilepsies have characteristic patterns of cortical thinning, and that these changes may be progressive over time. Leveraging a large multicenter cross-sectional cohort, we investigated whether regional morphometric changes occur in a sequential manner, and whether these changes in people with mesial temporal lobe epilepsy and hippocampal sclerosis (MTLE-HS) correlate with clinical features.

METHODS

We extracted regional measures of cortical thickness, surface area, and subcortical brain volumes from T1-weighted (T1W) magnetic resonance imaging (MRI) scans collected by the ENIGMA-Epilepsy consortium, comprising 804 people with MTLE-HS and 1625 healthy controls from 25 centers. Features with a moderate case-control effect size (Cohen d ≥ .5) were used to train an event-based model (EBM), which estimates a sequence of disease-specific biomarker changes from cross-sectional data and assigns a biomarker-based fine-grained disease stage to individual patients. We tested for associations between EBM disease stage and duration of epilepsy, age at onset, and antiseizure medicine (ASM) resistance.

RESULTS

In MTLE-HS, decrease in ipsilateral hippocampal volume along with increased asymmetry in hippocampal volume was followed by reduced thickness in neocortical regions, reduction in ipsilateral thalamus volume, and finally, increase in ipsilateral lateral ventricle volume. EBM stage was correlated with duration of illness (Spearman ρ = .293, p = 7.03 × 10-16 ), age at onset (ρ = -.18, p = 9.82 × 10-7 ), and ASM resistance (area under the curve = .59, p = .043, Mann-Whitney U test). However, associations were driven by cases assigned to EBM Stage 0, which represents MTLE-HS with mild or nondetectable abnormality on T1W MRI.

SIGNIFICANCE

From cross-sectional MRI, we reconstructed a disease progression model that highlights a sequence of MRI changes that aligns with previous longitudinal studies. This model could be used to stage MTLE-HS subjects in other cohorts and help establish connections between imaging-based progression staging and clinical features.

pySuStaIn: a Python implementation of the Subtype and Stage Inference algorithm

Progressive disorders are highly heterogeneous. Symptom-based clinical classification of these disorders may not reflect the underlying pathobiology. Data-driven subtyping and staging of patients has the potential to disentangle the complex spatiotemporal patterns of disease progression. Tools that enable this are in high demand from clinical and treatment-development communities. Here we describe the pySuStaIn software package, a Python-based implementation of the Subtype and Stage Inference (SuStaIn) algorithm. SuStaIn unravels the complexity of heterogeneous diseases by inferring multiple disease progression patterns (subtypes) and individual severity (stages) from cross-sectional data. The primary aims of pySuStaIn are to enable widespread application and translation of SuStaIn via an accessible Python package that supports simple extension and generalization to novel modelling situations within a single, consistent architecture.

Ordinal SuStaIn: Subtype and Stage Inference for Clinical Scores, Visual Ratings, and Other Ordinal Data

Subtype and Stage Inference (SuStaIn) is an unsupervised learning algorithm that uniquely enables the identification of subgroups of individuals with distinct pseudo-temporal disease progression patterns from cross-sectional datasets. SuStaIn has been used to identify data-driven subgroups and perform patient stratification in neurodegenerative diseases and in lung diseases from continuous biomarker measurements predominantly obtained from imaging. However, the SuStaIn algorithm is not currently applicable to discrete ordinal data, such as visual ratings of images, neuropathological ratings, and clinical and neuropsychological test scores, restricting the applicability of SuStaIn to a narrower range of settings. Here we propose 'Ordinal SuStaIn', an ordinal version of the SuStaIn algorithm that uses a scored events model of disease progression to enable the application of SuStaIn to ordinal data. We demonstrate the validity of Ordinal SuStaIn by benchmarking the performance of the algorithm on simulated data. We further demonstrate that Ordinal SuStaIn out-performs the existing continuous version of SuStaIn (Z-score SuStaIn) on discrete scored data, providing much more accurate subtype progression patterns, better subtyping and staging of individuals, and accurate uncertainty estimates. We then apply Ordinal SuStaIn to six different sub-scales of the Clinical Dementia Rating scale (CDR) using data from the Alzheimer's disease Neuroimaging Initiative (ADNI) study to identify individuals with distinct patterns of functional decline. Using data from 819 ADNI1 participants we identified three distinct CDR subtype progression patterns, which were independently verified using data from 790 ADNI2 participants. Our results provide insight into patterns of decline in daily activities in Alzheimer's disease and a mechanism for stratifying individuals into groups with difficulties in different domains. Ordinal SuStaIn is broadly applicable across different types of ratings data, including visual ratings from imaging, neuropathological ratings and clinical or behavioural ratings data.

Four distinct trajectories of tau deposition identified in Alzheimer's disease

Alzheimer's disease (AD) is characterized by the spread of tau pathology throughout the cerebral cortex. This spreading pattern was thought to be fairly consistent across individuals, although recent work has demonstrated substantial variability in the population with AD. Using tau-positron emission tomography scans from 1,612 individuals, we identified 4 distinct spatiotemporal trajectories of tau pathology, ranging in prevalence from 18 to 33%. We replicated previously described limbic-predominant and medial temporal lobe-sparing patterns, while also discovering posterior and lateral temporal patterns resembling atypical clinical variants of AD. These 'subtypes' were stable during longitudinal follow-up and were replicated in a separate sample using a different radiotracer. The subtypes presented with distinct demographic and cognitive profiles and differing longitudinal outcomes. Additionally, network diffusion models implied that pathology originates and spreads through distinct corticolimbic networks in the different subtypes. Together, our results suggest that variation in tau pathology is common and systematic, perhaps warranting a re-examination of the notion of 'typical AD' and a revisiting of tau pathological staging.

Sequence of clinical and neurodegeneration events in Parkinson's disease progression

Dementia is one of the most debilitating aspects of Parkinson's disease. There are no validated biomarkers that can track Parkinson's disease progression, nor accurately identify patients who will develop dementia and when. Understanding the sequence of observable changes in Parkinson's disease in people at elevated risk for developing dementia could provide an integrated biomarker for identifying and managing individuals who will develop Parkinson's dementia. We aimed to estimate the sequence of clinical and neurodegeneration events, and variability in this sequence, using data-driven statistical modelling in two separate Parkinson's cohorts, focusing on patients at elevated risk for dementia due to their age at symptom onset. We updated a novel version of an event-based model that has only recently been extended to cope naturally with clinical data, enabling its application in Parkinson's disease for the first time. The observational cohorts included healthy control subjects and patients with Parkinson's disease, of whom those diagnosed at age 65 or older were classified as having high risk of dementia. The model estimates that Parkinson's progression in patients at elevated risk for dementia starts with classic prodromal features of Parkinson's disease (olfaction, sleep), followed by early deficits in visual cognition and increased brain iron content, followed later by a less certain ordering of neurodegeneration in the substantia nigra and cortex, neuropsychological cognitive deficits, retinal thinning in dopamine layers, and further deficits in visual cognition. Importantly, we also characterize variation in the sequence. We found consistent, cross-validated results within cohorts, and agreement between cohorts on the subset of features available in both cohorts. Our sequencing results add powerful support to the increasing body of evidence suggesting that visual processing specifically is affected early in patients with Parkinson's disease at elevated risk of dementia. This opens a route to earlier and more precise detection, as well as a more detailed understanding of the pathological mechanisms underpinning Parkinson's dementia.

A comprehensive analysis of methods for assessing polygenic burden on Alzheimer's disease pathology and risk beyond APOE

Genome-wide association studies have identified dozens of loci that alter the risk to develop Alzheimer's disease. However, with the exception of the APOE-ε4 allele, most variants bear only little individual effect and have, therefore, limited diagnostic and prognostic value. Polygenic risk scores aim to collate the disease risk distributed across the genome in a single score. Recent works have demonstrated that polygenic risk scores designed for Alzheimer's disease are predictive of clinical diagnosis, pathology confirmed diagnosis and changes in imaging biomarkers. Methodological innovations in polygenic risk modelling include the polygenic hazard score, which derives effect estimates for individual single nucleotide polymorphisms from survival analysis, and methods that account for linkage disequilibrium between genomic loci. In this work, using data from the Alzheimer's disease neuroimaging initiative, we compared different approaches to quantify polygenic disease burden for Alzheimer's disease and their association (beyond the APOE locus) with a broad range of Alzheimer's disease-related traits: cross-sectional CSF biomarker levels, cross-sectional cortical amyloid burden, clinical diagnosis, clinical progression, longitudinal loss of grey matter and longitudinal decline in cognitive function. We found that polygenic scores were associated beyond APOE with clinical diagnosis, CSF-tau levels and, to a minor degree, with progressive atrophy. However, for many other tested traits such as clinical disease progression, CSF amyloid, cognitive decline and cortical amyloid load, the additional effects of polygenic burden beyond APOE were of minor nature. Overall, polygenic risk scores and the polygenic hazard score performed equally and given the ease with which polygenic risk scores can be derived; they constitute the more practical choice in comparison with polygenic hazard scores. Furthermore, our results demonstrate that incomplete adjustment for the APOE locus, i.e. only adjusting for APOE-ε4 carrier status, can lead to overestimated effects of polygenic scores due to APOE-ε4 homozygous participants. Lastly, on many of the tested traits, the major driving factor remained the APOE locus, with the exception of quantitative CSF-tau and p-tau measures.

Modeling longitudinal imaging biomarkers with parametric Bayesian multi-task learning

Longitudinal imaging biomarkers are invaluable for understanding the course of neurodegeneration, promising the ability to track disease progression and to detect disease earlier than cross-sectional biomarkers. To properly realize their potential, biomarker trajectory models must be robust to both under-sampling and measurement errors and should be able to integrate multi-modal information to improve trajectory inference and prediction. Here we present a parametric Bayesian multi-task learning based approach to modeling univariate trajectories across subjects that addresses these criteria. Our approach learns multiple subjects' trajectories within a single model that allows for different types of information sharing, that is, coupling, across subjects. It optimizes a combination of uncoupled, fully coupled and kernel coupled models. Kernel-based coupling allows linking subjects' trajectories based on one or more biomarker measures. We demonstrate this using Alzheimer's Disease Neuroimaging Initiative (ADNI) data, where we model longitudinal trajectories of MRI-derived cortical volumes in neurodegeneration, with coupling based on APOE genotype, cerebrospinal fluid (CSF) and amyloid PET-based biomarkers. In addition to detecting established disease effects, we detect disease related changes within the insula that have not received much attention within the literature. Due to its sensitivity in detecting disease effects, its competitive predictive performance and its ability to learn the optimal parameter covariance from data rather than choosing a specific set of random and fixed effects a priori, we propose that our model can be used in place of or in addition to linear mixed effects models when modeling biomarker trajectories. A software implementation of the method is publicly available.

Making use of longitudinal information in pattern recognition

Longitudinal designs are widely used in medical studies as a means of observing within-subject changes over time in groups of subjects, thereby aiming to improve sensitivity for detecting disease effects. Paralleling an increased use of such studies in neuroimaging has been the adoption of pattern recognition algorithms for making individualized predictions of disease. However, at present few pattern recognition methods exist to make full use of neuroimaging data that have been collected longitudinally, with most methods relying instead on cross-sectional style analysis. This article presents a principal component analysis-based feature construction method that uses longitudinal high-dimensional data to improve predictive performance of pattern recognition algorithms. The method can be applied to data from a wide range of longitudinal study designs and permits an arbitrary number of time-points per subject. We apply the method to two longitudinal datasets, one containing subjects with mild cognitive impairment along with healthy controls, the other with early dementia subjects and healthy controls. Across both datasets, we show improvements in predictive accuracy relative to cross-sectional classifiers for discriminating disease subjects from healthy controls on the basis of whole-brain structural magnetic resonance image-based voxels. In addition, we can transfer longitudinal information from one set of subjects to make disease predictions in another set of subjects. The proposed method is simple and, as a feature construction method, flexible with respect to the choice of classifier and image registration algorithm. Hum Brain Mapp 37:4385-4404, 2016. © 2016 Wiley Periodicals, Inc.