MRI predictors of amyloid pathology: results from the EMIF-AD Multimodal Biomarker Discovery study

Background

With the shift of research focus towards the pre-dementia stage of Alzheimer’s disease (AD), there is an urgent need for reliable, non-invasive biomarkers to predict amyloid pathology. The aim of this study was to assess whether easily obtainable measures from structural MRI, combined with demographic data, cognitive data and apolipoprotein E (APOE) ε4 genotype, can be used to predict amyloid pathology using machine-learning classification.

Methods

We examined 810 subjects with structural MRI data and amyloid markers from the European Medical Information Framework for Alzheimer’s Disease Multimodal Biomarker Discovery study, including subjects with normal cognition (CN, n = 337, age 66.5 ± 7.2, 50% female, 27% amyloid positive), mild cognitive impairment (MCI, n = 375, age 69.1 ± 7.5, 53% female, 63% amyloid positive) and AD dementia (n = 98, age 67.0 ± 7.7, 48% female, 97% amyloid positive). Structural MRI scans were visually assessed and Freesurfer was used to obtain subcortical volumes, cortical thickness and surface area measures. We first assessed univariate associations between MRI measures and amyloid pathology using mixed models. Next, we developed and tested an automated classifier using demographic, cognitive, MRI and APOE ε4 information to predict amyloid pathology. A support vector machine (SVM) with nested 10-fold cross-validation was applied to identify a set of markers best discriminating between amyloid positive and amyloid negative subjects.

Results

In univariate associations, amyloid pathology was associated with lower subcortical volumes and thinner cortex in AD-signature regions in CN and MCI. The multi-variable SVM classifier provided an area under the curve (AUC) of 0.81 ± 0.07 in MCI and an AUC of 0.74 ± 0.08 in CN. In CN, selected features for the classifier included APOE ε4, age, memory scores and several MRI measures such as hippocampus, amygdala and accumbens volumes and cortical thickness in temporal and parahippocampal regions. In MCI, the classifier including demographic and APOE ε4 information did not improve after additionally adding imaging measures.

Conclusions

Amyloid pathology is associated with changes in structural MRI measures in CN and MCI. An automated classifier based on clinical, imaging and APOE ε4 data can identify the presence of amyloid pathology with a moderate level of accuracy. These results could be used in clinical trials to pre-screen subjects for anti-amyloid therapies.

Alzheimer’s disease (AD) is characterized pathologically by beta-amyloid (Aβ) plaques and neurofibrillary tangles of misfolded tau protein [1]. As amyloid pathology may arise up to two decades before the onset of dementia, research focus has shifted towards the pre-dementia stage, which provides an opportunity for secondary prevention [2,3,4]. The design of clinical trials targeting the amyloid pathway in this early stage would be facilitated by the ability to recruit subjects with amyloid pathology. Amyloid pathology can be assessed in cerebrospinal fluid (CSF), obtainable by lumbar puncture, or on positron emission tomography (PET) scans. However, obtaining CSF is relatively invasive and PET scans are costly, invasive by exposing subjects to radiation and are not universally available. As the estimated prevalence of amyloid pathology between the ages of 60 and 80 ranges from 10 to 33% for cognitively normal (CN) subjects and from 37 to 60% for subjects with mild cognitive impairment (MCI) [5], assessing amyloid pathology with CSF or PET for screening purposes is likely inefficient. Finding minimally invasive biomarkers predicting amyloid pathology could reduce the number of invasive, costly and time-consuming measures in clinical trials.

Brain atrophy markers derived from structural magnetic resonance imaging (MRI) could serve as a potential biomarker for amyloid pathology [6,7,8,9,10,11,12]. In this study, we evaluate the use of easily obtainable MRI measures for the prediction of amyloid pathology. We included both visual rating scores, which can be easily performed in clinical settings, and quantitative measures of subcortical volumes, cortical thickness and surface area, which can be derived from freely available software and may be more sensitive than visual ratings. We first assessed univariate associations between MRI measures and amyloid pathology. Next, we used support vector machine (SVM) analysis to develop a multi-variable classifier for predicting brain amyloid pathology at a single subject level. Besides imaging measures, we also included other non-invasive measures relevant to AD in the classifier, including demographic information, cognitive testing and apolipoprotein E (APOE) ε4 genotype.

Participants

We included participants from the European Medical Information Framework for Alzheimer’s Disease Multimodal Biomarker Discovery (EMIF-AD MBD) study. The aim of this study was to discover novel diagnostic and prognostic markers for pre-dementia AD, by making use of existing data and samples [13]. The EMIF-AD MBD study pooled data of 494 CN, 526 MCI and 201 AD-dementia participants from three multicentre and eight single-centre studies. Inclusion criteria were: presence of normal cognition, MCI or a clinical diagnosis of AD-type dementia; availability of data on amyloid pathology, measured in CSF or on PET; age above 50 years; availability of MRI scans, plasma, DNA or CSF samples (at least two of the modalities); and absence of major neurological, psychiatric or somatic disorders that could cause cognitive impairment.

