ORIGINAL RESEARCH article

Front. Neurol., 08 September 2021

Sec. Applied Neuroimaging

Volume 12 - 2021 | https://doi.org/10.3389/fneur.2021.713520

Altered Coupling Between Resting-State Cerebral Blood Flow and Functional Connectivity Strength in Cervical Spondylotic Myelopathy Patients


Wuzeng Wei,&#x;Wuzeng Wei1,2Tao Wang&#x;Tao Wang2Tuersong Abulizi&#x;Tuersong Abulizi2Bing Li,Bing Li1,2Jun Liu,
Jun Liu1,2*
  • 1Department of Joints, Tianjin Hospital, Tianjin University, Tianjin, China
  • 2Clinical College of Orthopedics, Tianjin Medical University, Tianjin, China

Background: Changes in regional neural activity and functional connectivity in cervical spondylotic myelopathy (CSM) patients have been reported. However, resting-state cerebral blood flow (CBF) changes and coupling between CBF and functional connectivity in CSM patients are largely unknown.

Methods: Twenty-seven CSM patients and 24 sex/age-matched healthy participants underwent resting-state functional MRI and arterial spin labeling imaging to compare functional connectivity strength (FCS) and CBF between the two groups. The CBF–FCS coupling of the whole gray matter and specific regions of interest was also compared between the groups.

Results: Compared with healthy individuals, CBF–FCS coupling was significantly lower in CSM patients. The decrease in CBF–FCS coupling in CSM patients was observed in the superior frontal gyrus, bilateral thalamus, and right calcarine cortex, whereas the increase in CBF–FCS coupling was observed in the middle frontal gyrus. Moreover, low CBF and high FCS were observed in sensorimotor cortices and visual cortices, respectively.

Conclusion: In general, neurovascular decoupling at cortical level may be a potential neuropathological mechanism of CSM.

Introduction

Cervical spondylotic myelopathy (CSM) is a common degenerative spinal cord complication that results in significant mortalities (13). Lifestyle and technological advancements changes have paralleled the incidences of CSM (2, 4). Increasing evidence shows that CSM is linked to the altered structure (5) and function (611) of the brain. Functional MRI (fMRI), a blood-oxygen-level-dependent signal, measures the spontaneous neural activity of the brain. Assessment of functional alterations in CSM patients relies on regional homogeneity (ReHo), amplitude of low frequency fluctuation (ALFF), and functional connectivity (FC). These methods have revealed extensive functional reorganizations including regional alterations of specific brain regions (e.g., S1, M1, V1) (7, 1214), circuits (e.g., the thalamus-cortical circuit) (15, 16), subnetworks (e.g., the sensorimotor network and default mode network) (8, 1719), and the whole brain network (8) in CSM patients. Among the aforementioned neuroimaging metrics based on fMRI technique, FC is the most widely used for exploring brain functional reorganization in CSM patients. The inherent simplicity, sensitivity, and ease of interpretation favors the use of seed-based resting-state FC (20). However, this metric is based on the temporal correlation between time-courses extracted from priori-defined seed-regions; thus, it may be unsuitable if the disease pathology is unknown (21). Whole-brain functional connectivity strength (FCS) analysis is another data-driven method, which considers the temporal correlation between each voxel's time-course and all other voxels' time-courses in the brain (22). FCS is referred to as the voxel-level “degree centrality” in graph theory, and brain regions with high FCS are usually considered functional hubs, intricately and functionally connected with the rest of the brain (21, 23). FCS is applied in investigating connectivity changes in various psychological diseases including schizophrenia, ADHD, and AD, among others. (2427). Therefore, FCS may also be suitable for investigating cortical functional changes in CSM patients.

Cerebral blood flow (CBF) is the amount of blood delivered to a specific brain tissue over a given time (2830). Resting-state CBF is strongly related to the brain metabolism (glucose utilization, aerobic glycolysis, and oxygen consumption) (31, 32). Positron emission tomography (PET) and single photon emission computerized tomography (SPECT) have been used in assessing CBF changes in CSM (14). Corresponding assessments show that the spinal cord CBF reflects inflammatory processes in CSM pathology. CBF is also diagnostic and prognostic marker of CSM (3335). SPECT has revealed significantly low brain rs-CBF in the posterior cortical areas of CSM patients (36). Given that PET and SPECT are invasive procedures, they are cumbersome and have low spatial resolution. Contrarily, the arterial spin labeling (ASL) is a non-invasive MR perfusion technique for measuring rs-CBF using an endogenous contrast agent. Therefore, ASL-based CBF can reveal the metabolic state and inflammatory process in the cortices of CSM patients.

Neurovascular coupling (NVC) is an intrinsic brain function which reflects a strong connection between neuronal activity and blood supply (37). In the brain, the bulk of energy is used in driving spontaneous brain activity (38). Based on neurovascular coupling (NVC) hypothesis, brain regions with higher spontaneous regional neural activity tend to have a greater metabolic demand, which increases blood perfusion (39). Functional hubs in human brain display higher FCS and requires higher metabolic demand (4042). Numerous network analyses show that CBF reflects both anatomical and functional connections in the brain (FC, ICA, and FCS-based FC) are associated with (43, 44). Across-voxel CBF–FCS correlation (spatial correlation between FCS and CBF in the whole-brain or gray matter voxels) have been used in characterizing the coupling between the neural activity and hemodynamic response in the brain. This approach identifies altered neurovascular coupling in diseases such as schizophrenia, which cannot be detected using FCS or CBF alone (42, 4547). In recent years, it has been found that altered neurovascular coupling is associated with multiple neuropathology process in the cortex. For instance, abnormal astrocytes or GABA interneurons may disrupt neuronal activity and vascular response (42, 48). Other parameters such as nitric oxide and neuroinflammation levels can also affect neurovascular decoupling.

In the past decades, neuroinflammation, a critical factor affecting neurovascular coupling (49, 50), has been consistently reported in CSM patients (5154). Herein, we hypothesized that compared with healthy individuals, CSM patients would exhibit a significantly lower CBF–FCS coupling compared with healthy controls. As such, we separately explored FCS and CBF changes in specific brain regions of CSM patients. The resting-state fMRI signals and ASL-based cerebral blood flow analyses in 27 CSM patients, and 24 sex/age-matched healthy participants were collected and analyzed. The spatial correlation between CBF and FCS across voxels were then compared between CSM and healthy individuals. The associations between altered neurovascular coupling and the severity of clinical symptoms in patients with cervical compressive myelopathy based on Japanese Orthopedic Association (JOA) scores was also analyzed. CBF and FCS were also compared between groups.

Methods

Participants

