- Research
- Open access
- Published:
Amyloid-β misfolding as a plasma biomarker indicates risk for future clinical Alzheimer’s disease in individuals with subjective cognitive decline
Alzheimer's Research & Therapy volume 12, Article number: 169 (2020)
Abstract
Background
We evaluated Aβ misfolding in combination with Aβ42/40 ratio as a prognostic tool for future clinical progression to mild cognitive impairment (MCI) or dementia due to Alzheimer’s disease (AD) in individuals with subjective cognitive decline (SCD).
Methods
Baseline plasma samples (n = 203) from SCD subjects in the SCIENCe project and Amsterdam Dementia Cohort (age 61 ± 9 years; 57% male, mean follow-up time 2.7 years) were analyzed using immuno-infrared-sensor technology. Within 6 years of follow-up, 22 (11%) individuals progressed to MCI or dementia due to AD. Sensor readout values > 1646 cm− 1 reflected normal Aβ folding; readouts at ≤ 1646 cm− 1 reflected low and at < 1644 cm− 1 high misfolding. We used Cox proportional hazard models to quantify Aβ misfolding as a prognostic biomarker for progression to MCI and dementia due to AD. The accuracy of the predicted development of MCI/AD was determined by time-dependent receiver operating characteristic (t-ROC) curve analyses that take individual follow-up and conversion times into account. Statistical models were adjusted for age, sex, and APOEε4 status. Additionally, plasma Aβ42/40 data measured by SIMOA were statistically analyzed and compared.
Results
All 22 patients who converted to MCI or AD-dementia within 6 years exhibited Aβ misfolding at baseline. Cox analyses revealed a hazard ratio (HR) of 19 (95% confidence interval [CI] 2.2–157.8) for future conversion of SCD subjects with high misfolding and of 11 (95% CI 1.0–110.1) for those with low misfolding. T-ROC curve analyses yielded an area under the curve (AUC) of 0.94 (95% CI 0.86–1.00; 6-year follow-up) for Aβ misfolding in an age, sex, and APOEε4 model. A similar model with plasma Aβ42/40 ratio yielded an AUC of 0.92 (95% CI, 0.82–1.00). The AUC increased to 0.99 (95% CI, 0.99–1.00) after inclusion of both Aβ misfolding and the Aβ42/40 ratio.
Conclusions
A panel of structure- and concentration-based plasma amyloid biomarkers may predict conversion to clinical MCI and dementia due to AD in cognitively unimpaired subjects. These plasma biomarkers provide a noninvasive and cost-effective alternative for screening early AD pathological changes. Follow-up studies and external validation in larger cohorts are in progress for further validation of our findings.
Background
Biomarkers indicating Alzheimer’s disease (AD) in cognitively unimpaired individuals are essential for future therapeutic approaches [1]. In clinical trials, amyloid-β (Aβ) positron emission tomography (Aβ-PET) is used to visualize amyloid in the brain [2,3,4,5,6]. Aβ-PET targets fibrillary amyloid plaques; however, the technique is costly [3]. Since the increased deposition of Aβ42 in amyloid plaques is strongly associated with its decreased CSF concentration [7], the Aβ42 concentration in CSF is widely used as a fluid biomarker for AD [8,9,10,11,12]. Unfortunately, CSF necessitates invasive lumbar puncture, thereby limiting its use [13]. Therefore, there is an urgent need to identify blood-based biomarkers for AD. After several years of controversial results in cohort studies [14], the utility of the Aβ42/40 ratio as a noninvasive blood plasma biomarker has recently shown promising results, as it correlates with brain amyloid pathology [15,16,17]. Additionally, recent studies have suggested prognostic value for plasma Aβ42/40 [16, 18,19,20]. This is due to the novel methodological improvements that allow sensitive analysis of plasma Aβ concentrations by immunoassays such as SIMOA [16, 21] or mass spectrometry [15, 17, 22].
Complementary to these measurements of Aβ concentrations in plasma indicating indirectly the amyloid plaque formation, Aβ misfolding in blood plasma might serve as an additional biomarker for an early disease stage [23,24,25,26]. In vitro, it was shown that plaques are formed by initial misfolding of Aβ from a predominantly monomeric alpha-helical and disordered structure to a β-sheet-enriched secondary structure. This structural change triggers and initiates Aβ oligomerization and aggregation to much larger fibrils on the nanometer scale [27,28,29,30,31,32]. Since misfolding of Aβ causing peptide aggregation and plaque formation is believed to be one of the initial events in AD development and starts 15–20 years before clinical symptoms occur [4, 27, 33], it is conceivable that Aβ misfolding might be one of the earliest detectable events in AD pathogenesis. Therefore, Aβ misfolding might be a promising risk marker to identify high-risk individuals in the very early stages of the disease.
Aβ misfolding can be monitored using an immuno-infrared sensor (iRS), which measures the frequency of the C=O stretching vibration of the Aβ backbone [23, 34]. This vibration causes the amide I absorbance band, which in turn gives information about the secondary structure distribution of all Aβ isoforms [23].
We have recently validated Aβ misfolding as a structure biomarker in plasma for probable AD (prospective Essen cohort) [24], for prodromal AD (BioFINDER) [25], and for preclinical cognitively unimpaired AD subjects (population-based ESTHER cohort) [25]. Including APOEε4 as a risk factor, early AD pathology could be identified with an AUC of 0.84 as early as 14 years before the clinical diagnosis of AD [35]. Remarkably, consistent cutoffs have been observed across all studies with a general threshold of < 1644 cm− 1, indicating high misfolding in individuals with dementia due to AD. A second upper threshold of > 1646 cm− 1 was recently introduced to further differentiate between already low misfolding and a “normal” Aβ secondary structure distribution, as observed in individuals without dementia due to AD [26].
Here, we analyzed baseline plasma samples of individuals with subjective cognitive decline (SCD) from the memory clinic-based SCIENCe study and the Amsterdam Dementia Cohort, to further explore the potential of the structure-based plasma biomarker and the prognostic value in cognitively unimpaired individuals [16, 36,37,38]. However, SCD individuals have a higher risk of progression to dementia compared to cognitively unimpaired individuals without subjective complaints [39]. Furthermore, we explored, whether the combination of structure-based and concentration-based Aβ42/40 plasma biomarkers could increase the prognostic performance in a panel including other covariates to predict clinical progression to MCI or AD-dementia.
Methods
Study cohort
A total of 236 baseline EDTA plasma samples were received from individuals in the ongoing Amsterdam Dementia Cohort and SCIENCe project [36,37,38], which was a subset of the study population described in detail by Verberk et al. [16]. The inclusion criteria for the current study comprised a baseline diagnosis of SCD received within 0.5 years of plasma sampling, the plasma sample had to be nonhemolytic, at least one follow-up visit was available, and follow-up diagnosis confirmed either SCD or conversions to MCI or AD-dementia. Eight subjects progressed to non-AD-dementia and were excluded. Only the subjects with baseline data on possible covariates, such as age, sex, and APOEε4 status, were included. The total number of subjects diagnosed with SCD eligible for inclusion in the current study was 203.
All subjects were referred to the Alzheimer Center, Amsterdam between 2000 and 2016 because of cognitive complaints, where they were thoroughly screened for their neurological, physical, and neuropsychological functioning. Furthermore, brain magnetic resonance imaging, electroencephalography, Aβ1–42 CSF biomarker analysis by Innotest ELISAs (Fuijirebio, Ghent, Belgium), and APOEε4 genetic screening by polymerase chain reaction were performed. CSF Aβ42 levels were dichotomized into positive (< 813 pg/mL) or negative (> 813 pg/mL) [40]. In a multidisciplinary consensus meeting, the label “SCD” was assigned, since the subject’s report of cognitive worsening could not be objectified by any of the clinical or cognitive tests performed, and the criteria for MCI, dementia, or other medical conditions possibly explaining the perceived cognitive decline were not met [41]. Written informed consent to use medical data and biomaterials for research purposes was in place, in accordance with ethical approval from the VU University Medical Center based on the Helsinki Declaration (seventh revision).
Clinical follow-up
After the baseline visit, individuals were regularly examined, including repeated neurological, physical, and neuropsychological assessments. Clinical diagnosis was re-evaluated after each visit by clinical consensus using valid diagnostic criteria [16, 36,37,38]. Clinical progression was defined as the change in diagnosis to MCI [42, 43] or AD-dementia [44, 45], and the time point of progression was defined as the visit date corresponding to when the diagnosis first changed; when subjects first progressed to MCI and later to dementia, the time of the MCI diagnosis was used as the time point of clinical progression.
Blood plasma collection, processing, and preanalytics
Blood for preparation of EDTA plasma was drawn through venipuncture at the Amsterdam UMC, VU University Medical Center. Plasma tubes were centrifuged at 1800×g for 10 min within 1 h of collection, and plasma was aliquoted in 0.5-mL portions in polypropylene tubes and stored at − 80 °C. Samples were first thawed at the Amsterdam UMC for the quantitative determination of plasma concentrations of Aβ1–40 and Aβ1–42 with the SIMOA Human Neurology 3-Plex A assay kit (Quanterix, Lexington, MA) using the SIMOA HD-1 analyzer [16]. Thereafter, 200 μl of plasma per subject was shipped from Amsterdam to Bochum on dry ice.
Antibody and plasma preparation
For the iRS, we used the monoclonal antibody HJ5.1 as capture antibody, which specifically binds to all structural isoforms of Aβ (see supplementary Fig. S-1). The antibody HJ5.1 was developed, produced, and validated by Holtzman et al. from Washington University (Washington, USA). The antibody HJ5.1 has also been used for Aβ extraction in other assays by different groups [22]. We received the hybridoma cell line for in-house production of HJ5.1. Because of the low human anti-mouse antibody and interference of HJ5.1 with other human IgGs in our assay, we established a standardized sample preprocessing protocol for IgG removal from the plasma sample. In general, the depletion procedure comprised two different strategies for IgG removal. In strategy 1, the Pierce™ Spin column (0.5 ml volume, Thermo Scientific, Germany) was charged with 500 μl of Pierce™ Protein G Agarose resin (Thermo Scientific, Germany) and blocked with 150 μl of 1% (w/v) casein solution for 2 h. Afterwards, the column was washed with 3 ml of phosphate-buffered saline (PBS; pH 7.4) buffer. Subsequently, 100 μl of plasma sample was loaded on the prepared column and incubated overnight at 4 to 6 °C under light shaking conditions. As a final step, the remaining IgGs were eliminated from the plasma sample with magnetic beads. To this end, the plasma sample was eluted from the column and spiked with magnetic beads that are covalently linked with Protein A and Protein A/G and Protein G. The admit volume of the magnetic beads was 50 μl of each kind of bead. The plasma sample was incubated for at least 3 h under light shaking at 4–6 °C. Afterwards, the magnetic beads were completely removed from the sample via a magnetic stand, and the sample was immediately measured. To ensure that there is no antibody interference signal during sample analysis in our assay, all IgG-depleted plasma samples (strategy 1) simultaneously underwent quality control. Therefore, the surface of a separate sensor element was covered with antibody binding protein A employing silane chemistry. Such a sensor surface is highly specific to detect the lowest amounts of residual IgG interference within the plasma sample. The procedure of this sensor modification is described in detail by Budde et al. [46]. If no antibody interference was detected, the quality control confirmed the complete IgG removal from the plasma sample; otherwise, if strategy 1 was insufficient to completely remove IgGs from the sample and antibody interference was observed by two characteristic antibody absorbance peaks at 1638 cm− 1 and 1685 cm− 1, we extended IgG removal of this plasma sample by an alternative double column depletion strategy (strategy 2). Therefore, the first step of the depletion process was performed via Pierce™ Protein G Agarose resin as described above (strategy 1) followed by a second column charged with 500 μl of Pierce™ Protein A Agarose resin (Thermo Scientific, Germany). The protein A agarose column was incubated overnight at 4 to 6 °C under light shaking conditions.
To validate that the depletion process had no effect on the secondary structure of Aβ and on the diagnostic performance of the iRS, we compared the results of single- and double-depleted plasma samples measured with HJ 5.1 with each other and with previous measurements from these patient samples, in which the monoclonal antibody A8978 (Sigma Aldrich, Germany) was used [24, 25]. Since this antibody does not show any interference with human antibodies, we do not have to perform IgG depletion of the samples. Therefore, these measurements can serve as controls to indicate whether the depletion affects the Aβ secondary structure distribution. Our results with antibody HJ5.1 were highly reproducible (Fig. S-2), independent of the depletion strategy, and highly comparable to the diagnostic results obtained using antibody A8978 (Fig. S-3). Hence, the depletion procedure for IgG removal from the plasma samples had no effect on the secondary structure of Aβ in our assay. Both IgG removal procedures guaranteed the elimination of any antibody interference signal in our assay. Importantly, strategy 2 also had no influence on the Aβ secondary structure distribution within the required timeline (supplementary Fig. S-2). Samples were immediately measured after depletion to avoid further sample handling and additional freeze-thaw cycles. However, sample age in general showed no influence on the iRS readout (supplementary Fig. S-4). If sufficient, strategy 1 was preferred because of time and cost constraints as well as to avoid the hurdles of sample preparation. In the present study, strategy 1 was employed for 80% of the samples, and strategy 2 was employed for 20% of the samples.
Plasma analysis by the immuno-infrared sensor (iRS)
Baseline blood plasma samples were analyzed with our iRS [23] as schematically outlined in Fig. 1. The iRS has been validated in detail, including NHS silane generation and characterization, antibody batch-to-batch variation, antibody performance with synthetic Aβ, and standard reference CSF and blood plasma samples, matrix effects, lower and upper limits of quantification, assay selectivity, sample handling and documentation of zero background signals after Aβ immunodepletion. These procedures and protocols have been previously described in detail [23]. In short, the monoclonal antibody HJ5.1 (D. Holtzman, Washington) was covalently immobilized on the sensor surface by silane chemistry followed by surface blocking via a casein solution. By using a flow cuvette that housed the functionalized sensor surface, the plasma samples were circulated in a constant flow (1 ml/min) for 1 h to extract the total Aβ fraction. This step was followed by a washing step for 30 min with PBS buffer to remove unbound substances from the sensor surface. The binding signal of the sample stayed stable over circulation and washing time and did not result in any change in the Aβ structure during the analysis procedure (supplementary Fig. S-5). Importantly, the antibody HJ5.1 simultaneously binds to all structural isoforms of Aβ, including monomeric, oligomeric, and fibrillary Aβ. By difference spectroscopy, only the absorbance band of bound Aβ is elucidated, and the absorbance bands of all other components are subtracted. As a readout, the infrared frequency of the amide I absorbance band of the bound Aβ fraction was recorded. The amide I band records the absorbance of the secondary structure-sensitive C=O stretching vibration of the Aβ peptide backbone. Monomeric Aβ isoforms with alpha-helical and unstructured secondary structures show a maximum frequency at 1655 wavenumbers (cm− 1), whereas β-sheet enriched structures have a maximum at 1624 cm− 1 [23]. Misfolding from predominantly monomeric alpha-helical and disordered Aβ in healthy individuals to increased β-sheet enriched structures in individuals with AD shifts the frequency down. iRS readout values < 1644 cm− 1 were indicative of increased high misfolding, and values > 1646 cm− 1 reflected normal folding. Hence, values between ≥ 1644 and ≤ 1646 cm− 1 indicated a slightly increased misfolding level, defined as low misfolding. The threshold of 1644 cm− 1 (± 1 cm− 1) was empirically determined in [25]. The second upper threshold at 1646 cm− 1 (± 1 cm− 1) was recently introduced in Nabers et al. in the context of a two-step diagnostic workflow [26]. With this threshold, we identified individuals with a largely increased likelihood of disease onset. Plasma samples were analyzed at baseline, during which all individuals were clinically diagnosed as cognitively unimpaired.
Statistical analysis
In the statistical analyses, we did not differentiate between subjects who converted to either MCI or dementia due to AD and included them in one converter group. For the description of groups, mean ± SD were used for continuous variables and absolute number and percentage for categorical variables. Comparisons between groups were performed with nonparametric tests (Wilcoxon rank-sum test). p values less than 0.05 were considered significant. Since the patient-related follow-up period differed between the subjects within this longitudinal study, we applied statistical models that take event times (including censoring) into account. We used the Cox proportional hazard model, both unadjusted and adjusted for relevant covariates (sex, age, and APOEε4), to calculate the risk of clinical disease progression. Additionally, we used time-dependent receiver operating characteristic (t-ROC) curves to calculate the diagnostic accuracy for non-converted SCD vs. follow-up MCI or AD-dementia based on baseline biomarkers according to Uno et al. [47]. This method takes the interindividual differences in follow-up and conversion times into account. The t-ROC results could be interpreted in a way similar to that for usual ROC curves.
Survival curves were plotted to visualize the difference in conversion rates between subjects with normal folding and those low and high misfolding (unadjusted curves). For sensitivity analysis (for results, see appendix), similar Cox proportional hazard models were fitted for all available participants. For sensitivity analysis, patients with progression to other dementia were censored at the time of progression.
Analyses were performed with Origin 2017, MATLAB 2015, and R, version 3.5.1 using the packages “coin,, “survival,” “survminer,” and “timeROC.”
Results
An overview of our study population of n = 203 individuals with baseline SCD is listed in Table 1. Significant group differences between the non-converted SCD group and MCI/AD converters were observed for CSF Aβ42 (***p < 0.001, Wilcoxon rank-sum test) (Table 1), plasma Aβ42 (p value 0.003), and Aβ42/40 ratio (p value = 0.002). Furthermore, significant differences were also found for APOEε4, MMSE and age (p value < 0.001, 0.036 and < 0.001).
The structure-biomarker readout, the amide I band maxima, of all 203 participants is shown in Fig. 2 compared to the predefined cutoff values for low and high misfolding. Among the SCD subjects, who did not convert, 51 subjects displayed normal folding patterns (> 1646 cm− 1), the largest group (n = 77) showed low misfolding (≥1644 cm− 1 and ≤ 1646 cm− 1), and 52 subjects exhibited high misfolding (< 1644 cm− 1). In the group that had converted to MCI/AD within 6 years of follow-up, all 22 subjects showed already misfolding at baseline (indicated by red and yellow). A group of 16 of these 22 subjects had highly misfolded Aβ within a mean progression time of 2.2 ± 1.6 years, whereas the other 6 subjects who converted within 3.4 ± 2.2 years showed low Aβ misfolding. One subject that converted to AD was observed in the normal folding group above 1646 cm− 1 (green area). On post hoc inspection of this patient’s report, this individual converted not within the 6 years but after 9.15 years.
Furthermore, the mean progression time from SCD to MCI or AD-dementia was shorter when both the structure of Aβ and the decrease in CSF Aβ42 showed abnormal values. When Aβ misfolding and normal CSF Aβ42 were observed, the average progression time was 3.9 ± 2.8 years, but when Aβ misfolding and lower CSF Aβ42 were recorded, a much shorter progression time of 1.8 ± 1.2 years was observed.
The plasma Aβ structure biomarker correlated to the CSF Aβ42 biomarker, with p = 0.037 and r = 0.15, but not to the plasma Aβ biomarkers (supplementary Fig. S-6). In the converted group, 65% showed decreased CSF Aβ42, whereas in the non-converted SCD group, only 17% showed decreased CSF Aβ42 at baseline. The raw data shown in Fig. 2 largely overlap in the currently short mean follow-up time of 2.7 years. However, different individual follow-up times require an appropriate time-dependent statistical analysis.
To quantify the prognostic value provided by the structure biomarker to predict conversion to clinically diagnosed MCI or AD, we performed Cox proportional hazard analysis (Table 2), which takes into account the variability in follow-up times. If the biomarker frequency in cm− 1 was considered on a continuous scale within the covariate-adjusted Cox proportional hazard model, the risk for MCI or AD conversion increased by the factor of HR = 1.28 per cm− 1 decrease (95% CI 1.1–1.46, p = 0.0002). The HR for APOEε4, age, sex and the structure biomarker were also determined (Table 2). In the models adjusted for age, sex, and APOEε4, the group with high misfolding (< 1644 cm− 1) showed an HR of 19 (95% CI 2.249, 157.826), and the group with low misfolding (≥ 1644 cm− 1 and ≤ 1646 cm− 1) showed an HR of 11 (95% CI 1.048, 110.089) for future clinical conversion compared to the normal folding (> 1646 cm− 1) reference group.
In an extended Cox proportional hazard analysis, we also included those with other forms of dementia (supplementary Table 1). Other forms of dementias were diagnosed in 3.3% of n = 210 subjects. Of these, 1.9% had frontotemporal dementia, 0.5% had vascular dementia, and 1.0% were not further defined. Since misfolding is a specific marker for AD, but does not differentiate between healthy individuals and other forms of dementias, data from those with the other dementias were censored in the model when diagnosed [26]. The hazard ratios were also 11 and 19 when including those with other forms of dementia (supplementary Table 1).
The resulting Kaplan-Meier curves that visualize the probability of conversion within 6 years of follow-up for high misfolding, low misfolding, and normal folding are shown in Fig. 3 for the adjusted (a) and unadjusted (b) models. The three groups identified by the structure biomarker exhibited different levels of risk and differed in both models in contrast to the raw data shown in Fig. 2.
In the next step, we quantified the accuracy of the structure biomarker to predict the conversion to clinically diagnosed MCI or AD-dementia by t-ROC analysis. The converters at each year of follow-up are shown in detail in Table 2 in the supplementary material. In addition to the structure biomarker frequency, the covariates of sex, age, and APOEε4 status were also considered. In Fig. 4a, the ROC curves for 6 years of follow-up are shown. For a better overview of the performance during the follow-up period, the two biomarkers and the biomarker panel are summarized in Fig. 4b. In Fig. 4c, the detailed values are listed. When including the covariates, the AUC for the structure biomarker (in red) was well above 0.84 for all 6 years and increased to 0.94 in the sixth year (Fig. 4). The Aβ42/40 biomarker performance is shown in blue; in green, the performance of the biomarker panel, which includes the structure biomarker, Aβ42/40 and age, sex, and APOEε4, is shown. The structure biomarker and Aβ42/40 perform similarly in the first 3 years, but visually, a better performance of the structure biomarker was observed in the last 3 years. Most notably, with the biomarker panel, the AUC (Fig. 4, in green) increased up to 0.997 for 6 years with a narrower confidence interval than the single biomarker.
Discussion
In the present study, we validated the prognostic performance of Aβ misfolding as a structure biomarker by measuring the Aβ secondary structure distribution in baseline plasma samples from SCD subjects. The prognostic value of this biomarker for future conversion to clinical MCI and dementia due to AD was determined, adjusted for the covariates age, sex, and APOEε4. Additionally, we investigated whether the combination of different plasma biomarkers would increase the prognostic accuracy. For this purpose, we considered not only the misfolding of Aβ but also the plasma Aβ42/40 ratio level.
All 22 converters who developed MCI or dementia due to AD within 6 years of follow-up showed misfolding in the baseline plasma sample and thus could be correctly predicted. One converter who showed no misfolding at baseline converted later, that is, after 9.15 years. Moreover, this subject showed also no abnormal CSF Aβ42 at baseline. In this subject, increased Aβ misfolding in the plasma may have occurred after the baseline collection.
A large number of 52 SCD subjects with high Aβ misfolding had not yet converted and the converted and non-converted SCD subjects showed large overlap. These two groups might not have been differentiated in the raw data, as the average follow-up time in this cohort was only 2.7 years, which might be too short. Previous studies showed that our structure biomarker predicts conversion up to 14 years in advance [35]. In the ESTHER study, the average conversion time to AD-dementia was 8 years after baseline blood sample was taken [25]. In the current study, only those subjects with very rapid disease progression could be identified as converters. We hypothesize that subjects with high misfolding, but not subjects with normal folding (n = 51), will convert in the nearby future. The Amsterdam Dementia Cohort study is still ongoing, and it remains important for the future to repeat the statistical analysis with an extended follow-up period with the non-converted group.
To take this short follow-up time and the individual differences in follow-up and conversion times into account, an appropriate time-dependent statistical analysis was applied. This allowed us to determine the performance of the structure biomarker more accurately than by the raw data alone. Using this statistical analysis approach, the subjects with high misfolding showed a high risk, with an HR of 19, whereas those with low misfolding had a lower risk, with an HR of 11 compared to the normal folding group. Out of 22 converters, 16 showed high misfolding, whereas 6 converters showed low misfolding. The relative decrease in the readout indicated an increasing risk in the HR = 1.28 per cm− 1 decrease. Interestingly, converters with high misfolding of Aβ and an associated higher risk had a faster mean progression time (2.2 ± 1.6 years) than those with low Aβ misfolding (3.4 ± 2.2 years) and lower risk.
Kaplan-Meier survival curves in Fig. 3 demonstrated that the three groups (high misfolding, low misfolding, normal Aβ folding) diverged clearly from each other during the 6 years of follow-up (Fig. 3), which was in contrast to the raw data (Fig. 2). To quantify the accuracy of the prognosis by the structure biomarker, additional time-resolved ROC curves were determined, yielding high AUCs between 0.84 and 0.94 dependent on the length of follow-up time.
Furthermore, we investigated whether the combination of different plasma biomarkers could increase the prognostic precision. For this purpose, we considered not only misfolding of Aβ but also the plasma Aβ42/40 ratios measured by SIMOA technology [15]. The performance of the ratio biomarker was analyzed with the same advanced statistical approach as the structure biomarker and provided similar prognostic values. Most remarkably, the combination of both plasma biomarkers provided an added value with a significantly increased AUC to 0.99 (95% CI 0.99–1.00) for a follow-up period of 6 years. Additionally, it also provided a narrower confidence interval compared to models containing only one of the two markers. The biomarker panel can thus provide a very high prognostic accuracy to predict conversion to MCI or dementia due to AD using plasma biomarkers and covariates only.
The added value provided by the panel might be explained by different biological processes that the biomarkers reflect. Aβ42/40 primarily represents amyloid plaque formation, as this biomarker is highly associated with Aβ-PET scans and amyloidosis [15, 48]. The structure biomarker, on the other hand, monitors the misfolding into β-sheet-enriched structures, which is assumed to occur many years before plaque formation. In this study, Aβ misfolding correlated with the Aβ42 decrease in CSF, but did not correlate with plasma Aβ (see supplementary Fig. S-6). In our former studies, we already reported correlations between the structure-based biomarker in plasma and the Aβ42 decrease in CSF for prodromal to severe AD disease stages [24,25,26]. However, the structure-based biomarker does not correlate with the decrease of Aβ42 in plasma as detected with the iRS in the current study. This might be explained by the early disease state of these patients at baseline. In CSF, the Aβ42 concentration seems to be more sensitive to earlier changes in the brain than the concentration in plasma, which might explain the low correlation with the structure biomarker in plasma. These findings rather support the hypothesis that misfolding of Aβ is one of the initial events in AD development preceding Aβ accumulation in the brain (PET), which reduces the concentration of Aβ42 in CSF and plasma. The structure biomarker would therefore indicate an earlier stage of disease development than the Aβ concentration in plasma. In conclusion, it seems that by combining both plasma biomarkers, a broader time window can be monitored, starting with initial Aβ misfolding and the subsequent amyloid formation. Thus, improvements in prognostic accuracy were achieved when taking both plasma biomarkers into account. Interestingly, the mean conversion time from SCD to MCI or AD-dementia was shorter when both the structure biomarker and the concentration of Aβ42 in CSF showed abnormal values.
In addition to the potential for monitoring different disease-related biological processes, the use of multiple biomarkers might compensate for method-based biases as well. Changes in Aβ concentrations, e.g., due to the circadian cycle, due to an initial rise in liquid Aβ levels in presymptomatic individuals or due to other diseases, could affect the concentration-dependent measurement of SIMOA technology. iRS measurements, in contrast, should be more robust to concentration changes due to biological variation since SIMOA focuses on the ratio between individual conformational species of Aβ rather than on the absolute concentration of selected Aβ isoforms (e.g., Aβ1–42, Aβ1–40). However, misfolding of Aβ may be impaired by cross-seeding of other amyloidogenic proteins, e.g., in the presence of type II diabetes leading to AD-independent conformational changes, which in turn affects the iRS readout but not SIMOA results.
Among the strengths of our study is that we stratified risk groups for SCD patients to convert to clinical MCI or AD within the following 6 years by using only plasma biomarkers. While the previously measured ESTHER cohort is a prospective community-based cohort of older adults, the SCIENCe cohort measured here includes help-seeking individuals in a memory clinic setting. Since cognitive impairment can also be associated with other causes, such as mood states, stress or other neurological disorders, it is important to identify a biomarker that can identify SCD patients in a preclinical AD state. In the ESTHER cohort, we were able to show that the structural biomarker correlates with the future development of AD [25]. For SCD patients, a simple blood screening may be a great benefit. The earlier changes caused by AD can be detected, the more adequately the affected persons can be cared for and treated. In contrast to patients who already have objectively measurable cognitive deficits, SCD patients do not yet exhibit AD-associated irreversible brain damage. Early testing procedures, which can easily be performed on blood as routine and regular testing, can be used to initiate early countermeasures, including nontherapeutic countermeasures such as lifestyle changes, to help patients improve mentally and physically. On the other hand, patients who seem to have no risk can be reassured for the time being.
The Amsterdam Dementia Cohort study is still ongoing, and therefore, it remains important for the future to repeat the statistical analysis with an extended follow-up time with the non-converted group. To confirm the current analysis, it would be important to see that the groups significantly differ not only in the statistically analyzed data but also in the raw data. The broad overlap in the raw data is a weakness of the study, as the average follow-up time is currently too short. Furthermore, an additional external validation study with larger numbers must be performed to validate the findings here.
Conclusion
In summary, this study revealed that a panel of structure- and concentration-based Aβ plasma biomarkers precisely predicts conversion to clinical MCI and dementia due to AD 6 years in advance in individuals with SCD. It provides an earlier time window for screening high-risk symptom-free subjects for potential AD treatments. The proposed plasma biomarker panel including covariates offers a less invasive and cost-effective alternative to currently used CSF biomarkers and PET scanning.
Availability of data and materials
The datasets generated and analyzed during the current study are not publicly available due to local privacy regulations but are available from the corresponding author on reasonable request.
Change history
15 January 2021
An amendment to this paper has been published and can be accessed via the original article.
Abbreviations
- AD:
-
Alzheimer’s disease
- Aβ:
-
Amyloid-β
- AUC:
-
Area under curve
- CSF:
-
Cerebrospinal fluid
- SCD:
-
Subjective cognitive decline
- MCI:
-
Mild cognitive impairment
- APOE:
-
Apolipoprotein E
- SIMOA:
-
Single molecule array
- t-ROC:
-
Time-dependent receiver operating characteristic
- MMSE:
-
Mini-Mental State Examination
- HR:
-
Hazard ratio
- PBS:
-
Phosphate-buffered saline
- IR:
-
Infrared
- CI:
-
Confidence interval
- PET:
-
Positron emission tomography
- OD:
-
Other dementia
References
Bateman RJ, Blennow K, Doody R, Hendrix S, Lovestone S, Salloway S, et al. Plasma biomarkers of AD emerging as essential tools for drug development: an EU/US CTAD task force report. J Prev Alzheimers Dis. 2019;6:169–73.
Witte MM, Foster NL, Fleisher AS, Williams MM, Quaid K, Wasserman M, et al. Clinical use of amyloid-positron emission tomography neuroimaging: Practical and bioethical considerations. Alzheimers Dement Diagn Assess Dis Monit. 2015;1:358–67. Elsevier Inc.
O’Brien JT, Herholz K. Amyloid imaging for dementia in clinical practice. BMC Med. 2015;13:163.
Rowe CC, Bourgeat P, Ellis KA, Brown B, Lim YY, Mulligan R, et al. Predicting Alzheimer disease with β-amyloid imaging: results from the Australian imaging, biomarkers, and lifestyle study of ageing. Ann Neurol. 2013;74:905–13.
Blennow K, Mattsson N, Schöll M, Hansson O, Zetterberg H. Amyloid biomarkers in Alzheimer’s disease. Trends Pharmacol Sci. 2015;36:297–309.
Klunk WE, Engler H, Nordberg A, Wang Y, Blomqvist G, Holt DP, et al. Imaging brain amyloid in Alzheimer’s disease with Pittsburgh Compound-B. Ann Neurol. 2004;55:306–19.
Blennow K, Hampel H, Weiner M, Zetterberg H. Cerebrospinal fluid and plasma biomarkers in Alzheimer disease. Nat Rev Neurol. 2010;6:131–44. Nature Publishing Group.
Janelidze S, Pannee J, Mikulskis A, Chiao P, Zetterberg H, Blennow K, et al. Concordance between different amyloid immunoassays and visual amyloid positron emission tomographic assessment. JAMA Neurol. 2017;74:1492–501.
Pannee J, Portelius E, Minthon L, Gobom J, Andreasson U, Zetterberg H, et al. Reference measurement procedure for CSF amyloid beta (Aβ)1–42 and the CSF Aβ1–42/Aβ1–40 ratio – a cross-validation study against amyloid PET. J Neurochem. 2016;139:651–8.
Janelidze S, Zetterberg H, Mattsson N, Palmqvist S, Vanderstichele H, Lindberg O, et al. CSF Aβ42/Aβ40 and Aβ42/Aβ38 ratios: better diagnostic markers of Alzheimer disease. Ann Clin Transl Neurol. 2016;3:154–65.
Hansson O, Zetterberg H, Buchhave P, Andreasson U, Londos E, Minthon L, et al. Prediction of Alzheimer’s disease using the CSF Aβ42/Aβ40 ratio in patients with mild cognitive impairment. Dement Geriatr Cogn Disord. 2007;23:316–20.
Bateman RJ, Wen G, Morris JC, Holtzman DM. Fluctuations of CSF amyloid-ß levels: implications for a diagnostic and therapeutic biomarker. Neurology. 2007;69:1063–5.
Dubois B, Feldman HH, Jacova C, Hampel H, Molinuevo JL, Blennow K, et al. Advancing research diagnostic criteria for Alzheimer’s disease: the IWG-2 criteria. Lancet Neurol. 2014;13:614–29.
Lövheim H, Elgh F, Johansson A, Zetterberg H, Blennow K, Hallmans G, et al. Plasma concentrations of free amyloid β cannot predict the development of Alzheimer’s disease. Alzheimers Dement. 2017;13:778–82.
Schindler SE, Bollinger JG, Ovod V, Mawuenyega KG, Li Y, Gordon BA, et al. High-precision plasma β-amyloid 42/40 predicts current and future brain amyloidosis. Neurology. 2019;93:e1647–59.
Verberk IMW, Slot RE, Verfaillie SCJ, Heijst H, Prins ND, van Berckel BNM, et al. Plasma amyloid as Prescreener for the earliest Alzheimer pathological changes. Ann Neurol. 2018;84:648–58.
Nakamura A, Kaneko N, Villemagne VL, Kato T, Doecke J, Doré V, et al. High performance plasma amyloid-β biomarkers for Alzheimer’s disease. Nature. 2018;554:249–54.
Verberk IMW, Hendriksen HMA, van Harten AC, Wesselman LMP, Verfaillie SCJ, van den Bosch KA, et al. Plasma amyloid is associated with the rate of cognitive decline in cognitively normal elderly: the SCIENCe project. Neurobiol Aging. 2020;89:99–107. Elsevier Inc.
de Wolf F, Ghanbari M, Licher S, McRae-McKee K, Gras L, Weverling GJ, et al. Plasma tau, neurofilament light chain and amyloid-β levels and risk of dementia; a population-based cohort study. Brain. 2020;143:1220–32.
Albani D, Marizzoni M, Ferrari C, Fusco F, Boeri L, Raimondi I, et al. Plasma Aβ42 as a biomarker of prodromal Alzheimer’s disease progression in patients with amnestic mild cognitive impairment: evidence from the PharmaCog/E-ADNI study. Perry G, editor. J Alzheimers Dis. 2019;69:37–48.
Vergallo A, Mégret L, Lista S, Cavedo E, Zetterberg H, Blennow K, et al. Plasma amyloid β 40/42 ratio predicts cerebral amyloidosis in cognitively normal individuals at risk for Alzheimer’s disease. Alzheimers Dement. 2019;15:764–75.
Ovod V, Ramsey KN, Mawuenyega KG, Bollinger JG, Hicks T, Schneider T, et al. Amyloid β concentrations and stable isotope labeling kinetics of human plasma specific to central nervous system amyloidosis. Alzheimers Dement. 2017;13:841–9. Elsevier Inc.
Nabers A, Ollesch J, Schartner J, Kötting C, Genius J, Haußmann U, et al. An infrared sensor analysing label-free the secondary structure of the Abeta peptide in presence of complex fluids. J Biophotonics. 2016;9:224–34.
Nabers A, Ollesch J, Schartner J, Kötting C, Genius J, Hafermann H, et al. Amyloid-β-secondary structure distribution in cerebrospinal fluid and blood measured by an immuno-infrared-sensor: a biomarker candidate for Alzheimer’s disease. Anal Chem. 2016;88:2755–62.
Nabers A, Perna L, Lange J, Mons U, Schartner J, Güldenhaupt J, et al. Amyloid blood biomarker detects Alzheimer’s disease. EMBO Mol Med. 2018;10:e8763.
Nabers A, Hafermann H, Wiltfang J, Gerwert K. Aβ and tau structure-based biomarkers for a blood- and CSF-based two-step recruitment strategy to identify patients with dementia due to Alzheimer’s disease. Alzheimers Dement Diagn Assess Dis Monit. 2019;11:257–63.
Bateman RJ, Xiong C, Benzinger TLS, Fagan AM, Goate A, Fox NC, et al. Clinical and biomarker changes in dominantly inherited Alzheimer’s disease. N Engl J Med. 2012;367:795–804.
Sarroukh R, Cerf E, Derclaye S, Dufrêne YF, Goormaghtigh E, Ruysschaert J-M, et al. Transformation of amyloid β (1–40) oligomers into fibrils is characterized by a major change in secondary structure. Cell Mol Life Sci. 2011;68:1429–38.
Serpell LC. Alzheimer’s amyloid fibrils: structure and assembly. Biochim Biophys Acta Mol Basis Dis. 2000;1502:16–30.
Milanesi L, Sheynis T, Xue WF, Orlova EV, Hellewell AL, Jelinek R, et al. Direct three-dimensional visualization of membrane disruption by amyloid fibrils. Proc Natl Acad Sci U S A. 2012;109:20455–60.
Zawisza I, Rózga M, Bal W. Affinity of copper and zinc ions to proteins and peptides related to neurodegenerative conditions (Aβ, APP, α-synuclein, PrP). Coord Chem Rev. 2012;256:2297–307. Elsevier B.V.
Benilova I, Karran E, De Strooper B. The toxic Aβ oligomer and Alzheimer’s disease: an emperor in need of clothes. Nat Neurosci. 2012;15:349–57. Nature Publishing Group.
Villemagne VL, Burnham S, Bourgeat P, Brown B, Ellis KA, Salvado O, et al. Amyloid β deposition, neurodegeneration, and cognitive decline in sporadic Alzheimer’s disease: a prospective cohort study. Lancet Neurol. 2013;12:357–67. Elsevier Ltd.
Schartner J, Nabers A, Budde B, Lange J, Hoeck N, Wiltfang J, et al. An ATR–FTIR sensor unraveling the drug intervention of Methylene Blue, Congo Red, and Berberine on Human Tau and Aβ. ACS Med Chem Lett. 2017;8:710–4.
Stocker H, Nabers A, Perna L, Möllers T, Rujescu D, Hartmann A, et al. Prediction of Alzheimer’s disease diagnosis within 14 years through Aβ misfolding in blood plasma compared to APOE4 status, and other risk factors. Alzheimers Dement. 2020;16:283–91.
van der Flier WM, Scheltens P. Amsterdam dementia cohort: performing research to optimize care. J Alzheimers Dis. 2018;62:1091–111. Perry G, Avila J, Tabaton M, Zhu X, editors.
Slot RER, Verfaillie SCJ, Overbeek JM, Timmers T, Wesselman LMP, Teunissen CE, et al. Subjective Cognitive Impairment Cohort (SCIENCe): study design and first results. Alzheimers Res Ther. 2018;10:76.
van der Flier WM, Pijnenburg YAL, Prins N, Lemstra AW, Bouwman FH, Teunissen CE, et al. Optimizing patient care and research: the Amsterdam dementia cohort. J Alzheimers Dis. 2014;41:313–27.
Stewart R. Subjective cognitive impairment. Curr Opin Psychiatry. 2012;25:445–50.
Tijms BM, Willemse EAJ, Zwan MD, Mulder SD, Visser PJ, van Berckel BNM, et al. Unbiased approach to counteract upward drift in cerebrospinal fluid amyloid-β 1–42 analysis results. Clin Chem. 2018;64:576–85.
Jessen F, Amariglio RE, van Boxtel M, Breteler M, Ceccaldi M, Chételat G, et al. A conceptual framework for research on subjective cognitive decline in preclinical Alzheimer’s disease. Alzheimers Dement. 2014;10:844–52.
Petersen RC, Smith GE, Waring SC, Ivnik RJ, Tangalos EG, Kokmen E. Mild cognitive impairment. Arch Neurol. 1999;56:303–9.
Albert MS, DeKosky ST, Dickson D, Dubois B, Feldman HH, Fox NC, et al. The diagnosis of mild cognitive impairment due to Alzheimer’s disease: recommendations from the National Institute on Aging-Alzheimer’s Association workgroups on diagnostic guidelines for Alzheimer’s disease. Alzheimers Dement. 2011;7:270–9.
McKhann G, Drachman D, Folstein M, Katzman R, Price D, Stadlan EM. Clinical diagnosis of Alzheimer ’ s disease: report of the NINCDS-ADRDA work group under the auspices of department of health and human services task force on Alzheimer’s disease. Neurology. 1984;34:939.
McKhann GM, Knopman DS, Chertkow H, Hyman BT, Jack CR, Kawas CH, et al. The diagnosis of dementia due to Alzheimer’s disease: recommendations from the National Institute on Aging-Alzheimer’s association workgroups on diagnostic guidelines for Alzheimer’s disease. Alzheimers Dement. 2011;7:263–9.
Budde B, Schartner J, Tönges L, Kötting C, Nabers A, Gerwert K. Reversible Immuno-infrared sensor for the detection of Alzheimer’s disease related biomarkers. ACS Sensors. 2019;4:1851–6.
Uno H, Cai T, Tian L, Wei LJ. Evaluating prediction rules for t-year survivors with censored regression models. J Am Stat Assoc. 2007;102:527–37.
Jack CR, Bennett DA, Blennow K, Carrillo MC, Dunn B, Haeberlein SB, et al. NIA-AA research framework: toward a biological definition of Alzheimer’s disease. Alzheimers Dement. 2018;14:535–62.
Acknowledgements
The authors kindly thank David Holtzman and the members of Holtzman Lab for providing HJ5.1 hybridoma cells.
Funding
The presented research was funded by the Protein Research Unit within Europe (PURE), Ministry of Innovation, Science and Research of North-Rhine Westphalia, Germany, and the Center for Protein Diagnostics (PRODI), Ministry of Culture and Science of North-Rhine Westphalia, Germany. Patents regarding the iRS and the corresponding methods have been applied for (WO2015121339, EP16199792, EP16199805, EP17173322.3). Furthermore, funding was provided by Alzheimer Forschung Initiative e.V. (AFI, #19077CB).
Inge Verberk is supported by research grants from Gieskes-Strijbis Fonds and Alzheimer Nederland (NL-17004). Wiesje van der Flier holds the Pasman chair. The VUmc biobank is supported by VUmc. The Alzheimer Center Amsterdam is supported by Alzheimer Nederland and Stichting VUmc fonds, and the clinical database structure was developed with funding from Stichting Dioraphte. The SCIENCe prospective cohort study is supported by Gieske-Strijbis Fonds. Research of the Alzheimer Center Amsterdam is part of the neurodegeneration research program of Amsterdam Neuroscience. Open Access funding enabled and organized by Projekt DEAL.
Author information
Authors and Affiliations
Contributions
JS, IV, PS, WvdF, AN, CT, and KG conceived, designed, and supervised the study. JS, JL, and BB conducted the experiments, and JS, IV, PS, WvdF, AN, CT, and KG acquired the data. BB established the quality control procedures for iRS. JS established the procedure for capture antibody production. All authors analyzed the data and helped in preparing the figures. JL, IV, RD, NT, and AN performed the statistical analyses. All authors approved and critically evaluated the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Written informed consent to use medical data and biomaterials for research purposes was in place, in accordance with ethical approval from the VU University Medical Center based on the Helsinki Declaration (seventh revision).
Consent for publication
Not applicable.
Competing interests
The authors declare no conflicts of interest.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original online version of this article was revised: the authors noticed that the published figures have errors which was occurred during processing of the figures in production team
Supplementary Information
Additional file 1: Supplementary Table 1.
Cox proportional hazard regression models including other forms of dementia.
Additional file 2: Supplementary Table 2.
Risks and events by year with amide I maximum ≥1644 cm− 1, < 1644 cm− 1, and all subjects (n = 203).
Additional file 3: Supplementary Figure S-1.
HJ 5.1 binding characteristics. The antibody recognizes Aβ monomers derived from chicken telencephalonic cells with an amide I maximum of 1653 cm− 1 (A) and Aβ1–42 fibrils with a maximum of 1628 cm− 1 (B). CSF from healthy individuals yielded a maximum of 1645 cm− 1 (C) above the threshold and for AD subjects a maximum of 1640 cm− 1 below the threshold. (D).
Additional file 4: Supplementary Figure S-2.
Reproducibility of plasma measurements and specific binding of Aβ with HJ 5.1. (A) Repeated analysis of the same antibody-depleted plasma sample with HJ 5.1 results in the same read out at 1642 cm-1. (B) A control measurement of an Aβ-depleted sample shows no signal at all, indicating no unspecific binding of plasma proteins.
Additional file 5: Supplementary Figure S-3.
Depletion of plasma samples provide the same readout. (A) The same plasma sample shows the same spectra for both capture antibodies, HJ 5.1 (depleted sample, black) and A8978 (non-depleted sample, dotted line) indicating that both antibodies extract the same Aβ fractions out of the plasma sample and providing the same read out at 1640 cm-1. (B) Single depletion (dotted line) and double depletion (grey) of a sample containing few amounts of human antibodies results in the same read-out when measured with HJ 5.1 as with A8978 (black). (C) A different non-depleted plasma sample with a higher amount of antibodies is dominated by antibodies absorbing at 1638 cm-1 and a shoulder at 1685 cm-1 (grey). For comparison an antibody spectrum is shown in grey with a dotted line. Single depletion of a plasma sample containing a large amount of human antibodies still contains antibody contributions (blue). However, the double-depleted sample shows no longer antibody contribution (black, dotted line). The read-out is again the same as with the A8978 (black). All depleted samples underwent a quality control to ensure complete removal of human IgGs.
Additional file 6: Supplementary Figure S-4.
No correlation could be observed between amide I maximum and age of the samples (rs= 0.01, p-value = 0.87).
Additional file 7: Supplementary Figure S-5.
Binding of Aβ from plasma and CSF samples with HJ 5.1 as the capture antibody. (A) Amide I signal from a plasma sample of an SCD patient who did not convert (grey) and an SCD patient who did convert during the follow-up time (black). Amide I maxima varied from 1647-1651 cm-1 during circulation due to unbound plasma proteins that interfere with the signal. During washing, only the Aβ signal is left, showing constant and differentiable amide I maxima. (B) Amide I signal from a CSF healthy control sample (grey) and an AD subject (black). The signal stays stable during circulation and washing, indicating no change in the conformational structure of bound Aβ. Less interfering proteins in CSF samples did not influence the signal.
Additional file 8: Supplementary Figure S-6.
Correlation plots between CSF and plasma levels of Aβ and amide I maximum frequency. The appropriate biomarker is presented on the y-axis, the amide I maximum position on the x-axis and the spearman rank value (r) and p-value above the distribution points, equally. SCD subjects are plotted in blue, MCI subjects in orange and AD subjects in red. Biomarker cutoffs are shown as dashed lines. The structure-based biomarker showed weak correlation with CSF Aβ42 but not with the concentration biomarkers in plasma (Aβ40, Aβ42 and Aβ40/42).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Stockmann, J., Verberk, I.M.W., Timmesfeld, N. et al. Amyloid-β misfolding as a plasma biomarker indicates risk for future clinical Alzheimer’s disease in individuals with subjective cognitive decline. Alz Res Therapy 12, 169 (2020). https://doi.org/10.1186/s13195-020-00738-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13195-020-00738-8