From the 1221 subjects included in the EMIF-AD MBD study, MRI scans of 873 subjects were contributed by the different studies (Fig. 1). Based on visual assessment, 863 MRI scans were of sufficient quality for visual rating, consisting of 365 CN, 398 MCI and 100 AD-dementia participants. Data were obtained from the following cohorts: DESCRIPA [14], EDAR [15], PharmaCog [16] and single-centre studies at VU University Medical Centre [17], San Sebastian GAP [18], University of Antwerp [19], Leuven [20], University of Lausanne [21], University of Gothenburg [22] and Barcelona IDIBAPS [23]. Each study was approved by the local medical ethics committee. Subjects had provided written informed consent at the time of inclusion in the MBD study for sharing of data, fluid samples and scans.

Number of included subjects. EMIF-AD MBD European Medical Information Framework for Alzheimer’s Disease Multimodal Biomarker Discovery, MRI magnetic resonance imaging

Full size image

Clinical and cognitive data

From all parent cohorts, clinical information and neuropsychological tests were collected centrally, harmonized, pooled and stored in an online data platform as previously described [13]. In short, all parent cohorts administered the Mini-Mental State Examination (MMSE), and performed neuropsychological testing covering various cognitive domains, although the tests used varied across the different cohorts. For the cognitive domains memory, language, attention, executive functioning and visuo-construction, one priority test was selected from each cohort (Additional file 1: Table S1) and z-scores were computed based on local normative data when available, or published normative data from healthy controls otherwise.

APOE genotyping

For the entire EMIF-AD MBD cohort, APOE genotyping data from the local genetic analyses were available for 1121 (91%) individuals. Central genetic analyses were performed at Lübeck University, Germany for 805 DNA and 148 whole blood samples. From the blood samples, DNA was extracted using the QIAamp® DNA Blood Mini Kit (QIAGEN GmbH, Hilden, Germany) resulting in 953 DNA samples, of which 926 passed quality control. Genome-wide SNP genotyping was performed using the Infinium Global Screening Array (GSA) with Shared Custom Content (Illumina Inc.). APOE genotypes were determined either directly (rs7212) or by imputation (rs429358). For 80 samples for which no local APOE genotype was available, and for 45 mismatches between local and GSA-derived genotypes, the APOE genotype was determined using TaqMan assays (ThermoFisher Scientific, Foster City, CA, USA) on a QuantStudio-12 K-Flex system. TaqMan re-genotyping confirmed 23 GSA genotypes and 21 local genotypes. For one failed sample we retained the local genotype. We classified individuals as APOE ε4 carriers or non-carriers according to their genotype status at rs429358 (C-allele = ε4).

Amyloid classification

In the current selection (n = 863), amyloid status was defined by central analysis of CSF when available (n = 510), otherwise by local amyloid PET (n = 174) or local CSF (n = 179) measures. Central CSF analysis was performed at Gothenburg University, Sweden and included Aβ_1–40 and Aβ_1–42 measured using the V-PLEX Plus Aβ Peptide Panel 1 (6E10) Kit (Meso Scale Discovery, Rockville, MD, USA), as described by the manufacturer. The central cut-off value for Aβ positivity was an Aβ_42/40 ratio < 0.061. Amyloid PET was performed in one cohort using [¹⁸F]flutemetamol according to local standardized procedures, with a standardized uptake value ratio (SUVR) cut-off value > 1.38 used for abnormality [24]. In short, SUVR images were computed from spatially normalized summed images with cerebellar grey matter as the reference region. The cut-off value was derived from an independent dataset [25] and based on the statistical difference between AD dementia patients and cognitively normal subjects [24]. Local CSF amyloid was determined according to local protocols with local cut-off values. The number of amyloid positive subjects per diagnosis per cohort is presented in Additional file 1: Table S2.

MRI acquisition

At each site, imaging was acquired according to local protocols. From each parent cohort, we centrally collected the T1-weighted images, and if available also fluid-attenuated inversion recovery (FLAIR) and susceptibility weighted images (SWI) or T2*, at the VU University Medical Center, where a visual quality check was performed. The acquired sequences and acquisition parameters for the T1-weighted scans for each cohort are presented in Additional file 1: Table S3. Usually, MRI was assessed at baseline together with baseline cognitive and amyloid measures. For 104 subjects there was more than a 1-year difference between MRI acquisition and amyloid assessment. In cases where amyloid was abnormal and acquired before MRI, this subject was included in the analysis (n = 42). In cases where amyloid was normal and acquired after MRI, this subject was included in the analysis (n = 9). All other cases were excluded (n = 53). For 99 subjects there was more than a 1-year difference between baseline cognitive assessment and MRI. For these cases, we did not use the cognitive data in the multi-variable analysis. Demographic differences between subjects who were included and excluded for differences in time between MRI and amyloid or cognitive assessment are presented in Additional file 1: Tables S4 and S5.

