Yi Zhou
Jiaqi Shi*- Department of Orthopedics, Xiangyang Central Hospital, Affiliated Hospital of Hubei University of Arts and Science, Xiangyang, China
Background: Previous studies have shown the whole-brain global functional connectivity density (gFCD) and gray matter volume (GMV) alterations in patients with degenerative cervical myelopathy (DCM). However, no study aimed to investigate the associations between the spatial patterns of GMV and gFCD alterations in patients with DCM.
Methods: Structural data and resting-state functional MRI data of 35 DCM patients and 35 matched healthy controls were collected to assess their gFCD and GMV and investigate gFCD and GMV alterations in patients with DCM and their spatial pattern associations.
Results: In our current study, significant gFCD and GMV differences were observed in some regions of the visual system, sensorimotor cortices, and cerebellum between patients with DCM and healthy controls. In our findings, decreased gFCD was found in areas primarily located at the sensorimotor cortices, while increased gFCD was observed primarily within areas located at the visual system and cerebellum. Decreased GMV was seen in the left thalamus, bilateral supplementary motor area (SMA), and left inferior occipital cortices in patients with DCM, while increased GMV was observed in the cerebellum.
Conclusion: Our findings suggest that structural and functional alterations independently contributed to the neuropathology of DCM. However, longitudinal studies are still needed to further illustrate the associations between structural deficits and functional alterations underlying the onset of brain abnormalities as DCM develops.
Introduction
Degenerative cervical myelopathy (DCM), characterized by acquired cervical spine stenosis during the aging process (1), is becoming increasingly prevalent due to global aging and lifestyle changes. DCM has also been considered the most frequently reported non-traumatic spinal cord dysfunction in clinical practice (2), which requires timely diagnosis and spinal canal decompression to relieve myelopathy before irreparable damage to the spinal cord (3, 4). Although the surgical indications for DCM are quite clear (5), it is still difficult for spine surgeons to decide whether patients with long-term myelopathy (i.e., lasting for years or decades) or mild myelopathy would benefit from decompression surgery (6). Subsequently, simple and accurate markers are needed to monitor the myelopathy progression and inform surgical decisions (7). In the past decade, researchers have relied on conventional cervical spine MRI to monitor and determine spinal cord structural changes in DCM patients, e.g., high signal intensity on T2-weighted MR (8, 9). Its utility remains controversial as the information obtained from the spinal cord area is limited, i.e., the small cross-sectional area (10–12). Therefore, researchers set their sights on resting-state functional MRI for exploring brain alterations in patients with DCM (13–15). Several studies reported that DCM is associated with significant brain structural and functional changes (13–24). It has been shown that patients with DCM exhibit a significantly decreased gray matter volume (GMV) in various cortical and subcortical brain regions. Notably, the frontal, sensorimotor, and occipital cortices, thalami, and cerebellum are key regions with GMV decrease in DCM patients (13, 17, 18, 20–23). Further, results from fMRI studies indicated that functional alterations, which are characterized by reduced functional activity in subcortical regions such as the basal ganglia and some regions of the limbic system, and concurrent hyper-functional activity in cortical regions were observed in patients with DCM (15, 17–19, 25, 26). They have also been shown to exhibit functional alterations within frontal, temporal, and parietal lobes, sensorimotor cortex, cingulate gyrus, occipital gyrus, hippocampus, and thalami regions (26–32).
Whole-brain measurement of voxel-based morphometry (VBM) is an unbiased and fully automated method for evaluating GMV alterations (33, 34). These GMV changes can serve as structural biomarkers for clinical and research applications (33, 35). Functional connectivity density mapping (FCDM) is a voxel-wise data-driven method used to determine the number of functional connections between a given voxel and other voxels in the whole brain (36, 37). FCDM measures global functional connectivity density (gFCD) (38–41), an indicator of the number of resting-state functional connections of a given voxel with all other voxels in the entire brain, reflecting the one-to-many relationship (42). Traditional functional connectivity methods mainly involve measuring the connectivity strength between two voxels, regions, or networks, thus only reflecting a one-to-one relationship (42–44). In contrast, gFCD reflects the aberrant functional activity from several perspectives compared with traditional functional connectivity methods (42). A recent study reported that gFCD is a key biomarker for various diseases (45). Although previous studies reported the aberrations in GMV and gFCD as important characteristics of DCM (18, 46, 47), the associations between the spatial distribution patterns of GMV and gFCD alterations remain unexplored.
Therefore, we adopted VBM and FCDM approaches in the present study to explore the alterations in GMV and gFCD and their spatial distribution patterns in patients with DCM. Several previous studies reported that DCM is associated with significant anatomical and functional brain aberrations (47, 48). Anatomical aberrations are mainly located within sensorimotor regions, whereas functional aberrations are mainly observed in frontal, sensorimotor, and default-mode network regions (16–19). Notably, spatial and temporal distributions of anatomical and functional alterations are not correlated, and the aberrant patterns are inconsistent, indicating relatively “distinct patterns” (48). This study explores whether gFCD and GMV are altered in patients with DCM compared with matched healthy controls. An analysis was conducted to explore whether gFCD and GMV were simultaneously affected in particular brain regions and the relationships between the spatial distributions and aberrant patterns (increased or decreased).
Materials and Methods
Subjects
This retrospective study adheres to the STROBE guidelines (49). Ethical approval for this study was obtained from the local institutional review board. All participants provided written informed consent before participating in the current study. A total of 35 right-handed DCM patients were recruited consecutively from 2014 to 2020 according to the following inclusion criteria: (1) Patients with explicit cervical MR evidence of spinal cord compression (from C1 to C7); (2) patients clinically diagnosed with cervical myelopathy, with MRI findings consistent with the level of clinical signs. The common presenting symptoms include neck pain or stiffness, pain, weakness, or numbness (paraesthesias) in the upper limbs, loss of manual dexterity (clumsiness), stiffness, weakness, or numbness of the lower limbs, gait imbalance or unsteadiness, and falls. Objective physical signs of myelopathy include upper motor neuron signs in the upper and/or lower limbs (1); (3) patients with no history of surgery at the spinal canal; (4) patients capable of finishing the fMRI scan; (5) patients with no history of other diseases; (6) patients with no alcohol and substance abuse history. A total of 35 healthy subjects with matching age, gender, and education were recruited through advertisements according to the following inclusion criteria: (1) Individuals with no MR or clinical evidence of myelopathy; (2) individuals with no other spinal or brain neurological disorders or systemic diseases; (3) individuals capable of finishing the fMRI scan; (4) individuals with no history of other diseases; (5) Individuals with no alcohol and substance abuse history.
MRI Data Acquisition
MRI was performed using a 3.0-Tesla MR system (Discovery MR750, General Electric, Milwaukee, WI, USA). Tight but comfortable foam padding was used to minimize head motion, and earplugs were used to reduce scanner noise during MRI examination. Sagittal 3D T1-weighted images were acquired using a brain volume sequence with the following parameters: repetition time (TR) = 7.8 ms; echo time (TE) = 2.8 ms; inversion time (TI) = 440 ms; flip angle (FA) = 12°; field of view (FOV) = 256 × 256 mm; matrix = 256 × 256; slice thickness = 1 mm, no gap; and 188 sagittal slices. Resting-state fMRI data were acquired using a gradient-echo single-short Echo planar imaging sequence with the following parameters: TR/TE = 2,000/30 ms; FOV = 220 × 220 mm; matrix = 64 × 64; FA = 90°; slice thickness = 3 mm; gap = 0 mm; 48 interleaved transverse slices; and 180 volumes. All subjects were requested to keep their eyes closed, relax, minimize movement, think of nothing in particular, and not fall asleep during the scanning period. The images were evaluated right after the data acquisition by a senior radiologist, and subjects with low-quality structural or functional images were asked to retake the scan. Each patient was assessed by a senior orthopedic surgeon using the Japanese Orthopedic Association (JOA) score, which measures the sensorimotor function in patients with DCM.
fMRI Data Preprocessing
Resting-state fMRI data were preprocessed using the Statistical Parametric Mapping software (SPM8; http://www.fil.ion.ucl.ac.uk/spm). The first 10 volumes of each subject were discarded to obtain the equilibrated signal and allow the participants to adapt to the scanning noise. The remaining 170 volumes were corrected for the acquisition time delay between slices. Realignment was then performed to correct for motion between time points. The fMRI data of all subjects were within the defined motion thresholds, namely, translational or rotational motion parameters below 2 mm or <2°. The frame-wise displacement, which indexes volume-to-volume changes in head position, was also determined. Several nuisance covariates, including the six motion parameters, their first-time derivations, and average BOLD signals of the ventricular and white matter, were regressed from the data. A recent study reported that the signal spike caused by head motion significantly contaminates the final resting-state fMRI results even after adjusting for the six motion parameters (50). Therefore, spike volumes were adjusted when the frame-wise displacement (determined using the FD Jenkinson method) of the specific volume exceeded 0.5. The datasets were then band-pass filtered in a frequency range of 0.01 to 0.08 Hz (51). Individual structural images were linearly co-registered with the mean functional image in the normalization step and then linearly co-registered to the Montreal Neurological Institute (MNI) space. Each filtered functional volume was then spatially normalized to the MNI space using the co-registration parameters and resampled into a 3-mm cubic voxel.
GMV Calculation
The GMV of each voxel was determined using SPM12. Structural MR images were classified into gray matter (GM), white matter, and cerebrospinal fluid using the standard unified segmentation model. GM concentration images were non-linearly warped using the diffeomorphic anatomical registration through exponentiated lie algebra technique after the initial affine registration of the GM concentration map into the MNI space. The results were then resampled to a voxel size of 3 × 3 × 3 mm. The GMV of each voxel was obtained by multiplying the GM concentration map with the non-linear determinants derived from the spatial normalization step. The GMV images were further smoothed using a Gaussian kernel with 6 × 6 × 6 mm full-width at half maximum. The smoothed GMV maps were used for statistical analyses after spatial preprocessing.
Determination of FCD
The FCD of each voxel was determined using a Linux script developed in-house according to the method described by Tomasi and Volkow (45). Pearson's linear correlation analysis was conducted to evaluate the functional connectivity strength between voxels (52). Voxel pairs with a correlation coefficient of R > 0.6 were considered significantly connected and used for subsequent analysis (38–40, 42, 43). The FCD calculations were restricted to the cerebral gray matter mask regions. The gFCD at a given voxel, x0, was computed as the global number of functional connections, k(x0), between x0 and all other voxels. This calculation was performed for all x0 voxels in the brain. The FCD maps were spatially smoothed using a 6 × 6 × 6 mm Gaussian kernel to minimize differences in functional brain anatomy across subjects.
Statistical Analysis
Group differences in GMV and gFCD were compared in a voxel-wise manner using a two-sample t-test with age and gender as confounding variables. Multiple comparisons were corrected using the false discovery rate method with a corrected threshold of P < 0.05. Mean gFCD and GMV values of each cluster, with significant group differences in gFCD, were extracted for each subject and used for the region of interest (ROI)-based group comparisons. Group differences in GMV for each ROI were evaluated using two-sample t-tests after controlling for age and gender to explore if the regions with altered gFCD exhibited structural impairment. Pearson correlation analysis was conducted to explore the relationship between gFCD and GMV in each group and evaluate the associations between gFCD and clinical variables, including JOA scores and illness duration.
Results
Subject Demographics and Clinical Characteristics
The demographic and clinical characteristics of the subjects are summarized in Table 1. The results showed no significant gender or age differences between the two groups (t = 0.46, P = 0.53). A total of 35 patients underwent fMRI examinations, and the mean duration of illness was 32.4 ± 15.8 months. The mean preoperative JOA score was 11.2 ± 3.9.
Table 1. Demographic data of the two groups.
Distributed Specificity of gFCD Hubs in Patients and Controls
The results suggested that the strongest gFCD hubs in the two groups were mainly located in the precuneus, inferior parietal lobe, superior temporal gyrus, medial prefrontal cortex, and dorsal lateral prefrontal cortex. Most of the gFCD hubs were located in the default mode network and sensory cortices (Figure 1). The distribution of gFCD hubs was not significantly different between the two groups.
Figure 1. The spatial distributions of the averaged rescaled global functional connectivity density hubs between DCM patients and healthy controls. gFCD, global functional connectivity.
gFCD Changes in Patients With DCM
Voxel-wise analysis showed significantly increased gFCD (i.e., q-value < 0.05, as corrected by the FDR method) in bilateral inferior occipital cortices and cerebellum (Figure 2A, Table 2). Significantly decreased gFCD (i.e., q-value < 0.05, as corrected by the FDR method) was observed in the sensorimotor cortex, including bilateral supplementary motor area (SMA), left thalamus, and left superior frontal gyrus in patients with DCM (Figure 2B, Table 2).
Figure 2. The global functional connectivity density differences between degenerative cervical myelopathy patients and healthy controls. IOC, inferior occipital cortex; CV, cerebellum vermis; SMA, supplementary motor area; Tha, thalamus; SFG, superior frontal gyrus; l, left; r, right. *P < 0.05.
Table 2. Detailed information on intergroup global functional connectivity density (gFCD) differences.
Gray Matter Volume Changes in Patients With DCM
Decreased GMV was observed in the left thalamus, bilateral supplementary motor area (SMA), and cerebellum in patients with DCM (Figure 2, Table 3). In comparison, increased GMV was observed in the left inferior occipital cortices.
Table 3. Correlation coefficients among gFCD, GMV, and clinical measures.
Associations Between Spatial Distributions of gFCD and GMV
As previously mentioned, gFCD increases were primarily located in visual cortices and cerebellum. The gFCD decreases were primarily located in the sensorimotor regions, e.g., bilateral SMA, left thalamus, and left superior frontal gyrus. GMV decreases were observed in the left thalamus, bilateral SMA, and cerebellum in patients with DCM, while GMV increases were observed in the left inferior occipital cortices. Notably, the aberrant patterns of gFCD and GMV differ between regions. In the cerebellum, gFCD increased while GMV decreased, whereas in bilateral SMA and left thalamus, the decreased gFCD was accompanied by decreased GMV. In the left inferior occipital cortices, increased gFCD was accompanied by increased GMV. These findings implied a complex association between gFCD aberrations and GMV changes in DCM. We also performed correlation analysis to further investigate whether gFCD correlates with GMV changes. Multiple regions showed no statistical correlation between gFCD and GMV despite simultaneously being affected, either in patients with DCM or healthy controls.
Association Between gFCD/GMV and Clinical Variables
No statistical correlation was found between gFCD or GMV and illness duration or preoperative JOA scores. The detailed information can be found in Table 3. The reason was that despite several variables being significantly correlated, i.e., P < 0.05, these P-values could not survive the multiple comparison correction.
Discussion
In this study, we first applied the voxel-wise graph theory to investigate the spatial distribution patterns of gFCD and GMV alterations in patients with DCM. According to the results, gFCD and GMV were both affected in some regions of the visual system, sensorimotor cortices, and cerebellum, which are the key regions contributing to the pathology of DCM as indicated by several lines of evidence (15–17, 19, 24). More importantly, our findings showed that decreased gFCD was found in areas primarily located in the sensorimotor cortices, while increased gFCD was observed primarily in areas located in the visual system and cerebellum. Decreased GMV was observed in the left thalamus, bilateral SMA, and cerebellum in patients with DCM, while increased GMV was observed in the left inferior occipital cortices. Furthermore, gFCD and GMV were both affected in most of these regions, but the spatial distribution pattern of gFCD and GMV differed between regions. gFCD and GMV consistently decreased in regions other than the bilateral SMA and left thalamus and consistently increased in regions other than the left inferior occipital cortices. In the sensorimotor cortex, gFCD and GMV alterations were increased in DCM patients, with increased gFCD and increased GMV simultaneously observed in these regions.
Altered Brain Function and Structure Within the Sensorimotor Cortex Indicated Cortical Compensatory Changes Following Deficits of the Corticospinal Tract in DCM Patients
Our current findings are consistent with previous findings that patients with DCM exhibited gray matter atrophy within the sensorimotor cortex (47, 48, 53). Studies have shown that DCM patients exhibited decreased gray matter density within the primary motor cortex (M1), primary sensory cortex (S1), and supplementary motor area (SMA) (54). Moreover, several studies have also reported decreased regional functional activity and metabolic demand within these regions, indicating that the altered brain function in DCM patients could result from structural changes in the sensorimotor cortices (16, 20, 23). These findings showed that the neural deficits following chronic spinal cord injury further disrupt the signal transmission along the corticospinal tract, and the brain gradually adapted to these alterations and formed a series of compensatory changes at the cortical level.
The Cerebral Compensatory Changes Were Complex and Appeared in a Wide Range of Cerebral Regions Not Limited to the Sensorimotor Network
Compared with healthy controls, patients with DCM also exhibited altered regional function and functional connectivity within the visual network. The visual cortices, including the primary visual cortex (calcarine gyrus), and the secondary visual cortex (occipital lobe), exhibited significantly increased functional connectivity in patients with DCM (16, 17). These functional alterations were significantly correlated with visual acuity in patients with DCM (16). A resting-state fMRI study was conducted using the dynamic causal model to explore changes in the primary visual cortex, secondary visual cortex, and cerebellum. Their findings showed that patients with DCM exhibited higher bidirectional effective connectivity between the secondary visual cortex and cerebellum, and the increased effective connectivity was correlated with the prognosis of DCM patients (15). Therefore, the adaptive changes following myelopathy play a key role in the neuropathology of DCM. The gFCD increases within the occipital lobe and cerebellum of DCM patients further support previous findings (15–17, 19). In conclusion, the brain compensatory changes in DCM patients were complex and appeared in a wide range of brain regions of multiple brain networks.
The Brain Structural Alterations Exhibited Indirect Relationships With Functional Abnormalities in DCM Patients
According to the results of this study, the association between anatomical changes and functional activity aberrations is complex and has not been fully elucidated. Previous studies reported that structural deficits do not directly relate to functional abnormalities (55, 56). Longitudinal studies aiming to evaluate the onset of connectivity abnormalities as DCM develops should be conducted to explore the nature of associations between the structural deficits and functional connectivity alterations in DCM (56). Taken together, these results implied that aberrant gFCD and GMV are independent characteristics of DCM. Furthermore, no association between aberrant gFCD/GMV and clinical measures was detected in our current study, which is consistent with findings from previous studies (57). Thus, illness duration and symptom severity are not linearly correlated with brain functional or structural alterations (57, 58). In terms of clinical significance, imaging biomarkers are easily acquired for non-invasive measures in clinical practice, which could prove suitable for monitoring the development of DCM. However, this does not necessitate routine brain rs-fMRI for all DCM patients, especially those with progressive myelopathy, i.e., those needing immediate decompression surgery. However, patients with a long history of myelopathy or mild spinal cord compression but no clinical symptoms make it difficult for surgeons to decide whether surgical treatment is justified. Such brain rs-fMRI analyses may have potential utility for monitoring the progression of myelopathy, thus aiding clinical decision making.
Limitations
Our current study has several limitations. Firstly, no postoperative fMRI data were collected due to the possible artifacts and heating problems caused by surgical implants. We will collect these data when it is safe in the future. Secondly, all the patients included underwent long-term conservative intervention before enrolling in our study. Therefore, studies including drug-naïve patients with DCM are also warranted in the future. Thirdly, our current study only analyzed gFCD, while other resting-state fMRI metrics such as the amplitude of low-frequency fluctuations (ALFF) and regional homogeneity (ReHo) are also worth analyzing. Fourthly, our current study only included 35 patients, a relatively small sample. Future studies enrolling a larger sample are needed. Finally, our univariate approach could detect the differences in the amplitude of brain variables between patients and healthy controls but could not fully describe the spatial pattern differences between patients and controls. Therefore, multi-variate approaches are needed for interpreting these associations.
Conclusion
The present study provided new evidence on aberrant gFCD and decreased GMV in multiple brain regions associated with the DCM pathology. The results indicate that gFCD and GMV are simultaneously affected in multiple regions. However, their spatial distribution patterns in the intergroup difference maps were different. These findings suggest that structural and functional alterations may independently contribute to the neurobiology of DCM. However, the nature of the associations between the structural deficits and functional alterations should be explored through longitudinal studies to determine the onset of connectivity abnormalities as DCM develops.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by Xiangyang Central Hospital. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
YZ collected the data and wrote the manuscript. JS revised the manuscript and analyzed the data. Both authors contributed to the article and approved the submitted version.
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.
References
1. 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
2. Boogaarts HD, Bartels RH. Prevalence of cervical spondylotic myelopathy. Eur Spine J. (2015) 24(Suppl.2):139–41. doi: 10.1007/s00586-013-2781-x
3. 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
4. Karadimas SK, Gatzounis G, Fehlings MG. Pathobiology of cervical spondylotic myelopathy. Eur Spine J. (2015) 24(Suppl.2):132–8. doi: 10.1007/s00586-014-3264-4
5. Iyer A, Azad TD, Tharin S. Cervical spondylotic myelopathy. Clin Spine Surg. (2016) 29:408–14. doi: 10.1097/BSD.0000000000000397
6. Lebl DR, Bono CM. Update on the diagnosis and management of cervical spondylotic myelopathy. J Am Acad Orthop Surg. (2015) 23:648–60. doi: 10.5435/JAAOS-D-14-00250
7. Ganau M, Holly LT, Mizuno J, Fehlings MG. Future directions and new technologies for the management of degenerative cervical myelopathy. Neurosurg Clin N Am. (2018) 29:185–93. doi: 10.1016/j.nec.2017.09.006
8. Merali Z, Wang JZ, Badhiwala JH, Witiw CD, Wilson JR, Fehlings MG, et al. deep learning model for detection of cervical spinal cord compression in MRI scans. Sci Rep. (2021) 11:10473. doi: 10.1038/s41598-021-89848-3
9. Horáková M, Horák T, Valošek J, Rohan T, Koritáková E, Dostál M, et al. Semi-automated detection of cervical spinal cord compression with the Spinal Cord Toolbox. Quant Imaging Med Surg. (2022) 12:2261–79. doi: 10.21037/qims-21-782
10. Li XY, Lu SB, Sun XY, Kong C, Guo MC, Sun SY, et al. Clinical and magnetic resonance imaging predictors of the surgical outcomes of patients with cervical spondylotic myelopathy. Clin Neurol Neurosurg. (2018) 174:137–43. doi: 10.1016/j.clineuro.2018.09.003
11. Liu FJ, Sun YP, Shen Y, Ding WY, Wang LF. Prognostic value of magnetic resonance imaging combined with electromyography in the surgical management of cervical spondylotic myelopathy. Exp Ther Med. (2013) 5:1214–8. doi: 10.3892/etm.2013.934
12. Vallotton K, David G, Hupp M, Pfender N, Cohen-Adad J, Fehlings MG, et al. Tracking white and gray matter degeneration along the spinal cord axis in degenerative cervical myelopathy. J Neurotrauma. (2021) 38:2978–87. doi: 10.1089/neu.2021.0148
13. 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
14. 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
15. Zhao R, Song Y, Guo X, Yang X, Sun H, Chen X, et al. Enhanced information flow from cerebellum to secondary visual cortices leads to better surgery outcome in degenerative cervical myelopathy patients: a stochastic dynamic causal modeling study with functional magnetic resonance imaging. Front Hum Neurosci. (2021) 15:632829. doi: 10.3389/fnhum.2021.632829
16. 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
17. 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. (2020) 45:E272–9. doi: 10.1097/BRS.0000000000003245
18. Su Q, Zhao R, Wang S, Tu H, Guo X, Yang F. Identification and therapeutic outcome prediction of cervical spondylotic myelopathy based on the functional connectivity from resting-state functional MRI data: a preliminary machine learning study. Front Neurol. (2021) 12:711880. doi: 10.3389/fneur.2021.711880
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
20. 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
21. 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
22. Bharat JVK, Biswal B, Hyde S. Simultaneous assessment of flow and BOLD signals in resting-state functional connectivity maps. NMR Biomed. (1997) 10:165–70. doi: 10.1002/(SICI)1099-1492(199706/08)10:4/5<165::AID-NBM454>3.0.CO;2-7
23. 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
24. Zhao R, Guo X, Wang Y, Song Y, Su Q, Sun H, et al. Functional MRI evidence for primary motor cortex plasticity contributes to the disease's severity and prognosis of cervical spondylotic myelopathy patients. Eur Radiol. (2022) 2022:3. doi: 10.1007/s00330-021-08488-3
25. 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
26. 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
27. De Vico Fallani F, Astolfi L, Cincotti F, Mattia D, Marciani MG, Salinari S, et al. Cortical functional connectivity networks in normal and spinal cord injured patients: evaluation by graph analysis. Hum Brain Mapp. (2007) 28:1334–46. doi: 10.1002/hbm.20353
28. 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
29. Felix MS, Popa N, Djelloul M, Boucraut J, Gauthier P, Bauer S, et al. Alteration of forebrain neurogenesis after cervical spinal cord injury in the adult rat. Front Neurosci. (2012) 6:45. doi: 10.3389/fnins.2012.00045
30. Henderson LA, Gustin SM, Macey PM, Wrigley PJ, Siddall PJ. Functional reorganization of the brain in humans following spinal cord injury: evidence for underlying changes in cortical anatomy. J Neurosci. (2011) 31:2630–7. doi: 10.1523/JNEUROSCI.2717-10.2011
31. 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 NeurosurgSpine. (2018) 28:379–88. doi: 10.3171/2017.7.SPINE1768
32. Wrigley PJ, Gustin SM, Macey PM, Nash PG, Gandevia SC, Macefield VG, et al. Anatomical changes in human motor cortex and motor pathways following complete thoracic spinal cord injury. Cerebr Cortex. (2009) 19:224–32. doi: 10.1093/cercor/bhn072
33. Ashburner J, Friston KJ. Voxel-based morphometry–the methods. Neuroimage. (2000) 11:805–21. doi: 10.1006/nimg.2000.0582
34. Braggio D, Barbeito-Andrés B, Gonzalez P, Hallgrímsson B, Larrabide I. VBM sensitivity to localization and extent of mouse brain lesions: a simulation approach. Comput Methods Progr Biomed. (2020) 196:105636. doi: 10.1016/j.cmpb.2020.105636
35. Liu F, Guo W, Yu D, Gao Q, Gao K, Xue Z, et al. Classification of different therapeutic responses of major depressive disorder with multivariate pattern analysis method based on structural MR scans. PLoS ONE. (2012) 7:e40968. doi: 10.1371/journal.pone.0040968
36. Zhang J, Bi W, Zhang Y, Zhu M, Zhang Y, Feng H, et al. Abnormal functional connectivity density in Parkinson's disease. Behav Brain Res. (2015) 280:113–8. doi: 10.1016/j.bbr.2014.12.007
37. Zhuo C, Zhu J, Qin W, Qu H, Ma X, Tian H, et al. Functional connectivity density alterations in schizophrenia. Front Behav Neurosci. (2014) 8:404. doi: 10.3389/fnbeh.2014.00404
38. Guo J, Yang M, Biswal BB, Yang P, Liao W, Chen H. Abnormal functional connectivity density in post-stroke aphasia. Brain Topogr. (2019) 32:271–82. doi: 10.1007/s10548-018-0681-4
39. Chen L, Zhan Y, He F, Zhang S, Wu L, Gong H, et al. Altered functional connectivity density in young survivors of acute lymphoblastic leukemia using resting-state fMRI. Cancer Manag Res. (2020) 12:7033–41. doi: 10.2147/CMAR.S253202
40. Wang Y, Kwapong WR, Tu Y, Xia Y, Tang J, Miao H, et al. Altered resting-state functional connectivity density in patients with neuromyelitis optica-spectrum disorders. Mult Scler Relat Disord. (2020) 43:102187. doi: 10.1016/j.msard.2020.102187
41. Li H, Cao W, Zhang X, Sun B, Jiang S, Li J, et al. BOLD-fMRI reveals the association between renal oxygenation and functional connectivity in the aging brain. Neuroimage. (2019) 186:510–7. doi: 10.1016/j.neuroimage.2018.11.030
42. Cohen AD, Tomasi D, Shokri-Kojori E, Nencka AS, Wang Y. Functional connectivity density mapping: comparing multiband and conventional EPI protocols. Brain Imaging Behav. (2018) 12:848–59. doi: 10.1007/s11682-017-9742-7
43. Li W, Zhang J, Zhou C, Hou W, Hu J, Feng H, et al. Abnormal functional connectivity density in amyotrophic lateral sclerosis. Front Aging Neurosci. (2018) 10:215. doi: 10.3389/fnagi.2018.00215
44. Medvedev MG, Bushmarinov IS, Sun J, Perdew JP, Lyssenko KA. Density functional theory is straying from the path toward the exact functional. Science. (2017) 355:49–52. doi: 10.1126/science.aah5975
45. Tomasi D, Volkow ND. Mapping small-world properties through development in the human brain: disruption in schizophrenia. PLoS ONE. (2014) 9:e96176. doi: 10.1371/journal.pone.0096176
46. Wei W, Wang T, Abulizi T, Li B, Liu J. Altered coupling between resting-state cerebral blood flow and functional connectivity strength in cervical spondylotic myelopathy patients. Front Neurol. (2021) 12:713520. doi: 10.3389/fneur.2021.713520
47. 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
48. 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
49. Cuschieri S. The STROBE guidelines. Saudi J Anaesth. (2019) 13:S31–4. doi: 10.4103/sja.SJA_543_18
50. Power JD, Barnes KA, Snyder AZ, Schlaggar BL, Petersen SE. Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage. (2012) 59:2142–54. doi: 10.1016/j.neuroimage.2011.10.018
51. Liu F, Guo W, Liu L, Long Z, Ma C, Xue Z, et al. Abnormal amplitude low-frequency oscillations in medication-naive, first-episode patients with major depressive disorder: a resting-state fMRI study. J Affect Disord. (2013) 146:401–6. doi: 10.1016/j.jad.2012.10.001
52. Liu F, Guo W, Fouche JP, Wang Y, Wang W, Ding J, et al. Multivariate classification of social anxiety disorder using whole brain functional connectivity. Brain Struct Funct. (2015) 220:101–15. doi: 10.1007/s00429-013-0641-4
53. Grabher P, Mohammadi S, Trachsler A, Friedl S, David G, Sutter R, et al. Voxel-based analysis of grey and white matter degeneration in cervical spondylotic myelopathy. Sci Rep. (2016) 6:24636. doi: 10.1038/srep24636
54. Jütten K, Mainz V, Schubert GA, Fabian Gohmann R, Schmidt T, Ridwan H, et al. Cortical volume reductions as a sign of secondary cerebral and cerebellar impairment in patients with degenerative cervical myelopathy. NeuroImage Clin. (2021) 30:102624. doi: 10.1016/j.nicl.2021.102624
55. Ikuta T, Robinson DG, Gallego JA, Peters BD, Gruner P, Kane J, et al. Subcortical modulation of attentional control by second-generation antipsychotics in first-episode psychosis. Psychiatry Res. (2014) 221:127–34. doi: 10.1016/j.pscychresns.2013.09.010
56. Fornito A, Zalesky A, Pantelis C, Bullmore ET. Schizophrenia neuroimaging and connectomics. Neuroimage. (2012) 62:2296–314. doi: 10.1016/j.neuroimage.2011.12.090
57. 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
Keywords: fMRI, VBM, functional connectivity density, degenerative cervical myelopathy, gray matter atrophies
Citation: Zhou Y and Shi J (2022) Brain Structural and Functional Dissociated Patterns in Degenerative Cervical Myelopathy: A Case-Controlled Retrospective Resting-State fMRI Study. Front. Neurol. 13:895348. doi: 10.3389/fneur.2022.895348
Received: 13 March 2022; Accepted: 03 May 2022;
Published: 15 June 2022.
Edited by:
Patrick Stroman, Queen's University, CanadaReviewed by:
Mario Ganau, Oxford University Hospitals NHS Trust, United KingdomRui Zhao, Tianjin Medical University General Hospital, China
Mohamed Soliman, Cairo University, Egypt
Copyright © 2022 Zhou and Shi. 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: Jiaqi Shi, amlhcWlzaGkyMDIyJiN4MDAwNDA7MTYzLmNvbQ==