- 1Department of Rehabilitation, The Affiliated Hospital of Yunnan University, Kunming, China
- 2The Affiliated Hospital of Yunnan University, Kunming Medical University, Kunming, China
Introduction: The occurrence of Spinal cord injury (SCI) brings economic burden and social burden to individuals, families and society, and the complications after SCI greatly affect the rehabilitation and treatment of patients in the later stage.This study focused on the potential biomarkers that co-exist in SCI and sarcopenia, with the expectation to diagnose and prognose patients in the acute phase and rehabilitation phase using comprehensive data analysis.
Methods: The datasets used in this study were downloaded from Gene Expression Omnibus (GEO) database. Firstly, the datasets were analyzed with the “DEseq2” and “Limma” R package to identify differentially expressed genes (DEGs), which were then visualized using volcano plots. The SCI and sarcopenia DEGs that overlapped were used to construct a protein–protein interaction (PPI) network. Three algorithms were used to obtain a list of the top 10 hub genes. Next, validation of the hub genes was performed using three datasets. According to the results, the top hub genes were DCN, FSTL1, and COL12A1, which subsequently underwent were Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses. We also assessed immune cell infiltration with the CIBERSORT algorithm to explore the immune cell landscape. The correlations between the hub genes and age and body mass index were investigated. To illustrate the biological mechanisms of the hub genes more clearly, a single-cell RNA-seq dataset was assessed to determine gene expression when muscle injury occurred. According to our analysis and the role in muscle, we chose the fibro/adipogenic progenitors (FAPs) cluster in the next step of the analysis. In the sub cluster analysis, we use the “Monocle” package to perform the trajectory analysis in different injury time points and different cell states.
Results: A total of 144 overlapped genes were obtained from two datasets. Following PPI network analysis and validation, we finally identified three hub-genes (DCN, FSTL1, and COL12A1), which were significantly altered in sarcopenic SCI patients both before and after rehabilitation training. The three hub genes were also significantly expressed in the FAPs clusters. Furthermore, following injury, the expression of the hub genes changed with the time points, changing in FAPs cluster.
Discussion: Our study provides comprehensive insights into how muscle changes after SCI are associated with sarcopenia by moving from RNA-seq to RNA-SEQ, including Immune infiltration landscape, pesudotime change and so on. The three hub genes identified in this study could be used to distinguish the sarcopenia state at the genomic level. Additionally, they may also play a prognostic role in evaluating the efficiency of rehabilitation training.
1 Introduction
Both traumatic and non-traumatic spinal cord injury (SCI) associated with unfavorable consequences. According to data from American National Spinal Cord Injury Database, the proportion of SCI patients in the United States is 0.0054%. In addition, the age of patients with SCI is gradually getting younger, which deserves our attention (1). A study by a Chinese team showed that the number of SCI patients in China is increasing at a rate of approximately 10% per year, and will gradually comprise a relatively large population in the future (2). Morever, the latest hospital-based retrospective study from China, it shows that the number of people with spinal cord injury between 35 and 45 approximately equal to the number of people over 65 years old, and most people with spinal cord injury are over 45 years old (3). So, there was an aging trend in spinal cord injury patients, and the older population was the most common population for sarcopenia (4). Again, It is vital to focusing on these SCI patients with sarcopenia. SCI can have serious effects on physiology, psychology, and family. There is a strong correlation between the cost of medical treatment after SCI and the level of tetraplegia; the cost of a high level of tetraplegia is the highest regarding both the initial hospitalization and subsequent treatment plan compared to other levels of tetraplegia (1).
When SCI occurs, the motor system, central nervous system, including the sensory system with its proprioceptive component, as well as the sympathetic nervous system, and other related systems, will all exhibit varying degrees of change. Preventing the deterioration of these conditions is one of the core treatments of SCI. The restoration of spinal cord function and remodeling of neural pathways is one of the research hotspots at present; however, there are still many challenges from application and transformation to implementation (5). Therefore, recovery after SCI also remains a concern. This is particularly true in patients with SCI at the neck level and thoracic level, who will experience a longer rehabilitation process after SCI and face a more severe rehabilitation challenge. Among these sequelae, skeletal muscle is a major concern. After SCI, the paralyzed muscles have reduced metabolism and motor capability (6). This feature has been reported in many studies, in addition to skeletal muscle capability changes from different perspectives, including from the aspects of a decline in quality, performance, strength (7–9).
Sarcopenia is characterized by progressive and widespread skeletal muscle dysfunction with adverse consequences (10). Sarcopenia can be diagnosed with several tests, which include the grip strength test, assessment of the appendicular skeletal muscle mass by dual-energy X-ray absorptiometry, the timed-and-go test, and other methods (10). As mentioned above, patients with SCI can present with a decrease in muscle mass, strength, and performance; however, we still lack the awareness to recognize, prevent, and change it. One study on patients with SCI with an average age of 38 years found that sarcopenia was present in many of them (11). Another study also found that sarcopenia and sarcopenic obesity were common in patients with SCI (12). However, at present, unique criteria and biomarkers for the early sensitive identification and specific diagnosis of SCI are lacking, which not only affects the life and quality of life of patients, but also results in a poor prognosis and complications (11, 13).
Research on the molecular mechanisms of the association between sarcopenia and SCI remains sparse. Studying the molecular relationship between the two diseases could result in early diagnosis, later treatment, and rehabilitation of these patients. Furthermore, it also meets the requirements for precision medicine. The aim of this study was to explore the co-pathogenesis between sarcopenia and SCI. We identified hub genes associated with the two diseases and uses various bioinformatics methods to explore the biological mechanisms. The hub genes identified for sarcopenia and SCI in this study are expected to provide new insights and ideas for the diagnosis, treatment, and rehabilitation of both diseases.
2 Materials and methods
The framework of this work shown in Supplementary Figure S1, and all steps will clearly illustrate in subsequent sections.
2.1 Collecting datasets
We obtained our datasets from Gene Expression Omnibus (GEO): (1) GSE21497 (14) contains data on 20 SCI patients; (2) GSE111016 (15) contains data on 20 healthy subjects and 20 sarcopenia patients; (3) GSE111010 (15) contains data on 14 healthy subjects, nine sarcopenia patients, five low muscle mass patients, and 11 low muscle strength or function patients; (4) GSE111006 (15) contains data on 28 healthy subjects, four sarcopenia patients, five low muscle mass patients, and three low muscle strength or function patients; (5) GSE117525 (16) contains data on 53 young subjects, 73 healthy older subjects, and 61 frail older subjects; (6) GSE142426 (17) contains data on 15 healthy subjects and 15 SCI patients; and (7) GSE138826 (18) contains single-cell transcriptomics data on seven SCI patients. We have listed the detailed information of the datasets in Table 1.
2.2 Identification of differentially expressed genes (DEGs)
We processed the GSE21497 dataset using the “Limma” (19) to identify the DEGs and draw the volcano plot. Notably, we merged the datasets (GSE111016, GSE111006, and GSE111010), and used the “sva” (20) package to remove the batch effects and draw the volcano plots. In this merged dataset, we also compared the low mass group and low muscle function group, as its provided phenotype information. In the analysis, we set the threshold at logFC (log fold change) ≥1 and p value <0.05 to distinguish the significant DEGs by “DESeq2” (21) package.
2.3 Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses
To further explore the function behind DEGs, we used GO and KEGG enrichment analyses to investigate their biological significance using the “clusterProfiler” (22) package. In the GO analysis, we assessed the functional enrichment in three aspects, molecular function (MF), biological process (BP), and cellular component (CC). The p value thresholds of the GO and KEGG enrichment analyses were set as 0.01.
2.4 Protein–protein interaction (PPI) network construction
Firstly, we used the “venn” (23) package to obtain the genes that intersected in the GSE21497 and the merged datasets. Using the intersecting genes, we constructed the PPI network based on the STRING database (24). “Cytoscape” (Version 3.6) (25) was used to visualize the network. Finally, we used the “CytoHubba” (26) plugin to calculate the top genes based on the maximum neighborhood component algorithms, maximal correlation coefficient algorithms, and density of maximum neighborhood component algorithms.
2.5 Validation and evaluation of hub genes
We choose GSE117525, GSE21497, and GSE142426 as the validation datasets to verify the expression of hub genes. It is worth mentioning that GSE117525 and GSE14246 were both used to compare changes before and after rehabilitation training.
To provide the clinic value, we calculated the receiver operating characteristic (ROC) curves and area under the curve (AUC) with 95% confidence intervals (CI) using the “pROC” (27) and “plotROC” (28) packages to estimate the predictive accuracy of the hub genes. We the AUC to >0.60 to identify optimal biomarkers.
2.6 Immune cell infiltration landscape
We choose the GSE21497 and GSE111016 datasets to visualize the 28 types of immune cells by using gene expression using the CIBERSORT algorithm (29). A boxplot was used to exhibit the composition of immune cells in the datasets. We also used Spearman correlation analysis to evaluate the correlation between hub gene expression and immune cells between the different groups.
2.7 GO and KEGG enrichment analyses of hub genes
Based on the above analysis, we performed single gene analysis of the top hub genes. We also used the “clusterProfiler” package to conduct GO and KEGG enrichment analyses for the three genes to explore their underlying biological functions.
2.8 Correlation between body mass index (BMI) and age and hub genes
We further investigated the correlations between clinic information and hub genes using the “ggplot2” (30) package for the GSE21497 and GSE111016 datasets.
2.9 Single-cell data download and processing
The single-cell dataset GSE138826 was download from the GEO database. We filtered out the mitochondrial genes (15%) and hemoglobin genes (0.1%), identified cells with a total number of genes >300, and stipulated that genes must be expressed at least in three cells. All samples were eligible for further analysis.
We perform the SCTranform method to correct the batch effect using the “Seurat” (31) package. After finding the highly variable genes and scale data, the “RunPCA” package was used to reduce the dimension of the processed data. We choose dim = 15 to find cell clusters. Then, we used “UMAP” to perform the cell cluster visualization. We have listed our annotation genes for each cluster in Table 2. Furthermore, we used “Dotplot” and “Featureplot” to assess the expression level of each annotated gene to accurately mark these clusters.
2.10 Sub-cluster extraction and analysis
After visualizing the hub genes using “FeaturePlot,” a highly expressed subgroup (fibro/adipogenic progenitors (FAPs) cluster) of hub gene was extracted. After extracting the desired subpopulation of cells, we performed SCTranform and re-scaling of the raw counts, followed by dimensionality reduction and visualization using “UMAP.” In the FAPs, we listed our annotation genes for each cluster in Supplementary Table S1. The boxplot shows the expression difference in hub genes between different experimental time sequences.
We use “Monocle” (32) (Version 2) to perform the pseudo time trajectory analysis and identified the hub genes in the trajectory inference. The “Monocle” algorithm constructs the branches, and different branches represent cells in different trajectories.
3 Results
3.1 DEG identification
In the GSE21497, we obtain 639 DEGs using the “Limma” package. Next, we merged three datasets, GSE111016, GSE111006, and GSE111010, and used the “DESeq2” package (Figure 1A1). In this merged data, after comparing healthy subjects and sarcopenia subjects, we obtained 7,032 DEGs (Figure 1B1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1373605/fneur-15-1373605-HTML-r1/image_m/fneur-15-1373605-g001.jpg)
Figure 1. The differentially expressed genes (DEGs), and Kyoto Encyclopedia or Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses in GSE21497 and the merged datasets (including: GSE111016, GSE111010, and GSE111006). (A1,B1) Volcano map of DEGs in the GSE21497 and the merged dataset. (A2,A3) KEGG and GO enrichment analyses in GSE21497. (B2,B3) KEGG and GO enrichment analyses in the merged dataset.
3.2 Functional enrichment analysis
We conducted GO and KEGG analyses on the DEGs obtained from the above datasets, respectively, to discover the biological mechanism more clearly behind their respective DEGs.
Due to the excessive results obtained in the enrichment process, we selected the first few pathways with the smallest p value for display. In the GSE21497 dataset, the first few KEGG enrichment results were: carbohydrate digestion and absorption, hypertrophic cardiomyopathy, and viral protein interaction with cytokine and cytokine receptors. The top GO enrichment results were: apical plasma membrane, response to chemokines, and cellular response to chemokines (Figures 1A2,A3). Next, in the merged dataset (GSE11016, GSE111010, and GSE111006), the top KEGG enrichment results were: TGF-beta signaling pathway, protein digestion and absorption, and glutamatergic synapse. The top GO enrichment results were cytokine receptor binding, glycosaminoglycan binding, and receptor ligand activity (Figures 1B2,B3).
3.3 PPI network and identification of hub genes
Firstly, we used a Venn diagram to visualize the overlapping DEGs between the GSE21497 and merged dataset (Figure 2A); 141 DEGs overlapped. The STRING tool was used to construct the PPI network using the overlapping DEGs (Figure 2B). The PPI network contained 60 nodes and 116 edges; we set the interaction score at medium confidence. The “cytoHubba” plugin was applied in the PPI networks with the MNC algorithm, MCC algorithm, and DMNC algorithm (Figures 2C1–C3). Finally, 10 hub genes were selected for further analyses: DCN, FBN1, COL6A3, FSTL1, COL12A1, ADAMTS5, MFAP5, TMP4, GRIA1, and GRIK2 (Supplementary Table S1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1373605/fneur-15-1373605-HTML-r1/image_m/fneur-15-1373605-g002.jpg)
Figure 2. Protein–protein interaction (PPI) network construction and identification of hub genes. (A) The overlapped genes of the GSE21497 and the merged dataset. (B) The PPI network based on the STRING database. (C1–C3) Three different algorithms were used to obtain the hub genes, which included the MCC algorithm, MNC algorithm, and DMC algorithm.
3.4 Validation of hub genes
In the validation using the GSE117525 dataset, we compared the different training stage and different subjects’ statues to validate the hub genes (Figure 3C). In young vs. frail subjects and older vs. frail subjects, the expression of DCN remained significantly upregulated. The expression of COL12A1 significantly changed in the comparison of frail vs. older subjects and frail vs. young subjects. After training, the expression of DCN, FBN1, FSTL1, and COL6A3 were significant altered in both the frail subjects and older subjects’ groups. In the hub genes related to low muscle mass, the expression of TGFB1 and DCN changed significantly after training in both frail subjects and older subjects.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1373605/fneur-15-1373605-HTML-r1/image_m/fneur-15-1373605-g003.jpg)
Figure 3. Validation of the hub genes and the diagnostic value of the hub genes in the three datasets. (A) The significance of the expression level of the three hub genes in GSE21497 between day 2 and day 5. (B) The significance of the expression level of the three hub genes in GSE142426 between the control group and training group in patients with spinal cord injury (SCI). (C) The significance of the expression level of the three hub genes in sarcopenia in the GSE117525 dataset. *p < 0.01, **p < 0.05, ***p < 0.001.
In the validation using the GSE142426 dataset, we compared the SCI patients with healthy subjects after rehabilitation training in the lower limbs (Figure 3B). The change in expression of DCN and ADAMTS5 were observed. In the hub genes related to low muscle mass, the expression MSTN also changed significantly.
In the validation using the GSE21497 dataset, we compared the SCI patients at different times (Figure 3A). We observed that on day 5 there was changes in the expression of several hub genes, including DCN, FSTL1, COL12A1, MFAP5, KCNA1, GRIK2, FBN1, COL6A3, ADAMTS5.
3.5 Evaluation of the diagnostic value of the hub genes
After validation with the above datasets, we selected DCN, FSTL1, and COL12A1 for subsequent analyses and further evaluation. We evaluated the three hub genes against SCI and sarcopenia, respectively, using ROC curve analysis.
In the GSE21497 and GSE142426 datasets, the AUC values for the three hub genes were > 0.60 (Figures 4A,B). In the GSE117525 dataset, we performed different comparisons to show the diagnosis value. In the group of frail subjects and young subjects, DCN and FSTL1 showed AUC values >0.60; the AUC values of COL12A1 were lower, but it still >0.55 (Figure 4C1). When comparing before and after training exercise, the AUC values of three the hub genes were > 0.60 (Figure 4C2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1373605/fneur-15-1373605-HTML-r1/image_m/fneur-15-1373605-g004.jpg)
Figure 4. (A) Diagnostic value of the three hub genes in GSE21497. (B) Diagnostic values of the three hub genes in GSE142426. (C1,C2) Diagnostic value of the three hub genes in GSE117252.
3.6 The immune infiltration landscape of the hub genes
The results of the CIBERSORT algorithm revealed the immune cell infiltration landscape in both the SCI and sarcopenia after comparing the different groups. In the GSE21497 dataset, we assessed the proportion of immune cells and landscape (Figures 5A1,A2). FSTL1 had a significant correlation with M2 macrophages and neutrophils; DCN had a significant correlation with M1 macrophages, neutrophils, and CD4+ T cells (Figure 5A3). In the merged dataset (Figures 5B1,B2), FSTL1 had a significant correlation with neutrophils and plasma cells; DCN had a significant correlation with neutrophils (Figure 5B3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1373605/fneur-15-1373605-HTML-r1/image_m/fneur-15-1373605-g005.jpg)
Figure 5. Immune infiltration and correlation with immune cells. (A1,B1) The composition and proportion of immune cells in the GSE21497 and the merged dataset. (B2,B3) The landscape of 28 immune cells in the GSE21497 and the merged dataset. (A2,A3) The correlation of the three hub genes with immune cells in the GSE21497 and the merged dataset (HC: healthy control).
3.7 Enrichment analysis of the hub genes
DCN was markedly enriched in the TGF-beta signaling pathway in the KEGG analysis, and showed significant relationships with mitochondrial depolarization, regulation of mitochondrial depolarization, and negative regulation of cellular response to vascular endothelial growth factor stimulus. In the GO enrichment analysis, DCN was related to mitochondrial depolarization, regulation of mitochondrial depolarization, and negative regulation of cellular response to vascular endothelial growth factor stimulus (Figures 6A1,A2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1373605/fneur-15-1373605-HTML-r1/image_m/fneur-15-1373605-g006.jpg)
Figure 6. Kyoto Encyclopedia or Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment of the hub genes. (A1,A2) KEGG and GO enrichment of DCN. (B) GO enrichment of FSTL1. (C1,C2) KEGG and GO enrichment of COL12A1.
FSTL1 failed to show enrichment in the KEGG analysis. However, in the GO analysis, it showed significant relationships with the BMP signaling pathway, endothelium development, and other relevant biological processes (Figure 6B).
In the KEGG enrichment, COL12A1 only showed enrichment in protein digestion and absorption. In the GO enrichment, it showed strong relationships with endoderm development, collagen fibril organization, and endoderm formation (Figures 6C1,C2).
3.8 Correlations of hub genes with clinical features
In the GSE117525 dataset, we analyzed the relationships between age and BMI and hub genes. All three hub genes showed positive correlations with BMI and age (Figure 7).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1373605/fneur-15-1373605-HTML-r1/image_m/fneur-15-1373605-g007.jpg)
Figure 7. The relationships between the hub genes and age and body mass index. (A1,A2) The correlation between DCN and age and BMI. (B1,B2) The correlation between FSTL1 and age and BMI. (C1,C2) The correlation between COL12A1 and age and BMI.
3.9 Analysis of the single-cell RNA-seq dataset
After quality control, we obtained eight cell clusters, which were annotated with cell markers: CD14, PTPRC, CLEC12A, ADGRE1, CSF1R, H2-AB1, H2-EB1, CD34, HIC1, PDGFRA, PDGFRB, THY1, LY6A, ASB5, MYF5, S100A8, S100A9, CXCR4, MYH4, TNNC2, and KDR (Figure 8C). Finally, whole cells were allocated into eight cell clusters, including immune cells, FAP cells, MuSC cells, neutrophils, tenocytes, and endothelial Cells (Figure 8B). We also show the changes in cell clusters at different time points (Figure 8A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1373605/fneur-15-1373605-HTML-r1/image_m/fneur-15-1373605-g008.jpg)
Figure 8. Cell clusters landscape and expression level. (A) The change in cell populations at different time points. (B) The UMAP visualization. (C) The violine plot show the expression level of annotation markers (DPI: days post injury).
The cell expression level of the three hub genes can be visualized in the boxplot, which clearly demonstrates that the FAPs cluster expression is significant (Figures 9A–C). We also show the expression level of the three hub genes using “featureplot” (Figure 9D). Moreover, FAPs play a crucial role in maintaining muscle homeostasis and promoting regeneration (33, 34).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1373605/fneur-15-1373605-HTML-r1/image_m/fneur-15-1373605-g009.jpg)
Figure 9. The expression difference of hub genes in different cell clusters. (A–C) A boxplot was used to visualize the expression difference of hub genes in difference cell clusters. (D) “Featureplot” was used to visualize the expression difference of the hub genes in difference cell clusters. *p < 0.01, **p < 0.05, ***p < 0.001.
3.10 Pseudotime trajectory analysis of the FAPs cluster
After extracting the FAPs cluster, we obtained the 11 clusters and annotated them, including Ors1+ FAPs, Wisp1+ FAPs, Dlk1+ FAPs, Dpp4+ FAPs, Cxcl14+ FAPs, Bgn + FAPs, fibroblasts, Cxcl5+ FAPs, Csfr1+ FAPs, activated FAPs, and tenocytes (Figure 10A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1373605/fneur-15-1373605-HTML-r1/image_m/fneur-15-1373605-g010.jpg)
Figure 10. Identification of the fibro/adipogenic progenitors (FAPs) cluster and the regeneration trajectory-related analysis. (A) “UMAP” was used to visualize in the FAPs cluster. (B) The different trajectories: (1) trajectory inference analysis of FAPs by annotated markers; (2) trajectory inference analysis by different cell states; (3) psuedo time analysis colored by different shades of color; and (4) trajectory inference analysis by different time points. (C) The expression difference of the hub genes in different cells in trajectory interference. (D) The expression difference of the hub genes in the different groups. (E) The expression level of the hub genes in different cell states. *p < 0.01, **p < 0.05, ***p < 0.001.
We used “Monocle” to reduce the dimension with the “UMAP” method and to construct the pseudotime trajectory. After calculating the score of the pseudotime value, we obtained the pseudotime trajectory. Firstly, we performed the trajectory analysis according to the 11 cell types (Figure 10B). Through the evaluation of different states, we obtained seven differentiation states; red is the earliest differentiation state, and to the left is the subsequent differentiation state (Figure 10B). From the trajectory, we can also see that the darker the blue, the earlier the differentiation, and the lighter the blue, the later the differentiation (Figure 10B). Furthermore we show the pseudotime trajectory at different time points (Figure 10B).
According to the expression level of the three hub genes, we assessed the modulation change patterns. The expression of COL12A1 showed a decreasing trend. The expression of DCN and FSTL1 showed the same ascending trend, although the former one had a greater change rate. Furthermore, both peak at the end of trajectory (Figure 10C).
We also compared the expression level of the hub genes by group and differentiation states (Figures 10D,E).
4 Discussion
The correlation between SCI and sarcopenia has gained increasing attention in recent years. Following SCI, skeletal muscle atrophy is approximately 30–60%, and changes, such as decreased muscle endurance and increased fatigability, often occur (35). Patients with SCI are more likely to develop sarcopenia than the general population due to limb dysfunction, especially those with high spinal segmental injury, long-term bed state, or reduced nutritional intake (11, 13, 36, 37). In the 2019 European consensus, the criteria for sarcopenia included low muscle strength, low muscle quantity or quality, and low physical performance (10). In recent years, people have also become increasingly aware of the emergence of sarcopenia after SCI and have implemented many rehabilitation measures to improve muscle changes, such as electrical stimulation, exercise training, and other rehabilitation therapies (7, 38).
After the acute phase of treatment ends for many patients with SCI, most need to enter the rehabilitation phase. This link is extremely important for patients, not only to improve the remaining function, but also to replace the lost function learning. When a SCI occurs, muscle mass decreases significantly after 6 weeks; however, muscle mass can be increased by body weight supported treadmill training (BWSTT) or functional electrical stimulation (FES) (39). Low muscle quality and strength also decrease after SCI, and patients with a high injury plane showed a greater decrease (40, 41). Infiltration of intramuscular fat is also associated with muscle loss; however, FES and resistance training can improve these two symptoms to a certain extent (41). In addition to exercise training, the use of some drugs, such as testosterone and B2-adrenergic agonists, have been shown to improve muscle loss (42). In addition, in patients with SCI, the vastus lateralis show higher fibrosis compared with healthy controls, and muscle training, such as resistance training, could reduce muscle fibrosis (43, 44). However, there is still a lack of research on the effects of training at the genetic level.
Normal muscle tissue has a degree of regeneration that resists muscle loss, but beyond this level (especially after disease or aging), this function does not compensate for loss (45). The regeneration of muscle tissue is a coordinated result of multiple cells, among which satellite cells, FAP cells, have shown important regeneration potential in recent years (46, 47). FAPs plays an important role in muscle homeostasis and regeneration (33, 34). In a muscle where no injury has occurred, FAPs are in a quiescent condition. Furthermore, they are the predominant mononuclear cell population and an essential mesenchymal progenitor for supporting function (48). After muscle injury, FAPs will wake related cell populations to participant in regeneration. In an animal model experiment, researchers used Pdgfrα knock-out mice to simulate depletion of skeletal muscle, which showed a muscle regeneration delay, significant regenerative deficit, and a decrease in the muscle index (49). Furthermore, in sarcopenia-related research, when the depletion of FAPs occurs, it causes a reduction in muscle strength and weight (50). This is why FAPs was selected as a sub-cell cluster to conduct further analysis in this study.
In our study, we focus on the mechanism of SCI and sarcopenia in the muscle tissue. After validation, we obtained three hub genes: DCN, FSTL1, and COL12A1. According to the clinic information supplied by the datasets, we performed a correlation analysis between the hub genes and clinical information. The result showed age and BMI had a positive relationship with the hub genes.
DCN is an avid collagen-binding protein and, and in many skeleton and muscle diseases, DCN plays different roles due to its role in the pericellular matrix and extracellular matrix (51–53). DCN has been widely shown to regulate cartilage degeneration, bone matrix in osteoporosis, tendon architecture, and functional activity (54–56). In contrast, research on the association of DCN in SCI and sarcopenia is lacking. In SCI, few researchers focus on the muscular changes; most assess the transformation of the central nervous system. Because DCN can suppress scarring effects and inhibit fibrogenesis, researchers want to use it to reverse the scarring response after nervous system injury (57, 58). As previously mentioned, skeletal muscle loss in SCI patients can result in fibrosis; however, consistent rehabilitation training can alleviate this effect. It could be assumed that DCN may work in this training process and suppress fibrosis. In out results, DCN showed significant changes after training. Furthermore, in the KEGG enrichment analysis of DCN, it was found to play a crucial role in the TGF-beta signaling pathway, which may promote fibrosis, and if the pathway is suppressed, it will prevent fibrosis and improve regeneration (59–61). When the expression of DCN is elevated, the effects on the TGF-beta signaling pathway decline. Furthermore, DCN also participates in process of muscle regeneration. According to the single-cell RNA-seq analysis, the expression of DCN showed a significantly difference in the FAPs cell population, with an increasing trend in pseudotime trajectory. In the GO enrichment, analysis the result showed that DCN is related to mitochondrial function. There are few studies to show a correlation between DCN and mitochondrial function in muscle tissue; however, mitochondrial function, including maximal mitochondrial ATP production, dysregulation of mitochondrial dynamics, and mitochondrial proteolysis and mitophagy, may affect the muscle mass, function, and quality (62). In the skeletal muscle of patients with SCI, a mitochondrial oxygenation capability deficit has been observed, which is related to components such as succinate dehydrogenase decrease; rehabilitation training could increase mitochondrial function (63).
FSTL1 is a secreted factor from skeletal muscle; it participates in revascularization during ischemia and is related to the function of angiogenesis in skeletal muscle (64). For the capability of normal tissue, blood supply is critical, which is determined by vascular formation (65). In muscle regeneration, when mesenchymal or other stem cell’s function, the decisive factor is whether there is existing vasculature (66). In addition, after exercise, FSTL1 expression is elevate in the serum (67). However, evidence is lacking regarding a correlation between FSTL1 and indicies of muscle, including mass, quality, and performance. In the GO enrichment analysis, FSTL1 was shown to be related to the BMP signaling pathway, endothelial cell differentiation, and development. Research shows that BMP signaling controls muscle mass by the Smad family, and perturbed BMP signaling would lead to muscle loss (68, 69). Furthermore, BMP signaling relies on satellite cells to promote muscle regeneration. The links between FSLT1, BMP signaling, and satellite cells are elusive (70). There were no outcomes in the KEGG analysis of FSTL1; however, some studies have shown that it is related to Samd2/3 signaling and integrin β3/Wnt signaling (71, 72).
COL12A1 belongs to the family of fibril-associated collagens with interrupted triple helical domains and regulates cell communication, tissue repair, and regeneration (73, 74). COL12A1 is widely expressed in mesenchymal tissues in the embryo, which is consistent with our results in the GO enrichment analysis (75). In COL12A1-related myopathic Ehlers–Danlos syndrome, patients usually present with global muscle weakness and atrophy (76). In the KEGG analysis, COL12A1 was related to protein digestion and absorption. The skeletal muscle stores approximately 75% of the protein in the whole body, so when muscle loss or other sarcopenic situations begin, the balance between protein production and consumption is lost (77, 78). In patients with SCI, they may deplete more protein; therefore, it is vital for them to pay more attention to protein intake, and muscle loss is one of the recommendation standards for increasing protein intake (78). Thus, it is a worth exploring markers of skeletal muscle-related dysfunction.
In the immune infiltration landscape in SCI muscle tissue, there is no difference when comparing the time point of day 2 and day 5. We deem that the distance between the two time points is too short to make the difference significant. In sarcopenia datasets, compared with healthy control, the B cell naïve and B cell memory change significantly. Research shows that the capacity of B cells relates to skeletal regeneration and muscle strength (79). DCN and FSTL1 significantly correlated with M1 and M2 macrophages. Myogenesis and myotube fusion correlate with M2-biase macrophages and M1 macrophages convert to M2 macrophages in skeletal muscle regeneration (80). Furthermore, the infiltration of neutrophils will lead to fibrosis in muscle fiber regeneration, which significantly correlates with DCN and FSTL1 (79).
There are some limitations in our study should be acknowledged. Firstly, in microarray datasets or bulk-RNA seq datasets, patients with SCI are not formally diagnosed with sarcopenia. In that situation, we used the related sarcopenia datasets. Secondly, clinical information is lacking in many datasets, which made it hard for us to find the exact evidence and correlations between the hub genes and clinical features. Even although we still used limited clinical information to perform the correlation, a larger sample is needed to validate its value. In our future research, we will validate these hub genes in vitro and clearly explore the mechanism between SCI and sarcopenia in muscle tissue.
5 Conclusion
In our study, we identified three hub genes (DCN, FSTL1, and COL12A1) by bioinformatics analysis of GEO datasets. In patients with SCI, these hub genes show the underlying mechanism of the development of sarcopenia and could be a diagnostic and prognostic marker. We also provide the predictive values in a rehabilitation view to understand the clinical values of the hub genes, which will help in clinical decision-making. Furthermore, we present the immune change in SCI with sarcopenia. In general, the new sights provided by our study may promote the exploration of the comorbidity of SCI and sarcopenia. DCN, FSTL1, and COL12A1 are new candidate biomarkers for the comorbidity of SCI and sarcopenia.
Data availability statement
Information for existing publicly accessible datasets is contained within the article.
Ethics statement
Ethical approval was not required for the studies involving humans because all data comes from public databases and is not collected by us. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements because all data comes from public databases and is not collected by us. Ethical approval was not required for the study involving animals in accordance with the local legislation and institutional requirements because all data comes from public databases and is not collected by us.
Author contributions
BW: Formal analysis, Methodology, Software, Writing – original draft, Writing – review & editing. XY: Writing – review & editing. CL: Formal analysis, Writing – review & editing. RY: Writing – review & editing. TS: Writing – original draft. YY: Conceptualization, Funding acquisition, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the Zhou Mouwang Expert Workstation (202105AF150069).
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.2024.1373605/full#supplementary-material
SUPPLEMENTARY FIGURE S1 | The workflow of work.
References
1. Spinal cord injury (SCI) 2016 facts and figures at a glance. J Spinal Cord Med. (2016) 39:493–4. doi: 10.1080/10790268.2016.1210925
2. Jiang, B, Sun, D, Sun, H, Ru, X, Liu, H, Ge, S, et al. Prevalence, incidence, and external causes of traumatic spinal cord injury in China: a nationally representative cross-sectional survey. Front Neurol. (2021) 12:784647. doi: 10.3389/fneur.2021.784647
3. Zhou, H, Lou, Y, Chen, L, Kang, Y, Liu, L, Cai, Z, et al. Epidemiological and clinical features, treatment status, and economic burden of traumatic spinal cord injury in China: a hospital-based retrospective study. Neural Regen Res. (2024) 19:1126–32. doi: 10.4103/1673-5374.382257
4. Marzetti, E . Musculoskeletal aging and sarcopenia in the elderly. Int J Mol Sci. (2022) 23:2808. doi: 10.3390/ijms23052808
5. Hu, X, Xu, W, Ren, Y, Wang, Z, He, X, Huang, R, et al. Spinal cord injury: molecular mechanisms and therapeutic interventions. Signal Transduct Target Ther. (2023) 8:245. doi: 10.1038/s41392-023-01477-6
6. Fu, J, Wang, H, Deng, L, and Li, J. Exercise training promotes functional recovery after spinal cord injury. Neural Plast. (2016) 2016:1–7. doi: 10.1155/2016/4039580
7. Invernizzi, M, de Sire, A, and Fusco, N. Rethinking the clinical management of volumetric muscle loss in patients with spinal cord injury: synergy among nutritional supplementation, pharmacotherapy, and rehabilitation. Curr Opin Pharmacol. (2021) 57:132–9. doi: 10.1016/j.coph.2021.02.003
8. Gordon, T, and Tyreman, N. Electrical stimulation to promote muscle and motor unit force and endurance after spinal cord injury. J Physiol. (2023) 601:1449–66. doi: 10.1113/JP283972
9. Wiyanad, A, Amatachaya, P, Sooknuan, T, Somboonporn, C, Thaweewannakij, T, Saengsuwan, J, et al. The use of simple muscle strength tests to reflect body compositions among individuals with spinal cord injury. Spinal Cord. (2022) 60:99–105. doi: 10.1038/s41393-021-00650-4
10. Cruz-Jentoft, AJ, Bahat, G, Bauer, J, Boirie, Y, Bruyère, O, Cederholm, T, et al. Sarcopenia: revised European consensus on definition and diagnosis. Age Ageing. (2019) 48:16–31. doi: 10.1093/ageing/afy169
11. Dionyssiotis, Y, Skarantavos, G, Petropoulou, K, Galanos, A, Rapidi, CA, and Lyritis, GP. Application of current sarcopenia definitions in spinal cord injury. J Musculoskelet Neuronal Interact. (2019) 19:21–9.
12. Ishimoto, R, Mutsuzaki, H, Shimizu, Y, Kishimoto, H, Takeuchi, R, and Hada, Y. Prevalence of Sarcopenic obesity and factors influencing body composition in persons with spinal cord injury in Japan. Nutrients. (2023) 15:473. doi: 10.3390/nu15020473
13. Rodriguez, G, Berri, M, Lin, P, Kamdar, N, Mahmoudi, E, and Peterson, MD. Musculoskeletal morbidity following spinal cord injury: a longitudinal cohort study of privately-insured beneficiaries. Bone. (2021) 142:115700. doi: 10.1016/j.bone.2020.115700
14. Reich, KA, Chen, YW, Thompson, PD, Hoffman, EP, and Clarkson, PM. Forty-eight hours of unloading and 24 h of reloading lead to changes in global gene expression patterns related to ubiquitination and oxidative stress in humans. J Appl Physiol. (1985) 109:1404–15. doi: 10.1152/japplphysiol.00444.2010
15. Migliavacca, E, Tay, SKH, Patel, HP, Sonntag, T, Civiletto, G, McFarlane, C, et al. Mitochondrial oxidative capacity and NAD(+) biosynthesis are reduced in human sarcopenia across ethnicities. Nat Commun. (2019) 10:5808. doi: 10.1038/s41467-019-13694-1
16. Hangelbroek, RW, Fazelzadeh, P, Tieland, M, Boekschoten, MV, Hooiveld, GJ, van Duynhoven, JP, et al. Expression of protocadherin gamma in skeletal muscle tissue is associated with age and muscle weakness. J Cachexia Sarcopenia Muscle. (2016) 7:604–14. doi: 10.1002/jcsm.12099
17. Petrie, MA, Sharma, A, Taylor, EB, Suneja, M, and Shields, RK. Impact of short-and long-term electrically induced muscle exercise on gene signaling pathways, gene expression, and PGC1a methylation in men with spinal cord injury. Physiol Genomics. (2020) 52:71–80. doi: 10.1152/physiolgenomics.00064.2019
18. Oprescu, SN, Yue, F, Qiu, J, Brito, LF, and Kuang, S. Temporal dynamics and heterogeneity of cell populations during skeletal muscle regeneration. iScience. (2020) 23:100993. doi: 10.1016/j.isci.2020.100993
19. Ritchie, ME, Phipson, B, Wu, D, Hu, Y, Law, CW, Shi, W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. (2015) 43:e47. doi: 10.1093/nar/gkv007
20. Leek, JT, Johnson, WE, Parker, HS, Jaffe, AE, and Storey, JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. (2012) 28:882–3. doi: 10.1093/bioinformatics/bts034
21. Love, MI, Huber, W, and Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. (2014) 15:550. doi: 10.1186/s13059-014-0550-8
22. Yu, G, Wang, LG, Han, Y, and He, QY. Cluster profiler: an R package for comparing biological themes among gene clusters. OMICS. (2012) 16:284–7. doi: 10.1089/omi.2011.0118
23. Chen, H, and Boutros, PC. Venn diagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinformatics. (2011) 12:35. doi: 10.1186/1471-2105-12-35
24. Szklarczyk, D, Kirsch, R, Koutrouli, M, Nastou, K, Mehryary, F, Hachilif, R, et al. The STRING database in 2023: protein-protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Res. (2023) 51:D638–46. doi: 10.1093/nar/gkac1000
25. Shannon, P, Markiel, A, Ozier, O, Baliga, NS, Wang, JT, Ramage, D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498–504. doi: 10.1101/gr.1239303
26. Chin, CH, Chen, SH, Wu, HH, Ho, CW, Ko, MT, and Lin, CY. Cyto Hubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. (2014) 8:S11. doi: 10.1186/1752-0509-8-S4-S11
27. Robin, X, Turck, N, Hainard, A, Tiberti, N, Lisacek, F, Sanchez, JC, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinform. (2011) 12:77. doi: 10.1186/1471-2105-12-77
28. Sachs, MC . Plot ROC: a tool for plotting ROC curves. J Stat Softw. (2017) 79:2. doi: 10.18637/jss.v079.c02
29. Newman, AM, Liu, CL, Green, MR, Gentles, AJ, Feng, W, Xu, Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. (2015) 12:453–7. doi: 10.1038/nmeth.3337
31. Hao, Y, Hao, S, Andersen-Nissen, E, Mauck, WM 3rd, Zheng, S, Butler, A, et al. Integrated analysis of multimodal single-cell data. Cell. (2021) 184:3573–87.e29. doi: 10.1016/j.cell.2021.04.048
32. Trapnell, C, Cacchiarelli, D, Grimsby, J, Pokharel, P, Li, S, Morse, M, et al. The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells. Nat Biotechnol. (2014) 32:381–6. doi: 10.1038/nbt.2859
33. Theret, M, Rossi, FMV, and Contreras, O. Evolving roles of muscle-resident fibro-Adipogenic progenitors in health, regeneration, neuromuscular disorders, and aging. Front Physiol. (2021) 12:673404. doi: 10.3389/fphys.2021.673404
34. Uezumi, A, Fukada, S, Yamamoto, N, Takeda, S, and Tsuchida, K. Mesenchymal progenitors distinct from satellite cells contribute to ectopic fat cell formation in skeletal muscle. Nat Cell Biol. (2010) 12:143–52. doi: 10.1038/ncb2014
35. Qin, W, Bauman, WA, and Cardozo, C. Bone and muscle loss after spinal cord injury: organ interactions. Ann N Y Acad Sci. (2010) 1211:66–84. doi: 10.1111/j.1749-6632.2010.05806.x
36. Leone, GE, Shields, DC, Haque, A, and Banik, NL. Rehabilitation: neurogenic bone loss after spinal cord injury. Biomedicines. (2023) 11:2581. doi: 10.3390/biomedicines11092581
37. Cederholm, T, Barazzoni, R, Austin, P, Ballmer, P, Biolo, G, Bischoff, SC, et al. ESPEN guidelines on definitions and terminology of clinical nutrition. Clin Nutr. (2017) 36:49–64. doi: 10.1016/j.clnu.2016.09.004
38. Megía García, A, Serrano-Muñoz, D, Taylor, J, Avendaño-Coy, J, and Gómez-Soriano, J. Transcutaneous spinal cord stimulation and motor rehabilitation in spinal cord injury: a systematic review. Neurorehabil Neural Repair. (2020) 34:3–12. doi: 10.1177/1545968319893298
39. Castro, MJ, Apple, DF Jr, Staron, RS, Campos, GE, and Dudley, GA. Influence of complete spinal cord injury on skeletal muscle within 6 mo of injury. J Appl Physiol. (1985) 86:350–8. doi: 10.1152/jappl.1999.86.1.350
40. Butler, JE, and Thomas, CK. Effects of sustained stimulation on the excitability of motoneurons innervating paralyzed and control muscles. J Appl Physiol. (1985) 94:567–75. doi: 10.1152/japplphysiol.01176.2001
41. Atkins, KD, and Bickel, CS. Effects of functional electrical stimulation on muscle health after spinal cord injury. Curr Opin Pharmacol. (2021) 60:226–31. doi: 10.1016/j.coph.2021.07.025
42. Otzel, DM, Kok, HJ, Graham, ZA, Barton, ER, and Yarrow, JF. Pharmacologic approaches to prevent skeletal muscle atrophy after spinal cord injury. Curr Opin Pharmacol. (2021) 60:193–9. doi: 10.1016/j.coph.2021.07.023
43. Yarar-Fisher, C, Bickel, CS, Kelly, NA, Stec, MJ, Windham, ST, McLain, AB, et al. Heightened TWEAK-NF-κB signaling and inflammation-associated fibrosis in paralyzed muscles of men with chronic spinal cord injury. Am J Physiol Endocrinol Metab. (2016) 310:E754–61. doi: 10.1152/ajpendo.00240.2015
44. Lindsay, A, Larson, AA, Verma, M, Ervasti, JM, and Lowe, DA. Isometric resistance training increases strength and alters histopathology of dystrophin-deficient mouse skeletal muscle. J Appl Physiol. (1985) 126:363–75. doi: 10.1152/japplphysiol.00948.2018
45. Wallace, GQ, and McNally, EM. Mechanisms of muscle degeneration, regeneration, and repair in the muscular dystrophies. Annu Rev Physiol. (2009) 71:37–57. doi: 10.1146/annurev.physiol.010908.163216
46. Sousa-Victor, P, García-Prat, L, and Muñoz-Cánoves, P. Control of satellite cell function in muscle regeneration and its disruption in ageing. Nat Rev Mol Cell Biol. (2022) 23:204–26. doi: 10.1038/s41580-021-00421-2
47. Schmidt, M, Schüler, SC, Hüttner, SS, von Eyss, B, and von Maltzahn, J. Adult stem cells at work: regenerating skeletal muscle. Cell Mol Life Sci. (2019) 76:2559–70. doi: 10.1007/s00018-019-03093-6
48. Molina, T, Fabre, P, and Dumont, NA. Fibro-adipogenic progenitors in skeletal muscle homeostasis, regeneration and diseases. Open Biol. (2021) 11:210110. doi: 10.1098/rsob.210110
49. Wosczyna, MN, Konishi, CT, Perez Carbajal, EE, Wang, TT, Walsh, RA, Gan, Q, et al. Mesenchymal stromal cells are required for regeneration and homeostatic maintenance of skeletal muscle. Cell Rep. (2019) 27:2029–35.e5. doi: 10.1016/j.celrep.2019.04.074
50. Uezumi, A, Ikemoto-Uezumi, M, Zhou, H, Kurosawa, T, Yoshimoto, Y, Nakatani, M, et al. Mesenchymal bmp 3b expression maintains skeletal muscle integrity and decreases in age-related sarcopenia. J Clin Invest. (2021) 131:e139617. doi: 10.1172/JCI139617
51. Chen, S, and Birk, DE. The regulatory roles of small leucine-rich proteoglycans in extracellular matrix assembly. FEBS J. (2013) 280:2120–37. doi: 10.1111/febs.12136
52. Zhang, G, Chen, S, Goldoni, S, Calder, BW, Simpson, HC, Owens, RT, et al. Genetic evidence for the coordinated regulation of collagen fibrillogenesis in the cornea by decorin and biglycan. J Biol Chem. (2009) 284:8888–97. doi: 10.1074/jbc.M806590200
53. Chery, DR, Han, B, Zhou, Y, Wang, C, Adams, SM, Chandrasekaran, P, et al. Decorin regulates cartilage pericellular matrix micromechanobiology. Matrix Biol. (2021) 96:1–17. doi: 10.1016/j.matbio.2020.11.002
54. Han, B, Li, Q, Wang, C, Chandrasekaran, P, Zhou, Y, Qin, L, et al. Differentiated activities of decorin and biglycan in the progression of post-traumatic osteoarthritis. Osteoarthr Cartil. (2021) 29:1181–92. doi: 10.1016/j.joca.2021.03.019
55. Dunkman, AA, Buckley, MR, Mienaltowski, MJ, Adams, SM, Thomas, SJ, Satchell, L, et al. Decorin expression is important for age-related changes in tendon structure and mechanical properties. Matrix Biol. (2013) 32:3–13. doi: 10.1016/j.matbio.2012.11.005
56. Licini, C, Montalbano, G, Ciapetti, G, Cerqueni, G, Vitale-Brovarone, C, and Mattioli-Belmonte, M. Analysis of multiple protein detection methods in human osteoporotic bone extracellular matrix: from literature to practice. Bone. (2020) 137:115363. doi: 10.1016/j.bone.2020.115363
57. Esmaeili, M, Berry, M, Logan, A, and Ahmed, Z. Decorin treatment of spinal cord injury. Neural Regen Res. (2014) 9:1653–6. doi: 10.4103/1673-5374.141797
58. Logan, A, Baird, A, and Berry, M. Decorin attenuates gliotic scar formation in the rat cerebral hemisphere. Exp Neurol. (1999) 159:504–10. doi: 10.1006/exnr.1999.7180
59. Bernasconi, P, Di Blasi, C, Mora, M, Morandi, L, Galbiati, S, Confalonieri, P, et al. Transforming growth factor-beta 1 and fibrosis in congenital muscular dystrophies. Neuromuscul Disord. (1999) 9:28–33. doi: 10.1016/S0960-8966(98)00093-5
60. Cohn, RD, van Erp, C, Habashi, JP, Soleimani, AA, Klein, EC, Lisi, MT, et al. Angiotensin II type 1 receptor blockade attenuates TGF-beta-induced failure of muscle regeneration in multiple myopathic states. Nat Med. (2007) 13:204–10. doi: 10.1038/nm1536
61. Lefaucheur, JP, and Sébille, A. Muscle regeneration following injury can be modified in vivo by immune neutralization of basic fibroblast growth factor, transforming growth factor beta 1 or insulin-like growth factor I. J Neuroimmunol. (1995) 57:85–91. doi: 10.1016/0165-5728(94)00166-L
62. Gan, Z, Fu, T, Kelly, DP, and Vega, RB. Skeletal muscle mitochondrial remodeling in exercise and diseases. Cell Res. (2018) 28:969–80. doi: 10.1038/s41422-018-0078-7
63. O'Brien, LC, and Gorgey, AS. Skeletal muscle mitochondrial health and spinal cord injury. World J Orthop. (2016) 7:628–37. doi: 10.5312/wjo.v7.i10.628
64. Ouchi, N, Oshima, Y, Ohashi, K, Higuchi, A, Ikegami, C, Izumiya, Y, et al. Follistatin-like 1, a secreted muscle protein, promotes endothelial cell function and revascularization in ischemic tissue through a nitric-oxide synthase-dependent mechanism. J Biol Chem. (2008) 283:32802–11. doi: 10.1074/jbc.M803440200
65. Nomi, M, Atala, A, Coppi, PD, and Soker, S. Principals of neovascularization for tissue engineering. Mol Asp Med. (2002) 23:463–83. doi: 10.1016/S0098-2997(02)00008-0
66. Kasper, G, Dankert, N, Tuischer, J, Hoeft, M, Gaber, T, Glaeser, JD, et al. Mesenchymal stem cells regulate angiogenesis according to their mechanical environment. Stem Cells. (2007) 25:903–10. doi: 10.1634/stemcells.2006-0432
67. Görgens, SW, Raschke, S, Holven, KB, Jensen, J, Eckardt, K, and Eckel, J. Regulation of follistatin-like protein 1 expression and secretion in primary human skeletal muscle cells. Arch Physiol Biochem. (2013) 119:75–80. doi: 10.3109/13813455.2013.768270
68. Sartori, R, Schirwis, E, Blaauw, B, Bortolanza, S, Zhao, J, Enzo, E, et al. BMP signaling controls muscle mass. Nat Genet. (2013) 45:1309–18. doi: 10.1038/ng.2772
69. Sartori, R, Hagg, A, Zampieri, S, Armani, A, Winbanks, CE, Viana, LR, et al. Perturbed BMP signaling and denervation promote muscle wasting in cancer cachexia. Sci Transl Med. (2021) 13:eaay9592. doi: 10.1126/scitranslmed.aay9592
70. Stantzou, A, Schirwis, E, Swist, S, Alonso-Martin, S, Polydorou, I, Zarrouki, F, et al. BMP signaling regulates satellite cell-dependent postnatal muscle growth. Development. (2017) 144:2737–47. doi: 10.1242/dev.144089
71. Xi, Y, Hao, M, Liang, Q, Li, Y, Gong, DW, and Tian, Z. Dynamic resistance exercise increases skeletal muscle-derived FSTL1 inducing cardiac angiogenesis via DIP2A-Smad 2/3 in rats following myocardial infarction. J Sport Health Sci. (2021) 10:594–603. doi: 10.1016/j.jshs.2020.11.010
72. Cheng, S, Huang, Y, Lou, C, He, Y, Zhang, Y, and Zhang, Q. FSTL1 enhances chemoresistance and maintains stemness in breast cancer cells via integrin β3/Wnt signaling under mi R-137 regulation. Cancer Biol Ther. (2019) 20:328–37. doi: 10.1080/15384047.2018.1529101
73. Izu, Y, and Birk, DE. Collagen XII mediated cellular and extracellular mechanisms in development, regeneration, and disease. Front Cell Dev Biol. (2023) 11:1129000. doi: 10.3389/fcell.2023.1129000
74. Zou, Y, Zwolanek, D, Izu, Y, Gandhy, S, Schreiber, G, Brockmann, K, et al. Recessive and dominant mutations in COL12A1 cause a novel EDS/myopathy overlap syndrome in humans and mice. Hum Mol Genet. (2014) 23:2339–52. doi: 10.1093/hmg/ddt627
75. Chiquet, M, Birk, DE, Bönnemann, CG, and Koch, M. Collagen XII: protecting bone and muscle integrity by organizing collagen fibrils. Int J Biochem Cell Biol. (2014) 53:51–4. doi: 10.1016/j.biocel.2014.04.020
76. Mohassel, P, Liewluck, T, Hu, Y, Ezzo, D, Ogata, T, Saade, D, et al. Dominant collagen XII mutations cause a distal myopathy. Ann Clin Transl Neurol. (2019) 6:1980–8. doi: 10.1002/acn3.50882
77. Argilés, JM, Campos, N, Lopez-Pedrosa, JM, Rueda, R, and Rodriguez-Mañas, L. Skeletal muscle regulates metabolism via Interorgan crosstalk: roles in health and disease. J Am Med Dir Assoc. (2016) 17:789–96. doi: 10.1016/j.jamda.2016.04.019
78. Chen, L, Ming, J, Chen, T, Hébert, JR, Sun, P, Zhang, L, et al. Association between dietary inflammatory index score and muscle mass and strength in older adults: a study from National Health and nutrition examination survey (NHANES) 1999-2002. Eur J Nutr. (2022) 61:4077–89. doi: 10.1007/s00394-022-02941-9
79. Zhang, X, Li, H, He, M, Wang, J, Wu, Y, and Li, Y. Immune system and sarcopenia: presented relationship and future perspective. Exp Gerontol. (2022) 164:111823. doi: 10.1016/j.exger.2022.111823
Keywords: spinal cord injury, sarcopenia, rehabilitation training, skeletal muscle, genomics
Citation: Wang B, Yang X, Li C, Yang R, Sun T and Yin Y (2024) The shared molecular mechanism of spinal cord injury and sarcopenia: a comprehensive genomics analysis. Front. Neurol. 15:1373605. doi: 10.3389/fneur.2024.1373605
Edited by:
Mike Modo, University of Pittsburgh, United StatesReviewed by:
Syoichi Tashiro, Keio University, JapanIrma Ruslina Defi, Padjadjaran University, Indonesia
Copyright © 2024 Wang, Yang, Li, Yang, Sun and Yin. 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: Yong Yin, yyinpmr@126.com