Charting Disease Trajectories from Isolated REM Sleep Behavior Disorder to Parkinson's Disease

BACKGROUND

Clinical presentation and progression dynamics are variable in patients with Parkinson's disease (PD). Disease course mapping is an innovative disease modelling technique that summarizes the range of possible disease trajectories and estimates dimensions related to onset, sequence, and speed of progression of disease markers.

OBJECTIVE

To propose a disease course map for PD and investigate progression profiles in patients with or without rapid eye movement sleep behavioral disorders (RBD).

METHODS

Data of 919 PD patients and 88 isolated RBD patients from three independent longitudinal cohorts were analyzed (follow-up duration = 5.1; 95% confidence interval, 1.1-8.1] years). Disease course map was estimated by using eight clinical markers (motor and non-motor symptoms) and four imaging markers (dopaminergic denervation).

RESULTS

PD course map showed that the first changes occurred in the contralateral putamen 13 years before diagnosis, followed by changes in motor symptoms, dysautonomia, sleep-all before diagnosis-and finally cognitive decline at the time of diagnosis. The model showed earlier disease onset, earlier non-motor and later motor symptoms, more rapid progression of cognitive decline in PD patients with RBD than PD patients without RBD. This pattern was even more pronounced in patients with isolated RBD with early changes in sleep, followed by cognition and non-motor symptoms and later changes in motor symptoms.

CONCLUSIONS

Our findings are consistent with the presence of distinct patterns of progression between patients with and without RBD. Understanding heterogeneity of PD progression is key to decipher the underlying pathophysiology and select homogeneous subgroups of patients for precision medicine. © 2023 International Parkinson and Movement Disorder Society.

Association Between Diseases and Symptoms Diagnosed in Primary Care and the Subsequent Specific Risk of Multiple Sclerosis

BACKGROUND AND OBJECTIVES

Previous studies have reported a possible prodrome in multiple sclerosis (MS) defined by nonspecific symptoms including mood disorder or genitourinary symptoms and increased health care use detected several years before diagnosis. This study aimed to evaluate agnostically the associations between diseases and symptoms diagnosed in primary care and the risk of MS relative to controls and 2 other autoimmune inflammatory diseases with similar population characteristics, namely lupus and Crohn disease (CD).

METHODS

A case-control study was conducted using electronic health records from the Health Improvement Network database in the United Kingdom and France. We agnostically assessed the associations between 113 diseases and symptoms in the 5 years before and after diagnosis in patients with subsequent diagnosis of MS. Individuals with a diagnosis of MS were compared with individuals without MS and individuals with 2 other autoimmune diseases, CD and lupus.

RESULTS

The study population consisted of patients with MS (n = 20,174), patients without MS (n = 54,790), patients with CD (n = 30,477), and patients with lupus (n = 7,337). Twelve ICD-10 codes were significantly positively associated with the risk of MS compared with controls without MS. After considering ICD-10 codes suggestive of neurologic symptoms as the first diagnosis of MS, 5 ICD-10 codes remained significantly associated with MS: depression (UK: odds ratio 1.22, 95% CI 1.11-1.34), sexual dysfunction (1.47, 1.11-1.95), constipation (1.5, 1.27-1.78), cystitis (1.21, 1.05-1.39), and urinary tract infections of unspecified site (1.38, 1.18-1.61). However, none of these conditions was selectively associated with MS in comparisons with both lupus and CD. All 5 ICD-10 codes identified were still associated with MS during the 5 years after diagnosis.

DISCUSSION

We identified 5 health conditions associated with subsequent MS diagnosis, which may be considered not only prodromal but also early-stage symptoms. However, these health conditions overlap with prodrome of 2 other autoimmune diseases; hence, they lack specificity to MS.

Interaction of sex and onset site on the disease trajectory of amyotrophic lateral sclerosis

BACKGROUND

Studies showed the impact of sex and onset site (spinal or bulbar) on disease onset and survival in ALS. However, they mainly result from cross-sectional or survival analysis, and the interaction of sex and onset site on the different proxies of disease trajectory has not been fully investigated.

METHODS

We selected all patients with repeated observations in the PRO-ACT database. We divided them into four groups depending on their sex and onset site. We estimated a multivariate disease progression model, named ALS Course Map, to investigate the combined temporal changes of the four sub-scores of the revised ALS functional rating scale (ALSFRSr), the forced vital capacity (FVC), and the body mass index (BMI). We then compared the progression rate, the estimated age at onset, and the relative progression of the outcomes across each group.

RESULTS

We included 1438 patients from the PRO-ACT database. They were 51% men with spinal onset, 12% men with bulbar onset, 26% women with spinal onset, and 11% women with bulbar onset. We showed a significant influence of both sex and onset site on the ALSFRSr progression. The BMI decreased 8.9 months earlier (95% CI [3.9, 13.8]) in women than men, after correction for the onset site. Among patients with bulbar onset, FVC was impaired 2.6 months earlier (95% CI [0.6, 4.6]) in women.

CONCLUSION

Using a multivariable disease modelling approach, we showed that sex and onset site are important drivers of the progression of motor function, BMI, and FVC decline.

Impact of sex and APOE-ε4 genotype on patterns of regional brain atrophy in Alzheimer's disease and healthy aging