MRI visual rating

MRI scans with sufficient quality (n = 863) were visually rated by a single experienced rater, blinded to demographic information during rating. Medial temporal lobe atrophy (MTA) was assessed on coronal reconstructions of the T1-weighted images using a 5-point scale ranging from no atrophy (0) to end-stage atrophy (4) [26]. The MTA results from the left and right hemisphere were averaged. Global cortical atrophy (GCA) was assessed on transversal FLAIR or T1 images using a 4-point scale [27]. Posterior atrophy was assessed using a 4-point scale [28] and averaged over hemispheres. White matter hyperintensities were visually assessed on FLAIR images (n = 812) using the 4-point Fazekas scale (none, punctate, early confluent, confluent) [29]. Microbleeds were assessed on SWI and/or T2* images (n = 445) and defined as rounded hypointense homogeneous foci of up to 10 mm in diameter in the brain parenchyma. Microbleeds were dichotomized as present (≥ 1 microbleeds) or absent (0 microbleeds).

MRI quantitative analysis

Good quality 3D T1 images (n = 850) were uploaded on the N4U platform (https://neugrid4you.eu/) for automated quantitative processing. Subcortical volumes, cortical thickness and surface area measures were estimated from 3D T1 MRI using Freesurfer (version 5.3.0, https://surfer.nmr.mgh.harvard.edu) as previously described [30]. All segmentations were visually inspected. We excluded data from 20 subjects for subcortical volumes (five due to complete failure of the algorithm and 15 due to segmentation errors) and from 75 subjects for cortical thickness and surface area (five due to complete failure of the algorithm, 66 due to segmentations errors of the cortical ribbon and four for other failures). Subcortical volumes were normalized by total intracranial volume (TIV). Cortical thickness and surface area were available for 68 regions according to the Desikan–Killiany atlas implemented in Freesurfer. Additionally, we computed two AD-signature meta-ROI measures that have previously been presented in the literature: one by Dickerson et al. [10] consisting of the average cortical thickness in angular, precuneus, supramarginal, superior frontal, superior parietal, temporal pole, inferior temporal, medial temporal and inferior frontal cortex; and one by Jack et al. [31] consisting of the surface-area weighted average mean cortical thickness in entorhinal, inferior temporal, middle temporal and fusiform regions.

Statistical methods

Univariate analysis

Univariate statistical analyses were performed in R (version 3.3.1). Comparisons of clinical characteristics between amyloid positive and negative subjects within each diagnostic group were performed using independent t tests or Mann–Whitney U tests for continuous variables and chi-square tests for categorical variables. Baseline comparisons in quantitative MRI measures between groups were performed with linear mixed models (continuous outcome measures) (lme4 package, version 1.1–12; lmerTest package 2.0–36), mixed effects ordered logistic regressions (ordinal outcome measures) (ordinal package, version 2015.6–28) and mixed effects logistic regressions (dichotomous outcome measures) (lme4 package). In each model, we entered amyloid status (negative, positive) and diagnosis (CN, MCI and AD) and their interaction as fixed effects. Age (centred on mean), gender and APOE ε4 status were added as covariates. Cohort was added as a random intercept. The analyses were corrected within diagnostic group (in total 22 tests: five visual ratings, 14 subcortical volumes, three cortical thickness summary measures) for multiple hypothesis testing with the p.adjust() function using the false discovery rate, and indicated as p_FDR.

Multi-variable analysis

To find the best multi-variable predictor of amyloid pathology, we used a supervised machine-learning approach based on SVM analysis. In SVM, two classes are separated by finding a hyperplane that maximizes the margin of separation between data points of each class in a high-dimensional feature space. SVMs are used extensively in neuroimaging as they have been shown to predict outcomes with high accuracy and possess the ability to model diverse and high-dimensional data [32]. We built a classifier to separate amyloid positive from amyloid negative subjects separately in the CN and MCI subgroups and, for the sake of completeness, also in the whole sample (including CN, MCI and AD-dementia patients). To address the imbalance between the number of amyloid positive and amyloid negative subjects in each diagnostic group, we adopted the re-weighting strategy [33]. That means we adjusted weights of each SVM feature inversely proportional to amyloid positive versus negative frequencies.

Machine-learning approach

We used the python Scikit-learn library (version 0.19.1) to perform SVM classification [34]. To prevent overfitting (i.e. the classifier works perfectly on the training data, but is poorly generalizable to new data), we performed feature relevance evaluation and dimensionality reduction using a tree-based feature selection approach with a nested 10-fold cross-validation design [35, 36]. This was performed separately within each subgroup (CN, MCI and whole sample).

The nested cross-validation consists of an inner loop for model building and parameter estimation, and an outer loop for model testing. Consequently, the dataset was divided into two parts: a training plus validation subset and a test subset. In the inner loop, SVM models were trained with varying SVM hyper-parameters (i.e. cost parameters C and kernel function) based on a grid search, and a feature selection was performed using classification trees. The validation set was used to determine the SVM hyper-parameters over the grid of possible values. The performance of the resulting model, with optimized SVM hyper-parameters and features, was subsequently evaluated on the test set in the outer loop. For this outer loop, we used a 10-fold cross-validation scheme so that the data were divided into 10 equally sized parts. Nine of these were used as the training/validation set and one as the test set, and the 10 parts were permuted in each iteration of the outer loop so that each one was used for testing once. Finally, the SVM results were averaged over the 10 folds to estimate the predictive power of the proposed model on the whole dataset.

Feature selection

As the input for the classifier, we used demographic information, neuropsychological information, APOE ε4 genotype and MRI measures (visual ratings, subcortical volumes, regional cortical thickness and regional surface area measures). To combine information measured on different scales, continuous demographic and MRI measures were normalized to z-scores. In the adopted tree-based feature selection strategy, the Gini index was used to measure the relevance of each feature [37]. Features with a Gini index above the mean were kept, others were discarded. The complete list of features considered and selected, in the whole dataset and for CN and MCI separately, is reported in Additional file 1: Table S8.

Performance evaluation

To assess the performance of the classifier, we computed the averaged receiver operating characteristic (ROC) area under the curve (AUC), specificity, sensitivity and accuracy for the testing datasets. We initially maximized the Youden index, and then also explored the results when setting the sensitivity at 80%, 85%, 90%, 95% and 100%. To assess the added value of combining different sources of information, we also built classifiers including only demographic information and a single other biomarker type (neuropsychological tests, APOE ε4 genotype, MRI measures). Differences in AUC ROCs between classifiers were assessed with DeLong’s test.

Demographic and cognitive comparisons

We included 810 subjects divided over three diagnostic groups: CN (n = 337), MCI (n = 375) and AD dementia (n = 98). Within the CN group, 92 (27%) subjects were amyloid positive, in the MCI group 235 (63%) and in the AD-dementia group 95 (97%). Demographic and clinical data according to diagnosis and amyloid status are presented in Table 1. The amyloid positive MCI subjects were older and had lower cognitive scores compared to the amyloid negative MCI subjects. In CN, there were no differences in age or cognition between amyloid positive and amyloid negative subjects. Amyloid positive subjects were more often APOE ε4 carriers in both the MCI and CN groups.

Univariate association between MRI measures and amyloid pathology

Within the MCI group, subjects with amyloid pathology had higher visual rating scores of medial temporal lobe atrophy, global cortical atrophy and parietal atrophy compared to amyloid negative subjects (Table 2). There were no differences in Fazekas score or presence of microbleeds. Amyloid positive MCI subjects had statistically significantly lower bilateral hippocampus, amygdala, thalamus, left caudate and right putamen volumes, and a trend towards lower right caudate (p_uncorrected = 0.08) and bilateral accumbens (both p_uncorrected = 0.07) volumes compared to amyloid negative MCI subjects (Table 3). Amyloid positive MCI subjects also had lower whole brain average cortical thickness, as well as in the two AD-signature meta-ROIs, compared to amyloid negative MCI subjects.

In the CN group, amyloid positive subjects had statistically significantly lower right hippocampus, left amygdala, left thalamus and bilateral accumbens volumes compared to amyloid negative subjects. The effect of amyloid pathology on hippocampal volume was stronger in MCI subjects compared to CN subjects (significant interaction diagnosis × amyloid status). Amyloid positive CN subjects had lower values in the Jack AD-signature meta-ROI (p_uncorrected = 0.02), but not in the Dickerson AD-signature meta-ROI (p_uncorrected = 0.3) or whole brain average cortical thickness (p_uncorrected = 0.3) compared to amyloid negative CN subjects. There were no differences in visual rating scores between amyloid positive and amyloid negative CN subjects. All individual cortical thickness and surface area regions are presented in Additional file 1: Tables S6 and S7.

Compared to amyloid positive CN subjects, amyloid positive MCI subjects had lower bilateral hippocampal and amygdala volumes (all p < 0.001) and lower whole brain average cortical thickness (p = 0.001), as well as in the two AD-signature meta-ROIs (both p < 0.001).

Multi-variable classifier results

The features selected by the classifier in CN subjects, MCI subjects and the whole sample are presented in Additional file 1: Table S8. Across diagnoses, APOE ε4 genotype was the most important feature. Other relevant features selected across samples were age, the neuropsychological memory scores and various MRI measures such as hippocampus and amygdala volumes, as well as cortical thickness in temporal and parahippocampal regions (Fig. 2).

Freesurfer regions selected as features for the classifier in cognitively normal (top row) and mild cognitive impairment (bottom row). Colour bars represent averaged feature weight

Full size image

Combining the informative selected features in the SVM resulted in AUC = 0.81 ± 0.06 in MCI subjects, AUC = 0.74 ± 0.08 in CN subjects and AUC = 0.85 ± 0.05 in the whole sample to classify amyloid positive versus amyloid negative subjects (Fig. 3; Additional file 2: Figure S1). In MCI, the combined classifier including information from all modalities performed statistically significantly better than the classifiers based on demographic information combined with neuropsychology or imaging measures alone. The classifier including demographic variables and APOE ε4 genotype did not improve after additionally adding imaging and cognitive variables in MCI. In CN, the combined classifier including information from all modalities (demographics, cognitive, genetics and imaging) performed statistically significantly better than the classifiers including variables from only a subset of these modalities (Fig. 3; Additional file 2: Figure S2). The results from the SVM including only imaging variables are displayed in Additional file 2: Figure S3.

Classifier results. Receiver operating characteristic (ROC) curves of support vector machine classifier to predict amyloid pathology in cognitively normal (left panel) and mild cognitive impairment (right panel) subjects. Red: results from the combined classifier, including demographic information, neuropsychological tests, MRI measures and APOE ε4 genotype. Specific features selected presented in Additional file 1: Table S8. Classifier results from demographic information combined with only neuropsychology (green), or MRI measures (blue) or APOE ε4 genotype (yellow). ROC significant differences assessed with DeLong’s test. *p < 0.05, **p < 0.001, ns not significant. APOE, apolipoprotein E, MRI magnetic resonance imaging

Full size image

Table 4 presents the accuracy, sensitivity and specificity of the combined SVM in CN subjects, MCI subjects and the whole sample maximizing the Youden index, and at different levels of sensitivity. When optimizing both specificity and sensitivity, the positive predictive value for amyloid pathology was 0.84 in MCI (a 42% increase compared to the a priori probability (i.e. prevalence) of 0.59) and the negative predictive value was 0.62. In CN, the SVM obtained a positive predictive value of 0.41 (a 64% increase compared to the prevalence of 0.25) and a negative predictive value of 0.84.

In this study, we found that amyloid pathology is associated with brain atrophy in CN and MCI subjects. Using machine-learning techniques, we built a classifier based on a combination of demographic, cognitive, APOE ε4 genotype and MRI data that could predict amyloid status at single subject level with a moderate level of accuracy. The performance of the classifier was higher in MCI subjects than in CN subjects. These results are of interest for clinical trial designers who wish to recruit amyloid positive subjects for inclusion.

Our results on the association between amyloid pathology and MRI measures in MCI are in line with previous studies that also found more cortical and subcortical atrophy in amyloid positive compared to amyloid negative MCI subjects [6, 7]. In CN, amyloid pathology has previously been associated with cortical atrophy [9,10,11], and lower hippocampal volume in some studies [8, 9], but not in all [11, 38]. To capture cortical changes associated with AD, two different AD-signature meta-ROIs have been proposed in the literature [10, 31]. In MCI, both AD-signature measures were related to amyloid pathology. In CN, only the AD-signature meta-ROI by Jack et al. [31] was associated with amyloid pathology in our study, suggesting that this one is more sensitive in the early disease stage. We also found an effect of amyloid pathology on nucleus accumbens volume, which was most pronounced in CN subjects. Although nucleus accumbens volumes are not often measured in AD-related studies, it has been hypothesized that this structure could show secondary neurodegeneration in AD in response to reduced input from connections to medial-temporal lobe structures [39]. It should be noted, however, that the nucleus accumbens is a small structure, which is difficult to segment automatically. These results require further validation in future studies.

The optimal features selected in the SVM by the tree-based approach included some, but not all, of the variables that showed differences between amyloid positive and amyloid negative subjects in the univariate analyses. Similarly, some of the features selected did not show statistically significant univariate group differences, although for many a trend towards lower values in amyloid positive subjects compared to amyloid negative subjects was observed. By combining the selected features derived from demographic information, neuropsychological examination, MRI measures and APOE ε4 genotype, we were able to classify MCI and CN subjects as amyloid positive or negative with a moderate level of accuracy.

The AUC for prediction of amyloid pathology was slightly higher in the MCI group compared to the CN group, and in line with a previous study in MCI [40]. In that study, a SVM classifier to predict amyloid pathology in subjects with MCI was also developed. Using cognitive data, hippocampal volume, APOE ε4 genotype and peripheral blood protein markers from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) dataset, they obtained AUC = 0.80 for predicting amyloid pathology in subjects with MCI. In contrast to a previous study [12], we did not find that combining MRI markers with APOE ε4 genotype improved prediction of amyloid pathology in MCI over only including APOE ε4.

