Identification of gene coexpression modules, hub genes, and pathways related to spinal cord injury using integrated bioinformatics methods

Inhibition of IL1R1 or CASP4 attenuates spinal cord injury through ameliorating NLRP3 inflammasome-induced pyroptosis

Spinal cord injury (SCI) is a devastating trauma characterized by serious neuroinflammation and permanent neurological dysfunction. However, the molecular mechanism of SCI remains unclear, and few effective medical therapies are available at present. In this study, multiple bioinformatics methods were used to screen out novel targets for SCI, and the mechanism of these candidates during the progression of neuroinflammation as well as the therapeutic effects were both verified in a rat model of traumatic SCI. As a result, CASP4, IGSF6 and IL1R1 were identified as the potential diagnostic and therapeutic targets for SCI by computational analysis, which were enriched in NF-κB and IL6-JAK-STATA3 signaling pathways. In the injured spinal cord, these three signatures were up-regulated and closely correlated with NLRP3 inflammasome formation and gasdermin D (GSDMD) -induced pyroptosis. Intrathecal injection of inhibitors of IL1R1 or CASP4 improved the functional recovery of SCI rats and decreased the expression of these targets and inflammasome component proteins, such as NLRP3 and GSDMD. This treatment also inhibited the pp65 activation into the nucleus and apoptosis progression. In conclusion, our findings of the three targets shed new light on the pathogenesis of SCI, and the use of immunosuppressive agents targeting these proteins exerted anti-inflammatory effects against spinal cord inflammation by inhibiting NF-kB and NLRP3 inflammasome activation, thus blocking GSDMD -induced pyroptosis and immune activation.

Identification of Regeneration and Hub Genes and Pathways at Different Time Points after Spinal Cord Injury

Spinal cord injury (SCI) is a neurological injury that can cause neuronal loss around the lesion site and leads to locomotive and sensory deficits. However, the underlying molecular mechanisms remain unclear. This study aimed to verify differential gene time-course expression in SCI and provide new insights for gene-level studies. We downloaded two rat expression profiles (GSE464 and GSE45006) from the Gene Expression Omnibus database, including 1 day, 3 days, 7 days, and 14 days post-SCI, along with thoracic spinal cord data for analysis. At each time point, gene integration was performed using "batch normalization." The raw data were standardized, and differentially expressed genes at the different time points versus the control were analyzed by Gene Ontology enrichment analysis, the Kyoto Encyclopedia of Genes and Genomes pathway analysis, and gene set enrichment analysis. A protein-protein interaction network was then built and visualized. In addition, ten hub genes were identified at each time point. Among them, Gnb5, Gng8, Agt, Gnai1, and Psap lack correlation studies in SCI and deserve further investigation. Finally, we screened and analyzed genes for tissue repair, reconstruction, and regeneration and found that Anxa1, Snap25, and Spp1 were closely related to repair and regeneration after SCI. In conclusion, hub genes, signaling pathways, and regeneration genes involved in secondary SCI were identified in our study. These results may be useful for understanding SCI-related biological processes and the development of targeted intervention strategies.

Identification and analysis of spinal cord injury subtypes using weighted gene co-expression network analysis

Background

Spinal cord injury (SCI) has an immediate and devastating impact on the control over various movements and sensations. However, no effective therapies for SCI currently exist.

Methods

To identify and analyze SCI subtypes, we obtained the expression profile data of the 1,057 genes (889 intersection genes) in GSE45550 using weighted gene co-expression network analysis (WGCNA), and 14 co-expression gene modules were identified. Next, we filtered out the network degree top 10 (degree >80) genes, considered the final key SCI genes. A multifactor regulatory network (105 interaction pairs), consisting of messenger RNAs (mRNAs), long non-coding RNAs (lncRNAs), and transcription factors (TFs) was constructed. This network was involved in the co-expression of key genes. We selected the top 10 regulatory factors (degree >4) as core regulators in the multifactor regulatory network.

Results

The results of functional enrichment analysis of the target gene expressing the core regulatory factor [1,059] showed that these target genes were enriched in pathways for human cytomegalovirus infection, chronic myeloid leukemia, and pancreatic cancer. Further, we used the key genes in the co-expression network to categorize the SCI samples in GSE45550. The expression levels of the top 6 genes (CCNB2, CCNB1, CKS2, COL5A1, KIF20A, and RACGAP1) may act as potential marker genes for different SCI subtypes. On the basis of these different subtypes, 8 SCI core gene CDK1-associated drugs were also found to provide potential therapeutic options for SCI.

