Dataset Information

Weighted gene co-expression network analysis reveals specific modules and hub genes related to immune infiltration of osteoarthritis.

ABSTRACT:

Background

The incidence of osteoarthritis (OA), a chronic degenerative disease, is increasing every year. There is no effective clinical treatment for OA and the pathological mechanism remains unclear. Early diagnosis is an effective strategy to control the progress of OA. In this study, we aimed to identify potential early diagnostic biomarkers.

Methods

We downloaded the gene expression profile dataset, GSE51588 and GSE55235, from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) public database. The differentially expressed genes (DEGs) were screened out using the "limma" R package. Weighted gene co-expression network analysis (WGCNA) was utilized to build the co-expression network between the normal and OA samples. A Venn diagram was constructed to detect the hub genes. Potential molecular mechanisms and signaling pathways were enriched by gene set variation analysis (GSVA). Single sample gene set enrichment analysis (ssGSEA) was used to identify the immune infiltration of OA.

Results

We screened out three hub genes based on WGCNA and DEGs in this study. GSVA results showed that nuclear factor interleukin-3 (NFIL3) was related to tumor necrosis factor alpha (TNF-α) signaling via nuclear factor kappa-B (NF-κB), the reactive oxygen species pathway, and myelocytomatosis (MYC) targets v2. Highly-expressed ADM (adrenomedullin) pathways included TNF-α signaling via NF-κB, the reactive oxygen species pathway, and ultraviolet (UV) response up. OGN (osteoglycin)-enriched pathways included epithelial mesenchymal transition, coagulation, and peroxisome.

Conclusions

We identified three hub genes (NFIL3, ADM, and OGN) that were correlated to the development and progression of OA, which may provide new biomarkers for early diagnosis.

SUBMITTER: Cao J

PROVIDER: S-EPMC8576690 | biostudies-literature | 2021 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Weighted gene co-expression network analysis reveals specific modules and hub genes related to immune infiltration of osteoarthritis.

Cao Jiangang J Ding Han H Shang Jun J Ma Lei L Wang Qi Q Feng Shiqing S

Annals of translational medicine 20211001 20

<h4>Background</h4>The incidence of osteoarthritis (OA), a chronic degenerative disease, is increasing every year. There is no effective clinical treatment for OA and the pathological mechanism remains unclear. Early diagnosis is an effective strategy to control the progress of OA. In this study, we aimed to identify potential early diagnostic biomarkers.<h4>Methods</h4>We downloaded the gene expression profile dataset, GSE51588 and GSE55235, from the National Center for Biotechnology Informatio ...[more]

PMID: 34790731

Similar Datasets

Project description:BackgroundThe therapeutic management of obesity is challenging, hence further elucidating the underlying mechanisms of obesity development and identifying new diagnostic biomarkers and therapeutic targets are urgent and necessary. Here, we performed differential gene expression analysis and weighted gene co-expression network analysis (WGCNA) to identify significant genes and specific modules related to BMI based on gene expression profile data of 7 discordant monozygotic twins.ResultsIn the differential gene expression analysis, it appeared that 32 differentially expressed genes (DEGs) were with a trend of up-regulation in twins with higher BMI when compared to their siblings. Categories of positive regulation of nitric-oxide synthase biosynthetic process, positive regulation of NF-kappa B import into nucleus, and peroxidase activity were significantly enriched within GO database and NF-kappa B signaling pathway within KEGG database. DEGs of NAMPT, TLR9, PTGS2, HBD, and PCSK1N might be associated with obesity. In the WGCNA, among the total 20 distinct co-expression modules identified, coral1 module (68 genes) had the strongest positive correlation with BMI (r = 0.56, P = 0.04) and disease status (r = 0.56, P = 0.04). Categories of positive regulation of phospholipase activity, high-density lipoprotein particle clearance, chylomicron remnant clearance, reverse cholesterol transport, intermediate-density lipoprotein particle, chylomicron, low-density lipoprotein particle, very-low-density lipoprotein particle, voltage-gated potassium channel complex, cholesterol transporter activity, and neuropeptide hormone activity were significantly enriched within GO database for this module. And alcoholism and cell adhesion molecules pathways were significantly enriched within KEGG database. Several hub genes, such as GAL, ASB9, NPPB, TBX2, IL17C, APOE, ABCG4, and APOC2 were also identified. The module eigengene of saddlebrown module (212 genes) was also significantly correlated with BMI (r = 0.56, P = 0.04), and hub genes of KCNN1 and AQP10 were differentially expressed.ConclusionWe identified significant genes and specific modules potentially related to BMI based on the gene expression profile data of monozygotic twins. The findings may help further elucidate the underlying mechanisms of obesity development and provide novel insights to research potential gene biomarkers and signaling pathways for obesity treatment. Further analysis and validation of the findings reported here are important and necessary when more sample size is acquired.