Our results in CN are comparable to the result from a similar study using data from the ADNI and a monocentric cohort [41]. In that study, a machine-learning-based classifier including demographic variables, APOE ε4 genotype, cognitive testing and structural MRI data reached an AUC of around 0.6 in CN subjects to predict amyloid positivity. Other studies have used combinations of demographic information, APOE ε4 genotype and cognitive testing (without imaging measures) to predict amyloid positivity in CN [42, 43]. They obtained positive predictive values of 0.65 and 0.63 for amyloid positivity, which was a 43–59% increase compared to the baseline prevalence in the cohort (0.41 and 0.44 respectively). In comparison, in our study we obtained a positive predictive value of 0.41 for amyloid pathology in CN, with a baseline prevalence of 0.25 in our cohort, which is a 64% increase in predictive value. To recruit 1000 CN subjects with amyloid pathology, using the classifier could reduce the number of subjects needing to undergo amyloid assessment from 3925 to 2439, which is a 38% decrease. Assuming a cost of €850 for the pre-screening (including MRI, APOE genotyping and cognitive testing) and €3500 for an amyloid PET scan, using the classifier for pre-screening could reduce the total screening costs by nearly €2 million in this CN population. This example is based on an optimized sum of sensitivity and specificity (Youden index). For clinical trial design, it might be more interesting to optimize the sensitivity of the classifier, which would minimize the proportion of falsely excluded amyloid positive subjects, at the cost of the positive predictive value. As can be seen in Table 4, with increasing sensitivity (and higher negative predictive value), the positive predictive value of the classifier becomes lower, which would lead to increasing costs of pre-screening.

