Towards cascading genetic risk in Alzheimer's disease.

Alzheimer's disease typically progresses in stages, which have been defined by the presence of disease-specific biomarkers: amyloid (A), tau (T) and neurodegeneration (N). This progression of biomarkers has been condensed into the ATN framework, in which each of the biomarkers can be either positive (+) or negative (-). Over the past decades, genome-wide association studies have implicated ∼90 different loci involved with the development of late-onset Alzheimer's disease. Here, we investigate whether genetic risk for Alzheimer's disease contributes equally to the progression in different disease stages or whether it exhibits a stage-dependent effect. Amyloid (A) and tau (T) status was defined using a combination of available PET and CSF biomarkers in the Alzheimer's Disease Neuroimaging Initiative cohort. In 312 participants with biomarker-confirmed A-T- status, we used Cox proportional hazards models to estimate the contribution of APOE and polygenic risk scores (beyond APOE) to convert to A+T- status (65 conversions). Furthermore, we repeated the analysis in 290 participants with A+T- status and investigated the genetic contribution to conversion to A+T+ (45 conversions). Both survival analyses were adjusted for age, sex and years of education. For progression from A-T- to A+T-, APOE-e4 burden showed a significant effect [hazard ratio (HR) = 2.88; 95% confidence interval (CI): 1.70-4.89; P < 0.001], whereas polygenic risk did not (HR = 1.09; 95% CI: 0.84-1.42; P = 0.53). Conversely, for the transition from A+T- to A+T+, the contribution of APOE-e4 burden was reduced (HR = 1.62; 95% CI: 1.05-2.51; P = 0.031), whereas the polygenic risk showed an increased contribution (HR = 1.73; 95% CI: 1.27-2.36; P < 0.001). The marginal APOE effect was driven by e4 homozygotes (HR = 2.58; 95% CI: 1.05-6.35; P = 0.039) as opposed to e4 heterozygotes (HR = 1.74; 95% CI: 0.87-3.49; P = 0.12). The genetic risk for late-onset Alzheimer's disease unfolds in a disease stage-dependent fashion. A better understanding of the interplay between disease stage and genetic risk can lead to a more mechanistic understanding of the transition between ATN stages and a better understanding of the molecular processes leading to Alzheimer's disease, in addition to opening therapeutic windows for targeted interventions.

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.

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.

Robust Markers and Sample Sizes for Multicenter Trials of Huntington Disease.

OBJECTIVE: The identification of sensitive biomarkers is essential to validate therapeutics for Huntington disease (HD). We directly compare structural imaging markers across the largest collective imaging HD dataset to identify a set of imaging markers robust to multicenter variation and to derive upper estimates on sample sizes for clinical trials in HD. METHODS: We used 1 postprocessing pipeline to retrospectively analyze T1-weighted magnetic resonance imaging (MRI) scans from 624 participants at 3 time points, from the PREDICT-HD, TRACK-HD, and IMAGE-HD studies. We used mixed effects models to adjust regional brain volumes for covariates, calculate effect sizes, and simulate possible treatment effects in disease-affected anatomical regions. We used our model to estimate the statistical power of possible treatment effects for anatomical regions and clinical markers. RESULTS: We identified a set of common anatomical regions that have similarly large standardized effect sizes (>0.5) between healthy control and premanifest HD (PreHD) groups. These included subcortical, white matter, and cortical regions and nonventricular cerebrospinal fluid (CSF). We also observed a consistent spatial distribution of effect size by region across the whole brain. We found that multicenter studies were necessary to capture treatment effect variance; for a 20% treatment effect, power of >80% was achieved for the caudate (n = 661), pallidum (n = 687), and nonventricular CSF (n = 939), and, crucially, these imaging markers provided greater power than standard clinical markers. INTERPRETATION: Our findings provide the first cross-study validation of structural imaging markers in HD, supporting the use of these measurements as endpoints for both observational studies and clinical trials. ANN NEUROL 2020;87:751-762.

Serial CT analysis in idiopathic pulmonary fibrosis: comparison of visual features that determine patient outcome.

AIMS: Patients with idiopathic pulmonary fibrosis (IPF) receiving antifibrotic medication and patients with non-IPF fibrosing lung disease often demonstrate rates of annualised forced vital capacity (FVC) decline within the range of measurement variation (5.0%-9.9%). We examined whether change in visual CT variables could help confirm whether marginal FVC declines represented genuine clinical deterioration rather than measurement noise. METHODS: In two IPF cohorts (cohort 1: n=103, cohort 2: n=108), separate pairs of radiologists scored paired volumetric CTs (acquired between 6 and 24 months from baseline). Change in interstitial lung disease, honeycombing, reticulation, ground-glass opacity extents and traction bronchiectasis severity was evaluated using a 5-point scale, with mortality prediction analysed using univariable and multivariable Cox regression analyses. Both IPF populations were then combined to determine whether change in CT variables could predict mortality in patients with marginal FVC declines. RESULTS: On univariate analysis, change in all CT variables except ground-glass opacity predicted mortality in both cohorts. On multivariate analysis adjusted for patient age, gender, antifibrotic use and baseline disease severity (diffusing capacity for carbon monoxide), change in traction bronchiectasis severity predicted mortality independent of FVC decline. Change in traction bronchiectasis severity demonstrated good interobserver agreement among both scorer pairs. Across all study patients with marginal FVC declines, change in traction bronchiectasis severity independently predicted mortality and identified more patients with deterioration than change in honeycombing extent. CONCLUSIONS: Change in traction bronchiectasis severity is a measure of disease progression that could be used to help resolve the clinical importance of marginal FVC declines.

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.

