Identification of key genes associated with muscle atrophy after spinal cord injury and experimental verification in rats

Identification of key genes associated with muscle atrophy after spinal cord injury and experimental verification in rats | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article Identification of key genes associated with muscle atrophy after spinal cord injury and experimental verification in rats Xiang Wang, Jianzhong Huo, Yimin Gao This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-5504847/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Spinal cord injury (SCI) is a neurological disease with high morbidity and mortality. Post-SCI muscle atrophy is a cascade response to SCI, and failure to actively prevent its occurrence severely affects patients' mobility and quality of life. Therefore, deeply exploring the correlation between muscle atrophy after SCI and the molecular regulation mechanism is of great significance. Download GSE21497 expression profile data from the gene expression omnibus(GEO)database. Perform weighted gene co-expression network analysis (WGCNA) on the obtained differentially expressed genes (DEGs).Subsequently, we performed functional and pathway enrichment analyses of key modules. Construct a protein-protein interaction (PPI) network and screen core genes. Finally, the results were verified by real-time polymerase chain reaction(PCR). A total of 1007 DEGs were obtained, including 533 upregulated genes and 474 downregulated genes. WGCNA analysis identified 161 turquoise modules of DEGs as key modules related to SCI. Functional enrichment analysis showed that their involvement in multiple processes and pathways. FOS and CCL2 may be involved in the molecular pathophysiology of muscle atrophy after SCI, serving as potential targets for diagnosis or treatment of SCI-related muscle atrophy. Biological sciences/Computational biology and bioinformatics Biological sciences/Neuroscience Health sciences/Biomarkers Health sciences/Medical research Differentially expressed genes Spinal cord injury Muscle atrophy Hub genes Bioinformatics analysis Weighted gene co-expression network analysis Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Introduction Spinal cord injury (SCI) refers to a neurological disorder characterized by damage to the spinal cord, resulting from direct or indirect trauma, leading to dysfunction of motor, sensory, and autonomic functions below the level of injury within the neuroanatomical distribution 1 . According to statistics, approximately 100,000 to 800,000 individuals in China experience SCI annually, while in the United States, each SCI patient incurs lifetime costs ranging from $ 1.1 million to $ 4.7 million, imposing significant burdens on both families and society, severely affecting the physiological, psychological health, and quality of life of patients 2 , 3 . With the rapid development of modern transportation, construction, and industry, the incidence and disability rate of SCI have been steadily increasing each year 4 .Mechanical compression following SCI can cause direct mechanical damage to spinal cord tissue, including local hemorrhage, edema, and ischemic reactions, a process referred to as primary injury 5 .This type of injury triggers a cascade of cellular and molecular responses, including inflammation, lipid peroxidation, free radical generation, disruption of ion channels, axonal demyelination, glial scar formation, necrosis, and programmed cell death, collectively termed secondary injury 6 . One of the most apparent consequences following SCI is a significant decrease in muscle mass due to partial or complete interruption of descending motor pathways 7 . Post-SCI muscle atrophy is a cascade response to SCI, and failure to actively prevent its occurrence severely affects patients' mobility and quality of life 8 . Following SCI, changes occur in spinal cord motor neurons and nerve fibers, leading to dysfunction of motor units, which results in decreased muscle fatigue resistance, reduced strength, and slowed motor conduction velocity. These alterations exacerbate functional decline post-SCI and impede motor recovery 9 .Given the potential modifiability of post-SCI muscle atrophy 10 , studying its mechanisms and providing effective treatment strategies have become recent focal points. With the development of bioinformatics and the widespread use of high-throughput gene chip technology in exploring molecular changes in disease progression, there is now a theoretical basis and research direction for understanding the pathogenesis of post-SCI muscle atrophy 11 . This facilitates the targeted selection of genes with research value, potentially saving significant research time and effort. This study used the post SCI muscle atrophy expression profile data (GSE21497) from the gene expression omnibus database (GEO) as the research object, screened differentially expressed genes (DEGs), and performed weighted gene co-expression network analysis (WGCNA) and functional enrichment analysis to construct a protein protein interaction (PPI) network and screen for key genes. These analyses aim to provide a theoretical basis for exploring the molecular mechanisms underlying the development and maintenance of muscle atrophy post-SCI, as well as identifying important therapeutic targets for treating muscle atrophy. Finally, the results were verified by real-time polymerase chain reaction (PCR). Materials and Methods Microarray data analysis Microarray and gene expression profile dataset GSE21497 were obtained from the freely accessible public database NCBI-GEO ( https://www.ncbi.nlm.nih.gov/geo/ ). The gene sequencing platform used was GPL570 [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array. This database employed a self-controlled method to collect samples on the 2nd and 5th days after SCI occurrence from a total of 10 patients (9 males, 1 female, with an average age of 44 years, including 6 cases of tetraplegia and 4 cases of paraplegia). Muscle biopsy samples from the lateral vastus muscle were obtained. The online tool GEO2R ( https://www.ncbi.nlm.nih.gov/geo/geo2r/ ) was utilized to perform DEGs analysis between the two days and five days post-SCI group. Data processing and selection of DEGs Using the online Sangerbox 3.0 ( http://vip.sangerbox.com/home.html ), we performed ID conversion and gene annotation processing on dataset GSE21497 12 . Differential analysis was conducted using the R software package limma (version 3.40.6) to obtain DEGs between the two groups 13 . Based on the acquired expression profile dataset, the data were log2-transformed, and multivariate linear regression was performed using the lmFit function. Subsequently, the eBays function was employed to obtain the expression matrix, resulting in the identification of significant differences for each gene. DEGs were selected using a threshold of P 1.5. The online Sangerbox 3.0 software was then utilized to input the original data in TXT format, facilitating online construction of DEGs and the generation of volcano plots and heatmaps. Weighted Gene Co-Expression Network Analysis(WGCNA) Performing weighted gene co-expression network analysis (WGCNA) using the R software package on the obtained DEGs 14 .Inputting the DEGs data for sample selection, wherein sample and grouping information were used to set a cut-Height threshold for sample exclusion, removing outliers, and constructing sample clustering. After preparing the data for sample clustering observations, determining the optimal soft threshold (β) to ensure the constructed network conforms to a scale-free topology. Constructing a co-expression network based on the optimal soft threshold and partitioning genes into different modules, represented by branch colors on the clustering tree, and plotting gene clustering trees and module eigengene clustering diagrams. Based on these results, calculating the correlation between gene modules and phenotypes, identifying trait-associated modules, and plotting phenotype-gene correlation heatmaps. Functional enrichment analysis Gene Ontology (GO) and Kyoto Encyclopedia of Genes Genomes (KEGG) pathway enrichment analyses were conducted on DEGs utilizing the clusterprofiler in R package (version 3.14.3, http://bioconductor.org) 15,16 .Usin g the Database for Annotation, Visualization, and Integrated Discovery (DAVID, https://david.ncifcrf.gov/home.jsp ) for online gene function annotation and signal pathway enrichment analysis 17 , including biological processes (BP), molecular functions (MF), and cellular components (CC), as well as KEGG pathway analysis. Genes count ≥ 2 and P < 0.05 was used as the threshold of significant enrichment results. Gene Set enrichment analysis(GSEA) To observe whether there are changes in the distribution of pathway genes in all gene expression data between different time point groups of SCI patients. We obtained GSEA version3.0 from the GSEA website ( http://software.broadinstitute.org/gsea/ index.jsp) and divided the samples into two groups based on different observation time points 18 . Subsequently, we downloaded the subset from the Molecular Signatures Database ( http://www.gsea-msigdb.org/gsea/downloads.jsp ) to assess relevant pathways and molecular mechanisms 19 . Based on gene expression profiles and phenotype grouping, a minimum gene set of 5 and a maximum gene set of 5000 were established, with 1000 permutations. A significance threshold of P < 0.05 and FDR < 0.25 was considered statistically significant. Construction of protein-protein interaction (PPI) networks and identification of hub genes Using the STRING online database ( https://string-db.org/ ), the selected DEGs were inputted to analyze and construct PPI networks. The obtained data were imported into Cytoscape software to build protein co-expression network diagrams for comprehensive bioinformatics analysis 20 . Using the CytoHubba plugin in Cytoscape, the top 30 genes were selected, and subsequently, the top 10 hub genes were determined using various algorithms including Degree, maximum clique centrality (MCC), maximum neighborhood component (MNC), and edge percolated component (EPC). CytoHubba comprehensively uses four algorithms, namely Degree, MCC, MNC, and EPC, to evaluate the global connectivity, local compactness, topological importance of nodes, and edge permeability respectively. Selecting multiple algorithms can avoid the bias of a single algorithm, and cross-validation with multiple algorithms can improve the robustness of the results. To ensure the reliability of the screening results, we adopt a two-step weighting strategy: algorithm weight assignment and comprehensive score calculation.The raw data of the identified hub genes were inputted into the online Sangerbox 3.0 software to find intersections between sets and to generate Venn diagrams. Organize and analyze the differential expression of hub genes at different stages. Conduct statistical analysis on the expression differences of hub genes in the presence of injuries of varying degrees at the same stage. Establishment of a spinal cord injury model Forty healthy, adult, female Sprague-Dawley (SD) rats of specific pathogen-free (SPF) grade were obtained from the Animal Experiment Center of Shanxi Medical University (Animal Use License Number: SYXK (Jin) 2019 0008), with body weights ranging from 210 to 250 g. The experimental procedures complied with ethical standards for animal care and use. This study was approved by the Animal Ethics Committee of Shanxi Bethune Hospital (Approval Number: SBQKL-2021-029). All methods were performed in accordance with ARRIVE guidelines and relevant and regulations. The rats were randomly assigned into two groups (n = 20) using a random number table: the sham operation group (Sham group) and the spinal cord injury group (SCI group). In the Sham group, the rats underwent laminectomy and removal of the spinous processes at the T 8 − 9 vertebrae without spinal cord injury. In the SCI group, spinal cord injury was induced using the Allen method. Validation of hub genes via real-time PCR After the model was successfully established, 10 rats from each group were sacrificed at 2 and 5 days post-operation. All rats were euthanized by cervical dislocation under deep anesthesia with 2% sodium pentobarbital (100mg/kg), followed by transcardial perfusion normal saline (4°C). The lateral thigh muscles of rats was immediately harvested and stored at -80°C.Total RNA was extracted using TRIzol reagent (Invitrogen), followed by reverse transcription to synthesize cDNA. PCR amplification was performed, and electrophoresis was used for product detection and data analysis. The PCR conditions were as follows: pre-denaturation at 95°C for 30 s, denaturation at 94°C for 1 min, annealing at 56°C for 1 min, and extension at 72°C for 90 s, with a total of 35 cycles, followed by a final extension at 72°C for 10 min. In this study, the internationally commonly used 2 −ΔΔCT method was employed to calculate the relative mRNA expression levels. Statistical Analysis Data analysis was conducted using SPSS 26.0 statistical software(IBM Corporation, Armonk, NY, USA). The Shapiro-Wilk test was used to assess normality. For data that followed a normal distribution, measurements are presented as mean ± SD. To control the family-wise error rate (FWER) arising from multiple comparisons between time points, we applied the Benjamini-Hochberg false discovery rate (FDR) correction. Adjusted p-values were calculated using the GraphPad Prism (version 9.5.0), with a significance threshold set at FDR < 0.05. Group differences were analyzed by one-way ANOVA followed by Tukey’s post hoc test for pairwise comparisons. Figures were created using Adobe Photoshop CS6 and GraphPad Prism 9.5.0. Statistical significance was set at P < 0.05. Results Normalization of data and analysis of differential gene expression GSE21497 was downloaded from the GEO database and analyzed using the online tool GEO2R. In the box plot, the black lines are nearly at the same level, indicating minimal deviation between the two (Fig. 1 A). The analysis using the Uniform Manifold Approximation and Projection (UMAP) algorithm reveals distinct clustering between the SCI day 2 group and the SCI day 5 group, indicating significant differences between the two groups(Fig. 1 B). After filtering with P 1.5, a total of 1007 differentially expressed genes (DEGs) were identified on the 5th day post-SCI, including 533 upregulated genes and 474 downregulated genes. Volcano plots and heatmaps were generated using the online Sangerbox 3.0 tool(Figs. 1 C and D ). Data preparation and sample clustering The data matrix consists of a total of 20 samples, with 1007 DEGs. Further processing was performed using the WGCNA package in R software to conduct sample clustering. The clustering tree was examined to identify and filter out any outliers or anomalies. Based on the visual results of the sample dendrogram, the hierarchical clustering of the samples was obtained(Fig. 1 E). Selecting the optimal soft threshold After completing the sample clustering observations, the optimal soft threshold (soft thresholding or power) was determined to ensure the constructed network conforms to a scale-free topology. A scatter plot was generated for power values, and the results (Fig. 2 A) showed a clear trend throughout the range from 1, but the trend became less apparent after 18. Therefore, 18 was selected as the optimal power value (with a Power R-squared greater than 0.8). The mean connectivity plot (Fig. 2 A) depicts the network connectivity under different soft threshold values, and the relationship between average connectivity and power values was examined. The connectivity stabilizes at a value of 18, indicating its suitability for power selection. Constructing co-expression network Based on the optimal soft threshold, a co-expression network was constructed, and genes were divided into different modules. A gene clustering tree (Fig. 2 B) can be plotted, where the upper part represents the hierarchical clustering dendrogram of genes, and the lower part represents gene modules. Correspondence between the upper and lower parts indicates that genes with closer distances (clustered into the same branch) are assigned to the same module. Each color represents a module, with the grey module considered as a collection of genes that cannot be assigned to any module. Module correlation clustering and analysis of phenotype-gene correlation Explore the correlation between various modules, calculate module eigengenes, and plot a network diagram of module eigengenes. Darker colors indicate higher correlation coefficients (Fig. 2 C). Based on the above results, a heatmap of phenotype-gene correlation was plotted(Fig. 2 D). DEGs were categorized into four modules: blue, brown, turquoise, and grey. The color blocks on the right represent modules, while the color bar in the top right corner represents the range of correlations. Darker colors in the heatmap indicate higher correlation, with red indicating positive correlation and green indicating negative correlation. The numbers in each cell represent correlation and significance. Based on the absolute value of correlation, the most relevant modules were selected. The results indicate that the turquoise module (0.80, p < 0.0001) is the most negatively and positively correlated module with SCI. Therefore, the turquoise module, which includes 161 DEGs, is considered the key module associated with SCI. The grey module contains all genes that were not involved in clustering and is therefore considered an invalid module, not suitable for further analysis. Functional enrichment analysis The 161 DEGs in the turquoise module were subjected to GO and KEGG pathway analysis using the DAVID tool. The GO analysis revealed enrichment of the DEGs in 147 biological processes (BP), 146 cellular components (CC), and 145 molecular functions (MF) within the turquoise module. In the BP analysis, the DEGs are significantly enriched in terms such as negative regulation of cellular process, cytosol and response to organic substance(Fig. 3 A). In the CC analysis, the DEGs are significantly enriched in terms such as cytosol, endomembrane system and extracellular region(Fig. 3 B). In the MF analysis, the DEGs are significantly enriched in terms such as extracellular matrix structural, organic anion transmembrane transporter activity constituent and extracellular matrix binding(Figs. 3 C and D ). KEGG analysis indicates the involvement of the turquoise module in various signaling pathways including the TNF signaling pathway, cellular senescence and regulation of actin cytoskeleton(Figs. 3 E and F ). Gene set enrichment analysis (GSEA) The Gene Set Enrichment Analysis (GSEA) examined whether there were changes in the distribution of pathway gene sets among different groups of SCI patients. The larger the |NES| value, the smaller the FDR value, indicating higher reliability of the analysis results. The results indicate significant enrichment and upregulation of genes associated with the peroxisome proliferators-activated receptors (PPARs) signaling pathway (NES = 1.6916, FDR = 0.2060) and the adherens junction signaling pathway (NES = 1.6957, FDR = 0.2485) in both groups(Figs. 4 A and B ). Conversely, genes related to the DNA replication signaling pathway (NES=-1.6676, FDR = 0.2274) show significant enrichment and downregulation in both groups(Fig. 4 C). PPI network construction and selection and analysis on hub genes Using the STRING database ( https://cn.string-db.org ), a PPI network was constructed for the enriched DEGs in the turquoise module, with a threshold set to 0.4. The results were exported for further analysis. The file containing node and edge information was imported into Cytoscape using the appropriate plugin, resulting in the PPI network diagram (Fig. 5 A). Subsequently, the top 30 hub genes were selected (Fig. 5 B). Utilizing the CytoHubba plugin in Cytoscape, hub genes were identified using four different calculation methods: Edge Percolated Component (EPC), Maximum Clique Centrality (MCC), Maximum Neighborhood Component (MNC), and Degree. This process aided in identifying the top 10 critical genes (Figs. 6 A-D). The top 10 gene sets obtained from the four different algorithms were input into Sangerbox 3.0 software to determine the number of intersections in the datasets (Fig. 7 A). Finally, FOS and CCL2 were identified as the hub genes for our screening. Based on the data of 10 samples from the gene expression profile GSE21497 database in this study, it was found that there were differences in the expressions of FOS and CCL2 at 2 days and 5 days after the injury. Specifically, compared with the situation at 2 days after the injury, the protein expressions of FOS and CCL2 both decreased at 5 days after the injury (Fig. 7 B).We analyzed the expression levels of FOS and CCL2 in the lower limb muscle tissues of 6 patients with tetraplegia (severe SCI) and 4 patients with paraplegia (moderate SCI). The results showed that at 2 days, the expression levels of FOS and CCL2 in the lower limb muscle tissues of the 6 tetraplegic patients were higher. By the 5th day, although the expression levels of FOS and CCL2 in their lower limb muscle tissues decreased to some extent, they still remained at a relatively high level (Figs. 7 C and D ). This indicates that the higher the expression levels of FOS and CCL2, the more severe the muscle atrophy in patients. Validation of hub genes To further validate the reliability of WGCNA results, PCR was used to reconfirm two Hub genes.The primers for gene expression were designed as follows: FOS: forward 5'-AACAGATCCGAGCAGCTTCTA-3', reverse5'-TTTTGAGCTTCAACCGGCATC-3' CCL2: forward 5'-TCTTGCTGCACGTCCTTTATT-3', reverse 5'-GCATCCAGTTTTTGTATGTGCC-3'. The expression of FOS and CCL2 genes in the SCI model was analyzed via PCR. The results of our PCR experiments suggest that the expressions of FOS mRNA and CCL2 mRNA peaked at 2 days after SCI, and then slowly decreased at 5 days after SCI, but still remained at a relatively high level. (Figs. 7 E and F ) ( P ≤ 0.05). Discussion The secondary neuropathic changes following SCI are the primary cause of muscle atrophy, which not only impacts the care management and daily activities of SCI patients but also increases the risk of secondary complications (such as osteoporosis, diabetes, and cardiovascular diseases) 21 , 22 .Therefore, preventing muscle atrophy is crucial for maintaining metabolic health and normal daily activities after SCI. The pathogenesis of muscle atrophy following SCI involves various factors, including signal transduction, immune response, electrical conduction, stimulation, and metabolism 23 . SCI causes loss of central regulation function in peripheral nerves below the damaged segment, resulting in sensory, motor, and autonomic dysfunction, muscle paralysis, and decreased muscle loading 24 . Impairment of the neuromuscular system and reduced integrity of the musculoskeletal system are critical features of SCI and significant barriers to motor function recovery 25 . However, the etiology and precise underlying mechanisms of muscle atrophy after SCI remain poorly understood. The aim of our study is to identify relevant genes, biological processes, and pathway pathways that may play crucial roles in the development of muscle atrophy after SCI, thereby providing theoretical insights into the pathogenesis and therapeutic targets of muscle atrophy following SCI. In this study, we downloaded the gene dataset GSE21497. After data normalization and standardization, we used the online Sangerbox 3.0 software for computation, resulting in 1007 differentially expressed genes (DEGs) on day 5 after SCI, including 533 upregulated genes and 474 downregulated genes. WGCNA analysis was performed on the data matrix containing 20 samples and 1007 genes with expression differences. WGCNA can identify sets of genes with similar expression patterns, elucidate the relationships between gene sets and sample phenotypes, construct regulatory networks among genes within gene sets, and identify key regulatory genes 26 . All 1007 DEGs obtained through the “limma” algorithm were clustered into four modules. It was found that the 161 genes in the turquoise module were most closely related to SCI. GO enrichment analysis involves examining whether genes in a gene set are enriched at different levels to understand their functions in biological processes(BP), cellular components(CC), and molecular functions(MF) 27 . The results of this study show that DEGs are significantly enriched in biological processes such as “negative regulation of cellular process, cytosol and response to organic substance”. DEGs are also significantly enriched in cellular components including “cytosol, endomembrane system and extracellular region”, as well as molecular functions such as “extracellular matrix structural, organic anion transmembrane transporter activity constituent and extracellular matrix binding”. KEGG enrichment analysis is used to analyze the enrichment of gene pathways and functions in a given gene set 28 . KEGG analysis reveals that the turquoise module is involved in various signaling pathways such as the “TNF signaling pathway, ellular senescence and regulation of actin cytoskeleton”. GSEA is used to evaluate the distribution trend of genes from a predefined gene set within a ranked gene expression dataset associated with phenotypes, thereby determining their contribution to the phenotype 29 . We utilized the GSEA method to analyze the results, which indicated that following SCI, genes associated with the peroxisome proliferators-activated receptors (PPARs) signaling pathway and adherens junction signaling pathway were significantly increased in both groups, with upregulated expression, while genes related to the DNA replication signaling pathway were significantly decreased in both groups, with downregulated expression. PPARs pharmacological activation can be considered as a multifaceted therapeutic target, exerting anti-inflammatory, antioxidant, antiexcitotoxic, and energy-generating effects in some neurologic and inflammation-related diseases 30 . The Adherens junction signaling pathway also demonstrates significant enrichment in comprehensive proteomic and phosphoproteomic analyses following SCI 31 , reducing the loss of tight junction and adherens junction proteins in human brain microvascular endothelial cells (HBMECs). This reduction can attenuate blood-spinal cord barrier (BSCB) permeability and degradation of tight junction molecules (such as P120, β-catenin, Occludin, and Claudin5), ultimately improving functional recovery in a traumatic rat model 32 . Studies have shown that enhancement of the peroxisome proliferator-activated receptor-γ coactivator-1α(PGC-1α) mitochondrial biogenesis pathway upregulates nuclear respiratory factors (NRF-1,2) and mitochondrial transcription factor A (Tfam), stimulating mitochondrial biogenesis 33 .This leads to increased mitochondrial DNA replication and gene transcription, reducing oxidative damage and providing protective effects against atrophy-related skeletal muscle degeneration. The PPI network illustrates the relationships among proteins and even genes within a given gene set that play a protective role against atrophy-related skeletal muscle degeneration 34 . It is one of the commonly used methods for identifying core genes. Cytoscape ( https://cytoscape.org/ ) is an excellent software for graphical visualization and analysis of networks and can be equipped with various functional plugins, making it widely used in PPI analysis 35 . Based on the PPI network and the CytoHubba plugin with four different algorithms, the top 10 key genes were identified. These top 10 genes from each algorithm were input into the online Sangerbox 3.0 software for analysis, resulting in a Venn diagram showing FOS and CCL2 as the hub genes we selected. Moreover, the analysis of the data from the gene expression profile GSE21497 shows that the higher the expression levels of FOS and CCL2 are, the more severe the muscle atrophy of the patients will be. It was found that there were differences in the expressions of FOS and CCL2 at 2 days and 5 days after the injury. The results of our PCR experiments suggest that the expressions of FOS mRNA and CCL2 mRNA peaked at 2 days after SCI, and then slowly decreased at 5 days after SCI, but still remained at a relatively high level. Therefore, this dynamic change indicates that inhibiting the expression of CCL2 in the early stage of SCI and upregulating FOS in the chronic stage (such as local injection of FOS overexpression vectors) may be effective methods to alleviate muscle atrophy after SCI. FOS belongs to a class of genes known as immediate early genes (IEGs). It encodes a nuclear phosphoprotein that serves as a transcription factor 36 . FOS plays a significant role in regulating cell growth, division, proliferation, differentiation, and even programmed cell death 37 . As the protein product of the c-fos gene, under normal conditions, c-fos expression is minimal within cells. But, to some extent, the number of FOS-immunoreactive neurons correlates positively with the intensity of the stimulus received, making FOS a marker of injurious stimuli 38 . Research has shown that SCI significantly affects c-fos expression in the paraventricular nucleus (PVN) and supraoptic nucleus (SON) of the hypothalamus in rats. Following SCI at 1, 6, 12, and 24 hours, there is a significant increase in c-fos expression in the PVN, which persists for more than 24 hours. At 12 hours post-SCI, there is a significant increase in c-fos expression in the SON, reaching its peak at 24 hours post-SCI 39 . Chemokines are small-molecule peptides secreted by immune or non-immune cells, induced by inflammatory cytokines, growth factors, and pathological stimuli. They play an important role in inflammatory diseases, wound healing, spinal cord injury, and other processes 40 .Chemokine C-C motif ligand 2 (CCL2), also known as monocyte chemoattractant protein-1 (MCP-1), is a detrimental factor released after neural injury, primarily involved in the local inflammatory response following SCI 41 . After SCI, chronic pain, circulatory disturbances, and mechanical compression lead to changes in neurons of the dorsal root ganglion (DRG). After peripheral nerve or ganglion injury, CCL2 expression is upregulated in DRG neurons 42 . Research has shown that CCL2 plays a crucial role in recruiting macrophages during the early stages of muscle injury. Mice lacking CCL2 exhibit incomplete recovery of gastrocnemius muscle mass, muscle fiber cross-sectional area, and gastrocnemius muscle contraction characteristics during disuse atrophy recovery. Additionally, mice lacking CCL2 demonstrate reduced turnover of skeletal muscle collagen and significantly decreased recruitment of macrophages to the gastrocnemius muscle, resulting in poor recovery of muscle size and function as well as abnormal collagen remodeling 43 . We verified the differential expression of FOS and CCL2 in the SCI animal model through PCR, confirming the reliability of the bioinformatics analysis. The potential role of FOS: As an activator protein-1(AP-1) transcription factor, it may promote the activity of the ubiquitin-proteasome system by activating the Atrogin-1/MuRF1 pathway, thus facilitating the degradation of muscle fibers 44 . The potential role of CCL2: Recruiting macrophages through the CCR2 receptor, inducing the release of pro-inflammatory factors (such as TNF-α and IL-6), and accelerating the apoptosis of muscle cells 45 . These analyses link the changes in gene expression with known pathological mechanisms, pointing the way for subsequent functional studies. In terms of drug research and development, the already marketed JNK inhibitor (such as CC-930) can indirectly inhibit the activity of FOS by blocking the formation of the AP-1 complex 46 . The CCL2 antagonist (Bindarit, a CCL2 synthesis inhibitor) can reduce the infiltration density of macrophages and improve neuromotor function 47 .In terms of clinical translation, the serum CCL2 level can be used as a biomarker to predict the effective rate of Bindarit treatment and guide patient stratification. Clinical trials of the synergistic use of FOS inhibitors and rehabilitation training can also be carried out. Based on the existing research, the pro-inflammatory and pro-atrophic mechanisms of FOS and CCL2 make them ideal targets for the treatment of muscle atrophy following SCI. For example, the small molecule inhibitors T-5224 (targeting FOS/AP-1) and RS504393 (targeting CCL2/CCR2) have demonstrated anti-inflammatory and anti-fibrotic effects in other disease models 48 , 49 . Future studies could verify whether these drugs alleviate muscle atrophy in the SCI model by inhibiting the MuRF1/MAFbx pathway or regulating macrophage polarization. Additionally, combining gene silencing techniques or existing neurotrophic therapies may produce synergistic effects. In conclusion, as hub genes, FOS and CCL2, through the closed-loop path of "network screening-mechanism analysis-drug verification-clinical transformation", provide powerful targets for the development of precise drugs against muscle atrophy. Limitations Our study has several limitations. Firstly, the GSE21497 dataset used in this study is from the GEO database. Although it provides valuable information on the gene expression of muscle atrophy after SCI, there are still some limitations such as sample heterogeneity, and the lack of clinical information. Secondly,the samples included in this study are concentrated within 5 days after injury, which may not fully reflect the molecular characteristics of muscle atrophy in the chronic phase. Therefore, there may be a bias in temporal dynamics during the sample selection process. Finally, we have not directly verified the quantitative relationship between the expression of FOS and CCL2 and the severity of muscle atrophy in patients. In the later stage, we will continue to conduct clinical research related to this project. Conclusions In summary, we utilized WGCNA to analyze the expression profile data from GSE21497, identifying two hub genes, FOS and CCL2, in muscle atrophy following SCI. FOS and CCL2 may be involved in the molecular pathophysiology of muscle atrophy after SCI, serving as potential targets for diagnosis or treatment of SCI-related muscle atrophy. Abbreviations AP-1,activator protein-1 BP,biological processes BSCB,blood-spinal cord barrier CC,cellular components CCL2,Chemokine C-C motif ligand 2 DAVID, database for annotation, visualization, and integrated discovery DRG,dorsal root ganglion DEGs, differentially expressed genes EPC, edge percolated component GEO, gene expression omnibus GO, gene ontology GSEA,gene set enrichment analysis HBMECs,human brain microvascular endothelial cells KEGG, kyoto encyclopedia of genes genomes MCC, maximum clique centrality MCP-1,monocyte chemoattractant protein-1 MF,molecular functions MNC,maximum neighborhood component PCR, polymerase chain reaction PPARs,peroxisome proliferators-activated receptors PPI, protein protein interaction PVN,paraventricular nucleus SCI, spinal cord injury SON,supraoptic nucleus WGCNA, weighted gene co-expression network analysis Declarations Author Contributions X.W. conceived and designed the study; Y.G. analyzed and interpreted the data and drafted the work. X.W. assisted in revising the manuscript. J.H. reviewed the manuscript and supervised the whole project. All authors have approved the submission of the final version of the manuscript and agree to take responsibility for all aspects of the work. Funding This study was funded by the Taiyuan City Science and Technology Plan Project (Program NO.202201), the Key Medical Research Project of Shanxi Province (Program NO.2023XM034), and the Basic Research Plan of Shanxi Province(Program NO. 202103021223410,202303021222381 and 202303021221246 ). Data availability The data used in this study are available in public databases. Some data and all scripts for conducting the analyses in this study can be shared upon reasonable request to the corresponding author. Disclosure The authors have no conflicts of interest to declare. Ethics statement This study was approved by the Animal Ethics Committee of Shanxi Bethune Hospital (Approval Number: SBQKL-2021-029).The study is also reported in accordance with ARRIVE guidelines. References Liu, Z., Chen, X., Chen, P. & Wang, L. Transcranial direct current stimulation attenuates the chronic pain of osteoarthritis in rats via reducing NMDAR2B expressions in the spinal cord. J Musculoskelet Neuronal Interact . 23 (4),498-505(2023). Fallah, N., Noonan,V.K., & Sharwood,L.N. Editorial: Epidemiology, evidence-based care, and outcomes in spinal cord injury. Front Neurol . 4 (15),1383757(2024). Zhou, H. 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 . 19 (5),1126-1133(2024). Eli,I., Lerner,D.P., & Ghogawala,Z. Acute traumatic spinal cord injury. Neurol Clin . 39 (2):471-488(2021). Anjum, A. et al. Spinal cord injury: pathophysiology, multimolecular interactions, and underlying recovery mechanisms. Int J Mol Sci . 21 (20),7533 (2020). Widerström-Noga,E. Neuropathic pain and spinal cord injury: management, phenotypes, and biomarkers. Drugs . 83 (11),1001-1025(2023). Otzel,D.M., Kok,H.J., Graham,Z.A., Barton,E.R. & Yarrow,J.F. Pharmacologic approaches to prevent skeletal muscle atrophy after spinal cord injury. Curr Opin Pharmacol . 10 (60),193-199(2021). Chandrasekaran, S., Davis,J., Bersch,I., Goldberg,G. & Gorgey,A.S. Electrical stimulation and denervated muscles after spinal cord injury. Neural Regen Res . 15 (8),1397-1407(2020). Xu, X.et al. Mechanism of skeletal muscle atrophy after spinal cord injury: a narrative review. Front Nutr . 3 (10),1099143(2023). Gorgey, A.S.et al.Mitochondrial health and muscle plasticity after spinal cord injury. Eur J Appl Physiol . 119 (2),315-331(2019). Potter, L.A.et al.Assessing the impact of boldine on the gastrocnemius using multiomics profiling at 7 and 28 days post-complete spinal cord injury in young male mice. Physiol Genomic s. 55 (7),297-313(2023). Lin, X. et al. Pancancer analyses reveal genomics and clinical characteristics of the SETDB1 in human tumors. J Oncol . 2022,6115878(2022). Ritchie, M.E.et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res . 43 (7),e47(2015). Langfelder, P. & Horvath,S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics . 29 (9),559(2008). Altermann,E. & Klaenhammer, T.R. PathwayVoyager: pathway mapping using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. BMC Genomics . 3(6),60(2005). Yu,G., Wang,L.G., Han,Y. & He,Q.Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS . 16 (5),284-7(2012). Sherman,B.T. et al. DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res . 50 (W1),W216-W221(2022). Subramanian,A.et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA . 102 (43),15545-50(2005). Liberzon,A. et al. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst . 1 (6),417-425(2015). Doncheva,N.T., Morris,J.H., Gorodkin,J. &Jensen,L.J. Cytoscape StringApp: network analysis and visualization of proteomics data. J Proteome Res . 18 (2),623-632(2019). Nemeth,C., Banik,N.L. & Haque,A. Disruption of neuromuscular junction following spinal cord injury and motor neuron diseases. Int J Mol Sci . 25 (6),3520(2024). Rosales-Antequera,C., Viscor,G. & Araneda,O.F. Inflammation and oxidative stress as common mechanisms of pulmonary, autonomic and musculoskeletal dysfunction after spinal cord injury. Biology (Basel) . 11 (4),550(2022). Wei,L., Cai,G., Jiang,L., Gao,L., Yang,Z. & Zhang,W. Identification of key pathways and RNAs associated with skeletal muscle atrophy after spinal cord injury. J Musculoskelet Neuronal Interact . 21(4):550-559(2021). Zhang,C. et al. Early electrical field stimulation prevents the loss of spinal cord anterior horn motoneurons and muscle atrophy following spinal cord injury. Neural Regen Res . 13 (5),869-876(2018). Dionyssiotis,Y., Mavrogenis,A.,Trovas,G., Skarantavos,G., Papathanasiou,J. & Papagelopoulos,P.Bone and soft tissue changes in patients with spinal cord injury and multiple sclerosis. Folia Med (Plovdiv) . 56 (4),237-44(2014). Xue,J.et al. The identification and validation of hub genes associated with acute myocardial infarction using weighted gene co-expression network analysis. J Cardiovasc Dev Dis . 9 (1),30(2022). Gaudet,P., Logie,C., Lovering,R.C., Kuiper,M., Lægreid,A. & Thomas PD. Gene Ontology representation for transcription factor functions. Biochim Biophys Acta Gene Regul Mech . 1864 (11-12),194752(2021). Chen,L., Zhang,Y.H., Wang,S., Zhang,Y., Huang,T. & Cai,Y.D.Prediction and analysis of essential genes using the enrichments of gene ontology and KEGG pathways. PLoS One . 12 (9),e0184129(2017). Zito,A.et al. Gene set enrichment analysis of interaction networks weighted by node centrality. Front Genet . 24 (12),577623(2021). Esposito,E. & Cuzzocrea,S. Targeting the peroxisome proliferator-activated receptors (PPARs) in spinal cord injury. Expert Opin Ther Targets . 15 (8),943-59(2011). Li,Z.et al. A proteomic and phosphoproteomic landscape of spinal cord injury. Neurosci Lett . 25 (814),137449(2023). Zhou,Y.et al.Phenylbutyrate prevents disruption of blood-spinal cord barrier by inhibiting endoplasmic reticulum stress after spinal cord injury. Am J Transl Res . 8 (4),1864-75(2016). Kang,C. & Li,J.L. Role of PGC-1α signaling in skeletal muscle health and disease. Ann NYAcad Sci . 1271 (1),110-7(2012). Wang,J., Liu,Y., Zhang,Y., Liu,B. & Wei,Z. Identification of potential miRNA-mRNA regulatory network in denervated muscular atrophy by bioinformatic analysis. Biomed Res Int . 28;2022,6042591(2022). Doncheva,N.T.et al. Cytoscape stringApp 2.0: Analysis and visualization of heterogeneous biological networks. J Proteome Res . 2023; 22 (2):637-646. Cordier,F. & Creytens,D. New kids on the block: FOS and FOSB gene. J Clin Pathol . 22 (2),637-646(2023). Choi,H.et al. Amyloid β-induced elevation of O-GlcNAcylated c-Fos promotes neuronal cell death. Aging Cell . 22 (2),637-646(2023). Zhang,D., Liu,J., Zhu,T. & Zhou,C. Identifying c-fos expression as a strategy to investigate the actions of general anesthetics on the cCentral nervous system. Curr Neuropharmacol . 20 (1),55-71(2022). Xu,Y., Zheng,Z., Ho,K.P. & Qian,Z. Effects of spinal cord injury on c-fos expression in hypothalamic paraventricular nucleus and supraoptic nucleus in rats. Brain Res . 1087 (1):175-9(2006). Cambier,S. et al. The chemokines CXCL8 and CXCL12: molecular and functional properties, role in disease and efforts towards pharmacological intervention. Cell Mol Immunol. 20 (3),217-251(2023). Cambier,S., Gouwy,M. & Proost,P. Bone marrow-derived monocytes drive the inflammatory microenvironment in local and remote regions after thoracic spinal cord injury. J Neurotrauma . 36 (6),937-949(2019). Dansereau,M.A.et al. Mechanistic insights into the role of the chemokine CCL2/CCR2 axis in dorsal root ganglia to peripheral inflammation and pain hypersensitivity. J Neuroinflammation . 18 (1),79(2021). Ferrara,P.J.et al. Global deletion of CCL2 has adverse impacts on recovery of skeletal muscle fiber size and function and is muscle specific. J Appl Physiol . 134 (4),923-932(2023). Sartori,R. et al. Smad2 and 3 transcription factors control muscle mass in adulthood. Am J Physiol Cell Physiol . 296 (6):C1248-57(2009). Hofmann,C. et al. Transient inhibition of translation improves cardiac function after ischemia/reperfusion by attenuating the inflammatory response. Circulation . 150(16):1248-1267(2024). Plantevin,K.V.et al. Discovery of CC-930, an orally active anti-fibrotic JNK inhibitor. Bioorg Med Chem Lett . 22 (3):1433-8(2012). Shi,K. et al. tPA mobilizes immune cells that exacerbate hemorrhagic transformation in stroke. Circ Res . 128 (1):62-7(2021). Zhou,L., Xue,C., Chen,Z., Jiang,W., He,S. & Zhang,X. c-Fos is a mechanosensor that regulates inflammatory responses and lung barrier dysfunction during ventilator-induced acute lung injury. BMC Pulm Med . 22 (1):9(2022). Raghu,H. et al. CCL2/CCR2, but not CCL5/CCR5, mediates monocyte recruitment, inflammation and cartilage destruction in osteoarthritis. Ann Rheum Dis . 76 (5):914-92(2017). Additional Declarations No competing interests reported. Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-5504847","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":442552835,"identity":"a5ebc872-4bfc-478c-8a40-2d36e8b84cf0","order_by":0,"name":"Xiang Wang","email":"","orcid":"","institution":"Third Hospital of Shanxi Medical University, Shanxi Academy of Medical Sciences, Tongji Shanxi Hospital","correspondingAuthor":false,"prefix":"","firstName":"Xiang","middleName":"","lastName":"Wang","suffix":""},{"id":442552836,"identity":"b1e435c4-e24c-4275-b1dc-ac98f05de4fa","order_by":1,"name":"Jianzhong Huo","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA20lEQVRIiWNgGAWjYBACAxBOgDIYPhjYyJGmhXFGQZoxUVrgDGaeD4cTCWoxl0g+UPCg5o69OXvv4dc2BswJDOyHj27Ap8VyRlqCQcKxZ8yWPefSrHMM2PIYeNLSbuB12I0cA4MEtsNsQIaZcY4BTzGDBI8ZEVr+HeYBa7EwkEhsIEpLYtthCSDD+DEDkE1Yy5lnCQaJfYcNDM6cMWPsMUgwZiPol+PJxwx/fDtsb3C8x/jDjz//5fjZDx/DqwUI2GBxwyYBJgkoBwHmBzDGByJUj4JRMApGwQgEAEWbS1K+UuNjAAAAAElFTkSuQmCC","orcid":"","institution":"Shanxi Medical University Second Affiliated Hospital","correspondingAuthor":true,"prefix":"","firstName":"Jianzhong","middleName":"","lastName":"Huo","suffix":""},{"id":442552837,"identity":"9b9915df-ebe7-4eb1-83b7-c454ecb8b3f9","order_by":2,"name":"Yimin Gao","email":"","orcid":"","institution":"Third Hospital of Shanxi Medical University, Shanxi Academy of Medical Sciences, Tongji Shanxi Hospital","correspondingAuthor":false,"prefix":"","firstName":"Yimin","middleName":"","lastName":"Gao","suffix":""}],"badges":[],"createdAt":"2024-11-22 13:23:12","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-5504847/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-5504847/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":80786602,"identity":"5a17f5f2-a9a9-45ea-af58-ec6a355be4ab","added_by":"auto","created_at":"2025-04-17 06:01:33","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":5509302,"visible":true,"origin":"","legend":"\u003cp\u003e(\u003cstrong\u003eA\u003c/strong\u003e) Box plots of normalized differential gene expression show that the black lines representing the boxes are generally at similar levels, indicating comparability between the two groups of patient data with minimal deviation between them. (\u003cstrong\u003eB\u003c/strong\u003e) Uniform Manifold Approximation and Projection (UMAP) algorithm analysis was conducted on all included samples, revealing distinct clustering within both groups and significant differences between the two groups.(\u003cstrong\u003eC\u003c/strong\u003e) Volcano plot (pink triangles represent upregulated genes, blue-green triangles represent downregulated genes, black dots represent genes with no statistical significance). (\u003cstrong\u003eD\u003c/strong\u003e) Heatmap (representing differentially expressed proteins between the two groups of patients, where red indicates high expression levels and green indicates low expression levels). (\u003cstrong\u003eE\u003c/strong\u003e) Perform clustering of differentially expressed genes, where each line represents a gene, and similar genes are clustered into branches.\u003c/p\u003e","description":"","filename":"Figure1.png","url":"https://assets-eu.researchsquare.com/files/rs-5504847/v1/b7a9e02a0cf7e177c97d4d0a.png"},{"id":80784141,"identity":"f6a99c0b-8e17-4023-bd4c-0cef04c96a04","added_by":"auto","created_at":"2025-04-17 05:33:54","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":2002971,"visible":true,"origin":"","legend":"\u003cp\u003e(\u003cstrong\u003eA\u003c/strong\u003e)Estimation plot of soft threshold for scale-free co-expression networks. (\u003cstrong\u003eB\u003c/strong\u003e)Dendrogram of hierarchical clustering for all differentially expressed genes.(\u003cstrong\u003eC\u003c/strong\u003e) Hierarchical clustering dendrogram of module eigengenes.(\u003cstrong\u003eD\u003c/strong\u003e) Heatmap showing the correlation of each module with the two groups.\u003c/p\u003e","description":"","filename":"Figure2.png","url":"https://assets-eu.researchsquare.com/files/rs-5504847/v1/111618c5c934e86900a7de8d.png"},{"id":80784149,"identity":"fb6a44d3-cd2f-4e24-a87b-7b4f6166b2b7","added_by":"auto","created_at":"2025-04-17 05:33:55","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":3809624,"visible":true,"origin":"","legend":"\u003cp\u003e(\u003cstrong\u003eA\u003c/strong\u003e) The bubble plot to the BP enrichment analysis results for DEGs. (\u003cstrong\u003eB\u003c/strong\u003e) The bubble plot to the CC enrichment analysis results for DEGs. (\u003cstrong\u003eC\u003c/strong\u003e) The MF enrichment analysis results for DEGs(Circos plot). (\u003cstrong\u003eD\u003c/strong\u003e) The MF enrichment analysis results for DEGs (Bubble plot). (\u003cstrong\u003eE\u003c/strong\u003e)The KEGG enrichment analysis results for DEGs (Circos plot ).(\u003cstrong\u003eF\u003c/strong\u003e) The KEGG enrichment analysis results for DEGs(Bubble plot).\u003c/p\u003e","description":"","filename":"Figure3.png","url":"https://assets-eu.researchsquare.com/files/rs-5504847/v1/117ac22e4700a893724b5dce.png"},{"id":80785739,"identity":"930b44db-db0c-44ff-ab33-11f4a827db64","added_by":"auto","created_at":"2025-04-17 05:41:54","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":745208,"visible":true,"origin":"","legend":"\u003cp\u003eGene set enrichment analysis plot (\u003cstrong\u003eA\u003c/strong\u003e) Enrichment plot for the PPARs signaling pathway.(\u003cstrong\u003eB\u003c/strong\u003e) Enrichment plot for the Adherens Junction signaling pathway.(\u003cstrong\u003eC\u003c/strong\u003e) Enrichment plot for the DNA Replication signaling pathway.\u003c/p\u003e","description":"","filename":"Figure4.png","url":"https://assets-eu.researchsquare.com/files/rs-5504847/v1/e3493db71afacecbb54911e9.png"},{"id":80784142,"identity":"d71a23f8-2da3-41ff-b8de-f1b6905d532d","added_by":"auto","created_at":"2025-04-17 05:33:54","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":2053818,"visible":true,"origin":"","legend":"\u003cp\u003e(\u003cstrong\u003eA\u003c/strong\u003e) Construct a Protein-protein interaction (PPI) network using the Differentially Expressed Genes (DEGs) enriched in the turquoise module, and generate a PPI network diagram.. (\u003cstrong\u003eB\u003c/strong\u003e) Screen the top 30 hub genes and create a network diagram.\u003c/p\u003e","description":"","filename":"Figure5.png","url":"https://assets-eu.researchsquare.com/files/rs-5504847/v1/cd0757472b5bcf65f2dcff41.png"},{"id":80784136,"identity":"0d57f843-4a16-4db8-83fd-443da38cf657","added_by":"auto","created_at":"2025-04-17 05:33:54","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":1303472,"visible":true,"origin":"","legend":"\u003cp\u003eNetwork diagram of the top 10 critical genes identified by the four algorithms. (\u003cstrong\u003eA\u003c/strong\u003e) Degree algorithm: evaluates the global connectivity of nodes.(\u003cstrong\u003eB\u003c/strong\u003e) Maximum Neighborhood Component (MCC) algorithm: focuses on core regulatory nodes.(\u003cstrong\u003eC\u003c/strong\u003e) Maximum Neighborhood Component (MNC) algorithm: emphasizes the importance within the module. (\u003cstrong\u003eD\u003c/strong\u003e) Edge Percolated Component (EPC) algorithm: identifies key nodes.\u003c/p\u003e","description":"","filename":"Figure6.png","url":"https://assets-eu.researchsquare.com/files/rs-5504847/v1/1fbfd458f196f0e39dce3f8c.png"},{"id":80784170,"identity":"a223381c-8b1e-4a84-8bde-a8784e82bb82","added_by":"auto","created_at":"2025-04-17 05:33:55","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":4876301,"visible":true,"origin":"","legend":"\u003cp\u003e(\u003cstrong\u003eA\u003c/strong\u003e):Venn diagram depicting the overlapping hub genes from different algorithms. (B):The expression changes of FOS and CCL2 at different time points were compared,\u003cem\u003e *\u003c/em\u003e\u003csup\u003e\u003cem\u003e#\u003c/em\u003e\u003c/sup\u003e\u003cem\u003eP\u003c/em\u003e\u0026lt;0.05 vs 2 days after SCI.(C) and (D): Comparison of the expression levels of FOS and CCL2 in the lower limb muscle tissues of patients with SCI of different degrees,\u003cem\u003e *P\u003c/em\u003e\u0026lt;0.05 vs SCI(moderate), at 2 days,\u003csup\u003e\u003cem\u003e #\u003c/em\u003e\u003c/sup\u003e\u003cem\u003eP\u003c/em\u003e\u0026lt;0.05 vs SCI(moderate), at 5 days. (E) and (F):Validation of FOS mRNA and CCL2 mRNA expression in SCI rats via qRT-PCR, \u003cem\u003e*P\u003c/em\u003e\u0026lt;0.05 vs Sham group. \u003csup\u003e\u003cem\u003e#\u003c/em\u003e\u003c/sup\u003e\u003cem\u003eP\u003c/em\u003e\u0026lt;0.05 vs SCI group.\u003c/p\u003e","description":"","filename":"Figure7.png","url":"https://assets-eu.researchsquare.com/files/rs-5504847/v1/c44f0d51a0b31f5dbbf4c969.png"},{"id":86812994,"identity":"fdad09d7-c0d7-488e-a6f4-611ab07b56e6","added_by":"auto","created_at":"2025-07-15 22:01:35","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":20743003,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-5504847/v1/d70abbc4-bf18-4056-861e-a75d8360f2b9.pdf"}],"financialInterests":"No competing interests reported.","formattedTitle":"Identification of key genes associated with muscle atrophy after spinal cord injury and experimental verification in rats","fulltext":[{"header":"Introduction","content":"\u003cp\u003eSpinal cord injury (SCI) refers to a neurological disorder characterized by damage to the spinal cord, resulting from direct or indirect trauma, leading to dysfunction of motor, sensory, and autonomic functions below the level of injury within the neuroanatomical distribution\u003csup\u003e \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e \u003c/sup\u003e. According to statistics, approximately 100,000 to 800,000 individuals in China experience SCI annually, while in the United States, each SCI patient incurs lifetime costs ranging from \u003cspan\u003e$\u003c/span\u003e1.1\u0026nbsp;million to \u003cspan\u003e$\u003c/span\u003e4.7\u0026nbsp;million, imposing significant burdens on both families and society, severely affecting the physiological, psychological health, and quality of life of patients\u003csup\u003e \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e,\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e \u003c/sup\u003e. With the rapid development of modern transportation, construction, and industry, the incidence and disability rate of SCI have been steadily increasing each year\u003csup\u003e \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e \u003c/sup\u003e.Mechanical compression following SCI can cause direct mechanical damage to spinal cord tissue, including local hemorrhage, edema, and ischemic reactions, a process referred to as primary injury\u003csup\u003e \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e \u003c/sup\u003e.This type of injury triggers a cascade of cellular and molecular responses, including inflammation, lipid peroxidation, free radical generation, disruption of ion channels, axonal demyelination, glial scar formation, necrosis, and programmed cell death, collectively termed secondary injury\u003csup\u003e \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e \u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eOne of the most apparent consequences following SCI is a significant decrease in muscle mass due to partial or complete interruption of descending motor pathways\u003csup\u003e\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e\u003c/sup\u003e. Post-SCI muscle atrophy is a cascade response to SCI, and failure to actively prevent its occurrence severely affects patients' mobility and quality of life\u003csup\u003e\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u003c/sup\u003e. Following SCI, changes occur in spinal cord motor neurons and nerve fibers, leading to dysfunction of motor units, which results in decreased muscle fatigue resistance, reduced strength, and slowed motor conduction velocity. These alterations exacerbate functional decline post-SCI and impede motor recovery\u003csup\u003e\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e\u003c/sup\u003e.Given the potential modifiability of post-SCI muscle atrophy\u003csup\u003e\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e\u003c/sup\u003e, studying its mechanisms and providing effective treatment strategies have become recent focal points. With the development of bioinformatics and the widespread use of high-throughput gene chip technology in exploring molecular changes in disease progression, there is now a theoretical basis and research direction for understanding the pathogenesis of post-SCI muscle atrophy\u003csup\u003e\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e\u003c/sup\u003e. This facilitates the targeted selection of genes with research value, potentially saving significant research time and effort.\u003c/p\u003e \u003cp\u003eThis study used the post SCI muscle atrophy expression profile data (GSE21497) from the gene expression omnibus database (GEO) as the research object, screened differentially expressed genes (DEGs), and performed weighted gene co-expression network analysis (WGCNA) and functional enrichment analysis to construct a protein protein interaction (PPI) network and screen for key genes. These analyses aim to provide a theoretical basis for exploring the molecular mechanisms underlying the development and maintenance of muscle atrophy post-SCI, as well as identifying important therapeutic targets for treating muscle atrophy. Finally, the results were verified by real-time polymerase chain reaction (PCR).\u003c/p\u003e"},{"header":"Materials and Methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003eMicroarray data analysis\u003c/h2\u003e \u003cp\u003eMicroarray and gene expression profile dataset GSE21497 were obtained from the freely accessible public database NCBI-GEO (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.ncbi.nlm.nih.gov/geo/\u003c/span\u003e\u003cspan address=\"https://www.ncbi.nlm.nih.gov/geo/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). The gene sequencing platform used was GPL570 [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array. This database employed a self-controlled method to collect samples on the 2nd and 5th days after SCI occurrence from a total of 10 patients (9 males, 1 female, with an average age of 44 years, including 6 cases of tetraplegia and 4 cases of paraplegia). Muscle biopsy samples from the lateral vastus muscle were obtained. The online tool GEO2R (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.ncbi.nlm.nih.gov/geo/geo2r/\u003c/span\u003e\u003cspan address=\"https://www.ncbi.nlm.nih.gov/geo/geo2r/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) was utilized to perform DEGs analysis between the two days and five days post-SCI group.\u003c/p\u003e \u003c/div\u003e\n\u003ch3\u003eData processing and selection of DEGs\u003c/h3\u003e\n\u003cp\u003eUsing the online Sangerbox 3.0 (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://vip.sangerbox.com/home.html\u003c/span\u003e\u003cspan address=\"http://vip.sangerbox.com/home.html\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e), we performed ID conversion and gene annotation processing on dataset GSE21497\u003csup\u003e12\u003c/sup\u003e. Differential analysis was conducted using the R software package limma (version 3.40.6) to obtain DEGs between the two groups\u003csup\u003e\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e\u003c/sup\u003e. Based on the acquired expression profile dataset, the data were log2-transformed, and multivariate linear regression was performed using the lmFit function. Subsequently, the eBays function was employed to obtain the expression matrix, resulting in the identification of significant differences for each gene. DEGs were selected using a threshold of \u003cem\u003eP\u0026thinsp;\u0026lt;\u003c/em\u003e\u0026thinsp;0.05 and |logFC|\u003cem\u003e\u0026gt;\u003c/em\u003e1.5. The online Sangerbox 3.0 software was then utilized to input the original data in TXT format, facilitating online construction of DEGs and the generation of volcano plots and heatmaps.\u003c/p\u003e\n\u003ch3\u003eWeighted Gene Co-Expression Network Analysis(WGCNA)\u003c/h3\u003e\n\u003cp\u003ePerforming weighted gene co-expression network analysis (WGCNA) using the R software package on the obtained DEGs\u003csup\u003e\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e\u003c/sup\u003e.Inputting the DEGs data for sample selection, wherein sample and grouping information were used to set a cut-Height threshold for sample exclusion, removing outliers, and constructing sample clustering. After preparing the data for sample clustering observations, determining the optimal soft threshold (β) to ensure the constructed network conforms to a scale-free topology. Constructing a co-expression network based on the optimal soft threshold and partitioning genes into different modules, represented by branch colors on the clustering tree, and plotting gene clustering trees and module eigengene clustering diagrams. Based on these results, calculating the correlation between gene modules and phenotypes, identifying trait-associated modules, and plotting phenotype-gene correlation heatmaps.\u003c/p\u003e\n\u003ch3\u003eFunctional enrichment analysis\u003c/h3\u003e\n\u003cp\u003eGene Ontology (GO) and Kyoto Encyclopedia of Genes Genomes (KEGG) pathway enrichment analyses were conducted on DEGs utilizing the clusterprofiler in R package (version 3.14.3, \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://bioconductor.org)\u003c/span\u003e\u003cspan address=\"http://bioconductor.org)\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003csup\u003e15,16\u003c/sup\u003e\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e.Usin\u003c/span\u003e\u003cspan address=\"http://.Usin\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003eg the Database for Annotation, Visualization, and Integrated Discovery (DAVID, \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://david.ncifcrf.gov/home.jsp\u003c/span\u003e\u003cspan address=\"https://david.ncifcrf.gov/home.jsp\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) for online gene function annotation and signal pathway enrichment analysis\u003csup\u003e\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e\u003c/sup\u003e, including biological processes (BP), molecular functions (MF), and cellular components (CC), as well as KEGG pathway analysis. Genes count\u0026thinsp;\u003cem\u003e\u0026ge;\u003c/em\u003e\u0026thinsp;2 and \u003cem\u003eP\u0026thinsp;\u0026lt;\u003c/em\u003e\u0026thinsp;0.05 was used as the threshold of significant enrichment results.\u003c/p\u003e\n\u003ch3\u003eGene Set enrichment analysis(GSEA)\u003c/h3\u003e\n\u003cp\u003eTo observe whether there are changes in the distribution of pathway genes in all gene expression data between different time point groups of SCI patients. We obtained GSEA version3.0 from the GSEA website (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://software.broadinstitute.org/gsea/\u003c/span\u003e\u003cspan address=\"http://software.broadinstitute.org/gsea/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e index.jsp) and divided the samples into two groups based on different observation time points\u003csup\u003e\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e\u003c/sup\u003e. Subsequently, we downloaded the subset from the Molecular Signatures Database (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://www.gsea-msigdb.org/gsea/downloads.jsp\u003c/span\u003e\u003cspan address=\"http://www.gsea-msigdb.org/gsea/downloads.jsp\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) to assess relevant pathways and molecular mechanisms\u003csup\u003e\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e\u003c/sup\u003e. Based on gene expression profiles and phenotype grouping, a minimum gene set of 5 and a maximum gene set of 5000 were established, with 1000 permutations. A significance threshold of \u003cem\u003eP\u0026thinsp;\u0026lt;\u003c/em\u003e\u0026thinsp;0.05 and \u003cem\u003eFDR\u0026thinsp;\u0026lt;\u003c/em\u003e\u0026thinsp;0.25 was considered statistically significant.\u003c/p\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003eConstruction of protein-protein interaction (PPI) networks and identification of hub genes\u003c/h2\u003e \u003cp\u003eUsing the STRING online database (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://string-db.org/\u003c/span\u003e\u003cspan address=\"https://string-db.org/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e), the selected DEGs were inputted to analyze and construct PPI networks. The obtained data were imported into Cytoscape software to build protein co-expression network diagrams for comprehensive bioinformatics analysis\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e. Using the CytoHubba plugin in Cytoscape, the top 30 genes were selected, and subsequently, the top 10 hub genes were determined using various algorithms including Degree, maximum clique centrality (MCC), maximum neighborhood component (MNC), and edge percolated component (EPC). CytoHubba comprehensively uses four algorithms, namely Degree, MCC, MNC, and EPC, to evaluate the global connectivity, local compactness, topological importance of nodes, and edge permeability respectively. Selecting multiple algorithms can avoid the bias of a single algorithm, and cross-validation with multiple algorithms can improve the robustness of the results. To ensure the reliability of the screening results, we adopt a two-step weighting strategy: algorithm weight assignment and comprehensive score calculation.The raw data of the identified hub genes were inputted into the online Sangerbox 3.0 software to find intersections between sets and to generate Venn diagrams. Organize and analyze the differential expression of hub genes at different stages. Conduct statistical analysis on the expression differences of hub genes in the presence of injuries of varying degrees at the same stage.\u003c/p\u003e \u003c/div\u003e\n\u003ch3\u003eEstablishment of a spinal cord injury model\u003c/h3\u003e\n\u003cp\u003e Forty healthy, adult, female Sprague-Dawley (SD) rats of specific pathogen-free (SPF) grade were obtained from the Animal Experiment Center of Shanxi Medical University (Animal Use License Number: SYXK (Jin) 2019 0008), with body weights ranging from 210 to 250 g. The experimental procedures complied with ethical standards for animal care and use. This study was approved by the Animal Ethics Committee of Shanxi Bethune Hospital (Approval Number: SBQKL-2021-029). All methods were performed in accordance with ARRIVE guidelines and relevant and regulations. The rats were randomly assigned into two groups (n\u0026thinsp;=\u0026thinsp;20) using a random number table: the sham operation group (Sham group) and the spinal cord injury group (SCI group). In the Sham group, the rats underwent laminectomy and removal of the spinous processes at the T\u003csub\u003e8\u0026thinsp;\u0026minus;\u0026thinsp;9\u003c/sub\u003e vertebrae without spinal cord injury. In the SCI group, spinal cord injury was induced using the Allen method.\u003c/p\u003e\n\u003ch3\u003eValidation of hub genes via real-time PCR\u003c/h3\u003e\n\u003cp\u003eAfter the model was successfully established, 10 rats from each group were sacrificed at 2 and 5 days post-operation. All rats were euthanized by cervical dislocation under deep anesthesia with 2% sodium pentobarbital (100mg/kg), followed by transcardial perfusion normal saline (4\u0026deg;C). The lateral thigh muscles of rats was immediately harvested and stored at -80\u0026deg;C.Total RNA was extracted using TRIzol reagent (Invitrogen), followed by reverse transcription to synthesize cDNA. PCR amplification was performed, and electrophoresis was used for product detection and data analysis. The PCR conditions were as follows: pre-denaturation at 95\u0026deg;C for 30 s, denaturation at 94\u0026deg;C for 1 min, annealing at 56\u0026deg;C for 1 min, and extension at 72\u0026deg;C for 90 s, with a total of 35 cycles, followed by a final extension at 72\u0026deg;C for 10 min. In this study, the internationally commonly used 2\u003csup\u003e\u0026minus;ΔΔCT\u003c/sup\u003e method was employed to calculate the relative mRNA expression levels.\u003c/p\u003e \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003eStatistical Analysis\u003c/h2\u003e \u003cp\u003eData analysis was conducted using SPSS 26.0 statistical software(IBM Corporation, Armonk, NY, USA). The Shapiro-Wilk test was used to assess normality. For data that followed a normal distribution, measurements are presented as mean\u0026thinsp;\u0026plusmn;\u0026thinsp;SD. To control the family-wise error rate (FWER) arising from multiple comparisons between time points, we applied the Benjamini-Hochberg false discovery rate (FDR) correction. Adjusted p-values were calculated using the GraphPad Prism (version 9.5.0), with a significance threshold set at FDR\u0026thinsp;\u0026lt;\u0026thinsp;0.05. Group differences were analyzed by one-way ANOVA followed by Tukey\u0026rsquo;s post hoc test for pairwise comparisons. Figures were created using Adobe Photoshop CS6 and GraphPad Prism 9.5.0. Statistical significance was set at P\u0026thinsp;\u0026lt;\u0026thinsp;0.05.\u003c/p\u003e \u003c/div\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003ch2\u003eNormalization of data and analysis of differential gene expression\u003c/h2\u003e \u003cp\u003eGSE21497 was downloaded from the GEO database and analyzed using the online tool GEO2R. In the box plot, the black lines are nearly at the same level, indicating minimal deviation between the two (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eA). The analysis using the Uniform Manifold Approximation and Projection (UMAP) algorithm reveals distinct clustering between the SCI day 2 group and the SCI day 5 group, indicating significant differences between the two groups(Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eB). After filtering with \u003cem\u003eP\u0026thinsp;\u0026lt;\u003c/em\u003e\u0026thinsp;0.05 and |logFC|\u003cem\u003e\u0026gt;\u003c/em\u003e1.5, a total of 1007 differentially expressed genes (DEGs) were identified on the 5th day post-SCI, including 533 upregulated genes and 474 downregulated genes. Volcano plots and heatmaps were generated using the online Sangerbox 3.0 tool(Figs.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eC and \u003cb\u003eD\u003c/b\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec14\" class=\"Section2\"\u003e \u003ch2\u003eData preparation and sample clustering\u003c/h2\u003e \u003cp\u003eThe data matrix consists of a total of 20 samples, with 1007 DEGs. Further processing was performed using the WGCNA package in R software to conduct sample clustering. The clustering tree was examined to identify and filter out any outliers or anomalies. Based on the visual results of the sample dendrogram, the hierarchical clustering of the samples was obtained(Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eE).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec15\" class=\"Section2\"\u003e \u003ch2\u003eSelecting the optimal soft threshold\u003c/h2\u003e \u003cp\u003eAfter completing the sample clustering observations, the optimal soft threshold (soft thresholding or power) was determined to ensure the constructed network conforms to a scale-free topology. A scatter plot was generated for power values, and the results (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eA) showed a clear trend throughout the range from 1, but the trend became less apparent after 18. Therefore, 18 was selected as the optimal power value (with a Power R-squared greater than 0.8). The mean connectivity plot (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eA) depicts the network connectivity under different soft threshold values, and the relationship between average connectivity and power values was examined. The connectivity stabilizes at a value of 18, indicating its suitability for power selection.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec16\" class=\"Section2\"\u003e \u003ch2\u003eConstructing co-expression network\u003c/h2\u003e \u003cp\u003eBased on the optimal soft threshold, a co-expression network was constructed, and genes were divided into different modules. A gene clustering tree (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eB) can be plotted, where the upper part represents the hierarchical clustering dendrogram of genes, and the lower part represents gene modules. Correspondence between the upper and lower parts indicates that genes with closer distances (clustered into the same branch) are assigned to the same module. Each color represents a module, with the grey module considered as a collection of genes that cannot be assigned to any module.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec17\" class=\"Section2\"\u003e \u003ch2\u003eModule correlation clustering and analysis of phenotype-gene correlation\u003c/h2\u003e \u003cp\u003eExplore the correlation between various modules, calculate module eigengenes, and plot a network diagram of module eigengenes. Darker colors indicate higher correlation coefficients (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eC). Based on the above results, a heatmap of phenotype-gene correlation was plotted(Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eD). DEGs were categorized into four modules: blue, brown, turquoise, and grey. The color blocks on the right represent modules, while the color bar in the top right corner represents the range of correlations. Darker colors in the heatmap indicate higher correlation, with red indicating positive correlation and green indicating negative correlation. The numbers in each cell represent correlation and significance. Based on the absolute value of correlation, the most relevant modules were selected. The results indicate that the turquoise module (0.80, \u003cem\u003ep\u0026thinsp;\u0026lt;\u003c/em\u003e\u0026thinsp;0.0001) is the most negatively and positively correlated module with SCI. Therefore, the turquoise module, which includes 161 DEGs, is considered the key module associated with SCI. The grey module contains all genes that were not involved in clustering and is therefore considered an invalid module, not suitable for further analysis.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec18\" class=\"Section2\"\u003e \u003ch2\u003eFunctional enrichment analysis\u003c/h2\u003e \u003cp\u003eThe 161 DEGs in the turquoise module were subjected to GO and KEGG pathway analysis using the DAVID tool. The GO analysis revealed enrichment of the DEGs in 147 biological processes (BP), 146 cellular components (CC), and 145 molecular functions (MF) within the turquoise module. In the BP analysis, the DEGs are significantly enriched in terms such as negative regulation of cellular process, cytosol and response to organic substance(Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eA). In the CC analysis, the DEGs are significantly enriched in terms such as cytosol, endomembrane system and extracellular region(Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eB). In the MF analysis, the DEGs are significantly enriched in terms such as extracellular matrix structural, organic anion transmembrane transporter activity constituent and extracellular matrix binding(Figs.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eC and \u003cb\u003eD\u003c/b\u003e). KEGG analysis indicates the involvement of the turquoise module in various signaling pathways including the TNF signaling pathway, cellular senescence and regulation of actin cytoskeleton(Figs.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eE and \u003cb\u003eF\u003c/b\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec19\" class=\"Section2\"\u003e \u003ch2\u003eGene set enrichment analysis (GSEA)\u003c/h2\u003e \u003cp\u003eThe Gene Set Enrichment Analysis (GSEA) examined whether there were changes in the distribution of pathway gene sets among different groups of SCI patients. The larger the |NES| value, the smaller the FDR value, indicating higher reliability of the analysis results. The results indicate significant enrichment and upregulation of genes associated with the peroxisome proliferators-activated receptors (PPARs) signaling pathway (NES\u0026thinsp;=\u0026thinsp;1.6916, FDR\u0026thinsp;=\u0026thinsp;0.2060) and the adherens junction signaling pathway (NES\u0026thinsp;=\u0026thinsp;1.6957, FDR\u0026thinsp;=\u0026thinsp;0.2485) in both groups(Figs.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eA and \u003cb\u003eB\u003c/b\u003e). Conversely, genes related to the DNA replication signaling pathway (NES=-1.6676, FDR\u0026thinsp;=\u0026thinsp;0.2274) show significant enrichment and downregulation in both groups(Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eC).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec20\" class=\"Section2\"\u003e \u003ch2\u003ePPI network construction and selection and analysis on hub genes\u003c/h2\u003e \u003cp\u003eUsing the STRING database (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://cn.string-db.org\u003c/span\u003e\u003cspan address=\"https://cn.string-db.org\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e), a PPI network was constructed for the enriched DEGs in the turquoise module, with a threshold set to 0.4. The results were exported for further analysis. The file containing node and edge information was imported into Cytoscape using the appropriate plugin, resulting in the PPI network diagram (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eA). Subsequently, the top 30 hub genes were selected (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eB). Utilizing the CytoHubba plugin in Cytoscape, hub genes were identified using four different calculation methods: Edge Percolated Component (EPC), Maximum Clique Centrality (MCC), Maximum Neighborhood Component (MNC), and Degree. This process aided in identifying the top 10 critical genes (Figs.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eA-D). The top 10 gene sets obtained from the four different algorithms were input into Sangerbox 3.0 software to determine the number of intersections in the datasets (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eA). Finally, FOS and CCL2 were identified as the hub genes for our screening.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eBased on the data of 10 samples from the gene expression profile GSE21497 database in this study, it was found that there were differences in the expressions of FOS and CCL2 at 2 days and 5 days after the injury. Specifically, compared with the situation at 2 days after the injury, the protein expressions of FOS and CCL2 both decreased at 5 days after the injury (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eB).We analyzed the expression levels of FOS and CCL2 in the lower limb muscle tissues of 6 patients with tetraplegia (severe SCI) and 4 patients with paraplegia (moderate SCI). The results showed that at 2 days, the expression levels of FOS and CCL2 in the lower limb muscle tissues of the 6 tetraplegic patients were higher. By the 5th day, although the expression levels of FOS and CCL2 in their lower limb muscle tissues decreased to some extent, they still remained at a relatively high level (Figs.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eC and \u003cb\u003eD\u003c/b\u003e). This indicates that the higher the expression levels of FOS and CCL2, the more severe the muscle atrophy in patients.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec21\" class=\"Section2\"\u003e \u003ch2\u003eValidation of hub genes\u003c/h2\u003e \u003cp\u003eTo further validate the reliability of WGCNA results, PCR was used to reconfirm two Hub genes.The primers for gene expression were designed as follows:\u003c/p\u003e \u003cp\u003eFOS: forward 5'-AACAGATCCGAGCAGCTTCTA-3', reverse5'-TTTTGAGCTTCAACCGGCATC-3'\u003c/p\u003e \u003cp\u003eCCL2: forward 5'-TCTTGCTGCACGTCCTTTATT-3',\u003c/p\u003e \u003cp\u003ereverse 5'-GCATCCAGTTTTTGTATGTGCC-3'.\u003c/p\u003e \u003cp\u003eThe expression of FOS and CCL2 genes in the SCI model was analyzed via PCR. The results of our PCR experiments suggest that the expressions of FOS mRNA and CCL2 mRNA peaked at 2 days after SCI, and then slowly decreased at 5 days after SCI, but still remained at a relatively high level. (Figs.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eE and \u003cb\u003eF\u003c/b\u003e) (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026le;\u0026thinsp;0.05).\u003c/p\u003e \u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eThe secondary neuropathic changes following SCI are the primary cause of muscle atrophy, which not only impacts the care management and daily activities of SCI patients but also increases the risk of secondary complications (such as osteoporosis, diabetes, and cardiovascular diseases) \u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e,\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e\u003c/sup\u003e.Therefore, preventing muscle atrophy is crucial for maintaining metabolic health and normal daily activities after SCI. The pathogenesis of muscle atrophy following SCI involves various factors, including signal transduction, immune response, electrical conduction, stimulation, and metabolism\u003csup\u003e\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e\u003c/sup\u003e. SCI causes loss of central regulation function in peripheral nerves below the damaged segment, resulting in sensory, motor, and autonomic dysfunction, muscle paralysis, and decreased muscle loading\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e. Impairment of the neuromuscular system and reduced integrity of the musculoskeletal system are critical features of SCI and significant barriers to motor function recovery\u003csup\u003e\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e\u003c/sup\u003e. However, the etiology and precise underlying mechanisms of muscle atrophy after SCI remain poorly understood. The aim of our study is to identify relevant genes, biological processes, and pathway pathways that may play crucial roles in the development of muscle atrophy after SCI, thereby providing theoretical insights into the pathogenesis and therapeutic targets of muscle atrophy following SCI.\u003c/p\u003e \u003cp\u003eIn this study, we downloaded the gene dataset GSE21497. After data normalization and standardization, we used the online Sangerbox 3.0 software for computation, resulting in 1007 differentially expressed genes (DEGs) on day 5 after SCI, including 533 upregulated genes and 474 downregulated genes. WGCNA analysis was performed on the data matrix containing 20 samples and 1007 genes with expression differences. WGCNA can identify sets of genes with similar expression patterns, elucidate the relationships between gene sets and sample phenotypes, construct regulatory networks among genes within gene sets, and identify key regulatory genes\u003csup\u003e\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e\u003c/sup\u003e. All 1007 DEGs obtained through the \u0026ldquo;limma\u0026rdquo; algorithm were clustered into four modules. It was found that the 161 genes in the turquoise module were most closely related to SCI.\u003c/p\u003e \u003cp\u003eGO enrichment analysis involves examining whether genes in a gene set are enriched at different levels to understand their functions in biological processes(BP), cellular components(CC), and molecular functions(MF) \u003csup\u003e\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e\u003c/sup\u003e. The results of this study show that DEGs are significantly enriched in biological processes such as \u0026ldquo;negative regulation of cellular process, cytosol and response to organic substance\u0026rdquo;. DEGs are also significantly enriched in cellular components including \u0026ldquo;cytosol, endomembrane system and extracellular region\u0026rdquo;, as well as molecular functions such as \u0026ldquo;extracellular matrix structural, organic anion transmembrane transporter activity constituent and extracellular matrix binding\u0026rdquo;. KEGG enrichment analysis is used to analyze the enrichment of gene pathways and functions in a given gene set\u003csup\u003e\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e\u003c/sup\u003e. KEGG analysis reveals that the turquoise module is involved in various signaling pathways such as the \u0026ldquo;TNF signaling pathway, ellular senescence and regulation of actin cytoskeleton\u0026rdquo;.\u003c/p\u003e \u003cp\u003eGSEA is used to evaluate the distribution trend of genes from a predefined gene set within a ranked gene expression dataset associated with phenotypes, thereby determining their contribution to the phenotype\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. We utilized the GSEA method to analyze the results, which indicated that following SCI, genes associated with the peroxisome proliferators-activated receptors (PPARs) signaling pathway and adherens junction signaling pathway were significantly increased in both groups, with upregulated expression, while genes related to the DNA replication signaling pathway were significantly decreased in both groups, with downregulated expression. PPARs pharmacological activation can be considered as a multifaceted therapeutic target, exerting anti-inflammatory, antioxidant, antiexcitotoxic, and energy-generating effects in some neurologic and inflammation-related diseases\u003csup\u003e\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e\u003c/sup\u003e. The Adherens junction signaling pathway also demonstrates significant enrichment in comprehensive proteomic and phosphoproteomic analyses following SCI\u003csup\u003e\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e\u003c/sup\u003e, reducing the loss of tight junction and adherens junction proteins in human brain microvascular endothelial cells (HBMECs). This reduction can attenuate blood-spinal cord barrier (BSCB) permeability and degradation of tight junction molecules (such as P120, β-catenin, Occludin, and Claudin5), ultimately improving functional recovery in a traumatic rat model\u003csup\u003e\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e\u003c/sup\u003e. Studies have shown that enhancement of the peroxisome proliferator-activated receptor-γ coactivator-1α(PGC-1α) mitochondrial biogenesis pathway upregulates nuclear respiratory factors (NRF-1,2) and mitochondrial transcription factor A (Tfam), stimulating mitochondrial biogenesis\u003csup\u003e\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e\u003c/sup\u003e.This leads to increased mitochondrial DNA replication and gene transcription, reducing oxidative damage and providing protective effects against atrophy-related skeletal muscle degeneration.\u003c/p\u003e \u003cp\u003eThe PPI network illustrates the relationships among proteins and even genes within a given gene set that play a protective role against atrophy-related skeletal muscle degeneration\u003csup\u003e\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e\u003c/sup\u003e. It is one of the commonly used methods for identifying core genes. Cytoscape (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://cytoscape.org/\u003c/span\u003e\u003cspan address=\"https://cytoscape.org/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) is an excellent software for graphical visualization and analysis of networks and can be equipped with various functional plugins, making it widely used in PPI analysis\u003csup\u003e\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e\u003c/sup\u003e. Based on the PPI network and the CytoHubba plugin with four different algorithms, the top 10 key genes were identified. These top 10 genes from each algorithm were input into the online Sangerbox 3.0 software for analysis, resulting in a Venn diagram showing FOS and CCL2 as the hub genes we selected. Moreover, the analysis of the data from the gene expression profile GSE21497 shows that the higher the expression levels of FOS and CCL2 are, the more severe the muscle atrophy of the patients will be. It was found that there were differences in the expressions of FOS and CCL2 at 2 days and 5 days after the injury. The results of our PCR experiments suggest that the expressions of FOS mRNA and CCL2 mRNA peaked at 2 days after SCI, and then slowly decreased at 5 days after SCI, but still remained at a relatively high level. Therefore, this dynamic change indicates that inhibiting the expression of CCL2 in the early stage of SCI and upregulating FOS in the chronic stage (such as local injection of FOS overexpression vectors) may be effective methods to alleviate muscle atrophy after SCI.\u003c/p\u003e \u003cp\u003eFOS belongs to a class of genes known as immediate early genes (IEGs). It encodes a nuclear phosphoprotein that serves as a transcription factor\u003csup\u003e\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e\u003c/sup\u003e. FOS plays a significant role in regulating cell growth, division, proliferation, differentiation, and even programmed cell death\u003csup\u003e\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e\u003c/sup\u003e. As the protein product of the \u003cem\u003ec-fos\u003c/em\u003e gene, under normal conditions, \u003cem\u003ec-fos\u003c/em\u003e expression is minimal within cells. But, to some extent, the number of FOS-immunoreactive neurons correlates positively with the intensity of the stimulus received, making FOS a marker of injurious stimuli\u003csup\u003e\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e\u003c/sup\u003e. Research has shown that SCI significantly affects \u003cem\u003ec-fos\u003c/em\u003e expression in the paraventricular nucleus (PVN) and supraoptic nucleus (SON) of the hypothalamus in rats. Following SCI at 1, 6, 12, and 24 hours, there is a significant increase in \u003cem\u003ec-fos\u003c/em\u003e expression in the PVN, which persists for more than 24 hours. At 12 hours post-SCI, there is a significant increase in \u003cem\u003ec-fos\u003c/em\u003e expression in the SON, reaching its peak at 24 hours post-SCI\u003csup\u003e\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eChemokines are small-molecule peptides secreted by immune or non-immune cells, induced by inflammatory cytokines, growth factors, and pathological stimuli. They play an important role in inflammatory diseases, wound healing, spinal cord injury, and other processes\u003csup\u003e\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e\u003c/sup\u003e.Chemokine C-C motif ligand 2 (CCL2), also known as monocyte chemoattractant protein-1 (MCP-1), is a detrimental factor released after neural injury, primarily involved in the local inflammatory response following SCI\u003csup\u003e\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e\u003c/sup\u003e. After SCI, chronic pain, circulatory disturbances, and mechanical compression lead to changes in neurons of the dorsal root ganglion (DRG). After peripheral nerve or ganglion injury, CCL2 expression is upregulated in DRG neurons\u003csup\u003e\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e\u003c/sup\u003e. Research has shown that CCL2 plays a crucial role in recruiting macrophages during the early stages of muscle injury. Mice lacking CCL2 exhibit incomplete recovery of gastrocnemius muscle mass, muscle fiber cross-sectional area, and gastrocnemius muscle contraction characteristics during disuse atrophy recovery. Additionally, mice lacking CCL2 demonstrate reduced turnover of skeletal muscle collagen and significantly decreased recruitment of macrophages to the gastrocnemius muscle, resulting in poor recovery of muscle size and function as well as abnormal collagen remodeling\u003csup\u003e\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eWe verified the differential expression of FOS and CCL2 in the SCI animal model through PCR, confirming the reliability of the bioinformatics analysis. The potential role of FOS: As an activator protein-1(AP-1) transcription factor, it may promote the activity of the ubiquitin-proteasome system by activating the Atrogin-1/MuRF1 pathway, thus facilitating the degradation of muscle fibers\u003csup\u003e\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e\u003c/sup\u003e. The potential role of CCL2: Recruiting macrophages through the CCR2 receptor, inducing the release of pro-inflammatory factors (such as TNF-α and IL-6), and accelerating the apoptosis of muscle cells\u003csup\u003e\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e\u003c/sup\u003e. These analyses link the changes in gene expression with known pathological mechanisms, pointing the way for subsequent functional studies. In terms of drug research and development, the already marketed JNK inhibitor (such as CC-930) can indirectly inhibit the activity of FOS by blocking the formation of the AP-1 complex\u003csup\u003e\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e\u003c/sup\u003e. The CCL2 antagonist (Bindarit, a CCL2 synthesis inhibitor) can reduce the infiltration density of macrophages and improve neuromotor function\u003csup\u003e\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e\u003c/sup\u003e.In terms of clinical translation, the serum CCL2 level can be used as a biomarker to predict the effective rate of Bindarit treatment and guide patient stratification. Clinical trials of the synergistic use of FOS inhibitors and rehabilitation training can also be carried out.\u003c/p\u003e \u003cp\u003eBased on the existing research, the pro-inflammatory and pro-atrophic mechanisms of FOS and CCL2 make them ideal targets for the treatment of muscle atrophy following SCI. For example, the small molecule inhibitors T-5224 (targeting FOS/AP-1) and RS504393 (targeting CCL2/CCR2) have demonstrated anti-inflammatory and anti-fibrotic effects in other disease models\u003csup\u003e\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e,\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e\u003c/sup\u003e. Future studies could verify whether these drugs alleviate muscle atrophy in the SCI model by inhibiting the MuRF1/MAFbx pathway or regulating macrophage polarization. Additionally, combining gene silencing techniques or existing neurotrophic therapies may produce synergistic effects. In conclusion, as hub genes, FOS and CCL2, through the closed-loop path of \"network screening-mechanism analysis-drug verification-clinical transformation\", provide powerful targets for the development of precise drugs against muscle atrophy.\u003c/p\u003e \u003cdiv id=\"Sec23\" class=\"Section2\"\u003e \u003ch2\u003eLimitations\u003c/h2\u003e \u003cp\u003eOur study has several limitations. Firstly, the GSE21497 dataset used in this study is from the GEO database. Although it provides valuable information on the gene expression of muscle atrophy after SCI, there are still some limitations such as sample heterogeneity, and the lack of clinical information. Secondly,the samples included in this study are concentrated within 5 days after injury, which may not fully reflect the molecular characteristics of muscle atrophy in the chronic phase. Therefore, there may be a bias in temporal dynamics during the sample selection process. Finally, we have not directly verified the quantitative relationship between the expression of FOS and CCL2 and the severity of muscle atrophy in patients. In the later stage, we will continue to conduct clinical research related to this project.\u003c/p\u003e \u003c/div\u003e"},{"header":"Conclusions","content":"\u003cp\u003eIn summary, we utilized WGCNA to analyze the expression profile data from GSE21497, identifying two hub genes, FOS and CCL2, in muscle atrophy following SCI. FOS and CCL2 may be involved in the molecular pathophysiology of muscle atrophy after SCI, serving as potential targets for diagnosis or treatment of SCI-related muscle atrophy.\u003c/p\u003e"},{"header":"Abbreviations","content":"\u003cp\u003eAP-1,activator protein-1\u003c/p\u003e\n\u003cp\u003eBP,biological processes\u003c/p\u003e\n\u003cp\u003eBSCB,blood-spinal cord barrier\u003c/p\u003e\n\u003cp\u003eCC,cellular components\u003c/p\u003e\n\u003cp\u003eCCL2,Chemokine C-C motif ligand 2\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eDAVID, database for annotation, visualization, and integrated discovery\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eDRG,dorsal root ganglion\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eDEGs, differentially expressed genes\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eEPC, edge percolated component\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eGEO, gene expression omnibus\u003c/p\u003e\n\u003cp\u003eGO, gene ontology\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eGSEA,gene set enrichment analysis\u003c/p\u003e\n\u003cp\u003eHBMECs,human brain microvascular endothelial cells\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eKEGG, kyoto encyclopedia of genes genomes\u003c/p\u003e\n\u003cp\u003eMCC, maximum clique centrality\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eMCP-1,monocyte chemoattractant protein-1\u003c/p\u003e\n\u003cp\u003eMF,molecular functions\u003c/p\u003e\n\u003cp\u003eMNC,maximum neighborhood component\u0026nbsp;\u003c/p\u003e\n\u003cp\u003ePCR, polymerase chain reaction\u0026nbsp;\u003c/p\u003e\n\u003cp\u003ePPARs,peroxisome proliferators-activated receptors\u003c/p\u003e\n\u003cp\u003ePPI, protein protein interaction\u0026nbsp;\u003c/p\u003e\n\u003cp\u003ePVN,paraventricular nucleus\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eSCI, spinal cord injury\u003c/p\u003e\n\u003cp\u003eSON,supraoptic nucleus\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eWGCNA, weighted gene co-expression network analysis\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eAuthor Contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eX.W. conceived and designed the study; Y.G. analyzed and interpreted the data and drafted the work. X.W. assisted in revising the manuscript. J.H. reviewed the manuscript and supervised the whole project. All authors have approved the submission of the final version of the manuscript and agree to take responsibility for all aspects of the work.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis study was funded by the Taiyuan City Science and Technology Plan Project (Program NO.202201), the Key Medical Research Project of Shanxi Province (Program NO.2023XM034), and the Basic Research Plan of Shanxi Province(Program NO. 202103021223410,202303021222381 and 202303021221246 ).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eData availability\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe data used in this study are available in public databases. Some data and all scripts for conducting the analyses in this study can be shared upon reasonable request to the\u003c/p\u003e\n\u003cp\u003ecorresponding author.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDisclosure\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors have no conflicts of interest to declare.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eEthics statement\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis study was approved by the Animal Ethics Committee of Shanxi Bethune Hospital (Approval Number: SBQKL-2021-029).The study is also reported in accordance with ARRIVE guidelines.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eLiu, Z., Chen, X., Chen, P. \u0026amp; Wang, L. Transcranial direct current stimulation attenuates the chronic pain of osteoarthritis in rats via reducing NMDAR2B expressions in the spinal cord. \u003cem\u003eJ Musculoskelet Neuronal Interact\u003c/em\u003e. \u003cstrong\u003e23\u003c/strong\u003e(4),498-505(2023).\u003c/li\u003e\n\u003cli\u003eFallah, N., Noonan,V.K., \u0026amp; Sharwood,L.N. Editorial: Epidemiology, evidence-based care, and outcomes in spinal cord injury. \u003cem\u003eFront Neurol\u003c/em\u003e. \u003cstrong\u003e4\u003c/strong\u003e(15),1383757(2024).\u003c/li\u003e\n\u003cli\u003eZhou, H. et al. Epidemiological and clinical features, treatment status, and economic burden of traumatic spinal cord injury in China: a hospital-based retrospective study. \u003cem\u003eNeural Regen Res\u003c/em\u003e. \u003cstrong\u003e19\u003c/strong\u003e(5),1126-1133(2024).\u003c/li\u003e\n\u003cli\u003eEli,I., Lerner,D.P., \u0026amp; Ghogawala,Z. Acute traumatic spinal cord injury. \u003cem\u003eNeurol Clin\u003c/em\u003e. \u003cstrong\u003e39\u003c/strong\u003e(2):471-488(2021).\u003c/li\u003e\n\u003cli\u003eAnjum, A. et al. Spinal cord injury: pathophysiology, multimolecular interactions, and underlying recovery mechanisms. \u003cem\u003eInt J Mol Sci\u003c/em\u003e. \u003cstrong\u003e21\u003c/strong\u003e(20),7533 (2020).\u003c/li\u003e\n\u003cli\u003eWiderstr\u0026ouml;m-Noga,E. Neuropathic pain and spinal cord injury: management, phenotypes, and biomarkers. \u003cem\u003eDrugs\u003c/em\u003e. \u003cstrong\u003e83\u003c/strong\u003e(11),1001-1025(2023).\u003c/li\u003e\n\u003cli\u003eOtzel,D.M., Kok,H.J., Graham,Z.A., Barton,E.R. \u0026amp; Yarrow,J.F. Pharmacologic approaches to prevent skeletal muscle atrophy after spinal cord injury. \u003cem\u003eCurr Opin Pharmacol\u003c/em\u003e. \u003cstrong\u003e10\u003c/strong\u003e(60),193-199(2021).\u003c/li\u003e\n\u003cli\u003eChandrasekaran, S., Davis,J., Bersch,I., Goldberg,G. \u0026amp; Gorgey,A.S. Electrical stimulation and denervated muscles after spinal cord injury. \u003cem\u003eNeural Regen Res\u003c/em\u003e. \u003cstrong\u003e15\u003c/strong\u003e(8),1397-1407(2020).\u003c/li\u003e\n\u003cli\u003eXu, X.et al. Mechanism of skeletal muscle atrophy after spinal cord injury: a narrative review. \u003cem\u003eFront Nutr\u003c/em\u003e. \u003cstrong\u003e3\u003c/strong\u003e(10),1099143(2023).\u003c/li\u003e\n\u003cli\u003eGorgey, A.S.et al.Mitochondrial health and muscle plasticity after spinal cord injury. \u003cem\u003eEur J Appl Physiol\u003c/em\u003e. \u003cstrong\u003e119\u003c/strong\u003e(2),315-331(2019).\u003c/li\u003e\n\u003cli\u003ePotter, L.A.et al.Assessing the impact of boldine on the gastrocnemius using multiomics profiling at 7 and 28 days post-complete spinal cord injury in young male mice. \u003cem\u003ePhysiol Genomic\u003c/em\u003es. \u003cstrong\u003e55\u003c/strong\u003e(7),297-313(2023). \u003c/li\u003e\n\u003cli\u003eLin, X. et al. Pancancer analyses reveal genomics and clinical characteristics of the SETDB1 in human tumors. \u003cem\u003eJ Oncol\u003c/em\u003e. 2022,6115878(2022). \u003c/li\u003e\n\u003cli\u003eRitchie, M.E.et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. \u003cem\u003eNucleic Acids Res\u003c/em\u003e. \u003cstrong\u003e43\u003c/strong\u003e(7),e47(2015). \u003c/li\u003e\n\u003cli\u003eLangfelder, P. \u0026amp; Horvath,S. WGCNA: an R package for weighted correlation network analysis. \u003cem\u003eBMC Bioinformatics\u003c/em\u003e. \u003cstrong\u003e29\u003c/strong\u003e(9),559(2008).\u003c/li\u003e\n\u003cli\u003eAltermann,E. \u0026amp; Klaenhammer, T.R. PathwayVoyager: pathway mapping using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. \u003cem\u003eBMC Genomics\u003c/em\u003e. 3(6),60(2005).\u003c/li\u003e\n\u003cli\u003eYu,G., Wang,L.G., Han,Y. \u0026amp; He,Q.Y. clusterProfiler: an R package for comparing biological themes among gene clusters. \u003cem\u003eOMICS\u003c/em\u003e. \u003cstrong\u003e16\u003c/strong\u003e(5),284-7(2012).\u003c/li\u003e\n\u003cli\u003eSherman,B.T. et al. DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). \u003cem\u003eNucleic Acids Res\u003c/em\u003e. \u003cstrong\u003e50\u003c/strong\u003e(W1),W216-W221(2022).\u003c/li\u003e\n\u003cli\u003eSubramanian,A.et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. \u003cem\u003eProc Natl Acad Sci USA\u003c/em\u003e. \u003cstrong\u003e102\u003c/strong\u003e(43),15545-50(2005).\u003c/li\u003e\n\u003cli\u003eLiberzon,A. et al. The Molecular Signatures Database (MSigDB) hallmark gene set collection. \u003cem\u003eCell Syst\u003c/em\u003e. \u003cstrong\u003e1\u003c/strong\u003e(6),417-425(2015).\u003c/li\u003e\n\u003cli\u003eDoncheva,N.T., Morris,J.H., Gorodkin,J. \u0026amp;Jensen,L.J. Cytoscape StringApp: network analysis and visualization of proteomics data. \u003cem\u003eJ Proteome Res\u003c/em\u003e. \u003cstrong\u003e18\u003c/strong\u003e(2),623-632(2019).\u003c/li\u003e\n\u003cli\u003eNemeth,C., Banik,N.L. \u0026amp; Haque,A. Disruption of neuromuscular junction following spinal cord injury and motor neuron diseases. \u003cem\u003eInt J Mol Sci\u003c/em\u003e. \u003cstrong\u003e25\u003c/strong\u003e(6),3520(2024).\u003c/li\u003e\n\u003cli\u003eRosales-Antequera,C., Viscor,G. \u0026amp; Araneda,O.F. Inflammation and oxidative stress as common mechanisms of pulmonary, autonomic and musculoskeletal dysfunction after spinal cord injury. \u003cem\u003eBiology (Basel)\u003c/em\u003e. \u003cstrong\u003e11\u003c/strong\u003e(4),550(2022).\u003c/li\u003e\n\u003cli\u003eWei,L., Cai,G., Jiang,L., Gao,L., Yang,Z. \u0026amp; Zhang,W. Identification of key pathways and RNAs associated with skeletal muscle atrophy after spinal cord injury. \u003cem\u003eJ Musculoskelet Neuronal Interact\u003c/em\u003e. 21(4):550-559(2021).\u003c/li\u003e\n\u003cli\u003eZhang,C. et al. Early electrical field stimulation prevents the loss of spinal cord anterior horn motoneurons and muscle atrophy following spinal cord injury. \u003cem\u003eNeural Regen Res\u003c/em\u003e. \u003cstrong\u003e13\u003c/strong\u003e(5),869-876(2018).\u003c/li\u003e\n\u003cli\u003eDionyssiotis,Y., Mavrogenis,A.,Trovas,G., Skarantavos,G., Papathanasiou,J. \u0026amp; Papagelopoulos,P.Bone and soft tissue changes in patients with spinal cord injury and multiple sclerosis. \u003cem\u003eFolia Med (Plovdiv)\u003c/em\u003e. \u003cstrong\u003e56\u003c/strong\u003e(4),237-44(2014).\u003c/li\u003e\n\u003cli\u003eXue,J.et al. The identification and validation of hub genes associated with acute myocardial infarction using weighted gene co-expression network analysis. \u003cem\u003eJ Cardiovasc Dev Dis\u003c/em\u003e. \u003cstrong\u003e9\u003c/strong\u003e(1),30(2022). \u003c/li\u003e\n\u003cli\u003eGaudet,P., Logie,C., Lovering,R.C., Kuiper,M., L\u0026aelig;greid,A. \u0026amp; Thomas PD. Gene Ontology representation for transcription factor functions. \u003cem\u003eBiochim Biophys Acta Gene Regul Mech\u003c/em\u003e. \u003cstrong\u003e1864\u003c/strong\u003e(11-12),194752(2021).\u003c/li\u003e\n\u003cli\u003eChen,L., Zhang,Y.H., Wang,S., Zhang,Y., Huang,T. \u0026amp; Cai,Y.D.Prediction and analysis of essential genes using the enrichments of gene ontology and KEGG pathways. \u003cem\u003ePLoS One\u003c/em\u003e. \u003cstrong\u003e12\u003c/strong\u003e(9),e0184129(2017).\u003c/li\u003e\n\u003cli\u003eZito,A.et al. Gene set enrichment analysis of interaction networks weighted by node centrality. \u003cem\u003eFront Genet\u003c/em\u003e. \u003cstrong\u003e24\u003c/strong\u003e(12),577623(2021).\u003c/li\u003e\n\u003cli\u003eEsposito,E. \u0026amp; Cuzzocrea,S. Targeting the peroxisome proliferator-activated receptors (PPARs) in spinal cord injury. \u003cem\u003eExpert Opin Ther Targets\u003c/em\u003e. \u003cstrong\u003e15\u003c/strong\u003e(8),943-59(2011).\u003c/li\u003e\n\u003cli\u003eLi,Z.et al. A proteomic and phosphoproteomic landscape of spinal cord injury. \u003cem\u003eNeurosci Lett\u003c/em\u003e. \u003cstrong\u003e25\u003c/strong\u003e(814),137449(2023).\u003c/li\u003e\n\u003cli\u003eZhou,Y.et al.Phenylbutyrate prevents disruption of blood-spinal cord barrier by inhibiting endoplasmic reticulum stress after spinal cord injury. \u003cem\u003eAm J Transl Res\u003c/em\u003e. \u003cstrong\u003e8\u003c/strong\u003e(4),1864-75(2016). \u003c/li\u003e\n\u003cli\u003eKang,C. \u0026amp; Li,J.L. Role of PGC-1\u0026alpha; signaling in skeletal muscle health and disease. \u003cem\u003eAnn NYAcad Sci\u003c/em\u003e. \u003cstrong\u003e1271\u003c/strong\u003e(1),110-7(2012). \u003c/li\u003e\n\u003cli\u003eWang,J., Liu,Y., Zhang,Y., Liu,B. \u0026amp; Wei,Z. Identification of potential miRNA-mRNA regulatory network in denervated muscular atrophy by bioinformatic analysis. \u003cem\u003eBiomed Res Int\u003c/em\u003e. 28;2022,6042591(2022).\u003c/li\u003e\n\u003cli\u003eDoncheva,N.T.et al. Cytoscape stringApp 2.0: Analysis and visualization of heterogeneous biological networks. \u003cem\u003eJ Proteome Res\u003c/em\u003e. 2023;\u003cstrong\u003e22\u003c/strong\u003e(2):637-646. \u003c/li\u003e\n\u003cli\u003eCordier,F. \u0026amp; Creytens,D. New kids on the block: FOS and FOSB gene. \u003cem\u003eJ Clin Pathol\u003c/em\u003e. \u003cstrong\u003e22\u003c/strong\u003e(2),637-646(2023).\u003c/li\u003e\n\u003cli\u003eChoi,H.et al. Amyloid \u0026beta;-induced elevation of O-GlcNAcylated c-Fos promotes neuronal cell death. \u003cem\u003eAging Cell\u003c/em\u003e. \u003cstrong\u003e22\u003c/strong\u003e(2),637-646(2023).\u003c/li\u003e\n\u003cli\u003eZhang,D., Liu,J., Zhu,T. \u0026amp; Zhou,C. Identifying c-fos expression as a strategy to investigate the actions of general anesthetics on the cCentral nervous system. \u003cem\u003eCurr Neuropharmacol\u003c/em\u003e. \u003cstrong\u003e20\u003c/strong\u003e(1),55-71(2022).\u003c/li\u003e\n\u003cli\u003eXu,Y., Zheng,Z., Ho,K.P. \u0026amp; Qian,Z. Effects of spinal cord injury on c-fos expression in hypothalamic paraventricular nucleus and supraoptic nucleus in rats. \u003cem\u003eBrain Res\u003c/em\u003e. \u003cstrong\u003e1087\u003c/strong\u003e(1):175-9(2006).\u003c/li\u003e\n\u003cli\u003eCambier,S. et al. The chemokines CXCL8 and CXCL12: molecular and functional properties, role in disease and efforts towards pharmacological intervention. Cell Mol Immunol. \u003cstrong\u003e20\u003c/strong\u003e(3),217-251(2023).\u003c/li\u003e\n\u003cli\u003eCambier,S., Gouwy,M. \u0026amp; Proost,P. Bone marrow-derived monocytes drive the inflammatory microenvironment in local and remote regions after thoracic spinal cord injury. \u003cem\u003eJ Neurotrauma\u003c/em\u003e. \u003cstrong\u003e36\u003c/strong\u003e(6),937-949(2019).\u003c/li\u003e\n\u003cli\u003eDansereau,M.A.et al. Mechanistic insights into the role of the chemokine CCL2/CCR2 axis in dorsal root ganglia to peripheral inflammation and pain hypersensitivity. \u003cem\u003eJ Neuroinflammation\u003c/em\u003e. \u003cstrong\u003e18\u003c/strong\u003e(1),79(2021).\u003c/li\u003e\n\u003cli\u003eFerrara,P.J.et al. Global deletion of CCL2 has adverse impacts on recovery of skeletal muscle fiber size and function and is muscle specific. \u003cem\u003eJ Appl Physiol\u003c/em\u003e. \u003cstrong\u003e134\u003c/strong\u003e(4),923-932(2023).\u003c/li\u003e\n\u003cli\u003eSartori,R. et al. Smad2 and 3 transcription factors control muscle mass in adulthood. \u003cem\u003eAm J Physiol Cell Physiol\u003c/em\u003e. \u003cstrong\u003e296\u003c/strong\u003e(6):C1248-57(2009).\u003c/li\u003e\n\u003cli\u003eHofmann,C. et al. Transient inhibition of translation improves cardiac function after ischemia/reperfusion by attenuating the inflammatory response. \u003cem\u003eCirculation\u003c/em\u003e. 150(16):1248-1267(2024).\u003c/li\u003e\n\u003cli\u003ePlantevin,K.V.et al. Discovery of CC-930, an orally active anti-fibrotic JNK inhibitor. \u003cem\u003eBioorg Med Chem Lett\u003c/em\u003e.\u003cstrong\u003e 22\u003c/strong\u003e(3):1433-8(2012).\u003c/li\u003e\n\u003cli\u003eShi,K. et al. tPA mobilizes immune cells that exacerbate hemorrhagic transformation in stroke. \u003cem\u003eCirc Res\u003c/em\u003e. \u003cstrong\u003e128\u003c/strong\u003e(1):62-7(2021).\u003c/li\u003e\n\u003cli\u003eZhou,L., Xue,C., Chen,Z., Jiang,W., He,S. \u0026amp; Zhang,X. c-Fos is a mechanosensor that regulates inflammatory responses and lung barrier dysfunction during ventilator-induced acute lung injury. \u003cem\u003eBMC Pulm Med\u003c/em\u003e. \u003cstrong\u003e22\u003c/strong\u003e(1):9(2022).\u003c/li\u003e\n\u003cli\u003eRaghu,H. et al. CCL2/CCR2, but not CCL5/CCR5, mediates monocyte recruitment, inflammation and cartilage destruction in osteoarthritis. \u003cem\u003eAnn Rheum Dis\u003c/em\u003e. \u003cstrong\u003e76\u003c/strong\u003e(5):914-92(2017).\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Differentially expressed genes, Spinal cord injury, Muscle atrophy, Hub genes, Bioinformatics analysis, Weighted gene co-expression network analysis","lastPublishedDoi":"10.21203/rs.3.rs-5504847/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-5504847/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eSpinal cord injury (SCI) is a neurological disease with high morbidity and mortality. Post-SCI muscle atrophy is a cascade response to SCI, and failure to actively prevent its occurrence severely affects patients' mobility and quality of life. Therefore, deeply exploring the correlation between muscle atrophy after SCI and the molecular regulation mechanism is of great significance. Download GSE21497 expression profile data from the gene expression omnibus(GEO)database. Perform weighted gene co-expression network analysis (WGCNA) on the obtained differentially expressed genes (DEGs).Subsequently, we performed functional and pathway enrichment analyses of key modules. Construct a protein-protein interaction (PPI) network and screen core genes. Finally, the results were verified by real-time polymerase chain reaction(PCR). A total of 1007 DEGs were obtained, including 533 upregulated genes and 474 downregulated genes. WGCNA analysis identified 161 turquoise modules of DEGs as key modules related to SCI. Functional enrichment analysis showed that their involvement in multiple processes and pathways. FOS and CCL2 may be involved in the molecular pathophysiology of muscle atrophy after SCI, serving as potential targets for diagnosis or treatment of SCI-related muscle atrophy.\u003c/p\u003e","manuscriptTitle":"Identification of key genes associated with muscle atrophy after spinal cord injury and experimental verification in rats","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-04-17 05:33:49","doi":"10.21203/rs.3.rs-5504847/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"d4b61ca4-13b0-4693-b31e-f98481a3fded","owner":[],"postedDate":"April 17th, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":47111532,"name":"Biological sciences/Computational biology and bioinformatics"},{"id":47111533,"name":"Biological sciences/Neuroscience"},{"id":47111534,"name":"Health sciences/Biomarkers"},{"id":47111535,"name":"Health sciences/Medical research"}],"tags":[],"updatedAt":"2025-07-15T21:53:14+00:00","versionOfRecord":[],"versionCreatedAt":"2025-04-17 05:33:49","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-5504847","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-5504847","identity":"rs-5504847","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}