We chose SVM as a classification method for several reasons. First, it is based on a robust strategy (i.e. maximum-margin hyper-plane), which is considered to be one of the best to reduce the prediction error in a classification task [44, 45]. Second, only few parameters need to be tuned in order to make it fully operational, making SVM relatively easy to set up and use. Finally, it is particularly well suited for the separation of two classes (in this case, amyloid positive and amyloid negative).

A strength of our study is that, unlike previous studies [40,41,42], we performed our study in a heterogeneous cohort, in which data acquisition protocols were not standardized and different MR scanners and acquisition parameters were used. In this heterogeneous cohort, we showed a similar predictive accuracy compared to prospective research cohorts, which used standardized data acquisition protocols. This highlights the robustness of our approach and suggests that the results may also be generalizable to other cohorts. This will need to be tested in future studies. Our results may be of interest for studies recruiting subjects from parent cohorts to be included in (secondary) prevention studies targeting anti-amyloid therapeutics [4]. Our findings suggest that for individuals with MCI, screening for amyloid positivity can best be done by age and APOE ε4 genotype, with limited added value of MRI. In CN, MRI measures have an added value above the other markers.

This study has some limitations. First, we used data acquired at various centres, which had different inclusion criteria for subjects and used different protocols for data collection. However, as already discussed, this also increased generalizability. Second, not everyone had the same measure of amyloid pathology. When possible, we used centralized analysis of the CSF Aβ_42/40 ratio to identify amyloid positivity, which has been shown to correlate highly with PET measures of amyloid pathology [46, 47]. For data from one cohort, we only had amyloid PET data available. Although CSF and PET measures are usually in good agreement, some studies have suggested that CSF values might become abnormal earlier than PET [48, 49]. Finally, the same dataset was used to train and test the SVM classifier. Although nested k-fold cross-validation grants good generalizability of the SVM model [36], studies in independent datasets are needed to further validate our results.

