Hub Genes and Key Pathways of Intervertebral Disc Degeneration: Bioinformatics Analysis and Validation.
Ontology highlight
ABSTRACT: Objective
To identify significant pathways and genes in intervertebral disc degeneration (IDD) based on bioinformatics analysis.Design
The GEO database was used to download the GSE124272 dataset. Differentially expressed genes (DEGs) were analyzed using Limma package in R language. Then, gene ontologies (GO), Kyoto encyclopedia of genes and genomes (KEGG), and protein-protein interaction (PPI) networks were used to further identify hub genes. The mRNA expression levels of top six hub genes were verified.Results
We found 563 DEGs, of which 214 were upregulated and 349 were downregulated. The top 5 GO terms and pathways were shown including immune response, cell cycle, and p53 pathway. Based on the PPI analysis, we verified the mRNA expression levels of 6 hub genes. The mRNA levels of CHEK1, CDCA2, SKA3, and KIF20A were upregulated in degenerative NP tissue than in healthy NP tissue. However, the mRNA level of BUB1 and SPC25 was downregulated.Conclusions
This study may provide new biomarkers for the IDD and treatments to repair IDD related to CHEK1, CDCA2, SKA3, BUB1, KIF20A, and SPC25.
SUBMITTER: Zhang Z
PROVIDER: S-EPMC8449732 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Background: Intervertebral disc degeneration impairs the quality of patients lives. Even though there has been development of many therapeutic strategies, most of them remain unsatisfactory due to the limited understanding of the mechanisms that underlie the intervertebral disc degeneration. Questions/purposes: This study is meant to identify the key modules and hub genes related to the annulus fibrosus in intervertebral disc degeneration (IDD) through: (1) constructing a weighted gene co-expression network; (2) identifying key modules and hub genes; (3) verifying the relationships of key modules and hub genes with IDD; and (4) confirming the expression pattern of hub genes in clinical samples. Methods: The Gene Expression Omnibus provided 24 sets of annulus fibrosus microarray data. Differentially expressed genes between the annulus fibrosus of degenerative and non-degenerative intervertebral disc samples have gone through the Gene Ontology (GO) and pathway analysis. The construction of a gene network and classification of genes into different modules were conducted through performing Weighted Gene Co-expression Network Analysis. The identification of modules and hub genes that were most related to intervertebral disc degeneration was proceeded. In order to verify the relationships of the module and hub genes with intervertebral disc degeneration, Ingenuity Pathway Analysis was operated. Clinical samples were adopted to help verify the hub gene expression profile. Results: One thousand one hundred ninety differentially expressed genes were identified. Terms and pathways associated with intervertebral disc degeneration were presented by GO and pathway analysis. The construction of a Weighted Gene Coexpression Network was completed and clustering differentially expressed genes into four modules was also achieved. The module with the lowest P-value and the highest absolute correlation coefficient was selected and its relationship with intervertebral disc degeneration was confirmed by Ingenuity Pathway Analysis. The identification of hub genes and the confirmation of their expression profile were also realized. Conclusions: This study generated a comprehensive overview of the gene networks underlying annulus fibrosus in intervertebral disc degeneration. Clinical Relevance: Modules and hub genes identified in this study are highly associated with intervertebral disc degeneration, and may serve as potential therapeutic targets for intervertebral disc degeneration.
| S-EPMC7875098 | biostudies-literature
Project description:BackgroundIntervertebral disc degeneration (IDD) is widely known as the main contributor to low back pain which has a negative socioeconomic impact worldwide. However, the underlying mechanism remains unclear. This study aims to analyze the dataset GSE23130 using bioinformatics methods to identify the pivotal genes and pathways associated with IDD.Material/methodsThe gene expression data of GSE23130 was downloaded, and differentially expressed genes (DEGs) were extracted from 8 samples and 15 controls. GO and KEGG pathway enrichment analyses were performed. Also, protein-protein interaction (PPI) network was constructed and visualized, followed by identification of hub genes and key module.ResultsA total of 30 downregulated and 79 upregulated genes were identified. The DEGs were mainly enriched in the regulation of protein catabolic process, extracellular matrix organization, collagen fibril organization, and extracellular structure organization. Meanwhile, we found that most DEGs were primarily enriched in the PI3K-Akt signaling pathway. The top 10 hub genes were FN1, COL1A2, SPARC, COL3A1, CTGF, LUM, TIMP1, THBS2, COL5A2, and TGFB1.ConclusionsIn summary, key candidate genes and pathways were identified by using integrated bioinformatics analysis, which may provide insights into the underlying mechanisms and offer potential target genes for the treatment of IDD.
| S-EPMC7465721 | biostudies-literature
Project description:The exact mechanisms associated with secondary brain damage following traumatic brain injury (TBI) remain unclear; therefore, identifying the critical molecular mechanisms involved in TBI is essential. The mRNA expression microarray GSE2871 was downloaded from the Gene Expression Omnibus (GEO) repository. GSE2871 comprises a total of 31 cerebral cortex samples, including two post-TBI time points. The microarray features eight control and seven TBI samples, from 4 hours post-TBI, and eight control and eight TBI samples from 24 hours post-TBI. In this bioinformatics-based study, 109 and 66 differentially expressed genes (DEGs) were identified in a Sprague-Dawley (SD) rat TBI model, 4 and 24 hours post-TBI, respectively. Functional enrichment analysis showed that the identified DEGs were significantly enriched in several terms, such as positive regulation of nuclear factor-?B transcription factor activity, mitogen-activated protein kinase signaling pathway, negative regulation of apoptotic process, and tumor necrosis factor signaling pathway. Moreover, the hub genes with high connectivity degrees were primarily related to inflammatory mediators. To validate the top five hub genes, a rat model of TBI was established using the weight-drop method, and real-time quantitative polymerase chain reaction analysis of the cerebral cortex was performed. The results showed that compared with control rats, Tnf-?, c-Myc, Spp1, Cxcl10, Ptprc, Egf, Mmp9, and Lcn2 were upregulated, and Fn1 was downregulated in TBI rats. Among these hub genes, Fn1, c-Myc, and Ptprc may represent novel biomarkers or therapeutic targets for TBI. These identified pathways and key genes may provide insights into the molecular mechanisms of TBI and provide potential treatment targets for patients with TBI. This study was approved by the Experimental Animal Ethics Committee of the First Affiliated Hospital of Nanchang University, China (approval No. 003) in January 2016.
| S-EPMC7749465 | biostudies-literature
Project description:Although the pathology of sciatica has been studied extensively, the transcriptional changes in the peripheral blood caused by sciatica have not been characterized. An Agilent SurePrint G3 human gene expression microarray 8 ╳ 60 K was used to analyse gene expression in whole blood from patients with sciatica compared with healthy volunteers.
2021-01-22 | GSE150408 | GEO
Project description:This SuperSeries is composed of the SubSeries listed below.
2024-04-05 | GSE230809 | GEO
Project description:Background and objectiveIntervertebral disc degeneration (IDD) is a complex multifactorial and irreversible pathological process. In IDD, multiple competing endogenous RNAs (ceRNA, including mRNA, lncRNA, and pseudogenes) can compete to bind with miRNAs. However, the potential metabolic signatures in nucleus pulposus (NP) cells remain poorly understood. This study investigated key metabolic genes and the ceRNA regulatory mechanisms in the pathogenesis of IDD based on microarray datasets.MethodsWe retrieved and downloaded four independent IDD microarray datasets from the Gene Expression Omnibus. Combining the predicted interactions from online databases (miRcode, miRDB, miRTarBase, and TargetScan), differentially expressed lncRNAs (DElncRNAs), miRNAs (DEmiRNAs), and mRNAs (DEmRNAs) were identified. A ceRNA network was constructed and annotated using GO and KEGG pathway enrichment analyses. Moreover, we searched the online metabolic gene set and used support vector machine (SVM) to find the critical metabolic DEmRNA(s) and other DERNAs. Differential gene expression was validated with a merged dataset.ResultsA total of 45 DEmRNAs, 36 DElncRNAs, and only one DEmiRNA (miR-338-3p) were identified in the IDD microarray datasets. GO and KEGG pathway enrichment analyses revealed that the DEmRNAs were predominantly enriched in the PI3K-Akt signaling pathway, MAPK signaling pathway, IL-17 signaling pathway, apoptosis, and cellular response to oxidative stress. Based on SVM screening, 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase (PFK/FBPase) 2 is the critical metabolic gene with lower expression in IDD, and AC063977.6 is the key lncRNA with lower expression in IDD. The ceRNA hypothesis suggests that AC063977.6, miR-338-3p (high expression), and PFKFB2 are dysregulated as an axis in IDD.ConclusionsThe results suggest that lncRNA AC063977.6 correlate with PFKFB2, the vital metabolic signature gene, via targeting miR-338-3p during IDD pathogenesis. The current study may shed light on unraveling the pathogenesis of IDD.
| S-EPMC7937257 | biostudies-literature
Project description:Background:Basal-like breast cancer (BLBC) is the most aggressive subtype of breast cancer (BC) and links to poor outcomes. As the molecular mechanism of BLBC has not yet been completely discovered, identification of key pathways and hub genes of this disease is an important way for providing new insights into exploring the mechanisms of BLBC initiation and progression. Objective:The aim of this study was to identify potential gene signatures of the development and progression of the BLBC via bioinformatics analysis. Methods and results:The differential expressed genes (DEGs) including 40 up-regulated and 21 down-regulated DEGs were identified between GSE25066 and GSE21422 microarrays, and these DEGs were significantly enriched in the terms related to oncogenic or suppressive roles in BLBC progression. In addition, KEGG pathway and GSEA (Gene Set Enrichment Analysis) enrichment analyses were performed for DEGs between the basal type and non-basal-type breast cancer from GSE25066 microarray. These DEGs were enriched in pathways such as cell cycle, cytokine-cytokine receptor interaction, chemokine signaling pathway, central carbon metabolism signaling and TNF signaling pathway. Moreover, the protein-protein interaction (PPI) network was constructed with those 61 DEGs using the Cytoscape software, and the biological significance of putative modules was established using MCODE. The module 1 was found to be closely related with a term of mitosis regulation and enriched in cell cycle pathway, and thus confirmed the pathological characteristic of BLBC with a high mitotic index. Furthermore, prediction values of the top 10 hub genes such as CCNB2, BUB1, NDC80, CENPE, KIF2C, TOP2A, MELK, TPX2, CKS2 and KIF20A were validated using Oncomine and Kaplan-Meier plotter. Conclusion:Our results suggest the intriguing possibility that the hub genes and modules in the PPI network contributed to in-depth knowledge about the molecular mechanism of BLBC, paving a way for more accurate discovery of potential treatment targets for BLBC patients.
| S-EPMC6388944 | biostudies-literature
Project description:Non-small cell lung cancer (NSCLC) is the most common type of lung cancer, accounting for ~80% of all lung cancer cases. The aim of the present study was to identify key genes and pathways in NSCLC, in order to improve understanding of the mechanism of lung cancer. The GSE33532 gene expression dataset, containing 20 normal and 80 NSCLC samples, was used. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to obtain the enrichment data of differently expressed genes (DEGs). Disease modules within NSCLC were constructed by Cytoscape, using protein-protein interaction (PPI) from the Search Tool for the Retrieval of Interacting Genes database. In addition, the Kaplan Meier plotter KMplot was used to assess the top hub genes in the PPI network. As a result, 1,795 genes were identified in NSCLC; 729 were upregulated and 1,066 were downregulated. The results of the GO analysis indicated that the upregulated DEGs were significantly enriched in 'biological processes' (BP), including 'cell cycle and nuclear division'; the downregulated DEGs were also significantly enriched in BP, including 'response to wounding', 'anatomical structure morphogenesis' and 'response to stimulus'. Upregulated DEGs were also enriched in 'cell cycle', 'DNA replication' and the 'tumor protein 53 signaling pathway', while the downregulated DEGs were also enriched in 'complement and coagulation cascades', 'malaria' and 'cell adhesion molecules'. The top 9 hub genes were cyclin-dependent kinase 9 (CDK1), polo-like kinase 1, aurora kinase B, cell division cycle 20, baculoviral initiator of apoptosis repeat containing 5, mitotic checkpoint serine/threonine kinase B, proliferating cell nuclear antigen (PCNA), centromere protein A and MAD2 mitotic arrest deficient-like 1, and the KMplot results revealed that the high expression levels of these genes resulted in significantly low survival rates, compared with low expression samples (P<0.05), with the exception of PCNA and CDK1. In the pathway crosstalk analysis, 26 nodes and 41 interactions were divided into two groups: One module of the two groups primarily included 'metabolism of amino acid' and the other primarily contained 'tumor necrosis signaling' pathways. In conclusion, the present study assisted in improving the understanding of the molecular mechanisms underlying NSCLC development, and the results may help the understanding of the biological mechanism of NSCLC.
| S-EPMC6036325 | biostudies-literature
Project description:BackgroundThe competing endogenous RNA- (ceRNA-) mediated regulatory mechanisms are known to play a pivotal role in intervertebral disc degeneration (IDD). Our research intended to establish a ceRNA regulatory network related to IDD through bioinformatics analyses.MethodsThe expression profiles of circRNA, miRNA, and mRNA were obtained from the public Gene Expression Omnibus (GEO) datasets. Then, we use sequence-based bioinformatics methods to select differentially expressed mRNAs (DEmRNAs), microRNAs (DEmiRNAs), or circRNAs (DEcircRNAs) related to IDD. We used ChEA3 to verify the targets of transcription factors (TFs). Then, we used DAVID to annotate the DEmRNAs. Finally, we constructed a potentially circRNA-miRNA-mRNA network related to IDD by predicting in the database (ENCORI, TargetScan, miRecords, miRmap, and circBank).ResultsWe identified 31 common DEmRNAs by Venn analysis, of which MMP2 was regarded as the key hub genes. Simultaneously, miR-423-5p and miR-185-5p were predicted as the upstream molecules of MMP2. Furthermore, a total of six DEcircRNAs were predicted as the upstream circRNAs of miR-423-5p and miR-185-5p. Then, a potential circRNA-miRNA-mRNA network related to IDD was constructed by bioinformatics analysis.ConclusionA comprehensive ceRNA regulatory network was constructed, which was found to be significant in IDD progression.
| S-EPMC8357516 | biostudies-literature
Project description:The aim of the present study was to identify genes associated with and the underlying mechanisms in nasopharyngeal carcinoma (NPC) using microarray data. GSE12452 and GSE34573 gene expression profiles were obtained from the Gene Expression Omnibus (GEO) database. GEO2R was utilized to obtain differentially expressed genes (DEGs). In addition, the Database for Annotation, Visualization and Integrated Discovery was used to perform pathway enrichment analyses for DEGs using the Gene Ontology (GO) annotation along with the Kyoto Encyclopedia of Genes and Genomes (KEGG). Furthermore, Cytoscape was used to perform module analysis of the protein-protein interaction (PPI) network and pathways of the hub genes were studied. A total of 298 genes were ascertained as DEGs in the two datasets. To functionally categorize these DEGs, we obtained 82 supplemented GO terms along with 7 KEGG pathways. Subsequently, a PPI network consisting of 10 hub genes with high degrees of interaction was constructed. These hub genes included cyclin-dependent kinase (CDK) 1, structural maintenance of chromosomes (SMC) 4, kinetochore-associated (KNTC) 1, kinesin family member (KIF) 23, aurora kinase A (AURKA), ATAD (ATPase family AAA domain containing) 2, NDC80 kinetochore complex component, enhancer of zeste 2 polycomb repressive complex 2 subunit, BUB1 mitotic checkpoint serine/threonine kinase and protein regulator of cytokinesis 1. CDK1, SMC4, KNTC1, KIF23, AURKA and ATAD2 presented with high areas under the curve in receiver operator curves, suggesting that these genes may be diagnostic markers for nasopharyngeal carcinoma. In conclusion, it was proposed that CDK1, SMC4, KNTC1, KIF23, AURKA and ATAD2 may be involved in the tumorigenesis of NPC. Furthermore, they may be utilized as molecular biomarkers in early diagnosis of NPC.
| S-EPMC6732963 | biostudies-literature