Spatiotemporal Propagation of the Cortical Atrophy: Population and Individual Patterns

Beyond the usual suspects: multi-factorial computational models in the search for neurodegenerative disease mechanisms

From Alzheimer's disease to amyotrophic lateral sclerosis, the molecular cascades underlying neurodegenerative disorders remain poorly understood. The clinical view of neurodegeneration is confounded by symptomatic heterogeneity and mixed pathology in almost every patient. While the underlying physiological alterations originate, proliferate, and propagate potentially decades before symptomatic onset, the complexity and inaccessibility of the living brain limit direct observation over a patient's lifespan. Consequently, there is a critical need for robust computational methods to support the search for causal mechanisms of neurodegeneration by distinguishing pathogenic processes from consequential alterations, and inter-individual variability from intra-individual progression. Recently, promising advances have been made by data-driven spatiotemporal modeling of the brain, based on in vivo neuroimaging and biospecimen markers. These methods include disease progression models comparing the temporal evolution of various biomarkers, causal models linking interacting biological processes, network propagation models reproducing the spatial spreading of pathology, and biophysical models spanning cellular- to network-scale phenomena. In this review, we discuss various computational approaches for integrating cross-sectional, longitudinal, and multi-modal data, primarily from large observational neuroimaging studies, to understand (i) the temporal ordering of physiological alterations, i(i) their spatial relationships to the brain's molecular and cellular architecture, (iii) mechanistic interactions between biological processes, and (iv) the macroscopic effects of microscopic factors. We consider the extents to which computational models can evaluate mechanistic hypotheses, explore applications such as improving treatment selection, and discuss how model-informed insights can lay the groundwork for a pathobiological redefinition of neurodegenerative disorders.

Alzheimer’s Disease from the Amyloidogenic Theory to the Puzzling Crossroads between Vascular, Metabolic and Energetic Maladaptive Plasticity

Alzheimer’s disease (AD) is a progressive and degenerative disease producing the most common type of dementia worldwide. The main pathogenetic hypothesis in recent decades has been the well-known amyloidogenic hypothesis based on the involvement of two proteins in AD pathogenesis: amyloid β (Aβ) and tau. Amyloid deposition reported in all AD patients is nowadays considered an independent risk factor for cognitive decline. Vascular damage and blood–brain barrier (BBB) failure in AD is considered a pivotal mechanism for brain injury, with increased deposition of both immunoglobulins and fibrin. Furthermore, BBB dysfunction could be an early sign of cognitive decline and the early stages of clinical AD. Vascular damage generates hypoperfusion and relative hypoxia in areas with high energy demand. Long-term hypoxia and the accumulation within the brain parenchyma of neurotoxic molecules could be seeds of a self-sustaining pathological progression. Cellular dysfunction comprises all the elements of the neurovascular unit (NVU) and neuronal loss, which could be the result of energy failure and mitochondrial impairment. Brain glucose metabolism is compromised, showing a specific region distribution. This energy deficit worsens throughout aging. Mild cognitive impairment has been reported to be associated with a glucose deficit in the entorhinal cortex and in the parietal lobes. The current aim is to understand the complex interactions between amyloid β (Aβ) and tau and elements of the BBB and NVU in the brain. This new approach aimed at the study of metabolic mechanisms and energy insufficiency due to mitochondrial impairment would allow us to define therapies aimed at predicting and slowing down the progression of AD.

A multidimensional ODE-based model of Alzheimer's disease progression

Data-driven Alzheimer's disease (AD) progression models are useful for clinical prediction, disease mechanism understanding, and clinical trial design. Most dynamic models were inspired by the amyloid cascade hypothesis and described AD progression as a linear chain of pathological events. However, the heterogeneity observed in healthy and sporadic AD populations challenged the amyloid hypothesis, and there is a need for more flexible dynamical models that accompany this conceptual shift. We present a statistical model of the temporal evolution of biomarkers and cognitive tests that allows diverse biomarker paths throughout the disease. The model consists of two elements: a multivariate dynamic model of the joint evolution of biomarkers and cognitive tests; and a clinical prediction model. The dynamic model uses a system of ordinary differential equations to jointly model the rate of change of an individual's biomarkers and cognitive tests. The clinical prediction model is an ordinal logistic model of the diagnostic label. Prognosis and time-to-onset predictions are obtained by computing the clinical label probabilities throughout the forecasted biomarker trajectories. The proposed dynamical model is interpretable, free of one-dimensional progression hypotheses or disease staging paradigms, and can account for the heterogeneous dynamics observed in sporadic AD. We developed the model using longitudinal data from the Alzheimer's Disease Neuroimaging Initiative. We illustrate the patterns of biomarker rates of change and the model performance to predict the time to conversion from MCI to dementia.

