Kohei Matsubayashi1†
Munehisa Shinozaki2†Junichi Hata3Yuji Komaki4
Narihito Nagoshi1Osahiko Tsuji1Kanehiro Fujiyoshi5Masaya Nakamura1*
Hideyuki Okano2*- 1Department of Orthopaedic Surgery, School of Medicine, Keio University, Tokyo, Japan
- 2Department of Physiology, School of Medicine, Keio University, Tokyo, Japan
- 3Graduate School of Human Health Sciences, Tokyo Metropolitan University, Tokyo, Japan
- 4Live Animal Imaging Center, Central Institute for Experimental Animals, Kanagawa, Japan
- 5Department of Orthopaedic Surgery, Murayama Medical Center (NHO), Tokyo, Japan
Background: Spinal cord injury (SCI) causes severe sequelae and significant social loss, depending on the extent of the damage. Most previous studies have focused on the pathology of the spinal cord to develop treatments for SCI. However, it is now known that the brain, which is not directly damaged, also undergoes morphological changes after spinal cord injury, which could affect natural recovery and treatment. In recent years, magnetic resonance imaging (MRI) has been developed to analyze functional changes in the brain. Resting-state functional MRI (rsfMRI), which captures brain activity at rest, can calculate functional connections between brain areas and identify central hubs by network analysis.
Purpose: We aim to investigate functional connectivity in the brain using rsfMRI after SCI and to determine how brain-network main hubs change over time.
Methods: We evaluated rsfMRI in 10 mice of the contusional SCI model and calculated connectivity using graph theory. We evaluated “centrality,” a representative parameter of network analysis. The subtype of centrality was degree centrality, which indicates the hub function of a single area. The five times of rsfMRI were performed in each individual mouse: before injury and at 1, 3, 7, and 14 weeks post-injury.
Results: Before the injury, the degree centralities of the primary and secondary motor cortex were high, suggesting that these motor cortices served as main hubs for motor function. After SCI, the hub function of the motor cortices decreased by 14 weeks. In contrast, hub function in the external capsule and the putamen comparatively increased with time after injury, suggesting that the extrapyramidal/subcortical system, which runs the ventral side of the spinal cord and remains after injury in this model, becomes dominant.
Conclusion: We demonstrated the shift of the brain network hub after SCI. The results of this study provide basic information for understanding brain network changes after SCI and would be useful for treatment selection and evaluation of its efficacy in SCI patients.
Introduction
Spinal cord injury (SCI) causes dysfunction of the spinal cord at and below the level of injury, followed by impairments such as motor and sensory paralysis. Depending on the extent of the injury, spontaneous recovery is limited and often results in serious sequelae. Worldwide, one million people are affected by SCI each year, yielding a significant social loss (GBD 2016 Traumatic Brain Injury and Spinal Cord Injury Collaborators, 2019). To date, many researches have been conducted to develop treatments for SCI, and most of them have focused on the pathophysiology of the spinal cord itself. Meanwhile, it is known that there are morphological changes in the brain, which is not directly affected after SCI, and it is suggested that this hinders natural recovery and treatment (Jurkiewicz et al., 2006; Wrigley et al., 2009).
Neuromodulation and neurofeedback have been reported to treat brain or spinal cord disease by utilizing the plasticity of the brain. While some approaches to brain diseases are based on the reinforcement of residual brain activity on the affected side (Lorach et al., 2023), there are also treatments based on the theory of interhemispheric inhibition (Ferbert et al., 1992; Hummel and Cohen, 2006), in which interventions are performed on the unaffected side of the brain. Hence, it is necessary to correctly understand the changes in the brain in order to develop a brain-mediated treatment strategy.
Recent advances in magnetic resonance imaging (MRI) have established the methods to analyze brain activity and function. In particular, resting-state functional MRI (rs-fMRI), which captures brain activity at rest, has been developed to analyze functional relationships between different brain regions (Margulies et al., 2010; Lv et al., 2018). Changes in the relations between brain regions have been also reported in SCI studies (Huynh et al., 2019; Zhang et al., 2022), and several approaches are used in rsfMRI to calculate parameters. Among them, graph theory allows for the assessment of connections between hundreds of “nodes” based on functional brain localization. Graph theory analysis typically includes an overall measure such as small-worldness, or individual connectivity, which indicates the strength of connections between each “node.” However, by setting thresholds for each connectivity, traditional “network analysis” can be adapted to reveal the character of the brain network as an unweighted graph.
In this study, we evaluated rsfMRI using a mouse model of SCI to assess “centrality,” a typical network analysis parameter. We focused on degree centrality as a subtype of centrality, which represents the hub function of each node. Imaging was performed before and at 1, 3, 7, and 14 weeks after injury to observe how the major hubs of the brain changed after injury. Focusing on the primary motor cortex, we observed how motor-related hub function shifts after SCI. This is the first study to assess changes in degree centrality in a contusion injury model using longitudinal rsfMRI in mice. When targeting the brain for SCI treatments, brain activity is supposed to change before the symptoms change. Therefore, it is possible that rsfMRI reveals live brain activity which is more sensitive to treatment effects than physical phenotypes such as body function. The results of this study will provide basic information of brain network changes after SCI.
Materials and methods
Surgical procedures
Fourteen female C57BL/6J mice (8 weeks old, weighing 18–20 g, Clair Japan Inc) were used in this study. Surgical procedures are described previously (Matsubayashi et al., 2018). Mice were deeply anesthetized with a combination of ketamine (100 mg/kg) and xylazine (10 mg/kg) and fixed with a stereotaxic instrument. An incision was made to expose the surface of the skull. The surface was covered with a thin layer of dental cement using a brush-up technique (Superbond C&B, Sun Medical, Shiga, Japan). An acrylic head bar (3 × 3 × 27 mm3) was attached to the sagittal suture of the skull using dental cement. Mice were fed a high-energy diet. After the habituation and the first rs-fMRI imaging described below, mice were subjected to contusional incomplete spinal cord injury. Anesthetized with ketamine (100 mg/kg) and xylazine (10 mg/kg), the dura was exposed after laminectomy at the level of the T9/10 vertebrae. An IH impactor (Precision Systems and Instrumentation, United States, 60 kDyn) was used to induce incomplete contusional SCI at the level of the 10th thoracic vertebra, as described previously (Kamata et al., 2021; Shibata et al., 2021; Kitagawa et al., 2022). Mice for locomotor function underwent only the SCI treatment, and locomotor function was observed for 8 weeks. All animal experimental procedures were performed in accordance with the Laboratory Animal Welfare Act and the Guide for the Care and Use of Laboratory Animals (National Institutes of Health, Bethesda, MD, USA). All experiments were approved by the Animal Study/Committee of the Central Institute for Experimental Animals (approval number: 14078).
rs-fMRI in awake mice
Mice were habituated for 1 week each time before rs-fMRI, rs-fMRI was performed before and at 1, 3, 7, and 14 weeks post injury (wpi) as previously described (Ratering et al., 2008; Baltes et al., 2009; Basbaum et al., 2009; Bosshard et al., 2012; Komaki et al., 2016; Yoshida et al., 2016; Matsubayashi et al., 2018). rs-fMRI was performed on a 7.0 Tesla MRI system equipped with active shield gradient at a maximum intensity of 700 mT/m (Biospec; 70/16 Bruker BioSpin, Ettlingen, Germany) with a cryogenic orthogonal radiofrequency (RF) surface probe (CryoProbe; Bruker BioSpin AG, Fällanden, Switzerland) for increased sensitivity. The scan protocol used for rs-fMRI was as previously reported. High-resolution T2-weighted images (T2WI) of the whole brain were acquired using the Rapid Acquisition with Relaxation Enhancement (RARE) method with the following parameters: effective echo time (TE); 48 ms, repetition time (TR); 6100 ms, a RARE factor; 8, number of averages; 4, spatial resolution; 75 × 75 × 300 μm3, and number of slices; 52. The BOLD fMRI signals were acquired using gradient echo-planar imaging with the following parameters: TE; 20 ms, TR; 1,000 ms, fip angle; 55°, number of averages; 1, spatial resolution; 200 × 200 × 500 μm3, and number of slices; 16. Scans were repeated 600 times in 10 min.
Data analysis
rs-fMRI data analysis was performed as previously reported (Matsubayashi et al., 2018). Briefly, we used SPM12 (Wellcome Trust Centre for Neuroimaging, UCL Institute of Neurology, London, UK), and tailored software for MATLAB was used to adjust slice acquisition timing, motion correction, and registration of different mouse brains to C57BL/6J mouse 3D MRI brain atlas. Normalized functional images were smoothed with a 0.6 mm full width at half maximum filter. Functional connectivity analysis was performed using CONN, with temporal bandpass future applied in the range of 0.009 Hz to 0.1 Hz. Functional coupling between areas was defined as “edges” and analyzed using the regions defined in The Mouse Brain in Stereotaxic Coordinates (Paxinos and Franklin, 2019), where 244 functional areas “nodes” were set. Graph theory-measures were evaluated on raw and threshold networks. All graph theory indices were computed using R software and functions from the igrapgh library.1 “Strength” of the whole brain was obtained by adding all the functional connectivity. “Density,” defined as the percentage of possible connections present in a network, is a measure of the overall level of connectivity. To characterize centrality of network analysis, “degree centrality” of each node was calculated, where it corresponds to the number of incoming connections to a single node (Wheeler et al., 2015; Mareckova et al., 2022). To perform network analysis, a threshold value of 0.10 was used for all period data, and correlations that could withstand the threshold were considered connected. Data from a previously reported complete transection SCI animal model was introduced into the study for comparison with the extent of tissue damage (Matsubayashi et al., 2018). Since the previous data had been adjusted according to Allen’s atlas, the corresponding locations were identified according to previous literature (Chon et al., 2019). The centrality of each node was calculated. The three-dimensional configuration was visualized using Processing, version 4.2
Statistical analysis
All data except a supplement are presented as mean ± SEM. A supplement is presented as mean ± SD. Comparisons within the same model group (pre-injury vs. other time-points) were made by paired t-tests. Repeated measure ANOVA was used for comparisons between incomplete contusional SCI and transection-complete SCI.
Results
The overall network is temporarily increased after SCI
Graph theory analysis calculated 29,646 inter-node connectivities for 244 whole-brain (122 per hemisphere) regions (Supplementary Table 1). Each connectivity ranged from −1 to 1. The matrices of the connectivities in all time points are shown in Figure 1A. Connectivity was dominated by positive values. The sum of all connectivities was calculated as strength and divided by the number of connections between nodes to calculate density. Figure 1B shows their change over time. Strength tended to increase at 1w after SCI, suggesting that the brain network was disrupted and reorganized. It then decreased over time and was lower than before the injury at 14 wpi, suggesting that a limited network was gradually formed in 14 weeks. To perform classical network analysis, we set a threshold of 0.10 for connectivity and converted serial data to binary data (Figure 1C).
Figure 1. The matrix representing edges between nodes of the whole brain. (A) Matrix representing edges (connections) between 244 nodes of the whole brain. Those are of before injury (“pre”), 1, 3, 7, and 14 weeks post-injury (wpi) states. The lower left and upper right are mirror images of each other. Each connection takes a value between −1 (strong negative correlation, blue) and 1 (strong positive correlation, red). (B) Strength is the sum of the connectivity of all edges, and density is the average value of the connectivities. After injury, strength tended to increase from the baseline (“pre”), indicating that the brain network is being reorganized, and it decreases over time, indicating the reorganization is settled down. (C) Serial data from −1 to 1 were converted to binary data (0 or 1) using various threshold values for use in network analysis. We used 0.10 as the threshold value in subsequent network analysis with reference to the density value (Error bars are shown in SEM for all figures. n.s.: no significant differences).
The network change subsides over time
Serial data were used to evaluate the overall change over time in detail. Figure 2A shows a matrix of changes at each time period compared to pre-injury. The change of connectivity was different depending on the edges. Strength, the total sum of connectivity, increased rapidly in 1 wpi and remained higher than pre-injury in 3 wpi, while it decreased in 7 wpi and further decreased in 14 wpi (Figure 2B). Figure 2C shows the matrix of changes from the previous period: from 1 to 3 wpi showed the greatest decrease in strength, but from 3 to 7 wpi and 7 to 14 wpi the changes became smaller, indicating that the pathological condition of the brain network was settling down (Figure 2D).
Figure 2. Change in connectivity over time. (A) To observe the change of connectivities in detail, the difference from pre-injury (before injury) was calculated for the serial data and represented in matrices. Red indicates increased edges, and blue indicates decreased edges. Increased edges predominate at 1 wpi and decreasing edges predominate at 14 wpi. (B) A graph showing changes in strength compared to the pre-injury state. The strength tends to be higher than the pre-injury state at 1 wpi, and it gradually decreased to the lowest value at 14 wpi. (C) The difference from the previous time-point was calculated for the serial data and represented in matrices. (D) A graph showing the change in strength from the previous time point; it is negative (decreasing trend) at 3–7 wpi and 7–14 wpi, but not larger than the difference of 1–3 wpi, suggesting that the overall change is getting smaller over time (Error bars are shown in SEM for all figures. n.s.: no significant differences).
Hub functionality is reduced overall after SCI
Using the binary data, we calculated degree centrality, a representative network analysis. The degree centrality of each node represents the number of edges connecting the node and indicates hub function. Hence the potential value ranges from 0 to 243 in this study. Figure 3A shows oblique views of three-dimensional reconstruction based on the coordinates in the atlas, with nodes colored with an HSV lookup table from 0 to 243. See Supplementary Figure 1 for view directions and Supplementary Movie 1 for all view angles of reconstructed images. The actual mean value for each node ranged between 30 and 110. There was an overall decrease in centrality through the 14th week post-injury. Figure 3B shows the change in mean centrality for all nodes. In particular, there was a significant decrease in centrality at 14 weeks post-injury compared to each time point, including pre-injury. Figure 3C shows colored nodes with a ratio of 0.6 to 1.5, with the rate of change calculated with the pre-injury data as 1, for comparison with the pre-injury data. Figure 3D shows the statistical results compared to the pre-injury data demonstrating the number of nodes with significant differences. At 1w after injury, some nodes show a significant increase in centrality compared to pre-injury, while the number of nodes showing a significant decrease in centrality increases with time. Overall, it was found that the hub function of each node tended to decrease after SCI.
Figure 3. Degree centrality of each node. (A) Three-dimensionally reconstructed images of each node colored according to the degree centrality and located by the xyz coordinates. The value 0.10 was used as the threshold for calculating degree centrality from binary data. Degree centrality shows the number of edges of each node and represents the hub function of the network. The highest and lowest values were extracted from all time points and were fitted to a HSV lookup table from 0 to 243 and displayed in a continuous color from blue to red (red indicates nodes with higher degree centrality). (B) The change in mean centrality for all nodes. There was a significant decrease in centrality at 7 and 14 wpi compared to the pre-injury. (C) Three-dimensionally reconstructed images of each node colored according to the change from the pre-injury. (D) The statistical results compared to the pre-injury data. At 1w after injury, some nodes show a significant increase in centrality compared to pre-injury, while the number of nodes showing a significant decline in centrality increases with time (3D overlay image is from atlas.brain-map.org. Error bars are shown in SEM for all figures. ***p < 0.001).
Hubs of motor function shift from cortical to subcortical area
Next, we performed a network analysis focusing on motor function. We extracted nodes that had connectivities to either the left or right primary motor cortex before the injury as motor-related nodes. Only nodes whose connectivity with the primary motor cortex exceeded the threshold (0.1) in each individual were used to calculate centrality. Since no node was below the threshold for connectivity with the primary motor cortex in all individuals (nodes with subthreshold binding to the primary motor cortex were present in some individuals), 244 nodes were included in the analysis. The percentage change from before the injury was calculated. Figure 4A shows the nodes colored with ratios ranging from 0.6 to 1.5. Figure 4B shows the statistical results compared to pre-injury. In contrast to the overall analysis, many nodes show a significant decrease in centrality from 1w after injury, with few nodes showing an increase. Figure 4C shows the change in nodes that had high centrality before injury (approximately 2 SD above average). The highest centrality was seen in the primary and secondary motor cortices, which showed a predominant decrease at 14 weeks post-injury, suggesting that those hub functions were impaired. Figure 4D shows the changes in the nodes that showed high centrality at 14 weeks post-injury. The centrality of these nodes did not increase, but did not decrease, compared to pre-injury, suggesting that they comparatively took on a new main hub function in the overall decline. In particular, the centrality of the putamen and corpus callosum/external capsule was relatively high both before and after injury. Figure 4E shows the changes in centrality in the completely injured model for comparison with the current incomplete injury. Compared to pre-injury, the centrality of the sensorimotor cortex decreased in the incomplete injure model, whereas it increased predominantly in the complete injury model.
Figure 4. Change in degree centrality of nodes of motor function. (A) The percentage change from before the injury was calculated. The nodes colored with ratios ranging from 0.6 to 1.5. (B) Shows the statistical results compared to pre-injury. In contrast to the overall analysis, many nodes show a significant decrease in centrality from 1w after injury, with few nodes showing an increase. (C) The change in nodes that had high centrality before injury. The highest centrality was seen in the primary and secondary motor cortices, which showed a predominant decrease at 14 weeks post-injury. (D) The changes in the nodes that showed high centrality at 14 weeks post-injury. The centrality of these nodes did not increase, but did not decrease, compared to pre-injury, suggesting that they took on a new main hub function in the overall decline. In particular, the centrality of the putamen and corpus callosum/external capsule was relatively high both before and after injury. Please refer to the Supplementary Figure 2 for all time points. (E) The changes in centrality of sensorimotor cortex in the completely injured model for comparison with the current incomplete injury. The nodes of sensorimotor region were analyzed instead of primary motor cortex because the position was slightly different in the Paxson and Allen atlases. Compared to pre-injury, the centrality of the sensorimotor cortex decreased in the incomplete injure model, whereas it increased in the complete injury model (R, right; L, left; CA1-Or, CA1 oriens; CA1-Py, CA1 pyramidale; M1, primary motor cortex; M2, secondary motor cortex; S2, secondary sensorimotor cortex; cc/ec, corpus callosum/external capsule; Tu, olfactory tubercle; Cpu, putamen; A24b, cingulate cortex area A24b; DLEnt, Dorsolateral entorhinal cortex. Error bars are shown in SEM for all figures. *p < 0.05, **p < 0.01).
Discussion
In this study, we investigated longitudinal changes in brain network structure after SCI using rs-fMRI analysis in mice. Compared to intact mice (pre- injury), spinal cord-injured mice showed significant changes in the brain network. Graph theory and network analysis revealed that degree centrality, which represents hub function, was altered in motor-related brain regions. These results suggest that network analysis of rs-fMRI can reflect the pathophysiology of neural systems in SCI mice.
Recently, studies using rs-fMRI in SCI have increased. There are many reports in humans (Choe et al., 2013; Hou et al., 2014, 2016; Chisholm et al., 2015; Min et al., 2015a,b; Zhu et al., 2015; Oni-Orisan et al., 2016; Kaushal et al., 2017; Pan et al., 2017; Chen et al., 2018; Hawasli et al., 2018; Ge et al., 2019; Karunakaran et al., 2020; Li et al., 2020a,b, 2023; Park et al., 2020; Black et al., 2021; Huynh et al., 2021; Nakanishi et al., 2021; Zheng et al., 2021, 2022; Guo et al., 2022; Kim et al., 2022; Vallesi et al., 2022; Wang et al., 2022; Mandloi et al., 2023) and several reports in nonhuman primates (Rao et al., 2014, 2016, 2022; Wu et al., 2017) and rodents (Seminowicz et al., 2012; Matsubayashi et al., 2018; Wu et al., 2020; Sanganahalli et al., 2021, 2022). The analysis of rs-fMRI in SCI seems to be in its infancy due to differences in the various rs-fMRI analysis methods (seed-based, ICA, graph theory, ReHo, ALFF, etc.) (Margulies et al., 2010; Lv et al., 2018), the timing of imaging, and degree of injury, and the lack of sufficient evidence based on animal experiments. For example, while reports of intra-insular network involvement in pain after SCI show similar tendencies (Li et al., 2020a,b; Huynh et al., 2021; Mandloi et al., 2023), reports of changes in motor cortex connectivity in relation to motor function are inconsistent (Hou et al., 2014; Min et al., 2015a,b; Zhu et al., 2015; Hou et al., 2016; Oni-Orisan et al., 2016; Rao et al., 2016; Kaushal et al., 2017; Pan et al., 2017; Hawasli et al., 2018; Matsubayashi et al., 2018; Li et al., 2020a,b; Park et al., 2020; Nakanishi et al., 2021; Kim et al., 2022; Wang et al., 2022). In the present study, network analysis revealed a shift in main hub function from the pyramidal/cortical to the extrapyramidal/subcortical system. Graph theory can withstand network analysis to calculate centrality, etc. because it sets up a large number of nodes to which functional meaning is assigned (Min et al., 2015a,b; Kaushal et al., 2017; Matsubayashi et al., 2018). Recently, network analysis has become possible without using thresholded binary data (Gopalsamy et al., 2023), and there is a possibility that raw serial rs-fMRI data can be adapted. More research using rs-fMRI in SCI should be conducted, especially in animal models.
In this study, we used a thoracic SCI model in which the lower limbs are mainly affected and confirmed that dynamic network changes occur in the motor cortex. In particular, the shift of the motor hub from the cortical to the subcortical system is consistent with a previous human clinical study (Kim et al., 2022). Since cortical hub function was not reduced in the transection-complete SCI model, this may reflect a successful transition to compensatory neural circuits. However, unlike humans, the rodent motor cortex has limited function in voluntary movements and is often lumped together with the “sensorimotor” system. Particularly in the lower limb, the role of the motor cortex is thought to be limited to movements in coordination with the upper limb since there is sufficient movement of the lower limb even after the corticospinal tract is severed, and the axonal terminals are located in the dorsal horn of the gray matter. Why, then, have changes occurred in such a large network? First, the cells of the rodent motor cortex have fibers that extend to the lumbar spinal cord, but there are also branches in the brainstem and cervical spinal cord (Lu et al., 2022), suggesting that the axonal transection affected the larger network. Second, although we could not subdivide the motor cortex because there were only 244 nodes in the entire brain, analysis by the motor cortex subtype might have revealed that the upper limb region is relatively preserved. Third, besides the classical anatomical structures that have been revealed by tracers and electrical organization, new aspects of the brain map have been reported (Graziano, 2016; Gordon et al., 2023). Namely, brain regions could be divided according to their functions, in addition to corresponding body parts. As mentioned above, because rodents have upper and lower limb coordination during walking, the limitation of lower limb capacity may affect the overall function of controlling gait and cause significant changes in the motor cortex. Further studies are needed to elucidate the cause of the large changes in the motor cortex network in rodents after lower limb impairment.
rs-fMRI can calculate many parameters and may be useful in understanding the pathophysiology of SCI. In this study, significant differences were observed at 7 and 14 weeks after injury. In rodent SCI models, spontaneous recovery of motor function reaches a plateau at 6 weeks (Supplementary Figure 3), which is considered the “chronic phase” (Shinozaki et al., 2011, 2014). In the brain, however, further changes have been observed from 7 to 14 weeks, with continuous changes from 1 to 14 weeks in macroscopic measures such as mean order centrality, and more fine changes in individual degree centrality. There are reports of changes in rs-fMRI parameters with rehabilitation after SCI (Chisholm et al., 2015; Rao et al., 2022). Even when if changes in motor function are not detected, pathophysiological indices such as rs-fMRI may be useful in selecting treatment modalities and determining treatment efficacy.
There are several limitations with this paper. First, it mainly uses the same incomplete injured thoracic spinal cord model and does not represent the causal relationship between pathophysiology and brain networks. Second, the data were converted into binary data for network analysis, but there is no clear standard for setting the threshold value, unlike other studies (Gerhard et al., 2011). We chose 0.1 as the threshold because when we increased the threshold above 0.10, information centrality could not be calculated (data not shown). Third, because this study was performed only on the pathophysiological untreated model, it did not consider the choice of treatment or the evaluation of treatment effects, which would be revealed by the protocols with treatment interventions at different time points after injury. Fourth, the intervals of the time series are different, so it is not possible to examine the amount of change and the time series in detail. Fifth, although nodes with edges to primary motor cortices were extracted for the study of locomotion, there may be nodes more universally related to locomotion in rodents. Finally, due to the facility where the MRI was performed, behavioral analysis was not performed in the same individuals, and tissue evaluation was not performed. In particular, the rs-fMRI connectivity data would have been strongly supported if tissue assessment had been performed using tracers.
Our results suggest that network analysis of rs-fMRI data may reflect the pathophysiology of SCI. In the future, network analysis is expected to be used to select treatment modalities for SCI and to evaluate the efficacy of treatment modalities through further evidence of rs-fMRI in SCI.
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 authors.
Ethics statement
The animal study was approved by the Animal Study/Committee of the Central Institute for Experimental Animals (approval number: 14078). The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
OT, KF, MN, and HO contributed to conception and design of the study. KM, MS, JH, YK, and NN performed surgeries and experiments, collected the data, and contributed to figure generation. MS performed the statistical analysis and wrote the sections of the manuscript. All authors read and approved the submitted version.
Funding
This work was supported by the Japan Agency for Medical Research and Development (AMED: grant nos. JP20bk0104114 and JP15bm0204001 to HO and MN; AMED: grant no. 23bm1123037h0001 to HO and MS).
Acknowledgments
We appreciate K. Yasutake, M. Akizawa, T. Kobayashi, T. Harada, and all of whom are members of the Spinal Cord Research Team for their assistance with the experiments and animal care.
Conflict of interest
HO is a compensated scientific consultant at SanBio Co., Ltd. and K Pharma Inc. MN and MS are consultancy scientists at K Pharma Inc.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnmol.2023.1245902/full#supplementary-material
Supplementary Figure 1 | 3D schema of all nodes. Please refer supplementary Movie 1 for all directions.
Supplementary Figure 2 | Spontaneous recovery and stability of the contused incomplete spinal cord injury model used in this study. An open field score (Basso Mouse Scale) was observed. Scores reached a plateau at 3 weeks; N = 4. Error bars are shown as SD.
Supplementary Figure 3 | Changes over time of nodes with high centrality degree before or after injury. Error bars are shown as SEM.
Footnotes
References
Baltes, C., Radzwill, N., Bosshard, S., Marek, D., and Rudin, M. (2009). Micro MRI of the mouse brain using a novel 400 MHz cryogenic quadrature RF probe. NMR Biomed. 22, 834–842. doi: 10.1002/nbm.1396
Basbaum, A. I., Bautista, D. M., Scherrer, G., and Julius, D. (2009). Cellular and molecular mechanisms of pain. Cells 139, 267–284. doi: 10.1016/j.cell.2009.09.028
Black, S. R., King, J. B., Mahan, M. A., Anderson, J., and Butson, C. R. (2021). Functional hyperconnectivity and task-based activity changes associated with neuropathic pain after spinal cord injury: a pilot study. Front. Neurol. 12:613630. doi: 10.3389/fneur.2021.613630
Bosshard, S. C., Grandjean, J., Schroeter, A., Baltes, C., Zeilhofer, H. U., and Rudin, M. (2012). Hyperalgesia by low doses of the local anesthetic lidocaine involves cannabinoid signaling: an fMRI study in mice. Pain 153, 1450–1458. doi: 10.1016/j.pain.2012.04.001
Chen, Q., Zheng, W., Chen, X., Li, X., Wang, L., Qin, W., et al. (2018). Whether visual-related structural and functional changes occur in brain of patients with acute incomplete cervical cord injury: a multimodal based MRI study. Neuroscience 393, 284–294. doi: 10.1016/j.neuroscience.2018.10.014
Chisholm, A. E., Peters, S., Borich, M. R., Boyd, L. A., and Lam, T. (2015). Short-term cortical plasticity associated with feedback-error learning after locomotor training in a patient with incomplete spinal cord injury. Phys. Ther. 95, 257–266. doi: 10.2522/ptj.20130522
Choe, A. S., Belegu, V., Yoshida, S., Joel, S., Sadowsky, C. L., Smith, S. A., et al. (2013). Extensive neurological recovery from a complete spinal cord injury: a case report and hypothesis on the role of cortical plasticity. Front. Hum. Neurosci. 7:290. doi: 10.3389/fnhum.2013.00290
Chon, U., Vanselow, D. J., Cheng, K. C., and Kim, Y. (2019). Enhanced and unified anatomical labeling for a common mouse brain atlas. Nat. Commun. 10:5067. doi: 10.1038/s41467-019-13057-w
Ferbert, A., Priori, A., Rothwell, J. C., Day, B. L., Colebatch, J. G., and Marsden, C. D. (1992). Interhemispheric inhibition of the human motor cortex. J. Physiol. 453, 525–546. doi: 10.1113/jphysiol.1992.sp019243
GBD 2016 Traumatic Brain Injury and Spinal Cord Injury Collaborators (2019). Global, regional, and national burden of traumatic brain injury and spinal cord injury, 1990-2016: a systematic analysis for the global burden of disease study 2016. Lancet Neurol. 18, 56–87. doi: 10.1016/S1474-4422(18)30415-0
Ge, Y., Pan, Y., Wu, Q., and Dou, W. (2019). A distance-based neurorehabilitation evaluation method using linear SVM and resting-state fMRI. Front. Neurol. 10:1105. doi: 10.3389/fneur.2019.01105
Gerhard, F., Pipa, G., Lima, B., Neuenschwander, S., and Gerstner, W. (2011). Extraction of network topology from multi-electrode recordings: is there a small-world effect? Front. Comput. Neurosci. 5:4. doi: 10.3389/fncom.2011.00004
Gopalsamy, T., Thankappan, V., and Chandramohan, S. (2023). An efficient supply management in water flow network using graph spectral techniques. Environ. Sci. Pollut. Res. Int. 30, 2530–2543. doi: 10.1007/s11356-022-22335-y
Gordon, E. M., Chauvin, R. J., Van, A. N., Rajesh, A., Nielsen, A., Newbold, D. J., et al. (2023). A somato-cognitive action network alternates with effector regions in motor cortex. Nature 617, 351–359. doi: 10.1038/s41586-023-05964-2
Graziano, M. S. A. (2016). Ethological action maps: a paradigm shift for the motor cortex. Trends Cogn. Sci. 20, 121–132. doi: 10.1016/j.tics.2015.10.008
Guo, Y., Ge, Y., Li, J., Dou, W., and Pan, Y. (2022). Impact of injury duration on a sensorimotor functional network in complete spinal cord injury. J. Neurosci. Res. 100, 1765–1774. doi: 10.1002/jnr.25069
Hawasli, A. H., Rutlin, J., Roland, J. L., Murphy, R. K. J., Song, S. K., Leuthardt, E. C., et al. (2018). Spinal cord injury disrupts resting-state networks in the human brain. J. Neurotrauma 35, 864–873. doi: 10.1089/neu.2017.5212
Hou, J. M., Sun, T. S., Xiang, Z. M., Zhang, J. Z., Zhang, Z. C., Zhao, M., et al. (2014). Alterations of resting-state regional and network-level neural function after acute spinal cord injury. Neuroscience 277, 446–454. doi: 10.1016/j.neuroscience.2014.07.045
Hou, J., Xiang, Z., Yan, R., Zhao, M., Wu, Y., Zhong, J., et al. (2016). Motor recovery at 6 months after admission is related to structural and functional reorganization of the spine and brain in patients with spinal cord injury. Hum. Brain Mapp. 37, 2195–2209. doi: 10.1002/hbm.23163
Hummel, F. C., and Cohen, L. G. (2006). Non-invasive brain stimulation: a new strategy to improve neurorehabilitation after stroke? Lancet Neurol. 5, 708–712. doi: 10.1016/S1474-4422(06)70525-7
Huynh, V., Lutolf, R., Rosner, J., Luechinger, R., Curt, A., Kollias, S., et al. (2021). Supraspinal nociceptive networks in neuropathic pain after spinal cord injury. Hum. Brain Mapp. 42, 3733–3749. doi: 10.1002/hbm.25401
Huynh, V., Rosner, J., Curt, A., Kollias, S., Hubli, M., and Michels, L. (2019). Disentangling the effects of spinal cord injury and related neuropathic pain on supraspinal neuroplasticity: a systematic review on neuroimaging. Front. Neurol. 10:1413. doi: 10.3389/fneur.2019.01413
Jurkiewicz, M. T., Crawley, A. P., Verrier, M. C., Fehlings, M. G., and Mikulis, D. J. (2006). Somatosensory cortical atrophy after spinal cord injury: a voxel-based morphometry study. Neurology 66, 762–764. doi: 10.1212/01.wnl.0000201276.28141.40
Kamata, Y., Isoda, M., Sanosaka, T., Shibata, R., Ito, S., Okubo, T., et al. (2021). A robust culture system to generate neural progenitors with gliogenic competence from clinically relevant induced pluripotent stem cells for treatment of spinal cord injury. Stem Cells Transl. Med. 10, 398–413. doi: 10.1002/sctm.20-0269
Karunakaran, K. D., Yuan, R., He, J., Zhao, J., Cui, J. L., Zang, Y. F., et al. (2020). Resting-state functional connectivity of the thalamus in complete spinal cord injury. Neurorehabil. Neural Repair 34, 122–133. doi: 10.1177/1545968319893299
Kaushal, M., Oni-Orisan, A., Chen, G., Li, W., Leschke, J., Ward, B. D., et al. (2017). Evaluation of whole-brain resting-state functional connectivity in spinal cord injury: a large-scale network analysis using network-based statistic. J. Neurotrauma 34, 1278–1282. doi: 10.1089/neu.2016.4649
Kaushal, M., Oni-Orisan, A., Chen, G., Li, W., Leschke, J., Ward, D., et al. (2017). Large-scale network analysis of whole-brain resting-state functional connectivity in spinal cord injury: a comparative study. Brain Connect. 7, 413–423. doi: 10.1089/brain.2016.0468
Kim, A. R., Cha, H., Kim, E., Kim, S., Lee, H. J., Park, E., et al. (2022). Impact of fractional amplitude of low-frequency fluctuations in motor- and sensory-related brain networks on spinal cord injury severity. NMR Biomed. 35:e4612. doi: 10.1002/nbm.4612
Kitagawa, T., Nagoshi, N., Kamata, Y., Kawai, M., Ago, K., Kajikawa, K., et al. (2022). Modulation by DREADD reveals the therapeutic effect of human iPSC-derived neuronal activity on functional recovery after spinal cord injury. Stem Cell Rep. 17, 127–142. doi: 10.1016/j.stemcr.2021.12.005
Komaki, Y., Hikishima, K., Shibata, S., Konomi, T., Seki, F., Yamada, M., et al. (2016). Functional brain mapping using specific sensory-circuit stimulation and a theoretical graph network analysis in mice with neuropathic allodynia. Sci. Rep. 6:37802. doi: 10.1038/srep37802
Li, X., Chen, Q., Zheng, W., Chen, X., Wang, L., Qin, W., et al. (2020a). Inconsistency between cortical reorganization and functional connectivity alteration in the sensorimotor cortex following incomplete cervical spinal cord injury. Brain Imaging Behav. 14, 2367–2377. doi: 10.1007/s11682-019-00190-9
Li, X., Wang, L., Chen, Q., Hu, Y., Du, J., Chen, X., et al. (2020b). The reorganization of insular subregions in individuals with below-level neuropathic pain following incomplete spinal cord injury. Neural Plast. 2020, 1–9. doi: 10.1155/2020/2796571
Li, Y., Zhang, Y., Zhou, W., Li, R., Yu, J., Gong, L., et al. (2023). Depression mediated the relationships between precentral-subcortical causal links and motor recovery in spinal cord injury patients. Cereb. Cortex 33, 7250–7257. doi: 10.1093/cercor/bhad035
Lorach, H., Galvez, A., Spagnolo, V., Martel, F., Karakas, S., Intering, N., et al. (2023). Walking naturally after spinal cord injury using a brain-spine interface. Nature 618, 126–133. doi: 10.1038/s41586-023-06094-5
Lu, T., Shinozaki, M., Nagoshi, N., Nakamura, M., and Okano, H. (2022). 3D imaging of supraspinal inputs to the thoracic and lumbar spinal cord mapped by retrograde tracing and light-sheet microscopy. J. Neurochem. 162, 352–370. doi: 10.1111/jnc.15653
Lv, H., Wang, Z., Tong, E., Williams, L. M., Zaharchuk, G., Zeineh, M., et al. (2018). Resting-state functional MRI: everything that nonexperts have always wanted to know. AJNR Am. J. Neuroradiol. 39, 1390–1399. doi: 10.3174/ajnr.A5527
Mandloi, S., Syed, M., Shoraka, O., Ailes, I., Kang, K. C., Sathe, A., et al. (2023). The role of the insula in chronic pain following spinal cord injury: a resting-state fMRI study. J. Neuroimaging 33, 781–791. doi: 10.1111/jon.13117
Mareckova, K., Miles, A., Liao, Z., Andryskova, L., Brazdil, M., Paus, T., et al. (2022). Prenatal stress and its association with amygdala-related structural covariance patterns in youth. Neuroimage Clin. 34:102976. doi: 10.1016/j.nicl.2022.102976
Margulies, D. S., Bottger, J., Long, X., Lv, Y., Kelly, C., Schafer, A., et al. (2010). Resting developments: a review of fMRI post-processing methodologies for spontaneous brain activity. Magma 23, 289–307. doi: 10.1007/s10334-010-0228-5
Matsubayashi, K., Nagoshi, N., Komaki, Y., Kojima, K., Shinozaki, M., Tsuji, O., et al. (2018). Assessing cortical plasticity after spinal cord injury by using resting-state functional magnetic resonance imaging in awake adult mice. Sci. Rep. 8:14406. doi: 10.1038/s41598-018-32766-8
Min, Y. S., Chang, Y., Park, J. W., Lee, J. M., Cha, J., Yang, J. J., et al. (2015a). Change of brain functional connectivity in patients with spinal cord injury: graph theory based approach. Ann. Rehabil. Med. 39, 374–383. doi: 10.5535/arm.2015.39.3.374
Min, Y. S., Park, J. W., Jin, S. U., Jang, K. E., Nam, H. U., Lee, Y. S., et al. (2015b). Alteration of resting-state brain sensorimotor connectivity following spinal cord injury: a resting-state functional magnetic resonance imaging study. J. Neurotrauma 32, 1422–1427. doi: 10.1089/neu.2014.3661
Nakanishi, T., Nakagawa, K., Kobayashi, H., Kudo, K., and Nakazawa, K. (2021). Specific brain reorganization underlying superior upper limb motor function after spinal cord injury: a multimodal MRI study. Neurorehabil. Neural Repair 35, 220–232. doi: 10.1177/1545968321989347
Oni-Orisan, A., Kaushal, M., Li, W., Leschke, J., Ward, B. D., Vedantam, A., et al. (2016). Alterations in cortical sensorimotor connectivity following complete cervical spinal cord injury: a prospective resting-state fMRI study. PLoS One 11:e0150351. doi: 10.1371/journal.pone.0150351
Pan, Y., Dou, W. B., Wang, Y. H., Luo, H. W., Ge, Y. X., Yan, S. Y., et al. (2017). Non-concomitant cortical structural and functional alterations in sensorimotor areas following incomplete spinal cord injury. Neural Regen. Res. 12, 2059–2066. doi: 10.4103/1673-5374.221165
Park, E., Cha, H., Kim, E., Min, Y. S., Kim, A. R., Lee, H. J., et al. (2020). Alterations in power spectral density in motor- and pain-related networks on neuropathic pain after spinal cord injury. Neuroimage Clin. 28:102342. doi: 10.1016/j.nicl.2020.102342
Paxinos, G, and Franklin, KBJ. Paxinos and Franklin’s the mouse brain in stereotaxic coordinates. 5th. San Diego, CA: Academic Press/Elsevier. (2019).
Rao, J. S., Liu, Z., Zhao, C., Wei, R. H., Zhao, W., Yang, Z. Y., et al. (2016). Longitudinal evaluation of functional connectivity variation in the monkey sensorimotor network induced by spinal cord injury. Acta Physiol. 217, 164–173. doi: 10.1111/apha.12645
Rao, J. S., Ma, M., Zhao, C., Zhang, A. F., Yang, Z. Y., Liu, Z., et al. (2014). Fractional amplitude of low-frequency fluctuation changes in monkeys with spinal cord injury: a resting-state fMRI study. Magn. Reson. Imaging 32, 482–486. doi: 10.1016/j.mri.2014.02.001
Rao, J. S., Zhao, C., Wei, R. H., Feng, T., Bao, S. S., Zhao, W., et al. (2022). Neural regeneration therapy after spinal cord injury induces unique brain functional reorganizations in rhesus monkeys. Ann. Med. 54, 1867–1883. doi: 10.1080/07853890.2022.2089728
Ratering, D., Baltes, C., Nordmeyer-Massner, J., Marek, D., and Rudin, M. (2008). Performance of a 200-MHz cryogenic RF probe designed for MRI and MRS of the murine brain. Magn. Reson. Med. 59, 1440–1447. doi: 10.1002/mrm.21629
Sanganahalli, B. G., Chitturi, J., Herman, P., Elkabes, S., Heary, R., Hyder, F., et al. (2021). Supraspinal sensorimotor and pain-related reorganization after a hemicontusion rat cervical spinal cord injury. J. Neurotrauma 38, 3393–3405. doi: 10.1089/neu.2021.0190
Sanganahalli, B. G., Pavuluri, S., Chitturi, J., Herman, P., Elkabes, S., Heary, R., et al. (2022). Lateralized supraspinal functional connectivity correlate with pain and motor dysfunction in rat hemicontusion cervical spinal cord injury. Neurotrauma Rep. 3, 421–432. doi: 10.1089/neur.2022.0040
Seminowicz, D. A., Jiang, L., Ji, Y., Xu, S., Gullapalli, R. P., and Masri, R. (2012). Thalamocortical asynchrony in conditions of spinal cord injury pain in rats. J. Neurosci. 32, 15843–15848. doi: 10.1523/JNEUROSCI.2927-12.2012
Shibata, T., Tashiro, S., Shinozaki, M., Hashimoto, S., Matsumoto, M., Nakamura, M., et al. (2021). Treadmill training based on the overload principle promotes locomotor recovery in a mouse model of chronic spinal cord injury. Exp. Neurol. 345:113834. doi: 10.1016/j.expneurol.2021.113834
Shinozaki, M., Nakamura, M., Konomi, T., Kobayashi, Y., Takano, M., Saito, N., et al. (2014). Distinct roles of endogenous vascular endothelial factor receptor 1 and 2 in neural protection after spinal cord injury. Neurosci. Res. 78, 55–64. doi: 10.1016/j.neures.2013.09.011
Shinozaki, M., Takahashi, Y., Mukaino, M., Saito, N., Toyama, Y., Okano, H., et al. (2011). Novel concept of motor functional analysis for spinal cord injury in adult mice. J. Biomed. Biotechnol. 2011:157458. doi: 10.1155/2011/157458
Vallesi, V., Richter, J. K., Hunkeler, N., Abramovic, M., Hashagen, C., Christiaanse, E., et al. (2022). Functional connectivity and amplitude of low-frequency fluctuations changes in people with complete subacute and chronic spinal cord injury. Sci. Rep. 12:20874. doi: 10.1038/s41598-022-25345-5
Wang, L., Zheng, W., Yang, B., Chen, Q., Li, X., Chen, X., et al. (2022). Altered functional connectivity between primary motor cortex subregions and the whole brain in patients with incomplete cervical spinal cord injury. Front. Neurosci. 16:996325. doi: 10.3389/fnins.2022.996325
Wheeler, A. L., Wessa, M., Szeszko, P. R., Foussias, G., Chakravarty, M. M., Lerch, J. P., et al. (2015). Further neuroimaging evidence for the deficit subtype of schizophrenia: a cortical connectomics analysis. JAMA Psychiatry 72, 446–455. doi: 10.1001/jamapsychiatry.2014.3020
Wrigley, P. J., Gustin, S. M., Macey, P. M., Nash, P. G., Gandevia, S. C., Macefield, V. G., et al. (2009). Anatomical changes in human motor cortex and motor pathways following complete thoracic spinal cord injury. Cereb. Cortex 19, 224–232. doi: 10.1093/cercor/bhn072
Wu, T. L., Byun, N. E., Wang, F., Mishra, A., Janve, V. A., Chen, L. M., et al. (2020). Longitudinal assessment of recovery after spinal cord injury with behavioral measures and diffusion, quantitative magnetization transfer and functional magnetic resonance imaging. NMR Biomed. 33:e4216. doi: 10.1002/nbm.4216
Wu, R., Yang, P. F., and Chen, L. M. (2017). Correlated disruption of resting-state fMRI, LFP, and spike connectivity between area 3b and S2 following spinal cord injury in monkeys. J. Neurosci. 37, 11192–11203. doi: 10.1523/JNEUROSCI.2318-17.2017
Yoshida, K., Mimura, Y., Ishihara, R., Nishida, H., Komaki, Y., Minakuchi, T., et al. (2016). Physiological effects of a habituation procedure for functional MRI in awake mice using a cryogenic radiofrequency probe. J. Neurosci. Methods 274, 38–48. doi: 10.1016/j.jneumeth.2016.09.013
Zhang, L., Wang, L., Xia, H., Tan, Y., Li, C., and Fang, C. (2022). Connectomic mapping of brain-spinal cord neural networks: future directions in assessing spinal cord injury at rest. Neurosci. Res. 176, 9–17. doi: 10.1016/j.neures.2021.10.008
Zheng, W., Wang, L., Chen, Q., Li, X., Chen, X., Qin, W., et al. (2021). Functional reorganizations outside the sensorimotor regions following complete thoracolumbar spinal cord injury. J. Magn. Reson. Imaging 54, 1551–1559. doi: 10.1002/jmri.27764
Zheng, W., Wang, L., Yang, B., Chen, Q., Hu, Y., Du, J., et al. (2022). Cerebellum regulating cerebral functional cortex through multiple pathways in complete thoracolumbar spinal cord injury. Front. Neurosci. 16:914549. doi: 10.3389/fnins.2022.914549
Keywords: spinal cord injury, network analysis, corticospinal tract, motor function, brain network hubs, resting-state functional MRI, (rS-fMRI)
Citation: Matsubayashi K, Shinozaki M, Hata J, Komaki Y, Nagoshi N, Tsuji O, Fujiyoshi K, Nakamura M and Okano H (2023) A shift of brain network hub after spinal cord injury. Front. Mol. Neurosci. 16:1245902. doi: 10.3389/fnmol.2023.1245902
Edited by:
Aijun Wang, University of California, Davis, United StatesReviewed by:
Gerardo Rojas-Piloni, Universidad Nacional Autónoma de México, Campus UNAM-Juriquilla, MexicoXiaodong Sheldon Liu, Beijing University of Chinese Medicine, China
Copyright © 2023 Matsubayashi, Shinozaki, Hata, Komaki, Nagoshi, Tsuji, Fujiyoshi, Nakamura and Okano. 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: Masaya Nakamura, bWFzYUBrZWlvLmpw; Hideyuki Okano, aGlkb2thbm9Aa2Vpby5qcA==
†These authors have contributed equally to this work