Resting-state functional MRI signal fluctuations are correlated with brain amyloid-β deposition

Norman Scheel, Takashi Tarumi, Tsubasa Tomoto, C. Munro Cullum, Rong Zhang, David C. Zhu
doi: https://doi.org/10.1101/2021.04.22.21255924

Abstract

Mounting evidence suggests that amyloid-β (Aβ) and vascular etiologies are intertwined in the pathogenesis of Alzheimer’s disease. Spontaneous fluctuations of the brain blood-oxygen-level-dependent (BOLD) signal, as measured by resting-state functional MRI (rs-fMRI), have been shown to be associated with neuronal activities as well as cerebrovascular hemodynamics. Nevertheless, it is unclear if rs-fMRI BOLD fluctuations are associated with brain Aβ deposition in individuals with an elevated risk of Alzheimer’s disease.

We recruited 33 patients with amnestic mild cognitive impairment who underwent rs-fMRI and positron emission tomography (PET). The Aβ standardized uptake value ratio (SUVR) was calculated with cortical white matter as the reference region to improve sensitivity for cortical Aβ quantification. We calculated the amplitudes of low-frequency fluctuations (ALFF) of local BOLD signals in the frequency band of 0.01-0.08 Hz. Applying physiological/vascular signal regression in stepwise increasing levels on the rs-fMRI data, we examined whether local correlations between ALFF and brain Aβ deposition were driven by vascular hemodynamics, spontaneous neuronal activities, or both.

We found that ALFF and Aβ SUVR were negatively correlated in brain regions involving the default-mode and visual networks, with peak correlation at the precuneus, and angular, lingual, and fusiform gyri. Regions with higher ALFF had less Aβ accumulation. The correlated cluster sizes in MNI space were reduced from 3018 mm3 with no physiological/vascular regression to 1072 mm3 with strong physiological/vascular regression, with mean cluster r values at approximately -0.47.

Results demonstrate that both vascular hemodynamics and neuronal activities, as reflected by BOLD fluctuations, are negatively associated with brain Aβ deposition. These findings further imply that local brain blood fluctuations due to either vascular hemodynamics or neuronal activities can affect Aβ homeostasis.

Introduction

Alzheimer’s disease is a major neurodegenerative brain disorder and the most common type of dementia among older adults.1 Patients with Alzheimer’s disease exhibit memory loss, impaired decision-making capacity, disorientation, as well as changes in daily activity.2 These cognitive deficits and functional changes are accompanied and often preceded by brain pathophysiological changes including amyloid-β (Aβ) and tau protein depositions in the cerebral cortex, brain hypoperfusion, and neuronal degeneration.3 Mild cognitive impairment, particularly amnestic mild cognitive impairment, has been considered to represent a transitional phase between cognitive aging and Alzheimer’s disease, as many who develop amnestic mild cognitive impairment progress to Alzheimer’s disease.4

Two major hypotheses have been developed to understand the pathophysiological mechanisms of Alzheimer’s disease and to develop effective treatments. In the classic amyloid cascade hypothesis, the disruption of brain Aβ homeostasis has been proposed to be a primary driver of Alzheimer’s disease, leading to accumulations of Aβ plaques and neurofibrillary tangles, and consequently neurodegeneration and cognitive impairment.5 Conversely, the vascular hypothesis of Alzheimer’s disease proposes that cerebrovascular dysfunction is an important contributor to Alzheimer’s disease onset and progression.6–8 Mounting evidence indicates that Aβ and vascular etiology are intertwined in the pathogenesis of Alzheimer’s disease.9,10 In a recent study, we found that age-related carotid artery stiffness was associated positively with brain Aβ burden in patients with amnestic mild cognitive impairment, supporting the role of arterial aging in brain Aβ deposition.11 Recent animal and human studies also suggested that cerebral arterial oscillations and vasomotions influence brain Aβ clearance and deposition.12–14 Spontaneous cerebral blood flow (CBF) fluctuations are linked intrinsically to brain neuronal activities through neurovascular coupling.15–17 This phenomenon has been exploited by using resting-state functional MRI (rs-fMRI) to understand functional connectivities between brain regions and the disruption of functional connectivity in diseased states.16,18–26 Specifically, it has been proposed that spontaneous neuronal activities, which consume approximately 95% of the brain’s metabolism,16 are manifested as rs-fMRI blood-oxygen-level-dependent (BOLD) signal fluctuations.16,17,27

It has been widely recognized that spontaneous BOLD fluctuations are influenced by both systemic and cerebrovascular hemodynamics in addition to neuronal activities.28–32 In this regard, our previous work has shown that spontaneous fluctuations of BOLD signals and CBF velocity measured in the middle cerebral artery (upper-stream vascular signals) had similar spectral distributions and high correlation in the 0.01 - 0.08 Hz frequency band.33 These observations suggested that upper-stream cerebrovascular hemodynamic fluctuations can transmit downstream into the cerebral microcirculation and contribute importantly to regional BOLD fluctuations. On the other hand, local cerebral vasomotion, which has been linked to brain Aβ clearance, also may influence rs-fMRI BOLD fluctuations.12 A recent study also showed increased fMRI signal variability at cardiorespiratory frequencies in Alzheimer’s disease patients when compared to cognitively normal performing older adults.34 This study aimed to determine whether BOLD signal fluctuations as quantified by the amplitude of low-frequency fluctuations (ALFF) in the range of 0.01-0.08 Hz are related to brain Aβ deposition in patients with amnestic mild cognitive impairment.4,35–38 Furthermore, we examined whether the associations between local ALFF and brain Aβ depositions are influenced by the estimated upper-stream cerebrovascular hemodynamics, spontaneous neuronal activities, or both, using the well-established signal regression approaches.39–41 We hypothesized that local ALFF would be negatively correlated with brain Aβ deposition and that correlation strength between ALFF and Aβ deposition would decrease when the upper-stream vascular signals are regressed out from the raw BOLD signal under an assumption that upper-stream vascular effects drive brain Aβ clearance.

Materials and methods

Participants