Alzheimer's Disease (AD) is a heterogeneous disease that disproportionately affects women and people with the APOE-ε4 susceptibility gene. We aim to describe the not-well-understood influence of both risk factors on the dynamics of brain atrophy in AD and healthy aging. Regional cortical thinning and brain atrophy were modeled over time using non-linear mixed-effect models and the FreeSurfer software with t1-MRI scans from the Alzheimer's Disease Neuroimaging Initiative (N = 1,502 subjects, 6,728 images in total). Covariance analysis was used to disentangle the effect of sex and APOE genotype on the regional onset age and pace of atrophy, while correcting for educational level. A map of the regions mostly affected by neurodegeneration is provided. Results were confirmed on gray matter density data from the SPM software. Women experience faster atrophic rates in the temporal, frontal, parietal lobes and limbic system and earlier onset in the amygdalas, but slightly later onset in the postcentral and cingulate gyri as well as all regions of the basal ganglia and thalamus. APOE-ε4 genotypes leads to earlier and faster atrophy in the temporal, frontal, parietal lobes, and limbic system in AD patients, but not in healthy patients. Higher education was found to slightly delay atrophy in healthy patients, but not for AD patients. A cohort of amyloid positive patients with MCI showed a similar impact of sex as in the healthy cohort, while APOE-ε4 showed similar associations as in the AD cohort. Female sex is as strong a risk factor for AD as APOE-ε4 genotype regarding neurodegeneration. Women experience a sharper atrophy in the later stages of the disease, although not a significantly earlier onset. These findings may have important implications for the development of targeted intervention.

Multivariate disease progression modeling with longitudinal ordinal data

Disease modeling is an essential tool to describe disease progression and its heterogeneity across patients. Usual approaches use continuous data such as biomarkers to assess progression. Nevertheless, categorical or ordinal data such as item responses in questionnaires also provide insightful information about disease progression. In this work, we propose a disease progression model for ordinal and categorical data. We built it on the principles of disease course mapping, a technique that uniquely describes the variability in both the dynamics of progression and disease heterogeneity from multivariate longitudinal data. This extension can also be seen as an attempt to bridge the gap between longitudinal multivariate models and the field of item response theory. Application to the Parkinson's progression markers initiative cohort illustrates the benefits of our approach: a fine-grained description of disease progression at the item level, as compared to the aggregated total score, together with improved predictions of the patient's future visits. The analysis of the heterogeneity across individual trajectories highlights known disease trends such as tremor dominant or postural instability and gait difficulties subtypes of Parkinson's disease.

A comparison between early presentation of dementia with Lewy Bodies, Alzheimer's disease and Parkinson's disease: evidence from routine primary care and UK Biobank data

OBJECTIVE

To simultaneously contrast prediagnostic clinical characteristics of individuals with a final diagnosis of dementia with Lewy Bodies, Parkinson's disease, Alzheimer's disease compared to controls without neurodegenerative disorders.

METHODS

Using the longitudinal THIN database in the UK, we tested the association of each neurodegenerative disorder with a selected list of symptoms and broad families of treatments, and compared the associations between disorders to detect disease-specific effects. We replicated the main findings in the UK Biobank.

RESULTS

We used data of 28,222 patients with PD, 20,214 with AD, 4,682 with DLB and 20,214 controls. All neurodegenerative disorders were significantly associated with the presence of multiple clinical characteristics before their diagnosis including sleep disorders, falls, psychiatric symptoms and autonomic dysfunctions. When comparing DLB patients with patients with PD and AD patients, falls, psychiatric symptoms and autonomic dysfunction were all more strongly associated with DLB in the five years preceding the first neurodegenerative diagnosis. The use of statins was lower in patients who developed PD and higher in patients who developed DLB compared to AD. In PD patients, the use of statins was associated with the development of dementia in the five years following PD diagnosis.

INTERPRETATION

Prediagnostic presentations of falls, psychiatric symptoms and autonomic dysfunctions were more strongly associated with DLB than PD and AD. This study also suggests that whilst several associations with medications are similar in neurodegenerative disorders, statin usage is negatively associated with Parkinson's Disease but positively with DLB and AD as well as development of dementia in PD. This article is protected by copyright. All rights reserved.

Forecasting individual progression trajectories in Alzheimer's disease

The anticipation of progression of Alzheimer's disease (AD) is crucial for evaluations of secondary prevention measures thought to modify the disease trajectory. However, it is difficult to forecast the natural progression of AD, notably because several functions decline at different ages and different rates in different patients. We evaluate here AD Course Map, a statistical model predicting the progression of neuropsychological assessments and imaging biomarkers for a patient from current medical and radiological data at early disease stages. We tested the method on more than 96,000 cases, with a pool of more than 4,600 patients from four continents. We measured the accuracy of the method for selecting participants displaying a progression of clinical endpoints during a hypothetical trial. We show that enriching the population with the predicted progressors decreases the required sample size by 38% to 50%, depending on trial duration, outcome, and targeted disease stage, from asymptomatic individuals at risk of AD to subjects with early and mild AD. We show that the method introduces no biases regarding sex or geographic locations and is robust to missing data. It performs best at the earliest stages of disease and is therefore highly suitable for use in prevention trials.

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.

Asymptotic Analysis of a Matrix Latent Decomposition Model

Matrix data sets arise in network analysis for medical applications, where each network belongs to a subject and represents a measurable phenotype. These large dimensional data are often modeled using lower-dimensional latent variables, which explain most of the observed variability and can be used for predictive purposes. In this paper, we provide asymptotic convergence guarantees for the estimation of a hierarchical statistical model for matrix data sets. It captures the variability of matrices by modeling a truncation of their eigendecomposition. We show that this model is identifiable, and that consistent Maximum A Posteriori (MAP) estimation can be performed to estimate the distribution of eigenvalues and eigenvectors. The MAP estimator is shown to be asymptotically normal for a restricted version of the model.

Benchmarking off-the-shelf statistical shape modeling tools in clinical applications

