Skip to main content

Inter-network connectivity and amyloid-beta linked to cognitive decline in preclinical Alzheimer’s disease: a longitudinal cohort study

Abstract

Background

Amyloid-beta (Aβ) has a dose-response relationship with cognition in healthy adults. Additionally, the levels of functional connectivity within and between brain networks have been associated with cognitive performance in healthy adults. Aiming to explore potential synergistic effects, we investigated the relationship of inter-network functional connectivity, Aβ burden, and memory decline among healthy individuals and individuals with preclinical, prodromal, or clinical Alzheimer’s disease.

Methods

In this longitudinal cohort study (ADNI2), participants (55–88 years) were followed for a maximum of 5 years. We included cognitively healthy participants and patients with mild cognitive impairment (with or without elevated Aβ) or Alzheimer’s disease. Associations between memory decline, Aβ burden, and connectivity between networks across the groups were investigated using linear and curvilinear mixed-effects models.

Results

We found a synergistic relationships between inter-network functional connectivity and Aβ burden on memory decline. Dose-response relationships between Aβ and memory decline varied as a function of directionality of inter-network connectivity across groups. When inter-network correlations were negative, the curvilinear mixed-effects models revealed that higher Aβ burden was associated with greater memory decline in cognitively normal participants, but when inter-network correlations were positive, there was no association between the magnitude of Aβ burden and memory decline. Opposite patterns were observed in patients with mild cognitive impairment. Combining negative inter-network correlations with Aβ burden can reduce the required sample size by 88% for clinical trials aiming to slow down memory decline.

Conclusions

The direction of inter-network connectivity provides additional information about Aβ burden on the rate of expected memory decline, especially in the preclinical phase. These results may be valuable for optimizing patient selection and decreasing study times to assess efficacy in clinical trials.

Background

Elevated levels of amyloid-beta (Aβ), a neuropathological hallmark of Alzheimer’s disease (AD) [1,2,3,4,5,6], are crucial in identifying the earliest stages of AD. While Aβ does not relate well to cognition cross-sectionally, it has a dose-response relationship with cognition in healthy adults [7]. However, dementia-related pathologies, including Aβ, only explain 41% of variation in cognitive decline [8] and a low signal-to-noise ratio of biomarkers in the asymptomatic stages of the disease for entry criteria in clinical trials has been reported as a reason for trial failure [9].

Healthy aging has previously been associated with reduced activity within the default mode network (DMN), a network often associated with memory-related processes such as thinking about the future, episodic memory, and autobiographical memory [10]. Reduced connectivity within the DMN, as measured with resting-state functional magnetic resonance imaging, is also associated with impaired cognitive performance in old age [11]. Additionally, levels of functional connectivity within cognition-related intrinsic brain networks, such as the DMN, and Aβ burden have a synergistic effect on memory decline in clinically normal older individuals [12]. Various networks predict cognitive decline which may indicate that Aβ burden impacts the interaction between networks. Negative correlations between the DMN and task-positive networks has been positively associated with cognitive performance in young individuals [13, 14]. These negative correlation patterns are also referred to as “anti-correlations” in the literature and have been described as intrinsically organized antagonistic activation patterns between networks in the brain [15]. These patterns are reported to be part of a mechanism that facilitates cognition, possibly by reinforcing connections between two loci in the brain dedicated to cognitive functions [14, 16].

Previous studies have reported attenuated effects of aging with negative correlations between resting-state networks, with negative correlations further decreasing in patients with mild cognitive impairment (MCI) and AD [17,18,19]. While detrimental effects of Aβ on cognition and on functional connectivity within networks have been shown [2, 20,21,22], it remains unknown whether Aβ modulates functional connectivity between networks. Given the close relationship between inter-network connectivity and cognitive decline and these first reported associations with Aβ, combining information from inter-network connectivity with Aβ may reduce noise when selecting asymptomatic individuals at risk for AD in clinical trials aimed at slowing down cognitive decline.

To that end, we investigated whether functional connectivity between the DMN and task-positive networks predicts memory decline differently among cognitively normal individuals, MCI patients with and without elevated Aβ, or AD patients. Additionally, we investigated whether the relationship between inter-network connectivity and memory decline depends on Aβ levels in a dose-response type relationship. We expected that, as the magnitude of negative correlations decreases, memory performance would be lower, and that the strength of these associations would show a dose-response relationship with Aβ burden. To investigate domain-specificity of our findings, we have also investigated our hypotheses using executive functions as a control outcome measure.

Methods

Data used in this article were obtained from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) database (adni.loni.usc.edu). The ADNI was launched in 2003, led by Principal Investigator Michael W. Weiner, MD. The main goal of the ADNI is to test whether magnetic resonance imaging, positron emission tomography, other biological markers, and clinical and neuropsychological assessment can be combined to measure the progression of MCI and AD. For up-to-date information, see www.adni-info.org.

Participants

