Lili Huang1,2,3,4†Xin Chen1,2,3,4†Wenshan Sun1,2,3,4
Haifeng Chen1,2,3,4
Qing Ye1,2,3,4Dan Yang1,2,3,4Mengchun Li1,2,3,4Caimei Luo1,2,3,4Junyi Ma1,2,3,4Pengfei Shao1,2,3,4Hengheng Xu1,2,3,4
Bing Zhang5Xiaolei Zhu1,2,3,4*
Yun Xu1,2,3,4*- 1Department of Neurology, Drum Tower Hospital, Medical School and the State Key Laboratory of Pharmaceutical Biotechnology, Institute of Brain Science, Nanjing University, Nanjing, China
- 2Jiangsu Key Laboratory for Molecular Medicine, Medical School of Nanjing University, Nanjing, China
- 3Jiangsu Province Stroke Center for Diagnosis and Therapy, Nanjing, China
- 4Nanjing Neurological Medical Center, Nanjing, China
- 5Department of Radiology, Affiliated Drum Tower Hospital of Nanjing University Medical School, Nanjing, China
Objective: To characterize earlier damage pattern of white matter (WM) microstructure in cerebral small vessel disease (CSVD) and its relationship with cognitive domain dysfunction.
Methods: A total of 144 CSVD patients and 100 healthy controls who underwent neuropsychological measurements and diffusion tensor imaging (DTI) examination were recruited. Cognitive function, emotion, and gait were assessed in each participant. The automated fiber quantification (AFQ) technique was used to extract different fiber properties between groups, and partial correlation and general linear regression analyses were performed to assess the relationship between position-specific WM microstructure and cognitive function.
Results: Specific segments in the association fibers, commissural WM regions of interest (ROIs), and projection fibers were damaged in the CSVD group [P < 0.05, family-wise error (FWE) correction], and these damaged segments showed interhemispheric symmetry. In addition, the damage to specific tract profiles [including the posteromedial component of the right cingulum cingulate (CC), the occipital lobe portion of the callosum forceps major, the posterior portion of the left superior longitudinal fasciculus (SLF), and the bilateral anterior thalamic radiation (ATR)] was related to the dysfunction in specific cognitive domains. Among these tracts, we found the ATR to be the key set of tracts whose profiles were most associated with cognitive dysfunction. The left ATR was a specific fiber bundle associated with episode memory and language function, whereas the fractional anisotropy (FA) values of the intermediate component of the right ATR were negatively correlated with executive function and gait evaluation. It should be noted that the abovementioned relationships could not survive the Bonferroni correction (p < 0.05/27), so we chose more liberal uncorrected statistical thresholds.
Conclusions: Damage to the WM fiber bundles showed extensive interhemispheric symmetry and was limited to particular segments in CSVD patients. Disruption of strategically located fibers was associated with different cognitive deficits, especially the bilateral ATR.
Introduction
Most cerebral small vessel diseases (CSVDs) are associated with age and hypertension. One of the clinical manifestations of the CSVD is cognitive impairment, which may ultimately progress to dementia (Dichgans et al., 2016). The prevalence of dementia in association with CSVD is as high as 36–67% (Roh and Lee, 2014). Deficits in general cognitive function and other cognitive domains, especially information processing speed and executive function, are common in CSVD patients (Prins et al., 2005; Lawrence et al., 2014). The MRI findings in CSVD include white matter hyperintensities (WMH), lacunar infarction (LI), microbleeds, perivascular spaces, and brain atrophy (Wardlaw et al., 2013). A quantitatively dependent relationship was found between the number of lacunas/WMH and the degree of cognitive impairment (CI; Prins and Scheltens, 2015), but variations in clinical symptoms did not completely correspond to conventional MRI markers (Ter Telgte et al., 2018).
Predicting the progression of cognitive decline and its mechanism in the context of CSVD needs to be studied. Increasing evidence has shown a significant association between CSVD's radiological manifestations and dementia pathology. Magnetic resonance-based diffusion tensor imaging (DTI) provides a measure of the diffusion of water in white matter (WM) tracts, allowing researchers to examine the integrity of axonal pathways in patients with CSVD (Beaulieu, 2002). Damage to the microstructural integrity of the WM tract visible in DTI always precedes macroscopic lesions and the destruction of fibers connecting key brain regions (De Groot et al., 2013). There are two parameters to evaluate microstructural integrity of the WM tract, including mean diffusivity (MD), a measure of water diffusion averaged in all spatial directions, and fractional anisotropy (FA), which provides information about the directionality of water diffusion. Loss of microstructural integrity is typically accompanied by a decrease in FA and/or an increase in MD. Increasing evidence has indicated that DTI parameters are early markers of cognitive decline than volume measures (Nir et al., 2013). Voxel-based analysis (VBA) and tract-based spatial statistics (TBSS) are two measures of DTI analysis that have been used to investigate abnormalities in WM diffusion properties in CSVD studies (Della Nave et al., 2007; Papma et al., 2014; Tuladhar et al., 2015). These studies showed that WMH was associated with an increased MD value, when compared with healthy controls, in multiple brain regions, including the corpus callosum, internal and external capsule, cingulum bundle, and other fibers (Papma et al., 2014; Tuladhar et al., 2015). However, these two analysis techniques have drawbacks. It is difficult for VBA to co-register all subjects accurately in a common space due to the huge variations in the shape of the WM fiber bundles among subjects (Wassermann et al., 2011; Pasi et al., 2016). Although TBSS is more accurate and specific than VBA, it fails to precisely define the actual tract location of fibers in an individual (Bach et al., 2014; Pasi et al., 2016). In addition, spurious effects in multiple skeleton points make the results of TBSS less accurate (Pasi et al., 2016).
The new, automated fiber quantification (AFQ) method could capitalize on the precision of tractography for automatically identifying and localizing fiber tracts in individual brains (Yeatman et al., 2012; Chen et al., 2020). In this study, AFQ was applied to evaluate tract profiles of WM integrity in CSVD patients compared with normal controls. FA and MD were extracted at 100 locations along 20 fiber tracts to find damaged fiber tract segments that cause disruption of neural circuits and the locations that are responsible for certain clinical manifestations. These tract profiles visible on DTI may be early imaging markers for CSVD patients with cognitive decline and may provide possible anatomical mechanisms. In addition, the identification of these location-specific properties will offer us a new target for treatment.
Materials and Methods
Study Population
One hundred and forty-four patients with CSVD on MRI were consecutively enrolled between January 2017 and January 2019 from the Department of Neurology of the Affiliated Drum Tower Hospital of Nanjing University Medical School, China. Based on the established research criteria, CSVD is defined as lesions of moderate-to-severe WMH (Prins and Scheltens, 2015; Fazekas scores of 2 or higher) and/or anatomically appropriate lacunar infarcts on neuroimaging, with or without perivascular spaces, microbleeds, and brain atrophy (Wardlaw et al., 2013; Lawrence et al., 2014). The symptoms of CSVD include acute symptoms (e.g., lacunar syndromes and transient ischemic attack) or subacute manifestations (e.g., cognitive impairment, motor disturbances, and emotional disorders). The inclusion criteria were being patients aged between 45 and 85 years and a diagnosis of CSVD based on neuroimaging (WMH and/or lacunes). The exclusion criteria included the presence of any of the following: (1) cerebral infarctions larger than 15 mm in diameter; (2) non-small vascular disease WMH mimics (e.g., multiple sclerosis); (3) intracranial or extracranial large artery stenosis >50%; (4) intracranial hemorrhage; (5) other diseases interfering with neuropsychological tests such as Alzheimer's disease, Parkinson's disease or other psychiatric disorders; (6) MRI contraindications or claustrophobia; and (7) left-handedness. We also recruited 100 healthy individuals from a free clinic or the community during the same time for the healthy control (HC) group. The inclusion criteria comprised: (1)without clinical complaints of cognitive impairment; (2)the scores of clinical scales were within the normal range; (3) without CSVD lesions in MRI. This study was approved by the Ethics Committee of Nanjing Drum Tower Hospital, and written informed consent was obtained from all participants.
Magnetic Resonance Imaging Scanning
All MRI scans were conducted using a Philips 3.0-T scanner (Achieva 3.0 T TX, Philips Medical Systems, The Netherlands) equipped with an eight-channel head coil to diminish scanner noise and restrict head motion. A three-dimensional, high-resolution, T1-weighted turbo gradient echo sequence was performed with the following parameters: repetition time (TR) = 9.8 ms, echo time (TE) = 4.6 ms, flip angle (FA) = 8°, slices = 192, field of view (FOV) = 250 × 250 mm2, acquisition matrix = 256 × 256, and thickness = 1.0 mm. Fluid-attenuated inversion recovery (FLAIR) images were obtained with TR/TE/inversion time (TI) of 4,500/333/1,600 ms, number of slices was 200, voxel size was 0.95 × 0.95 × 0.95 mm3, and acquisition matrix was 270 × 260. Diffusion tensor images were acquired by spin-echo echo planar imaging (TR = 9154 ms, TE = 55 ms, FOV = 224 ×224 mm2, matrix = 112 ×112, slice thickness = 2.5 mm, voxel size = 2 ×2 ×2.5 mm3, with 32 directional diffusion encoding (b = 1,000 s/mm2), and an acquisition without diffusion encoding (b = 0 s/mm2).
Conventional Magnetic Resonance Imaging Analysis
The number of LIs was counted on T1-weighted images and FLAIR sequences by two expert neurologists. The WMH volume, gray matter volume, and intracranial volume (ICV) for each subject was obtained from AccuBrain, a multi-atlas-based anatomical segmentation tool (Wang et al., 2019). The volumes of WMHs and gray matter were divided by the ICV for each participant as a process of normalization (brain volume/ICV ×100 %).
Magnetic Resonance Image Preprocessing
For the three-dimensional (3D) T1 anatomical scans, brain extraction was first performed to delete non-brain structures by a brain extraction tool (BET, a tool in the FSL software). We then used the FMRIB Software Library (FSL version 5.0.9, 1; https://www.fmrib.ox.ac.uk/fsl/) to average and rotate T1 images to align with the AC–PC plane. Raw DTI images were also initially preprocessed by BET to remove non-brain tissue, and then eddy current distortions and head motion were corrected by an affine image registration and eddy current program (tools in the FSL software). Whole-brain diffusion metrics were calculated at each voxel by the DTIFIT program (command in the FSL software), producing FA and MD values.
Automated Fiber Quantification
We used the AFQ software to reconstruct 20 WM fiber tracts in the brain of each subject (Yeatman et al., 2012). The automated calculation was performed by MATLAB R2012b (MathWorks, Inc., Natick, Massachusetts, USA). The AFQ pipeline mainly includes six steps, which were described in Yeatman's paper (Yeatman et al., 2012). The first step is whole-brain deterministic tractography, which is performed on all voxels whose FA > 0.3 in a WM mask. Termination of the tracking occurs when the FA value is < 0.2 or the minimum angle between two consecutive tractography steps is more than 30 degrees. Second, the fiber tract is segmented using waypoint regions of interest (ROIs). Third, the identified fiber tract is refined based on the fiber tract probability maps (Hua et al., 2008). The fourth step is fiber tract cleaning, which filters out fibers more than four standard deviations (SDs) above the mean fiber length or those deviating >5 SDs from the core of the fiber bundle. Fifth, the fiber bundles are clipped at each waypoint ROI to determine the core of the tract. Finally, in the sixth step, 100 equidistant nodes along each tract core are segmented, and diffusivity metrics are extracted using spline interpolation. The 20 identified WM tracts are listed in Table 1. However, we failed to track some fibers, especially the bilateral cingulum hippocampus (CH) and bilateral arcuate fasciculus (AF), because the threshold setting in fiber tracking makes it difficult to identify the tracts adjacent to the gray matter (Johnson et al., 2014; Chen et al., 2020). Therefore, we excluded those four fiber tracts from the analyses.
Table 1. Successful identification rate for each of 20 fiber tracts in HC and CSVD groups.
Neuropsychological Measurements
The neuropsychological evaluation protocol for each subject was completed by a professional neuropsychologist within 1 day of the MRI scans. General cognitive function was assessed using the Mini-Mental State Examination (MMSE) and Beijing version of the Montreal Cognitive Assessment (MoCA). Information processing speed was evaluated by the Stroop Color and Word Tests B (SCWT-B) and the Trail Making Test-A (TMT-A). The SCWT-C and TMT-B are measures of executive function. Episodic memory was evaluated by the visual reproduction–long-delayed recall portion of the Wechsler Memory Scale (WMS-VR-DR) and the Huashan version of the Auditory Verbal Learning Test-delayed recall (AVLT-DR). Language was assessed using the Category Verbal Fluency (CVF) and the Boston Naming Test (BNT). Visuospatial function was evaluated using the Clock Drawing Test (CDT) and the Visual Reproduction-copy (VR-C). It is worth noting that the results from the TMT and SCWT represent time, so they were inversely transformed to maintain consistency. The raw test scores were transformed into Z-scores and averaged to evaluate general cognitive function and function in various cognitive domains (Chen et al., 2019a). In addition, anxiety and depressive disorders were represented by the Hamilton Anxiety Scale (HAMA) and Hamilton Depression Scale (HAMD), respectively. The Timed Up and Go (TUG) test and the Tinetti Balance and Gait Analysis evaluated the fall risk and gait disorders of each subject, respectively. The Z-scores were shown in Table 3, and the raw scores were shown in Supplementary Table 2.
Subjects with CI were distinguished by the MMSE and MoCA scores and their educational experience. The cut-offs of CI for MMSE and MoCA were made based on the studies conducted among the Chinese population. The subject was defined as cognitive impaired when his scores were equal to or lower than the cut-off values of MMSE or MoCA (Gu et al., 2019). (The cut-offs are shown in Supplementary Table 1).
Statistical Analysis
Demographics, clinical information, and neuropsychological data from the HC and CSVD groups were analyzed using SPSS 22.0 (Chicago, IL, USA). The comparisons between normally distributed continuous variables were performed by two sample t-tests, and Mann-Whitney U-tests were applied to compare the non-normally distributed neuropsychological and demographic data. In addition, chi-squared (χ2) tests were used in the comparisons of categorical data, such as sex and the existence of hypertension. The threshold for statistical significance was set at p < 0.05.
To compare the group differences in tract profiles (FA and MD), pointwise analyses were performed based on the “Randomize” program (tool in the FSL software). Age, sex, years of education, and other demographic data that varied between the groups, such as the existence of hypertension and history of LI/ transient ischemic attack (TIA), were included as covariates in the general linear model (GLM). The between-group difference at each node was tested using a permutation-based statistical analysis with 5,000 random permutations. Significance was set at p < 0.05 with a family-wise error (FWE) correction (Sun et al., 2015), and only differences that included three or more adjacent nodes along a tract are reported (Banfi and Koschutnig, 2019).
Then, we calculated the average FA and MD values of 100 nodes along the fiber tracts with differentiated fiber segments to explore different WM fibers between the groups from the level of fiber bundles.
Next, the partial correlations controlling for age, sex, years of education, existence of hypertension, and history of LI/TIA were applied in the CSVD group to assess the relationships between neurofunctional assessments, such as cognitive domain measurements, and all other variables of interest, including WMH volume, gray matter volume, LI count, and tract properties that differed across the groups. Bonferroni correction was conducted to adjust the false-positive rate (p < 0.05/27). To determine whether the integrity of the segmental fiber microstructure was related to the multidimensional clinical symptoms independent of macrostructural burden and to build predictive models of cognitive functions, a multiple linear regression analysis was performed in the CSVD patients. The cognitive domain function was taken as the dependent variable, and the significant factors in the correlation analysis were taken as the independent variables controlling for age, sex, years of education, existence of hypertension, and history of LI/TIA. All tests were two tailed, and the threshold was p < 0.05.
Results
Demographic and Neuropsychological Characteristics
The demographic and clinical data of the HC and CSVD groups are presented in Table 2. There was no significant difference in sex (p = 0.742) or years of education (p = 0.059). The individuals in the CSVD group were significantly older than those in the HC group. When vascular risk factors were compared between the two groups, the CSVD group showed a significantly increased history of hypertension (p = 0.001) and LI/TIA (p < 0.001). Thus, age, history of hypertension, and history of LI/TIA were employed as the covariates in the subsequent analyses.
Table 2. Demographic and clinical characteristics.
LI lesions were present in 59% of the enrolled CSVD patients and nearly 80.6% of the patients presented with WMH. As expected, compared with the HCs, the CSVD patients had more extensive WMH volumes (p < 0.001), a greater number of LIs (p < 0.001), and a significantly reduced volume of gray matter (p < 0.001). Among the CSVD group, 66 subjects were cognitively normal, accounting for 45.8%, and nearly 54.2% were those with CI. Compared with the HC group, the CSVD group also showed poor performance in the general cognitive function and all of the cognitive domains, including information processing speed, executive function, episodic memory, language, and visuospatial function. In addition, they obtained higher scores on the HAMA and HAMD, and significant differences in the two gait evaluations were shown between the two groups. The detailed neuropsychological characteristics are shown in Table 3.
Table 3. Z-scores of neuropsychological measurement.
Fiber Tracking and Group Comparisons at Pointwise Levels
Using the AFQ, we successfully identified 16 fiber tracts in both of the groups, including the bilateral anterior thalamic radiation (ATR), bilateral corticospinal tract (CST), bilateral cingulum cingulate (CC), callosum forceps major, callosum forceps minor, bilateral inferior fronto-occipital fasciculus (IFOF), bilateral inferior longitudinal fasciculus (ILF), bilateral superior longitudinal fasciculus (SLF), and bilateral uncinated fasciculus (UF). Pointwise alterations in FA and MD were analyzed between the groups.
As shown in Figure 1, compared with HCs, the CSVD individuals demonstrated significantly decreased FA in the anterior component of the left ATR (nodes 8–32), the anterior and intermediate components of the right ATR (nodes 1–15; nodes 20–29; and nodes 43–61), the superior portion of the right CST (nodes 65–73), and the occipital lobe portions of the callosum forceps major (nodes 6–14 and nodes 88–91).
Figure 1. Significantly damaged locations of fiber tracts between the HC and CSVD groups in the pointwise comparison of FA profiles (FWE correction, p < 0.05). The blue line indicates the HC group, the orange line indicates the CSVD group, the shadows represent the confidence interval, and the pink fiber segments are the position of the fiber segment that is different between the groups. (A) Significantly decreased FA values in the anterior component of the left ATR (nodes 8–32) in CSVD group. (B) Significantly decreased FA values in the anterior and intermediate components of the right ATR (nodes 1–15; nodes 20–29; and nodes 43–61) in CSVD group. (C) Significantly decreased FA values in the superior portion of the right CST (nodes 65–73) in CSVD group. (D) Significantly decreased FA values in the occipital lobe portions of the callosum forceps major (nodes 6–14 and nodes 88–91) in CSVD group. HC, health control; CSVD, cerebral small vessel disease; FA, fractional anisotropy; ATR, anterior thalamic radiation; CST, corticospinal tract.
In the pointwise comparison of the MD properties, the significantly damaged locations in the fiber tracts (FWE correction, p < 0.05) are shown: (1) the anterior half of the left ATR (nodes 1–56); (2) the anterior half of the right ATR (nodes 1–54); (3) the superior portion of the left CST (nodes 63–100); (4) the superior portion of the right CST (nodes 59–100); (5) the posteromedial component of the right CC (nodes 64–68); (6) the frontal lobe portions of the callosum forceps minor (nodes 10–20 and nodes 73–100); (7) the anterior and intermediate portions of the left IFOF (nodes 27-40 and nodes 80–100); (8) the anterior and intermediate components of the right IFOF (nodes 15–49 and nodes 72–100); (9) the anterior, intermediate and posterior portions of the right ILF (nodes 12–21; nodes 37–45; nodes 59–71; and nodes 98–100); (10) the anterior and posterior portions of the left SLF (nodes 1–43; nodes 78–82; nodes 87–93; and nodes 95–100); and (11) the entire fiber bundle of the right SLF (nodes 1–100). The segments with damaged MD values are shown in Figure 2.
Figure 2. Significantly damaged locations of fiber tracts between the HC and CSVD groups in the pointwise comparison of MD profiles (FWE correction, p < 0.05). The blue line indicates the HC group, the orange line indicates the CSVD group, the shadows represent the confidence interval, and the pink fiber segments are the position of the fiber segment that is different between the groups. The significantly increased MD values in CSVD group (FWE correction, p < 0.05) were located in (A) the anterior half of the left ATR (nodes 1–56); (B) the anterior half of the right ATR (nodes 1–54); (C) the superior portion of the left CST (nodes 63–100); (D) the superior portion of the right CST (nodes 59–100); (E) the posteromedial component of the right CC (nodes 64–68); (F) the frontal lobe portions of the callosum forceps minor (nodes 10–20 and nodes 73–100); (G) the anterior and intermediate portions of the left IFOF (nodes 27–40 and nodes 80–100); (H) the anterior and intermediate components of the right IFOF (nodes 15–49 and nodes 72–100); (I) the anterior, intermediate and posterior portions of the right ILF (nodes 12–21; nodes 37–45; nodes 59–71; and nodes 98–100); (J) the anterior and posterior portions of the left SLF (nodes 1–43; nodes 78–82; nodes 87–93; and nodes 95–100); and (K) the entire fiber bundle of the right SLF (nodes 1–100). HC, health control; CSVD, cerebral small vessel disease; MD, mean diffusivity; ATR, anterior thalamic radiation; CST, corticospinal tract; CC, cingulum cingulate; IFOF, inferior fronto-occipital fasciculus; ILF, inferior longitudinal fasciculus; SLF, superior longitudinal fasciculus.
Group Difference in WM Tract Level
Then, when we compared the diffusion metrics in the WM tract level, the mean FA values of the right CST were no longer different between the groups, whereas other tracts remained significantly different (details in Supplementary Figures 1, 2).
Relationships Between Conventional MRI Markers of CSVD and Neuropsychological Measurements in CSVD Group
Partial correlations were performed between the cognitive domains tests/emotional assessments/gait evaluations and the conventional MRI markers like corrected WMH volumes, gray matter volumes, and LI counts. We only found that general cognitive function was negatively associated with WMH (r = −0.194, p = 0.023; details in Figure 3).
Figure 3. Correlations between traditional CSVD image markers and neuropsychological examinations. Partial correlation was conducted controlling for age, gender, years of education, history of hypertension, and history of LI/TIA in patients with CSVD. General cognitive function was negatively associated with WMH (r = −0.194, p = 0.023). CSVD, cerebral small vessel disease; LI, lacunar infarction; TIA, transient ischemic attack;WMH, white matter hyperintensities.
Relationships Between Diffusion Metrics and Neuropsychological Measurements in CSVD Group
Partial correlations controlling for age, sex, years of education, history of hypertension, and history of LI/TIA were performed between the cognitive domains tests/emotional assessments/gait evaluations and mean diffusion metrics (FA or MD) in each fiber cluster where group differences emerged.
No significant correlations existed between the visuospatial function/HAMD scores and any tract profile. We found significant negative correlations between the general cognitive function scores and the mean MD values of the posteromedial component of the right CC (nodes 64–68, r = −0.186, p = 0.029; Figure 5A). Information processing speeds were positively correlated with the mean FA values of the occipital lobe portions of the callosum forceps major (nodes 6–14, r = 0.228, p = 0.007; nodes 88–91, r = 0.171, p = 0.045; Figures 4A,B); and processing speeds were also negatively correlated with the mean MD values of the frontal lobe portion of the callosum forceps minor (nodes 73–100, r = −0.168, p = 0.049; Figure 5B). Significant positive correlations were observed between the executive function scores and the FA values of the occipital lobe portion of the callosum forceps major (nodes 6–14, r = 0.174, p = 0.042; Figure 4D). It was beyond our expectation that executive function scores were negatively related to the FA values of the intermediate component of the right ATR (nodes 43–61, r = −0.185, p = 0.030; Figure 4C). Higher cortical MD values in the anterior half of the left ATR (nodes 1–56) were correlated with impairments in the episodic memory (r = −0.2, p = 0.019; Figure 5C). Lower FA values in the anterior component of the left ATR (nodes 8–32, r = 0.216, p = 0.011; Figure 4E) and higher MD values in the anterior half of the left ATR (nodes 1–56, r = −0.198, p = 0.020; Figure 5D), frontal lobe portion of the callosum forceps minor (nodes 10–20, r = −0.178, p = 0.037; Figure 5E), and anterior component of the right IFOF (nodes 72–100, r = −0.222, p = 0.009; Figure 5F) showed associations with impairments in the language domain. TUG scores were positively correlated with the FA values of the intermediate component of the right ATR (nodes 43–61, r = 0.187, p = 0.028; Figure 4F). When focusing on anxiety, HAMA scores were negatively related to the MD values of the anterior and posterior portions of the left SLF (nodes 1–43, r = −0.195, p = 0.025; nodes 78–82, r = −0.191, p = 0.026; Figures 5G,H). It is worth noting that the correlations listed above were not able to survive Bonferroni correction (p < 0.05/27), so we chose more liberal uncorrected statistical thresholds.
Figure 4. Correlations between FA profiles and neuropsychological examinations. Partial correlation was conducted controlling for age, gender, years of education, history of hypertension, and history of LI/TIA in patients with CSVD. (A,B) Information processing speed was positively correlated with the mean FA values of the occipital lobe portion of the callosum forceps major (nodes 6–14, r = 0.228, p = 0.007; nodes 88–91, r = 0.171, p = 0.045). (C,D) Executive function was negatively related to the FA value of the intermediate component of the right ATR (nodes 43–61, r = −0.185, p = 0.030), whereas was positively correlated with the FA of the occipital lobe portion of the callosum forceps major (nodes 6–14, r = 0.174, p = 0.042). (E) Language function was positively correlated with the FA of anterior component of the left ATR (nodes 8–32, r = 0.216, p = 0.011). (F) TUG showed positive correlated with the FA value of the intermediate component of the right ATR (nodes 43–61, r = 0.187, p = 0.028). CSVD, cerebral small vessel disease; LI, lacunar infarction; TIA, transient ischemic attack; WMH, white matter hyperintensities; ATR, anterior thalamic radiation.
Figure 5. Correlations between MD profiles and neuropsychological examinations. Partial correlation was conducted controlling for age, gender, years of education, history of hypertension, and history of LI/TIA in patients with CSVD. (A) The mean MD value of the posteromedial component of the right CC (nodes 64–68) correlated negatively with general cognitive function (r = −0.186, p = 0.029). (B) The mean MD value of the frontal lobe portion of the callosum forceps minor (nodes 73–100) correlated negatively with information processing speed (r = −0.168, p = 0.037). (C) Episode memory was negatively associated with the MD value of the anterior half of the left ATR (nodes 1–56, r = −0.2, p = 0.019). (D–F) Higher MD in the anterior half of the left ATR (nodes 1–56, r = −0.198, p = 0.020), frontal lobe portion of the callosum forceps minor (nodes 10–20, r = −0.178, p = 0.037), and anterior component of the right IFOF (nodes 72–100, r = −0.222, p = 0.009) showed significant impairment in the language domain. (G,H) HAMA was negatively related to the MD values of the anterior and posterior portion of left SLF (nodes 1–43, r = −0.195, p = 0.025; nodes 78–82, r = −0.191, p = 0.026). CSVD, cerebral small vessel disease; LI, lacunar infarction; TIA, transient ischemic attack; WMH, white matter hyperintensities; MD, mean diffusivity; CC, cingulum cingulate; ATR, anterior thalamic radiation; IFOF, inferior fronto-occipital fasciculus; SLF, superior longitudinal fasciculus; HAMA, Hamilton Anxiety Scale.
Also, we summarized the correlations between the mean values of different WM fibers and the cognition assessments in Supplementary Figure 3 and detected fewer strategic WM fibers than from the pointwise level analysis. We did not find any significantly damaged WM fibers, which were related to general cognitive function, episodic memory, visuospatial function, TUG, or HAMD. We found the mean FA values of forceps major were significantly related to both executive function (r = 0.181, p = 0.034) and information processing speed (r = 0.259, p = 0.002), and the mean MD values of right CC correlated negatively with information processing speed (r = −0.200, p = 0.0.018). In addition, the mean MD values of left ATR showed an impact on the language function. Besides, only the damage to the left SLF showed an effect on the impairment of anxiety disorder (r = −0.192, p = 0.026).
Multiple Linear Regression Analysis
The multiple regression analysis revealed that the FA value of the occipital lobe portion of the callosum forceps major, the MD value of the anterior half of the left ATR, and the FA value of the anterior component of the left ATR were independent predictors of information processing speed, episodic memory, and language domain scores, respectively. Interestingly, we also found that the FA values of the intermediate component of the right ATR were negatively related to the executive function scores but positively correlated to the TUG scores. HAMA scores were negatively related to the MD values of the posterior portion of the left SLF (details in Table 4).
Table 4. Multiple linear regression analysis.
Discussion
In this study, for the first time, we found specific changes in the microstructural integrity in patients with CSVD by AFQ: (1) WM microstructural damage presented in an extensive interhemispheric symmetrical pattern; (2) the damage did not occur along the entire trajectory of a WM fascicle but along some vulnerable fibrous segments. The pointwise level analysis is more sensitive than the analysis in the WM fiber level where more damaged fiber segments can be found, and more microstructural damages affecting clinical cognition can be detected; and (3) impairments in the different WM portions might mediate different cognitive domain dysfunction, but the anterior middle portion of the ATR causes impairments in multiple recognition domains. These results suggested that the segmental WM fascicle microstructural damage detected by the AFQ may be a marker for the early prediction of CSVD-related cognitive decline, revealing the pathological mechanism of CSVD-related CI from an anatomical perspective.
We first studied the WM microstructural integrity between the CSVD patients and the HCs using the AFQ, which included both the WMH and the normal appearing WM on MRI. The AFQ is a more sensitive tractography method that allows the identification of a complete map of the distribution of cerebral damage in those with CSVD at the individual level. As diffusion properties always vary along the tract trajectory, the AFQ provides a framework for quantifying diffusion metrics at multiple nodes along the trajectory of each fiber tract, generating the “tract profile” (Yeatman et al., 2012). In this study, we observed that the CSVD patients exhibited a significantly damaged tract profile along the multiple bundles, with FA decreased in 25% (4/16) and MD increased in 68.75% (11/16) of the fiber tracts examined. These tracts belong to three fiber systems, namely, the association WM fibers (CC, IFOF, ILF, and SLF), commissural WM ROIs (callosum forceps major and callosum forceps minor), and projection WM fibers (CST and ATR). When compared from the WM fiber level, fewer damaged fibers were found, and that is one of the advantages of AFQ that exploration at the node level can reduce the coverage of some small damage clusters (Chen et al., 2020). As Figures 1, 2 show, most of the microstructural damage showed interhemispheric symmetry and heterogeneity, and this pattern was partly in accordance with macroscopic WMH on T2-weighted images in the CSVD patients (Pantoni, 2010). Mascalchi et al. (2017) focused on CADASIL, an inherited CSVD, and declared an extensive symmetric pattern of fiber damage by DTI analysis. These similar results may be due to extensive neuropathological changes in the leptomeningeal and small penetrating cerebral arteries, such as the loss of smooth muscle cells in the tunica media, deposits of fibrohyaline materials, and thickening of the vessel wall with incomplete lumen occlusion (Pantoni, 2010; Mascalchi et al., 2017). In addition, we found that the diffusion properties varied along the WM fascicle trajectory in the CSVD patients and that there existed some vulnerable nodes or segments along a fiber, consistent with the findings on Alzheimer's disease (Zhang et al., 2019).
Exhilaratingly, this study found that specific WM fascicle segments were related to particular brain functions, such as cognition, emotion, and gait. And compared with strategic WM fibers, more strategic segments were detected, which may have an impact on clinical manifestations. In particular, these segmental correspondences to particular cognitive domains were identified through the multiple linear regression analysis (Table 4). In other words, the specific WM fascicle segment might mediate the corresponding brain dysfunction if the WM microstructure was damaged in CSVD patients.
First, this study found that with the damage of the microstructure of nodes 64–68 in the right CC, the general cognitive function of CSVD patients tended to decline. The cingulum bundle has been implicated in episodic memory, executive control, emotion, and pain. Numerous studies have revealed that cingulum abnormalities could lead to Alzheimer's disease, schizophrenia, depression, posttraumatic stress disorder, and mild cognitive impairment (MCI; Bubb et al., 2018). Our results are supported by some brain perfusion and metabolism studies by PET, that declared reduced metabolism in a network of limbic structures such as the posterior cingulate in MCI patients (Nestor et al., 2003). A resting-state functional MRI study also found that the posterior cingulate cortex is the main connecting hub of the main cortical networks (Buckner et al., 2009).
Second, this study showed that there was a positive relationship between the FA values of nodes 6–14 of the forceps major and information processing speed or executive function. FA was most commonly reported to be associated with processing speed in the genu and body of the corpus callosum in elderly adults (Kerchner et al., 2012). Our findings highlight the occipital lobe portion nodes 6–14 of the callosum forceps major as an anatomical structure potentially impacting processing speed. This is partly supported by Duering et al., who identified this area as strategic fibers affecting processing speed in patients with CADASIL (Duering et al., 2011). We further proposed that segmental damage of the key fiber bundles might be sufficient to reduce information processing speeds due to the disconnection of interhemispheric information transmission.
Notably, the present study first showed that different WM fascicle segments of the ATR might be related to different cognitive domains. The MD values of nodes 1–56 (anterior half) of the left ATR were an independent predictor of episodic memory scores. The FA values of nodes 8–32 (anterior portion) of the left ATR were independently related to language function, and the FA values of nodes 43–61 (medium portion) of the right ATR were independently related to both executive and gait functions.
The ATR begins in the anterior and mediodorsal nuclei of the thalamus and projects to the anterior cingulate cortex and dorsolateral frontal regions (Mamiya et al., 2018), being the core structure of frontothalamic circuitry. Studies have shown that lesions to the anterior nuclei or their afferent WM tracts result in deficits in the encoding of new stimuli and thus affect the declarative memory process (Van Der Werf et al., 2003). The anterior nuclei are involved in the selection of memory materials, and the mediodorsal nucleus is responsible for coordinating and selecting strategies to retrieve declarative memory (Van Der Werf et al., 2003). A PET analysis in those with vascular dementia exhibited a hypometabolic pattern in the frontothalamic circuitry, which was related to memory disruptions (Longarzo et al., 2019). In addition, hypometabolism in patients with amnestic-MCI occurs mainly in the left hemisphere (Seo et al., 2009). Asymmetry in cortical activation has also been suggested, for example, language-based working memory tasks preferentially activated the left hemisphere ( O'Sullivan et al., 2005), which is consistent with our results that only the left ATR was related to memory. Our results may also serve as a support for the existence of multiple memory systems. The concept of an “extended hippocampal system” composed of the hippocampus, mammillary bodies, fornix, and anterior thalamic nuclei was proposed by Aggleton and Saunders (1997). The disruption of these systems leads to deficits in the recall of episodic events, and connections between the thalamus and cortical structures were proposed to be involved in memory processing (Alexinsky, 2001). Disconnections between the prefrontal cortex and the mediodorsal or anterior thalamic nuclei might represent the pathophysiological substrates for aberrant WM microstructure in CSVD-related cognitive decline. Notably, our results emphasized the importance of the integrity of prefrontal-subcortical loop in memory management in patients with CSVD.
In addition, recent research has established that widely distributed, isolated, and shared cerebral regions are involved in both phonological and semantic fluency (Li et al., 2017). These regions are located in the left thalamus, frontal lobe, and parietal lobe (Li et al., 2017). When one produces words in response to particular cues, the left thalamus is activated. The left prefrontal lobe was reported to be involved in category fluency processing (Birn et al., 2010). Therefore, it is not surprising to find that the left ATR, which connects these key brain regions, contributes to language dysfunction (Han et al., 2013). It is also worth noting that the verbal memory and attention/executive control function are the basis of semantic and phonological fluencies (Li et al., 2017). Given the role of the ATR in memory and executive functions, the disconnection of the anterior component of the left ATR could lead to disruptions in verbal fluency (Li et al., 2017).
The dorsolateral prefrontal cortex (DLPFC) is a hub in the executive network, and functional and structural abnormalities in this network contribute to both executive impairments and decreased processing speed (Cummings, 1993; Poole et al., 2018). Ping et al. reported that the integrity of WM fibers connecting the anterior cingulate cortex and the frontal regions was related to test scores in attention; this correlation only existed in the right hemisphere (Mamiya et al., 2018). The involvement of frontothalamic projection fiber tracts in the executive function has been well-established (Tekin and Cummings, 2002) and is in line with our result of the impact of ATR on the executive function and gait evaluation. However, the FA values in the right ATR were negatively related to the executive function and gait. In other words, poor executive function or gait was related to more intact fibers in the right ATR. Some possible explanations were available to expound this seemingly conflicting trend. First, we speculated that the reason for this phenomenon was the plasticity and compensatory effects of WM fibers. Well-organized brain networks have the capacity to adapt to local deficits, and reconstruction of cellular elements within the affected networks after stroke has been reported. For example, the remodeling of functional connections in the photothrombotic infarcts in the murine cortex is often accompanied by dendritic spine turnover and the formation of new structural connections (Brown et al., 2009). The plasticity of fiber microstructure can be compensated temporarily, and we hypothesized that decompensation may occur with the further deterioration of patients' cognitive function, and this definitely needs to be confirmed by our follow-up studies in a larger sample. Next, structural damage may be accompanied by enhanced functional connectivity, which can enhance one's cognition (Chen et al., 2019b). This of course needs further research including functional magnetic resonance. Besides, cognitive reconfiguration could be another explanation. For example, the reconfiguration of the memory process by shifting from the damaged fornix to the left parahippocampal cingulum (PHC) is another confirmation of the compensatory effects of microstructural integrity (Ray and Metzler-Baddeley, 2015). In addition, since the results cannot survive rigorous multiple comparisons, other explanations cannot be ruled out, such as an increase in the FA values caused by the degeneration of crossing-fibers in these regions. In short, a larger sample is needed to conduct a further longitudinal study to validate the abovementioned speculations.
Few fiber tracts other than the posterior portion of the left SLF were associated with emotional disorders. To our surprise, with the emotional deterioration associated with anxiety, the microstructural integrity of the left SLF increased. A study by Lai and Wu involving patients with major depressive disorders reported lower FA in the left SLF of those who were severely depressed (Lai and Wu, 2014), whereas Frodl et al. (2012) found that patients at risk for major depression showed high FA in the left SLF, in agreement with our findings. We hypothesized that the reason for this contrasting result was that the mood disorders of the CSVD subjects in our study were mild or just at risk, and the subjects were resilient. Lower MD has been associated with resilience against anxiety. Resilience to stress was related to the capacity to release corticotrophin-releasing hormone and corticosterone mediated by the expression of brain-derived neurotrophic factor (BDNF) in a negative feedback system (Charney, 2004; Taliaz et al., 2011; Frodl et al., 2012). BDNF is associated with the development of structural neural circuits that control stress adaptations (Taliaz et al., 2011). This may explain the adaptation and plasticity of the left SLF in patients at high risk for mood disorders.
We also analyzed demographic variables in the multiple linear regression model and found that the executive function, information processing speed, HAMA, and TUG scores yielded significant correlations with age, while the episodic memory, language, information processing speed, and executive function scores yielded significant correlations with education. Additionally, sex is also an independent influencing factor for language and anxiety disorders.
Age-related mechanisms of neurodegenerative diseases are impairments in the connectivity of neural networks, which include both primary cortical neuronal degeneration and slowly progressing myelin-degeneration processes such as microvascular dysfunction and β-amyloid deposits (Fellgiebel et al., 2005). As our results suggested, such functionally relevant WM impairments can be sensitively detected and quantified by the AFQ. Education is efficacious on global cognition, select cognitive domains, and psychosocial functioning in people with CI (Kessels et al., 2017; Arce Renteria and Vonk, 2019). The positive effect of education on cognition shows a kind of cognitive reserve function. The cognitive reserve obtained from education cannot protect the brain from pathological changes caused by advanced age but can regulate the relationship between pathology and pathological expression. In other words, cognitive reserve is a kind of accumulation of experience and ability, which may have a positive impact on brain growth, compliance, and compensatory flexibility in the face of injury (Brickman et al., 2011). Studies have confirmed that women have a higher lifetime prevalence rate of emotional and anxiety disorders, which may be related to physiological factors such as sex hormones and, of course, social factors (Riecher-Rossler, 2017). The sex differences in patients with CSVD needs further stratified analyses.
The presence of many macrostructural lesions, measured as WMH volume and LI counts, was not significantly correlated with cognitive domains such as executive function. It should be noted that our study included the spectrum of CSVD patients, and the global dysfunction of the sample was mild. Thus, DTI measures were more sensitive in detecting WM changes and displayed close associations with cognitive function. Microstructural damage detected by DTI precedes the development of WMH on T2-W FLAIR MRI. DTI parameters are more sensitive and serve as early markers of microstructural neurodegenerative changes in CSVD patients (Ye and Bai, 2018). It should also be noted that we included the spectrum of CSVD patients, rather than limiting the study to patients with specific diagnoses such as vascular cognitive dysfunction, with the purpose of generalizing the results to all patients with at least minimal CSVD load, regardless of the final diagnosis.
However, several limitations in this study need to be noted. First, we failed to track some fibers, especially the bilateral CH and bilateral AF, because the threshold setting in the fiber tracking makes it difficult to identify tracts adjacent to the gray matter. The failure to analyze these fibers may have resulted in the loss of important findings. Second, LI is another important clinical outcome of the CSVD, and WM damage was found to be closely related to the incidence of LI (Xu et al., 2018). We did not focus on the effects of strategically located WM fibers on the incidence and recurrence of lacunar infarction. We will follow up with these patients for future studies. Third, this was a cross-sectional study, and longitudinal studies are needed to confirm these findings, especially the changes in the WM tract profiles with the development of CSVD. Last but not the least, in this exploratory work, our microstructural damages were explored from a whole-brain perspective, which involved both macrostructural lesions and normal-appearing white matter. So we cannot confirm if the relationships between cognition and DTI metrics occur in regions outside of lesions. We are now trying to identify and delineate the lesions manually to construct a normal-appearing white matter mask, and more convincing relationships between microstructural damage and cognition will be achieved in future work.
Conclusion
The current study demonstrated that the damage to WM fiber bundles caused by the CSVD showed extensive interhemispheric symmetry and was limited to particular segments. The intermediate component of the right ATR and the anterior half of the left ATR were strategically located fibers affecting CSVD-related CI. In addition, a conflicting trend between the integrity of right ATR and the executive function or gait needs further research, and we speculated that multiple forms of compensation may exist in the brain. These findings contribute to the understanding of the neural mechanisms underlying CI in those with CSVD and shed light on the development of novel markers of CSVD-related CI.
Data Availability Statement
The datasets generated for this study are available on request to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Nanjing Drum Tower Hospital. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
LH, XC, XZ, and YX designed the study, conducted the statistical analysis, and drafted the manuscript. WS, HC, and QY helped to conduct the statistical analysis and conduct the data preprocessing. DY, ML, CL, JM, PS, HX, and BZ helped to collect the data and conduct the data preprocessing. XZ and YX revised the study and manuscript and provided the critical and intellectual content for the study. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by the National Key Research and Development Program of China (2016YFC1300504 and 2016YFC0901004), the National Natural Science Foundation of China (81920108017 and 81630028), the Key Research and Development Program of Jiangsu Province of China (BE2016610), and Jiangsu Province Key Medical Discipline (ZDXKA2016020).
Conflict of Interest
The 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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2020.598242/full#supplementary-material
Supplementary Figure 1. Significantly different WM fibers between the HC and CSVD groups in comparison of mean FA profiles. (A) The mean FA values were significantly decreased in left ATR in CSVD group. (B) The mean FA values were significantly decreased in right ATR in CSVD group. (C) The mean FA values were significantly decreased in Callosum Forcepts Major in CSVD group. HC, health control; CSVD, cerebral small vessel disease; FA, fractional anisotropy; ATR, anterior thalamic radiation. *indicates a statistical difference between groups.
Supplementary Figure 2. Significantly different WM fibers between the HC and CSVD groups in comparison of mean MD profiles. The significant increased mean MD values were located in (A) left ATR, (B) right ATR, (C) left CST, (D) right CST, (E) right CC, (F) Callosum Forcepts Minor, (G) left IFOF, (H) right IFOF, (I) right ILF, (J) left SLF, (K) right SLF. HC, health control; CSVD, cerebral small vessel disease; MD, mean diffusivity; ATR, anterior thalamic radiation; CST, corticospinal tract; CC, cingulum cingulate; IFOF, inferior fronto-occipital fasciculus; ILF, inferior longitudinal fasciculus; SLF, superior longitudinal fasciculus. *indicates a statistical difference between groups.
Supplementary Figure 3. The correlations between the mean FA or MD values and the cognition assessment in the CSVD group. (A) The mean FA values of callosum forceps major were significantly related to executive function (r = 0.181, p = 0.034); (B) The mean FA values of callosum forceps major were significantly related to information processing speed (r = 0.259, p = 0.002); (C) The mean MD values of right CC correlated negatively with information processing speed (r = −0.200, p = 0.0.018). (D) The mean MD values of the left ATR showed an impact on the language function (r = 0.210, P = 0.014); (E) The mean MD values of the left SLF showed a significant negative correlation with HAMA (r = −0.192, p = 0.026).
References
Aggleton, J. P., and Saunders, R. C. (1997). The relationships between temporal lobe and diencephalic structures implicated in anterograde amnesia. Memory 5, 49–71. doi: 10.1080/741941143
Alexinsky, T. (2001). Differential effect of thalamic and cortical lesions on memory systems in the rat. Behav. Brain Res. 122, 175–191. doi: 10.1016/S0166-4328(01)00182-6
Arce Renteria, M., and Vonk, J. M. J. (2019). Illiteracy, dementia risk, and cognitive trajectories among older adults with low education. Neurology 93, e2247–e2256. doi: 10.1212/WNL.0000000000008587
Bach, M., Laun, F. B., Leemans, A., Tax, C. M., Biessels, G. J., Stieltjes, B., et al. (2014). Methodological considerations on tract-based spatial statistics (TBSS). Neuroimage 100, 358–369. doi: 10.1016/j.neuroimage.2014.06.021
Banfi, C., and Koschutnig, K. (2019). White matter alterations and tract lateralization in children with dyslexia and isolated spelling deficits. Hum. Brain Mapp. 40, 765–776. doi: 10.1002/hbm.24410
Beaulieu, C. (2002). The basis of anisotropic water diffusion in the nervous system - a technical review. NMR Biomed. 15, 435–455. doi: 10.1002/nbm.782
Birn, R. M., Kenworthy, L., Case, L., Caravella, R., Jones, T. B., Bandettini, P. A., et al. (2010). Neural systems supporting lexical search guided by letter and semantic category cues: a self-paced overt response fMRI study of verbal fluency. Neuroimage 49, 1099–1107. doi: 10.1016/j.neuroimage.2009.07.036
Brickman, A. M., Siedlecki, K. L., Muraskin, J., Manly, J. J., Luchsinger, J. A., Yeung, L. K., et al. (2011). White matter hyperintensities and cognition: testing the reserve hypothesis. Neurobiol. Aging 32, 1588–1598. doi: 10.1016/j.neurobiolaging.2009.10.013
Brown, C. E., Aminoltejari, K., Erb, H., Winship, I. R., and Murphy, T. H. (2009). In vivo voltage-sensitive dye imaging in adult mice reveals that somatosensory maps lost to stroke are replaced over weeks by new structural and functional circuits with prolonged modes of activation within both the peri-infarct zone and distant sites. J. Neurosci. 29, 1719–1734. doi: 10.1523/JNEUROSCI.4249-08.2009
Bubb, E. J., Metzler-Baddeley, C., and Aggleton, J. P. (2018). The cingulum bundle: anatomy, function, and dysfunction. Neurosci. Biobehav. Rev. 92, 104–127. doi: 10.1016/j.neubiorev.2018.05.008
Buckner, R. L., Sepulcre, J., Talukdar, T., Krienen, F. M., Liu, H., and Hedden, T. (2009). Cortical hubs revealed by intrinsic functional connectivity: mapping, assessment of stability, and relation to Alzheimer's disease. J. Neurosci. 29, 1860–1873. doi: 10.1523/JNEUROSCI.5062-08.2009
Charney, D. S. (2004). Psychobiological mechanisms of resilience and vulnerability: implications for successful adaptation to extreme stress. Am. J. Psychiatry 161, 195–216. doi: 10.1176/appi.ajp.161.2.195
Chen, H., Huang, L., Yang, D., Ye, Q., Guo, M., Qin, R., et al. (2019a). Nodal global efficiency in front-parietal lobe mediated periventricular white matter hyperintensity (PWMH)-related cognitive impairment. Front. Aging Neurosci. 11:347. doi: 10.3389/fnagi.2019.00347
Chen, H. F., Huang, L. L., Li, H. Y., Qian, Y., Yang, D., Qing, Z., et al. (2020). Microstructural disruption of the right inferior fronto-occipital and inferior longitudinal fasciculus contributes to WMH-related cognitive impairment. CNS Neurosci. Ther. 26,576–588. doi: 10.1111/cns.13283
Chen, X., Huang, L., Ye, Q., Yang, D., Qin, R., Luo, C., et al. (2019b). Disrupted functional and structural connectivity within default mode network contribute to WMH-related cognitive impairment. Neuroimage Clin. 24:102088. doi: 10.1016/j.nicl.2019.102088
Cummings, J. L. (1993). Frontal-subcortical circuits and human behavior. Arch. Neurol. 50, 873–880. doi: 10.1001/archneur.1993.00540080076020
De Groot, M., Verhaaren, B. F., De Boer, R., Klein, S., Hofman, A., Van Der Lugt, A., et al. (2013). Changes in normal-appearing white matter precede development of white matter lesions. Stroke 44, 1037–1042. doi: 10.1161/STROKEAHA.112.680223
Della Nave, R., Foresti, S., Pratesi, A., Ginestroni, A., Inzitari, M., Salvadori, E., et al. (2007). Whole-brain histogram and voxel-based analyses of diffusion tensor imaging in patients with leukoaraiosis: correlation with motor and cognitive impairment. AJNR Am. J. Neuroradiol. 28, 1313–1319. doi: 10.3174/ajnr.A0555
Dichgans, M., Wardlaw, J., Smith, E. E., Zietmann, V., Seshadri, S., Sachdev, P., et al. (2016). METACOHORTS for the study of vascular disease and its contribution to cognitive decline and neurodegeneration: an initiative of the Joint Programme for Neurodegenerative Disease Research. Alzheimers Dement. 12, 1235–1249. doi: 10.1016/j.jalz.2016.06.004
Duering, M., Zieren, N., Herve, D., Jouvent, E., Reyes, S., Peters, N., et al. (2011). Strategic role of frontal white matter tracts in vascular cognitive impairment: a voxel-based lesion-symptom mapping study in CADASIL. Brain 134, 2366–2375. doi: 10.1093/brain/awr169
Fellgiebel, A., Muller, M. J., Wille, P., Dellani, P. R., Scheurich, A., Schmidt, L. G., et al. (2005). Color-coded diffusion-tensor-imaging of posterior cingulate fiber tracts in mild cognitive impairment. Neurobiol. Aging 26, 1193–1198. doi: 10.1016/j.neurobiolaging.2004.11.006
Frodl, T., Carballedo, A., Fagan, A. J., Lisiecka, D., Ferguson, Y., and Meaney, J. F. (2012). Effects of early-life adversity on white matter diffusivity changes in patients at risk for major depression. J. Psychiatry Neurosci. 37, 37–45. doi: 10.1503/jpn.110028
Gu, Y., Liu, R., Qin, R., Chen, X., Zou, J., Jiang, Y., et al. (2019). Characteristic changes in the default mode network in hypertensive patients with cognitive impairment. Hypertens. Res. 42, 530–540. doi: 10.1038/s41440-018-0176-4
Han, Z., Ma, Y., Gong, G., He, Y., Caramazza, A., and Bi, Y. (2013). White matter structural connectivity underlying semantic processing: evidence from brain damaged patients. Brain 136, 2952–2965. doi: 10.1093/brain/awt205
Hua, K., Zhang, J., Wakana, S., Jiang, H., Li, X., Reich, D. S., et al. (2008). Tract probability maps in stereotaxic spaces: analyses of white matter anatomy and tract-specific quantification. Neuroimage 39, 336–347. doi: 10.1016/j.neuroimage.2007.07.053
Johnson, R. T., Yeatman, J. D., Wandell, B. A., Buonocore, M. H., Amaral, D. G., and Nordahl, C. W. (2014). Diffusion properties of major white matter tracts in young, typically developing children. Neuroimage 88, 143–154. doi: 10.1016/j.neuroimage.2013.11.025
Kerchner, G. A., Racine, C. A., Hale, S., Wilheim, R., Laluz, V., Miller, B. L., et al. (2012). Cognitive processing speed in older adults: relationship with white matter integrity. PLoS ONE 7:e50425. doi: 10.1371/journal.pone.0050425
Kessels, R. P., Eikelboom, W. S., Schaapsmeerders, P., Maaijwee, N. A., Arntz, R. M., Van Dijk, E. J., et al. (2017). Effect of formal education on vascular cognitive impairment after stroke: a meta-analysis and study in young-stroke patients. J. Int. Neuropsychol. Soc. 23, 223–238. doi: 10.1017/S1355617716001016
Lai, C. H., and Wu, Y. T. (2014). Alterations in white matter micro-integrity of the superior longitudinal fasciculus and anterior thalamic radiation of young adult patients with depression. Psychol. Med. 44, 2825–2832. doi: 10.1017/S0033291714000440
Lawrence, A. J., Chung, A. W., Morris, R. G., Markus, H. S., and Barrick, T. R. (2014). Structural network efficiency is associated with cognitive impairment in small-vessel disease. Neurology 83, 304–311. doi: 10.1212/WNL.0000000000000612
Li, M., Zhang, Y., Song, L., Huang, R., Ding, J., Fang, Y., et al. (2017). Structural connectivity subserving verbal fluency revealed by lesion-behavior mapping in stroke patients. Neuropsychologia 101, 85–96. doi: 10.1016/j.neuropsychologia.2017.05.008
Longarzo, M., Cavaliere, C., Orsini, M., Tramontano, L., Aiello, M., Salvatore, M., et al. (2019). A multimodal imaging study in a case of bilateral thalamic damage with multidomain cognitive impairment. Front. Neurol. 10:1048. doi: 10.3389/fneur.2019.01048
Mamiya, P. C., Richards, T. L., and Kuhl, P. K. (2018). Right forceps minor and anterior thalamic radiation predict executive function skills in young bilingual adults. Front. Psychol. 9:118. doi: 10.3389/fpsyg.2018.00118
Mascalchi, M., Pantoni, L., Giannelli, M., Valenti, R., Bianchi, A., Pracucci, G., et al. (2017). Diffusion tensor imaging to map brain microstructural changes in CADASIL. J. Neuroimaging 27, 85–91. doi: 10.1111/jon.12374
Nestor, P. J., Fryer, T. D., Smielewski, P., and Hodges, J. R. (2003). Limbic hypometabolism in Alzheimer's disease and mild cognitive impairment. Ann. Neurol. 54, 343–351. doi: 10.1002/ana.10669
Nir, T. M., Jahanshad, N., Villalon-Reina, J. E., Toga, A. W., Jack, C. R., Weiner, M. W., et al. (2013). Effectiveness of regional DTI measures in distinguishing Alzheimer's disease, MCI, and normal aging. Neuroimage Clin. 3, 180–195. doi: 10.1016/j.nicl.2013.07.006
O'Sullivan, M., Barrick, T. R., Morris, R. G., Clark, C. A., and Markus, H. S. (2005). Damage within a network of white matter regions underlies executive dysfunction in CADASIL. Neurology 65, 1584–1590. doi: 10.1212/01.wnl.0000184480.07394.fb
Pantoni, L. (2010). Cerebral small vessel disease: from pathogenesis and clinical characteristics to therapeutic challenges. Lancet Neurol. 9, 689–701. doi: 10.1016/S1474-4422(10)70104-6
Papma, J. M., De Groot, M., De Koning, I., Mattace-Raso, F. U., Van Der Lugt, A., Vernooij, M. W., et al. (2014). Cerebral small vessel disease affects white matter microstructure in mild cognitive impairment. Hum. Brain Mapp. 35, 2836–2851. doi: 10.1002/hbm.22370
Pasi, M., Van Uden, I. W., Tuladhar, A. M., De Leeuw, F. E., and Pantoni, L. (2016). White matter microstructural damage on diffusion tensor imaging in cerebral small vessel disease: clinical consequences. Stroke 47, 1679–1684. doi: 10.1161/STROKEAHA.115.012065
Poole, V. N., Wooten, T., Iloputaife, I., Milberg, W., Esterman, M., and Lipsitz, L. A. (2018). Compromised prefrontal structure and function are associated with slower walking in older adults. Neuroimage Clin. 20, 620–626. doi: 10.1016/j.nicl.2018.08.017
Prins, N. D., and Scheltens, P. (2015). White matter hyperintensities, cognitive impairment and dementia: an update. Nat. Rev. Neurol. 11, 157–165. doi: 10.1038/nrneurol.2015.10
Prins, N. D., Van Dijk, E. J., Den Heijer, T., Vermeer, S. E., Jolles, J., Koudstaal, P. J., et al. (2005). Cerebral small-vessel disease and decline in information processing speed, executive function and memory. Brain 128, 2034–2041. doi: 10.1093/brain/awh553
Ray, N. J., and Metzler-Baddeley, C. (2015). Cholinergic basal forebrain structure influences the reconfiguration of white matter connections to support residual memory in mild cognitive impairment. J. Neurosci. 35, 739–747. doi: 10.1523/JNEUROSCI.3617-14.2015
Riecher-Rossler, A. (2017). Sex and gender differences in mental disorders. Lancet Psychiatry 4, 8–9. doi: 10.1016/S2215-0366(16)30348-0
Roh, J. H., and Lee, J. H. (2014). Recent updates on subcortical ischemic vascular dementia. J. Stroke 16, 18–26. doi: 10.5853/jos.2014.16.1.18
Seo, S. W., Cho, S. S., Park, A., Chin, J., and Na, D. L. (2009). Subcortical vascular versus amnestic mild cognitive impairment: comparison of cerebral glucose metabolism. J. Neuroimaging 19, 213–219. doi: 10.1111/j.1552-6569.2008.00292.x
Sun, H., Lui, S., Yao, L., Deng, W., Xiao, Y., Zhang, W., et al. (2015). Two patterns of white matter abnormalities in medication-naive patients with first-episode schizophrenia revealed by diffusion tensor imaging and cluster analysis. JAMA Psychiatry 72, 678–686. doi: 10.1001/jamapsychiatry.2015.0505
Taliaz, D., Loya, A., Gersner, R., Haramati, S., Chen, A., and Zangen, A. (2011). Resilience to chronic stress is mediated by hippocampal brain-derived neurotrophic factor. J. Neurosci. 31, 4475–4483. doi: 10.1523/JNEUROSCI.5725-10.2011
Tekin, S., and Cummings, J. L. (2002). Frontal-subcortical neuronal circuits and clinical neuropsychiatry: an update. J. Psychosom. Res. 53, 647–654. doi: 10.1016/S0022-3999(02)00428-2
Ter Telgte, A., van Leijsen, E. M., Wiegertjes, K., Klijn, C. J., Tuladhar, A. M., and de Leeuw, F. E. (2018). Cerebral small vessel disease: from a focal to a global perspective. Nat. Rev. Neurol. 14, 387–398. doi: 10.1038/s41582-018-0014-y
Tuladhar, A. M., Van Norden, A. G., De Laat, K. F., Zwiers, M. P., Van Dijk, E. J., Norris, D. G., et al. (2015). White matter integrity in small vessel disease is related to cognition. Neuroimage Clin. 7, 518–524. doi: 10.1016/j.nicl.2015.02.003
Van Der Werf, Y. D., Jolles, J., Witter, M. P., and Uylings, H. B. (2003). Contributions of thalamic nuclei to declarative memory functioning. Cortex 39, 1047–1062. doi: 10.1016/S0010-9452(08)70877-3
Wang, C., Zhao, L., Luo, Y., Liu, J., Miao, P., Wei, S., et al. (2019). Structural covariance in subcortical stroke patients measured by automated MRI-based volumetry. Neuroimage Clin. 22:101682. doi: 10.1016/j.nicl.2019.101682
Wardlaw, J. M., Smith, E. E., Biessels, G. J., Cordonnier, C., Fazekas, F., Frayne, R., et al. (2013). Neuroimaging standards for research into small vessel disease and its contribution to ageing and neurodegeneration. Lancet Neurol. 12, 822–838. doi: 10.1016/S1474-4422(13)70124-8
Wassermann, D., Rathi, Y., Bouix, S., Kubicki, M., Kikinis, R., Shenton, M., et al. (2011). White matter bundle registration and population analysis based on Gaussian processes. Inf. Process. Med. Imaging 22, 320–332. doi: 10.1007/978-3-642-22092-0_27
Xu, X., Gao, Y., Liu, R., Qian, L., Chen, Y., Wang, X., et al. (2018). Progression of white matter hyperintensities contributes to lacunar infarction. Aging Dis. 9, 444–452. doi: 10.14336/AD.2017.0808
Ye, Q., and Bai, F. (2018). Contribution of diffusion, perfusion and functional MRI to the disconnection hypothesis in subcortical vascular cognitive impairment. Stroke Vasc. Neurol. 3, 131–139. doi: 10.1136/svn-2017-000080
Yeatman, J. D., Dougherty, R. F., Myall, N. J., Wandell, B. A., and Feldman, H. M. (2012). Tract profiles of white matter properties: automating fiber-tract quantification. PLoS ONE 7:e49790. doi: 10.1371/journal.pone.0049790
Keywords: automated fiber quantification, cerebral small vessel disease, cognitive decline, damage mode, white matter microstructure damage
Citation: Huang L, Chen X, Sun W, Chen H, Ye Q, Yang D, Li M, Luo C, Ma J, Shao P, Xu H, Zhang B, Zhu X and Xu Y (2021) Early Segmental White Matter Fascicle Microstructural Damage Predicts the Corresponding Cognitive Domain Impairment in Cerebral Small Vessel Disease Patients by Automated Fiber Quantification. Front. Aging Neurosci. 12:598242. doi: 10.3389/fnagi.2020.598242
Received: 24 August 2020; Accepted: 07 December 2020;
Published: 11 January 2021.
Edited by:
Boon-Seng Wong, Singapore Institute of Technology, SingaporeReviewed by:
Jeremy Andrew Elman, University of California, San Diego, United StatesDantao Peng, China-Japan Friendship Hospital, China
Copyright © 2021 Huang, Chen, Sun, Chen, Ye, Yang, Li, Luo, Ma, Shao, Xu, Zhang, Zhu and Xu. 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: Yun Xu, xuyun20042001@aliyun.com; Xiaolei Zhu, zhuquelee@126.com
†These authors have contributed equally to this work and share first authorship