Amyloid pathology is associated with structural MRI changes in AD typical regions in CN subjects and in subjects with MCI. We developed a classifier that can predict amyloid pathology at a single subject level using a combination of easily obtainable, non-invasive measures. Our results are of interest for trial designers who intend to recruit a large number of amyloid positive subjects. Implementing pre-screening procedures consisting of simple, non-invasive tests could substantially reduce screening failure rates. In future studies, the classifier might be improved by adding data from other minimally invasive tests, such as blood proteins and genetic markers [40]. In the EMIF-AD MBD study, plasma proteomics and metabolomics, and genomics and epigenomics, will also be analysed.

AD:

Alzheimer’s disease

APOE :

Apolipoprotein E

CN:

Cognitively normal

CSF:

Cerebrospinal fluid

EMIF-AD MBD:

European Medical Information Framework for Alzheimer’s Disease Multimodal Biomarker Discovery

GCA:

Global cortical atrophy visual rating

MCI:

Mild cognitive impairment

MMSE:

Mini-Mental State Examination

MRI:

Magnetic resonance imaging

MTA:

Medial temporal lobe atrophy visual rating

PET:

Positron emission tomography

SVM:

Support vector machine

TIV:

Total intracranial volume

The authors acknowledge the contribution of the personnel of the Genomic Service Facility at the VIB-U Antwerp Center for Molecular Neurology.

Funding