Temporal Dynamics of the Scale for the Assessment and Rating of Ataxia in Spinocerebellar Ataxias

BACKGROUND

The Scale for the Assessment and Rating of Ataxia (SARA) is the reference clinical scale to assess the severity of cerebellar ataxia. In the context of upcoming therapeutic trials, a reliable clinical outcome is needed to assess the efficiency of treatments.

OBJECTIVE

The aim is to precisely assess and compare temporal dynamics of SARA and a new f-SARA.

METHODS

We analyzed data from four cohorts (EUROSCA, RISCA, CRC-SCA, and SPATAX) comprising 1210 participants and 4092 visits. The linearity of the progression and the variability were assessed using an ordinal Bayesian mixed-effect model (Leaspy). We performed sample size calculations for therapeutic trials with different scenarios to improve the responsiveness of the scale.

RESULTS

Seven of the eight different items had a nonlinear progression. The speed of progression was different between most of the items, with an average time for a one-point increase from 3.5 years [3.4; 3.6] (median, 95% credible interval) for the fastest item to 11.4 [10.9; 12.0] years. The total SARA score had a linear progression with an average time for a one-point increase of 0.95 [0.92; 0.98] years. After removing the four last items and rescaling all items from 0 to 4, variability increased and progression was slower and thus would require a larger sample size in a future therapeutic trial.

CONCLUSION

Despite a heterogeneous temporal dynamics at the item level, the global progression of SARA was linear. Changing the initial scale deteriorates the responsiveness. This new information about the temporal dynamics of the scale should help design the outcome of future clinical trials. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Forecasting individual progression trajectories in Huntington disease enables more powered clinical trials

Variability in neurodegenerative disease progression poses great challenges for the evaluation of potential treatments. Identifying the persons who will experience significant progression in the short term is key for the implementation of trials with smaller sample sizes. We apply here disease course mapping to forecast biomarker progression for individual carriers of the pathological CAG repeat expansions responsible for Huntington disease. We used data from two longitudinal studies (TRACK-HD and TRACK-ON) to synchronize temporal progression of 15 clinical and imaging biomarkers from 290 participants with Huntington disease. We used then the resulting HD COURSE MAP to forecast clinical endpoints from the baseline data of 11,510 participants from ENROLL-HD, an external validation cohort. We used such forecasts to select participants at risk for progression and compute the power of trials for such an enriched population. HD COURSE MAP forecasts biomarkers 5 years after the baseline measures with a maximum mean absolute error of 10 points for the total motor score and 2.15 for the total functional capacity. This allowed reducing sample sizes in trial up to 50% including participants with a higher risk for progression ensuring a more homogeneous group of participants.

Multi-task longitudinal forecasting with missing values on Alzheimer's disease

BACKGROUND AND OBJECTIVE

Machine learning techniques typically used in dementia assessment are not able to learn multiple tasks jointly and deal with time-dependent heterogeneous data containing missing values. In this paper, we reformulate SSHIBA, a recently introduced Bayesian multi-view latent variable model, for jointly learning diagnosis, ventricle volume, and ADAS score in dementia on longitudinal data with missing values.

METHODS

We propose a novel Bayesian Variational inference framework capable of simultaneously imputing missing values and combining information from several views. This way, we can combine different data views from different time-points in a common latent space and learn the relationships between each time-point, using the semi-supervised formulation to fully exploit the temporal structure of the data and handle missing values. In turn, the model can combine all the available information to simultaneously model and predict multiple output variables.

RESULTS

We applied the proposed model to jointly predict diagnosis, ventricle volume, and ADAS score in dementia. The comparison of imputation strategies demonstrated the superior performance of the semi-supervised formulation of the model, improving the best baseline methods. Moreover, the performance in simultaneous prediction of diagnosis, ventricle volume, and ADAS score led to an improved prediction performance over the best baseline method.

