- 1First Central Clinical College, Tianjin Medical University, Tianjin, China
- 2Department of Radiology and Tianjin Key Laboratory of Functional Imaging, Tianjin Medical University General Hospital, Tianjin, China
- 3Department of Radiology, Affiliated Hospital of Hebei University, Baoding, China
- 4Tianjin International Joint Research Center for Neural Engineering, Academy of Medical Engineering and Translational Medicine, Tianjin University, Tianjin, China
- 5MR Collaboration, Siemens Healthcare Ltd., Beijing, China
- 6Department of Radiology, Tianjin First Central Hospital, Tianjin, China
Objective: To investigate the relationship between changes in cerebral blood flow (CBF) and gray matter (GM) microstructure in Alzheimer’s disease (AD) and mild cognitive impairment (MCI).
Methods: A recruited cohort of 23 AD patients, 40 MCI patients, and 37 normal controls (NCs) underwent diffusional kurtosis imaging (DKI) for microstructure evaluation and pseudo-continuous arterial spin labeling (pCASL) for CBF assessment. We investigated the differences in diffusion- and perfusion-related parameters across the three groups, including CBF, mean diffusivity (MD), mean kurtosis (MK), and fractional anisotropy (FA). These quantitative parameters were compared using volume-based analyses for the deep GM and surface-based analyses for the cortical GM. The correlation between CBF, diffusion parameters, and cognitive scores was assessed using Spearman coefficients, respectively. The diagnostic performance of different parameters was investigated with k-nearest neighbor (KNN) analysis, using fivefold cross-validation to generate the mean accuracy (mAcc), mean precision (mPre), and mean area under the curve (mAuc).
Results: In the cortical GM, CBF reduction primarily occurred in the parietal and temporal lobes. Microstructural abnormalities were predominantly noted in the parietal, temporal, and frontal lobes. In the deep GM, more regions showed DKI and CBF parametric changes at the MCI stage. MD showed most of the significant abnormalities among all the DKI metrics. The MD, FA, MK, and CBF values of many GM regions were significantly correlated with cognitive scores. In the whole sample, the MD, FA, and MK were associated with CBF in most evaluated regions, with lower CBF values associated with higher MD, lower FA, or lower MK values in the left occipital lobe, left frontal lobe, and right parietal lobe. CBF values performed best (mAuc = 0.876) for distinguishing the MCI from the NC group. Last, MD values performed best (mAuc = 0.939) for distinguishing the AD from the NC group.
Conclusion: Gray matter microstructure and CBF are closely related in AD. Increased MD, decreased FA, and MK are accompanied by decreased blood perfusion throughout the AD course. Furthermore, CBF values are valuable for the predictive diagnosis of MCI and AD. GM microstructural changes are promising as novel neuroimaging biomarkers of AD.
Introduction
Alzheimer’s disease (AD) is a neurodegenerative disorder characterized by extracellular accumulation of amyloid β protein (Aβ) plaques and intracellular neurofibrillary tau tangles (Guo et al., 2020). Studies within the last few decades provide growing evidence that cerebral vascular dysfunction in various molecular pathways plays a central role in AD pathogenesis. Due to blood–brain barrier (BBB) breakdown, the accumulation of excess neurotoxic Aβ occurs, resulting in impaired BBB function and cerebral hypoperfusion (Chakraborty et al., 2017; Yamazaki and Kanekiyo, 2017). In addition, restricted Aβ clearance further damages surrounding neurons and neurites, leading to microstructural gray matter (GM) loss, as in the hippocampus with associated cognitive changes (Wang et al., 2018). Given that brain microstructural and cerebral blood flow (CBF) changes occur years before the onset of clinical symptoms, efforts to detect these alterations early could significantly enhance our ability to diagnose AD sooner.
Diffusion kurtosis imaging (DKI) is an MRI technique that extends conventional diffusion tensor imaging (DTI) by estimating the kurtosis of the water diffusion probability distribution function, allowing for subtle pathology quantification in gray and white matter architecture (Jensen et al., 2005). DKI provides kurtosis metrics such as mean kurtosis (MK) related to the non-Gaussianity of water diffusion for reflecting microstructure complexity. In addition, diffusion parameters such as fractional anisotropy (FA) and mean diffusivity (MD) are also acquired from DKI, symbols of neuronal status, including axonal, dendritic, and synaptic integrity. Although DKI applications in AD and mild cognitive impairment (MCI) have been largely performed to probe alterations in white matter, they have great potential for detecting GM microstructural changes.
One study showed that microstructural abnormalities in cortical GM occurred primarily in the parietal and frontal lobes in AD (Gong et al., 2017). Moreover, some research suggested that cortical MD was sensitive enough to predict progression to MCI (Wang et al., 2018; Rodriguez-Vieitez et al., 2021). Arterial spin labeling (ASL) is a non-invasive technique that labels arterial blood water as an endogenous contrast tracer to measure CBF at the tissue level (Hernandez-Garcia et al., 2019). The diagnostic performance detected by ASL was similar to that identified by single photon emission computed tomography (SPECT) (Johnson et al., 2005). Hence, ASL has become a fundamental CBF imaging technique in brain research. In addition, some studies suggested that CBF can serve as a reliable hallmark in the early diagnosis of AD and MCI (Hays et al., 2016; Korte et al., 2020; Camargo and Wang, 2021).
With combined cerebral microstructure and CBF, Alisch et al. (2021) found that lower microstructural integrity was associated with lower CBF in the choroid plexus. This founding suggested that adequate CBF is paramount for nutrient and oxygen delivery and the maintenance of tissue integrity. Moreover, Bouhrara et al. (2022) found that FA and MD were associated with CBF in most regions evaluated, including the frontal and parietal white matter lobes in healthy adults. All these findings indicated that the CBF changes might affect cerebral tissue health. However, Bouhrara’s study only focused on cohorts of cognitively unimpaired subjects, and their results about association between CBF and GM microstructure were mostly not significant in most regions, likely due to the limitation of DTI to capture changes in GM as compared to DKI technique. Furthermore, the association between GM microstructural integrity and perfusion in individuals with cognitive impairment has currently received little consideration. So it is necessary that we combined DKI and pseudo-continuous arterial spin labeling (pCASL) to investigate further associations established in cognitive impairments.
Given the impairment of neurovascular coupling in AD pathogenesis, there may have a potential link between them. Therefore, the specific aims of this study were to analyze abnormalities of the GM microstructure and CBF and determine the association of these changes with cognitive function, to investigate the relationship between microstructural and perfusion parameters in combined DKI and ASL assessment in early AD patients, and to explore the diagnostic performance of MD, FA, MK, and CBF values in AD and MCI patients.
Materials and methods
Subjects
The study was approved by the hospital ethics committee, and all subjects signed informed consent. Twenty-three AD patients, 40 MCI patients, and 37 healthy age- and sex-matched controls were recruited in our hospital from September 2021 to December 2022. All subjects underwent MRI examinations and cognitive evaluations with psychological tests such as the mini-mental state examination (MMSE), Montreal cognitive assessment (MoCA), Clinical Dementia Rating (CDR) Scale, activities of daily living (ADL) assessment, Hachinski Ischemic Scale (HIS), and the Hamilton Depression Rating Scale (HAMD).
Patients met the criteria for AD based on the National Institute on Aging and Alzheimer’s Association (NIA-AA) (McKhann et al., 2011) and for MCI using Petersen’s criteria (Petersen et al., 2018). The exclusion criteria were as follows: (1) patients with brain diseases other than AD or MCI; (2) severe cerebrovascular diseases and HIS ≥4; (3) patients with a psychiatric disorder or severe head trauma history; (4) poor image quality in post-processing; or (5) left-handedness.
MRI acquisition
All the participants underwent MR scans on a 3T MAGNETOM Prisma MR scanner (Siemens Healthcare, Erlangen, Germany) equipped with a 64-channel head-neck coil. In addition, 3D-pCASL images were acquired using the following parameters: TR/TE = 4,100/36.68 ms, FOV = 240 mm × 240 mm, matrix = 96 × 96, slice thickness = 3 mm, multiple post-labeling delay time = 500/1,000/1,500/2,000/2,500 ms, flip angle = 120°, and 12 pairs of label and control images with a scan time of 7 min and 11 s. Furthermore, diffusion-weighted images were obtained with a spin-echo echo-planar imaging (SE-EPI) sequence. The DKI parameters were TR/TE = 3,800/72 ms, FOV = 220 mm × 220 mm, voxel = 2.0 mm × 2.0 mm × 2.2 mm, matrix = 110 × 110, slices = 60, simultaneous multi-slice acquisition (SMS) acceleration factor = 2, b = 0/2,000/,3000 s/mm2, and directions = 64. The parameters for 3D T1-weighted imaging acquisition were TR/TE = 1,550/2.98 ms, flip angle = 9°, FOV = 256 mm × 256 mm, voxel = 1.0 mm × 1.0 mm × 1.0 mm, and slices = 176.
Image analysis
ASL image processing
Raw ASL data were analyzed with CereFlow software programs.1 Motion correction was performed on all ASL images. Pairwise subtraction was used between the label and control images, followed by averaging to generate the mean image. Quantitative CBF maps were calculated using a previously reported model (Hernandez-Garcia et al., 2019). Then, registration to the AAL atlas was performed on CBF maps to calculate the CBF values of every ROI.
DKI image processing
First, all raw diffusion images were converted to the Neuroimaging Informatics Technology Initiative (NIFTI) format. Then, a rigid body transformation between the b = 0 image and all the diffusion-weighted acquisition was applied to mitigate motion effects. After correcting for movement, non-brain tissue was removed using the Brain Extraction Tool with the FSL toolbox (Smith et al., 2004) (Oxford Centre for Functional MRI of the brain).2 Finally, Diffusional Kurtosis Estimator (DKE) software (Tabesh et al., 2011)3 was used to calculate the MK, FA, and MD maps by loading 4D NIFTI diffusion images.
FreeSurfer processing
The segmentation of T1-weighted images was performed with the FreeSurfer package4 using the recon-all pipeline that has been described in a previous publication (Montal et al., 2018). This processing comprises motion correction and a skull strip, talairach transform computation, automated segmentation of subcortical white matter, and deep GM structures.
We used a surface-based approach to process cortical diffusion and perfusion images because the commonly used voxel-based morphometry approach has defects in GM analyses, where partial volume effects may make cortical parameter measurements that are biased due to the contribution from the CSF and white matter signal on the GM voxels. The detailed surface-based DKI steps were as follows: (1) the skull-stripped b0 diffusion images were co-registered to the segmented structural T1-weighted volume using the bbregister command, a boundary-based registration algorithm in FreeSurfer (Greve and Fischl, 2009). All registrations were visually inspected to detect and correct processing errors in the registration; (2) DKI volumes, including MD, FA, and MK, were projected to each individual’s surface space generated during the cortical segmentation. At each vertex, cortical DKI maps were sampled at the midpoint along the normal vector between white and pial surfaces using the mri_vol2surf command (Greve et al., 2014); and (3) All maps were normalized to a standard surface template (fsaverage) and smoothed with a 15-mm full-width half-maximum (FWHM) Gaussian kernel.
The process of surface-based perfusion images was similar to the described approach. First, a boundary-based registration of the CBF and M0 volumes to the T1-weighted volume was performed, followed by projection on an averaged surface. Then, all maps were smoothed with a 15-mm FWHM Gaussian kernel. In addition, cortical Desikan-Killiany atlases were used to extract DKI and CBF GM parameters.
All the maps analysis were applied a partial volume effect correction in order to improve the accuracy of measurements. After motion correction, maps estimation and registration, the contribution of the CSF in each voxel, as segmented by FreeSurfer’s gtmseg, was computed to correct the diffusion – and perfusion – maps in the cortical GM (Greve et al., 2016). Then, mri_gtmpvc command was used to carry out partial volume correction. The maps without contamination of CSF were performed further process.
Statistical analyses
The demographic data and cognitive scores of all subjects were analyzed with SPSS 26.0 statistical software (IBM Corp., Armonk, NY, USA). Normal distribution and homogeneity of variance were applied to all quantitative data. The subjects’ age, education, MMSE, and MoCA differences were tested with Kruskal–Wallis, and sex differences were tested with the chi-squared test in the AD, MCI, and normal control (NC) groups.
We first performed group analyses for MD, FA, MK, and CBF of the cortical GM with general linear models, as implemented in FreeSurfer [family-wise error (FWE), p < 0.05] between the NC, MCI, and AD stages. In this work, we used the Monte Carlo simulation with 5,000 repeats to generate a distribution of maximum cluster size. Then, ROI parameters in the deep GM were compared using ANOVA, and the Bonferroni correction was performed in post hoc pairwise comparisons between the sub-groups. The correlation between DKI parameters, CBF, and cognitive scores was analyzed with Spearman’s correlation. Simple linear regression adjusted age and sex was applied to analyze the association between GM CBF and diffusion indices calculated in eight defined GM ROIs.
Finally, we used k-nearest neighbor (KNN) with the Sklearn library in Python to establish classification model and investigate the diagnostic power of DKI and CBF parameters for MCI and AD dementia. For this analysis, we first selected ROIs with f_classif method as predicators for the parameters of each metric, f-value was higher, the diagnostic value was higher. For the classification of AD and NC, the most important ROIs were the left inferior temporal, the left isthmus cingulate and the left lateral occipital for CBF; the most important ROIs were bilateral hippocampus and the right inferior temporal for MD; the most important ROIs were the left entorhinal, the right precentral and the right rostral middle frontal for FA; the most important ROIs were the left caudal anterior cingulate, the left posterior cingulate and the right caudal middle frontal for MK. For the classification of MCI and NC, the most important ROIs were the left entorhinal, the left inferior parietal and the left inferior temporal for CBF; the most important ROIs were the left entorhinal, the right middle temporal and the right superior parietal for MD; the most important ROIs were the left lateral orbital frontal, the right precentral and the left putamen for FA; the most important ROIs were the right fusiform, the right lateral occipital and the left putamen for MK. We then generated mean area under the curve (mAuc), mean accuracy (mAcc), and mean precision (mPre) by running cross_val_score command (cv = 5) for each diagnostic group.
Results
Demographics
Demographic data and neuropsychological score results of the three groups are shown in Table 1. There were no significant differences in sex (p = 0.304), age (p = 0.754), or education duration (p = 0.135) in the AD, MCI, or NC groups. However, there were significant differences in the MMSE (p < 0.001) and MoCA scores (p < 0.001) in the three groups.
Abnormalities in DKI metrics of the cortical and deep gray matter
Compared with the NC, the MCI patients showed areas of increased MD in the precuneus, entorhinal, superior parietal, isthmus cingulate, pars opercularis, and caudal middle frontal lobe in the left hemisphere, and middle temporal, fusiform, precuneus, postcentral, precentral, rostral middle frontal, superior frontal, medial orbitofrontal, lateral orbitofrontal, superior parietal, inferior parietal, and pars opercularis in the right hemisphere (Figure 1A, second row and Supplementary Table 1). The AD patients revealed a widespread increased pattern that included the precuneus, insula, precentral, caudal middle frontal, rostral middle frontal, superior parietal, inferior parietal, middle temporal, and supramarginal regions in the left hemisphere, and medial orbitofrontal, inferior temporal, superior parietal, and lateral occipital regions in the right hemisphere (Figure 1B, second row and Supplementary Table 2). Compared with the MCI, the MD values of the bilateral fusiform, bilateral rostral middle frontal, left posterior cingulate, left lingual, left superior frontal, right insula, right precuneus, right inferior temporal, right caudal middle frontal, and right inferior parietal regions were significantly higher in the AD group (Figure 1C, second row and Supplementary Table 3).
Figure 1. Surface-based statistical map showing the clusters with significant abnormalities in CBF and DKI metrics on the AD continuum. The differences in the CBF and DKI metrics in the NC vs. MCI (A), NC vs. AD (B), and MCI vs. AD (C). Clusters survived correction for multiple comparisons implemented in FreeSurfer by Monte Carlo simulation with 5,000 repeats, with the family-wise error correction settled at p < 0.05. The blue color represents decreased values, and the red color represents increased values.
The MCI patients displayed regions of decreased FA compared to the NC in the left lateral orbitofrontal, right precentral, and right superior temporal areas (Figure 1A, third row and Supplementary Table 1). The AD patients showed a decreased FA in the superior parietal, inferior parietal, precuneus, entorhinal, and posterior cingulate regions in the left hemisphere and in the pars triangularis, precentral, middle temporal, superior frontal, and rostral middle frontal areas in the right hemisphere (Figure 1B, third row and Supplementary Table 2). Compared with the MCI patients, the FA values of the bilateral rostral middle frontal, left temporal pole, left supramarginal, left superior parietal, right precuneus, right postcentral, right caudal middle frontal, and right lateral occipital regions were significantly lower in the AD group (Figure 1C, third row and Supplementary Table 3).
Significant decreases in MK were shown between the MCI and NC groups in the right precuneus and fusiform regions (Figure 1A, fourth row and Supplementary Table 1). As in the AD versus NC comparisons, decreases in MK were noted in the temporal pole, superior temporal, pars opercularis, and caudal anterior cingulate regions in the left hemisphere and in the precuneus, precentral, insula, caudal middle frontal, rostral middle frontal, superior temporal, lateral occipital, inferior parietal, and pericalcarine areas in the right hemisphere (Figure 1B, fourth row and Supplementary Table 2). In addition, compared with the MCI group, the MK value of the left posterior cingulate was significantly lower in the AD group (Figure 1C, fourth row and Supplementary Table 3).
In the deep GM, MD showed the most regions with significant abnormalities in the bilateral amygdala, putamen, and left hippocampus at the MCI stage compared to the NC. FA also revealed changes in the left putamen. In addition to the bilateral amygdala, putamen, and hippocampus, significant differences in MD were shown between the AD and NC groups in the bilateral globus pallidus and left thalamus. Abnormalities of FA occurred in the bilateral putamen, right amygdala, and right thalamus at the AD stage. Similarly, compared with the MCI group, the MD value was significantly increased in AD patients, including the bilateral hippocampus, amygdala, putamen, globus pallidus, and the left thalamus. The FA value was significantly decreased in AD, including in the right putamen. However, no significant changes were found in MK (Table 2 and Figure 2).
Table 2. Group differences of regional CBF and DKI metrics in deep gray matter between NC and MCI/AD.
Figure 2. Box plots of fractional anisotropy (FA) in the bilateral putamen, cerebral blood flow (CBF) in the left caudate and right putamen, and mean diffusivity (MD) in the bilateral hippocampus. Significantly lower FA was observed in the bilateral putamen in the AD group and in the left putamen in the MCI group. In the left caudate and right putamen, significantly lower CBF was observed in the AD group. Significantly higher MD was observed in the bilateral hippocampus in the AD group and in the left hippocampus in the MCI group.
Abnormalities of CBF in the cortical and deep gray matter
Compared with the NC, the MCI patients showed areas of decreased CBF in the bilateral inferior temporal, isthmus cingulate, caudal middle frontal, left lateral occipital, superior parietal, and the right inferior parietal regions (Figure 1A, first row and Supplementary Table 1). As in the AD vs. NC comparisons, decreases in CBF were noted in the bilateral inferior temporal, lateral orbital frontal, medial orbitofrontal, parahippocampus, the left caudal middle frontal, isthmus cingulate, precentral, caudal anterior cingulate, entorhinal, superior temporal, lateral occipital, right precuneus, rostral middle frontal, inferior parietal, temporal pole, inferior parietal, and pericalcarine regions (Figure 1B, first row and Supplementary Table 2). Compared with the MCI, the CBF values of the left lateral occipital, superior temporal, inferior temporal, lingual, supramarginal, inferior parietal, precentral, caudal middle frontal, right fusiform, middle temporal, and lingual regions were significantly lower in the AD group (Figure 1C, first row and Supplementary Table 3).
In the deep GM, CBF showed significant decreases in the bilateral caudate, putamen, and right hippocampus at the AD stage compared to the NC. Compared with the MCI, the CBF values of the left caudate and right putamen were significantly lower in the AD group. There is no significant difference between the MCI and NC groups (Table 2 and Figure 2).
Correlations between imaging parameters and neuropsychological scores
As showed in Figure 3, The FA and MK values of the left putamen (r = 0.490, p < 0.001; r = 0.498, p < 0.001), the MD value of the right middle temporal lobe (r = −0.622, p < 0.001), and the CBF value of the inferior temporal lobe (r = 0.574, p < 0.001) were strongly correlated with the MMSE scores. Moreover, the FA and MK values of the left putamen (r = 0.474, p < 0.001; r = 0.472, p < 0.001), the MD value of the right inferior temporal lobe (r = −0.609, p < 0.001), and the CBF value of the left inferior temporal lobe (r = 0.588, p < 0.001) were strongly correlated with the MoCA scores.
Figure 3. Correlation analyses with maximal correlation coefficients between DKI-derived parameters (MD, FA, and MK), CBF, and neuropsychological testing scores (MMSE scores, MoCA scores). MD, mean diffusivity; FA, fractional anisotropy; MK, mean kurtosis; MMSE, mini-mental state examination; MoCA, Montreal cognitive assessment.
Correlations between DKI parameters and CBF values
As illustrated in Table 3, a higher CBF was associated with a lower MD, higher FA, or higher MK. These associations were statistically significant in most ROIs, including the left frontal lobe [for CBF and FA: β(SE) = 0.9(0.3) × 10–3, p = 0.014], left parietal lobe [for CBF and MD: β(SE) = −3.3(2.0) × 10–3, p = 0.036; for CBF and FA: β(SE) = 0.6(0.2) × 10–3, p = 0.025], right parietal lobe [for CBF and MD: β(SE) = −3.2(0.9) × 10–3, p = 0.019; for CBF and MK: β(SE) = 1.3(0.6) × 10–3, p = 0.041], left temporal lobe [for CBF and MD: β(SE) = −5.5(2.0) × 10–3, p = 0.001], right temporal lobe [for CBF and MD: β(SE) = −5.1(0.7) × 10–3, p < 0.001; for CBF and FA: β(SE) = −0.7(0.2) × 10–3, p = 0.018], left occipital lobe [for CBF and MD: β(SE) = −4.9(0.9) × 10–3, p < 0.001], and right occipital lobe [for CBF and MD: β(SE) = −3.9(1.0) × 10–3, p = 0.004].
Table 3. Regression coefficient and significance in simple linear regression models between CBF and DKI metrics.
KNN analyses for diagnostic performance
An exploratory KNN analysis was undertaken using different MD, FA, MK, and CBF metrics to test the diagnostic performance of these diffusion and perfusion parameters. For distinguishing the AD group from the NC group, MD values performed best (mAuc = 0.939, mAcc = 0.867, mPre = 0.96), MK values performed worst (mAuc = 0.694, mAcc = 0.767, mPre = 0.883). For distinguishing the MCI group from the NC group, CBF values performed best (mAuc = 0.876, mAcc = 0.761, mPre = 0.782), FA values performed worst (mAuc = 0.574, mAcc = 0.523, mPre = 0.550) (Figure 4 and Supplementary Table 4).
Figure 4. A bar graph of mAuc, mAcc, and mPre in the CBF and DKI metrics between the three groups. The mAuc, mAcc, and mPre data of the CBF and DKI metrics to distinguish the AD from the NC group (A) and the MCI from the NC group (B). mAuc, mean area under the curve; mAcc, mean accuracy; mPre, mean precision.
Discussion
This study assessed cortical microstructure and perfusion changes and their relationship in the AD continuum. In addition, the diagnostic efficacy of DKI and CBF parameters was explored. We found the following: (1) the MD values increased significantly, and the FA, MK, and CBF decreased significantly with AD progression, suggesting GM microstructural and perfusion abnormalities in both AD and MCI patients; (2) GM microstructural and perfusion changes are intimately related in AD, providing evidence of associations between regional changes in blood supply and cerebral tissue integrity; (3) the MD, FA, MK, and CBF values of many regions were significantly correlated with the MMSE and MoCA scores, demonstrating that cortical microstructural and blood flow abnormalities associate with cognitive decline; and (4) the CBF values had the best diagnostic power for distinguishing MCI from NC; MD values had the highest diagnostic ability to distinguish AD from NC.
Previous cross-sectional and longitudinal studies on CBF in AD demonstrated that reduced CBF served as an essential imaging biomarker of early disease detection and progression (Asllani et al., 2008; Weijs et al., 2022), and differences in CBF abnormalities were more sensitive than those noted with structural atrophy. In our study, the CBF in numerous regions decreased in AD and MCI patients compared to the NC group, as expected. Furthermore, our study confirmed that the CBF deficit was linked to cognitive impairment. These results were consistent with other AD perfusion reports (Binnewijzend et al., 2013; Ding et al., 2014), showing that CBF decreased gradually with the progression of AD. In addition, one study confirmed that Aβ constricted brain capillaries at pericyte locations, resulting in cerebral perfusion reduction (Nortley et al., 2019). Hence, decreased CBF is closely related to brain vasculature damage. CBF should be considered a potential target for prevention and treatment from a cardiovascular point of view in future research.
We found a significant elevation of MD in the cortical and deep GM. Consistent with previous studies (Rose et al., 2008; Gong et al., 2017), the changes in the cortical GM primarily occurred in the frontal and parietal cortices, including the superior frontal, caudal middle frontal, rostral middle frontal, medial orbital frontal, lateral orbital frontal, precentral gyrus, superior parietal, and inferior parietal lobule in the MCI group; also, the impairment regions were more extensive at the AD stage. MD can be used to examine differences in brain structural integrity. In the present study, increased MD in the cortex suggested early microstructure breakdown due to synaptic and dendritic loss and damage to cellular membranes that occur prior to the development of clinical dementia (Le Bihan, 2003).
As fewer barriers to diffusion appeared, molecules were able to diffuse more freely, and MD was generally expected to increase, implying a decline in FA that was positively associated with neuron density (Preziosa et al., 2019). Notably, we found MD was the most sensitive metric for capturing the loss of microstructural integrity in AD and MCI among the three DKI metrics, as well as serving as meaningful imaging biomarkers for assessing cognitive status. Rodriguez-Vieitez et al. (2021) suggested that higher cortical MD predicted faster hippocampal atrophy and clinical progression to MCI. Furthermore, the MK metric can complement conventional MD and FA for detecting microstructural changes and serving as diagnostic imaging biomarkers.
As with previous studies (Kantarci et al., 2005; Muller et al., 2005; Fellgiebel et al., 2006), the most significant microstructural abnormality reflected by MD in the deep GM were in the hippocampus, which was consistent with the neuropathological features of AD. Douaud et al. (2013) found that left hippocampal MD was the best predictor of future progression to AD. Interestingly, we found MK exhibited poor power at capturing microstructural differences in the deep GM, while previous studies have reported the advantages of diffusional kurtosis metrics in the deep GM (Pelizzari et al., 2019). However, some research has also shown that alterations in kurtosis parameters were limited in AD (Tu et al., 2021; Raj et al., 2022). Wang et al. (2018) found that alterations of MK in subcortical nuclei were weak in AD. We hypothesized that the microstructural changes were covered up by the volume loss. So, in order to comprehensively study the change pattern of microstructure in deep GM in cognitive impairments, future studies need to involve GM volume and microstructure simultaneously.
Robust correlations between DKI parameters, CBF, and cognitive scores were found in many regions, particularly in the MD value of the right inferior temporal lobe and the left hippocampus, verifying that cognition is related to changes in the GM microstructure. Evidence confirmed the association between the microstructural changes of hippocampus and cognitive scores (Shahid et al., 2022). A previous study showed inferior temporal tau burden were associated with activities of daily living impairment in AD (Halawa et al., 2019). Unlike our study, Gong et al. (2013) found an association between GM MD and the MMSE score in the frontal and parietal lobes. While our results were consistent with the sequence of pathology progression in AD (Berron et al., 2021).
We found that a lower CBF was associated with a higher MD, lower FA, and lower MK. These data supported an intimate relationship between tissue integrity and cerebral perfusion changes. In addition, this finding is consistent with the association between CBF and DTI metrics in a cross sectional study reported by Kiely et al. (2022). Meanwhile, a study also found that a positive association between hypoperfusion and white matter damage, and myelin integrity declines with CBF across cognitively normal subjects (Bouhrara et al., 2020). This pattern occurs whether regions are in GM or white matter, to some extent, it reveals that the blood supply offers some protection to cerebral structures. For the results of the article, it is possible that Aβ accumulation leads to vascular injuries, such as vascular wall destruction, pericyte contraction, or neurovascular unit uncoupling (Chakraborty et al., 2017; Sweeney et al., 2018; Apátiga-Pérez et al., 2022), and vascular damage leads to a further CBF reduction. Due to disruptive homeostasis between the blood supply and neuronal activity, damage to the synapses and neurons occurs. These results are biologically plausible, as described previously. In addition, we also found that CBF, MD, and FA changes coexisted in the parietal and temporal lobes. This progression pattern indicated that synergy effects between CBF and brain microstructure may occur early in AD. Therefore, it is essential to explore the extent of this potential association to understand the role of the blood supply and brain tissue integrity in AD.
By analyzing the diagnostic performance of DKI and CBF parameters using the KNN model, we found that the diagnostic value of MD was highest for AD and NC diagnoses, while the diagnostic value of CBF was best for MCI and NC diagnoses. A machine learning study (Collij et al., 2016) showed that CBF can be used to classify the diagnosis of AD disease stages with good accuracy. Furthermore, in the present study, we noted that GM microstructure had great diagnosis value, a finding yet to be extensively discussed in the literature. Some papers explored the predictive performance of hippocampal diffusivity. In a 3-year longitudinal study, Douaud et al. (2013) found that hippocampal diffusivity was more sensitive than white matter tracts and was also a better predictor than CSF Aβ and tau measurements. Similarly, a study suggested that hippocampal MD had high diagnostic accuracy in detecting dementia (Brueggen et al., 2018). The importance of cerebral microstructure should be considered in early disease detection.
This study had some limitations. A limited sample size and unbalanced groups resulted in some bias. Future research needs large-scale, multi-center studies to explore the underlying relationship between brain perfusion and diffusion. Then, we hope to conduct a longitudinal study to trace the trajectory of the cerebral microstructure and blood flow. Only with longitudinal comparisons can we confirm the dynamic performance of these changes with AD progression. Moreover, patients at the preclinical stage (i.e., subjective cognitive decline, SCD) should be enrolled in future studies. In this way, we would be able to confirm that these changes occur earlier and they could become a latent target for early diagnosis and treatment.
Conclusion
In conclusion, this study showed that GM microstructure and CBF were closely related in AD. Increased MD, decreased FA, and MK were accompanied by decreased blood perfusion throughout the AD course. CBF and MD values help predict the AD and MCI diagnoses. These results should be considered in future clinical trials.
Data availability statement
The original contributions presented in this study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving human participants were reviewed and approved by the Tianjin First Central Hospital Medical Ethics Committee. The patients/participants provided their written informed consent to participate in this study.
Author contributions
XN and YG designed the study, collected the data, and wrote the manuscript. ZC contributed to the image processing. TL and YC contributed to the data analysis and interpretation. XZ contributed to the revision of the manuscript. HN participated in the design of the experiment and provided the main idea for the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by the National Natural Scientific Foundation of China (Grant number 30870713, PI: HN), the Tianjin Natural Science Foundation Project (Grant number 21JCZDJC01210, PI: HN), and the Tianjin Health Science and Technology Project (Grant number ZC20161, PI: HN).
Conflict of interest
XZ was employed by Siemens Healthcare Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2023.1205838/full#supplementary-material
Footnotes
- ^ http://www.cereflow.cn
- ^ http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/
- ^ https://www.nitrc.org/projects/dke
- ^ http://surfer.nmr.mgh.harvard.edu/
References
Alisch, J., Kiely, M., Triebswetter, C., Alsameen, M. H., Gong, Z., Khattar, N., et al. (2021). Characterization of age-related differences in the human choroid plexus volume, microstructural integrity, and blood perfusion using multiparameter magnetic resonance imaging. Front. Aging Neurosci. 13:734992. doi: 10.3389/fnagi.2021.734992
Apátiga-Pérez, R., Soto-Rojas, L. O., Campa-Córdoba, B. B., Luna-Viramontes, N. I., Cuevas, E., Villanueva-Fierro, I., et al. (2022). Neurovascular dysfunction and vascular amyloid accumulation as early events in Alzheimer’s disease. Metab. Brain Dis. 37, 39–50. doi: 10.1007/s11011-021-00814-4
Asllani, I., Habeck, C., Scarmeas, N., Borogovac, A., Brown, T. R., and Stern, Y. (2008). Multivariate and univariate analysis of continuous arterial spin labeling perfusion MRI in Alzheimer’s disease. J. Cereb. Blood Flow Metab. 28, 725–736. doi: 10.1038/sj.jcbfm.9600570
Berron, D., Vogel, J. W., Insel, P. S., Pereira, J. B., Xie, L., Wisse, L., et al. (2021). Early stages of tau pathology and its associations with functional connectivity, atrophy and memory. Brain 144, 2771–2783. doi: 10.1093/brain/awab114
Binnewijzend, M. A., Kuijer, J. P., Benedictus, M. R., van der Flier, W. M., Wink, A. M., Wattjes, M. P., et al. (2013). Cerebral blood flow measured with 3D pseudocontinuous arterial spin-labeling MR imaging in Alzheimer disease and mild cognitive impairment: a marker for disease severity. Radiology 267, 221–230. doi: 10.1148/radiol.12120928
Bouhrara, M., Alisch, J., Khattar, N., Kim, R. W., Rejimon, A. C., Cortina, L. E., et al. (2020). Association of cerebral blood flow with myelin content in cognitively unimpaired adults. BMJ Neurol. Open 2:e000053. doi: 10.1136/bmjno-2020-000053
Bouhrara, M., Triebswetter, C., Kiely, M., Bilgel, M., Dolui, S., Erus, G., et al. (2022). Association of cerebral blood flow with longitudinal changes in cerebral microstructural integrity in the Coronary Artery Risk Development in Young Adults (CARDIA) Study. JAMA Netw. Open 5:e2231189. doi: 10.1001/jamanetworkopen.2022.31189
Brueggen, K., Dyrba, M., Kilimann, I., Henf, J., Hoffmann, W., Thyrian, J. R., et al. (2018). Hippocampal mean diffusivity for the diagnosis of dementia and mild cognitive impairment in primary care. Curr. Alzheimer Res. 15, 1005–1012. doi: 10.2174/1567205015666180613114829
Camargo, A., and Wang, Z. (2021). Longitudinal cerebral blood flow changes in normal aging and the Alzheimer’s disease continuum identified by arterial spin labeling MRI. J. Alzheimers Dis. 81, 1727–1735. doi: 10.3233/JAD-210116
Chakraborty, A., de Wit, N. M., van der Flier, W. M., and de Vries, H. E. (2017). The blood brain barrier in Alzheimer’s disease. Vasc. Pharmacol. 89, 12–18. doi: 10.1016/j.vph.2016.11.008
Collij, L. E., Heeman, F., Kuijer, J. P., Ossenkoppele, R., Benedictus, M. R., Moller, C., et al. (2016). Application of machine learning to arterial spin labeling in mild cognitive impairment and alzheimer disease. Radiology 281, 865–875. doi: 10.1148/radiol.2016152703
Ding, B., Ling, H. W., Zhang, Y., Huang, J., Zhang, H., Wang, T., et al. (2014). Pattern of cerebral hyperperfusion in Alzheimer’s disease and amnestic mild cognitive impairment using voxel-based analysis of 3D arterial spin-labeling imaging: initial experience. Clin. Interv. Aging 9, 493–500. doi: 10.2147/CIA.S58879
Douaud, G., Menke, R. A., Gass, A., Monsch, A. U., Rao, A., Whitcher, B., et al. (2013). Brain microstructure reveals early abnormalities more than two years prior to clinical progression from mild cognitive impairment to Alzheimer’s disease. J. Neurosci. 33, 2147–2155. doi: 10.1523/JNEUROSCI.4437-12.2013
Fellgiebel, A., Dellani, P. R., Greverus, D., Scheurich, A., Stoeter, P., and Muller, M. J. (2006). Predicting conversion to dementia in mild cognitive impairment by volumetric and diffusivity measurements of the hippocampus. Psychiatry Res. 146, 283–287. doi: 10.1016/j.pscychresns.2006.01.006
Gong, N., Wong, C., Chan, C., Leung, L., and Chu, Y. (2013). Correlations between microstructural alterations and severity of cognitive deficiency in Alzheimer’s disease and mild cognitive impairment: a diffusional kurtosis imaging study. Magn. Reson. Imaging 31, 688–694. doi: 10.1016/j.mri.2012.10.027
Gong, N. J., Chan, C. C., Leung, L. M., Wong, C. S., Dibb, R., and Liu, C. (2017). Differential microstructural and morphological abnormalities in mild cognitive impairment and Alzheimer’s disease: Evidence from cortical and deep gray matter. Hum. Brain Mapp. 38, 2495–2508. doi: 10.1002/hbm.23535
Greve, D. N., and Fischl, B. (2009). Accurate and robust brain image alignment using boundary-based registration. Neuroimage 48, 63–72. doi: 10.1016/j.neuroimage.2009.06.060
Greve, D. N., Salat, D. H., Bowen, S. L., Izquierdo-Garcia, D., Schultz, A. P., Catana, C., et al. (2016). Different partial volume correction methods lead to different conclusions: An 18F-FDG-PET study of aging. Neuroimage 132, 334–343. doi: 10.1016/j.neuroimage.2016.02.042
Greve, D. N., Svarer, C., Fisher, P. M., Feng, L., Hansen, A. E., Baare, W., et al. (2014). Cortical surface-based analysis reduces bias and variance in kinetic modeling of brain PET data. Neuroimage 92, 225–236. doi: 10.1016/j.neuroimage.2013.12.021
Guo, T., Zhang, D., Zeng, Y., Huang, T. Y., Xu, H., and Zhao, Y. (2020). Molecular and cellular mechanisms underlying the pathogenesis of Alzheimer’s disease. Mol. Neurodegener. 15:40. doi: 10.1186/s13024-020-00391-7
Halawa, O. A., Gatchel, J. R., Amariglio, R. E., Rentz, D. M., Sperling, R. A., Johnson, K. A., et al. (2019). Inferior and medial temporal tau and cortical amyloid are associated with daily functional impairment in Alzheimer’s disease. Alzheimers Res. Ther. 11:14. doi: 10.1186/s13195-019-0471-6
Hays, C. C., Zlatar, Z. Z., and Wierenga, C. E. (2016). The utility of cerebral blood flow as a biomarker of preclinical Alzheimer’s Disease. Cell Mol. Neurobiol. 36, 167–179. doi: 10.1007/s10571-015-0261-z
Hernandez-Garcia, L., Lahiri, A., and Schollenberger, J. (2019). Recent progress in ASL. Neuroimage 187, 3–16. doi: 10.1016/j.neuroimage.2017.12.095
Jensen, J. H., Helpern, J. A., Ramani, A., Lu, H., and Kaczynski, K. (2005). Diffusional kurtosis imaging: the quantification of non-gaussian water diffusion by means of magnetic resonance imaging. Magn. Reson. Med. 53, 1432–1440. doi: 10.1002/mrm.20508
Johnson, N. A., Jahng, G. H., Weiner, M. W., Miller, B. L., Chui, H. C., Jagust, W. J., et al. (2005). Pattern of cerebral hypoperfusion in Alzheimer disease and mild cognitive impairment measured with arterial spin-labeling MR imaging: initial experience. Radiology 234, 851–859. doi: 10.1148/radiol.2343040197
Kantarci, K., Petersen, R. C., Boeve, B. F., Knopman, D. S., Weigand, S. D., O’Brien, P. C., et al. (2005). DWI predicts future progression to Alzheimer disease in amnestic mild cognitive impairment. Neurology 64, 902–904. doi: 10.1212/01.WNL.0000153076.46126.E9
Kiely, M., Triebswetter, C., Gong, Z., Laporte, J. P., Faulkner, M. E., Akhonda, M., et al. (2022). Evidence of an association between cerebral blood flow and microstructural integrity in normative aging using a holistic MRI Approach. J. Magn. Reson. Imaging doi: 10.1002/jmri.28508 [Epub ahead of print].
Korte, N., Nortley, R., and Attwell, D. (2020). Cerebral blood flow decrease as an early pathological mechanism in Alzheimer’s disease. Acta Neuropathol. 140, 793–810. doi: 10.1007/s00401-020-02215-w
Le Bihan, D. (2003). Looking into the functional architecture of the brain with diffusion MRI. Nat. Rev. Neurosci. 4, 469–480. doi: 10.1038/nrn1119
McKhann, G. M., Knopman, D. S., Chertkow, H., Hyman, B. T., Jack, C. J., Kawas, C. H., et al. (2011). 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. 7, 263–269. doi: 10.1016/j.jalz.2011.03.005
Montal, V., Vilaplana, E., Alcolea, D., Pegueroles, J., Pasternak, O., González Ortiz, S., et al. (2018). Cortical microstructural changes along the Alzheimer’s disease continuum. Alzheimers Dement. 14, 340–351. doi: 10.1016/j.jalz.2017.09.013
Muller, M. J., Greverus, D., Dellani, P. R., Weibrich, C., Wille, P. R., Scheurich, A., et al. (2005). Functional implications of hippocampal volume and diffusivity in mild cognitive impairment. Neuroimage 28, 1033–1042.
Nortley, R., Korte, N., Izquierdo, P., Hirunpattarasilp, C., Mishra, A., Jaunmuktane, Z., et al. (2019). Amyloid beta oligomers constrict human capillaries in Alzheimer’s disease via signaling to pericytes. Science 365:aav9518. doi: 10.1126/science.aav9518
Pelizzari, L., Laganà, M. M., Di Tella, S., Rossetto, F., Bergsland, N., Nemni, R., et al. (2019). Combined assessment of diffusion parameters and cerebral blood flow within basal ganglia in early Parkinson’s Disease. Front. Aging Neurosci. 11:134. doi: 10.3389/fnagi.2019.00134
Petersen, R. C., Lopez, O., Armstrong, M. J., Getchius, T., Ganguli, M., Gloss, D., et al. (2018). Practice guideline update summary: Mild cognitive impairment: Report of the Guideline Development, Dissemination, and Implementation Subcommittee of the American Academy of Neurology. Neurology 90, 126–135. doi: 10.1212/WNL.0000000000004826
Preziosa, P., Kiljan, S., Steenwijk, M. D., Meani, A., van de Berg, W., Schenk, G. J., et al. (2019). Axonal degeneration as substrate of fractional anisotropy abnormalities in multiple sclerosis cortex. Brain 142, 1921–1937. doi: 10.1093/brain/awz143
Raj, S., Vyas, S., Modi, M., Garg, G., Singh, P., Kumar, A., et al. (2022). Comparative evaluation of diffusion kurtosis imaging and diffusion tensor imaging in detecting cerebral microstructural changes in alzheimer disease. Acad Radiol. 29(Suppl. 3), S63–S70. doi: 10.1016/j.acra.2021.01.018
Rodriguez-Vieitez, E., Montal, V., Sepulcre, J., Lois, C., Hanseeuw, B., Vilaplana, E., et al. (2021). Association of cortical microstructure with amyloid-beta and tau: Impact on cognitive decline, neurodegeneration, and clinical progression in older adults. Mol. Psychiatry 26, 7813–7822. doi: 10.1038/s41380-021-01290-z
Rose, S. E., Janke, A. L., and Chalk, J. B. (2008). Gray and white matter changes in Alzheimer’s disease: a diffusion tensor imaging study. J. Magn. Reson. Imaging 27, 20–26. doi: 10.1002/jmri.21231
Shahid, S. S., Wen, Q., Risacher, S. L., Farlow, M. R., Unverzagt, F. W., Apostolova, L. G., et al. (2022). Hippocampal-subfield microstructures and their relation to plasma biomarkers in Alzheimer’s disease. Brain 145, 2149–2160. doi: 10.1093/brain/awac138
Smith, S. M., Jenkinson, M., Woolrich, M. W., Beckmann, C. F., Behrens, T. E., Johansen-Berg, H., et al. (2004). Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 23(Suppl. 1), S208–S219. doi: 10.1016/j.neuroimage.2004.07.051
Sweeney, M. D., Sagare, A. P., and Zlokovic, B. V. (2018). Blood-brain barrier breakdown in Alzheimer disease and other neurodegenerative disorders. Nat. Rev. Neurol. 14, 133–150. doi: 10.1038/nrneurol.2017.188
Tabesh, A., Jensen, J. H., Ardekani, B. A., and Helpern, J. A. (2011). Estimation of tensors and tensor-derived measures in diffusional kurtosis imaging. Magn. Reson. Med. 65, 823–836. doi: 10.1002/mrm.22655
Tu, M. C., Huang, S. M., Hsu, Y. H., Yang, J. J., Lin, C. Y., and Kuo, L. W. (2021). Discriminating subcortical ischemic vascular disease and Alzheimer’s disease by diffusion kurtosis imaging in segregated thalamic regions. Hum. Brain Mapp. 42, 2018–2031. doi: 10.1002/hbm.25342
Wang, M., Wei, X., Fu, J., Li, W., Yu, M., Li, P., et al. (2018). Subcortical nuclei in Alzheimer’s disease: a volumetric and diffusion kurtosis imaging study. Acta Radiol. 59, 1365–1371. doi: 10.1177/0284185118758122
Weijs, R. W. J., Shkredova, D. A., Brekelmans, A. C. M., Thijssen, D. H. J., and Claassen, J. A. H. R. (2022). Longitudinal changes in cerebral blood flow and their relation with cognitive decline in patients with dementia: Current knowledge and future directions. Alzheimer Dement. doi: 10.1002/alz.12666
Keywords: microstructure, cerebral blood flow (CBF), diffusion kurtosis imaging (DKI), arterial spin labeling (ASL), Alzheimer’s disease, mild cognitive impairment (MCI), gray matter
Citation: Niu X, Guo Y, Chang Z, Li T, Chen Y, Zhang X and Ni H (2023) The correlation between changes in gray matter microstructure and cerebral blood flow in Alzheimer’s disease. Front. Aging Neurosci. 15:1205838. doi: 10.3389/fnagi.2023.1205838
Received: 14 April 2023; Accepted: 17 May 2023;
Published: 02 June 2023.
Edited by:
Rashid Deane, University of Rochester, United StatesReviewed by:
Imad Rashid Khan, University of Rochester Medical Center, United StatesMustapha Bouhrara, National Institute on Aging (NIH), United States
Copyright © 2023 Niu, Guo, Chang, Li, Chen, Zhang and Ni. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hongyan Ni, hongyni@163.com
†These authors have contributed equally to this work and share first authorship