Network Analyses of Differentially Expressed Genes in Osteoarthritis to Identify Hub Genes.
Ontology highlight
ABSTRACT:
SUBMITTER: Li Z
PROVIDER: S-EPMC6500622 | biostudies-literature | 2019
REPOSITORIES: biostudies-literature
Publications
Network Analyses of Differentially Expressed Genes in Osteoarthritis to Identify Hub Genes.
Li Zhaoyan Z Zhong Lei L Du Zhenwu Z Chen Gaoyang G Shang Jing J Yang Qiwei Q Wang Qingyu Q Song Yang Y Zhang Guizhen G
BioMed research international 20190416
<h4>Background</h4>Osteoarthritis (OA) is the most common degenerative disease in orthopedics. However, the cause and underlying molecular mechanisms are not clear. This study aims to identify the hub genes and pathways involved in the occurrence of osteoarthritis.<h4>Methods</h4>The raw data of GSE89408 were downloaded from the Gene Expression Omnibus (GEO) database, and the differentially expressed genes (DEGs) were identified by R software. The DAVID database was used for pathway and gene ont ...[more]
Similar Datasets
Project description:Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer. Identifying characteristic genes of PTC are of great importance to reveal its potential genetic mechanisms. In this paper, we proposed a framework, as well as a measure named Normalized Centrality Measure (NCM), to identify characteristic genes of PTC. The framework consisted of four steps. First, both up-regulated genes and down-regulated genes, collectively called differentially expressed genes (DEGs), were screened and integrated together from four datasets, that is, GSE3467, GSE3678, GSE33630, and GSE58545; second, an interaction network of DEGs was constructed, where each node represented a gene and each edge represented an interaction between linking nodes; third, both traditional measures and the NCM measure were used to analyze the topological properties of each node in the network. Compared with traditional measures, more genes related to PTC were identified by the NCM measure; fourth, by mining the high-density subgraphs of this network and performing Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, several meaningful results were captured, most of which were demonstrated to be associated with PTC. The experimental results proved that this network framework and the NCM measure are useful for identifying more characteristic genes of PTC.
| S-EPMC6356810 | biostudies-literature
Project description:ObjectiveMultiple sclerosis (MS) is an autoimmune neuroinflammatory disease of the nervous system. However, the precise molecular mechanisms underlying MS have yet to be fully elucidated. In this study, our aim was to provide novel insight into the pathogenesis of MS and provide a resource for identifying new biomarkers and therapeutics for MS.MethodsIn this study, we analyzed the gene expression profiles (GSE21942) and miRNA expression profiles (GSE61741) of MS patient samples that were downloaded from the GEO database and identified differentially expressed mRNAs and miRNAs (DEmRNAs, DEmiRNAs). Next, we constructed a protein-protein interaction (PPI) network and a MS-specific ceRNA network (MCEN) by integrating expression profiles, interaction pairs of mRNA-miRNAs and lncRNA-miRNAs. Then, according to the modular structure of the PPI network, we identified hub DEmRNAs and generated a ceRNA subnetwork so that we could analyze the key lncRNAs that were associated with MS.ResultsWe first identified 4 modules by constructing a PPI network using DEmRNAs. Functional enrichment analysis showed these modules were enriched in immune-related pathways. Then, we constructed the MCEN and the hub gene-associated ceRNA subnetwork using a comprehensive computational approach. We identified three key lncRNAs (LINC00649, TP73-AS1 and MALAT1) and further identified key lncRNA-mediated ceRNAs within the subnetwork. Finally, by analyzing LINC00649-miR-1275-CD20, we identified 6 drugs that may represent novel drugs for MS.ConclusionCollectively, our results provide novel insight for the discovery of biomarkers and therapeutics for MS and provide a suitable foundation from which to design future investigations of the pathogenic mechanisms associated with MS.
| S-EPMC8902536 | biostudies-literature
Project description:IntroductionOsteoarthritis (OA) refers to a commonly seen degenerative joint disorder and a major global public health burden. According to the existing literature, osteoarthritis is related to epigenetic changes, which are important for diagnosing and treating the disease early. Through early targeted treatment, costly treatments and poor prognosis caused by advanced osteoarthritis can be avoided.MethodsThis study combined gene differential expression analysis and weighted gene co-expression network analysis (WGCNA) of the transcriptome with epigenome microarray data to discover the hub gene of OA. We obtained 2 microarray datasets (GSE114007, GSE73626) in Gene Expression Omnibus (GEO). The R software was utilized for identifying differentially expressed genes (DEGs) and differentially methylated genes (DMGs). By using WGCNA to analyze the relationships between modules and phenotypes, it was discovered that the blue module (MEBlue) has the strongest phenotypic connection with OA (cor = 0.92, p = 4e-16). The hub genes for OA, also known as the hub methylated differentially expressed genes, were identified by matching the MEblue module to differentially methylated differentially expressed genes. Furthermore, this study used Gene set variation analysis (GSVA) to identify specific signal pathways associated with hub genes. qRT-PCR and western blotting assays were used to confirm the expression levels of the hub genes in OA patients and healthy controls.ResultsThree hub genes were discovered: HTRA1, P2RY6, and RCAN1. GSVA analysis showed that high HTRA1 expression was mainly enriched in epithelial-mesenchymal transition and apical junction; high expression of P2RY6 was mainly enriched in the peroxisome, coagulation, and epithelial-mesenchymal transition; and high expression of RCAN1 was mainly enriched in epithelial-mesenchymal-transition, TGF-β-signaling, and glycolysis. The results of the RT-qPCR and WB assay were consistent with the findings.DiscussionThe three genes tested may cause articular cartilage degeneration by inducing chondrocyte hypertrophy, regulating extracellular matrix accumulation, and improving macrophage pro-inflammatory response, resulting in the onset and progression of osteoarthritis. They can provide new ideas for targeted treatment of osteoarthritis.
| S-EPMC10351907 | biostudies-literature
Project description:The present study was performed to explore the underlying molecular mechanisms and screen hub genes of osteoarthritis (OA) via bioinformatics analysis. In total, twenty-five OA synovial tissue samples and 25 normal synovial tissue samples were derived from three datasets, namely, GSE55457, GSE55235, and GSE1919, and were used to identify the differentially expressed genes (DEGs) of OA by R language. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of DEGs were conducted using the Database for Annotation, Visualization, and Integrated Discovery (DAVID). A Venn diagram was built to show the potential hub genes identified in all three datasets. The STRING database was used for constructing the protein-protein interaction (PPI) networks and submodules of DEGs. We identified 507 upregulated and 620 downregulated genes. Upregulated DEGs were significantly involved in immune response, MHC class II receptor activity, and presented in the extracellular region, while downregulated DEGs were mainly enriched in response to organic substances, extracellular region parts, and cadmium ion binding. Results of KEGG analysis indicated that the upregulated DEGs mainly existed in cell adhesion molecules (CAMs), while downregulated DEGs were significantly involved in the MAPK signaling pathway. A total of eighteen intersection genes were identified across the three datasets. These include Nell-1, ATF3, RhoB, STC1, and VEGFA. In addition, 10 hub genes including CXCL12, CXCL8, CCL20, and CCL4 were found in the PPI network and module construction. Identification of DEGs and hub genes associated with OA may be helpful for revealing the molecular mechanisms of OA and further promotes the development of relevant biomarkers and drug targets.
| S-EPMC9207392 | biostudies-literature
Project description:BackgroundInadequate immunity caused by poor immune surveillance leads to tumorigenesis, while excessive immunity due to breakdown of immune tolerance causes autoimmune genesis. Although the function of immunity during the onset of these two processes appears to be distinct, the underlying mechanism is shared. To date, gene expression data for large bodies of clinical samples are available, but the resemblances of tumorigenesis and autoimmune genesis in terms of immune responses remains to be summed up.MethodsConsidering the high disease prevalence, we chose invasive ductal carcinoma (IDC) and systemic lupus erythematosus (SLE) to study the potential commonalities of immune responses. We obtained gene expression data of IDC/SLE patients and normal controls from five IDC databases (GSE29044, GSE21422, GSE22840, GSE15852, and GSE9309) and five SLE databases (GSE154851, GSE99967, GSE61635, GSE50635, and GSE17755). We intended to identify genes differentially expressed in both IDC and SLE by using three bioinformatics tools including GEO2R, the limma R package, and Weighted Gene Co-expression Network Analysis (WGCNA) to perform function enrichment, protein-protein network, and signaling pathway analyses.ResultsThe mRNA levels of signal transducer and activator of transcription 1 (STAT1), 2'-5'-oligoadenylate synthetase 1 (OAS1), 2'-5'-oligoadenylate synthetase like (OASL), and PML nuclear body scaffold (PML) were found to be differentially expressed in both IDC and SLE by using three different bioinformatics tools of GEO2R, the limma R package and WGCNA. From the combined databases in this study, the mRNA levels of STAT1 and OAS1 were increased in IDC while reduced in SLE. And the mRNA levels of OASL and PML were elevated in both IDC and SLE. Based on Kyoto Encyclopedia of Genes and Genomes pathway analysis and QIAGEN Ingenuity Pathway Analysis, both IDC and SLE were correlated with the changes of multiple components involved in the Interferon (IFN)-Janus kinase (JAK)-signal transducer and activator of transcription (STAT) signaling pathway.ConclusionThe expression levels of STAT1 and OAS1 manifest the opposite expression tendency across cancer and autoimmune disease. They are components in the IFN-JAK-STAT signaling pathway related to both tumorigenesis and autoimmune genesis. STAT1 and OAS1-associated IFN-JAK-STAT signaling could explain the commonalities during tumorigenesis and autoimmune genesis and render significant information for more precise treatment from the point of immune homeostasis.
| S-EPMC9921081 | biostudies-literature
Project description:Muscle, a multinucleate syncytium formed by the fusion of mononuclear myoblasts, arises from quiescent progenitors (satellite cells) via activation of muscle-specific transcription factors (MyoD, Myf5, myogenin: MYOG, and MRF4). Subsequent to a decline in Pax7, induction in the expression of MYOG is a hallmark of myoblasts that have entered the differentiation phase following cell cycle withdrawal. It is evident that MYOG function cannot be compensated by any other myogenic regulatory factors (MRFs). Despite a plethora of information available regarding MYOG, the mechanism by which MYOG regulates muscle cell differentiation has not yet been identified. Using an RNA-Seq approach, analysis of MYOG knock-down muscle satellite cells (MSCs) have shown that genes associated with cell cycle and division, DNA replication, and phosphate metabolism are differentially expressed. By constructing an interaction network of differentially expressed genes (DEGs) using GeneMANIA, cadherin-associated protein (CTNNA2) was identified as the main hub gene in the network with highest node degree. Four functional clusters (modules or communities) were identified in the network and the functional enrichment analysis revealed that genes included in these clusters significantly contribute to skeletal muscle development. To confirm this finding, in vitro studies revealed increased expression of CTNNA2 in MSCs on day 12 compared to day 10. Expression of CTNNA2 was decreased in MYOG knock-down cells. However, knocking down CTNNA2, which leads to increased expression of extracellular matrix (ECM) genes (type I collagen α1 and type I collagen α2) along with myostatin (MSTN), was not found significantly affecting the expression of MYOG in C2C12 cells. We therefore propose that MYOG exerts its regulatory effects by acting upstream of CTNNA2, which in turn regulates the differentiation of C2C12 cells via interaction with ECM genes. Taken together, these findings highlight a new mechanism by which MYOG interacts with CTNNA2 in order to promote myoblast differentiation.
| S-EPMC4511796 | biostudies-literature
Project description:BackgroundOsteoarthritis (OA) is the most common chronic disorder of joints; however, the key genes and transcription factors (TFs) associated with OA are still unclear. Through bioinformatics tools, the study aimed to understand the mechanism of genes associated with the development of OA.MethodsFour gene expression profiling datasets were used to identify differentially expressed genes (DEGs) between OA and healthy control samples by a meta-analysis. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed with Multifaceted Analysis Tool for Human Transcriptome (MATHT). Subsequently, a protein-protein interaction (PPI) network was constructed for these DEGs. Significant network modules were identified using ReactomeFIViz, and the pathway of each module was enriched using MATHT. In addition, TFs in the DEGs were identified.ResultsIn total, 690 DEGs were identified between OA and healthy control samples, including 449 upregulated and 241 downregulated DEGs. Additionally, 622 nodes and 2752 interactions constituted the PPI network, including 401 upregulated and 221 downregulated DEGs. Among them, FOS, TWIST1, POU2F1, SMARCA4, and CREBBP were also identified as TFs. RT-PCR results showed that the expression levels of Fos, Twist1, Pou2f1, Smarca4, and Crebbp decreased in mice with OA. In addition, FOS, TWIST1, SMARCA4, and CREBBP were involved in the positive regulation of transcription from the RNA polymerase II promoter.ConclusionsTWIST1, POU2F1, SMARCA4, and CREBBP may play an important role in OA pathology.
| S-EPMC5842536 | biostudies-literature
Project description:Increasing evidence has indicated that modulation of epigenetic mechanisms, especially methylation and long-non-coding RNA (lncRNA) regulation, plays a pivotal role in the process of atherosclerosis; however, few studies focused on revealing the epigenetic-related subgroups during atherosclerotic progression using unsupervised clustering analysis. Hence, we aimed to identify the epigenetics-related differentially expressed genes associated with atherosclerosis subtypes and characterize their clinical utility in atherosclerosis. Eighty samples with expression data (GSE40231) and 49 samples with methylation data (GSE46394) from a large artery plaque were downloaded from the GEO database, and aberrantly methylated-differentially expressed (AMDE) genes were identified based on the relationship between methylation and expression. Furthermore, we conducted weighted correlation network analysis (WGCNA) and co-expression analysis to identify the core AMDE genes strongly involved in atherosclerosis. K-means clustering was used to characterize two subtypes of atherosclerosis in GSE40231, and then 29 samples were recognized as validation dataset (GSE28829). In a blood sample cohort (GSE90074), chi-square test and logistic analysis were performed to explore the clinical implication of the K-means clusters. Furthermore, significance analysis of microarrays and prediction analysis of microarrays (PAM) were applied to identify the signature AMDE genes. Moreover, the classification performance of signature AMDE gene-based classifier from PAM was validated in another blood sample cohort (GSE34822). A total of 1,569 AMDE mRNAs and eight AMDE long non-coding RNAs (lncRNAs) were identified by differential analysis. Through the WGCNA and co-expression analysis, 32 AMDE mRNAs and seven AMDE lncRNAs were identified as the core genes involved in atherosclerosis development. Functional analysis revealed that AMDE genes were strongly related to inflammation and axon guidance. In the clinical analysis, the atherosclerotic subtypes were associated with the severity of coronary artery disease and risk of adverse events. Eight genes, including PARP15, SERGEF, PDGFD, MRPL45, UBR1, STAU1, WIZ, and LSM4, were selected as the signature AMDE genes that most significantly differentiated between atherosclerotic subtypes. Ultimately, the area under the curve of signature AMDE gene-based classifier for atherosclerotic subtypes was 0.858 and 0.812 in GSE90074 and GSE34822, respectively. This study identified the AMDE genes (lncRNAs and mRNAs) that could be implemented in clinical clustering to recognize high-risk atherosclerotic patients.
| S-EPMC7767999 | biostudies-literature
Project description:The aim of the present study was to screen the differentially expressed genes (DEGs) associated with glaucoma and investigate the changing patterns of the expression of these genes. The GSE2378 gene microarray data of glaucoma was downloaded from the Gene Expression Omnibus database, which included seven normal samples and eight glaucoma astrocyte samples. Taking into account the corresponding associations between the probe ID and gene symbols, the DEGs were identified prior to and subsequent to the summation of probe level values using the Limma package in R language, followed by Gene Ontology (GO) and pathway enrichment analyses. Interaction networks of the DEGs were constructed using the Biomolecular Interaction Network Database, and cluster analysis of the genes in the networks was performed using ClusterONE. Subsequent to the summation of probe value, a total of 223 genes were identified as DEGs between the normal and glaucoma samples, including 74 downregulated and 149 upregulated genes. In addition, the DEGs were found to be associated with several functions, including response to wounding, extracellular region part and calcium ion binding. The most significantly enriched pathways were complement and coagulation cascades, arrhythmogenic right ventricular cardiomyopathy and extracellular matrix (ECM)‑receptor interaction. Furthermore, interaction networks were constructed of the DEGs prior to and subsequent to the summation of probe values, and HNF4A and CEBPD were identified as hub genes. Additionally, 37 and 31 GO terms were identified to be enriched in the two DEGs of the networks prior to and subsequent to summation, respectively. The results indicated the identified genes associated with ECM as important, and the CEBPD gene was considered to be a critical gene in glaucoma. The findings of the present study offer a potential reference value in further investigations of glaucoma at the gene level.
| S-EPMC4581750 | biostudies-literature
Project description:Introduction: Abnormal status of gene expression plays an important role in tumorigenesis, progression and metastasis of breast cancer. Mechanisms of gene silence or activation were varied. Methylation of genes may contribute to alteration of gene expression. This study aimed to identify differentially expressed hub genes which may be regulated by DNA methylation and evaluate their prognostic value in breast cancer by bioinformatic analysis. Methods: GEO2R was used to obtain expression microarray data from GSE54002, GSE65194 and methylation microarray data from GSE20713, GSE32393. Differentially expressed-aberrantly methylated genes were identified by FunRich. Biological function and pathway enrichment analysis were conducted by DAVID. PPI network was constructed by STRING and hub genes was sorted by Cytoscape. Expression and DNA methylation of hub genes was validated by UALCAN and MethHC. Clinical outcome analysis of hub genes was performed by Kaplan Meier-plotter database for breast cancer. IHC was performed to analyze protein levels of EXO1 and Kaplan-Meier was used for survival analysis. Results: 677 upregulated-hypomethylated and 361 downregulated-hypermethylated genes were obtained from GSE54002, GSE65194, GSE20713 and GSE32393 by GEO2R and FunRich. The most significant biological process, cellular component, molecular function enriched and pathway for upregulated-hypomethylated genes were viral process, cytoplasm, protein binding and cell cycle respectively. For downregulated-hypermethylated genes, the result was peptidyl-tyrosine phosphorylation, plasma membrane, transmembrane receptor protein tyrosine kinase activity and Rap1 signaling pathway (All p< 0.05). 12 hub genes (TOP2A, MAD2L1, FEN1, EPRS, EXO1, MCM4, PTTG1, RRM2, PSMD14, CDKN3, H2AFZ, CCNE2) were sorted from 677 upregulated-hypomethylated genes. 4 hub genes (EGFR, FGF2, BCL2, PIK3R1) were sorted from 361 downregulated-hypermethylated genes. Differential expression of 16 hub genes was validated in UALCAN database (p<0.05). 7 in 12 upregulated-hypomethylated and 2 in 4 downregulated-hypermethylated hub genes were confirmed to be significantly hypomethylated or hypermethylated in breast cancer using MethHC database (p<0.05). Finally, 12 upregulated hub genes (TOP2A, MAD2L1, FEN1, EPRS, EXO1, MCM4, PTTG1, RRM2, PSMD14, CDKN3, H2AFZ, CCNE2) and 3 downregulated genes (FGF2, BCL2, PIK3R1) contributed to significant unfavorable clinical outcome in breast cancer (p<0.05). High expression level of EXO1 protein was significantly associated with poor OS in breast cancer patients (p=0.03). Conclusion: Overexpression of TOP2A, MAD2L1, FEN1, EPRS, EXO1, MCM4, PTTG1, RRM2, PSMD14, CDKN3, H2AFZ, CCNE2 and downregulation of FGF2, BCL2, PIK3R1 might serve as diagnosis and poor prognosis biomarkers in breast cancer by more research validation. EXO1 was identified as an individual unfavorable prognostic factor. Methylation might be one of the major causes leading to abnormal expression of those genes. Functional analysis and pathway enrichment analysis of those genes would provide novel ideas for breast cancer research.
| S-EPMC6856906 | biostudies-literature