Statistical shape modeling (SSM) is widely used in biology and medicine as a new generation of morphometric approaches for the quantitative analysis of anatomical shapes. Technological advancements of in vivo imaging have led to the development of open-source computational tools that automate the modeling of anatomical shapes and their population-level variability. However, little work has been done on the evaluation and validation of such tools in clinical applications that rely on morphometric quantifications(e.g., implant design and lesion screening). Here, we systematically assess the outcome of widely used, state-of-the-art SSM tools, namely ShapeWorks, Deformetrica, and SPHARM-PDM. We use both quantitative and qualitative metrics to evaluate shape models from different tools. We propose validation frameworks for anatomical landmark/measurement inference and lesion screening. We also present a lesion screening method to objectively characterize subtle abnormal shape changes with respect to learned population-level statistics of controls. Results demonstrate that SSM tools display different levels of consistencies, where ShapeWorks and Deformetrica models are more consistent compared to models from SPHARM-PDM due to the groupwise approach of estimating surface correspondences. Furthermore, ShapeWorks and Deformetrica shape models are found to capture clinically relevant population-level variability compared to SPHARM-PDM models.

Changes in the use of psychotropic drugs during the course of Alzheimer's disease: A large-scale longitudinal study of French medical records

INTRODUCTION

We aim to understand how patients with Alzheimer's disease (AD) are treated by identifying in a longitudinal fashion the late-life changes in patients' medical history that precede and follow AD diagnosis.

METHODS

We use prescription history of 34,782 patients followed between 1996 and 2019 by French general practitioners. We compare patients with an AD diagnosis, patients with mild cognitive impairment (MCI), and patients free of mental disorders. We use a generalized mixed-effects model to study the longitudinal changes in the prescription of eight drug types for a period 15 years before diagnosis and 10 years after.

RESULTS

In the decades preceding diagnosis, we find that future AD patients are treated significantly more than MCI patients with most psychotropic drugs and that most studied drugs are increasingly prescribed with age. At the time of diagnosis, all psychotropic drugs except benzodiazepines show a significant increase in prescription, while other drugs are significantly less prescribed. In the 10 years after diagnosis, nearly all categories of drugs are less and less prescribed including antidementia drugs.

DISCUSSION

Pre-diagnosis differences between future AD patients and MCI patients may indicate that subtle cognitive changes are recognized and treated as psychiatric symptoms. The disclosure of AD diagnosis drastically changes patients' care, priority being given to the management of psychiatric symptoms. The decrease of all prescriptions in the late stages may reflect treatment discontinuation and simplification of therapeutic procedures. This study therefore provides new insights into the medical practices for management of AD.

Gaussian Graphical Model Exploration and Selection in High Dimension Low Sample Size Setting

Gaussian graphical models (GGM) are often used to describe the conditional correlations between the components of a random vector. In this article, we compare two families of GGM inference methods: the nodewise approach and the penalised likelihood maximisation. We demonstrate on synthetic data that, when the sample size is small, the two methods produce graphs with either too few or too many edges when compared to the real one. As a result, we propose a composite procedure that explores a family of graphs with a nodewise numerical scheme and selects a candidate among them with an overall likelihood criterion. We demonstrate that, when the number of observations is small, this selection method yields graphs closer to the truth and corresponding to distributions with better KL divergence with regards to the real distribution than the other two. Finally, we show the interest of our algorithm on two concrete cases: first on brain imaging data, then on biological nephrology data. In both cases our results are more in line with current knowledge in each field.

Clinica: An Open-Source Software Platform for Reproducible Clinical Neuroscience Studies

We present Clinica (www.clinica.run), an open-source software platform designed to make clinical neuroscience studies easier and more reproducible. Clinica aims for researchers to (i) spend less time on data management and processing, (ii) perform reproducible evaluations of their methods, and (iii) easily share data and results within their institution and with external collaborators. The core of Clinica is a set of automatic pipelines for processing and analysis of multimodal neuroimaging data (currently, T1-weighted MRI, diffusion MRI, and PET data), as well as tools for statistics, machine learning, and deep learning. It relies on the brain imaging data structure (BIDS) for the organization of raw neuroimaging datasets and on established tools written by the community to build its pipelines. It also provides converters of public neuroimaging datasets to BIDS (currently ADNI, AIBL, OASIS, and NIFD). Processed data include image-valued scalar fields (e.g., tissue probability maps), meshes, surface-based scalar fields (e.g., cortical thickness maps), or scalar outputs (e.g., regional averages). These data follow the ClinicA Processed Structure (CAPS) format which shares the same philosophy as BIDS. Consistent organization of raw and processed neuroimaging files facilitates the execution of single pipelines and of sequences of pipelines, as well as the integration of processed data into statistics or machine learning frameworks. The target audience of Clinica is neuroscientists or clinicians conducting clinical neuroscience studies involving multimodal imaging, and researchers developing advanced machine learning algorithms applied to neuroimaging data.

Understanding the Variability in Graph Data Sets through Statistical Modeling on the Stiefel Manifold

Network analysis provides a rich framework to model complex phenomena, such as human brain connectivity. It has proven efficient to understand their natural properties and design predictive models. In this paper, we study the variability within groups of networks, i.e., the structure of connection similarities and differences across a set of networks. We propose a statistical framework to model these variations based on manifold-valued latent factors. Each network adjacency matrix is decomposed as a weighted sum of matrix patterns with rank one. Each pattern is described as a random perturbation of a dictionary element. As a hierarchical statistical model, it enables the analysis of heterogeneous populations of adjacency matrices using mixtures. Our framework can also be used to infer the weight of missing edges. We estimate the parameters of the model using an Expectation-Maximization-based algorithm. Experimenting on synthetic data, we show that the algorithm is able to accurately estimate the latent structure in both low and high dimensions. We apply our model on a large data set of functional brain connectivity matrices from the UK Biobank. Our results suggest that the proposed model accurately describes the complex variability in the data set with a small number of degrees of freedom.

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.

