Dataset Information

Identification of Hub Genes Associated with Diabetes Mellitus and Tuberculosis Using Bioinformatic Analysis.

ABSTRACT:

Purpose

To investigate the potential pathophysiological association between tuberculosis (TB) and diabetes mellitus (DM) using bioinformatic analyses.

Patients and methods

Gene expression datasets for healthy controls (HCs), TB patients, DM patients, TB+DM patients (GSE114192), and metformin-treated cells (GSE102677) were obtained from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified from pairwise dataset comparisons TB vs HCs and DM vs HCs. DEGs were verified by comparing them to DEGs for TB+DM vs HCs. Enrichment analysis of DEGs common to all three dataset comparisons was conducted using DAVID. The protein-protein interaction (PPI) network was established via STRING and visualised in Cytoscape. Hub genes were identified using the Cytoscape plug-in cytoHubba and then were verified using reverse transcription-quantitative polymerase chain reaction (RT-qPCR) analysis. Targeted miRNA prediction analysis identified metformin treatment-induced gene expression changes in peripheral blood mononuclear cells.

Results

A total of 422 DEGs were common to all three dataset comparisons. Functional enrichment analysis revealed these DEGs were enriched for functional terms of type I interferon signaling pathway, innate immune response, inflammatory response, and infectious diseases. Ten hub genes identified using PPI network analysis were screened for interactions with metformin target gene INS using cytoHubba based on maximal clique centrality (MCC) score. Subsequently, five hub genes were predicted to functionally interact with INS, including STAT1, IFIT3, RSAD2, IFI44L, and XAF1, as verified by RT-qPCR. Meanwhile, seven miRNAs (miR-3680-3p, miR-3059-5p, miR-629-3p, miR-29b-2-5p, miR-514b-5p, miR-4755-5p, miR-4691-3p) were associated with regulation of hub genes. Notably, six hub genes (STAT1, IFIT3, RSAD2, ISG15, IFI44, IFI6) were down-regulated in cells exposed to both metformin and Mycobacterium tuberculosis antigens.

Conclusion

Network hub genes hold promise as disease status biomarkers and as metformin treatment targets for alleviating TB and DM. This study describes a strategy for exploring pathogenic mechanisms of diseases such as TB and DM.

SUBMITTER: Liu S

PROVIDER: S-EPMC8331204 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundBladder cancer (BLCA) is the ninth most common cancer worldwide, with high mortality and recurrence rates. Studies have increasingly reported that molecular diagnosis contributes to the early diagnosis and prognostic assessment of diseases. Thus, this study aims to find new biomarkers for the diagnosis and prognosis of BLCA.MethodsThe microarray datasets GSE147983 and The Cancer Genome Atlas (TCGA)-BLCA mRNA were obtained from the Gene Expression Omnibus (GEO) and TCGA. Differentially expressed genes (DEGs) were screened using the R "Limma" package. The "ClusterProfiler" package was used to conduct Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the DEGs. A DEG protein-protein interaction (PPI) network was constructed using the Search Tool for the Retrieval of Interacting Genes (STRING) database and visualized using Cytoscape. The functional module was reanalyzed using Cytoscape's Molecular Complex Detection ("MCODE") plugin, and key genes related to BLCA were identified via the "cytoHubba" plugin. Gene Expression Profiling Interactive Analysis 2 (GEPIA2) and the Tumor Immune Estimation Resource (TIMER) were used to verify the correlation between hub gene expression and immunity. A survival analysis of hub genes was performed using the Kaplan-Meier Plotter online tool.ResultsA total of 355 DEGs were screened out, including 236 upregulated and 119 downregulated DEGs. Some of the GO terms and pathways, such as chromosome separation, cell cycle, and cell senescence, were found to be significantly enriched in the DEGs. The key genes were kinesin family member 11 (KIF11), DLG associated protein 5 (DLGAP5), non-SMC condensin I complex subunit G (NCAPG), cell division cycle 20 (CDC20), cyclin B2 (CCNB2), BUB1 mitotic checkpoint serine (BUB1B), TPX2 microtubule nucleation factor (TPX2), NUF2 component of NDC80 kinetochore complex (NUF2), kinesin family member 2C (KIF2C), and cyclin B1 (CCNB1). Nine of them were immune-related, including KIF11, DLGAP5, NCAPG, CDC20, CCNB2, BUB1B, NUF2, KIF2C, and CCNB1. Survival analysis showed that the overexpression of BUB1B, CCNB1, CDC20, and DLGAP5 significantly reduced overall survival (OS) in patients with BLCA.ConclusionsThis study provided a theoretical basis for elucidating the pathogenesis and evaluating the prognosis of BLCA by screening potential biomarkers of BLCA.