CONCLUSIONS

The results demonstrate that the proposed SSHIBA framework can learn an excellent imputation of the missing values and outperforming the baselines while simultaneously predicting three different tasks.

Italian, European, and international neuroinformatics efforts: An overview

Neuroinformatics is a research field that focusses on software tools capable of identifying, analysing, modelling, organising and sharing multiscale neuroscience data. Neuroinformatics has exploded in the last two decades with the emergence of the Big Data phenomenon, characterised by the so-called 3Vs (volume, velocity and variety), which provided neuroscientists with an improved ability to acquire and process data faster and more cheaply thanks to technical improvements in clinical, genomic and radiological technologies. This situation has led to a 'data deluge', as neuroscientists can routinely collect more study data in a few days than they could in a year just a decade ago. To address this phenomenon, several neuroimaging-focussed neuroinformatics platforms have emerged, funded by national or transnational agencies, with the following goals: (i) development of tools for archiving and organising analytical data (XNAT, REDCap and LabKey); (ii) development of data-driven models evolving from reductionist approaches to multidimensional models (RIN, IVN, HBD, EuroPOND, E-DADS and GAAIN BRAIN); and (iii) development of e-infrastructures to provide sufficient computational power and storage resources (neuGRID, HBP-EBRAINS, LONI and CONP). Although the scenario is still fragmented, there are technological and economical attempts at both national and international levels to introduce high standards for open and Findable, Accessible, Interoperable and Reusable (FAIR) neuroscience worldwide.

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.

Progression along data-driven disease timelines is predictive of Alzheimer's disease in a population-based cohort

Data-driven disease progression models have provided important insight into the timeline of brain changes in AD phenotypes. However, their utility in predicting the progression of pre-symptomatic AD in a population-based setting has not yet been investigated. In this study, we investigated if the disease timelines constructed in a case-controlled setting, with subjects stratified according to APOE status, are generalizable to a population-based cohort, and if progression along these disease timelines is predictive of AD. Seven volumetric biomarkers derived from structural MRI were considered. We estimated APOE-specific disease timelines of changes in these biomarkers using a recently proposed method called co-initialized discriminative event-based modeling (co-init DEBM). This method can also estimate a disease stage for new subjects by calculating their position along the disease timelines. The model was trained and cross-validated on the Alzheimer's Disease Neuroimaging Initiative (ADNI) dataset, and tested on the population-based Rotterdam Study (RS) cohort. We compared the diagnostic and prognostic value of the disease stage in the two cohorts. Furthermore, we investigated if the rate of change of disease stage in RS participants with longitudinal MRI data was predictive of AD. In ADNI, the estimated disease timeslines for ϵ4 non-carriers and carriers were found to be significantly different from one another (p<0.001). The estimate disease stage along the respective timelines distinguished AD subjects from controls with an AUC of 0.83 in both APOEϵ4 non-carriers and carriers. In the RS cohort, we obtained an AUC of 0.83 and 0.85 in ϵ4 non-carriers and carriers, respectively. Progression along the disease timelines as estimated by the rate of change of disease stage showed a significant difference (p<0.005) for subjects with pre-symptomatic AD as compared to the general aging population in RS. It distinguished pre-symptomatic AD subjects with an AUC of 0.81 in APOEϵ4 non-carriers and 0.88 in carriers, which was better than any individual volumetric biomarker, or its rate of change, could achieve. Our results suggest that co-init DEBM trained on case-controlled data is generalizable to a population-based cohort setting and that progression along the disease timelines is predictive of the development of AD in the general population. We expect that this approach can help to identify at-risk individuals from the general population for targeted clinical trials as well as to provide biomarker based objective assessment in such trials.

Inter-Cohort Validation of SuStaIn Model for Alzheimer's Disease