Anomaly detection for the individual analysis of brain PET images

Purpose: In clinical practice, positron emission tomography (PET) images are mostly analyzed visually, but the sensitivity and specificity of this approach greatly depend on the observer's experience. Quantitative analysis of PET images would alleviate this problem by helping define an objective limit between normal and pathological findings. We present an anomaly detection framework for the individual analysis of PET images. Approach: We created subject-specific abnormality maps that summarize the pathology's topographical distribution in the brain by comparing the subject's PET image to a model of healthy PET appearance that is specific to the subject under investigation. This model was generated from demographically and morphologically matched PET scans from a control dataset. Results: We generated abnormality maps for healthy controls, patients at different stages of Alzheimer's disease and with different frontotemporal dementia syndromes. We showed that no anomalies were detected for the healthy controls and that the anomalies detected from the patients with dementia coincided with the regions where abnormal uptake was expected. We also validated the proposed framework using the abnormality maps as inputs of a classifier and obtained higher classification accuracies than when using the PET images themselves as inputs. Conclusions: The proposed method was able to automatically locate and characterize the areas characteristic of dementia from PET images. The abnormality maps are expected to (i) help clinicians in their diagnosis by highlighting, in a data-driven fashion, the pathological areas, and (ii) improve the interpretability of subsequent analyses, such as computer-aided diagnosis or spatiotemporal modeling.

Are endocasts reliable proxies for brains? A 3D quantitative comparison of the extant human brain and endocast

Endocasts (i.e., replicas of the inner surface of the bony braincase) constitute a critical proxy for qualifying and quantifying variations in brain shape and organization in extinct taxa. In the absence of brain tissues preserved in the fossil record, endocasts provide the only direct evidence of brain evolution. However, debates on whether or not information inferred from the study of endocasts reflects brain shape and organization have polarized discussions in paleoneurology since the earliest descriptions of cerebral imprints in fossil hominin crania. By means of imaging techniques (i.e., MRIs and CT scans) and 3D modelling methods (i.e., surface-based comparisons), we collected consistent morphological (i.e., shape) and structural (i.e., sulci) information on the variation patterns between the brain and the endocast based on a sample of extant human individuals (N = 5) from the 3D clinical image database of the Steve Biko Academic Hospital in Pretoria (South Africa) and the Hôpitaux Universitaires Pitié Salpêtrière in Paris (France). Surfaces of the brain and endocast of the same individual were segmented from the 3D MRIs and CT images, respectively. Sulcal imprints were automatically detected. We performed a deformation-based shape analysis to compare both the shape and the sulcal pattern of the brain and the endocast. We demonstrated that there is close correspondence in terms of morphology and organization between the brain and the corresponding endocast with the exception of the superior region. By comparatively quantifying the shape and organization of the brain and endocast, this work represents an important reference for paleoneurological studies.

A Diffeomorphic Vector Field Approach to Analyze the Thickness of the Hippocampus From 7 T MRI

OBJECTIVE

7-Tesla MRI of the hippocampus enhances the visualization of its internal substructures. Among these substructures, the cornu Ammonis and subiculum form a contiguous folded ribbon of gray matter. Here, we propose a method to analyze local thickness measurements of this ribbon.

METHODS

We introduce an original approach based upon the estimation of a diffeomorphic vector field that traverses the ribbon. The method is designed to handle specificities of the hippocampus and corresponding 7-Tesla acquisitions: highly convoluted surface, non-closed ribbon, incompletely defined inner/outer boundaries, anisotropic acquisitions. We furthermore propose to conduct group comparisons using a population template built from the central surfaces of individual subjects.

RESULTS

We first assessed the robustness of our approach to anisotropy, as well as to inter-rater variability, on a post-mortem scan and on in vivo acquisitions respectively. We then conducted a group study on a dataset of in vivo MRI from temporal lobe epilepsy (TLE) patients and healthy controls. The method detected local thinning patterns in patients, predominantly ipsilaterally to the seizure focus, which is consistent with medical knowledge.

CONCLUSION

This new technique allows measuring the thickness of the hippocampus from 7-Tesla MRI. It shows good robustness with respect to anisotropy and inter-rater variability and has the potential to detect local atrophy in patients.

SIGNIFICANCE

As 7-Tesla MRI is increasingly available, this new method may become a useful tool to study local alterations of the hippocampus in brain disorders. It is made freely available to the community (code: https://github.com/aramis-lab/hiplay7-thickness, postmortem segmentation: https://doi.org/10.5281/zenodo.3533264).

Predicting the progression of mild cognitive impairment using machine learning: A systematic, quantitative and critical review

We performed a systematic review of studies focusing on the automatic prediction of the progression of mild cognitive impairment to Alzheimer's disease (AD) dementia, and a quantitative analysis of the methodological choices impacting performance. This review included 172 articles, from which 234 experiments were extracted. For each of them, we reported the used data set, the feature types, the algorithm type, performance and potential methodological issues. The impact of these characteristics on the performance was evaluated using a multivariate mixed effect linear regressions. We found that using cognitive, fluorodeoxyglucose-positron emission tomography or potentially electroencephalography and magnetoencephalography variables significantly improved predictive performance compared to not including them, whereas including other modalities, in particular T1 magnetic resonance imaging, did not show a significant effect. The good performance of cognitive assessments questions the wide use of imaging for predicting the progression to AD and advocates for exploring further fine domain-specific cognitive assessments. We also identified several methodological issues, including the absence of a test set, or its use for feature selection or parameter tuning in nearly a fourth of the papers. Other issues, found in 15% of the studies, cast doubts on the relevance of the method to clinical practice. We also highlight that short-term predictions are likely not to be better than predicting that subjects stay stable over time. These issues highlight the importance of adhering to good practices for the use of machine learning as a decision support system for the clinical practice.

Predicting PET-derived myelin content from multisequence MRI for individual longitudinal analysis in multiple sclerosis

Multiple sclerosis (MS) is a demyelinating and inflammatory disease of the central nervous system (CNS). The demyelination process can be repaired by the generation of a new sheath of myelin around the axon, a process termed remyelination. In MS patients, the demyelination-remyelination cycles are highly dynamic. Over the years, magnetic resonance imaging (MRI) has been increasingly used in the diagnosis of MS and it is currently the most useful paraclinical tool to assess this diagnosis. However, conventional MRI pulse sequences are not specific for pathological mechanisms such as demyelination and remyelination. Recently, positron emission tomography (PET) with radiotracer [¹¹C]PIB has become a promising tool to measure in-vivo myelin content changes which is essential to push forward our understanding of mechanisms involved in the pathology of MS, and to monitor individual patients in the context of clinical trials focused on repair therapies. However, PET imaging is invasive due to the injection of a radioactive tracer. Moreover, it is an expensive imaging test and not offered in the majority of medical centers in the world. In this work, by using multisequence MRI, we thus propose a method to predict the parametric map of [¹¹C]PIB PET, from which we derived the myelin content changes in a longitudinal analysis of patients with MS. The method is based on the proposed conditional flexible self-attention GAN (CF-SAGAN) which is specifically adjusted for high-dimensional medical images and able to capture the relationships between the spatially separated lesional regions during the image synthesis process. Jointly applying the sketch-refinement process and the proposed attention regularization that focuses on the MS lesions, our approach is shown to outperform the state-of-the-art methods qualitatively and quantitatively. Specifically, our method demonstrated a superior performance for the prediction of myelin content at voxel-wise level. More important, our method for the prediction of myelin content changes in patients with MS shows similar clinical correlations to the PET-derived gold standard indicating the potential for clinical management of patients with MS.

Convolutional neural networks for classification of Alzheimer's disease: Overview and reproducible evaluation

Numerous machine learning (ML) approaches have been proposed for automatic classification of Alzheimer's disease (AD) from brain imaging data. In particular, over 30 papers have proposed to use convolutional neural networks (CNN) for AD classification from anatomical MRI. However, the classification performance is difficult to compare across studies due to variations in components such as participant selection, image preprocessing or validation procedure. Moreover, these studies are hardly reproducible because their frameworks are not publicly accessible and because implementation details are lacking. Lastly, some of these papers may report a biased performance due to inadequate or unclear validation or model selection procedures. In the present work, we aim to address these limitations through three main contributions. First, we performed a systematic literature review. We identified four main types of approaches: i) 2D slice-level, ii) 3D patch-level, iii) ROI-based and iv) 3D subject-level CNN. Moreover, we found that more than half of the surveyed papers may have suffered from data leakage and thus reported biased performance. Our second contribution is the extension of our open-source framework for classification of AD using CNN and T1-weighted MRI. The framework comprises previously developed tools to automatically convert ADNI, AIBL and OASIS data into the BIDS standard, and a modular set of image preprocessing procedures, classification architectures and evaluation procedures dedicated to deep learning. Finally, we used this framework to rigorously compare different CNN architectures. The data was split into training/validation/test sets at the very beginning and only the training/validation sets were used for model selection. To avoid any overfitting, the test sets were left untouched until the end of the peer-review process. Overall, the different 3D approaches (3D-subject, 3D-ROI, 3D-patch) achieved similar performances while that of the 2D slice approach was lower. Of note, the different CNN approaches did not perform better than a SVM with voxel-based features. The different approaches generalized well to similar populations but not to datasets with different inclusion criteria or demographical characteristics. All the code of the framework and the experiments is publicly available: general-purpose tools have been integrated into the Clinica software (www.clinica.run) and the paper-specific code is available at: https://github.com/aramis-lab/AD-DL.

Ensemble Learning of Convolutional Neural Network, Support Vector Machine, and Best Linear Unbiased Predictor for Brain Age Prediction: ARAMIS Contribution to the Predictive Analytics Competition 2019 Challenge

We ranked third in the Predictive Analytics Competition (PAC) 2019 challenge by achieving a mean absolute error (MAE) of 3.33 years in predicting age from T1-weighted MRI brain images. Our approach combined seven algorithms that allow generating predictions when the number of features exceeds the number of observations, in particular, two versions of best linear unbiased predictor (BLUP), support vector machine (SVM), two shallow convolutional neural networks (CNNs), and the famous ResNet and Inception V1. Ensemble learning was derived from estimating weights via linear regression in a hold-out subset of the training sample. We further evaluated and identified factors that could influence prediction accuracy: choice of algorithm, ensemble learning, and features used as input/MRI image processing. Our prediction error was correlated with age, and absolute error was greater for older participants, suggesting to increase the training sample for this subgroup. Our results may be used to guide researchers to build age predictors on healthy individuals, which can be used in research and in the clinics as non-specific predictors of disease status.

Early Cognitive, Structural, and Microstructural Changes in Presymptomatic C9orf72 Carriers Younger Than 40 Years

Importance

Presymptomatic carriers of chromosome 9 open reading frame 72 (C9orf72) mutation, the most frequent genetic cause of frontotemporal lobar degeneration and amyotrophic lateral sclerosis, represent the optimal target population for the development of disease-modifying drugs. Preclinical biomarkers are needed to monitor the effect of therapeutic interventions in this population.

Objectives

To assess the occurrence of cognitive, structural, and microstructural changes in presymptomatic C9orf72 carriers.

Design, Setting, and Participants

The PREV-DEMALS study is a prospective, multicenter, observational study of first-degree relatives of individuals carrying the C9orf72 mutation. Eighty-four participants entered the study between October 2015 and April 2017; 80 (95%) were included in cross-sectional analyses of baseline data. All participants underwent neuropsychological testing and magnetic resonance imaging; 63 (79%) underwent diffusion tensor magnetic resonance imaging. Gray matter volumes and diffusion tensor imaging metrics were calculated within regions of interest. Anatomical and microstructural differences between individuals who carried the C9orf72 mutation (C9+) and those who did not carry the C9orf72 mutation (C9-) were assessed using linear mixed-effects models. Data were analyzed from October 2015 to April 2017.

Main Outcomes and Measures

Differences in neuropsychological scores, gray matter volume, and white matter integrity between C9+ and C9- individuals.

Results

Of the 80 included participants, there were 41 C9+ individuals (24 [59%] female; mean [SD] age, 39.8 [11.1] years) and 39 C9- individuals (24 [62%] female; mean [SD] age, 45.2 [13.9] years). Compared with C9- individuals, C9+ individuals had lower mean (SD) praxis scores (163.4 [6.1] vs 165.3 [5.9]; P = .01) and intransitive gesture scores (34.9 [1.6] vs 35.7 [1.5]; P = .004), atrophy in 8 cortical regions of interest and in the right thalamus, and white matter alterations in 8 tracts. When restricting the analyses to participants younger than 40 years, compared with C9- individuals, C9+ individuals had lower praxis scores and intransitive gesture scores, atrophy in 4 cortical regions of interest and in the right thalamus, and white matter alterations in 2 tracts.

Conclusions and Relevance

Cognitive, structural, and microstructural alterations are detectable in young C9+ individuals. Early and subtle praxis alterations, underpinned by focal atrophy of the left supramarginal gyrus, may represent an early and nonevolving phenotype related to neurodevelopmental effects of C9orf72 mutation. White matter alterations reflect the future phenotype of frontotemporal lobar degeneration/amyotrophic lateral sclerosis, while atrophy appears more diffuse. Our results contribute to a better understanding of the preclinical phase of C9orf72 disease and of the respective contribution of magnetic resonance biomarkers.

Trial Registration

clinicaltrials.gov Identifier: NCT02590276.

Neurite density is reduced in the presymptomatic phase of C9orf72 disease

OBJECTIVE

To assess the added value of neurite orientation dispersion and density imaging (NODDI) compared with conventional diffusion tensor imaging (DTI) and anatomical MRI to detect changes in presymptomatic carriers of chromosome 9 open reading frame 72 (C9orf72) mutation.

METHODS

The PREV-DEMALS (Predict to Prevent Frontotemporal Lobar Degeneration and Amyotrophic Lateral Sclerosis) study is a prospective, multicentre, observational study of first-degree relatives of individuals carrying the C9orf72 mutation. Sixty-seven participants (38 presymptomatic C9orf72 mutation carriers (C9+) and 29 non-carriers (C9-)) were included in the present cross-sectional study. Each participant underwent one single-shell, multishell diffusion MRI and three-dimensional T1-weighted MRI. Volumetric measures, DTI and NODDI metrics were calculated within regions of interest. Differences in white matter integrity, grey matter volume and free water fraction between C9+ and C9- individuals were assessed using linear mixed-effects models.

RESULTS

Compared with C9-, C9+ demonstrated white matter abnormalities in 10 tracts with neurite density index and only 5 tracts with DTI metrics. Effect size was significantly higher for the neurite density index than for DTI metrics in two tracts. No tract had a significantly higher effect size for DTI than for NODDI. For grey matter cortical analysis, free water fraction was increased in 13 regions in C9+, whereas 11 regions displayed volumetric atrophy.

CONCLUSIONS

NODDI provides higher sensitivity and greater tissue specificity compared with conventional DTI for identifying white matter abnormalities in the presymptomatic C9orf72 carriers. Our results encourage the use of neurite density as a biomarker of the preclinical phase.

TRIAL REGISTRATION NUMBER

NCT02590276.

Spatiotemporal analysis for detection of pre-symptomatic shape changes in neurodegenerative diseases: Initial application to the GENFI cohort

Brain atrophy as measured from structural MR images, is one of the primary imaging biomarkers used to track neurodegenerative disease progression. In diseases such as frontotemporal dementia or Alzheimer's disease, atrophy can be observed in key brain structures years before any clinical symptoms are present. Atrophy is most commonly captured as volume change of key structures and the shape changes of these structures are typically not analysed despite being potentially more sensitive than summary volume statistics over the entire structure. In this paper we propose a spatiotemporal analysis pipeline based on Large Diffeomorphic Deformation Metric Mapping (LDDMM) to detect shape changes from volumetric MRI scans. We applied our framework to a cohort of individuals with genetic variants of frontotemporal dementia and healthy controls from the Genetic FTD Initiative (GENFI) study. Our method, take full advantage of the LDDMM framework, and relies on the creation of a population specific average spatiotemporal trajectory of a relevant brain structure of interest, the thalamus in our case. The residuals from each patient data to the average spatiotemporal trajectory are then clustered and studied to assess when presymptomatic mutation carriers differ from healthy control subjects. We found statistical differences in shape in the anterior region of the thalamus at least five years before the mutation carrier subjects develop any clinical symptoms. This region of the thalamus has been shown to be predominantly connected to the frontal lobe, consistent with the pattern of cortical atrophy seen in the disease.

Fluid-attenuated inversion recovery MRI synthesis from multisequence MRI using three-dimensional fully convolutional networks for multiple sclerosis

Multiple sclerosis (MS) is a white matter (WM) disease characterized by the formation of WM lesions, which can be visualized by magnetic resonance imaging (MRI). The fluid-attenuated inversion recovery (FLAIR) MRI pulse sequence is used clinically and in research for the detection of WM lesions. However, in clinical settings, some MRI pulse sequences could be missed because of various constraints. The use of the three-dimensional fully convolutional neural networks is proposed to predict FLAIR pulse sequences from other MRI pulse sequences. In addition, the contribution of each input pulse sequence is evaluated with a pulse sequence-specific saliency map. This approach is tested on a real MS image dataset and evaluated by comparing this approach with other methods and by assessing the lesion contrast in the synthetic FLAIR pulse sequence. Both the qualitative and quantitative results show that this method is competitive for FLAIR synthesis.

Reduction of recruitment costs in preclinical AD trials: validation of automatic pre-screening algorithm for brain amyloidosis

We propose a method for recruiting asymptomatic Amyloid positive individuals in clinical trials, using a two-step process. We first select during a pre-screening phase a subset of individuals which are more likely to be amyloid positive based on the automatic analysis of data acquired during routine clinical practice, before doing a confirmatory PET-scan to these selected individuals only. This method leads to an increased number of recruitments and to a reduced number of PET-scans, resulting in a decrease in overall recruitment costs. We validate our method on three different cohorts, and consider five different classification algorithms for the pre-screening phase. We show that the best results are obtained using solely cognitive, genetic and socio-demographic features, as the slight increased performance when using MRI or longitudinal data is balanced by the cost increase they induce. We show that the proposed method generalizes well when tested on an independent cohort, and that the characteristics of the selected set of individuals are identical to the characteristics of a population selected in a standard way. The proposed approach shows how Machine Learning can be used effectively in practice to optimize recruitment costs in clinical trials.

An Automated Pipeline for the Analysis of PET Data on the Cortical Surface

We present a fully automatic pipeline for the analysis of PET data on the cortical surface. Our pipeline combines tools from FreeSurfer and PETPVC, and consists of (i) co-registration of PET and T1-w MRI (T1) images, (ii) intensity normalization, (iii) partial volume correction, (iv) robust projection of the PET signal onto the subject's cortical surface, (v) spatial normalization to a template, and (vi) atlas statistics. We evaluated the performance of the proposed workflow by performing group comparisons and showed that the approach was able to identify the areas of hypometabolism characteristic of different dementia syndromes: Alzheimer's disease (AD) and both the semantic and logopenic variants of primary progressive aphasia. We also showed that these results were comparable to those obtained with a standard volume-based approach. We then performed individual classifications and showed that vertices can be used as features to differentiate cognitively normal and AD subjects. This pipeline is integrated into Clinica, an open-source software platform for neuroscience studies available at www.clinica.run.

Double Diffeomorphism: Combining Morphometry and Structural Connectivity Analysis

The brain is composed of several neural circuits which may be seen as anatomical complexes composed of grey matter structures interconnected by white matter tracts. Grey and white matter components may be modeled as 3-D surfaces and curves, respectively. Neurodevelopmental disorders involve morphological and organizational alterations which cannot be jointly captured by usual shape analysis techniques based on single diffeomorphisms. We propose a new deformation scheme, called double diffeomorphism, which is a combination of two diffeomorphisms. The first one captures changes in structural connectivity, whereas the second one recovers the global morphological variations of both grey and white matter structures. This deformation model is integrated into a Bayesian framework for atlas construction. We evaluate it on a data-set of 3-D structures representing the neural circuits of patients with Gilles de la Tourette syndrome (GTS). We show that this approach makes it possible to localise, quantify, and easily visualise the pathological anomalies altering the morphology and organization of the neural circuits. Furthermore, results also indicate that the proposed deformation model better discriminates between controls and GTS patients than a single diffeomorphism.

Reproducible evaluation of classification methods in Alzheimer's disease: Framework and application to MRI and PET data

A large number of papers have introduced novel machine learning and feature extraction methods for automatic classification of Alzheimer's disease (AD). However, while the vast majority of these works use the public dataset ADNI for evaluation, they are difficult to reproduce because different key components of the validation are often not readily available. These components include selected participants and input data, image preprocessing and cross-validation procedures. The performance of the different approaches is also difficult to compare objectively. In particular, it is often difficult to assess which part of the method (e.g. preprocessing, feature extraction or classification algorithms) provides a real improvement, if any. In the present paper, we propose a framework for reproducible and objective classification experiments in AD using three publicly available datasets (ADNI, AIBL and OASIS). The framework comprises: i) automatic conversion of the three datasets into a standard format (BIDS); ii) a modular set of preprocessing pipelines, feature extraction and classification methods, together with an evaluation framework, that provide a baseline for benchmarking the different components. We demonstrate the use of the framework for a large-scale evaluation on 1960 participants using T1 MRI and FDG PET data. In this evaluation, we assess the influence of different modalities, preprocessing, feature types (regional or voxel-based features), classifiers, training set sizes and datasets. Performances were in line with the state-of-the-art. FDG PET outperformed T1 MRI for all classification tasks. No difference in performance was found for the use of different atlases, image smoothing, partial volume correction of FDG PET images, or feature type. Linear SVM and L2-logistic regression resulted in similar performance and both outperformed random forests. The classification performance increased along with the number of subjects used for training. Classifiers trained on ADNI generalized well to AIBL and OASIS. All the code of the framework and the experiments is publicly available: general-purpose tools have been integrated into the Clinica software (www.clinica.run) and the paper-specific code is available at: https://gitlab.icm-institute.org/aramislab/AD-ML.

Spatiotemporal Propagation of the Cortical Atrophy: Population and Individual Patterns

Repeated failures in clinical trials for Alzheimer’s disease (AD) have raised a strong interest for the prodromal phase of the disease. A better understanding of the brain alterations during this early phase is crucial to diagnose patients sooner, to estimate an accurate disease stage, and to give a reliable prognosis. According to recent evidence, structural alterations in the brain are likely to be sensitive markers of the disease progression. Neuronal loss translates in specific spatiotemporal patterns of cortical atrophy, starting in the enthorinal cortex and spreading over other cortical regions according to specific propagation pathways. We developed a digital model of the cortical atrophy in the left hemisphere from prodromal to diseased phases, which is built on the temporal alignment and combination of several short-term observation data to reconstruct the long-term history of the disease. The model not only provides a description of the spatiotemporal patterns of cortical atrophy at the group level but also shows the variability of these patterns at the individual level in terms of difference in propagation pathways, speed of propagation, and age at propagation onset. Longitudinal MRI datasets of patients with mild cognitive impairments who converted to AD are used to reconstruct the cortical atrophy propagation across all disease stages. Each observation is considered as a signal spatially distributed on a network, such as the cortical mesh, each cortex location being associated to a node. We consider how the temporal profile of the signal varies across the network nodes. We introduce a statistical mixed-effect model to describe the evolution of the cortex alterations. To ensure a spatiotemporal smooth propagation of the alterations, we introduce a constrain on the propagation signal in the model such that neighboring nodes have similar profiles of the signal changes. Our generative model enables the reconstruction of personalized patterns of the neurodegenerative spread, providing a way to estimate disease progression stages and predict the age at which the disease will be diagnosed. The model shows that, for instance, APOE carriers have a significantly higher pace of cortical atrophy but not earlier atrophy onset.

Revolution of Alzheimer Precision Neurology. Passageway of Systems Biology and Neurophysiology

The Precision Neurology development process implements systems theory with system biology and neurophysiology in a parallel, bidirectional research path: a combined hypothesis-driven investigation of systems dysfunction within distinct molecular, cellular, and large-scale neural network systems in both animal models as well as through tests for the usefulness of these candidate dynamic systems biomarkers in different diseases and subgroups at different stages of pathophysiological progression. This translational research path is paralleled by an "omics"-based, hypothesis-free, exploratory research pathway, which will collect multimodal data from progressing asymptomatic, preclinical, and clinical neurodegenerative disease (ND) populations, within the wide continuous biological and clinical spectrum of ND, applying high-throughput and high-content technologies combined with powerful computational and statistical modeling tools, aimed at identifying novel dysfunctional systems and predictive marker signatures associated with ND. The goals are to identify common biological denominators or differentiating classifiers across the continuum of ND during detectable stages of pathophysiological progression, characterize systems-based intermediate endophenotypes, validate multi-modal novel diagnostic systems biomarkers, and advance clinical intervention trial designs by utilizing systems-based intermediate endophenotypes and candidate surrogate markers. Achieving these goals is key to the ultimate development of early and effective individualized treatment of ND, such as Alzheimer's disease. The Alzheimer Precision Medicine Initiative (APMI) and cohort program (APMI-CP), as well as the Paris based core of the Sorbonne University Clinical Research Group "Alzheimer Precision Medicine" (GRC-APM) were recently launched to facilitate the passageway from conventional clinical diagnostic and drug development toward breakthrough innovation based on the investigation of the comprehensive biological nature of aging individuals. The APMI movement is gaining momentum to systematically apply both systems neurophysiology and systems biology in exploratory translational neuroscience research on ND.

The endocranial shape of Australopithecus africanus: surface analysis of the endocasts of Sts 5 and Sts 60

Assessment of global endocranial morphology and regional neuroanatomical changes in early hominins is critical for the reconstruction of evolutionary trajectories of cerebral regions in the human lineage. Early evidence of cortical reorganization in specific local areas (e.g. visual cortex, inferior frontal gyrus) is perceptible in the non-human South African hominin fossil record. However, to date, little information is available regarding potential global changes in the early hominin brain. The introduction of non-invasive imaging techniques opens up new perspectives for the study of hominin brain evolution. In this context, our primary aim in this study is to explore the organization of the Australopithecus africanus endocasts, and highlight the nature and extent of the differences distinguishing A. africanus from the extant hominids at both local and global scales. By means of X-ray-based imaging techniques, we investigate two A. africanus specimens from Sterkfontein Member 4, catalogued as Sts 5 and Sts 60, respectively a complete cranium and a partial cranial endocast. Endocrania were virtually reconstructed and compared by using a landmark-free registration method based on smooth and invertible surface deformation. Both local and global information provided by our deformation-based approach are used to perform statistical analyses and topological mapping of inter-specific variation. Statistical analyses indicate that the endocranial shape of Sts 5 and Sts 60 approximates the Pan condition. Furthermore, our study reveals substantial differences with respect to the extant human condition, particularly in the parietal regions. Compared with Pan, the endocranial shape of the fossil specimens differs in the anterior part of the frontal gyri.