The protocol for this study was approved by the ethical review board of Tianjin Hospital, China. All participants consented to participate in this study. In general, 29 CSM patients and 24 healthy individuals were recruited in this study between 2019 and 2020. All CSM patients were right-handed. To be included, CSM patients must have fulfilled the following: (a) with cord compression on cervical spine based on MRI-analysis; (b) with sensorimotor deficit of extremities or bladder/bowel dysfunction; (c) with spinal cord compression; (d) with no history of cervical spine surgery; (e) available for the entire study period; (f) with no stenosis of extracranial vertebral artery and the carotid artery after Doppler ultrasound examination; (g) with no clinical or history of other neurologic, psychiatric, ocular, or systemic diseases like hypertension and diabetes; and (h) with no history of alcohol and substance abuse. The 24 healthy individuals were matched to the CMS patients with regard to age, gender, and education. Specific inclusion criteria included: (a) absence of spine compression, (b) other spinal or brain neurological disorders or systemic disease, and (c) availability for the entire study period. The demographic data of the study participants at baseline are shown in Table 1. Two CSM patients with excessive head motion (more than five times the average scrubbing) were excluded from the study. Therefore, the final analyses were based on 27 CSM patients.

TABLE 1
www.frontiersin.org

Table 1. The demographic data study participants at baseline.

Clinical Assessments

Cognitive aspects of the CSM patients were evaluated using the Montreal Cognitive Assessment (MOCA) and the Mini-Mental State Examination (MMSE). CSM patients underwent a thorough assessment for the severity of cervical compressive myelopathy based on Japanese Orthopedic Association (JOA) scores. Details of clinical assessment for each patient are shown in Supplementary Table 1.

Acquisition of Clinical Data

The 3T fMRI data were acquired using a MAGNETOM Prisma 3T MR scanner (Siemens, Erlangen, Germany) equipped with a 64-channel phase-array head-neck coil. Participants kept their heads still during the scanning process. The head was supported with a sponge pad to minimize unconscious movement. The participants closed their eyes but remained awake, and restrained from specific and strong thoughts. BOLD signals were captured simultaneously using a prototype multi-slice gradient echo-planar imaging (EPI) sequence. The specific test parameters included 30 ms echo time (TE), 800 ms repetition time (TR), 222 × 222 mm field of view (FOV), 74 × 74 matrix, 3 × 3 mm in-plane resolution, flip angle (FA) of 54 degrees, slice thickness of 3 mm, 0 section thickness, 18 slices, transversal slice orientation, bandwidth of 1,690 Hz/pixel, slice acceleration factor of 4, and phase encoding acceleration factor of 2 under parallel acquisition technique (PAT) mode. Overall, 450 images were captured in 6 min.

A high-resolution 3D T1 structural image (two inversion contrast magnetization prepared rapid gradient echo sequences, MP2RAGE) was also captured under the following parameters: 4,000 ms/3.41 ms TR/TE, inversion time (TI1/TI2) of 700 ms/2,110 ms, FA1/FA2 of 4/5°, 256 × 240 matrix, FOV of 256 × 240 mm, 192 slices, 1 × 1 mm in-plane resolution, 1 mm slice thickness, and sagittal slice orientation. The process lasted 6 min and 42 s. Participants closed their eyes, sat in relaxed position and moved minimally, thought about nothing in particular, and stayed awake during the entire process.

Resting-state perfusion imaging was performed using a pseudo-continuous ASL (pcASL) sequence equipped with a 3D fast spin-echo acquisition and background suppression platform (TR/TE of 3,600/21.8 ms; post-label delay of 2,025 ms; spiral in readout of 8 arms with 512 sample points; FA of 180°; FOV of 220 × 220 mm; reconstruction matrix of 128 × 128; slice thickness of 3 mm, no gap; 44 axial slices; eight excitations; and 284 s of acquisition time). The label and control whole-brain image volumes were eight TRs. In total, eight pairs of label and control volumes were captured.

fMRI Data Preprocessing