Project description:Purpose:Despite advances in characterizing the neurobiology of emotional disorders, there is still a significant lack of scientific understanding of the pathophysiological mechanisms governing major depressive disorder (MDD). This study attempted to elucidate the molecular circuitry of MDD and to identify more potential genes associated with the pathogenesis of the disease. Patients and Methods:Microarray data from the GSE98793 dataset were downloaded from the NCBI Gene Expression Omnibus (GEO) database, including 128 patients with MDD and 64 healthy controls. Weighted gene coexpression network analysis (WGCNA) was performed to find modules of differentially expressed genes (DEGs) with high correlations followed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses to obtain further biological insight into the top three key modules. The protein-protein interaction (PPI) network, the modules from the PPI network, and the gene annotation enrichment of modules were analyzed, as well. Results:We filtered 3276 genes that were considered significant DEGs for further WGCNA analysis. By performing WGCNA, we found that the turquoise, blue and brown functional modules were all strongly correlated with MDD development, including immune response, neutrophil degranulation, ribosome biogenesis, T cell activation, glycosaminoglycan biosynthetic process, and protein serine/threonine kinase activator activity. Hub genes were identified in the key functional modules that might have a role in the progression of MDD. Functional annotation showed that these modules primarily enriched such KEGG pathways as the TNF signaling pathway, T cell receptor signaling pathway, primary immunodeficiency, Th1, Th2 and Th17 cell differentiation, autophagy and RNA degradation and oxidative phosphorylation. These results suggest that these genes are closely related to autophagy and cellular immune function. Conclusion:The results of this study may help to elucidate the pathophysiology of MDD development at the molecular level and explore the potential molecular mechanisms for new interventional strategies.

Project description:BackgroundWhile large-scale genomic analyses symbolize a precious attempt to decipher the molecular foundation of uterine leiomyosarcoma (ULMS), bioinformatics results associated with the occurrence of ULMS based totally on WGCNA and CIBERSORT have not yet been reported. This study aimed to screen the hub genes and the immune cell infiltration pattern in ULMS by bioinformatics methods.MethodsFirstly, the GSE67463 dataset, including 25 ULMS tissues and 29 normal myometrium (NL) tissues, was downloaded from the public database. The differentially expressed genes (DEGs) were screened by the 'limma' package and hub modules were identified by weighted gene co-expression network analysis (WGCNA). Subsequently, gene function annotations were performed to investigate the biological role of the genes from the intersection of two groups (hub module and DEGs). The above genes were calculated in the protein-protein interaction (PPI) network to select the hub genes further. The hub genes were validated using external data (GSE764 and GSE68295). In addition, the differential immune cell infiltration between UL and ULMS tissues was investigated using the CIBERSORT algorithm. Finally, we used western blot to preliminarily detect the hub genes in cell lines.ResultsWGCNA analysis revealed a green-yellow module possessed the highest correlation with ULMS, including 1063 genes. A total of 172 DEGs were selected by thresholds set in the 'limma' package. The above two groups of genes were intersected to obtain 72 genes for functional annotation analysis. Interestingly, it indicated that 72 genes were mainly involved in immune processes and the Neddylation pathway. We found a higher infiltration of five types of cells (memory B cells, M0-type macrophages, mast cells activated, M1-type macrophages, and T cells follicular helper) in ULMS tissues than NL tissues, while the infiltration of two types of cells (NK cells activated and mast cells resting) was lower than in NL tissues. In addition, a total of five genes (KDR, CCL21, SELP, DPT, and DCN) were identified as the hub genes. Internal and external validation demonstrated that the five genes were over-expressed in NL tissues compared with USML tissues. Finally, the correlation analysis results indicate that NK cells activated and mast cells activated positively correlated with the hub genes. However, M1-type macrophages had a negative correlation with the hub genes. Moreover, only the DCN may be associated with the Neddylation pathway.ConclusionA series of evidence confirm that the five hub genes and the infiltration of seven types of immune cells are related to USML occurrence. These hub genes may affect the occurrence of USML through immune-related and Neddylation pathways, providing molecular evidence for the treatment of USML in the future.

Dataset Information

Weighted gene co-expression network analysis reveals specific modules and hub genes related to immune infiltration of osteoarthritis.

Background

Methods

Results

Conclusions

Publications

Weighted gene co-expression network analysis reveals specific modules and hub genes related to immune infiltration of osteoarthritis.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets