GLRX5 as a Differential Diagnosis Marker for Kashin-Beck Disease and Osteoarthritis and Its Clinical Significance

Abstract Background: Kashin-Beck disease (KBD) is currently an endemic form of osteoarthritis. In this study, we explored novel KBD diagnostic biomarkers.Methods: The GSE59446 dataset was used to conduct Weighted Gene Co-expression Network Analysis (WGCNA) and differentially expressed genes (DEGs) analysis with peripheral blood samples of 100 healthy individuals and 100 KBD patients. As part of the gene ontology pathway enrichment analysis, genes related to SONFH and DEGs were selected from the extraction module. Then, central DEGs were selected for LASSO analysis, and, based on SVM-RFE and DEG results, overlapping genes were identified as key KBD genes. Next, we analyzed the correlations between the selected genes and age, gender, and other factors to eliminate their influences on gene expression. Finally, we evaluated the diagnostic value of key KBD genes using case information collected by us.Results: Seven gene co-expression modules were created using WGCNA. The turquoise module was identified as a KBD key module since it showed the highest correlation to KBD. The functional enrichment analysis revealed that the genes associated with this key module were mainly involved in mitochondrial reactions, protein heterooligomerization, and negatively regulating cysteine-type endopeptidase-dependent apoptotic processes. Additionally, 12 key genes were identified using the LASSO analysis, 5 major genes using SVM-RFE analysis, and 36 DEGs were screened through the "limma" R package. The GLRX5 gene - pivotal in DEGs, LASSO, and SVM-RFE - was further aggregated as the key KBD gene. Correlation analyses confirmed the GLRX5 diagnostic value for KBD and that it was not related to age, gender, and other factors. Finally, data from our patients demonstrated that GLRX5 can be a KBD diagnostic biomarker.Conclusions: We demonstrated that the target gene GLRX5 can be a KBD non-invasive diagnosis biomarker.

Download Full-text

HEMGN and SLC2A1 might be potential diagnostic biomarkers of steroid-induced osteonecrosis of femoral head: study based on WGCNA and DEGs screening

BMC Musculoskeletal Disorders ◽

10.1186/s12891-021-03958-7 ◽

2021 ◽

Vol 22 (1) ◽

Author(s):

Zhixin Wu ◽

Yinxian Wen ◽

Guanlan Fan ◽

Hangyuan He ◽

Siqi Zhou ◽

...

Keyword(s):

Femoral Head ◽

Roc Curve ◽

Enrichment Analysis ◽

Diagnostic Value ◽

Functional Enrichment ◽

Pathway Enrichment Analysis ◽

Overlapping Genes ◽

Hub Genes ◽

Diagnostic Biomarkers ◽

Key Genes

Abstract Background Steroid-induced osteonecrosis of the femoral head (SONFH) is a chronic and crippling bone disease. This study aims to reveal novel diagnostic biomarkers of SONFH. Methods The GSE123568 dataset based on peripheral blood samples from 10 healthy individuals and 30 SONFH patients was used for weighted gene co-expression network analysis (WGCNA) and differentially expressed genes (DEGs) screening. The genes in the module related to SONFH and the DEGs were extracted for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Genes with |gene significance| > 0.7 and |module membership| > 0.8 were selected as hub genes in modules. The DEGs with the degree of connectivity ≥5 were chosen as hub genes in DEGs. Subsequently, the overlapping genes of hub genes in modules and hub genes in DEGs were selected as key genes for SONFH. And then, the key genes were verified in another dataset, and the diagnostic value of key genes was evaluated by receiver operating characteristic (ROC) curve. Results Nine gene co-expression modules were constructed via WGCNA. The brown module with 1258 genes was most significantly correlated with SONFH and was identified as the key module for SONFH. The results of functional enrichment analysis showed that the genes in the key module were mainly enriched in the inflammatory response, apoptotic process and osteoclast differentiation. A total of 91 genes were identified as hub genes in the key module. Besides, 145 DEGs were identified by DEGs screening and 26 genes were identified as hub genes of DEGs. Overlapping genes of hub genes in the key module and hub genes in DEGs, including RHAG, RNF14, HEMGN, and SLC2A1, were further selected as key genes for SONFH. The diagnostic value of these key genes for SONFH was confirmed by ROC curve. The validation results of these key genes in GSE26316 dataset showed that only HEMGN and SLC2A1 were downregulated in the SONFH group, suggesting that they were more likely to be diagnostic biomarkers of SOFNH than RHAG and RNF14. Conclusions Our study identified that two key genes, HEMGN and SLC2A1, might be potential diagnostic biomarkers of SONFH.

