Abstract
Background:
Stroke is an important cause of cognitive impairment. Rich club organization, a highly interconnected network brain core region, is closely related to cognition. We hypothesized that the disturbance of rich club organization exists in patients with post-stroke cognitive impairment (PSCI).
Methods:
We collected data on resting-state functional magnetic resonance imaging (rs-fMRI) with 21 healthy controls (HC), 16 hemorrhagic stroke (hPSCI), and 21 infarct stroke (iPSCI). 3D shape features and first-order statistics of stroke lesions were extracted using 3D slicer software. Additionally, we assessed cognitive function using the Montreal Cognitive Assessment (MoCA) and Mini-Mental State Examination (MMSE).
Results:
Normalized rich club coefficients were higher in hPSCI and iPSCI than HC at low-degree k-levels (k = 1–8 in iPSCI, k = 2–8 in hPSCI). Feeder and local connections were significantly decreased in PSCI patients versus HC, mainly distributed in salience network (SN), default-mode network (DMN), cerebellum network (CN), and orbitofrontal cortex (ORB), especially involving the right and left caudate with changed nodal efficiency. The feeder and local connections of significantly between-group difference were positively related to MMSE and MoCA scores, primarily distributed in the sensorimotor network (SMN) and visual network (VN) in hPSCI, SN, and DMN in iPSCI. Additionally, decreased local connections and low-degree ϕnorm(k) were correlated to 3D shape features and first-order statistics of stroke lesions.
Conclusion:
This study reveals the disrupted low-degree level rich club organization and relatively preserved functional core network in PSCI patients. Decreased feeder and local connections in cognition-related networks (DMN, SN, CN, and ORB), particularly involving the caudate nucleus, may offer insight into pathological mechanism of PSCI patients. The shape and signal features of stroke lesions may provide an essential clue for the damage of functional connectivity and the whole brain networks.
Introduction
The stroke prevalence in China was 1,115 cases per 100,000 people in 2013, and it shows the tendency of marked increase (). Most survivors after stroke are accompanied by different types of dysfunction. Among them, post-stroke cognitive impairment (PSCI) is a series of syndromes that meets the diagnostic criteria of cognitive impairment after stroke, containing the deficits of various cognitive domains such as memory, attention, language, executive function, and visual-spatial ability (). About 30% of stroke survivors had cognitive impairment (), which can influence the quality of life with stroke patients and increases the financial pressure of the family ().
Resting-state functional magnetic resonance imaging (fMRI) was used to analyze the dynamic interactions of blood oxygen level-dependent (BOLD) signal, reflected the spontaneous functional activity of neurons among brain regions (; ), which could evaluate cognitive function (). Recently, some studies showed that the alterations in functional connectivity (FC) were closely associated with PSCI. The reduction of FCs across the cerebral hemisphere after acute ischemic stroke was correlated with impaired cognitive performance (; ). The FC analyses of the default-mode network (DMN) were decreased after stroke (; ), which might reveal the pathogenesis of cognitive impairment (). A study investigated the altered FC of the hippocampal subfields belonging to DMN in stroke patients and found the decline of FC between hippocampal and inferior parietal lobule (). Stroke patients had significant impairment in the FCs between the DMN nodes by seed-based connectivity analysis (). Our team found that the fractional amplitude of low-frequency fluctuations and seed-based FC values distributed in DMN in hPSCI and salience network (SN) in iPSCI were decreased compared with healthy volunteers (). Interestingly, after receiving cognitive therapy, the FC in hippocampal increased in PSCI patients (). also found that restored DMN connectivity was relevant to the recovery of cognitive function. These studies mainly paid attention to the DMN connectivity and seed-based FC in PSCI patients. Cognitive impairment caused by stroke might be not only limited to localized brain area damage but also related to the disturbance of neural information transfer in the whole brain network (). However, a few studies investigated the changes in whole-brain network topology based on functional connectivity in PSCI patients.
Graph-theoretic analyses in resting-state fMRI could describe brain network topology and the characteristic of clinical syndromes and assess rich club organization of brain networks (). Rich club organization was densely interconnected in the central hub of a network, promoting interregional neural signaling and integrated information (). We segmented the cortical surface and cerebellum into 116 regions using the Automated Anatomical Labeling (AAL), and each brain region was defined as a node. Some core or hub brain regions are named rich club nodes, and the remaining nodes are peripheral nodes. Rich club connections, highly interconnected networks in brain core regions, could most efficiently transmit neural information in the brain network, which is vital for cognition due to directly regulating brain function across multiple isolated areas and optimizing cognitive processes (). Feeder connections were a core–periphery structure, integrating and transferring information from rich club and local connections. Local connections are connected between periphery regions and created a periphery–periphery structure. Rich club analysis might be helpful in understanding stroke pathophysiology, underlying diagnosis biomarker, and functional recovery (). Rich club analysis has been extensively applied to cognitive impairment caused by Alzheimer’s disease (AD) (; ), cerebral small vessel disease (), and frontal and temporal gliomas (). However, a few studies of rich club analysis on PSCI patients had been conducted.
We analyzed the changes in rich club organization using graph theory and assumed that the disturbance of rich club organization existed in hemorrhagic post-stroke cognitive impairment (hPSCI) and infarct post-stroke cognitive impairment (iPSCI). In addition, we explored the correlations between the clinical performance of PSCI, radiomic features of stroke lesions, and the significant alterations of FCs in patients. We hypothesized that those relationships might be different in various types of stroke.
Materials and Methods
Participants
Fifty-eight Chinese subjects were recruited into this study between January 2019 and December 2020. Among them, 37 cognitive impairment patients were recruited from the Department of Rehabilitation Medicine in Zhongnan Hospital, including 16 hPSCI patients and 21 iPSCI patients. Twenty-one healthy controls (HC) with matched age, sex, and education level were enrolled from the community by advertisements. In conformity with the Helsinki Declaration, the study was approved by the Medical Research Ethics Committee and Institutional Review Board of Zhongnan Hospital (Approval Number: 2019012). All participants provided written informed consent.
All subjects underwent two extensive neuropsychological tests assessed by a professional therapist, including the Beijing version of the Montreal Cognitive Assessment (MoCA) () and the Chinese version of the Mini-Mental State Examination (MMSE) (). The inclusion criteria applied to all PSCI patients included the following: (1) the diagnosis of stroke under the Chinese National Criteria in Diagnostic Essentials of Various Cerebrovascular Diseases revised at the Fourth National Cerebrovascular Diseases Academic Conference in 1995, confirmed by CT or MRI examination (); (2) diagnosis of first-time stroke; (3) onset time ranging 7 days from 3 months in the early subacute period (); (4) age range from 40 to 80 (); (5) reporting mild or moderate impairment, 10 ≤ MMSE score ≤ 27 (>6 years of education) (); (6) right-handed subjects; (7) no severe aphasia and completing the entire experiment; and (8) all participants signed an informed consent form. All patients and HC were required to meet the following exclusion criteria: (1) vital signs were unstable; (2) after craniotomy or skull defect; (3) cognitive decline caused by AD, encephalitis, traumatic brain injury, and Parkinson’s disease; (4) unable to complete neuropsychological tests for severe depression or affective disorders; and (5) contraindications of MRI scanning.
Image Acquisition
All subjects were scanned on a MAGNETOM Trio 3.0 Tesla MR scanner (Siemens, Germany). The scanning protocol contained T1-weighted images and rs-fMRI data using the standard head coil. It was important for subjects that we used headphones and cushions to reduce the noise in scanning. All participants were required to keep their eyes closed, reduce movement, and stay awake. All T1-weighted images were acquired by a 3D magnetization-prepared rapid gradient echo (3D-MPRAGE) sequence for anatomical reference. The scanning parameters were the following: TR/TE = 2,000/2.3 ms, thickness = 1.0 mm, FA = 8°, FOV = 225 mm × 240 mm, and voxel size = 1 mm × 1 mm × 1 mm. All rs-fMRI datasets were obtained with a gradient echo-planar imaging sequence (EPI) with the following parameters: TR/TE = 2,000/30 ms, FOV = 240 mm × 240 mm, flip angle (FA) = 78°, matrix = 64 × 64, thickness = 4.0 mm, number of slices = 35, and voxel size = 3.75 mm × 3.75 mm × 4 mm. Each subject was scanned 480 s and acquired 240 time points.
Image Segmentation
A professional radiologist manually segmented all lesion regions with the T1-weighted images in hPSCI and iPSCI patients using the 3D-Slicer software version 4.13.01. 3D shape features and first-order statistics were extracted with radiomics module. Besides, another radiologist validated the above information. Intraclass correlation coefficient (ICC) was used to assess the consistency of the two radiologists. The reliability was great when the ICC value was higher than 0.75.
Image Preprocessing
Image preprocessing was performed with the DPABI software2 implemented in SPM123, which was run in the environment of MATLAB (Mathworks, Natick, MA, United States). First, we converted the DICOM files to NIFTI images. Second, considering that subjects need to adapt to the scanning environment for a bit of time, the first 10 volumes of BOLD sequences were discarded. Third, the effects of head motion and slice timing were corrected in the remaining volumes of BOLD sequences. Fourth, the functional images were normalized to the Montreal Neurological Institute (MNI) space by using the information of T1 structural images with a resampled voxel size of 3 mm × 3 mm × 3 mm. Fifth, the images underwent spatial smoothing using the full width at half maximum (FWHM) Gaussian kernel of 4 mm × 4 mm × 4 mm to reduce the spatial noise. Finally, the high-frequency noises of respiratory and heartbeat were eliminated by removing the linear trend of time courses and temporally band-pass filtering (0.01-0.08 Hz). Notably, the range of head motion in all subjects was no more than 2° of rotation in the x, y, or z directions, or 2 mm of displacement.
Network Construction
This study used the AAL 116 template to segment the cortical surface and cerebellum into 116 regions (58 regions per hemisphere) as network nodes (). We extracted the averaged time series of the BOLD signal of each node. Then, we computed the Pearson’s correlation coefficients between a node, and the other remaining nodes acted as the edges. In this way, a 116 × 116 functional connection matrix was successfully constructed for each participant. Sparsity threshold values ranged from 0.05 to 0.5 at intervals of 0.01, which removed spurious edges (; Zhu et al., 2021).
Detection of Rich Club Organization
The process of network analysis was performed using the GRETNA toolbox4 in the environment of MATLAB.
The degree (k) was defined as the number of edges that the node shared with other nodes. Rich club coefficients (ϕ(k)) were the density of connections between rich club nodes. The higher the rich club coefficients, the closer the connections between rich club nodes. The formula for calculating ϕ(k) is as follows (; ):
Among them, N represents the number of all nodes in a parcelation. E represents the number of edges. Here, ϕ(k) is defined as the ratio between the total actual number of edges (E>k) and the possible number of maximum edges that connect nodes of degree k or higher over a range of k-values.
As higher-degree nodes are likely to tend to be interconnected to each other, ϕ(k) is typically normalized. Normalized rich club coefficients (ϕnorm(k)) were calculated using ϕ(k) normalized relative to a set of 1,000 comparable random networks (ϕrandom(k)) (; ). ϕnorm(k) was higher than 1 to indicate the existence of a rich club organization in the brain network for each participant (). The formula for calculating ϕnorm(k) is as follows ():
According to the classification of network nodes into the rich club and peripheral regions, the network edges were classified into local connections, connecting between peripheral nodes; feeder connections, connecting one rich club node with one peripheral node, and rich club connections, connecting between rich club nodes.
Statistical Analysis
IBM SPSS performed the statistical analysis for version 235. Mean ± standard deviation represents the numerical variables. We used the Shapiro–Wilk (S–W) test to analyze the normal distribution. Group differences were tested using the chi-square test in sex, lesion side, number of lesions and lesion location, and one-way analysis of variance (ANOVA) in age, years of education, MMSE, and MoCA score. Two-sample t-test was used to test the difference in stroke duration and lesion volume in hPSCI and iPSCI.
We used ANOVA to evaluate group differences in continuous variables of connectivity strength, nodal efficiency, and ϕnorm(k). Then two-sample t-tests (HC group compared with patient groups), acting as post hoc tests, were conducted, with FDR corrected to the p-values to correct for multiple comparisons (). A value of p < 0.05 represented significant difference.
Pearson’s correlations were used to assess the relationship between the clinical performance of cognitive impairment, radiomic features, and FC in each group.
Results
Demographic and Clinical Data
This study recruited 16 hPSCI, 21 iPSCI, and 21 HC. Table 1 shows the demographic and clinical information for all participants. There were no significant group differences for sex, age, and years of education. The stroke duration, lesion side, and number of lesions were not significantly different in hPSCI and iPSCI patients. Significant group differences were found in MMSE and MoCA scores. The patients with hPSCI and iPSCI performed significantly worse than HC on cognitive function, and there was no significant difference between hPSCI and iPSCI. The lesion location in hPSCI patients was principally concentrated on the thalamus and basal ganglia, while iPSCI patients occurred with cortical infarction (p = 0.003). The lesion volume in iPSCI patients was higher than that in hPSCI patients (p = 0.004).
TABLE 1
| Group | HC | hPSCI | iPSCI | p |
| N = 21 | N = 16 | N = 21 | ||
| Sex (male/female) | 17/4 | 12/4 | 18/3 | 0.568 |
| Age (years) | 60.10 ± 6.59 | 60.38 ± 9.78 | 55.81 ± 10.61 | 0.216 |
| Lesion side (left/right) | - | 10/6 | 16/5 | 0.475 |
| Number of lesions (single/multiple) | - | 12/4 | 9/12 | 0.093 |
| Lesion location (cortex/subcortex) | - | 3/13 | 15/6 ‡ | 0.003 |
| Lesion volume (voxel) | - | 25,027 ± 32,775 | 83,149 ± 69,981 ‡ | 0.004 |
| Education (years) | 11.81 ± 2.80 | 11.31 ± 3.22 | 12.24 ± 3.03 | 0.652 |
| Stroke duration (days) | - | 49.93 ± 22.80 | 47.57 ± 50.00 | 0.866 |
| MoCA | 29.10 ± 1.14 | 16.62 ± 4.57 * | 15.48 ± 5.22 † | < 0.001 |
| MMSE | 29.33 ± 0.73 | 20.69 ± 4.66 * | 18.57 ± 5.58 † | < 0.001 |
Baseline of demographic and clinical characteristics.
Values for age, education, MoCA, and MMSE were derived from ANOVA with post hoc test. Sex, lesion side, number of lesions, and lesion location were derived from chi-square test. Two-sample t-tests were used for lesion volume and stroke duration. *Compared with the HC group, p < 0.001 in the hPSCI group; †compared with the HC group, p < 0.001 in the iPSCI group. ‡Compared with the hPSCI group, p < 0.05 in the iPSCI group. HC, healthy controls; hPSCI, hemorrhagic post-stroke cognitive impairment; iPSCI. infarct post-stroke cognitive impairment; MMSE, Mini-Mental State Examination; MoCA, Montreal Cognitive Assessment; ANOVA, one-way analysis of variance.
Rich Club Nodes of All Groups
Rich club organization existed in the brain network in all groups when ϕnorm(k) was greater than 1 (Figure 1A). Most subjects (90%) were found to have the rich club effects at degree levels (k) (; ). We selected the top 17 (15%) highest-degree nodes to represent rich club nodes based on the averaged nodal degree across all groups (Figure 1B). The identified rich club nodes comprised the following: the left and right precuneus, left and right fusiform, left and right lingual, right inferior temporal gyrus, left and right cerebellum_6, left and right middle cingulum, right superior frontal gyrus, left cerebellum_1, right supplementary motor area, left calcarine, and vermis_4, 5, right thalamus. The rest of the nodes were defined as peripheral nodes.
FIGURE 1
The Disturbance of Rich Club Organization
The ϕnorm(k) was higher than 1 at k value ranging from 1 to 36. ϕnorm(k) was significantly higher in hPSCI and iPSCI patients than HC at low-degree k-levels, detailedly k = 1–8 in the iPSCI group and k = 2–8 in the hPSCI group (Figure 1A).
We found that rich club connections had no significant differences in PSCI patients relative to HC. FCs were lower in hPSCI patients than HC with nine feeder connections (Figure 2A and Supplementary Table 1) and 71 local connections (Figure 2B and Supplementary Table 2), which was mainly distributed in the SN, DMN, cerebellum network (CN), and orbitofrontal cortex (ORB). FCs were significantly decreased in iPSCI patients versus HC with 20 feeder connections (Figure 2C and Supplementary Table 3) and 117 local connections (Figure 2D and Supplementary Table 4), primarily distributed in the SN, DMN, CN, and ORB.
FIGURE 2
Whole-Brain Functional Connectivity of Abnormal Nodes
Some nodes of significantly decreased FCs existed in hPSCI and iPSCI versus HC, which might be essential nodes leading to reduced communication. We identified two abnormal nodes with highest edges of aberrant connections according to plotting the frequency of aberrant connections linked by the 116 nodes in hPSCI and iPSCI groups.
Compared with HC, two peripheral nodes with the most declined connections were found in hPSCI, including left caudate (CAU_L) with 31 abnormal edges and right caudate (CAU_R) with 4 abnormal edges. Two peripheral nodes with the most declined connections were found in iPSCI, including CAU_L with 43 abnormal edges and CAU_R with 5 abnormal edges. Group differences were seen in CAU_L efficiency [F (2,55) = 32.51, p < 0.001] and CAU_R efficiency [F (2,55) = 19.94, p < 0.001]. A significant decrease in CAU_L efficiency was observed in hPSCI versus HC (p < 0.001) (Figure 3A). CAU_L efficiency was lower in hPSCI than in iPSCI (p < 0.001) (Figure 3A). CAU_L efficiency was significantly increased in iPSCI patients versus HC (p = 0.0382) (Figure 3A). Compared with HC, the significant declines in CAU_R efficiency were seen in hPSCI (p = 0.0003) and iPSCI groups (p < 0.001) (Figure 3B).
FIGURE 3
Behavioral Correlation Analysis
We first found that some FCs had significant differences between the two patient groups and HC. Then these FCs were positively correlated with the MMSE and MoCA scores, chiefly distributed in the sensorimotor network (SMN) and visual network (VN) in hPSCI (Figure 3C and Supplementary Table 5), SN, and DMN in iPSCI (Figure 3D and Supplementary Table 6).
Radiomic Feature Correlation Analysis
The ICC value was [(0.81 ± 0.037), p < 0.05], representing great consistency for the segmentation of hemorrhage and infraction lesions. The FCs linked to CAU_L and radiomic features were performed with correlation analysis. We found the negative correlation between decreased local connections and most 3D shape features in hPSCI patients, and a positive correlation in most first-order statistics (Figure 4A and Supplementary Table 7). There was a positive correlation between local connections and most 3D shape features, and a negative correlation in most first-order statistics in iPSCI patients (Figure 5 and Supplementary Table 9). Besides, ϕnorm(k) in hPSCI patients (k = 3) and iPSCI patients (k = 2–8) were positively correlated to most 3D shape features and first-order statistics (Figure 4B and Supplementary Table 8).
FIGURE 4
FIGURE 5
Discussion
In this study, we used a graph-theoretical approach to analyze the FCs of whole brain in healthy volunteers, the patients of iPSCI and hPSCI. Our study had some findings, which are as follows: (1) ϕnorm(k) was higher at low-degree k-levels in patients than HC. (2) The undisturbed rich club connections were relatively preserved in patients with PSCI. PSCI tended to affect the feeder and local connections of the SN, DMN, CN, and ORB, especially involving the caudate nucleus (CAU) belonging to the SN. (3) Aberrant CAU efficiencies were observed in the hPSCI and iPSCI groups. (4) Decreased connections involving SMN and VN in hPSCI, and the ones involving SN and DMN in iPSCI, were positively related to the MMSE and MoCA scores. (5) 3D shape features and first-order statistics were correlated with decreased local connections and low-degree ϕnorm(k) in patients.
The Changes in Normalized Rich Club Coefficients
We observed that ϕnorm(k) increased in the lower-degree regions (k = 1-8 in iPSCI and k = 2-8 in hPSCI) compared with the HC group. Our result indicated that rich club nodes in low-degree regions were more likely to be connected to peripheral nodes, which suggested a communication relay to integrate global information by forming locally clustered and segregated communities of the connectome (). Besides, the higher-degree regions found no significant difference in patients, implying that the high-degree k-level regions remain as highly interconnected nodes (). With the value of k degree increased, the main core regions of the brain network were maintained by peeling off low-degree nodes. In our study, most patients with PSCI were mild and moderate, which might only damage some regions at lower-degree (peripheral nodes). considered that ϕnorm(k) was a potential metric to detect the network differences, and found that ϕnorm(k) was increased within low-degree k-levels rather than high-degree in AD and mild cognitive impairment (MCI) patients versus HC, which was consistent with our result. The brain connectivity of peripheral regions may be disrupted rather than rich club connections in patients with mild and moderate cognitive impairment after stroke.
Changes of Rich Club, Feeder, and Local Connections
We found no significant difference in rich club connections, which indicated that rich club connections were relatively stable in hPSCI and iPSCI, which was in accordance with the finding of ϕnorm(k). Rich club connections played an essential role in the integrative brain processes and formed a central backbone for global brain communication, whose interruption severely impacted the whole brain (). A study found that rich club connections kept stable in patients with subjective cognitive decline (). Our study showed that the relatively stable rich club connections in patients with mild and moderate PSCI suggested the retention of partial cognitive function, similar to the damaged pattern of subjective cognitive decline patients.
Our result showed that feeder and local connections were declined in hPSCI and iPSCI patients. Feeder connections could create a core–periphery structure, transferring information transmission based on network among these nodes, integrating and receiving information from various peripheral regions (). Our result indicated that the ability of received and integrated information was declined. Local connections were the periphery–periphery structure by linking non-rich club nodes to non-rich club nodes (). We observed the descend of the transferred information in local connections. Notably, the number of damaged local connections was much greater than that of feeder connections, which manifested that mild and moderate PSCI might mainly damage local connections. Because peripheral regions reduced persistence and lower levels in the hierarchical network, local and feeder connections were easier to damage (). found that the decline in feeder and local connections were main factors leading to the lower efficiency of brain network in MCI and AD, which was in accord with this study. Thus, we speculated that the stroke events might disturb the ability to receive and integrate information of feeder and local connections, eventually causing cognitive impairment.
Additionally, we detected that descended feeder and local connections were majorly distributed in the SN, DMN, CN, and ORB in both patient groups. When individuals were concentrated on the internal tasks rather than the external environment, DMN preferentially was active, and SN was an essential network in the dynamic switching between DMN and central executive network (CEN) (). The SN and DMN played crucial roles in various cognitive domains (; ). confirmed that the damaged interaction between DMN, CEN, and SN was closely related to MCI. Interestingly, we also found that CN connectivity was declined in patients. The output channel from the cerebellum comprises substantial projections to relevant cerebral regions, such as the prefrontal cortex (). The cerebellum could form extensive neural circuits between the cerebrum and cerebellum to modulate cognitive function. showed that cognitive impairment after cerebellar stroke might be related to interrupted brain network. Besides, the connections of neurons in ORB were able to integrate some signals from the external world and internal states to output the thalamus and prefrontal cortex, ongoing cognitive processing (). found that the alteration in sulcogyral patterns in the orbitofrontal cortex occurred in MCI patients. These previous studies were similar to our results. The disturbance of feeder and local connections was primarily distributed in the SN, DMN, CN, and ORB, which may be an important factor leading to cognitive impairment after stroke.
Furthermore, although there was no statistical difference in FCs and clinical performance between the two patient groups, MMSE and MoCA scores were higher in hPSCI than in iPSCI patients. The number of damaged feeder and local connections was greater in iPSCI than in hPSCI, which might be associated with the severity of cognitive impairment, consistent with the result of neuropsychological tests.
Aberrant Connections of Caudate Nucleus
In this study, the CAU was a primary peripheral region of decreased FCs in both patient groups. Besides, the connections of CAU with other nodes were widely disturbed. CAU was an essential part of SN, which was involved in integrating global brain information () and cooperated with the hippocampus to enhance memory performance (). The damage in CAU, through cortical–caudate–cortical loops, may influence cognitive functions (). found that CAU_L and CAU_R contained widely aberrant connections in subjective cognitive decline patients, which was supported by our study.
Our study showed that CAU_L and CAU_R efficiencies were significantly declined in hPSCI patients versus HC. In this study, the lesion location in hPSCI patients was principally the thalamus and basal ganglia. CAU was a part of the basal ganglia. A study showed that the alteration in effective connectivity network in the prefrontal–basal ganglia circuit influenced cognitive function in patients with subcortical stoke (). Remarkably, PSCI patients with subcortical injury had less extra path from the left ventral anterior nucleus in the thalamus to CAU connection (). Therefore, we inferred that the decrease in CAU efficiencies might be related to the damage of the prefrontal–basal ganglia circuit, which was similar to the above previous study.
We found that CAU_L efficiency was increased, and CAU_R efficiency was decreased in iPSCI. In this study, iPSCI patients mainly had cortical infarction, and the CAU was not damaged. We inferred that iPSCI patients might have the compensatory phenomenon, improving nodal efficiency by reducing path length to relieve the onset of clinical symptoms. However, due to the severe brain damage, the clinical symptoms of cognitive impairment still appeared. The functional reorganization of CAU_L might underlie compensatory plasticity mechanisms (; ). However, CAU_R efficiency decreased in iPSCI. The hemispheric lateralization might be an essential factor to modulate cognitive specialization (). A possible reason was that CAU was a region of leftward asymmetry according to AAL atlas (), which might explain the differences in nodal efficiency.
Behavioral Correlation
Our study found the FCs of significant between-group differences in SN and DMN in iPSCI patients, which was positively related to MMSE and MoCA scores. Most iPSCI patients had multiple cortical infarctions. The interruption of SN connectivity was related to the decline in cognitive function in elderly people (). A study found that the FC involved in DMN was decreased in infarct stroke patients to account for the occurrence of PSCI (). We showed that the impairment of DMN and SN connectivity might be closely related to overall cognition in iPSCI, which corresponds to the above previous studies.
Our result showed that the FCs of significant between-group differences in SMN and VN in hPSCI were positively correlated with MMSE and MoCA scores. In this study, the thalamus and basal ganglia were damaged in most hPSCI patients. The thalamus and basal ganglia played roles in motor and cognition through cortico-basal ganglia-thalamo-cortical loops (). The output of basal ganglia targeted the motor cortex, prefrontal cortex, and the frontal eye field, implying that damaged basal ganglia produced not only the alteration in motor function but also the aberration of cognitive impairment (). As an important relayed nucleus, the thalamus was involved in the transmission of motor pathways and cognitive processes (). Thalamus was a first-order subcortical nucleus along the retina–cortical pathway, transmitting peripheral signals to the visual cortex (). The disturbance of the brain network caused by stroke not only existed in the vicinity of the lesion but also occurred in remotely unaffected brain areas (). Our results suggested that the SMN and VN related to motor and cognition might be impaired. A study found that the difference of FCs in early sensory networks, especially involving SMN and VN, was relevant to the difference in age-related cognitive performance (), which was the same as our results. The cognitive impairment caused by the thalamus and basal ganglia stroke was not only related to the damage of the local brain structure but also might be affected by the interference of the brain network mainly in SMN and VN.
Radiomic Feature Correlation
In this study, we found the negative correlation between decreased local connections (CN linked to CAU_L) and 3D shape features, suggesting that the size of hemorrhagic lesions might affect the interruption of the brain network. Most of the infarct lesions were located in the blood supply area of the middle cerebral artery, consistent with the affected cortical network (distributed in DMN, ExN, and ORB). Decreased FCs were positively related to 3D shape features in iPSCI patients, which implied that the larger infarct lesions might be beneficial for repairing decreased FCs because the larger the edema and penumbra, more active neurons may remain.
First-order statistics was derived from statistical moments of the image intensity histogram, which described the distribution of voxel intensity (). The decreased local connections (distributed in CN, SN, ExN, and DMN) were positively correlated to first-order statistics in hPSCI patients, which indicated that the signal loss caused by hemosiderin deposition might damage FCs in the brain network. Infarct regions finally became glial scars with low signal through necrosis, liquefaction, absorption, and repair. We found a negative correlation between decreased local connections (distributed in DMN, ExN, and ORB) and first-order statistics in iPSCI patients, which suggested that the repair of brain tissue might restore the part of FCs in the brain network. Additionally, ϕnorm(k) was positively correlated with 3D shape features and first-order statistics in hPSCI (k = 3) and iPSCI (k = 2–7) patients, indicating that the features of lesions might influence the FCs in low-degree k-level.
Limitations
There were some limitations in the current study. First, this study lacked a control group of non-cognitive impairment after stroke. Second, we only collected the data on the overall cognitive function of two cognitive scales, lacking the analysis of specific cognitive domains, such as attention, memory, and executive function. Third, considering that the sample size may be small, we did not compare the difference between mild and severe cognitive impairment in PSCI patients. Finally, future study should investigate the topological properties in the varying degrees of cognitive impairment and analyze the difference of cognitive impairment caused by the subcortical and cortical lesions.
Conclusion
We find that rich club connections are not damaged in PSCI patients. However, the disturbances in feeder and local connections are more severe in hPSCI and iPSCI than HC, primarily distributed in cognition-related networks (SN, DMN, CN, and ORB). Rich club analysis reveals that the decrease in FCs principally focuses on CAU belonging to SN. The abnormal FCs in brain networks (SMN and VN in hPSCI, SN, and DMN in iPSCI) may be related to cognitive performance. The radiomic features of stroke lesions may imply the damage of FCs and low-degree rich club organization. These findings are helpful for further understanding of the pathogenesis with PSCI.
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.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Ethics statement
The studies involving human participants were reviewed and approved by the Medical Research Ethics Committee and Institutional Review Board of Zhongnan Hospital. The patients/participants provided their written informed consent to participate in this study.
Author contributions
GM contributed to the writing of the original draft, conceptualization, methodology, and performed the formal analysis. BR contributed to the conceptualization, performed formal analysis, and wrote, reviewed, and edited the draft. SW wrote, reviewed, and edited the draft and conducted the formal analysis. PF, ZC, XZ, and JZ conducted the investigation. LC handled the visualization. HX and WL handled the supervision and funding acquisition. All authors contributed to the article and approved the submitted version.
Funding
This study was supported in part by the National Key Technology Research and Development Program of China (2018YFC2002300) and the National Natural Science Foundation of China (Grant No. 81771819).
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/fnins.2021.796530/full#supplementary-material
References
1
BallG.AljabarP.ZebariS.TusorN.ArichiT.MerchantN.et al (2014). Rich-club organization of the newborn human brain.Proc. Natl. Acad. Sci.1117456–7461. 10.1073/pnas.1324118111
2
BernhardtJ.HaywardK. S.KwakkelG.WardN. S.WolfS. L.BorschmannK.et al (2017). Agreed definitions and a shared vision for new standards in stroke recovery research: the Stroke Recovery and Rehabilitation Roundtable taskforce.Int. J. Stroke12444–450. 10.1177/1747493017711816
3
BiglerE. D. (2014). Magnetic resonance imaging in the evaluation of cognitive function.Pediatr. Blood Cancer611724–1728. 10.1002/pbc.25110
4
BucknerR. L. (2013). The cerebellum and cognitive function: 25 years of insight from anatomy and neuroimaging.Neuron80807–815. 10.1016/j.neuron.2013.10.044
5
BurtonL.TysonS. (2015). Screening for cognitive impairment after stroke: a systematic review of psychometric properties and clinical utility.J. Rehabil. Med.47193–203. 10.2340/16501977-1930
6
CaeyenberghsK.LeemansA. (2014). Hemispheric lateralization of topological organization in structural brain networks.Hum. Brain Mapp.354944–4957. 10.1002/hbm.22524
7
CaoR.WangX.GaoY.LiT.ZhangH.HussainW.et al (2020). Abnormal Anatomical Rich-Club Organization and Structural–Functional Coupling in Mild Cognitive Impairment and Alzheimer’s Disease.Front. Neurol.11:53. 10.3389/fneur.2020.00053
8
CarterA. R.AstafievS. V.LangC. E.ConnorL. T.RengacharyJ.StrubeM. J.et al (2010). Resting interhemispheric functional magnetic resonance imaging connectivity predicts performance after stroke.Ann. Neurol.67365–375. 10.1002/ana.21905
9
ChandG. B.WuJ.HajjarI.QiuD. (2017). Interactions of the Salience Network and Its Subsystems with the Default-Mode and the Central-Executive Networks in Normal Aging and Mild Cognitive Impairment.Brain Connect.7401–412. 10.1089/brain.2017.0509
10
ChuaX. Y.ChooR. W. M.HaN. H. L.CheongC. Y.WeeS. L.YapP. L. K. (2019). Mapping modified Mini-Mental State Examination (MMSE) scores to dementia stages in a multi-ethnic Asian population.Int. Psychogeriatr.31147–151. 10.1017/S1041610218000704
11
Dacosta AguayoR.GrañaM.Iturria MedinaY.Fernández AndújarM.López CancioE.CáceresC.et al (2015). Impairment of functional integration of the default mode network correlates with cognitive outcome at three months after stroke.Hum. Brain Mapp.36577–590. 10.1002/hbm.22648
12
DaianuM.JahanshadN.NirT. M.JackC. R.WeinerM. W.BernsteinM. A.et al (2015). Rich club analysis in the Alzheimer’s disease connectome reveals a relatively undisturbed structural core network.Hum. Brain Mapp.363087–3103. 10.1002/hbm.22830
13
DaianuM.MezherA.MendezM. F.JahanshadN.JimenezE. E.ThompsonP. M. (2016). Disrupted rich club network in behavioral variant frontotemporal dementia and early-onset Alzheimer’s disease.Hum. Brain Mapp.37868–883. 10.1002/hbm.23069
14
DingX.LiC. Y.WangQ. S.DuF. Z.KeZ. W.PengF.et al (2014). Patterns in default-mode network connectivity for determining outcomes in cognitive function in acute stroke patients.Neuroscience277637–646. 10.1016/j.neuroscience.2014.07.060
15
FuhJ. L.WangS. J. (1995). Caudate hemorrhage: clinical features, neuropsychological assessments and radiological findings.Clin. Neurol. Neurosurg.97296–299. 10.1016/0303-8467(95)00059-s
16
GouldenN.KhusnulinaA.DavisN. J.BracewellR. M.BokdeA. L.McNultyJ. P.et al (2014). The salience network is responsible for switching between the default mode network and the central executive network: replication from DCM.NeuroImage99180–190. 10.1016/j.neuroimage.2014.05.052
17
HeX.QinW.LiuY.ZhangX.DuanY.SongJ.et al (2014). Abnormal salience network in normal aging and in amnestic mild cognitive impairment and Alzheimer’s disease.Hum. Brain Mapp.353446–3464. 10.1002/hbm.22414
18
JiangL.GengW.ChenH.ZhangH.BoF.MaoC.et al (2018). Decreased functional connectivity within the default-mode network in acute brainstem ischemic stroke.Eur. J. Radiol.105221–226. 10.1016/j.ejrad.2018.06.018
19
JokinenH.MelkasS.YlikoskiR.PohjasvaaraT.KasteM.ErkinjunttiT.et al (2015). Post-stroke cognitive impairment is common even after successful clinical recovery.Eur. J. Neurol.221288–1294. 10.1111/ene.12743
20
JungJ.LaverickR.NaderK.BrownT.MorrisH.WilsonM.et al (2021). Altered hippocampal functional connectivity patterns in patients with cognitive impairments following ischaemic stroke: a resting-state fMRI study.Neuroimage Clin.32:102742. 10.1016/j.nicl.2021.102742
21
KimD.MinB. (2020). Rich-club in the brain’s macrostructure: insights from graph theoretical analysis.Comput. Struct. Biotec.181761–1773. 10.1016/j.csbj.2020.06.039
22
KwonH. S.LeeD.LeeM. H.YuS.LimJ.YuK.et al (2020). Post-stroke cognitive impairment as an independent predictor of ischemic stroke recurrence: Picasso sub-study.J. Neurol.267688–693. 10.1007/s00415-019-09630-4
23
LanciegoJ. L.LuquinN.ObesoJ. A. (2012). Functional Neuroanatomy of the Basal Ganglia.Cold Spring Harb. Perspect. Med.2:a009621. 10.1101/cshperspect.a009621
24
LimJ.KangD. (2015). Stroke Connectome and Its Implications for Cognitive and Behavioral Sequela of Stroke.J. Stroke17256–267. 10.5853/jos.2015.17.3.256
25
LiuX.HeC.FanD.ZhuY.ZangF.WangQ.et al (2021). Disrupted rich-club network organization and individualized identification of patients with major depressive disorder.Prog. Neuropsychopharmacol. Biol. Psychiatry108:110074. 10.1016/j.pnpbp.2020.110074
26
LiuY.YangK.HuX.XiaoC.RaoJ.LiZ.et al (2020). Altered Rich-Club Organization and Regional Topology Are Associated With Cognitive Decline in Patients With Frontal and Temporal Gliomas.Front. Hum. Neurosci.14:23. 10.3389/fnhum.2020.00023
27
MedagliaJ. D. (2017). Graph Theoretic Analysis of Resting State Functional MR Imaging.Neuroimag. Clin. N Am.27593–607. 10.1016/j.nic.2017.06.008
28
MiddletonF. A.StrickP. L. (2000). Basal ganglia and cerebellar loops: motor and cognitive circuits 1.Brain Res. Rev.31236–250. 10.1016/s0165-0173(99)00040-5
29
MüllerN. C. J.KonradB. N.KohnN.Müoz-LópezM.CzischM.FernándezG.et al (2018). Hippocampal–caudate nucleus interactions support exceptional memory performance.Brain Struct. Funct.2231379–1389. 10.1007/s00429-017-1556-2
30
NeurologyC.NeurosurgeryC. (1996). Classification, diagnostic criteria and evaluation of neurological impairment for stroke patients.Zhonghua Shen Jing Ge Za Zhi29379–380.
31
ParkJ.KimY.ChangW. H.ParkC.ShinY.KimS. T.et al (2014). Significance of longitudinal changes in the default-mode network for cognitive recovery after stroke.Eur. J. Neurosci.402715–2722. 10.1111/ejn.12640
32
PiniL.YoussovK.SambataroF.Bachoud LeviA. C.VallesiA.JacquemotC. (2020). Striatal connectivity in pre-manifest Huntington’s disease is differentially affected by disease burden.Eur. J. Neurol.272147–2157. 10.1111/ene.14423
33
PuigJ.BlascoG.Alberich-BayarriA.SchlaugG.DecoG.BiarnesC.et al (2018). Resting-State Functional Connectivity Magnetic Resonance Imaging and Outcome After Acute Stroke.Stroke492353–2360. 10.1161/STROKEAHA.118.021319
34
RehmeA. K.GrefkesC. (2013). Cerebral network disorders after stroke: evidence from imaging-based connectivity analyses of active and resting brain states in humans.J. Physiol.59117–31. 10.1113/jphysiol.2012.243469
35
RikhyeR. V.WimmerR. D.HalassaM. M. (2018). Toward an Integrative Theory of Thalamic Function.Annu. Rev. Neurosci.41163–183. 10.1146/annurev-neuro-080317-062144
36
RubinovM.SpornsO. (2010). Complex network measures of brain connectivity: uses and interpretations.Neuroimage521059–1069. 10.1016/j.neuroimage.2009.10.003
37
RudebeckP. H.RichE. L. (2018). Orbitofrontal cortex.Curr. Biol.28R1083–R1088. 10.1016/j.cub.2018.07.018
38
SaalmannY. B.KastnerS. (2011). Cognitive and Perceptual Functions of the Visual Thalamus.Neuron71209–223. 10.1016/j.neuron.2011.06.027
39
SpornsO. (2013). Structure and function of complex brain networks.Dialogues Clin. Neurosci.15247–262.
40
StoccoA.LebiereC.O’ReillyR. C.AndersonJ. R. (2012). Distinct contributions of the caudate nucleus, rostral prefrontal cortex, and parietal cortex to the execution of instructed tasks.Cogn. Affect. Behav. Neurosci.12611–628. 10.3758/s13415-012-0117-7
41
StummeJ.JockwitzC.HoffstaedterF.AmuntsK.CaspersS. (2020). Functional network reorganization in older adults: graph-theoretical analyses of age, cognition and sex.Neuroimage214:116756. 10.1016/j.neuroimage.2020.116756
42
SunJ.TanL.YuJ. (2014). Post-stroke cognitive impairment: epidemiology, mechanisms and management.Ann. Transl. Med.2:80. 10.3978/j.issn.2305-5839.2014.08.05
43
TuladharA. M.SnaphaanL.ShumskayaE.RijpkemaM.FernandezG.NorrisD. G.et al (2013). Default Mode Network Connectivity in Stroke Patients.PLoS One8:e66556. 10.1371/journal.pone.0066556
44
Tzourio-MazoyerN.LandeauB.PapathanassiouD.CrivelloF.EtardO.DelcroixN.et al (2002). Automated Anatomical Labeling of Activations in SPM Using a Macroscopic Anatomical Parcellation of the MNI MRI Single-Subject Brain.Neuroimage15273–289. 10.1006/nimg.2001.0978
45
van den HeuvelM. P.SpornsO. (2011). Rich-Club Organization of the Human Connectome.J. Neurosci.3115775–15786. 10.1523/JNEUROSCI.3539-11.2011
46
van den HeuvelM. P.SpornsO.CollinG.ScheeweT.MandlR. C. W.CahnW.et al (2013). Abnormal Rich Club Organization and Functional Brain Dynamics in Schizophrenia.JAMA Psychiat.70:783. 10.1001/jamapsychiatry.2013.1328
47
van LeijsenE.van UdenI.BergkampM. I.van der HolstH. M.NorrisD. G.ClaassenJ.et al (2019). Longitudinal changes in rich club organization and cognition in cerebral small vessel disease.Neuroimage Clin.24:102048. 10.1016/j.nicl.2019.102048
48
WangD.YaoQ.YuM.XiaoC.FanL.LinX.et al (2019). Topological Disruption of Structural Brain Networks in Patients With Cognitive Impairment Following Cerebellar Infarction.Front. Neurol.10:759. 10.3389/fneur.2019.00759
49
WangS.RaoB.ChenL.ChenZ.FangP.MiaoG.et al (2021). Using Fractional Amplitude of Low-Frequency Fluctuations and Functional Connectivity in Patients With Post-stroke Cognitive Impairment for a Simulated Stimulation Program.Front. Aging Neurosci.13:724267. 10.3389/fnagi.2021.724267
50
WangZ.ZhangX.LiuR.WangY.QingZ.LuJ.et al (2020). Altered sulcogyral patterns of orbitofrontal cortex in patients with mild cognitive impairment.Psychiatry Res. Neuroimag.302:111108. 10.1016/j.pscychresns.2020.111108
51
WuS.WuB.LiuM.ChenZ.WangW.AndersonC. S.et al (2019). Stroke in China: advances and challenges in epidemiology, prevention, and management.Lancet Neurol.18394–405. 10.1016/S1474-4422(18)30500-3
52
XiaoF.LuC.ZhaoD.ZouQ.XuL.LiJ.et al (2019). Independent Component Analysis and Graph Theoretical Analysis in Patients with Narcolepsy.Neurosci. Bull.35743–755. 10.1007/s12264-018-0307-6
53
YanT.WangW.YangL.ChenK.ChenR.HanY. (2018). Rich club disturbances of the human connectome from subjective cognitive decline to Alzheimer’s disease.Theranostics83237–3255. 10.7150/thno.23772
54
YangS.JiangC.YeH.TaoJ.HuangJ.GaoY.et al (2014). Effect of Integrated Cognitive Therapy on Hippocampal Functional Connectivity Patterns in Stroke Patients with Cognitive Dysfunction: a Resting-State fMRI Study.Evid. Based Compl. Alt.2014962304–962309. 10.1155/2014/962304
55
YipS. S. F.AertsH. J. W. L. (2016). Applications and limitations of radiomics.Phys. Med. Biol.61R150–R166. 10.1088/0031-9155/61/13/R150
56
ZhangJ.LiZ.CaoX.ZuoL.WenW.ZhuW.et al (2020). Altered Prefrontal–Basal Ganglia Effective Connectivity in Patients With Poststroke Cognitive Impairment.Front. Neurol.11:577482. 10.3389/fneur.2020.577482
57
ZhangT.ZhangY.RenJ.YangC.ZhouH.LiL.et al (2021). Aberrant basal ganglia-thalamo-cortical network topology in juvenile absence epilepsy: a resting-state EEG-fMRI study.Seizure8478–83. 10.1016/j.seizure.2020.11.015
58
ZhuJ.ZengQ.ShiQ.LiJ.DongS.LaiC.et al (2021). Altered Brain Functional Network in Subtypes of Parkinson’s Disease: a Dynamic Perspective.Front. Aging Neurosci.13:710735. 10.3389/fnagi.2021.710735
Summary
Keywords
post-stroke cognitive impairment, resting-state functional magnetic resonance imaging, functional connectivity, rich club organization, graph theory
Citation
Miao G, Rao B, Wang S, Fang P, Chen Z, Chen L, Zhang X, Zheng J, Xu H and Liao W (2022) Decreased Functional Connectivities of Low-Degree Level Rich Club Organization and Caudate in Post-stroke Cognitive Impairment Based on Resting-State fMRI and Radiomics Features. Front. Neurosci. 15:796530. doi: 10.3389/fnins.2021.796530
Received
17 October 2021
Accepted
31 December 2021
Published
16 February 2022
Volume
15 - 2021
Edited by
Noorbakhsh Amiri Golilarz, Southern Illinois University Carbondale, United States
Reviewed by
Zhongxiang Ding, Zhejiang University, China; Mahdieh Ghasemi, University of Neyshabur, Iran
Updates
Copyright
© 2022 Miao, Rao, Wang, Fang, Chen, Chen, Zhang, Zheng, Xu and Liao.
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: Haibo Xu, xuhaibo@whu.edu.cnWeijing Liao, weijingliaorehab@163.com
†These authors have contributed equally to this work and share first authorship
This article was submitted to Brain Imaging Methods, a section of the journal Frontiers in Neuroscience
Disclaimer
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.