Alzheimer’s disease (AD) is a neurodegenerative disorder which spans several years from preclinical manifestations to dementia. In recent years, interest in the application of machine learning (ML) algorithms to personalized medicine has grown considerably, and a major challenge that such models face is the transferability from the research settings to clinical practice. The objective of this work was to demonstrate the transferability of the Subtype and Stage Inference (SuStaIn) model from well-characterized research data set, employed as training set, to independent less-structured and heterogeneous test sets representative of the clinical setting. The training set was composed of MRI data of 1043 subjects from the Alzheimer’s disease Neuroimaging Initiative (ADNI), and the test set was composed of data from 767 subjects from OASIS, Pharma-Cog, and ViTA clinical datasets. Both sets included subjects covering the entire spectrum of AD, and for both sets volumes of relevant brain regions were derived from T1-3D MRI scans processed with Freesurfer v5.3 cross-sectional stream. In order to assess the predictive value of the model, subpopulations of subjects with stable mild cognitive impairment (MCI) and MCIs that progressed to AD dementia (pMCI) were identified in both sets. SuStaIn identified three disease subtypes, of which the most prevalent corresponded to the typical atrophy pattern of AD. The other SuStaIn subtypes exhibited similarities with the previously defined hippocampal sparing and limbic predominant atrophy patterns of AD. Subject subtyping proved to be consistent in time for all cohorts and the staging provided by the model was correlated with cognitive performance. Classification of subjects on the basis of a combination of SuStaIn subtype and stage, mini mental state examination and amyloid-β_1-42 cerebrospinal fluid concentration was proven to predict conversion from MCI to AD dementia on par with other novel statistical algorithms, with ROC curves that were not statistically different for the training and test sets and with area under curve respectively equal to 0.77 and 0.76. This study proves the transferability of a SuStaIn model for AD from research data to less-structured clinical cohorts, and indicates transferability to the clinical setting.

Simulating the outcome of amyloid treatments in Alzheimer's disease from imaging and clinical data

In this study, we investigate SimulAD, a novel quantitative instrument for the development of intervention strategies for disease-modifying drugs in Alzheimer's disease. SimulAD is based on the modeling of the spatio-temporal dynamics governing the joint evolution of imaging and clinical biomarkers along the history of the disease, and allows the simulation of the effect of intervention time and drug dosage on the biomarkers' progression. When applied to multi-modal imaging and clinical data from the Alzheimer's Disease Neuroimaging Initiative the method enables to generate hypothetical scenarios of amyloid lowering interventions. The results quantify the crucial role of intervention time, and provide a theoretical justification for testing amyloid modifying drugs in the pre-clinical stage. Our experimental simulations are compatible with the outcomes observed in past clinical trials, and suggest that anti-amyloid treatments should be administered at least 7 years earlier than what is currently being done in order to obtain statistically powered improvement of clinical endpoints.

AD Course Map charts Alzheimer's disease progression

Alzheimer's disease (AD) is characterized by the progressive alterations seen in brain images which give rise to the onset of various sets of symptoms. The variability in the dynamics of changes in both brain images and cognitive impairments remains poorly understood. This paper introduces AD Course Map a spatiotemporal atlas of Alzheimer's disease progression. It summarizes the variability in the progression of a series of neuropsychological assessments, the propagation of hypometabolism and cortical thinning across brain regions and the deformation of the shape of the hippocampus. The analysis of these variations highlights strong genetic determinants for the progression, like possible compensatory mechanisms at play during disease progression. AD Course Map also predicts the patient's cognitive decline with a better accuracy than the 56 methods benchmarked in the open challenge TADPOLE. Finally, AD Course Map is used to simulate cohorts of virtual patients developing Alzheimer's disease. AD Course Map offers therefore new tools for exploring the progression of AD and personalizing patients care.

Medical Informatics Platform (MIP): A Pilot Study Across Clinical Italian Cohorts

Introduction: With the shift of research focus to personalized medicine in Alzheimer's Dementia (AD), there is an urgent need for tools that are capable of quantifying a patient's risk using diagnostic biomarkers. The Medical Informatics Platform (MIP) is a distributed e-infrastructure federating large amounts of data coupled with machine-learning (ML) algorithms and statistical models to define the biological signature of the disease. The present study assessed (i) the accuracy of two ML algorithms, i.e., supervised Gradient Boosting (GB) and semi-unsupervised 3C strategy (Categorize, Cluster, Classify-CCC) implemented in the MIP and (ii) their contribution over the standard diagnostic workup. Methods: We examined individuals coming from the MIP installed across 3 Italian memory clinics, including subjects with Normal Cognition (CN, n = 432), Mild Cognitive Impairment (MCI, n = 456), and AD (n = 451). The GB classifier was applied to best discriminate the three diagnostic classes in 1,339 subjects, and the CCC strategy was used to refine the classical disease categories. Four dementia experts provided their diagnostic confidence (DC) of MCI conversion on an independent cohort of 38 patients. DC was based on clinical, neuropsychological, CSF, and structural MRI information and again with addition of the outcome from the MIP tools. Results: The GB algorithm provided a classification accuracy of 85% in a nested 10-fold cross-validation for CN vs. MCI vs. AD discrimination. Accuracy increased to 95% in the holdout validation, with the omission of each Italian clinical cohort out in turn. CCC identified five homogeneous clusters of subjects and 36 biomarkers that represented the disease fingerprint. In the DC assessment, CCC defined six clusters in the MCI population used to train the algorithm and 29 biomarkers to improve patients staging. GB and CCC showed a significant impact, evaluated as +5.99% of increment on physicians' DC. The influence of MIP on DC was rated from "slight" to "significant" in 80% of the cases. Discussion: GB provided fair results in classification of CN, MCI, and AD. CCC identified homogeneous and promising classes of subjects via its semi-unsupervised approach. We measured the effect of the MIP on the physician's DC. Our results pave the way for the establishment of a new paradigm for ML discrimination of patients who will or will not convert to AD, a clinical priority for neurology.

Statistical Disease Progression Modeling in Alzheimer Disease

Background: The characterizing symptom of Alzheimer disease (AD) is cognitive deterioration. While much recent work has focused on defining AD as a biological construct, most patients are still diagnosed, staged, and treated based on their cognitive symptoms. But the cognitive capability of a patient at any time throughout this deterioration reflects not only the disease state, but also the effect of the cognitive decline on the patient's pre-disease cognitive capability. Patients with high pre-disease cognitive capabilities tend to score better on cognitive tests that are sensitive early in disease relative to patients with low pre-disease cognitive capabilities at a similar disease stage. Thus, a single assessment with a cognitive test is often not adequate for determining the stage of an AD patient. Repeated evaluation of patients' cognition over time may improve the ability to stage AD patients, and such longitudinal assessments in combinations with biomarker assessments can help elucidate the time dynamics of biomarkers. In turn, this can potentially lead to identification of markers that are predictive of disease stage and future cognitive decline, possibly before any cognitive deficit is measurable.

Methods and Findings: This article presents a class of statistical disease progression models and applies them to longitudinal cognitive scores. These non-linear mixed-effects disease progression models explicitly model disease stage, baseline cognition, and the patients' individual changes in cognitive ability as latent variables. Maximum-likelihood estimation in these models induces a data-driven criterion for separating disease progression and baseline cognition. Applied to data from the Alzheimer's Disease Neuroimaging Initiative, the model estimated a timeline of cognitive decline that spans ~15 years from the earliest subjective cognitive deficits to severe AD dementia. Subsequent analyses demonstrated how direct modeling of latent factors that modify the observed data patterns provides a scaffold for understanding disease progression, biomarkers, and treatment effects along the continuous time progression of disease.

Conclusions: The presented framework enables direct interpretations of factors that modify cognitive decline. The results give new insights to the value of biomarkers for staging patients and suggest alternative explanations for previous findings related to accelerated cognitive decline among highly educated patients and patients on symptomatic treatments.

A survey on machine and statistical learning for longitudinal analysis of neuroimaging data in Alzheimer's disease

BACKGROUND AND OBJECTIVES

Recently, longitudinal studies of Alzheimer's disease have gathered a substantial amount of neuroimaging data. New methods are needed to successfully leverage and distill meaningful information on the progression of the disease from the deluge of available data. Machine learning has been used successfully for many different tasks, including neuroimaging related problems. In this paper, we review recent statistical and machine learning applications in Alzheimer's disease using longitudinal neuroimaging.

METHODS

We search for papers using longitudinal imaging data, focused on Alzheimer's Disease and published between 2007 and 2019 on four different search engines.

RESULTS

After the search, we obtain 104 relevant papers. We analyze their approach to typical challenges in longitudinal data analysis, such as missing data and variability in the number and extent of acquisitions.

CONCLUSIONS

Reviewed works show that machine learning methods using longitudinal data have potential for disease progression modelling and computer-aided diagnosis. We compare results and models, and propose future research directions in the field.

Predicting time to dementia using a quantitative template of disease progression

INTRODUCTION