Download Full-text

gprofiler2 -- an R package for gene list functional enrichment analysis and namespace conversion toolset g:Profiler

F1000Research ◽

10.12688/f1000research.24956.2 ◽

2020 ◽

Vol 9 ◽

pp. 709 ◽

Cited By ~ 1

Author(s):

Liis Kolberg ◽

Uku Raudvere ◽

Ivan Kuzmin ◽

Jaak Vilo ◽

Hedi Peterson

Keyword(s):

Gene List ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Automated Analysis ◽

R Package ◽

Biological Data ◽

Functional Enrichment ◽

Link Type ◽

Functional Profiling ◽

Rest Api

g:Profiler (https://biit.cs.ut.ee/gprofiler) is a widely used gene list functional profiling and namespace conversion toolset that has been contributing to reproducible biological data analysis already since 2007. Here we introduce the accompanying R package, gprofiler2, developed to facilitate programmatic access to g:Profiler computations and databases via REST API. The gprofiler2 package provides an easy-to-use functionality that enables researchers to incorporate functional enrichment analysis into automated analysis pipelines written in R. The package also implements interactive visualisation methods to help to interpret the enrichment results and to illustrate them for publications. In addition, gprofiler2 gives access to the versatile gene/protein identifier conversion functionality in g:Profiler enabling to map between hundreds of different identifier types or orthologous species. The gprofiler2 package is freely available at the CRAN repository.

Download Full-text

Identification of novel genes and pathways in carotid atheroma using integrated bioinformatic methods

Scientific Reports ◽

10.1038/srep18764 ◽

2016 ◽

Vol 6 (1) ◽

Cited By ~ 10

Author(s):

Wenqing Nai ◽

Diane Threapleton ◽

Jingbo Lu ◽

Kewei Zhang ◽

Hongyuan Wu ◽

...

Keyword(s):

Signaling Pathway ◽

Transforming Growth Factor ◽

Interaction Network ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Transforming Growth Factor Β ◽

Functional Enrichment ◽

Differentially Expressed ◽

Pathway Enrichment Analysis ◽

Atheromatous Plaques

Abstract Atherosclerosis is the primary cause of cardiovascular events and its molecular mechanism urgently needs to be clarified. In our study, atheromatous plaques (ATH) and macroscopically intact tissue (MIT) sampled from 32 patients were compared and an integrated series of bioinformatic microarray analyses were used to identify altered genes and pathways. Our work showed 816 genes were differentially expressed between ATH and MIT, including 443 that were up-regulated and 373 that were down-regulated in ATH tissues. GO functional-enrichment analysis for differentially expressed genes (DEGs) indicated that genes related to the “immune response” and “muscle contraction” were altered in ATHs. KEGG pathway-enrichment analysis showed that up-regulated DEGs were significantly enriched in the “FcεRI-mediated signaling pathway”, while down-regulated genes were significantly enriched in the “transforming growth factor-β signaling pathway”. Protein-protein interaction network and module analysis demonstrated that VAV1, SYK, LYN and PTPN6 may play critical roles in the network. Additionally, similar observations were seen in a validation study where SYK, LYN and PTPN6 were markedly elevated in ATH. All in all, identification of these genes and pathways not only provides new insights into the pathogenesis of atherosclerosis, but may also aid in the development of prognostic and therapeutic biomarkers for advanced atheroma.

Download Full-text

Identification of Latent Diagnostic Biomarkers and Biological Pathways in Dermatomyositis Based on WGCNA

Journal of Oncology ◽

10.1155/2021/1920111 ◽

2021 ◽

Vol 2021 ◽

pp. 1-12

Author(s):

Shaxi Ouyang ◽

Yifang Liu ◽

Changjuan Xiao ◽

Qinghua Zeng ◽

Xun Luo ◽

...

Keyword(s):

Cell Lines ◽

Differentially Expressed Gene ◽

Interaction Network ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Receptor Interaction ◽

Interferon Α ◽

Pathway Enrichment Analysis ◽

Oligoadenylate Synthetase

Introduction. Dermatomyositis (DM) is a chronic autoimmune disease of predominantly lymphocytic infiltration mainly involving the transverse muscle. Its pathogenesis is remaining unknown. This research is designed to probe the latent pathogenesis of dermatomyositis, identify potential biomarkers, and reveal the pathogenesis of dermatomyositis through information biology analysis of gene chips. Methods. In this study, we utilised the GSE14287 and GSE11971 datasets rooted in the Gene Expression Omnibus (GEO) databank, which included a total of 62 DM samples and 9 normal samples. The datasets were combined, and the differentially expressed gene sets were subjected to weighted gene coexpression network analysis, and the hub gene was screened using a protein interaction network from genes in modules highly correlated with dermatomyositis progression. Results. A total of 3 key genes—myxovirus resistance-2 (MX2), oligoadenylate synthetase 1 (OAS1), and oligoadenylate synthetase 2 (OAS2)—were identified in combination with cell line samples, and the expressions of the 3 genes were verified separately. The results showed that MX2, OAS1, and OAS2 were highly expressed in LPS-treated cell lines compared to normal cell lines. The results of pathway enrichment analysis of the genes indicated that all 3 genes were enriched in the cytosolic DNA signalling and cytokine and cytokine receptor interaction signalling pathways; the results of functional enrichment analysis showed that all 3 were enriched in interferon-α response and interferon-γ response functions. Conclusions. This is important for the study of the pathogenesis and objective treatment of dermatomyositis and provides important reference information for the targeted therapy of dermatomyositis.

Download Full-text

Co-expression network analysis identifies specific hub genes in association with developmental neuronal remodeling in Drosophila melanogaster

10.7287/peerj.preprints.27586v1 ◽

2019 ◽

Author(s):

Yunze Liu ◽

Xiaojie Sun ◽

Aijun Qu

Keyword(s):

Drosophila Melanogaster ◽

Interaction Network ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Pathway Enrichment Analysis ◽

Hub Genes ◽

Protein Protein Interaction ◽

Neuronal Remodeling ◽

Gene Modules

As an evolutionarily conserved mechanism, developmental neuronal remodeling is needed for the proper wiring of the nervous system and is critical for understanding the neurodevelopment mechanisms. Previous studies have shown that during metamorphosis lots of Drosophila melanogaster mushroom body neurons experience stereotypic remodeling. However, the related regulators and downstream executors of pathways are yet unclear, especially studies of transcriptional gene co-expression analysis of nervous systems remain insufficient. In this study, we develop a weighted gene co-expression network (WGCNA) to classify gene modules associated with neuronal remodeling. Moreover, functional and pathway enrichment analysis with protein-protein network construction is applied to detect high informative hub genes in the targeted gene module. Thus, we select a total of five hub genes that play critical roles in neuronal remodeling and identify them with functional enrichment analysis and protein-protein interaction network. Overall, this study provides insight into the underlying molecular mechanism of developmental neuronal remodeling in Drosophila melanogaster.

Download Full-text

Co-expression network analysis identifies specific hub genes in association with developmental neuronal remodeling in Drosophila melanogaster

10.7287/peerj.preprints.27586 ◽

2019 ◽

Author(s):

Yunze Liu ◽

Xiaojie Sun ◽

Aijun Qu

Keyword(s):

Drosophila Melanogaster ◽

Interaction Network ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Pathway Enrichment Analysis ◽

Hub Genes ◽

Protein Protein Interaction ◽

Neuronal Remodeling ◽

Gene Modules

Download Full-text

Multiomic Analysis of Methylation and Transcriptome Reveals a Novel Signature in Esophageal Cancer

Dose-Response ◽

10.1177/1559325820942075 ◽

2020 ◽

Vol 18 (3) ◽

pp. 155932582094207

Author(s):

Yi-qi Jin ◽

Dong-liu Miao

Keyword(s):

Esophageal Cancer ◽

Messenger Rna ◽

Cox Model ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

R Package ◽

The Cancer Genome Atlas ◽

Functional Enrichment ◽

Cox Analysis ◽

Clinicopathological Variables

Background: Epigenetic alterations have been shown to lead to human carcinogenesis. The aim of this study was to perform an integrative analysis to develop an epigenetic signature to predict overall survival (OS) of esophageal cancer. Methods: DNA methylation and messenger RNA expression data of esophageal cancer samples were downloaded from The Cancer Genome Atlas database and were incorporated and analyzed using an R package MethylMix. Functional enrichment analysis of the methylation-related differentially expressed genes (DEGs) was performed. Epigenetic signature and nomogram associated with the OS of esophageal cancer were established by the multivariate Cox model. Results: A total of 71 methylation-related DEGs were identified. Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that these genes were involved in the biological process related to the initiation and progression of esophageal cancer. Two-gene (FAM24B and FAM200A) risk signature for OS was developed by multivariate Cox analysis, of which had high accuracy. The signature is independent of clinicopathological variables and indicated better predictive power than other clinicopathological variables. Moreover, we developed a novel prognostic nomogram based on risk score and 3 clinicopathological factors. Conclusions: Our study indicated possible methylation-related DEGs and established an epigenetic signature, which may provide novel insights for understanding the pathogenesis of esophageal cancer.

Download Full-text

Systems Level Analysis and Identification of Pathways and Key Genes Associated with Delirium

Genes ◽

10.3390/genes11101225 ◽

2020 ◽

Vol 11 (10) ◽

pp. 1225 ◽

Cited By ~ 1

Author(s):

Yukiko Takahashi ◽

Tomoyoshi Terada ◽

Yoshinori Muto

Keyword(s):

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

High Tendency ◽

Diagnostic Biomarkers ◽

Protein Protein Interaction ◽

Gene Sets ◽

Network Modules ◽

Complex Detection ◽

Human Accelerated Regions

Delirium is a complex pathophysiological process, and multiple contributing mechanisms have been identified. However, it is largely unclear how the genes associated with delirium contribute and which of them play key roles. In this study, the genes associated with delirium were retrieved from the Comparative Toxicogenomics Database (CTD) and integrated through a protein–protein interaction (PPI) network. Delirium-associated genes formed a highly interconnected PPI subnetwork, indicating a high tendency to interact and agglomerate. Using the Molecular Complex Detection (MCODE) algorithm, we identified the top two delirium-relevant network modules, M1 and M5, that have the most significant enrichments for the delirium-related gene sets. Functional enrichment analysis showed that genes related to neurotransmitter receptor activity were enriched in both modules. Moreover, analyses with genes located in human accelerated regions (HARs) provided evidence that HAR-Brain genes were overrepresented in the delirium-relevant network modules. We found that four of the HAR-Brain genes, namely APP, PLCB1, NPY, and HTR2A, in the M1 module were highly connected and appeared to exhibit hub properties, which might play vital roles in delirium development. Further understanding of the function of the identified modules and member genes could help to identify therapeutic intervention targets and diagnostic biomarkers for delirium.

Download Full-text

LncRNA SNHG8 is identified as a key regulator of acute myocardial infarction by RNA-seq analysis

Lipids in Health and Disease ◽

10.1186/s12944-019-1142-0 ◽

2019 ◽

Vol 18 (1) ◽

Cited By ~ 1

Author(s):

Liu-An Zhuo ◽

Yi-Tao Wen ◽

Yong Wang ◽

Zhi-Fang Liang ◽

Gang Wu ◽

...

Keyword(s):

Myocardial Infarction ◽

Acute Myocardial Infarction ◽

Roc Curve ◽

Regulatory Networks ◽

Pairwise Comparison ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Diagnostic Value ◽

Functional Enrichment ◽

Rna Seq

Abstract Background Long noncoding RNAs (lncRNAs) are involved in numerous physiological functions. However, their mechanisms in acute myocardial infarction (AMI) are not well understood. Methods We performed an RNA-seq analysis to explore the molecular mechanism of AMI by constructing a lncRNA-miRNA-mRNA axis based on the ceRNA hypothesis. The target microRNA data were used to design a global AMI triple network. Thereafter, a functional enrichment analysis and clustering topological analyses were conducted by using the triple network. The expression of lncRNA SNHG8, SOCS3 and ICAM1 was measured by qRT-PCR. The prognostic values of lncRNA SNHG8, SOCS3 and ICAM1 were evaluated using a receiver operating characteristic (ROC) curve. Results An AMI lncRNA-miRNA-mRNA network was constructed that included two mRNAs, one miRNA and one lncRNA. After RT-PCR validation of lncRNA SNHG8, SOCS3 and ICAM1 between the AMI and normal samples, only lncRNA SNHG8 had significant diagnostic value for further analysis. The ROC curve showed that SNHG8 presented an AUC of 0.850, while the AUC of SOCS3 was 0.633 and that of ICAM1 was 0.594. After a pairwise comparison, we found that SNHG8 was statistically significant (PSNHG8-ICAM1 = 0.002; PSNHG8-SOCS3 = 0.031). The results of a functional enrichment analysis of the interacting genes and microRNAs showed that the shared lncRNA SNHG8 may be a new factor in AMI. Conclusions Our investigation of the lncRNA-miRNA-mRNA regulatory networks in AMI revealed a novel lncRNA, lncRNA SNHG8, as a risk factor for AMI and expanded our understanding of the mechanisms involved in the pathogenesis of AMI.

Download Full-text

Identifying Synergistic Mechanisms of Multiple Ingredients in Shuangbai Tablets against Proteinuria by Virtual Screening and a Network Pharmacology Approach

Evidence-based Complementary and Alternative Medicine ◽

10.1155/2020/1027271 ◽

2020 ◽

Vol 2020 ◽

pp. 1-15

Author(s):

Xiaoqin Ma ◽

Meixiang Yu ◽

Chenxia Hao ◽

Wanhua Yang

Keyword(s):

Molecular Docking ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Network Pharmacology ◽

Pathway Enrichment Analysis ◽

Active Ingredients ◽

Human Phenotype ◽

Herbal Mixture ◽

Relationship Of

Shuangbai Tablets (SBT), a traditional herbal mixture, has shown substantial clinical efficacy. However, a systematic mechanism of its active ingredients and pharmacological mechanisms of action against proteinuria continues being lacking. A network pharmacology approach was effectual in discovering the relationship of multiple ingredients and targets of the herbal mixture. This study aimed to identify key targets, major active ingredients, and pathways of SBT against proteinuria by network pharmacology approach combined with thin layer chromatography (TLC). Human phenotype (HP) disease analysis, gene ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, and molecular docking were used in this study. To this end, a total of 48 candidate targets of 118 active ingredients of SBT were identified. Network analysis showed PTGS2, ESR1, and NOS2 to be the three key targets, and beta-sitosterol, quercetin, and berberine were the three major active ingredients; among them one of the major active ingredients, quercetin, was discriminated by TLC. These results of the functional enrichment analysis indicated that the most relevant disease including these 48 candidate proteins is proteinuria, SBT treated proteinuria by sympathetically regulating multiple biological pathways, such as the HIF-1, RAS, AGE-RAGE, and VEGF signaling pathways. Additionally, molecular docking validation suggested that major active ingredients of SBT were capable of binding to HIF-1A and VEGFA of the main pathways. Consequently, key targets, major active ingredients, and pathways based on data analysis of SBT against proteinuria were systematically identified confirming its utility and providing a new drug against proteinuria.

Download Full-text