MR data were preprocessed using the toolbox Data Processing Assistant for rs-fMRI (DPARSF) (http://www.restfmri.net/forum/DPARSF) pipeline. Briefly, 450 volumes were acquired for functional scanning. The first 10 volumes of each functional scans were used for acclimatization of scanning and magnetization stabilization and were excluded. Motion correction was performed to remove the influences of head movement (due to significantly shortened TR, slice-timing correction was not applied). Liner-drift, Friston-24 parameters, mean global signal, white matter signal, and CSF signal covariates were regressed out to minimize non-neural signals. Subsequently, a scrubbing step for high motion time-points was also performed. The standard FD Jenkinson value was set at 0.5. Time-points exceeding the threshold were scrubbed using the cubic spine method (scrubbing time-points before bad time-points: 2; scrubbing time-points after bad time-point). Effects of high-frequency noise were removed using a band pass filter (0.01–0.08 Hz). Functional images were co-registered to structural images and spatially normalized using the Montreal Neurological Institute template. Each voxel was resampled to 3 × 3 × 3 mm3, and then smoothed with a 6 mm full-width-half-maximum isotropic Gaussian kernel before data analysis.

CBF Analysis

To quantify CBF, the ASL images were analyzed using SPM12 (Statistical Parametric Mapping) (http://www.fil.ion.ucl.ac.uk/spm/) and ASL Data Processing Toolbox (ASLtbx) (http://www.cfn.upenn.edu). Pulsed ASL data were processed as follows: first, the ASL images were realigned to correct for head motion. Second, the perfusion difference of images was calculated through sinc subtraction of label/control pairs. Third, relative CBF images of each subject were normalized using the SPM12 software based on the standard Montreal Neurological Institute (MNI) template. Finally, spatial smoothing was performed again using an FWHM of 6 mm. The intervening effect of head-motion, education, age and gender was removed using multiple liner regression analysis.

Functional Connectivity Strength (FCS) Analysis

In this study, we performed voxel-wise FCS analysis. For a given voxel, we (1) first computed Pearson correlation coefficients between the BOLD time-course of this voxel and the rest of voxels within the gray matter mask; (2) the FCS of this voxel is calculated as the number of the coefficients that surpass the threshold of 0.2. In this study, because of the removal of the global signal, there is potential possibility for inducing controversial negative correlations; thus, we restricted our analysis to positive correlations only. Then, this procedure was repeated once for each voxel within the gray matter mask. Therefore, for each subject, we obtained an FCS map and it was spatially smoothed with a 6 × 6 × 6 mm FWHM Gaussian kernel. To remove potential influences of head-motion, education, age, and gender, we set these variables as covariates and regress out these effects by multiple liner regression.

Whole Gray Matter CBF–FCS Coupling Analysis

For subject-level analysis, both CBF and FCS maps were normalized into z-scores by dividing the SD of global values within the gray matter mask with the difference between the mean and individual global values within the gray matter mask. The coupling between FCS and CBF was based on the spatial correlation between CBF and FCS maps across all voxels within the gray matter mask.

For group-level analysis, the average of both CBF and FCS maps of each group were first obtained. The CBF and FCS maps were then converted into z-scores by subtracting the mean of individual global values and dividing the results by the SD of global values within the gray matter mask. The correlation between CBF and FCS maps of each group were then calculated. To reveal group-level CBF–FCS coupling difference between CSM and HC, the correlation coefficients of the two groups were compared. This analysis was performed using an online platform (http://comparingcorrelations.org) as previously described (55).

Subsequently, the spatial correlation between the CBF and FCS map for each subject was calculated to obtain a subject-level whole-brain CBF–FCS coupling. The difference in subject-level CBF–FCS coupling between CSM and HC was analyzed using a two-sample t-test.

Region-Wise CBF–FCS Coupling Analysis

For region-wise analysis, both CBF and FCS maps were converted into z-scores by subtracting the mean of a group from global values within the gray matter mask and dividing the difference by the corresponding SD of the group. For a given subject, the coupling between FCS and CBF was the spatial correlation between CBF and FCS maps across all voxels within each subregion, based on the Anatomical Automatic Labeling (AAL) template. This generated 116 coefficients representing the CBF–FCS coupling for a given subregion. The difference in CBF and FCS between CSM patients and healthy individuals for a subregion were analyzed using two sample t-tests. The resultant p-values were corrected by the False Discovery Rate (FDR) method.

Region-Wise CBF and FCS Analysis

The CBF and FCS maps were first normalized into z-scores, then the mean values for each sub-region defined by AAL template were obtained and compared between the two groups. The resultant p-values were corrected by FDR for multiple comparison correction.

Voxel-Wise Comparisons in CBF and FCS

To better explore the differences between CSM patients and healthy individuals, we compared the differences in CBF and FCS between the two groups in a voxel-wise manner while controlling for head motion, education, age, and sex. Multiple voxel-wise comparisons were also performed by FDR at statistical significance of corrected p < 0.05.

ROI-Wise Analysis

Although the multiple comparison correction method (FDR correction) was used in our previous region-wise analyses, this method may overlook important differences within brain regions. Therefore, to better explore functional or structural alterations in CSM patients, 12 regions of interest (ROIs) were selected for further analysis. They included the left and right precentral gyrus, the left and right postcentral gyrus, the left and right supplementary motor area, the left and right calcarine gyrus, the left and right precuneus, and the left and right thalamus (7, 12, 19, 56, 57). The ROIs were described based on the corresponding brain areas defined in the AAL. The association between CBF and FCS between voxels in a given region was also analyzed. The mean CBF and FCS between CSM patients and healthy individuals along age, gender, education years, and head motion were also compared.

The Relationship Between Altered Brain Function and Specific Clinical Parameters

The possible relationship between the whole-brain level CBF and FCS in the whole gray matter voxels was analyzed using Cox regression analysis. The relationship between CBF as well as FCS and several clinical parameters including preoperative/postoperative JOA scores, JOA recovery, as well as MoCA and MMSE scores was also analyzed. The relationship between CBF as well as FCS and the aforementioned clinical parameters in a given brain region was also evaluated. Finally, the relationship between CBF as well as FCS in ROI and the afore-mentioned clinical parameters was also performed.

Results

Whole Gray Matter CBF–FCS Coupling Changes in CSM

Compared with healthy controls, CSM patients exhibited significant decreased CBF–FCS coupling for both subject- and group-level analysis. In group-level analysis, the mean CBF and FCS of each group (i.e., mean maps across subjects within each group) were obtained. The correlation coefficient between group-mean FCS map and group-mean CBF map in CSM was 0.28; the correlation coefficient between group-mean FCS map and group-mean CBF map in HC was 0.33. These two correlation coefficients were compared, and the resultant z-value was −46.0489, p < 0.0001 [Dunn and Clark's method (58)], z-value was −45.6988, p < 0.0001 [Steiger's method (59)], and z-value was −46.0489, p < 0.0001 [Raghunathan's method (60)]. In subject-level analysis, the correlation coefficients between the CBF map and FCS map across whole gray matter voxels were compared between CSM patients and healthy controls. We observed significantly decreased CBF–FCS coupling in CSM patients (t-value = 2.109, df = 49, p = 0.0401) (Figure 1).

FIGURE 1
www.frontiersin.org

Figure 1. The CBF–FCS coupling changes in the whole gray matter of cervical spondylotic myelopathy (CSM) patients. Scatter plots show the spatial across voxels correlations between CBF and FCS in CSM patients (red) and healthy controls (blue). Although there is significant correlation between CBF and FCS in both CSM patients and healthy individuals, compared to healthy individuals CSM patients display significantly weak CBF–FCS coupling. Error bars represent the SD. CBF, cerebral blood flow; FCS, functional connectivity strength.

Region-Wise CBF–FCS Coupling, CBF and FCS Changes in CSM Patients

Compared with healthy controls, CSM patients exhibited significantly low CBF–FCS coupling in the superior frontal gyrus (t-value = −3.51, q-value = 0.04, FDR corrected), but substantially high CBF–FCS coupling in the middle frontal gyrus (t-value = 3.56, q-value = 0.03, FDR corrected) (Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2. Region-wise CBF–FCS coupling changes based on anatomical automatic labeling (AAL) in cervical spondylotic myelopathy (CSM) patients. Significant CBF-FCS correlation across voxels within the middle frontal gyrus and superior frontal gyrus were observed between CSM patients and healthy individuals. Error bars represent the SD. FDR-corrected p-values are shown above the bars. CBF, cerebral blood flow; FCS, functional connectivity strength.

Compared with healthy controls, the FCS of CSM patients was significantly low in the left precentral gyrus (t-value = −3.54, q-value = 0.03, FDR corrected) and left postcentral gyrus (t-value = −3.50, q-value = 0.04, FDR corrected). Also, the FCS of CSM patients was substantially high in the right hippocampus (t-value = 3.62, q-value = 0.02, FDR corrected), right parahippocampus (t-value = 3.51, q-value = 0.04, FDR corrected), right amygdala (t-value = 3.49, q-value = 0.04, FDR corrected), right pallidum (t-value = 3.61, q-value = 0.02, FDR corrected), right thalamus (t-value = 3.55, q-value = 0.03, FDR corrected), and left thalamus (t-value = 3.63, q-value = 0.02, FDR corrected) (Figure 3).

FIGURE 3
www.frontiersin.org

Figure 3. Region-wise CBF–FCS coupling changes in cervical spondylotic myelopathy (CSM) patients based on AAL. FCS was significantly different between CSM patients and healthy individuals. Error bars represent the SD, FDR-corrected p-values are shown above bars. CBF, cerebral blood flow; FCS, functional connectivity strength; Precentral_L, Left precentral gyrus; Hippocampus_R, Right hippocampus; ParaHippocampal_R, Right parahippocampus; Amygdala_R, Right Amygdala; Postcentral_L, Left postcentral gyrus; Pallidum_R, Right pallidum; Thalamus_L, Left Thalamus; Thalamus_R, Right thalamus.

However, there were no significant differences in region-wise CBF between CSM patients and healthy individuals, possibility because of the small sample size used in this study and the relatively high strict multiple comparison correction method used.

Voxel-Wise CBF and FCS Changes in CSM Patients

Compared with healthy individuals, the FCS of CSM patients was significantly high in the left thalamus, right thalamus, right anterior cingulate gyrus, and right hippocampus. Contrarily, the FCS of CSM patients was significantly low in the left postcentral gyrus. Details of these clusters are shown in Table 2. The spatial distribution of the clusters is shown in Figure 4.

TABLE 2
www.frontiersin.org

Table 2. Difference in voxel-wise FCS between CSM patients and healthy individuals.

FIGURE 4
www.frontiersin.org

Figure 4. Voxel-wise functional connectivity strength (FCS) changes in cervical spondylotic myelopathy (CSM) patients and healthy patients. Significant difference in FCS was found between the two groups. Scatter plots for the relationship between altered FCS and clinical measures including Montreal cognitive assessments (right), and preoperative Japanese Orthopedics Association scores (Left) across voxels. FCS, functional connectivity strength; ACC, anterior cingulate gyrus.

ROI-Wise CBF–FCS Coupling, CBF and FCS Changes in CSM Patients

Voxel-wise and region-wise analysis is highly robust and stringent, given the multiple comparison corrections. Therefore, based on functional alterations in CSM patients, several analyses in ROIs were performed.

Compared with healthy individuals, the CBF of CSM patients was significantly low in the left precentral gyrus (t-value = −2.65, p = 0.01, uncorrected), right calcarine gyrus (t-value = −2.05, p = 0.04, uncorrected), left postcentral gyrus (t-value = −2.27, p = 0.02, uncorrected), left precuneus (t-value = −2.93, p = 0.005, uncorrected), and right precuneus (t-value = −2.47, p = 0.01, uncorrected) (Figure 5).

FIGURE 5
www.frontiersin.org

Figure 5. CBF–FCS coupling changes in ROI in cervical spondylotic myelopathy (CSM) patients based on anatomical automatic labeling. We found significant difference in CBF–FCS coupling, as well as CBF rate and FCS. Error bars represent the SD. P-values (uncorrected) are shown above bars. CBF, cerebral blood flow; FCS, functional connectivity strength; ROI, region of interests; Precentral_L, Left precentral gyrus; Precentral_R, Right precentral gyrus; Supp_Motor_Area_L, Left supplementary motor area; Supp_Motor_Area_R, Right supplementary motor area; Calcarine_L, Left calcarine cortex; Calcarine_R, Right calcarine cortex; Postcentral_L, Left postcentral gyrus; Postcentral_R, Right postcentral gyrus; Precuneus_L, Left precuneus cortex; Precuneus_R, Right precuneus cortex; Thalamus_L, Left Thalamus; Thalamus_R, Right thalamus.

The FCS of CSM patients was also significantly low in the left precentral gyrus (t-value = −3.54, p < 0.001, uncorrected), right precentral gyrus (t-value = −2.35, p = 0.02, uncorrected), and left postcentral gyrus (t-value = −3.41, p < 0.001, uncorrected) but substantially high in the left (t-value = 3.63, p < 0.001, uncorrected), and right thalamus (t-value = 3.55, p < 0.001, uncorrected) (Figure 5).

CBF–FCS coupling was also low in the right calcarine gyrus (t-value = −2.04, p = 0.04, uncorrected), left thalamus (t-value = −2.20, p = 0.03, uncorrected), and right thalamus (t-value = −2.23, p = 0.03, uncorrected) (Figure 5).

Correlation Between Brain Alterations and Clinical Measures

Subject-level analysis revealed no correlation between CBF and FCS levels in the whole gray matter of CSM patients (Supplementary Table 2).

Region-wise analysis revealed a strong positive correlation between CBF as well as FCS in the left middle frontal gyrus and postoperative JOA scores (severity of clinical symptoms in patients with cervical compressive myelopathy) in CSM patients (r = 0.50, p < 0.05, uncorrected) (Supplementary Table 3). The FCS in the left precentral gyrus also positively correlated with preoperative JOA scores (r = 0.53, p < 0.05) (Supplementary Table 4).

Voxel-wise analysis revealed that the FCS in the left precentral gyrus positively correlated with the preoperative JOA scores (r = 0.63, p < 0.001), whereas the FCS in the right hippocampus positively correlated with MoCA scores (r = 0.45, p < 0.05) (Figure 4).

ROI-wise analysis on its part revealed a strong negative correlation between postoperative JOA scores and CBF as well as FCS and in the left calcarine gyrus (r = −0.37, p < 0.05). A negative correlation was also observed between preoperative JOA scores and CBF as well as FCS in the right calcarine gyrus (r = −0.37, p < 0.05). A strong negative correlation was observed between preoperative JOA scores and CBF as well as FCS in the left thalamus (r = −0.48, p < 0.05). We also found a negative correlation between postoperative JOA scores and CBF as well as FCS in the right thalamus (r = −0.38, p < 0.05) (Supplementary Table 5). Preoperative JOA score positively correlated with FCS in the left precentral gyrus. We observed a negative correlation between the FCS in the left calcarine and preoperative JOA scores (r = −0.37, p < 0.05) as well as JOA recovery (r = 0.40, p < 0.05) following decompression surgery. The FCS in the left precuneus negatively correlated with preoperative JOA scores (r = −0.54, p < 0.05) as well as JOA recovery (r = 0.50, p < 0.05), whereas the FCS in the right precuneus negatively correlated with preoperative JOA scores (r = −0.61, p < 0.05) as well as JOA recovery (r = 0.53, p < 0.05) (Supplementary Table 6). In this section, no multiple comparison correction method was performed because the relative strict multiple comparison method may cover up the mild association between brain variables and clinical measures within the selected ROIs.

Discussion

To the best of our knowledge, this is the first study investigating CBF–FCS coupling changes in CSM patients using a combination of BOLD and ASL techniques. We found that, compared with healthy individuals, the CBF–FCS coupling in the whole gray matter was significantly low. CSM patients also exhibited altered CBF–FCS coupling in the superior/mid frontal gyrus. Moreover, the FCS in the pre/postcentral gyrus was substantially low in CSM patients, in contrast with rCBF in the thalamus, hippocampus, and visual cortices, which were all substantially high. These findings deepen our understanding on the neural pathological characteristics of CSM. We used blood oxygenation level-dependent (BOLD) signals in exploring functional alterations in the brain under certain disease complications such as Alzheimer's disease (AD), Parkinson's disease (PD), epilepsy, attention deficit hyperactivity disorder (ADHD), and mood disorders.

We observed a significant whole-brain wise across-voxel correlation between rCBF and FCS in healthy adults, consistent with previous findings (27, 45). Proper neurovascular functioning is necessary for normal coupling between rCBF and FCS (both neuronal and vascular components). Astrocytes maintain normal coupling between neuronal and vascular components by coordinating the activity and responses mediated by these components (48). Although we observed a significant across-voxel correlation between rCBF and FCS in CSM patients, it was substantially lower than that of healthy patients, implying decoupling between rCBF and FCS in CSM patients. This neurovascular decoupling may have resulted from prolonged neuroinflammation caused by myelopathy (6163), in which the abnormally functioning astrocytes in CSM patients may have reduced the regional neural activity and blood supply (62, 64, 65). Neuroinflammation would have also directly affected the vascular component, leading to neurovascular decoupling in CSM patients. However, we did not find a significant difference in rCBF between CSM patients and healthy individuals. Therefore, neurovascular decoupling in CSM is caused by abnormal functioning of the astrocyte.

We observed low across-voxel correlation between rCBF and FCS in the superior frontal gyrus of CSM patients. The superior frontal gyrus regulates motor functions. Previous neuroimaging studies have demonstrated that CSM patients display altered neural activity in the superior frontal gyrus. Takenaka et al. observed significantly high functional connectivity between superior frontal and lingue gyrus in (healthy individuals), and surgery substantially improved this phenomenon in CSM patients (66). In their subsequent study using the same cohort, Takenaka et al. found that decompression surgery increased ALFF in the superior frontal gyrus of CSM patients (9). Herein, we found a high across-voxel correlation between rCBF and FCS in the middle frontal gyrus of (CSM patients), a brain region that participates in processing chronic pain (67, 68). The high neurovascular coupling may have resulted from the chronic pain in CSM patients. Taken together, CSM patients display altered superior frontal and middle frontal gyrus function. Meanwhile, the neurovascular coupling (measured by cross-voxel correlation between rCBF and FCS) provides complementary information on the pathological changes in CSM.

FCS levels reflect the contribution of each voxel in information transmission in the whole brain. Information transmission efficiency is directly proportional to HFCS levels. Low FCS implies local brain damage (reduction in gray matter and alterations in white matter) (69, 70). Herein, we found FCS within sensorimotor cortices including pre and postcentral gyrus, and the FCS within left precentral gyrus was significantly low in CSM patients, and positively correlated with preoperative JOA scores. This implies that regional impairments are directly proportional to the severity of clinical symptoms. We speculate the functional impairments result from the structural damage in sensorimotor cortices in CSM patients. Converging evidence shows that CSM patients have little gray matter (7173) and low metabolism in the sensorimotor cortices (14, 74). Therefore, the structural changes may result from the FCS changes in the precentral and postcentral gyrus. Moreover, we also observed a significant increase in FCS within the thalamus and hippocampus of CSM patients, and high FCS in the hippocampus positively correlated with MoCA scores. High global efficiency of the default mode network positively correlates with MoCA scores in CSM patients. This implies that the high global efficiency in the default mode network compensates for the local impairment (measured by signal variability) in the default mode network (19). As the hippocampus has been shown to be an important brain region participating in many cognitive functions, the high FCS in the hippocampus was also associated with compensatory or adaptive changes in CSM, necessary for maintaining normal cognitive function in these individuals.

What causes functional changes in the brain of CSM patients is largely unknown. However, it is thought that CSM patients develop stenosis of transverse foramen that compresses vertebral artery, altering cortical blood supply. However, in our current study, altered CBF in CSM patients was only observed in a few regions. We found no significant voxel-wise or region-wise difference in CBF between CSM patients and healthy individuals. However, we observed altered CBF and FCS coupling in CSM patients, consistent with the idea that functional impairments in CSM patients is caused by neuroinflammation.

Regarding limitations, first, the CSM patients included in our study received a long conservative treatment before surgery (NSAID), which could have influenced the observed outcomes. Future studies using treatment-naïve CSM patients are necessary to validate our findings. Second, both CBF and FCS are indirect measures of vascular and neuronal function; thus, our findings may not be conclusive. Third, postoperative fMRI assessment was not performed due to the possible artifacts of the surgical implants and possible heating of these materials. Means of conducting safe postoperative fMRI are necessary to allow this process in patients with surgical implants. Lastly, the sample size of our current study is relatively small. Future studies with a larger sample size are needed to further confirm our conclusion.

In conclusion, we observed an altered neurovascular coupling in CSM by combining BOLD and ASL techniques. Specifically, we found a decreased across-voxel correlation between rCBF and FCS in the superior frontal gyrus involved in motor control and motor planning, and an increased across-voxel correlation between rCBF and FCS in the middle frontal gyrus implicated in cognitive function. These findings suggest that the neurovascular decoupling in the brain may be a potential neural mechanism involved in the pathophysiology of CSM.

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 TianJin Hospital. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

JL designed the study. WW wrote the manuscript. TW, BL, and TA collected the data. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the Tianjin Health Commission Foundation (Grant Number: ZC20192), Tianjin Biomedical Engineering Technology Key Project fund (Grant Number: 18ZXSGSY00010), Tianjin Science Technology Key Project fund (Grant Number: 18YFZCSY00890), and Natural Science Foundation of Tianjin fund (Grant Number: 20JCQNJC01170).

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.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2021.713520/full#supplementary-material

References

1. Akter F, Kotter M. Pathobiology of degenerative cervical myelopathy. Neurosurg Clin N Am. (2018) 29:13–9. doi: 10.1016/j.nec.2017.09.015

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Badhiwala JH, Ahuja CS, Akbar MA, Witiw CD, Nassiri F, Furlan JC, et al. Degenerative cervical myelopathy - update and future directions. Nat Rev Neurol. (2020) 16:108–24. doi: 10.1038/s41582-019-0303-0

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Davies BM, Mowforth OD, Smith EK, Kotter MR. Degenerative cervical myelopathy. BMJ. (2018) 360:k186. doi: 10.1136/bmj.k186

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Lv Y, Tian W, Chen D, Liu Y, Wang L, Duan F. The prevalence and associated factors of symptomatic cervical Spondylosis in Chinese adults: a community-based cross-sectional study. BMC Musculoskelet Disord. (2018) 19:325. doi: 10.1186/s12891-018-2234-0

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Wen CY, Cui JL, Mak KC, Luk KD, Hu Y. Diffusion tensor imaging of somatosensory tract in cervical spondylotic myelopathy and its link with electrophysiological evaluation. Spine J. (2014) 14:1493–500. doi: 10.1016/j.spinee.2013.08.052

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Bhagavatula ID, Shukla D, Sadashiva N, Saligoudar P, Prasad C, Bhat DI. Functional cortical reorganization in cases of cervical spondylotic myelopathy and changes associated with surgery. Neurosurg Focus. (2016) 40:E2. doi: 10.3171/2016.3.FOCUS1635

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Chen Z, Zhao R, Wang Q, Yu C, Li F, Liang M, et al. Functional connectivity changes of the visual cortex in the cervical spondylotic myelopathy patients: a resting-state fMRI study. Spine. (2019) 45:E272–9. doi: 10.1097/BRS.0000000000003245

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Kuang C, Zha Y. Abnormal intrinsic functional activity in patients with cervical spondylotic myelopathy: a resting-state fMRI study. Neuropsychiatr Dis Treat. (2019) 15:2371–83. doi: 10.2147/NDT.S209952

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Takenaka S, Kan S, Seymour B, Makino T, Sakai Y, Kushioka J, et al. Resting-state amplitude of low-frequency fluctuation is a potentially useful prognostic functional biomarker in cervical myelopathy. Clin Orthop Relat Res. (2020) 478:1667–80. doi: 10.1097/CORR.0000000000001157

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Zhou FQ, Tan YM, Wu L, Zhuang Y, He LC, Gong HH. Intrinsic functional plasticity of the sensory-motor network in patients with cervical spondylotic myelopathy. Sci Rep. (2015) 5:9975. doi: 10.1038/srep09975

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Liu X, Qian W, Jin R, Li X, Luk KD, Wu EX, et al. Amplitude of low frequency fluctuation (ALFF) in the cervical spinal cord with stenosis: a resting state fMRI study. PLoS ONE. (2016) 11:e0167279. doi: 10.1371/journal.pone.0167279

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Chen Z, Wang Q, Liang M, Zhao R, Zhu J, Xiong W, et al. Visual cortex neural activity alteration in cervical spondylotic myelopathy patients: a resting-state fMRI study. Neuroradiology. (2018) 60:921–32. doi: 10.1007/s00234-018-2061-x

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Ryan K, Goncalves S, Bartha R, Duggal N. Motor network recovery in patients with chronic spinal cord compression: a longitudinal study following decompression surgery. J Neurosurg Spine. (2018) 28:379–88. doi: 10.3171/2017.7.SPINE1768

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Aleksanderek I, McGregor SMK, Stevens TK, Goncalves S, Bartha R, Duggal N. Cervical spondylotic myelopathy: metabolite changes in the primary motor cortex after surgery. Radiology. (2017) 282:817–25. doi: 10.1148/radiol.2016152083

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Peng X, Tan Y, He L, Ou Y. Alterations of functional connectivity between thalamus and cortex before and after decompression in cervical spondylotic myelopathy patients: a resting-state functional MRI study. Neuroreport. (2020) 31:365–71. doi: 10.1097/WNR.0000000000001346

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Zhou F, Wu L, Liu X, Gong H, Luk KD, Hu Y. Characterizing thalamocortical disturbances in cervical spondylotic myelopathy: revealed by functional connectivity under two slow frequency bands. PLoS ONE. (2015) 10:e0125913. doi: 10.1371/journal.pone.0125913

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Dong Y, Holly LT, Albistegui-Dubois R, Yan X, Marehbian J, Newton JM, et al. Compensatory cerebral adaptations before and evolving changes after surgical decompression in cervical spondylotic myelopathy. J Neurosurg Spine. (2008) 9:538–51. doi: 10.3171/SPI.2008.10.0831

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Tan Y, Zhou F, Wu L, Liu Z, Zeng X, Gong H, et al. Alteration of regional homogeneity within the sensorimotor network after spinal cord decompression in cervical spondylotic myelopathy: a resting-state fMRI study. Biomed Res Int. (2015) 2015:647958. doi: 10.1155/2015/647958

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Zhao R, Su Q, Chen Z, Sun H, Liang M, Xue Y. Neural correlates of cognitive dysfunctions in cervical spondylotic myelopathy patients: a resting-state fMRI study. Front Neurol. (2020) 11:596795. doi: 10.3389/fneur.2020.596795

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Fox MD, Raichle ME. Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nat Rev Neurosci. (2007) 8:700–11. doi: 10.1038/nrn2201

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Nair A, Keown CL, Datko M, Shih P, Keehn B, Müller RA. Impact of methodological variables on functional connectivity findings in autism spectrum disorders. Hum Brain Map. (2014) 35:4035–48. doi: 10.1002/hbm.22456

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Buckner RL, Sepulcre J, Talukdar T, Krienen FM, Liu H, Hedden T, et al. Cortical hubs revealed by intrinsic functional connectivity: mapping, assessment of stability, and relation to Alzheimer's disease. J Neurosci. (2009) 29:1860–73. doi: 10.1523/JNEUROSCI.5062-08.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Wang L, Xia M, Li K, Zeng Y, Su Y, Dai W, et al. The effects of antidepressant treatment on resting-state functional brain networks in patients with major depressive disorder. Hum Brain Map. (2015) 36:768–78. doi: 10.1002/hbm.22663

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Pijnenburg R, Scholtens LH, Mantini D, Vanduffel W, Barrett LF, van den Heuvel MP. Biological characteristics of connection-wise resting-state functional connectivity strength. Cerebral Cortex. (2019) 29:4646–53. doi: 10.1093/cercor/bhy342

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Wang S, Zhan Y, Zhang Y, Lv L, Wu R, Zhao J, et al. Abnormal functional connectivity strength in patients with adolescent-onset schizophrenia: a resting-state fMRI study. Eur Child Adolesc Psychiatry. (2017) 26:839–45. doi: 10.1007/s00787-017-0958-2

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Yang M, He H. The effects of music intervention on functional connectivity strength of the brain in schizophrenia. Neural Plast. (2018) 2018:2821832. doi: 10.1155/2018/2821832

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Zuo L, Zhou Y, Wang S, Wang B, Gu H, Chen J. Abnormal brain functional connectivity strength in the overactive bladder syndrome: a resting-state fMRI study. Neural Plast. (2019) 131:64–70. doi: 10.1016/j.urology.2019.05.019

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Golay X, Hendrikse J, Lim TC. Perfusion imaging using arterial spin labeling. Top Magn Reson Imaging. (2004) 15:10–27. doi: 10.1097/00002142-200402000-00003

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Kindler J, Jann K, Homan P, Hauf M, Walther S, Strik W, et al. Static and dynamic characteristics of cerebral blood flow during the resting state in schizophrenia. Schizophr Bull. (2015) 41:163–70. doi: 10.1093/schbul/sbt180

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Williams DS, Detre JA, Leigh JS, Koretsky AP. Magnetic resonance imaging of perfusion using spin inversion of arterial water. Proc Natl Acad Sci USA. (1992) 89:212–6. doi: 10.1073/pnas.89.1.212

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Raichle ME, MacLeod AM, Snyder AZ, Powers WJ, Gusnard DA, Shulman GL, et al. default mode of brain function. Proc Natl Acad Sci USA. (2001) 98:676–82. doi: 10.1073/pnas.98.2.676

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Vaishnavi SN, Vlassenko AG, Rundle MM, Snyder AZ, Mintun MA, Raichle ME. Regional aerobic glycolysis in the human brain. Proc Natl Acad Sci USA. (2010) 107:17757–62. doi: 10.1073/pnas.1010459107

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Eicker SO, Langen KJ, Galldiks N, Stoffels G, Herdmann J, Steiger HJ, et al. Clinical value of 2-deoxy-[18F]fluoro-D-glucose positron emission tomography in patients with cervical spondylotic myelopathy. Neurosurg Focus. (2013) 35:E2. doi: 10.3171/2013.3.FOCUS1379

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Floeth FW, Stoffels G, Herdmann J, Eicker S, Galldiks N, Rhee S, et al. Prognostic value of 18F-FDG PET in monosegmental stenosis and myelopathy of the cervical spinal cord. J Nucl Med. (2011) 52:1385–91. doi: 10.2967/jnumed.111.091801

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Sakushima K, Yabe I, Shiga T, Yashima-Yamada M, Tsuji-Akimoto S, Terae S, et al. FDG-PET SUV can distinguish between spinal sarcoidosis and myelopathy with canal stenosis. J Neurol. (2011) 258:227–30. doi: 10.1007/s00415-010-5729-7

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Hoshimaru M. Neuropsychological improvement in patients with cervical spondylotic myelopathy after posterior decompression surgery. Neurol Med Chir. (2010) 50:554–9. doi: 10.2176/nmc.50.554

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Segarra M, Aburto MR, Hefendehl J, Acker-Palmer A. Neurovascular interactions in the nervous system. Ann Rev Cell Dev Biol. (2019) 35:615–35. doi: 10.1146/annurev-cellbio-100818-125142

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Phillips AA, Chan FH, Zheng MM, Krassioukov AV, Ainslie PN. Neurovascular coupling in humans: physiology, methodological advances and clinical implications. J Cereb Blood Flow Metab. (2016) 36:647–64. doi: 10.1177/0271678X15617954

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Masamoto K, Hirase H, Yamada K, Kanno I. Neurovascular coupling-What next? Prog Brain Res. (2016) 225:269–72. doi: 10.1016/bs.pbr.2016.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Jin M, Wang L, Wang H, Han X, Diao Z, Guo W, et al. Disturbed neurovascular coupling in hemodialysis patients. PeerJ. (2020) 8:e8989. doi: 10.7717/peerj.8989

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Junejo RT, Braz ID, Lucas SJ, van Lieshout JJ, Phillips AA, Lip GY, et al. Neurovascular coupling and cerebral autoregulation in atrial fibrillation. J Cerebral Blood Flow Metab. (2020) 40:1647–57. doi: 10.1177/0271678X19870770

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Guo X, Zhu J, Zhang N, Zhang L, Qi Y, Cai H, et al. Altered neurovascular coupling in neuromyelitis optica. PeerJ. (2019) 40:976–86. doi: 10.1002/hbm.24426

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Hosford PS, Gourine AV. What is the key mediator of the neurovascular coupling response? Neurosci Biobehav Rev. (2019) 96:174–81. doi: 10.1016/j.neubiorev.2018.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Iadecola C, Gottesman RF. Neurovascular and cognitive dysfunction in hypertension. Hum Brain Mapp. (2019) 124:1025–44. doi: 10.1161/CIRCRESAHA.118.313260

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Zhu J, Zhuo C, Xu L, Liu F, Qin W, Yu C. Altered coupling between resting-state cerebral blood flow and functional connectivity in schizophrenia. Schizophr Bull. (2017) 43:1363–74. doi: 10.1093/schbul/sbx051

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Zhu J, Zhang Y, Zhang B, Yang Y, Wang Y, Zhang C, et al. Abnormal coupling among spontaneous brain activity metrics and cognitive deficits in major depressive disorder. J Affect Disord. (2019) 252:74–83. doi: 10.1016/j.jad.2019.04.030

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Zhu J, Zhu DM, Qian Y, Li X, Yu Y. Altered spatial and temporal concordance among intrinsic brain activity measures in schizophrenia. J Psychiatr Res. (2018) 106:91–8. doi: 10.1016/j.jpsychires.2018.09.015

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Liu LR, Liu JC, Bao JS, Bai QQ, Wang GQ. Interaction of microglia and astrocytes in the neurovascular unit. Front Immunol. (2020) 11:1024. doi: 10.3389/fimmu.2020.01024

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Nedeljkovic N. Complex regulation of ecto-5'-nucleotidase/CD73 and A(2A)R-mediated adenosine signaling at neurovascular unit: a link between acute and chronic neuroinflammation. Pharmacol Res. (2019) 144:99–115. doi: 10.1016/j.phrs.2019.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Price BR, Johnson LA, Norris CM. Reactive astrocytes: the nexus of pathological and clinical hallmarks of Alzheimer's disease. Ageing Res Rev. (2021) 68:101335. doi: 10.1016/j.arr.2021.101335

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Gohmann RF, Blume C, Zvyagintsev M, Mainz V, Clusmann H, Wiesmann M, et al. Cervical spondylotic myelopathy: changes of fractional anisotropy in the spinal cord and magnetic resonance spectroscopy of the primary motor cortex in relation to clinical symptoms and their duration. Eur J Radiol. (2019) 116:55–60. doi: 10.1016/j.ejrad.2019.04.009

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Hirai T, Uchida K, Nakajima H, Guerrero AR, Takeura N, Watanabe S, et al. The prevalence and phenotype of activated microglia/macrophages within the spinal cord of the hyperostotic mouse (twy/twy) changes in response to chronic progressive spinal cord compression: implications for human cervical compressive myelopathy. PLoS ONE. (2013) 8:e64528. doi: 10.1371/journal.pone.0064528

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Laliberte AM, Karadimas SK, Vidal PM, Satkunendrarajah K, Fehlings MG. Mir21 modulates inflammation and sensorimotor deficits in cervical myelopathy: data from humans and animal models. Brain Commun. (2021) 3:fcaa234. doi: 10.1093/braincomms/fcaa234

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Liu W, Nair G, Vuolo L, Bakshi A, Massoud R, Reich DS, et al. In vivo imaging of spinal cord atrophy in neuroinflammatory diseases. Ann Neurol. (2014) 76:370–8. doi: 10.1002/ana.24213

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Diedenhofen B, Musch J. cocor: a comprehensive solution for the statistical comparison of correlations. PLoS ONE. (2015) 10:e0121945. doi: 10.1371/journal.pone.0121945

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Sun LQ, Li M, Li YM. Predictors for surgical outcome of laminoplasty for cervical spondylotic myelopathy. World Neurosurg. (2016) 94:89–96. doi: 10.1016/j.wneu.2016.06.092

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Zhou F, Gong H, Liu X, Wu L, Luk KD, Hu Y. Increased low-frequency oscillation amplitude of sensorimotor cortex associated with the severity of structural impairment in cervical myelopathy. PLoS ONE. (2014) 9:e104442. doi: 10.1371/journal.pone.0104442

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Dunn OJ, Clark V. Comparison of tests of the equality of dependent correlation coefficients. J Am Stat Assoc. (1971) 66:904–8. doi: 10.1080/01621459.1971.10482369

CrossRef Full Text | Google Scholar

59. Steiger JH. Tests for comparing elements of a correlation matrix. Psychol Bull. (1980) 87:245–51. doi: 10.1037/0033-2909.87.2.245

CrossRef Full Text | Google Scholar

60. Raghunathan TE, Rosenthal R, Rubin DB. Comparing correlated but nonoverlapping correlations. Psychol Methods. (1996) 1:178–83. doi: 10.1037/1082-989X.1.2.178

CrossRef Full Text | Google Scholar

61. Hashimoto K, Kobatake Y, Asahina R, Yamato O, Islam MS, Sakai H, et al. Up-regulated inflammatory signatures of the spinal cord in canine degenerative myelopathy. Res Vet Sci. (2021) 135:442–9. doi: 10.1016/j.rvsc.2020.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Lien BV, Tuszynski MH, Lu P. Astrocytes migrate from human neural stem cell grafts and functionally integrate into the injured rat spinal cord. Exp Neurol. (2019) 314:46–57. doi: 10.1016/j.expneurol.2019.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Takahashi H, Aoki Y, Nakajima A, Sonobe M, Terajima F, Saito M, et al. Axonal damage is remarkable in patients with acutely worsening symptoms of compression myelopathy: biomarkers in cerebrospinal fluid samples. Eur Spine J. (2018) 27:1824–30. doi: 10.1007/s00586-018-5549-5

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Yovich JV, Gould DH, LeCouteur RA. Chronic cervical compressive myelopathy in horses: patterns of astrocytosis in the spinal cord. Aust Vet J. (1991) 68:334–7. doi: 10.1111/j.1751-0813.1991.tb03092.x

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Yovich JV, le Couteur RA, Gould DH. Chronic cervical compressive myelopathy in horses: clinical correlations with spinal cord alterations. Austr Vet J. (1991) 68:326–34. doi: 10.1111/j.1751-0813.1991.tb03091.x

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Takenaka S, Kan S. Towards prognostic functional brain biomarkers for cervical myelopathy: a resting-state fMRI study. Sci Rep. (2019) 9:10456. doi: 10.1038/s41598-019-46859-5

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Huang X, Chen J. Impaired frontal-parietal control network in chronic prostatitis/chronic pelvic pain syndrome revealed by graph theoretical analysis: a DTI study. Eur J Neurosci. (2021) 53:1060–71. doi: 10.1111/ejn.14962

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Liu S, Gong Q, Liu T, Lu C, Qin Z, Cui H, et al. The neural correlates of pain-related fear: a meta-analysis comparing fear conditioning studies using painful and non-painful stimuli. Eur J Neurosci. (2020) 119:52–65. doi: 10.1016/j.neubiorev.2020.09.016

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Fitzsimmons J, Kubicki M, Shenton ME. Review of functional and anatomical brain connectivity findings in schizophrenia. Curr Opin Psychiatry. (2013) 26:172–87. doi: 10.1097/YCO.0b013e32835d9e6a

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Fornito A, Harrison BJ, Goodby E, Dean A, Ooi C, Nathan PJ, et al. Functional dysconnectivity of corticostriatal circuitry as a risk phenotype for psychosis. J Am Med Assoc Psychiatry. (2013) 70:1143–51. doi: 10.1001/jamapsychiatry.2013.1976

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Jones JG, Cen SY, Lebel RM, Hsieh PC, Law M. Diffusion tensor imaging correlates with the clinical assessment of disease severity in cervical spondylotic myelopathy and predicts outcome following surgery. Am J Neuroradiol. (2013) 34:471–8. doi: 10.3174/ajnr.A3199

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Kuang C, Zha Y, Liu C, Chen J. Altered topological properties of brain structural covariance networks in patients with cervical spondylotic myelopathy. Front Hum Neurosci. (2020) 14:364. doi: 10.3389/fnhum.2020.00364

PubMed Abstract | CrossRef Full Text | Google Scholar

73. Liu M, Tan Y, Zhang C, He L. Cortical anatomy plasticity in cases of cervical spondylotic myelopathy associated with decompression surgery: a strobe-compliant study of structural magnetic resonance imaging. Medicine. (2021) 100:e24190. doi: 10.1097/MD.0000000000024190

PubMed Abstract | CrossRef Full Text | Google Scholar

74. Aleksanderek I, Stevens TK, Goncalves S, Bartha R, Duggal N. Metabolite and functional profile of patients with cervical spondylotic myelopathy. J Neurosurg Spine. (2017) 26:547–53. doi: 10.3171/2016.9.SPINE151507

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: arterial spin labeling, cerebral blood flow, functional magnetic resonance imaging, functional connectivity, cervical spondylotic myelopathy, resting-state fMRI

Citation: Wei W, Wang T, Abulizi T, Li B and Liu J (2021) Altered Coupling Between Resting-State Cerebral Blood Flow and Functional Connectivity Strength in Cervical Spondylotic Myelopathy Patients. Front. Neurol. 12:713520. doi: 10.3389/fneur.2021.713520

Received: 23 May 2021; Accepted: 20 July 2021;
Published: 08 September 2021.

Edited by:

Wei Zhang, Peking University First Hospital, China

Reviewed by:

Jiajia Zhu, First Affiliated Hospital of Anhui Medical University, China
Jie Yang, Central South University, China

Copyright © 2021 Wei, Wang, Abulizi, Li and Liu. 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: Jun Liu, bGl1anVuNTEwOTk5ODEyQDE2My5jb20=

These authors have contributed equally to this work

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.