Thirty-three amnestic mild cognitive impairment subjects (64.4 ± 6.4 years of age, 19 females) were recruited who represent a subgroup of participants enrolled in a proof of concept investigation aimed to identify effects of exercise training on neurocognitive function in amnestic mild cognitive impairment (ClinicalTrials.gov, NCT01146717).42 Data presented here were obtained at baseline before any intervention. The diagnosis of amnestic mild cognitive impairment was based on the Petersen criteria as modified by the Alzheimer’s Disease Neuroimaging Initiative (ADNI) project (http://adni-info.org).43 Diagnostic assessments included the Clinical Dementia Rating (CDR) scale,44 the Mini-Mental State Examination (MMSE),45 and the Wechsler Memory Scale-Revised Logical Memory (LM) immediate and delayed recall trials.46 Heart rate (HR) and brachial cuff blood pressure (BP) were measured >3 times with an ECG-gated electro-sphygmomanometer (Suntech, Morrisville, NC, USA). Obtained values were averaged to obtain HR, systolic (SBP), and diastolic BP (DBP). Participants with major psychiatric disorders, major or unstable medical conditions, uncontrolled hypertension, diabetes mellitus, or chronic inflammatory diseases were excluded, as were participants with a cardiac pacemaker or any metal plates or pins in their body preventing them from undergoing MRI scanning (detailed inclusion and exclusion criteria are provided in ClinicalTrials.gov NCT01146717). All participants gave informed consent. This study was approved by the Institutional Review Boards of the University of Texas Southwestern Medical Center and Texas Health Presbyterian Hospital of Dallas and performed by the guidelines of the Declaration of Helsinki and Belmont Report.

MRI Measurements

Rs-fMRI data were collected on a Philips Achieva 3T scanner (Philips Healthcare, Best, the Netherlands) under an “eye-closed” condition with the following parameters: gradient recalled echo planner imaging (EPI), 29 contiguous 5-mm axial slices, 30-ms time of echo (TE), 1500-ms time of repetition (TR), 60° flip angle, a 24-cm field of view, 80 × 80 matrix size and 200 time-points. High-resolution T1-weighted 3D MPRAGE (Magnetization-Prepared Rapid Acquisition Gradient-Echo) images were also collected to cover the whole brain with the following parameters: TE/TR = 3.7/8.1 ms, flip angle = 12°, field of view = 256 mm × 204 mm, 160 1-mm slices, resolution = 1 mm × 1 mm × 1 mm, SENSE factor = 2 and scan duration = 4 minutes.

PET image acquisition

After an intravenous bolus injection of 10 mCi 18F-florbetapir (also called AV45), participants were positioned in a Siemens (Munich, Germany) ECAT HR PET scanner for data acquisition, using laser guidance for precise head positioning. Velcro straps and foam wedges were used to secure the participant’s head. To ensure the brain was completely in the field of view and absent of rotation in either the transverse or sagittal planes, a 2-minute scout scan was acquired. At 50 minutes post-injection, 2 frames of 5-minute PET emission scan and a 7-minute transmission scan were acquired in 3D mode using the following parameters: matrix size = 128 × 128, resolution = 5 mm × 5 mm, slice thickness = 2.42 mm, and field of view = 58.3 cm. Emission images were processed by iterative reconstruction, 4 iterations, and 16 subsets with a 3-mm full width at half maximum (FWHM) ramp filter. The transmission image was reconstructed using back-projection and a 6-mm FWHM Gaussian filter for attenuation correction 47,48.

RS-fMRI BOLD signal fluctuation quantification

The quantification of rs-fMRI BOLD local signal fluctuation using ALFF was first proposed by Zang et al..38 It was first used to investigate regional rs-fMRI signal fluctuations in attention deficit hyperactivity disorder. To calculate ALFF, rs-fMRI recordings are first pre-processed according to current standards, including slice-timing correction, motion correction, and spatial blurring (different strategies used here are detailed later). Importantly, no temporal filtering is applied before ALFF calculation. Fast Fourier transformation (FFT) is carried out on each preprocessed voxel time course. The ALFF index is calculated as two times the amplitude of the FFT output, divided by the sample length (number of data points in the time course), then averaged across the specified frequency band of interest.

A frequency band of 0.01 - 0.08 Hz was used in this study. Figure 1(a) shows a typical pre-processed rs-fMRI time course at a brain voxel. Figure 1(b) shows the corresponding Fourier spectrum and the frequency band for ALFF computation. ALFF at each voxel is subsequently normalized by the global mean ALFF from all voxels within the brain to generate a standardized ALFF index (sALFF) that allows comparisons across participants. ALFF computations were carried out using the ALFF function included in the DPARSFA toolbox from the “DPABI V4.3_200401” software package49,50 in MATLAB R2019a (MathWorks, Natick, MA, USA).

Figure 1
Figure 1 Time series and spectral analysis of spontaneous BOLD signal fluctuations

(a) time course at a voxel in the occipital region of one study participant; (b) corresponding Fourier spectrum. Across the applied frequency band of 0.01-0.08 Hz, the amplitude of low-frequency fluctuations (ALFF) yields 605.7 as of this voxel’s fluctuation amplitude (plotted in the green shade).

Anatomical image pre-processing and MNI-registration

One common challenge in brain imaging studies of older adults is the high degree of anatomical variability due to atrophy, leading to difficulties in performing group-wise image analyses. FreeSurfer51, an automated data processing software for cortical surface reconstruction and anatomical segmentation of brain MRI scans, provides a robust brain segmentation in comparable study cohorts, as demonstrated in prior studies.26 The high-resolution T1-weighted volumetric MR images of each participant were processed using the default FreeSurfer “recon-all” pipeline. This pipeline creates a robust non-linear transformation matrix that warps images from a participant’s native space to MNI standard space, according to the MNI305 template resolution (1 mm × 1 mm × 1 mm with a matrix size of 256 × 256 × 256). This transformation matrix was applied to warp the results of the sALFF and AV45-SUVR calculations to the MNI305 template. While the transformation to the standard template performed well overall for brains with atrophy, artifacts in subcortical regions were noticeable due to non-linear stretching, especially for brains with enlarged ventricles. These artifacts generated spatial noise in our data analyses. Subsequently, for computational efficiency, all images in MNI305 template resolution were resampled to fit the MNI152 template with a resolution of 2 mm × 2 mm × 2 mm at a matrix size of 91 × 109 × 91.

fMRI preprocessing and physiological/vascular signal regression

Resting-state fMRI pre-processing was first carried using AFNI software 52 in native space. The “afni_proc.py” routine in AFNI was used to generate the script to pre-process the rs-fMRI data. For each participant, any signal spikes in the signal intensity time courses were first detected and removed. The acquisition timing difference was then corrected for different slice locations. With the third functional volume as the reference, rigid-body motion correction was carried out in three translational and three rotational directions. The displacements due to motion in these six directions as well as the corresponding motion derivatives at each time point were estimated and then modeled as motion regressors. The data points with excess motion (normalized motion derivative□>□0.5 or voxel outliers□>□10%) were identified to generate another motion regressor. For each participant, spatial blurring with a full width half maximum (FWHM) of 4 mm was used to reduce random noise. At each voxel, motion-related signal changes, baseline, linear, 2nd-order, and 3rd-order system-induced signal trends were modeled as linear regressors. The “3dDeconvolve” function in AFNI was applied to remove the motion and system introduced noises from each voxel’s time course, as modeled in the regressors described above. Up to this step, physiological signals, likely due to upper-stream vascular effects, have not yet been regressed from the signal. The resulting fMRI dataset up to this step will be called “No Physiological Regression” to emphasize that no attempt has been made to remove upper-stream vascular effects.

The standardized amplitude of the BOLD signal fluctuations (sALFF) within the frequency range of 0.01 Hz to 0.08 Hz was then calculated for each voxel within the brain. For group analyses, subject-space sALFF maps were non-linearly transformed into MNI305 standard space, using the transformation matrix created by the FreeSurfer processing pipeline as discussed above. MNI305 sALFF images were subsequently resampled to fit the MNI152 standard template for improved computational efficiency. The correlation between sALFF and brain amyloid β accumulation, calculated as AV45 PET amyloid β standardized uptake value ratio (SUVR), was then calculated in MNI152 standard space as discussed further in the following section. To understand the upper-stream physiological/vascular effects, the “no physiological regression” dataset was further processed with three progressively more aggressive physiological signal regression procedures. For each procedure, sALFF was calculated again followed by the aforementioned standard-space transformations and correlation analyses. The three regression procedures, in addition to the “no physiological regression”, are listed below:

  1. “No physiological regression:” as described above, this procedure has no special treatment to remove upper-stream physiological/vascular effects before sALFF computation.

  2. “WM/CSF regression:” White matter (WM) and cerebrospinal fluid (CSF) do not contain neuronal brain activity. Consequently, fluctuations of mean BOLD signal time courses at these regions can be assumed to mostly represent upper-stream vascular effects.53 Hence, CSF and WM signals were modeled and used as additional regressors in the “3dDeconvolve” step to reduce the upper-stream vascular influence, before sALFF computation.

  3. “WM/CSF/GS regression:” The mean global signal (GS) was assumed to be primarily driven by an overall upper-stream vascular effect and added as an additional regressor to this “WM/CSF regression” procedure for further physiological signal removal before sALFF computation. With this procedure, the resulting signal time course should closely reflect neuronal brain activity.53 While global signal removal can potentially remove some neuronal signals from fMRI data, current research suggests that the upper-stream vascular components most likely dominate the global signal.54– 56 Recently, Xifra-Porxas and colleagues reported that physiological signals accounted for an r of about 0.6 within the global signal.57 Here, global signal regression was used to further investigate the effect of the upper-stream vascular signals on Aβ deposition in the brain and its influence in this cascaded setup.

  4. “Aggressive AROMA:” A rs-fMRI noise removing technique called “ICA-AROMA” was developed by Pruim and colleagues39 to aggressively remove motion artifacts and physiological noise. The authors demonstrated that this technique could effectively identify motion and physiological signals as independent components and then remove them from rs-fMRI data.58 “ICA-AROMA” was implemented using FSL’s MELODIC routine 59. In this technique, single-subject spatial independent component analysis (ICA) is performed on a motion-corrected whole-brain rs-fMRI dataset in native space. Using a pre-trained classifier, the resulting components are sorted into noise and non-noise components. The time courses of components identified as noise are subsequently regressed from the data using FSL’s “regfilt” function, where partial regression is called “non-aggressive AROMA,” and full regression is called “aggressive AROMA.” The “aggressive AROMA” approach was employed here to maximize the effect of removing upper-stream vascular signals.

The overall signal processing steps described above are also presented in supplementary Figure S1.

PET image processing

Using “flirt” from the FSL software package,59 a six-degree linear transformation was applied to align PET AV45 SUVR images to the 3D MPRAGE anatomical MRI images. All alignments were visually inspected and improved with manual adjustments if necessary.

Using the segmentation of the high-resolution 3D MPRAGE image from the FreeSurfer analysis, white-matter (WM) and cerebellar regions were isolated. An additional erosion of the white-matter map, using FSL’s “fslmaths” with a box kernel of 4 mm × 4 mm × 4 mm ensured that only white-matter regions were obtained. The AV45 uptake images were then normalized by the mean uptake at the white matter and also the cerebellar regions to generate two versions of SUVR images. Previous studies have shown that SUVR calculation using WM signal can improve discriminatory power for detecting Alzheimer’s disease and mild cognitive impairment.60 In the following results, SUVR will be reported for WM normalization, yet posthoc analyses with the AV45-SUVR images normalized by cerebellar signal were also carried out to confirm observations. Congruent to the spatial normalization of sALFF maps, AV45-SUVR images were non-linearly transformed to match the MNI305 standard template, using the transformation matrix created by the FreeSurfer processing pipeline. AV45-SUVR images in MNI305 template space were subsequently resampled to fit the MNI152 standard template for better computational efficiency.

Statistical analysis

Correlations between sALFF and AV45-SUVR across participants were carried out for each brain voxel in MNI152 template space. Monte-Carlo simulations, using 10,000 random permutations for each correlation value r were applied to approximate the corresponding p values, through constructing voxel-specific empirical cumulative distributions.61 These computations were realized in MATLAB R2019a (MathWorks, Natick, MA, USA). Cluster-level multiple-comparison corrections were computed using AFNI’s52 “3dClustSim” Monte-Carlo simulation software. In these simulations, we estimated the minimum cluster sizes required to achieve a corrected p ≤ 0.05 at a voxel-level p < 0.001. As the estimated minimum cluster size is highly dependent on data spatial smoothness, which in turn depends on data preprocessing, cluster size thresholds were estimated independently for each physiological/vascular signal removal procedure, described earlier. Autocorrelation function parameters (-acf) employed by “3dClustSim” were estimated using “3dFWHMx” (AFNI) on the preprocessed fMRI data (MNI transformed and resampled to 91 x 109 x 91 isotropic 2-mm voxels). This procedure yields mean parameters (across participants) for each of the four different autocorrelation functions, specific to each fMRI preprocessing procedure. We also estimated the spatial autocorrelation of the AV45-SUVR maps using “3dFWHMx”. Our pre-processed fMRI images showed less auto-correlation than the AV45-SUVR maps. To be conservative, we used the autocorrelation estimates of the fMRI images to calculate the required cluster sizes to achieve statistical significance. Subsequently, “3dClusterize” (AFNI) was used to generate correlation cluster maps that survive the estimated cluster correction thresholds for each preprocessing procedure.

The preprocessing-specific significant clusters were then used to assess whether the correlations between BOLD signal fluctuations and brain amyloid deposition were rather influenced by upper-stream physiological/vascular effects or neuronal activations. For this purpose, we illustrated the distributions of the Fisher z transformed r values of the AV45-SUVR-sALFF correlations, using the voxel populations identified as significant clusters resulting from the analyses based on the “no physiological regression” and the “aggressive AROMA” procedures.

Data availability

The raw data were generated at the University of Texas Southwestern Medical Center and UTSW and Institute for Exercise and Environmental Medicine Texas Health Presbyterian Hospital, Dallas, TX, USA. Derived data supporting the findings of this study are available from the corresponding authors upon reasonable request.

Results

Table 1 shows the participant characteristics. Overall, our amnestic mild cognitive impairment patients were well educated and showed early amnestic mild cognitive impairment symptoms, as shown by CDR scores of 0.5, a mean MMSE score of 29.5 ± 0.8, and mean LM immediate recall of 11.4 ± 2.3 and delayed recall of 9.7 ± 1.8. These older adults had an average heart rate of 62 bpm and mean systolic blood pressure of 118 mmHg but included hypertension (≥140 mmHg). The presence of apolipoprotein E4 (APOE4) was assessed in 28 participants and 25% were APOE4 positive. The participants had a mean cortical AV45-SUVR of 0.57 ± 0.05 when normalized with the white matter as a reference and 1.15 ± 0.08 with the whole cerebellum as a reference.

View this table:
Table 1 Participant characteristics

Figures 2 and 3 show graphical representations of select sagittal and axial slices of mean AV45-SUVR, mean sALFF for each preprocessing procedure, and the significant correlation clusters between them. Table 2 provides a detailed list of significant correlation clusters with associated brain regions for each preprocessing procedure. Using the “no physiological regression” procedure, cluster level significance was estimated with voxel-level p < 0.001 and a corrected p ≤ 0.05, yielding a minimum cluster size of 493 voxels (2-mm isotropic). Using this threshold, four significant negatively correlated clusters between sALFF and AV45-SUVR were found spanning across the lingual gyrus, precuneus, cuneus, and fusiform gyrus, with the most prominent cluster reaching a size of approximately 3018 mm3 (see Table 2, Figures 2b & 3b).

View this table:
Table 2 Correlation cluster analysis of sALFF and AV45-SUVR
Figure 2
Figure 2 Voxel-wise correlation analysis between brain amyloid-β deposition and the amplitude of low-frequency BOLD signal fluctuations (sALFF) - sagittal.

Images represent slices near the mid-sagittal region (right and left 9.25mm in MNI space). (a) mean AV45-SUVR (standardized uptake value ratio with the white matter as a scaling reference) color maps of 33 amnestic mild cognitive impairment participants are shown along with the corresponding mean sALFF (standardized amplitude of low-frequency fluctuations) color maps after the physiological signal cleaning procedures of (b) “no physiological regression”, (c) “WM/CSF regression”, (d) “WM/CSF/GS regression” and (e) “aggressive AROMA”. Significant correlation clusters between AV45-SUVR and sALFF for each regression procedure (b)-(e) are shown in the row below the mean sALFF images. Clusters are depicted in blue, as r values are negative (see Table 2 for detailed statistical analysis data). AV45-SUVR maps are nearly identical if the cerebellum, instead of white matter, is used as a reference, except the color scale range becomes 0 to 2.

Figure 3
Figure 3 Voxel-wise correlation analysis between brain amyloid-β deposition and the amplitude of low-frequency BOLD signal fluctuations (sALFF) - axial.

Images represent an axial slice (top half) above the corpus callosum (42mm in inferior-superior plane MNI) and an axial slice (bottom half) across the thalamus (0.75mm in inferior-superior plane MNI). (a) mean AV45-SUVR (standardized uptake value ratio with white matter scaling) color maps of 33 amnestic mild cognitive impairment participants are shown, along with the corresponding mean sALFF (standardized amplitude of low-frequency fluctuations) color maps after the physiological signal cleaning procedures of (b) “no physiological regression”, (c) “WM/CSF regression”, (d) “WM/CSF/GS regression”, and (e) “aggressive AROMA”. Significant correlation clusters between AV45-SUVR and sALFF for each regression procedure (b)-(e) are shown in the row below the mean sALFF images. Clusters are depicted in blue, as r values are negative (see Table 2 for detailed statistical analysis data). AV45-SUVR maps are nearly identical if the cerebellum, instead of white matter, is used as a reference, except the color scale range becomes 0 to 2.

After applying the “WM/CSF regression” procedure using the same statistical criteria as above, the corrected p ≤ 0.05 was achieved with a minimum cluster size of 434 voxels. At this threshold, four clusters were identified as significant. The significant brain regions are consistent with the “no physiological regression” methods, spanning the lingual gyrus, fusiform gyri, cuneus, and precuneus, yet at slightly smaller cluster sizes for all except the lingual gyrus cluster which slightly increases in size (see Table 2, Figures 2c & 3c).

After applying the “WM/CSF/GS regression” procedure using the same statistical criteria as above, the corrected p ≤ 0.05 was achieved with a minimum cluster size of 344 voxels. Again, four clusters were identified as significant at brain regions consistent with the prior procedures. While the cluster at the right lingual gyrus, fusiform gyrus, and cuneus increases in size once more, clusters at regions involving precuneus and left fusiform gyrus show smaller cluster sizes (see Table 2, Figures 2d & 3d).

After applying the “aggressive AROMA” procedure using the same statistical criteria as above, the corrected p ≤ 0.05 was achieved with a minimum cluster size of 384 voxels. Two clusters spanning precuneus, cuneus, angular gyrus, fusiform gyrus, and lingual gyrus were identified as significant, with the most prominent cluster reaching a size of approximately 1072 mm3 (see Table 2, Figures 2e & 3e).

Overall, significant negative correlations were found between sALFF and AV45-SUVR with all four regression procedures, decreasing in cluster sizes in a stepwise manner based on the aggressiveness of physiological/vascular signal removal. Scatterplots, demonstrating the negative correlation between the mean sALFF and mean AV45-SUVR at the location of the most significant cluster, found using the “no physiological regression” procedure, are depicted in Figure 4, across the four regression procedures. The corresponding scatterplots for all significant voxel populations across all regression procedures are presented in Supplementary Figures S2-S5.

Figure 4
Figure 4 Correlation of brain amyloid-β deposition and the amplitude of low-frequency BOLD signal fluctuations (sALFF) for “No physiological regression” cluster #1 (Table 2).

Depicted are scatterplots for each physiological/vascular regression procedure, left to right: (a) “no physiological regression”, (b) “WM/CSF regression”, (c) “WM/CSF/GS regression”, and (e) “aggressive AROMA”. Blue crosses represent the 33 subjects, using mean sALFF and mean AV45 within this cluster. The plots show reduced correlation strength from (a) to (d). For scatterplots of the remaining significant clusters see Supplementary Figures 2-5.

To further understand the effect of the different physiological/vascular signal regression procedures on sALFF and AV45-SUVR correlations, the Fisher z transformed r-value distributions at significant voxel locations were illustrated in Figure 5. Using all significant voxel locations for either the “no physiological regression” procedure (4563 voxels) or the “aggressive AROMA” procedure (944 voxels), r-value distributions are depicted for each regression procedure at these specific voxel locations. Figure 5 (a) shows the violin plots of the r-value distributions of significant voxels from the no physiological/vascular regression procedure with mean and standard deviations as follows: - 0.47 ± 0.06 (no physiologic regression), -0.45 ± 0.07 (WM/CSF regression), -0.43 ± 0.08 (WM/CSF/GS regression) and -0.32 ± 0.13 (aggr. AROMA). Figure 5 (b) shows the opposite effect, depicting the violin plots with r-value distributions of the significant voxel population from the “aggressive AROMA” procedure. The r-value means and standard deviations at these voxel locations corresponding to each regression procedure are: -0.37 ± 0.14 (no physiologic regression), -0.36 ± 0.13 (WM/CSF regression), -0.38 ± 0.13 (WM/CSF/GS regression) and -0.46 ± 0.06 (aggr. AROMA).

Figure 5
Figure 5 Impact of stepwise regression of upper-stream vascular effects on the correlations between brain amyloid-β deposition and the amplitude of low-frequency BOLD signal fluctuations (sALFF) in the significant clusters identified with “no physiological regression” (a) and “aggressive AROMA” (b). Fisher z transformed r-value distributions shown in violin plots:

(a) Using the r-values of all voxels within the significant clusters reported in Table 2 under “no physiological regression”. This comparison shows a stepwise decrease of negative correlation magnitude as more aggressive physiological signal removal procedures were applied, suggesting that BOLD-fluctuations at these voxel locations are more driven by upper-stream vascular effects; (b) Using the r-values of all voxels within the significant clusters reported in Table 2 under “aggressive AROMA”. This comparison shows a stepwise increase of negative correlation magnitude as more aggressive physiological signal removal procedures were applied, suggesting that BOLD fluctuations at these voxel locations are more likely driven by neuronal activity.

Discussion

The main finding of this study is the observation of significant negative correlations between the sALFF of rs-fMRI BOLD signals and brain amyloid deposition measured with PET AV45-SUVR at regions associated with default mode and visual networks in patients with amnestic mild cognitive impairment. We have also shown that stepwise removal of the upper-stream physiological/vascular signals reduced the correlation strength between sALFF and AV45-SUVR at these regions. However, even the most intensive upper-stream physiological/vascular signal regression procedures implemented in this study did not completely abolish the negative correlations between sALFF and AV45-SUVR in regions associated with default-mode and visual networks. These findings suggest low-frequency upper-stream cerebral arterial pressure, as well as blood flow fluctuations, are transmitted into the cerebral microcirculation, where they may play an important role in brian Aβ deposition.33 Conversely, our data also support the hypothesis that local neuronal activity fluctuations contribute to Aβ homeostasis, independent of the upper-stream vascular effect. Below, we discuss the potential mechanisms and methodological considerations of this study.

Potential mechanisms

Brain metabolic waste products, including Aβ, are cleared partly via the recently identified brain lymphatic system.13,62 This system is a highly organized brain metabolite transportation system between brain interstitial fluid (ISF), cerebrospinal fluid (CSF), and the cerebral vasculature.13,62 One theory proposed that CSF enters the perivascular spaces from the subarachnoid space and is propelled deep into the brain by arterial pulsatility. The CSF subsequently enters the neuropil via the astrocytes AQP4 water channels. It then mixes with brain ISF in the extracellular space and leaves the brain, along with brain waste products via the perivenous space into the systemic lymphatic system.13 Conversely, other studies have shown that brain waste products also can be cleared via backflow through the periarterial basement space.62 Regardless of specific vascular clearance pathways, cerebral arterial pulsatility has been proposed as a key driver of brain waste clearance through these pathways.14 Furthermore, a recent study showed that spontaneous cerebral arteriolar vasomotion at a frequency of ∼0.1 Hz was correlated with paravascular clearance of Aβ in the awake mouse brain (van Veluw et al. Neuron 2020), while MRI studies also suggest that low-frequency waves in the CSF of ∼ 0.02 and ∼ 0.05 Hz are likely related to brain Aβ clearance.63,64 In this regard, our previous study has demonstrated that brain BOLD signal fluctuations at the regions associated with the default-mode and the visual networks have similar spectral distribution to the fluctuations of systemic blood pressure and CBF velocity measured from the middle cerebral artery, suggesting these upstream vascular signals may transmit downstream into the cerebral microcirculation and impact Aβ clearance.33

The present study extended these previous studies by showing that significant negative correlations exist between spontaneous low-frequency BOLD fluctuations and brain amyloid deposition, measured using PET AV45-SUVR, in default-mode and visual networks in patients with amnestic mild cognitive impairment. These findings taken together with the previous studies discussed above suggest that hemodynamic fluctuations transmitted from the upper-stream cardiovascular and cerebrovascular systems into the cerebral microcirculation may influence Aβ homeostasis in individuals with elevated risk for Alzheimer’s disease dementia. Consistent with this hypothesis, we have found that age-related carotid artery stiffness was associated positively with brain Aβ burden in patients with amnestic mild cognitive impairment, suggesting cerebral arterial stiffening may attenuate the transmission of arterial pressure and blood flow oscillations into the brain, leading to a reduction of brain amyloid clearance.11,65

Interestingly, even with the most intensive upper-stream physiological/vascular signal removal procedure implemented in this study, negative correlations between sALFF and AV45-SUVR were still present in default-mode and visual regions. These findings suggest that local neuronal activity fluctuations might also prevent brain Aβ deposition.

Methodological considerations

For calculating the PET AV45-SUVR, we chose white matter instead of the more traditional cerebellum as the reference for normalization. The comparison between reference regions has been comprehensively discussed by Brendel et al..60 In their study, they performed a discrimination analysis of healthy controls (HC) vs. mild cognitive impairment vs. Alzheimer’s disease of roughly 1000 subjects from the ADNI dataset. They compared the PET AV45-SUVR discriminative power using different reference regions, including the cerebellum, brainstem, and cortical white matter. Their receiver operating characteristics (ROC) analyses found that using white matter as the reference had the best performance in discriminating the diagnosis groups. Specifically, when compared to the more typical cerebellar reference scaling, a higher discriminatory power between healthy controls and Alzheimer’s disease was found when the white matter was used as a reference. Inter-subject variability was lowest when the white matter was used as a reference as well. Assessment of the longitudinal amyloid deposition was also more reliable and consistent when using white matter as the reference. Using white matter as the reference for the PET AV45-SUVR calculation, Brendel et al. showed that additional partial volume effect correction could lead to 0-2% higher scores in sensitivity and specificity, thus reaching the highest discriminative power of all comparisons. Since the effect of partial volume correction appeared small and its benefits have not been consistent,66 we analyzed the AV45 data without partial volume correction. We repeated our analyses with the cerebellum as the reference for AV45-SUVR calculation, finding the same trend of negative correlation between local sALFF and AV45-SUVR (data not shown).

The mean images of the well-established standardized amplitude of low-frequency fluctuations (sALFF) of fMRI recordings in Figures 2 and 3 show a strong effect of preprocessing strategy. As fMRI data is generally suffering from a multitude of confounds and noise signals, the goal of advanced preprocessing strategies is to achieve cleaned BOLD response signals stemming only from neuronal activity. In the present study, we employed different noise regression models, to remove physiological/vascular signals. All regression strategies had inherent motion confound modeling, yet the amount of physiological/vascular regression varied. While white matter (WM) and cerebrospinal fluid (CSF) signals are widely accepted to represent physiological/vascular signal components,53 the much-debated global signal has recently been shown to also represent a substantial amount of physiological/vascular signals.57,67 Using the data-driven “ICA-AROMA” approach warrants noise-regressors that are closely tailored to the data and promise an effective and robust noise regression,58 making it the most demanding regression model in this study.

Limitations and strengths

Only 33 participants were available in this study and thus the statistical power for voxel-wise analyses was limited. However, our voxel-wise analyses have clearly shown a correlation between sALFF and PET AV45 SUVR. Changes in the correlation strength due to different regression procedures were further presented with cluster-based correlation analyses. However, a causal relationship between brain Aβ deposition and BOLD fluctuation quantified by ALFF cannot be drawn from this cross-sectional study.

The discriminative power of PET AV45 imaging in Alzheimer’s disease has been proven in many studies, nonetheless, there are limitations.68 PET images inherently contain a high level of spatial blurring, in addition to a relatively low spatial resolution provided by signal detection. Additionally, the high level of non-specific binding to white matter leads to segmentation and registration inaccuracies at the border between gray- and white matter.69 Especially in the situation of grey matter atrophy, commonly observed in Alzheimer’s disease brains, the signal at the cortical grey matter is inevitably contaminated by white matter signal, resulting in reduced specificity. Using the advanced freesurfer segmentation from the anatomical MRI recordings, paired with the manually supervised alignment of AV45 PET images to the anatomical images, we can reduce these segmentation issues in AV45 PET images.

Conclusions

This study has shown that local spontaneous BOLD signal fluctuations are associated negatively with brain Aβ accumulation measured with amyloid PET in individuals with amnestic mild cognitive impairment who are at a higher risk for Alzheimer’s disease. Furthermore, we found that physiological signal regression decreased the cluster size of negative correlations between low-frequency BOLD signal fluctuation and PET AV45 SUVR. Nevertheless, the most intensive physiological/vascular regression procedures did not completely abolish the correlations between BOLD signal fluctuations and PET AV45 SUVR. Therefore, these findings collectively suggest that both cerebral hemodynamic fluctuations and local neuronal activity play important roles in brain Aβ homeostasis in older adults with amnestic mild cognitive impairment.

Data Availability

The raw data were generated at the University of Texas Southwestern Medical Center and UTSW and Institute for Exercise and Environmental Medicine Texas Health Presbyterian Hospital, Dallas, TX, USA. Derived data supporting the findings of this study are available from the corresponding authors upon reasonable request.

Funding

Funding towards this work was received from the National Institutes of Health (R01AG033106, R01HL102457, and R01AG057571).

Competing interests

The authors report no competing interests.

Supplementary material

Supplementary material is available at medRxiv online.

Acknowledgements

We thank each of the study participants for their effort and time contributing to the study. The 18F-florbetapir PET radiotracer was provided to the study by Avid Radiopharmaceuticals.

Abbreviations

ADNI
Alzheimer’s Disease Neuroimaging Initiative
AFNI
Analysis of Functional NeuroImages
ALFF
Amplitude of Low-Frequency Fluctuations
APOE4
Apolipoprotein E4
AROMA
ICA-based Automatic Removal Of Motion Artifacts
AV45
18F-Florbetapir
Aβ
Amyloid-β
BOLD
Blood-Oxygen-Level-Dependent
CBF
Cerebral Blood Flow
CDR
Clinical Dementia Rating
DPABI
Toolbox for Data Processing & Analysis of Brain Imaging
DPARSFA
Data Processing Assistant for Resting-State fMRI advanced edition
EPI
Echo Planar Imaging
FFT
Fast Fourier Transformation
flirt
FMRIB’s Linear Image Registration Tool
fMRI
functional MRI
FMRIB
Functional Magnetic Resonance Imaging of the Brain
FSL
FMRIB Software Library
FWHM
Full Width at Half Maximum
GS
Global Signal
ICA
Independent Component Analysis
ISF
Interstitial Fluid
MELODIC
Multivariate Exploratory Linear Optimized Decomposition into Independent Components
MMSE
Mini-Mental State Examination
MNI
Montreal Neurological Institute
MPRAGE
Magnetization-Prepared Rapid Acquisition Gradient-Echo
MRI
Magnetic Resonance Imaging
ROC
Receiver Operating Characteristics
rs-fMRI
resting-state fMRI
sALFF
standardized ALFF
SUVR
Standardized Uptake Value Ratio
TE
Time of Echo
TR
Time of Repetition
WM
White Matter
WMS-LM
Wechsler Memory Scale – Revised Logical Memory

References

  1. 1.↵
    Matthews KA, Xu W, Gaglioti AH, et al. Racial and ethnic estimates of Alzheimer’s disease and related dementias in the United States (2015–2060) in adults aged ≥65 years. Alzheimer’s Dement. 2019;15(1):17–24. doi:10.1016/j.jalz.2018.06.3063
    OpenUrlCrossRefPubMed
  2. 2.↵
    McKhann GM, Knopman DS, Chertkow H, 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. Alzheimer’s Dement. 2011;7(3):263–269. doi:10.1016/j.jalz.2011.03.005
    OpenUrlCrossRefPubMedWeb of Science
  3. 3.↵
    Querfurth HW, LaFerla FM. Alzheimer’s disease. N Engl J Med. 2010;362(4):329– 344. doi:10.1056/NEJMra0909142
    OpenUrlCrossRefPubMedWeb of Science
  4. 4.↵
    Petersen RC, Lopez O, Armstrong MJ, et al. Practice guideline update summary: Mild cognitive impairment report of theguideline development, dissemination, and implementation. Neurology. 2018;90(3):126–135. doi:10.1212/WNL.0000000000004826
    OpenUrlCrossRefPubMed
  5. 5.↵
    Hardy JA, Higgins GA. Alzheimer’s disease: the amyloid cascade hypothesis. Science. 1992;256(5054):184–185. doi:10.1126/science.1566067
    OpenUrlFREE Full Text
  6. 6.↵
    De la Torre JC, Mussivand T. Can disturbed brain microcirculation cause Alzheimer’s disease? Neurol Res. 1993;15(3):146–153. doi:10.1080/01616412.1993.11740127
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.
    de la Torre J. The Vascular Hypothesis of Alzheimer’s Disease: A Key to Preclinical Prediction of Dementia Using Neuroimaging. J Alzheimers Dis. 2018;63(1):35–52. doi:10.3233/JAD-180004
    OpenUrlCrossRefPubMed
  8. 8.↵
    Iadecola C, Duering M, Hachinski V, et al. Vascular Cognitive Impairment and Dementia: JACC Scientific Expert Panel. J Am Coll Cardiol. 2019;73(25):3326–3344. doi:10.1016/j.jacc.2019.04.034
    OpenUrlFREE Full Text
  9. 9.↵
    Gaugler J, James B, Johnson T, Scholz K, Weuve J. 2016 Alzheimer’s disease facts and figures. Alzheimer’s Dement. 2016;12(4):459–509. doi:10.1016/j.jalz.2016.03.001
    OpenUrlCrossRefPubMed
  10. 10.↵
    Jack CR, Bennett DA, Blennow K, et al. NIA-AA Research Framework: Toward a biological definition of Alzheimer’s disease. Alzheimer’s Dement. 2018;14(4):535– 562. doi:10.1016/j.jalz.2018.02.018
    OpenUrlCrossRefPubMed
  11. 11.↵
    Pasha EP, Rutjes E, Tomoto T, et al. Carotid Stiffness is Associated with Brain Amyloid-β Burden in Amnestic Mild Cognitive Impairment. J Alzheimer’s Dis. 2020;74(3):925–935. doi:10.3233/JAD-191073
    OpenUrlCrossRef
  12. 12.↵
    van Veluw SJ, Hou SS, Calvo-Rodriguez M, et al. Vasomotion as a Driving Force for Paravascular Clearance in the Awake Mouse Brain. Neuron. 2020;105(3):549-561.e5. doi:10.1016/j.neuron.2019.10.033
    OpenUrlCrossRef
  13. 13.↵
    Nedergaard M, Goldman SA. Glymphatic failure as a final common pathway to dementia. Science (80-). 2020;370(6512):50–56. doi:10.1126/science.abb8739
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    Iliff JJ, Wang M, Zeppenfeld DM, et al. Cerebral arterial pulsation drives paravascular CSF-Interstitial fluid exchange in the murine brain. J Neurosci. 2013;33(46):18190–18199. doi:10.1523/JNEUROSCI.1592-13.2013
    OpenUrlAbstract/FREE Full Text
  15. 15.↵
    Biswal B, Yetkin FZ, Haughton VM, Hyde JS. Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn Reson Med. 1995;34(4):537–541. http://www.ncbi.nlm.nih.gov/pubmed/8524021
    OpenUrlCrossRefPubMedWeb of Science
  16. 16.↵
    Fox MD, Raichle ME. Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nat Rev Neurosci. Published online 2007. doi:10.1038/nrn2201
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    Heeger DJ, Ress D. What does fMRI tell us about neuronal activity? Nat Rev Neurosci. 2002;3(2):142–151. doi:10.1038/nrn730
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    Greicius MD, Krasnow B, Reiss AL, Menon V. Functional connectivity in the resting brain: A network analysis of the default mode hypothesis. Proc Natl Acad Sci U S A. 2003;100(1):253–258. doi:10.1073/pnas.0135058100
    OpenUrlAbstract/FREE Full Text
  19. 19.
    Greicius MD, Supekar K, Menon V, Dougherty RF. Resting-state functional connectivity reflects structural connectivity in the default mode network. Cereb Cortex. 2009;19(1):72–78. doi:10.1093/cercor/bhn059
    OpenUrlCrossRefPubMedWeb of Science
  20. 20.
    Johnston JM, Vaishnavi SN, Smyth MD, et al. Loss of resting interhemispheric functional connectivity after complete section of the corpus callosum. J Neurosci. 2008;28(25):6453–6458. doi:10.1523/JNEUROSCI.0573-08.2008
    OpenUrlAbstract/FREE Full Text
  21. 21.
    Ebisch SJH, Gallese V, Willems RM, et al. Altered intrinsic functional connectivity of anterior and posterior insula regions in high-functioning participants with autism spectrum disorder. Hum Brain Mapp. 2011;32(7):1013–1028. doi:10.1002/hbm.21085
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.
    Fair DA, Posner J, Nagel BJ, et al. Atypical default network connectivity in youth with attention-deficit/ hyperactivity disorder. Biol Psychiatry. 2010;68(12):1084–1091. doi:10.1016/j.biopsych.2010.07.003
    OpenUrlCrossRefPubMedWeb of Science
  23. 23.
    Baudrexel S, Witte T, Seifried C, et al. Resting state fMRI reveals increased subthalamic nucleus-motor cortex connectivity in Parkinson’s disease. Neuroimage. 2011;55(4):1728–1738. doi:10.1016/j.neuroimage.2011.01.017
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.
    Chen G, Ward BD, Xie C, et al. Classification of AD, MCI, and NC Status with Large- Scale Network Analysis Based on rs-fMRI. Radiology. 2011;259(1):213–221. doi:10.1148/radiol.10100734/-/DC1
    OpenUrlCrossRefPubMedWeb of Science
  25. 25.
    Park CH, Chang WH, Ohn SH, et al. Longitudinal changes of resting-state functional connectivity during motor recovery after stroke. Stroke. 2011;42(5):1357–1362. doi:10.1161/STROKEAHA.110.596155
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    Zhu DC, Majumdar S, Korolev IO, Berger KL, Bozoki AC. Alzheimer’s disease and amnestic mild cognitive impairment weaken connections within the default-mode network: A multi-modal imaging study. J Alzheimer’s Dis. 2013;34(4):969–984. doi:10.3233/JAD-121879
    OpenUrlCrossRefPubMedWeb of Science
  27. 27.↵
    Raichle ME, Mintun MA. Brain Work and Brain Imaging. Annu Rev Neurosci. 2006;29(1):449–476. doi:10.1146/annurev.neuro.29.051605.112819
    OpenUrlCrossRefPubMedWeb of Science
  28. 28.↵
    Shmueli K, van Gelderen P, de Zwart J a, et al. Low-frequency fluctuations in the cardiac rate as a source of variance in the resting-state fMRI BOLD signal. Neuroimage. 2007;38(2):306–320. doi:10.1016/j.neuroimage.2007.07.037
    OpenUrlCrossRefPubMedWeb of Science
  29. 29.
    Chang C, Cunningham JP, Glover GH. Influence of heart rate on the BOLD signal: the cardiac response function. Neuroimage. 2009;44(3):857–869. doi:10.1016/j.neuroimage.2008.09.029
    OpenUrlCrossRefPubMedWeb of Science
  30. 30.
    Chang C, Glover GH. Time-frequency dynamics of resting-state brain connectivity measured with fMRI. Neuroimage. 2010;50(1):81–98. doi:10.1016/j.neuroimage.2009.12.011
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.
    Chang C, Metzger CD, Glover GH, Duyn JH, Heinze H-J, Walter M. Association between heart rate variability and fluctuations in resting-state functional connectivity. Neuroimage. 2013;68:93–104. doi:10.1016/j.neuroimage.2012.11.038
    OpenUrlCrossRefPubMedWeb of Science
  32. 32.↵
    Tsvetanov KA, Henson RNA, Jones PS, et al. The effects of age on resting-state BOLD signal variability is explained by cardiovascular and cerebrovascular factors. Psychophysiology. 2020;(July):1–20. doi:10.1111/psyp.13714
    OpenUrlCrossRef
  33. 33.↵
    Zhu DC, Tarumi T, Khan MA, Zhang R. Vascular coupling in resting-state fMRI: Evidence from multiple modalities. J Cereb Blood Flow Metab. 2015;35(12):1910– 1920. doi:10.1038/jcbfm.2015.166
    OpenUrlCrossRefPubMed
  34. 34.↵
    Tuovinen T, Kananen J, Rajna Z, Lieslehto J, Korhonen V. The variability of functional MRI brain signal increases in Alzheimer ’ s disease at cardiorespiratory frequencies. Sci Rep. Published online 2020:1–11. doi:10.1038/s41598-020-77984-1
    OpenUrlCrossRef
  35. 35.↵
    Skidmore FM, Yang M, Baxter L, et al. Apathy, depression, and motor symptoms have distinct and separable resting activity patterns in idiopathic Parkinson disease. Neuroimage. Published online July 14, 2011. doi:10.1016/j.neuroimage.2011.07.012
    OpenUrlCrossRefPubMed
  36. 36.
    Skidmore FM, Yang M, Baxter L, et al. Reliability analysis of the resting state can sensitively and specifically identify the presence of Parkinson disease. Neuroimage. Published online September 6, 2011. doi:10.1016/j.neuroimage.2011.06.056
    OpenUrlCrossRefPubMed
  37. 37.
    Sato JR, Hoexter MQ, Fujita A, Rohde LA. Evaluation of pattern recognition and feature extraction methods in ADHD prediction. Front Syst Neurosci. 2012;6(September):68. doi:10.3389/fnsys.2012.00068
    OpenUrlCrossRefPubMed
  38. 38.↵
    Zang YF, Yong H, Chao-Zhe Z, et al. Altered baseline brain activity in children with ADHD revealed by resting-state functional MRI. Brain Dev. 2007;29(2):83–91. doi:10.1016/j.braindev.2006.07.002
    OpenUrlCrossRefPubMedWeb of Science
  39. 39.↵
    Pruim RHR, Mennes M, van Rooij D, Llera Arenas A, Buitelaar JK, Beckmann CF. ICA-AROMA: A robust ICA-based strategy for removing motion artifact from fMRI data. Neuroimage. Published online 2015. doi:10.1016/j.neuroimage.2015.02.064
    OpenUrlCrossRefPubMed
  40. 40.
    Caballero-Gaudes C, Reynolds RC. Methods for cleaning the BOLD fMRI signal. Neuroimage. 2016;(December):0–1. doi:10.1016/j.neuroimage.2016.12.018
    OpenUrlCrossRefPubMed
  41. 41.↵
    Poline J-B, Brett M. The general linear model and fMRI: Does love last forever? Neuroimage. 2012;62(2):871–880. doi:10.1016/j.neuroimage.2012.01.133
    OpenUrlCrossRefPubMed
  42. 42.↵
    Tarumi T, Rossetti H, Thomas BP, et al. Exercise Training in Amnestic Mild Cognitive Impairment: A One-Year Randomized Controlled Trial. J Alzheimer’s Dis. 2019;71(2):421–433. doi:10.3233/JAD-181175
    OpenUrlCrossRef
  43. 43.↵
    Petersen RC, Smith GE, Waring SC, Ivnik RJ, Tangalos EG, Kokmen E. Mild cognitive impairment: clinical characterization and outcome. Arch Neurol. 1999;56(3):303–308. doi:http://dx.doi.org/10.1212/01.CON.0000429175.29601.97
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    Morris JC. The clinical dementia rating (CDR): Current version and scoring rules. Neurology. 1993;43(11):2412–2414. doi:10.1212/wnl.43.11.2412-a
    OpenUrlFREE Full Text
  45. 45.↵
    Folstein MF, Folstein SE, McHugh PR. “Mini-mental state”: a practical method for grading the cognitive state of patients for the clinician. J psychiat Res. 1975;12(3):189– 198.
    OpenUrlCrossRefPubMedWeb of Science
  46. 46.↵
    Wechsler D. Wechsler Memory Scale-Revised: Manual. San Antonio, TX: Psychological Corporation 1987; 1987.
  47. 47.↵
    Clark CM, Schneider JA, Bedell BJ, et al. Use of florbetapir-PET for imaging β- amyloid pathology. JAMA - J Am Med Assoc. 2011;305(3):275–283. doi:10.1001/jama.2010.2008
    OpenUrlCrossRefPubMedWeb of Science
  48. 48.↵
    Rodrigue KM, Kennedy KM, Devous MD, et al. β-amyloid burden in healthy aging: Regional distribution and cognitive consequences. Neurology. 2012;78(6):387–395. doi:10.1212/WNL.0b013e318245d295
    OpenUrlAbstract/FREE Full Text
  49. 49.↵
    Chao-Gan Y, Yu-Feng Z. DPARSF: A MATLAB Toolbox for Pipeline Data Analysis of Resting-State fMRI. Front Syst Neurosci. 2010;4(May):13. doi:10.3389/fnsys.2010.00013
    OpenUrlCrossRefPubMed
  50. 50.↵
    Yan CG, Wang X Di, Zuo XN, Zang YF. DPABI: Data Processing & Analysis for (Resting-State) Brain Imaging. Neuroinformatics. 2016;14(3):339–351. doi:10.1007/s12021-016-9299-4
    OpenUrlCrossRefPubMed
  51. 51.↵
    Fischl B, Salat DH, Busa E, et al. Whole brain segmentation: Automated labeling of neuroanatomical structures in the human brain. Neuron. 2002;33(3):341–355. doi:10.1016/S0896-6273(02)00569-X
    OpenUrlCrossRefPubMedWeb of Science
  52. 52.↵
  53. 53.↵
    Fox MD, Snyder AZ, Vincent JL, Corbetta M, Essen DC Van, Raichle ME. The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc Natl Acad Sci U S A. 2005;102(27):9673–9678. doi:10.1073/pnas.0504136102
    OpenUrlAbstract/FREE Full Text
  54. 54.↵
    Xifra-Porxas A, Kassinopoulos M, Mitsis GD. Physiological and head motion signatures in static and time-varying functional connectivity and their subject discriminability. bioRxiv. Published online 2020:2020.02.04.934554. doi:10.1101/2020.02.04.934554
    OpenUrlAbstract/FREE Full Text
  55. 55.
    Kassinopoulos M, Mitsis GD. Identification of physiological response functions to correct for fluctuations in resting-state fMRI related to heart rate and respiration. Neuroimage. 2019;202(August):116150. doi:10.1016/j.neuroimage.2019.116150
    OpenUrlCrossRefPubMed
  56. 56.↵
    Kassinopoulos M, Mitsis GD. White Matter Denoising Improves the Identifiability of Large-Scale Networks and Reduces the Effects of Motion in fMRI Functional Connectivity. bioRxiv. Published online 2020:837609. doi:10.1101/837609
    OpenUrlAbstract/FREE Full Text
  57. 57.↵
    Xifra-Porxas A, Kassinopoulos M, Prokopiou P, Boudrias M-H, Mitsis GD. Does global signal regression remove alpha power fluctuations? An EEG-fMRI study in humans at rest. In: 26TH ANNUAL MEETING OF THE ORGANIZATION FOR HUMAN BRAIN MAPPING. ; 2020:1498.
  58. 58.↵
    Pruim RHR, Mennes M, Buitelaar JK, Beckmann CF. Evaluation of ICA-AROMA and alternative strategies for motion artifact removal in resting-state fMRI. Neuroimage. Published online 2015. doi:10.1016/j.neuroimage.2015.02.063
    OpenUrlCrossRefPubMed
  59. 59.↵
    Jenkinson M, Beckmann CF, Behrens TEJ, Woolrich MW, Smith SM. FSL. Neuroimage. 2012;62:782–790. doi:10.1016/j.neuroimage.2011.09.015
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    Brendel M, Högenauer M, Delker A, et al. Improved longitudinal [18F]-AV45 amyloid PET by white matter reference and VOI-based partial volume effect correction. Neuroimage. 2015;108:450–459. doi:10.1016/j.neuroimage.2014.11.055
    OpenUrlCrossRefPubMed
  61. 61.↵
    Golland P, Liang F, Mukherjee S, Panchenko D. Permutation tests for classification. Learn Theory. 2005;1:1-22. Accessed February 27, 2013. http://www.springerlink.com/index/WT2P7419VNDKX4HL.pdf
    OpenUrl
  62. 62.↵
    Weller RO, Hawkes CA, Kalaria RN, Werring DJ, Carare RO. White matter changes in dementia: Role of impaired drainage of interstitial fluid. Brain Pathol. 2015;25(1):63–78. doi:10.1111/bpa.12218
    OpenUrlCrossRef
  63. 63.↵
    Kiviniemi V, Wang X, Korhonen V, et al. Ultra-fast magnetic resonance encephalography of physiological brain activity-Glymphatic pulsation mechanisms? J Cereb Blood Flow Metab. 2016;36(6):1033–1045. doi:10.1177/0271678X15622047
    OpenUrlCrossRefPubMed
  64. 64.↵
    Fultz NE, Bonmassar G, Setsompop K, et al. Coupled electrophysiological, hemodynamic, and cerebrospinal fluid oscillations in human sleep. Science (80-). 2019;366(6465):628–631. doi:10.1126/science.aax5440
    OpenUrlAbstract/FREE Full Text
  65. 65.↵
    Sugawara J, Tarumi T, Xing C, et al. Older Age and Male Sex are Associated with Higher Cerebrovascular Impedance. J Appl Physiol. 2020;(214). doi:10.1152/japplphysiol.00396.2020
    OpenUrlCrossRef
  66. 66.↵
    Schwarz CG, Gunter JL, Lowe VJ, et al. A Comparison of Partial Volume Correction Techniques for Measuring Change in Serial Amyloid PET SUVR. J Alzheimer’s Dis. 2019;67(1):181–195. doi:10.3233/JAD-180749
    OpenUrlCrossRef
  67. 67.↵
    Murphy K, Birn RM, Handwerker D a, Jones TB, Bandettini P a. The impact of global signal regression on resting state correlations: are anti-correlated networks introduced? Neuroimage. 2009;44(3):893–905. doi:10.1016/j.neuroimage.2008.09.036
    OpenUrlCrossRefPubMedWeb of Science
  68. 68.↵
    Saint-Aubert L, Barbeau EJ, Péran P, et al. Cortical florbetapir-PET amyloid load in prodromal Alzheimer’s disease patients. EJNMMI Res. 2013;3(1):1–22. doi:10.1186/2191-219X-3-43
    OpenUrlCrossRefPubMed
  69. 69.↵
    Wolk DA, Zhang Z, Boudhar S, Clark CM, Pontecorvo MJ, Arnold SE. Amyloid imaging in Alzheimer’s disease: Comparison of florbetapir and Pittsburgh compound- B positron emission tomography. J Neurol Neurosurg Psychiatry. 2012;83(9):923–926. doi:10.1136/jnnp-2012-302548
    OpenUrlAbstract/FREE Full Text
Back to top
Posted April 26, 2021.
Download PDF

Supplementary Material

Data/Code
Email
Resting-state functional MRI signal fluctuations are correlated with brain amyloid-β deposition
Norman Scheel, Takashi Tarumi, Tsubasa Tomoto, C. Munro Cullum, Rong Zhang, David C. Zhu
medRxiv 2021.04.22.21255924; doi: https://doi.org/10.1101/2021.04.22.21255924
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area