Conclusions

These results may provide a novel therapeutic strategy for the treatment of SCI.

Integrated Bioinformatics Analysis of Hub Genes and Pathways Associated with a Compression Model of Spinal Cord Injury in Rats

Background

Spinal cord injury (SCI) is a serious nervous system condition that can cause lifelong disability. The aim of this study was to identify potential molecular mechanisms and therapeutic targets for SCI.

Material/Methods

We constructed a weighted gene coexpression network and predicted which hub genes are involved in SCI. A compression model of SCI was established in 45 Sprague-Dawley rats, which were divided into 5 groups (n=9 per group): a sham operation group, and 1, 3, 5, and 7 days post-SCI groups. The spinal cord tissue on the injured site was harvested on 1, 3, 5, and 7 days after SCI and 3 days after surgery in the sham operation group. High-throughput sequencing was applied to investigate the expression profile of the mRNA in all samples. Differentially expressed genes were screened and included in weighted gene coexpression network analysis (WGCNA). Co-expressed modules and hub genes were identified by WGCNA. The biological functions of each module were investigated using the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases.

Results

According to the RNA-seq data, a total of 1965 differentially expressed genes were screened, and WGCNA identified 10 coexpression modules and 5 hub genes. Module function analysis revealed that SCI was associated with immune response, cell division, neuron projection development, and collagen fibril organization.

Conclusions

Our study revealed dynamic changes in a variety of biological processes following SCI and identified 5 hub genes via WGCNA. These results provide insights into the molecular mechanisms and therapeutic targets of SCI.

Identification of Key Modules, Hub Genes, and Noncoding RNAs in Chronic Rhinosinusitis with Nasal Polyps by Weighted Gene Coexpression Network Analysis

Chronic rhinosinusitis with nasal polyps (CRSwNP) is a chronic inflammatory disease with relatively easy recurrence. However, the precise molecular mechanisms of this disease are poorly known. Based on gene sequencing data obtained from the Gene Expression Omnibus (GEO) database, we constructed coexpression networks by weighted gene coexpression network analysis (WGCNA). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed by the Database for Annotation, Visualization, and Integrated Discovery (DAVID). The core gene of pathogenesis, CRSwNP, was screened by protein-protein interaction data (PPI) from the HPRD database. Unsupervised clustering was applied to screen hub genes related to the phenotype of CRSwNP. Blue and turquoise modules were found to be most significantly related to the pathogenicity of CRSwNP. Functional enrichment analysis showed that cell proliferation in the blue modules, the apoptotic process in the turquoise module, and the cancer pathway in both modules were mostly significantly correlated with the development of CRSwNP. The noncoding RNAs (long noncoding RNA and microRNA) and the top 10 core genes in each module were found to be associated with the pathogenesis of CRSwNP. A total of nine hub genes were identified to be related to the CRSwNP phenotype. By qRT-PCR analysis, AKT1, CDH1, PIK3R1, CBL, LRP1, MALAT1, and XIST were proven to be associated with the pathogenesis of CRSwNP. AGR2, FAM3D, PIP, DSE, and TMC were identified to be related to the CRSwNP phenotype. Further exploration of these genes will reveal more important information about the mechanisms of CRSwNP.

Identification of RAD54 homolog B as a promising therapeutic target for breast cancer

Breast cancer is a recognized threat to the health of women globally. Due to the lack of the knowledge about the molecular pathogenesis of breast cancer, therapeutic strategies remain inadequate, especially for aggressive breast cancer. In the present study, sequential bioinformatics analysis was performed using data from the GSE20711 dataset, and the results demonstrated that three genes may impact the survival of patients with breast cancer. One of these genes, RAD54 homolog B (RAD54B), may be a potential prognostic factor for breast cancer. A signature was established that could evaluate the overall survival for patients with breast cancer based on the risk score calculated from RAD54B expression and the Tumor-Node-Metastasis (TNM) stage [risk score=expRAD54B × 0.236 + TNM stage (I/II=0 or III/IV=1) ×1.025]. In addition, based on the GSE85871 dataset and inhibitory assay, the study identified a natural compound, Japonicone A, which may reduce the proliferation of breast cancer cells by inhibiting the expression of RAD54B. Overall, the present study identified a novel candidate gene and a candidate compound as promising therapeutic targets for the treatment of breast cancer.