The present study was conducted as part of the EMIF-AD project which has received support from the Innovative Medicines Initiative Joint Undertaking under EMIF grant agreement n° 115372, resources of which are composed of a financial contribution from the European Union’s Seventh Framework Programme (FP7/2007–2013) and an EFPIA companies’ in-kind contribution. The DESCRIPA study was funded by the European Commission within the Fifth Framework Programme (QLRT-2001-2455). The EDAR study was funded by the European Commission within the Fifth Framework Programme (contract # 37670). The VUmc Alzheimer Center is supported by Stichting Alzheimer Nederland and Stichting VUmc fonds, and the clinical database structure was developed with funding from Stichting Dioraphte. The Leuven cohort was funded by the Stichting voor Alzheimer Onderzoek (grant numbers #11020, #13007 and #15005). The GAP study is supported by grants from the Department of Economic Promotion, Rural Areas and Territorial Balance of the Provincial Government of Gipuzkoa (124/16), the Department of Health of the Basque Government (2016111096), the Carlos III Institute of Health (PI15/00919, PN de I + D + I 2013–2016), Obra Social Kutxa-Fundazioa and anonymous private donors. The Gothenburg MCI study was supported by the Sahlgrenska University Hospital, Gothenburg, Sweden. The Lausanne cohort study was supported by a grant from the Swiss National Research Foundation to JP (SNF 320030_141179). The research at VIB-CMN is funded in part by the University of Antwerp Research Fund. RV is a senior clinical investigator of the Flemish Research Foundation (FWO). CVB is partly supported by the Flemish government-initiated Flanders Impulse Program on Networks for Dementia Research (VIND) and the Methusalem Excellence Program, the Research Foundation Flanders (FWO) and the University of Antwerp Research Fund, Belgium. FB is supported by the NIHR UCLH Biomedical Research Centre. HZ is supported by the Dementia Research Institute at UCL and is a Wallenberg Academy Fellow. SJV receives research support from ZonMw. VW has received funding from the European Union’s Horizon 2020 Research and Innovation Programme under grant agreement no. 666992.

Availability of data and materials

The raw imaging data used in this study are stored on the online NeuGrid platform (https://neugrid4you.eu/). Clinical data used in this study and the derived imaging measures are stored on the online tranSMART platform from EMIF-AD. Access to the data is available upon request and after approval of each data provider, via the EMIF-AD platform.

1. Mara ten Kate and Alberto Redolfi contributed equally to this work.

Affiliations

1. Alzheimer Center & Department of Neurology, VU University Medical Center, PO Box 7057, 1007 MB, Amsterdam, the Netherlands
   • Mara ten Kate
   • , Philip Scheltens
   •  & Pieter Jelle Visser
2. Laboratory of Epidemiology & Neuroimaging, IRCCS San Giovanni di Dio Fatebenefratelli, Brescia, Italy
   • Alberto Redolfi
   • , Enrico Peira
   • , Silvia Bianchetti
   •  & Giovanni Frisoni
3. Alzheimer Centrum Limburg, Department of Psychiatry and Neuropsychology, Maastricht University, Maastricht, the Netherlands
   • Isabelle Bos
   • , Stephanie J. Vos
   • , Frans R. J. Verhey
   •  & Pieter Jelle Visser
4. University Hospital Leuven, Leuven, Belgium
   • Rik Vandenberghe
   • , Silvy Gabel
   •  & Jolien Schaeverbeke
5. Laboratory for Cognitive Neurology, Department of Neurosciences, KU Leuven, Leuven, Belgium
   • Rik Vandenberghe
   • , Silvy Gabel
   •  & Jolien Schaeverbeke
6. AP-HM, CHU Timone, CIC CPCET, Service de Pharmacologie Clinique et Pharmacovigilance, Marseille, France
   • Olivier Blin
7. Neurosciences Therapeutic Area Unit, GlaxoSmithKline R&D, Stevenage, UK
   • Jill C. Richardson
8. U1171 Inserm, CHU Lille, Degenerative and Vascular Cognitive Disorders, University of Lille, Lille, France
   • Regis Bordet
9. Sahlgrenska Academy, Institute of Neuroscience and Physiology, Section for Psychiatry and Neurochemistry, University of Gothenburg, Gothenburg, Sweden
   • Anders Wallin
   •  & Carl Eckerstrom
10. Barcelona βeta Brain Research Center, Pasqual Maragall Foundation, Barcelona, Spain
    • José Luis Molinuevo
11. Reference Center for Biological Markers of Dementia (BIODEM), Institute Born-Bunge, University of Antwerp, Antwerp, Belgium
    • Sebastiaan Engelborghs
    •  & Johannes Streffer
12. Department of Neurology and Memory Clinic, Hospital Network Antwerp (ZNA) Middelheim and Hoge Beuken, Antwerp, Belgium
    • Sebastiaan Engelborghs
13. Neurodegenerative Brain Diseases, Center for Molecular Neurology, VIB, Antwerp, Belgium
    • Christine Van Broeckhoven
14. Laboratory of Neurogenetics, Institute Born-Bunge, University of Antwerp, Antwerp, Belgium
    • Christine Van Broeckhoven
15. Department of Neurology, Center for Research and Advanced Therapies, CITA-Alzheimer Foundation, San Sebastian, Spain
    • Pablo Martinez-Lage
16. Department of Psychiatry, University Hospital of Lausanne, Lausanne, Switzerland
    • Julius Popp
17. Geriatric Psychiatry, Department of Mental Health and Psychiatry, Geneva University Hospitals, Geneva, Switzerland
    • Julius Popp
18. Memory and Dementia Center, 3rd Department of Neurology, “G Papanicolau” General Hospital, Aristotle University of Thessaloniki, Thessaloniki, Greece
    • Magdalini Tsolaki
19. University of Oxford, Oxford, UK
    • Alison L. Baird
    •  & Simon Lovestone
20. King’s College London, London, UK
    • Cristina Legido-Quigley
21. Lübeck Interdisciplinary Platform for Genome Analytics, University of Lübeck, Lubeck, Germany
    • Lars Bertram
    •  & Valerija Dobricic
22. School of Public Health, Imperial College London, London, UK
    • Lars Bertram
23. Department of Psychology, University of Oslo, Oslo, Norway
    • Lars Bertram
24. Department of Psychiatry and Neurochemistry, University of Gothenburg, Mölndal, Sweden
    • Henrik Zetterberg
25. Department of Molecular Neuroscience, UCL Institute of Neurology, Queen Square, London, UK
    • Henrik Zetterberg
26. UK Dementia Research Institute at UCL, London, UK
    • Henrik Zetterberg
27. Clinical Neurochemistry Laboratory, Sahlgrenska University Hospital, Mölndal, Sweden
    • Henrik Zetterberg
28. UCB Biopharma SPRL, Braine-l’Alleud, Belgium
    • Johannes Streffer
29. Janssen Pharmaceutical Research and Development, Titusville, NJ, USA
    • Gerald P. Novak
30. MAAT, Archamps, France
    • Jerome Revillard
31. Teva Pharmaceuticals, Inc., Malvern, PA, USA
    • Mark F. Gordon
32. Boehringer Ingelheim Pharmaceuticals, Inc., Ridgefield, CT, USA
    • Mark F. Gordon
33. Worldwide Research and Development, Pfizer Inc, Cambridge, MA, USA
    • Zhiyong Xie
34. Department of Radiology and Nuclear Medicine, VUMC, Amsterdam, the Netherlands
    • Viktor Wottschel
    •  & Frederik Barkhof
35. University of Geneva, Geneva, Switzerland
    • Giovanni Frisoni
36. Institutes of Neurology and Healthcare Engineering, UCL, London, UK
    • Frederik Barkhof

Contributions

MtK and AR provided data analyses, statistical analysis and data interpretation, and wrote the manuscript. EP provided data analyses. IB and SJV provided data management. RV, SG, JS, PS OB, JCR, RB, AW, CE, JLM, SE, CVB, PM-L, JP, MT and FRJV provided data and sample collection. ALB was responsible for plasma proteomics in EMIF-AD MBD. CL-Q was responsible for plasma metabolomics in EMIF-AD MBD. LB was responsible for genetic analyses in EMIF-AD MBD. VD was responsible for APOE analyses in EMIF-AD MBD. HZ was responsible for CSF analysis in EMIF-AD MBD. SL and JS provided study design and study coordination. SB provided data analyses. GPN is a member of the EMIF-AD MBD imaging workgroup and provided data interpretation. JR provided MRI data management. MFG and ZX are members of the EMIF-AD MBD imaging workgroup and provided data interpretation. VW provided critical revision of the manuscript. GF provided data and sample collection, and supervision of the project. PJV provided study design, study coordination and supervision of the project, and is a member of the EMIF-AD MBD imaging workgroup. FB provided study design, study coordination and supervision of the project, and is head of the EMIF-AD MBD imaging workgroup. All authors revised the manuscript and approved the final version.

Corresponding author

Correspondence to Mara ten Kate.

Ethics approval and consent to participate

The local medical ethics committee in each centre approved the study. Subjects had already provided written informed consent at the time of inclusion in the cohort for use of data, samples and scans.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests. HZ has served on scientific advisory boards of Eli Lilly and Roche Diagnostics, has received travel support from Teva and is a co-founder of Brain Biomarker Solutions in Gothenburg AB, a GU Ventures-based platform company at the University of Gothenburg. MFG’s current employer is Teva Pharmaceuticals, Inc., Malvern, PA, USA; his former employer was Boehringer Ingelheim Pharmaceuticals, Inc., Ridgefield, CT, USA. Any views expressed in this publication represent the personal opinions of the authors and not those of their respective employer. JCR is a full-time employee of GlaxoSmithkline. PM-L reports personal fees from Lilly, Axon, General Electric and Nutricia for advisory boards, and lecturing fees from Lilly, Nutricia, Piramal. RV was principal investigator of the phase 1 and 2 [¹⁸F]flutemetamol trials. RV’s institution has clinical trial agreements (RV as PI) with AbbVie, Biogen, EliLilly, Merck and Novartis, and consultancy agreements (RV as PI) with Novartis and Cytox Ltd. SL has done consultancy for Eaisi, EIP Pharma, SomaLogic, Merck and Optum Labs.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Reprints and Permissions