We included data from a total of 122 eligible participants from the ADNI2 study (as of March 2017), of which seven were removed due to low imaging data quality, resulting in a total of 115 participants. Qualified clinicians working for the ADNI categorized participants into four groups (cognitively normal, early MCI, late MCI, or AD) based on diagnostic procedures from the ADNI protocol [23]. However, for the purposes of our study, we grouped all patients with MCI based on their Aβ levels being below or above a pathological cutoff point of 1.11 18F-AV-45 florbetapir positron emission tomography standardized uptake value ratio [24], respectively (referred to as MCI and MCI+ in this manuscript). MCI diagnoses were based on the Petersen criteria for MCI [25, 26], while patients with AD met the NINCDS-ADRDA criteria for probable Alzheimer’s disease [27]. For a detailed overview of exclusion, inclusion, and diagnostic criteria and procedures, please refer to the ADNI2 protocol (http://adni.loni.usc.edu/wp-content/uploads/2008/07/adni2-procedures-manual.pdf). Additionally, only participants with complete resting-state functional magnetic resonance imaging and Aβ data at baseline were eligible for inclusion, and patients with AD were only eligible if their Aβ levels were above the cutoff point to ensure only AD-related pathology was represented in this group. A complete listing of all participant IDs that were included in the final analyses can be found in Additional file 1 (Table S1).

Materials and equipment

Test batteries

The main outcome measures we used were the composite memory score (ADNI-Mem) and the composite executive functions score (ADNI-EF), which have good validity and are ideally suited to track changes over time. These scores were derived by combining scores related to memory performance and executive functions. ADNI-Mem score is derived from several test batteries, including the Rey Auditory Verbal Learning Test, Alzheimer’s Disease Assessment Schedule-Cognition, Mini-Mental State Examination, and Wechsler Memory Scale-Revised [28]. ADNI-EF score is derived from tests including WAIS-R Digit Symbol Substitution, Digit Span Backwards, Trails A and B, Category Fluency, and Clock Drawing [29].

Biomarker assessment

18F-AV-45 florbetapir positron emission tomography measures were used to quantify levels of neocortical Aβ at baseline. The duration of positron emission tomography imaging was 20 min and started 50 min after injection of tracer fluid. The neocortical standardized uptake value ratio is the mean uptake in an aggregate of the frontal lobe, cingulate cortex, lateral parietal, and lateral temporal regions relative to mean uptake in the whole cerebellum, including white and gray matter. Further processing of positron emission tomography images was performed as described in a previous report [24]. Participants were characterized as Aβ-positive if they exceeded the cutoff value of 1.11 standardized uptake value ratio, as previously determined in ADNI cohorts [24].

Imaging equipment and acquisition

Imaging data were acquired using Philips Medical Systems 3.0-Tesla magnetic resonance systems. Structural T1-weighted gradient echo pulse sequence data with dimensions of 170 × 256 × 256 mm with a voxel resolution of 1.2 × 1 × 1 mm were acquired in sagittal orientation with a repetition time of 6.8 ms, echo time of 3.1 ms, flip angle of 9°, and slice thickness of 1.2 mm. Additionally, resting-state functional magnetic resonance imaging scans of 7 min were obtained, consisting of 140 volumes of T2*-weighted data with 48 slices per volume, dimensions of 64 × 64 × 48 mm, and a voxel resolution of 3.3 × 3.3 × 3.3 mm. Functional images were acquired in transverse orientation with repetition time of 3000 ms, echo time of 30 ms, flip angle of 80°, and a slice thickness of 3.3 mm.

Image preprocessing and de-noising

For preprocessing we used the default preprocessing pipeline for volume-based analyses within the NITRC CONN toolbox (version 17.a; http://www.nitrc.org/projects/conn/) [30]. Images were realigned and unwarped, centered, slice-time corrected, segmented (gray/white/cerebrospinal fluid), normalized to MNI space, and outliers were identified with the Artifact Detection Toolbox. Outliers were regressed out (scrubbing) using conservative settings (95th percentile in normative sample, global signal Z value = 2, motion = 0.5 mm, and we discarded the first three volumes). Functional data were smoothed with a Gaussian kernel of 6 mm full width at half maximum. As the reliability and validity of negative correlation measures has been the subject of debate [31], de-noising of data involved regressing out principal components of the signal from white matter and cerebrospinal fluid following the CompCor method [32] which reduces the presence of artificial negative correlations. Additionally, linear detrending and an after-regression BOLD signal band-pass filter (0.008 < f <  0.09 Hz) were applied. Data quality was ensured by visual inspection of histograms of functional connectivity values before and after de-noising for each participant. Histograms that did not show a normal distribution of functional connectivity values after de-noising were indicative of a suboptimal de-noising process and these participants were excluded from data analysis to preserve a high level of data quality (n = 7). Motion parameters were entered as regressor in our analyses. Volumes with motion above 0.5 mm were entered as regressors in our general linear models (“scrubbing”). Furthermore, all images were inspected for irregularities during the entire preprocessing and analysis process.

Functional connectivity analyses

A bivariate correlation, hemodynamic response function-weighted ROI-to-ROI analysis was performed within CONN, using four independent network regions of interest (ROI) as sources, with a total of 19 structurally defined hub regions. The regions of interest include the DMN and regions representing task-positive networks, including the dorsal attention network (DAN), salience network (SN), and frontoparietal network (FPN). These regions of interest, which are part of the CONN software package, were generated by an independent component analysis on 497 healthy control participants (293 females) as part of the Human Connectome Project (http://www.humanconnectome.org) [30]. As a reference for future replication efforts, peak coordinates of the regions of interest used for this functional connectivity analysis are provided in Additional file 1 (Table S2) and the networks used are visualized in Fig. 1. To assess inter-network connectivity, the average time course signal from each network was extracted and Fisher r-to-z transformed inter-network correlation values were exported for statistical analyses.

Fig. 1
figure 1

Linear effects of functional connectivity on memory decline over time per group. The effect of baseline functional connectivity between networks on memory performance over time was investigated in three network pairs in all groups. Functional connectivity values were standardized to ensure that the mean reflects a z-value of 0. Red, green, and blue lines indicate the estimated marginal means for the moderation by negative (–1 standard deviation (SD)), positive (+1 SD), and no (mean) correlation between networks, respectively, but the analyses were performed using functional connectivity measures continuously. Network combinations from top to bottom: DMN-DAN, DMN-SN, and DMN-FPN. The DMN is shown in dark blue, the DAN in green, the SN in yellow, and the FPN in light blue. The brain images give a superior viewpoint of the brain; top = anterior, L = left, R = right, bottom = posterior. The AD group only had a maximum follow-up time of 2 years (Table 1). All p values are corrected for multiple comparisons using FDR. Significant effects are indicated by an asterisk and bold font. a DMN-DAN: significant effects were found in the cognitively normal and MCI+ groups. b DMN-SN: significant effects were found in the cognitively normal and MCI+ groups. c DMN-FPN: no significant effects were found for this network pair

Statistical analysis

All data were analyzed using R version 3.3.1 (https://www.r-project.org/) and MATLAB R2016a. Baseline participant characteristics were compared with analysis of variance for continuous variables and chi-squared test for categorical variables. Linear regression analyses were used for cross-sectional data. Because of their ability to deal with missing data in longitudinal studies without being subject to complete-case bias [33, 34], longitudinal analyses were performed with linear and curvilinear mixed-effect regression analyses using maximum likelihood, utilizing nlme version 3.1–128 [35]. Fixed effects were predictor of interest, a random intercept for each participant, and random slope for time. Time was calculated as the number of years since baseline assessment. During the construction of statistical models, we first investigated effects in the whole sample. In case of a significant effect, the sample was split based on clinical groups and the level of amyloid burden, for which interactive analyses were performed with the cognitively normal group as reference. We also performed within-group analyses to better understand the patterns. For all models, we compared the Akaike Information Criteria between models with a random intercept and random slope or a random intercept only using the Log-likelihood ratio test and selected the most parsimonious model. In all models, age, sex, education, and their interaction with time were included as covariates if p <  0.10 (using the Wald t statistic) [36]. To improve longitudinal data quality [37], variability in head size and intracranial volume, as measured with FreeSurfer version 5.1 [38], was added as a covariate [39]. Additionally, intracranial volume can be used as a measure of head size, which correlates with total gray matter atrophy. It is important to control for atrophy as atrophy can reduce blood flow [40], and AD patients may show more atrophy than patients with MCI or healthy controls. The most complex model constructed for analyses is described in detail in Additional file 1 (Box 1). Residual plots and Q-Q plots were examined for all models. Significance was set at p <  0.050 (two-sided) and results were corrected for multiple comparisons using the false-discovery rate (FDR). Power calculations were performed using MATLAB scripts for the mixed-effects model using the slope of memory decline and residual variance [40] to estimate the number of participants a clinical trial would need to enroll to detect slowing of memory decline of 30% (two arms for 4-year annual assessments, 80% power, alpha = 0.05) using high amyloid levels as an inclusion criterion. These mixed-effects models included the covariates, random intercept, and random slope.

Results

Demographics

The total sample consisted of 115 participants, including 53 females (46%) with a mean age of 72.91 (standard deviation (SD) = 6.68) at baseline. Of these participants, 28 were in the cognitively normal group, 27 were in the MCI group, 36 were in the MCI+ group, and 24 were in the AD group. Other participant characteristics are summarized in Table 1.

Table 1 Summary table of participant characteristics at baseline

Analyses of variance indicated group differences in 18F-AV-45 florbetapir standardized uptake value ratio, ADNI-Mem, ADNI-EF, clinical dementia rating scale scores, and follow-up times. The proportion of cases at each time point were as follows: baseline, 115 (100%); month 6, 110 (96%); month 12, 98 (85%); month 24, 74 (64%); month 36, 6 (5%); month 48, 29 (25%); and month 60, 1 (< 1%). Results of Tukey tests for post-hoc differences between group characteristics are available in Additional file 1 (Table S3). Additionally, we found no differences in image preprocessing parameters between groups (Additional file 1: Table S4).

Effects of functional connectivity and group on memory at baseline

There were no significant associations between functional connectivity between networks and memory performance at baseline across the entire sample or within the groups. An association was found between functional connectivity between the DMN and FPN and executive functions across the entire sample at baseline. Further analyses of the association between functional connectivity and executive functions within the diagnosis groups reveal an association in the amyloid-negative MCI group (Additional file 1: Table S5).

Longitudinal linear effects of functional connectivity and group on memory

Linear mixed-effect models were performed to investigate effects of functional connectivity and group on memory decline over time. The associations between functional connectivity between the networks and time for the whole sample were not significant (Additional file 1: Table S6). Using the cognitively normal group as a reference level, a significant three-way interaction effect of time, group, and functional connectivity was found in comparison with the MCI+ group for the DMN-DAN and DMN-SN correlations, but not for the DMN-FPN correlations (Table 2). Adding inter-network functional connectivity to the models with the covariates contributed significantly to the explained variance of memory decline for the DMN-DAN (R2 difference = 0.03, 95% confidence interval (CI) 0.01 to 0.07; p = 0.04), but not for the DMN-SN (R2 difference = 0.001, 95% CI −0.001 to 0.01; p = 0.21) or the DMN-FPN (R2 difference < 0.001, 95% CI −0.001 to 0.01; p = 0.48).

Table 2 Associations between group, inter-network functional connectivity, and memory decline

Figure 1 and Table 2 show significant associations between DMN-DAN and DMN-SN functional connectivity and memory decline between the cognitively normal and MCI+ groups. Post-hoc linear mixed models within each group confirmed that, for the cognitively normal group, positive correlations between these network pairs are positively associated with memory performance over time, whereas in the MCI+ group positive correlations were negatively associated with memory performance (Additional file 1: Table S7). No significant associations were found for executive functions (Additional file 1: Table S8).

As the cognitively normal group consisted of 11 individuals with Aβ levels above the threshold, we also examined associations between functional inter-network connectivity, Aβ as a continuous variable, and memory decline in the cognitively normal group. These results showed that this interaction was significant for the DMN-DAN and the DMN-SN, but not for the DMN-FPN. This shows that Aβ moderates the association between functional connectivity and memory decline also in the cognitively normal group in such a way that negative correlations were associated with memory decline with elevated Aβ (Table 3). Similar findings are observed when grouping cognitively normal individuals based on the Aβ cutoff value (Additional file 1: Table S9).

Table 3 Three-way interaction of functional connectivity, amyloid-beta, and time in the cognitively normal group

Longitudinal curvilinear effects of functional connectivity and amyloid-beta on memory

Since our linear mixed model results showed opposite moderations by Aβ on the association between functional connectivity and memory performance across the diagnostic groups, we examined curvilinear (quadratic) mixed-effect models for these network interactions using Aβ as a continuous variable to investigate a possible dose-response relationship with memory decline. To ensure that individuals with the highest Aβ levels would not drive associations, we excluded the AD group from these analyses. In the whole sample, excluding the AD group, a quadratic three-way interaction effect of functional connectivity, Aβ, and time was found for the DMN-DAN, but marginally not for the DMN-SN correlations (Table 4). To visualize this interaction, we plotted simple slopes for different values of 18F-AV-45 florbetapir standardized uptake value ratio in both groups, broken down by positive and negative inter-network correlations and holding all other fixed effects constant (Fig. 2). Values for Aβ were chosen to reflect Aβ negativity (0.9), the cutoff (1.1), slightly elevated Aβ (1.3), moderately elevated Aβ (1.5), and high Aβ (1.7).

Table 4 Curvilinear associations between functional inter-network connectivity, amyloid-beta burden, and memory decline
Fig. 2
figure 2

Curvilinear three-way interaction of functional connectivity, amyloid-beta and time and its relationship with memory performance. To visualize the three-way interaction, the whole sample, excluding the AD group, was split based on whether functional connectivity between DMN-DAN was below or above zero, allowing us to visualize the dose-response relationship of Aβ burden (measured with 18F-AV-45 florbetapir standardized uptake value ratio) on memory decline in both subgroups. 18F-AV-45 florbetapir standardized uptake value ratio values over 1.11 indicate Aβ positivity. Linear or curvilinear graphs were drawn based on which type of association showed the best fit for each model. All p values are FDR-corrected. Significant effects are indicated by an asterisk and bold font. a Effects within the subgroup of cognitively normal participants and MCI patients showing positive correlations between networks. b Curvilinear effects within the subgroup of cognitively normal participants and MCI patients showing negative correlations between networks. c Effects within the cognitively normal subgroup with positive correlations between networks; no significant dose-response relationship of Aβ on memory (practice effects). d Effects within the MCI subgroup with positive correlations between networks showing a dose-response relationship of Aβ on memory. e Effects within the cognitively normal subgroup with negative correlations between networks showing a curvilinear dose-response relationship of Aβ on memory. f Effects within the MCI subgroup with negative correlations between networks showing no (borderline) significant effect of Aβ burden on memory

Figure 2 shows that the dose-response relationship between Aβ and memory decline is modulated by both diagnosis and the direction of DMN-DAN correlations. Within the cognitively normal group, negative DMN-DAN correlations have a curvilinear dose-response relationship with memory decline. However, when networks are positively correlated, we found no significant association between Aβ and memory decline. Conversely, in the MCI group, positive DMN-DAN correlations were associated with a linear dose-response relationship between Aβ and rate of memory. No moderation effects of Aβ were found for the MCI group with negative correlations, although linear effects approached significance (Table 4).

Power estimations were performed in the cognitively normal group to investigate the added value of adding functional inter-network connectivity to higher levels of Aβ burden as an inclusion criterion on the required sample size per arm in clinical trials (Fig. 3).

Fig. 3
figure 3

Power analysis in the cognitively normal group. To assess the effect of including inter-network connectivity as an inclusion criterion in clinical trials, we performed a power analysis in the cognitively normal group. The x axis describes the memory slope reduction in percentages and the y axis describes the number of participants needed per arm to detect this reduction with 80% power and α = 0.05 in a 4-year trial with annual assessments. The green line shows the sample size needed when only amyloid is used as an inclusion criterion and inter-network connectivity between the DMN and DAN is not considered. The orange line shows the required sample size when positive inter-network correlations between the DMN and DAN and amyloid are used as inclusion criteria. The blue line shows the required sample size when negative inter-network correlations between the DMN and DAN and amyloid are used as inclusion criteria

These results show that selecting cognitively normal participants with higher levels of Aβ burden and who also have a negative DMN-DAN correlation may greatly reduce the sample size needed to detect amyloid-related changes in memory decline in clinical trials as compared with only including preclinical AD individuals or those with positive correlations. For instance, the required sample size per arm to detect an effect size of 0.3 (30%) can be reduced by 88% when using both amyloid and negative inter-network correlations (n = 37) as inclusion criteria in clinical trials towards asymptomatic AD, as opposed to considering only amyloid (n = 318) in the inclusion process.

While functional inter-network connectivity is a relatively affordable and easy measure to collect in participants, this additional requirement can increase the rate of screen failures at inclusion.

Discussion

This study provides evidence for a dose-response relationship between Aβ burden and inter-network connectivity on memory decline over 4 years in both healthy and patient populations. Previous studies suggested that negative correlations between task-positive and task-negative networks might have a beneficial effect on memory in young individuals [13, 14]. We now show that in older adults these associations depend on disease stage and the burden of Aβ. In cognitively normal older individuals, negative correlations were associated with memory decline, and this association became stronger in individuals with higher levels of Aβ. Conversely, in patients with MCI, the magnitude of Aβ burden predicted the rate of memory decline associated with positive inter-network correlations.

These results suggest that the direction of baseline inter-network connectivity can provide additional information when combined with baseline Aβ burden about the rate of expected memory decline and can add important information to the complexity of factors contributing to cognitive decline. Additionally, these findings can have prognostic implications when using functional connectivity as a potential biomarker and show the importance of including both the direction of inter-network connectivity and Aβ burden for projecting cognitive trajectories as outcome, especially in preclinical AD.

Interestingly, the associations we observed involved mainly the DMN-DAN correlations. The DMN is one of the most investigated networks and is among the first networks where Aβ accumulates significantly early in the disease [41]. It is therefore likely that other functional networks, especially networks tightly coupled to the DMN, are indirectly susceptible to Aβ-related alterations. Animal and electrophysiological studies have provided evidence that Aβ burden leads to synaptic dysfunction and network disorganization [42].

Widespread inter-network reorganization under the influence of Aβ have indeed been reported for the coupling between DMN and DAN [43], DMN and SN [19], and DMN and FPN [43], suggesting that topographical closeness to certain hubs, such as the association cortices of the DMN, may put networks at increased risk for Aβ toxicity. Our findings now show that functional implications of the interactions between large-scale brain networks and amyloid in relation to memory decline were mainly observed for the DMN-DAN correlations, with borderline significant associations for the DMN-SN correlations. This is in line with previous research suggesting that an optimal interaction between these specific networks is important for cognitive functioning [13, 14, 44]. The fact that we observed associations between greater positive inter-network correlations and better memory performance over time in clinically normal individuals suggests the presence of compensation mechanisms involving possible reconfiguration of large-scale networks [45, 46]. In the context of preclinical AD, this compensation may provide the architectural basis for maintaining optimal memory performance in old age [45, 47]. These compensation mechanisms may be moderated by cognitive reserve [46]. Individuals with preclinical AD may be able to compensate for AD pathology in the form of positive inter-network correlations. Furthermore, previous studies reported differential connectivity patterns for the anterior versus posterior DMN [48]. Increased within-network functional connectivity in the anterior DMN was associated with higher levels of amyloid pathology in preclinical AD [21], indicative of regional compensatory mechanisms. However, as the disease progresses, compensatory mechanisms may be dependent on cognitive reserve to stave off memory decline. Attenuated negative correlations between the DMN-DAN may signal impending neuronal breakdown resulting from Aβ toxicity in the prodromal phase, a finding consistent with previous work [46].

The DMN and DAN may operate in a negative correlation pattern to facilitate attention-related processes during cognitive tasks. The DAN is thought to coordinate cross-talk between brain networks, especially during memory tasks. The coordination between brain networks may become especially important in MCI, where medial temporal lobe structures become functionally isolated from other regions or networks. In MCI, a coordinated antiphase functional connection between the DMN and DAN may be fundamental in maintaining optimal levels of memory performance [46].

Finally, functional connectivity deteriorates as the disease progresses to more advanced stages [48]. This may partially explain the lack of connectivity-related findings in the AD group, as coherent connectivity patterns between networks may dissipate leaving no association between inter-network connectivity and cognition.

Our dose-response findings relating Aβ burden to memory decline, depending on functional connectivity status, showed different patterns between cognitively normal individuals and patients with MCI. This has important implications for early detection of individuals at risk of cognitive decline and the selection of preclinical and prodromal AD individuals for trials using a combination of functional connectivity and Aβ. Patients with MCI showed a dose-response relationship between Aβ and memory decline for positive DMN-DAN correlations. In clinically normal individuals, Aβ has a dose-response relationship with memory decline when DMN-DAN correlations are negative. Selecting preclinical AD individuals with positive DMN-DAN correlations could lead to the selection of individuals where individuals with higher levels of Aβ have similar practice effects as individuals with lower levels of Aβ burden. Thus, the direction of inter-network connectivity in combination with Aβ burden can have important implications for participant selection in clinical trials. This notion is further outlined by the results from our power analysis showing that the selection of cognitively normal individuals with elevated Aβ levels and negative inter-network correlations reduced the required sample size by 88% to slow down memory decline by 30%.

Furthermore, we did not observe independent or synergistic effects of inter-network connectivity or Aβ on cognition in patients with AD in the 2-year follow-up period. This could be related to widespread deposits of both Aβ and tau that may have impacted the integrity of the entire brain. We also did not observe between-group differences in inter-network connectivity of the DMN-DAN at baseline, which is in contrast to previous reports in similar populations [17,18,19]. These discrepancies may be related to sample size, methodological differences (regional versus network-based analyses), and the larger age range of the participants in our study.

Interestingly, the results found in our study seem to be domain-specific for memory. This may be explained by the DMN being specifically associated with memory functions [10, 46]. Since we found the most convincing results in the DMN-DAN correlation, this may emphasize the domain-specific nature of the interaction between DMN-DAN correlation, amyloid burden, and memory decline in preclinical AD.

Limitations

Due to the observational nature of our study, our data do not allow for any causal inferences between changes in functional connectivity between networks or Aβ deposition. Determining the temporal direction of the three-way association between amyloid, memory decline, and inter-network functional connectivity requires further experimental investigation. Furthermore, longitudinal data of the AD group are limited to a follow-up period of 2 years. This may partially explain the lack of significant results found in this group, although it is also possible that the disease process and associated pathological accumulations are too widespread and may have affected functional networks in multiple ways, all associated with cognitive decline.

Future directions

With the recent development of tau positron emission tomography tracers, future studies can investigate how tau pathology may impact inter-network connectivity since research has suggested a dynamic influence of tau on connectivity within networks, depending on Aβ levels [49,50,51]. Results from such a study could ultimately culminate in clinical trials where participants can be chosen using multi-modal selective criteria for memory decline due to AD.

Additionally, future studies should investigate whether carriers of the apolipoprotein ε4 allele show stonger synergistic effects between AD pathology [52, 53] and inter-network correlations on cognitive decline. Such investigations may further refine the selection criteria for trials and may have implications for determining response to treatment.

Finally, previous research has shown that AD variants, such as posterior cortical atrophy and early-onset AD, may have unique networks that are preferentially affected in the disease process [54]. Thus, future research may want to investigate how Aβ, inter-network correlations, and cognitive decline are associated in these different variants of AD and whether inter-network correlations may be associated with different cognitive domains in these variants.

Conclusions

In conclusion, our results show that the direction of inter-network connectivity provides additional information to baseline Aβ burden about the rate of expected memory decline. These findings add important information for the understanding of factors contributing to cognitive decline. These results also suggest that when using functional connectivity as a biomarker or selection criterion for trials in preclinical populations, the directionality of inter-network connectivity might aid in selecting individuals that are more likely to be on an Aβ-related negative memory trajectory. Including information about both Aβ burden and inter-network functional connectivity can improve recruitment strategies and decrease the time to determine the efficacy of clinical trials in the asymptomatic phase of the disease.

Abbreviations

Aβ:

Amyloid-beta

AD:

Alzheimer’s disease

ADNI:

Alzheimer’s Disease Neuroimaging Initiative

DAN:

Dorsal attention network

DMN:

Default mode network

FDR:

False-discovery rate

FPN:

Frontoparietal network

MCI:

Mild cognitive impairment

ROI:

Region of interest

SN:

Salience network

References

  1. McKhann G, Knopman DS, Chertkow H, Hyman BT, Jack CR Jr, 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.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Jagust WJ, Mormino EC. Lifespan brain activity, beta-amyloid, and Alzheimer’s disease. Trends Cogn Sci. 2011;15:520–6.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Jack CR, Knopman DS, Jagust WJ, Petersen RC, Weiner MW, Aisen PS, et al. Tracking pathophysiological processes in Alzheimer's disease: an updated hypothetical model of dynamic biomarkers. Lancet Neurol. 2013;12:207–16.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  4. Price JL, Davis P, Morris J, White D. The distribution of tangles, plaques and related immunohistochemical markers in healthy aging and Alzheimer’s disease. Neurobiol Aging. 1991;12:295–312.

    Article  PubMed  CAS  Google Scholar 

  5. Mulder C, Verwey NA, van der Flier WM, Bouwman FH, Kok A, van Elk EJ, et al. Amyloid-β (1–42), total tau, and phosphorylated tau as cerebrospinal fluid biomarkers for the diagnosis of Alzheimer disease. Clin Chem. 2010;56:248–53.

    Article  PubMed  CAS  Google Scholar 

  6. Braak H, Braak E. Neuropathological staging of Alzheimer-related changes. Acta Neuropathol. 1991;82:239–59.

    Article  PubMed  CAS  Google Scholar 

  7. Farrell ME, Kennedy KM, Rodrigue KM, et al. Association of longitudinal cognitive decline with amyloid burden in middle-aged and older adults: evidence for a dose-response relationship. JAMA Neurol. 2017;74:830–8.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Boyle PA, Wilson RS, Yu L, Barr AM, Honer WG, Schneider JA, et al. Much of late life cognitive decline is not due to common neurodegenerative pathologies. Ann Neurol. 2013;74:478–89.

    Article  PubMed  Google Scholar 

  9. Gauthier S, Albert M, Fox N, Goedert M, Kivipelto M, Mestre-Ferrandiz J, et al. Why has therapy development for dementia failed in the last two decades? Alzheimers Dement. 2016;12:60–4.

    Article  PubMed  Google Scholar 

  10. Damoiseaux JS, Beckmann CF, Arigita EJS, Barkhof F, Scheltens P, Stam CJ, et al. Reduced resting-state brain activity in the “default network” in normal aging. Cereb Cortex. 2008;18:1856–64.

    Article  PubMed  CAS  Google Scholar 

  11. Andrews-Hanna JR, Snyder AZ, Vincent JL, Lustig C, Head D, Raichle Marcus E, et al. Disruption of large-scale brain systems in advanced aging. Neuron. 2007;56:924–35.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Buckley RF, Schultz AP, Hedden T, Papp KV, Hanseeuw BJ, Marshall G, et al. Functional network integrity presages cognitive decline in preclinical Alzheimer disease. Neurology. 2017;89(1):29–37.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Fornito A, Harrison BJ, Zalesky A, Simons JS. Competitive and cooperative dynamics of large-scale brain functional networks supporting recollection. Proc Natl Acad Sci U S A. 2012;109:12788–93.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Kelly AM, Uddin LQ, Biswal BB, Castellanos FX, Milham MP. Competition between functional brain networks mediates behavioral variability. NeuroImage. 2008;39:527–37.

    Article  PubMed  Google Scholar 

  15. Fox MD, Snyder AZ, Vincent JL, Corbetta M, Van Essen DC, Raichle ME. The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc Natl Acad Sci U S A. 2005;102:9673–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Pinsk MA, Kastner S. Unconscious networking. Nature. 2007;447:46.

    Article  PubMed  CAS  Google Scholar 

  17. Wang K, Liang M, Wang L, Tian L, Zhang X, Li K, et al. Altered functional connectivity in early Alzheimer’s disease: a resting-state fMRI study. Hum Brain Mapp. 2007;28:967–78.

    Article  PubMed  Google Scholar 

  18. Esposito R, Cieri F, Chiacchiaretta P, Cera N, Lauriola M, Di Giannantonio M, et al. Modifications in resting state functional anticorrelation between default mode network and dorsal attention network: comparison among young adults, healthy elders and mild cognitive impairment patients. Brain Imaging Behav. 2018;12:127–41.

    Article  PubMed  Google Scholar 

  19. Brier MR, Thomas JB, Snyder AZ, Benzinger TL, Zhang D, Raichle ME, et al. Loss of intranetwork and internetwork resting state functional connections with Alzheimer’s disease progression. J Neurosci. 2012;32:8890–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. Mormino EC, Betensky RA, Hedden T, Schultz AP, Amariglio RE, Rentz DM, et al. Synergistic effect of beta-amyloid and neurodegeneration on cognitive decline in clinically normal individuals. JAMA Neurol. 2014;71:1379–85.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Mormino EC, Smiljic A, Hayenga AO, Onami SH, Greicius MD, Rabinovici GD, et al. Relationships between beta-amyloid and functional connectivity in different components of the default mode network in aging. Cereb Cortex. 2011;21:2399–407.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Sheline YI, Raichle ME, Snyder AZ, Morris JC, Head D, Wang S, et al. Amyloid plaques disrupt resting state default mode network connectivity in cognitively normal elderly. Biol Psychiatry. 2010;67:584–7.

    Article  PubMed  CAS  Google Scholar 

  23. Mueller SG, Weiner MW, Thal LJ, Petersen RC, Jack C, Jagust W, et al. The Alzheimer’s disease neuroimaging initiative. Neuroimaging Clin N Am. 2005;15:869–77.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Landau SM, Mintun MA, Joshi AD, Koeppe RA, Petersen RC, Aisen PS, et al. Amyloid deposition, hypometabolism, and longitudinal cognitive decline. Ann Neurol. 2012;72:578–86.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Petersen RC. Mild cognitive impairment as a diagnostic entity. J Intern Med. 2004;256:183–94.

    Article  PubMed  CAS  Google Scholar 

  26. Petersen RC, Aisen P, Beckett LA, Donohue M, Gamst A, Harvey DJ, et al. Alzheimer’s disease neuroimaging initiative (ADNI) clinical characterization. Neurology. 2010;74:201–9.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 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.

    Article  PubMed  CAS  Google Scholar 

  28. Crane PK, Carle A, Gibbons LE, Insel P, Mackin RS, Gross A, et al. Development and assessment of a composite score for memory in the Alzheimer’s Disease Neuroimaging Initiative (ADNI). Brain Imaging Behav. 2012;6:502–16.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Gibbons LE, Carle AC, Mackin RS, Harvey D, Mukherjee S, Insel P, et al. A composite score for executive functioning, validated in Alzheimer’s Disease Neuroimaging Initiative (ADNI) participants with baseline mild cognitive impairment. Brain Imaging Behav. 2012;6:517–27.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Whitfield-Gabrieli S, Nieto-Castanon A. Conn: a functional connectivity toolbox for correlated and anticorrelated brain networks. Brain connectivity. 2012;2:125–41.

    Article  PubMed  Google Scholar 

  31. Murphy K, Birn RM, Handwerker DA, Jones TB, Bandettini PA. The impact of global signal regression on resting state correlations: are anti-correlated networks introduced? NeuroImage. 2009;44:893–905.

    Article  PubMed  Google Scholar 

  32. Behzadi Y, Restom K, Liau J, Liu TT. A component-based noise correction method (CompCor) for BOLD and perfusion-based fMRI. NeuroImage. 2007;37:90–101.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Field A. Discovering statistics using IBM SPSS statistics. London: SAGE Publications Ltd; 2013.

  34. Ibrahim JG, Molenberghs G. Missing data methods in longitudinal studies: a review. Test (Madr). 2009;18:1–43.

    Article  Google Scholar 

  35. Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team. nlme: linear and nonlinear mixed effects models. R package version 3.1–128. 2016. http://CRAN.R-project.org/package=nlme.

  36. Bischof GN, Rodrigue KM, Kennedy KM, Devous MD, Park DC. Amyloid deposition in younger adults is linked to episodic memory performance. Neurology. 2016;87:2562–6.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Whitwell JL, Crum WR, Watt HC, Fox NC. Normalization of cerebral volumes by use of intracranial volume: implications for longitudinal quantitative MR imaging. Am J Neuroradiol. 2001;22:1483–9.

    PubMed  CAS  Google Scholar 

  38. Fischl B. FreeSurfer. Neuroimage. 2012;62:774–81.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Buckner RL, Head D, Parker J, Fotenos AF, Marcus D, Morris JC, et al. A unified approach for morphometric and functional data analysis in young, old, and demented adults using automated atlas-based head size normalization: reliability and validation against manual measurement of total intracranial volume. NeuroImage. 2004;23:724–38.

    Article  PubMed  Google Scholar 

  40. Hanseeuw BJ, Betensky RA, Mormino EC, Schultz AP, Sepulcre J, Becker JA, et al. PET staging of amyloidosis using striatum. Alzheimers Dement. 2018. https://doi.org/10.1016/j.jalz.2018.04.011.

  41. Elman JA, Madison CM, Baker SL, Vogel JW, Marks SM, Crowley S, et al. Effects of beta-amyloid on resting state functional connectivity within and between networks reflect known patterns of regional vulnerability. Cereb Cortex. 2016;26:695–707.

    PubMed  Google Scholar 

  42. Mucke L, Selkoe DJ. Neurotoxicity of amyloid β-protein: synaptic and network dysfunction. Cold Spring Harb Perspect Med. 2012;2:a006338.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Grothe MJ, Teipel SJ. Spatial patterns of atrophy, hypometabolism, and amyloid deposition in Alzheimer’s disease correspond to dissociable functional brain networks. Hum Brain Mapp. 2016;37:35–53.

    Article  PubMed  Google Scholar 

  44. Keller JB, Hedden T, Thompson TW, Anteraper SA, Gabrieli JD, Whitfield-Gabrieli S. Resting-state anticorrelations between medial and lateral prefrontal cortex: association with working memory, aging, and individual differences. Cortex. 2015;64:271–80.

    Article  PubMed  Google Scholar 

  45. Cabeza R, Anderson ND, Locantore JK, McIntosh AR. Aging gracefully: compensatory brain activity in high-performing older adults. NeuroImage. 2002;17:1394–402.

    Article  PubMed  Google Scholar 

  46. Franzmeier N, Buerger K, Teipel S, Stern Y, Dichgans M, Ewers M, et al. Cognitive reserve moderates the association between functional network anti-correlations and memory in MCI. Neurobiol Aging. 2017;50:152–62.

    Article  PubMed  Google Scholar 

  47. Spreng RN, Stevens WD, Viviano JD, Schacter DL. Attenuated anticorrelation between the default and dorsal attention networks with aging: evidence from task and rest. Neurobiol Aging. 2016;45:149–60.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Damoiseaux JS, Prater KE, Miller BL, Greicius MD. Functional connectivity tracks clinical deterioration in Alzheimer's disease. Neurobiol Aging. 2012;33:828. e19-.e30

    Article  PubMed  Google Scholar 

  49. Schultz AP, Chhatwal JP, Hedden T, Mormino EC, Hanseeuw BJ, Sepulcre J, et al. Phases of hyperconnectivity and hypoconnectivity in the default mode and salience networks track with amyloid and tau in clinically normal individuals. J Neurosci. 2017;37:4323–31.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  50. Sepulcre J, Sabuncu MR, Li Q, El Fakhri G, Sperling R, Johnson KA. Tau and amyloid β proteins distinctively associate to functional network changes in the aging brain. Alzheimers Dement. 2017;13:1261–9.

    Article  PubMed  Google Scholar 

  51. Jacobs HI, Hedden T, Schultz AP, Sepulcre J, Perea RD, Amariglio RE, et al. Structural tract alterations predict downstream tau accumulation in amyloid-positive older individuals. Nat Neurosci. 2018;21:424.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  52. Corder EH, Saunders AM, Strittmatter WJ, Schmechel DE, Gaskell PC, Small G, et al. Gene dose of apolipoprotein E type 4 allele and the risk of Alzheimer's disease in late onset families. Science. 1993;261:921–3.

    Article  PubMed  CAS  Google Scholar 

  53. Strittmatter WJ, Saunders AM, Schmechel D, Pericak-Vance M, Enghild J, Salvesen GS, et al. Apolipoprotein E: high-avidity binding to beta-amyloid and increased frequency of type 4 allele in late-onset familial Alzheimer disease. Proc Natl Acad Sci U S A. 1993;90:1977–81.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  54. Lehmann M, Madison CM, Ghosh PM, Seeley WW, Mormino E, Greicius MD, et al. Intrinsic connectivity networks in healthy subjects explain clinical variability in Alzheimer’s disease. Proc Natl Acad Sci U S A. 2013;110:11606–11.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to thank Ed H.B.M. Gronenschild, PhD (Faculty of Health, Medicine and Life Sciences, School for Mental Health and Neuroscience, Department of Psychiatry and Neuropsychology, Alzheimer Center Limburg, Maastricht University, Maastricht, The Netherlands) for his invaluable assistance in the processing of imaging data.

Data used in the preparation of this article were obtained from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) database (adni.loni.usc.edu). As such, the investigators within the ADNI contributed to the design and implementation of ADNI and/or provided data but did not participate in the analysis or writing of this report. A complete listing of ADNI investigators can be found at https://adni.loni.usc.edu/wp-content/uploads/how_to_apply/ADNI_Acknowledgement_List.pdf.

Funding

HILJ is supported by funding from NWO (Veni—451–14-035), the European Union’s Horizon 2020 Research and Innovation Programme under the Marie Sklodowska-Curie Grant agreement (IF-2015-GF, 706714), a standard grant of the International Stichting voor Alzheimer Onderzoek (#15007) and a research grant by the Deutsche Forschungsgemeinschaft (JA/2336/1–1). ADNI is funded through generous contributions from the following: AbbVie, Alzheimer’s Association; Alzheimer’s Drug Discovery Foundation; Araclon Biotech; BioClinica, Inc.; Biogen; Bristol-Myers Squibb Company; CereSpir, Inc.; Cogstate; Eisai Inc.; Elan Pharmaceuticals, Inc.; Eli Lilly and Company; EuroImmun; F. Hoffmann-La Roche Ltd. and its affiliated company Genentech, Inc.; Fujirebio; GE Healthcare; IXICO Ltd.; Janssen Alzheimer Immunotherapy Research & Development, LLC.; Johnson & Johnson Pharmaceutical Research & Development LLC.; Lumosity; Lundbeck; Merck & Co., Inc.; Meso Scale Diagnostics, LLC.; NeuroRx Research; Neurotrack Technologies; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Piramal Imaging; Servier; Takeda Pharmaceutical Company; and Transition Therapeutics. The Canadian Institutes of Health Research is providing funds to support ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the National Institutes of Health (www.fnih.org). The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimer’s Therapeutic Research Institute at the University of Southern California. ADNI data are disseminated by the Laboratory for Neuro Imaging at the University of Southern California and funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering. Data collection and sharing for this project was funded by the Alzheimer’s Disease Neuroimaging Initiative (ADNI) (National Institutes of Health Grant U01 AG024904) and DOD ADNI (Department of Defense award number W81XWH-12-2-0012). The funding sources had no involvement related to this manuscript in any way.

Availability of data and materials

The datasets generated and/or analyzed during the current study are available in the ADNI repository, http://adni.loni.usc.edu/data-samples/access-data/.

Author information

Authors and Affiliations

Authors

Consortia

Contributions

All authors contributed to the conception and design of the study. RWEVH and HILJ performed analyses on the data. All authors interpreted the data. RWEVH drafted the initial version of the manuscript. JMR and HILJ revised the manuscript for important intellectual content, and all authors read and approved the final manuscript.

Corresponding author

Correspondence to Roy W. E. Van Hooren.

Ethics declarations

Ethics approval and consent to participate

ADNI is a multisite investigation. All ADNI sites obtained Institutional Review Board approval for the use of humans for research from their respective Institutional Review Board within 5 months of receiving the final protocol and all participants provided written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional file

Additional file 1:

Table S1. List of participant IDs included in final analyses. Table S2. Peak coordinates of regions of interest in the MNI space. Box 1. Regression model construction. Table S3. Tukey tests for differences of baseline characteristics. Table S4. Summary of maximum motion parameter values and scrubbed volumes across groups. Table S5. Baseline results. Table S6. Associations between functional inter-network connectivity and memory decline in the whole sample. Table S7. Longitudinal linear effects of inter-network functional connectivity between each network pair on memory, per group. Table S8. Associations between group, inter-network functional connectivity, and decline in executive functions. Table S9. Associations between functional connectivity, time, amyloid status, memory, and executive functioning in the cognitively normal group. (DOCX 68 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Van Hooren, R.W.E., Riphagen, J.M., Jacobs, H.I.L. et al. Inter-network connectivity and amyloid-beta linked to cognitive decline in preclinical Alzheimer’s disease: a longitudinal cohort study. Alz Res Therapy 10, 88 (2018). https://doi.org/10.1186/s13195-018-0420-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13195-018-0420-9

Keywords