Project description:The chronic metabolic disease named type 2 diabetes (T2D) accounts for over 90% of diabetes mellitus. An increasing number of evidences have revealed that hypoxia has a significantly suppressive effect on cell-mediated immunity, as well as the utilization of glucose in diabetics. Therefore, we aimed to screen and identify hypoxia-immune-related hub genes in T2D through bioinformatic analysis. The Gene Expression Omnibus (GEO) database was used to get T2D gene expression profile data in the peripheral blood samples (GSE184050), and hypoxia-related genes were acquired from Molecular Signatures Database (MSigDB). Differentially expressed mRNAs (DEGs) and lncRNAs (DELs) between T2D and normal samples were identified by DeSeq2 package. The clusterProfiler package was used to perform enrichment analyses for the overlapped genes of DEGs and hypoxia-related genes. Further, Hypoxia-related hub genes were discovered using two machine learning algorithms. Next, the compositional patterns of immune and stromal cells in T2D and healthy groups were estimated by using xCell algorithm. Moreover, we used the weighted correlation network analysis (WGCNA) to examine the connection between genes and immune cells to screen immune-related genes. Gene Set Enrichment Analysis (GSEA) was used to investigate the functions of the hypoxia-immune-related hub genes. Finally, two peripheral blood cohorts of T2D (GSE184050 and GSE95849) as well as the quantitative real-time PCR (qRT-PCR) experiments for clicinal peripheral blood samples with T2D were used for verification analyses of hub genes. And meanwhile, a lncRNA-TF-mRNA network was constructed. Following the differentially expressed analysis, 38 out of 3822 DEGs were screened as hypoxia-related DEGs, and 493 DELs were found. These hypoxia-related DEGs were mainly enriched in the GO terms of pyruvate metabolic process, cytoplasmic vesicle lumen and monosaccharide binding, and the KEGG pathways of glycolysis/gluconeogenesis, pentose phosphate pathway and biosynthesis of nucleotide sugars. Moreover, 7 out of hypoxia-related DEGs were identified as hub genes. There were six differentially expressed immune cell types between T2D and healthy samples, which were further used as the clinical traits for WGCNA to identify AMPD3 and IER3 as the hypoxia-immune-related hub genes. The results of the KEGG pathways of genes in high-expression groups of AMPD3 and IER3 were mainly concentrated in glycosaminoglycan degradation and vasopressin-regulated water reabsorption, while the low-expression groups of AMPD3 and IER3 were mainly associated with RNA degradation and nucleotide excision repair. Finally, when compared to normal samples, both the AMPD3 and IER3 were highly expressed in the T2D groups in the GSE184050 and GSE95849 datasets. The result of lncRNA-TF-mRNA regulatory network showed that lncRNAs such as BACH1-IT1 and SNHG15 might induce the expression of the corresponding TFs such as TFAM and THAP12 and upregulate the expression of AMPD3. This study identified AMPD3 and IER3 as hypoxia-immune-related hub genes and potential regulatory mechanism for T2D, which provided a new perspective for elucidating the upstream molecular regulatory mechanism of diabetes mellitus.