Characterization of longitudinal trajectories of biomarkers implicated in sporadic Alzheimer's disease (AD) in decades before clinical diagnosis is important for disease prevention and monitoring.

METHODS

We used a multivariate Bayesian model to temporally align 1369 Alzheimer's disease Neuroimaging Initiative participants based on the similarity of their longitudinal biomarker measures and estimated a quantitative template of the temporal evolution of cerebrospinal fluid Aβ 1 - 42 , p- ta u 181 p , and t-tau and hippocampal volume, brain glucose metabolism, and cognitive measurements. We computed biomarker trajectories as a function of time to AD dementia and predicted AD dementia onset age in a disjoint sample.

RESULTS

Quantitative template showed early changes in verbal memory, cerebrospinal fluid Aβ1-42 and p-tau181p, and hippocampal volume. Mean error in predicted AD dementia onset age was < 1.5 years.

DISCUSSION

Our method provides a quantitative approach for characterizing the natural history of AD starting at preclinical stages despite the lack of individual-level longitudinal data spanning the entire disease timeline.

Multi-study validation of data-driven disease progression models to characterize evolution of biomarkers in Alzheimer's disease

Understanding the sequence of biological and clinical events along the course of Alzheimer's disease provides insights into dementia pathophysiology and can help participant selection in clinical trials. Our objective is to train two data-driven computational models for sequencing these events, the Event Based Model (EBM) and discriminative-EBM (DEBM), on the basis of well-characterized research data, then validate the trained models on subjects from clinical cohorts characterized by less-structured data-acquisition protocols.

Seven independent data cohorts were considered totalling 2389 cognitively normal (CN), 1424 mild cognitive impairment (MCI) and 743 Alzheimer's disease (AD) patients. The Alzheimer's Disease Neuroimaging Initiative (ADNI) data set was used as training set for the constriction of disease models while a collection of multi-centric data cohorts was used as test set for validation. Cross-sectional information related to clinical, cognitive, imaging and cerebrospinal fluid (CSF) biomarkers was used.

Event sequences obtained with EBM and DEBM showed differences in the ordering of single biomarkers but according to both the first biomarkers to become abnormal were those related to CSF, followed by cognitive scores, while structural imaging showed significant volumetric decreases at later stages of the disease progression. Staging of test set subjects based on sequences obtained with both models showed good linear correlation with the Mini Mental State Examination score (R²_EBM = 0.866; R²_DEBM = 0.906). In discriminant analyses, significant differences (p-value ≤ 0.05) between the staging of subjects from training and test sets were observed in both models. No significant difference between the staging of subjects from the training and test was observed (p-value > 0.05) when considering a subset composed by 562 subjects for which all biomarker families (cognitive, imaging and CSF) are available.

Event sequence obtained with DEBM recapitulates the heuristic models in a data-driven fashion and is clinically plausible. We demonstrated inter-cohort transferability of two disease progression models and their robustness in detecting AD phases. This is an important step towards the adoption of data-driven statistical models into clinical domain.

•

Data-driven event sequences describe evolution of relevant biomarkers in AD.

•

Agreement between event sequences and heuristic AD progression models

•

Accuracy in classifying subjects from clinical cohorts up to 91%

•

Staging of subjects and MMSE scores of individuals show linear relation.

•

Transferability of AD progression models based on research data to clinical cohorts

Flexible Locally Weighted Penalized Regression With Applications on Prediction of Alzheimer's Disease Neuroimaging Initiative's Clinical Scores

In recent years, we have witnessed the explosion of large-scale data in various fields. Classical statistical methodologies, such as linear regression or generalized linear regression, often show inadequate performance on heterogeneous data because the key homogeneity assumption fails. In this paper, we present a flexible framework to handle heterogeneous populations that can be naturally grouped into several ordered subtypes. A local model technique utilizing ordinal class labels during the training stage is proposed. We define a new "progression score" that captures the progression of ordinal classes, and use a truncated Gaussian kernel to construct the weight function in a local regression framework. Furthermore, given the weights, we apply sparse shrinkage on the local fitting to handle high dimensionality. In this way, our local model is able to conduct variable selection on each query point. Numerical studies show the superiority of our proposed method over several existing ones. Our method is also applied to the Alzheimer's Disease Neuroimaging Initiative data to make predictions on the longitudinal clinical scores based on different modalities of baseline brain image features.