Effect of Genetic Risk on the Relationship Between rs-fMRI Complexity and Tau and Amyloid PET in Alzheimer's Disease.

Reduced functional magnetic resonance imaging (fMRI)-complexity in Alzheimer's disease (AD) progression has been demonstrated and found to be associated with tauopathy and cognition. However, association of fMRI-complexity with amyloid and influence of genetic risk (APOEɛ4) remain unknown. Here we investigate the association between fMRI-complexity, tau-PET, and amyloid-PET as well as influence of APOE genotype using multivariate generalized linear models. We show that fMRI-complexity has a strong association with tau but not amyloid deposition and that the presence of an APOEɛ4 allele enhances this effect. Thus fMRI-complexity provides a surrogate marker of impaired brain functionality in AD progression.

Investigating the factors that explain white matter hyperintensity load in older Indians.

White matter hyperintensities are areas of hyperintense signal on MRI that typically represent cerebrovascular pathology. While focal white matter hyperintensities are common among older individuals, extensive white matter hyperintensities have been found to accelerate the progression of dementia. However, little is currently known about how various socioeconomic, health, lifestyle and environmental factors affect the severity of these lesions, particularly in low- and middle-income countries such as India. We investigated this question using cross-sectional MRI data (n = 126) from a pilot neuroimaging sub-study of an ongoing, nationally representative epidemiological study of late-life cognition in India. As a screening step, we estimated white matter hyperintensity load from fluid-attenuated inversion recovery MRI using a fully automated technique and tested for associations with each factor separately, controlling for age, sex and estimated total intracranial volume in each case. A combined model of white matter hyperintensity load included five factors which were significant after multiple comparisons correction: systolic blood pressure, body mass index, urbanicity status (urban versus rural living), daily chore hours and the frequency of store trips. This model explained an additional 27% of the variance in white matter hyperintensity load (54 versus 27% for the baseline model with only age, sex and estimated total intracranial volume). We accounted for the possibility of reverse causality by additionally controlling for concurrent markers of neurodegeneration and cognitive impairment, with no substantial change in our findings. Overall, our findings suggest that controlling high blood pressure and maintaining both a healthy body mass index and high levels of physical activity may reduce white matter hyperintensity load in older Indian adults, helping to prevent or delay dementia.

The temporal event-based model: Learning event timelines in progressive diseases.

Timelines of events, such as symptom appearance or a change in biomarker value, provide powerful signatures that characterise progressive diseases. Understanding and predicting the timing of events is important for clinical trials targeting individuals early in the disease course when putative treatments are likely to have the strongest effect. However, previous models of disease progression cannot estimate the time between events and provide only an ordering in which they change. Here, we introduce the temporal event-based model (TEBM), a new probabilistic model for inferring timelines of biomarker events from sparse and irregularly sampled datasets. We demonstrate the power of the TEBM in two neurodegenerative conditions: Alzheimer's disease (AD) and Huntington's disease (HD). In both diseases, the TEBM not only recapitulates current understanding of event orderings but also provides unique new ranges of timescales between consecutive events. We reproduce and validate these findings using external datasets in both diseases. We also demonstrate that the TEBM improves over current models; provides unique stratification capabilities; and enriches simulated clinical trials to achieve a power of 80% with less than half the cohort size compared with random selection. The application of the TEBM naturally extends to a wide range of progressive conditions.

Deep phenotyping and genomic data from a nationally representative study on dementia in India.

The Harmonized Diagnostic Assessment of Dementia for the Longitudinal Aging Study in India (LASI-DAD) is a nationally representative in-depth study of cognitive aging and dementia. We present a publicly available dataset of harmonized cognitive measures of 4,096 adults 60 years of age and older in India, collected across 18 states and union territories. Blood samples were obtained to carry out whole blood and serum-based assays. Results are included in a venous blood specimen datafile that can be linked to the Harmonized LASI-DAD dataset. A global screening array of 960 LASI-DAD respondents is also publicly available for download, in addition to neuroimaging data on 137 LASI-DAD participants. Altogether, these datasets provide comprehensive information on older adults in India that allow researchers to further understand risk factors associated with cognitive impairment and dementia.

Nomograms of human hippocampal volume shifted by polygenic scores.

Nomograms are important clinical tools applied widely in both developing and aging populations. They are generally constructed as normative models identifying cases as outliers to a distribution of healthy controls. Currently used normative models do not account for genetic heterogeneity. Hippocampal volume (HV) is a key endophenotype for many brain disorders. Here, we examine the impact of genetic adjustment on HV nomograms and the translational ability to detect dementia patients. Using imaging data from 35,686 healthy subjects aged 44-82 from the UK Biobank (UKB), we built HV nomograms using Gaussian process regression (GPR), which - compared to a previous method - extended the application age by 20 years, including dementia critical age ranges. Using HV polygenic scores (HV-PGS), we built genetically adjusted nomograms from participants stratified into the top and bottom 30% of HV-PGS. This shifted the nomograms in the expected directions by ~100 mm3 (2.3% of the average HV), which equates to 3 years of normal aging for a person aged ~65. Clinical impact of genetically adjusted nomograms was investigated by comparing 818 subjects from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database diagnosed as either cognitively normal (CN), having mild cognitive impairment (MCI) or Alzheimer's disease (AD) patients. While no significant change in the survival analysis was found for MCI-to-AD conversion, an average of 68% relative decrease was found in intra-diagnostic-group variance, highlighting the importance of genetic adjustment in untangling phenotypic heterogeneity.

Association of Amygdala Development With Different Forms of Anxiety in Autism Spectrum Disorder.

BACKGROUND: The amygdala is widely implicated in both anxiety and autism spectrum disorder. However, no studies have investigated the relationship between co-occurring anxiety and longitudinal amygdala development in autism. Here, the authors characterize amygdala development across childhood in autistic children with and without traditional DSM forms of anxiety and anxieties distinctly related to autism. METHODS: Longitudinal magnetic resonance imaging scans were acquired at up to four time points for 71 autistic and 55 typically developing (TD) children (∼2.5-12 years, 411 time points). Traditional DSM anxiety and anxieties distinctly related to autism were assessed at study time 4 (∼8-12 years) using a diagnostic interview tailored to autism: the Anxiety Disorders Interview Schedule-IV with the Autism Spectrum Addendum. Mixed-effects models were used to test group differences at study time 1 (3.18 years) and time 4 (11.36 years) and developmental differences (age-by-group interactions) in right and left amygdala volume between autistic children with and without DSM or autism-distinct anxieties and TD children. RESULTS: Autistic children with DSM anxiety had significantly larger right amygdala volumes than TD children at both study time 1 (5.10% increase) and time 4 (6.11% increase). Autistic children with autism-distinct anxieties had significantly slower right amygdala growth than TD, autism-no anxiety, and autism-DSM anxiety groups and smaller right amygdala volumes at time 4 than the autism-no anxiety (-8.13% decrease) and autism-DSM anxiety (-12.05% decrease) groups. CONCLUSIONS: Disparate amygdala volumes and developmental trajectories between DSM and autism-distinct forms of anxiety suggest different biological underpinnings for these common, co-occurring conditions in autism.

Spatial-Temporal Patterns of ß-Amyloid Accumulation: A Subtype and Stage Inference Model Analysis.

BACKGROUND AND OBJECTIVES: ß-amyloid (Aß) staging models assume a single spatial-temporal progression of amyloid accumulation. We assessed evidence for Aß accumulation subtypes by applying the data-driven Subtype and Stage Inference (SuStaIn) model to amyloid-PET data. METHODS: Amyloid-PET data of 3,010 participants were pooled from 6 cohorts (ALFA+, EMIF-AD, ABIDE, OASIS, and ADNI). Standardized uptake value ratios were calculated for 17 regions. We applied the SuStaIn algorithm to identify consistent subtypes in the pooled dataset based on the cross-validation information criterion and the most probable subtype/stage classification per scan. The effects of demographics and risk factors on subtype assignment were assessed using multinomial logistic regression. RESULTS: Participants were mostly cognitively unimpaired (n = 1890 [62.8%]), had a mean age of 68.72 (SD 9.1) years, 42.1% were APOE ε4 carriers, and 51.8% were female. A 1-subtype model recovered the traditional amyloid accumulation trajectory, but SuStaIn identified 3 optimal subtypes, referred to as frontal, parietal, and occipital based on the first regions to show abnormality. Of the 788 (26.2%) with strong subtype assignment (>50% probability), the majority was assigned to frontal (n = 415 [52.5%]), followed by parietal (n = 199 [25.3%]) and occipital subtypes (n = 175 [22.2%]). Significant differences across subtypes included distinct proportions of APOE ε4 carriers (frontal 61.8%, parietal 57.1%, occipital 49.4%), participants with dementia (frontal 19.7%, parietal 19.1%, occipital 31.0%), and lower age for the parietal subtype (frontal/occipital 72.1 years, parietal 69.3 years). Higher amyloid (Centiloid) and CSF p-tau burden was observed for the frontal subtype; parietal and occipital subtypes did not differ. At follow-up, most participants (81.1%) maintained baseline subtype assignment and 25.6% progressed to a later stage. DISCUSSION: Whereas a 1-trajectory model recovers the established pattern of amyloid accumulation, SuStaIn determined that 3 subtypes were optimal, showing distinct associations with Alzheimer disease risk factors. Further analyses to determine clinical utility are warranted.

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.

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.