Project description:The underlying causes of esophageal cancer (EC) are unknown. To explore the molecular mechanisms that lead to EC, gene expression profiles of large cohorts of patients with EC were obtained from The Cancer Genome Atlas and the Gene Expression Omnibus (GEO) databases (GSE5364, GSE20347 and GSE23400). The present study identified 83 upregulated and 22 downregulated genes between EC and normal tissue using R statistical software and the GEO2R web tool. The Database for Annotation, Visualization and Integrated Discovery was used to identify the associated pathways, and for functional annotation of the differentially expressed genes (DEGs). Protein-protein interactions of these DEGs were analyzed based on the Search Tool for the Retrieval of Interacting Genes database, and hub genes were visualized using Cytoscape software. An online Kaplan-Meier plotter survival analysis tool was utilized to evaluate the prognostic value of hub gene expression in patients with EC. Further analysis of an additional dataset from GEO (GSE21293) revealed that these genes were associated with infiltration and metastasis in EC. In addition, the Gene Expression Profiling Interactive Analysis tool was used to evaluate expression levels of hub genes in patients with EC for different pathological stages. The Ualcan analysis tool was used to evaluate the expression levels of hub genes for different histological types. Overall, ubiquitin conjugating enzyme E2 C, cyclin dependent kinase inhibitor 3, CDC28 protein kinase regulatory subunit 2, kinesin family member 20A (KIF20A) and RAD51 associated protein 1 (RAD51AP1) were upregulated in EC tissues compared with normal tissues, and upregulation of these genes was a poor prognostic factor for patients with EC, indicating that these genes may mediate EC cell infiltration and metastasis. Among the hub genes, KIF-20A had potential value for predicting the pathological stage of EC. KIF20A and RAD51AP1 were more informative biomarkers of esophageal squamous cell carcinoma. Further studies are required to explore the value of these genes in the treatment of EC.

Project description:BackgroundPulmonary arterial hypertension (PAH) is a disease or pathophysiological syndrome which has a low survival rate with abnormally elevated pulmonary artery pressure caused by known or unknown reasons. In addition, the pathogenesis of PAH is not fully understood. Therefore, it has become an urgent matter to search for clinical molecular markers of PAH, study the pathogenesis of PAH, and contribute to the development of new science-based PAH diagnosis and targeted treatment methods.MethodsIn this study, the Gene Expression Omnibus (GEO) database was used to downloaded a microarray dataset about PAH, and the differentially expressed genes (DEGs) between PAH and normal control were screened out. Moreover, we performed the functional enrichment analyses and protein-protein interaction (PPI) network analyses of the DEGs. In addition, the prediction of miRNA and transcriptional factor (TF) of hub genes and construction miRNA-TF-hub gene network were performed. Besides, the ROC curve was used to evaluate the diagnostic value of hub genes. Finally, the potential drug targets for the 5 identified hub genes were screened out.Results69 DEGs were identified between PAH samples and normal samples. GO and KEGG pathway analyses revealed that these DEGs were mostly enriched in the inflammatory response and cytokine-cytokine receptor interaction, respectively. The miRNA-hub genes network was conducted subsequently with 131 miRNAs, 7 TFs, and 5 hub genes (CCL5, CXCL12, VCAM1, CXCR1, and SPP1) which screened out via constructing the PPI network. 17 drugs interacted with 5 hub genes were identified.ConclusionsThrough bioinformatic analysis of microarray data sets, 5 hub genes (CCL5, CXCL12, VCAM1, CXCR1, and SPP1) were identified from DEGs between control samples and PAH samples. Studies showed that the five hub genes might play an important role in the development of PAH. These 5 hub genes might be potential biomarkers for diagnosis or targets for the treatment of PAH. In addition, our work also indicated that paying more attention on studies based on these 5 hub genes might help to understand the molecular mechanism of the development of PAH.