Identification of hub genes related to the progression of type 1 diabetes by computational analysis.
Ontology highlight
ABSTRACT: Background
Type 1 diabetes (T1D) is a serious threat to childhood life and has fairly complicated pathogenesis. Profound attempts have been made to enlighten the pathogenesis, but the molecular mechanisms of T1D are still not well known.Methods
To identify the candidate genes in the progression of T1D, expression profiling by high throughput sequencing dataset GSE123658 was downloaded from Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were identified, and gene ontology (GO) and pathway enrichment analyses were performed. The protein-protein interaction network (PPI), modules, target gene - miRNA regulatory network and target gene - TF regulatory network analysis were constructed and analyzed using HIPPIE, miRNet, NetworkAnalyst and Cytoscape. Finally, validation of hub genes was conducted by using ROC (Receiver operating characteristic) curve and RT-PCR analysis. A molecular docking study was performed.Results
A total of 284 DEGs were identified, consisting of 142 up regulated genes and 142 down regulated genes. The gene ontology (GO) and pathways of the DEGs include cell-cell signaling, vesicle fusion, plasma membrane, signaling receptor activity, lipid binding, signaling by GPCR and innate immune system. Four hub genes were identified and biological process analysis revealed that these genes were mainly enriched in cell-cell signaling, cytokine signaling in immune system, signaling by GPCR and innate immune system. ROC curve and RT-PCR analysis showed that EGFR, GRIN2B, GJA1, CAP2, MIF, POLR2A, PRKACA, GABARAP, TLN1 and PXN might be involved in the advancement of T1D. Molecular docking studies showed high docking score.Conclusions
DEGs and hub genes identified in the present investigation help us understand the molecular mechanisms underlying the advancement of T1D, and provide candidate targets for diagnosis and treatment of T1D.
SUBMITTER: Prashanth G
PROVIDER: S-EPMC8028841 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundType 2 diabetes mellitus (T2DM) is one of the most common chronic diseases in the world with complicated pathogenesis. This study aimed to identify differentially expressed genes (DEGs) and molecular pathways in T2DM using bioinformatics analysis.Materials and methodsTo explore potential therapeutic targets for T2DM, we analyzed three microarray datasets (GSE50397, GSE38642, and GSE44035) acquired from the Gene Expression Omnibus (GEO) database. DEGs between T2DM islet and normal islet were picked out by the GEO2R tool and Venn diagram software. Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) to identify the pathways and functional annotation of DEGs. Then, protein-protein interaction (PPI) of these DEGs was visualized by Cytoscape with Search Tool for the Retrieval of Interacting Genes/Proteins (STRING).ResultsIn total, we identified 36 DEGs in the three datasets, including 32 up-regulated genes and four down-regulated genes. The improved functions and pathways of the DEGs enriched in cytokine-cytokine receptor interaction, pathways in cancer, PI3K-Akt signaling pathway, and Rheumatoid arthritis. Among them, ten hub genes with a high degree of connectivity were selected. Furthermore, via the re-analysis of DAVID, four genes (IL6, MMP3, MMP1, and IL11) were significantly enriched in the Rheumatoid arthritis pathway.ConclusionOur study, based on the GEO database, identified four significant up-regulated DEGs and provided novel targets for diagnosis and treatment of T2DM.
| S-EPMC7250707 | biostudies-literature
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.
| S-EPMC10160629 | biostudies-literature
Project description:BackgroundColorectal cancer (CRC) is the third most common cancer in the world. The present study is aimed at identifying hub genes associated with the progression of CRC.MethodThe data of the patients with CRC were obtained from the Gene Expression Omnibus (GEO) database and assessed by weighted gene co-expression network analysis (WGCNA), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses performed in R by WGCNA, several hub genes that regulate the mechanism of tumorigenesis in CRC were identified. Differentially expressed genes in the data sets GSE28000 and GSE42284 were used to construct a co-expression network for WGCNA. The yellow, black and blue modules associated with CRC level were filtered. Combining the co-expression network and the PPI network, 15 candidate hub genes were screened.ResultsAfter validation using the TCGA-COAD dataset, a total of 10 hub genes (MT1X, MT1G, MT2A, CXCL8, IL1B, CXCL5, CXCL11, IL10RA, GZMB, KIT) closely related to the progression of CRC were identified. The expressions of MT1G, CXCL8, IL1B, CXCL5, CXCL11 and GZMB in CRC tissues were higher than normal tissues (p-value < 0.05). The expressions of MT1X, MT2A, IL10RA and KIT in CRC tissues were lower than normal tissues (p-value < 0.05).ConclusionsBy combinating with a series of methods including GO enrichment analysis, KEGG pathway analysis, PPI network analysis and gene co-expression network analysis, we identified 10 hub genes that were associated with the progression of CRC.
| S-EPMC7664468 | biostudies-literature
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.
| S-EPMC7491038 | biostudies-literature
Project description:BackgroundType 2 Diabetes Mellitus (T2DM) is an independent risk factor of hepatocellular carcinoma (HCC). However, the related genes and modules to hepatocarcinogenesis and progression in T2DM remain unclear.MethodsThe microarray data from Gene Expression Omnibus (GEO) were analyzed to screen differentially expressed genes (DEGs) of T2DM and HCC dataset. Then, weighted gene co-expression network analysis (WGCNA) was performed on these DEGs to detect the modules and genes, respectively. Common genes in modules with clinical interests of T2DM and HCC were obtained and annotated via GOSemSim package and Metascape. Genes related to late-stage HCC and high glycated haemoglobin (HbA1c) were also identified. These genes were validated by UALCAN analysis and univariate cox regression based on The Cancer Genome Atlas (TCGA). Finally, another two independent datasets were applied to confirm the results of our study.ResultsA total of 1288 and 1559 DEGs of T2DM and HCC were screened, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment revealed several shared pathways in two diseases, such as pathways in cancer and metabolism. A total of 37 common genes correlated with T2DM and HCC were then identified with WGCNA. Furthermore, 12 genes from modules associated with late-stage HCC and high HbA1c were regarded as hub genes. Among these genes, 8 genes associated with tumor invasion and metastasis were validated by UALCAN analysis. Moreover, downregulations of ACAT1, SLC2A2, PCK1 and ABAT were significantly associated with poorer prognosis in HCC patients with elevated HbA1c. Additionally, the expressions of PCK1 and ABAT were raised in HepG2 cells pre-treated with metformin and phenformin.ConclusionsThe present study confirmed several metabolic genes related to hyperglycemia and malignant tumor, which may provide not only new insights into the pathogenesis of hepatocarcinogenesis and progression in T2DM, but also novel therapeutic targets for T2DM patients with HCC in the future.
| S-EPMC8053303 | biostudies-literature
Project description:Hepatocellular carcinoma (HCC) is the fourth leading cause of cancer‑related deaths among cancer patients. Genes correlated with the progression and prognosis of HCC are critically needed to be identified. In the present study, 3 Gene Expression Omnibus (GEO) datasets (GSE46408, GSE65372 and GSE84402) were used to analyze the differentially expressed genes (DEGs) between HCC and non‑tumor liver tissues. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were conducted to clarify the functional roles of DEGs. A protein‑protein interaction network was established to screen the hub genes associated with HCC. The prognostic values of hub genes in HCC patients were analyzed using The Cancer Genome Atlas (TCGA) database. The expression levels of hub genes were validated based on ONCOMINE, TCGA and Human Protein Atlas (HPA) databases. Notably, 56 upregulated and 33 downregulated DEGs were markedly enriched under various GO terms and four KEGG terms. Among these DEGs, 10 hub genes with high connectivity degree were identified, including cyclin B1, cyclin A2, cyclin B2, condensin complex subunit 3, PDZ binding kinase, nucleolar and spindle‑associated protein 1, aurora kinase A, ZW10 interacting kinetochore protein, protein regulator of cytokinesis 1 and kinesin family member 4A. The upregulated expression levels of these hub genes in HCC tissues were further confirmed by ONCOMINE, TCGA, and HPA databases. Additionally, the increased mRNA expression of each hub gene was related to the unfavorable disease‑free survival and overall survival of HCC patients. The present study identified ten genes associated with HCC, which may help to provide candidate targets for the diagnosis and treatment of HCC.
| S-EPMC6908929 | biostudies-literature
Project description:Spinal cord injury (SCI) is a disturbance of peripheral and central nerve conduction that causes disability in sensory and motor function. Currently, there is no effective treatment for SCI. Mitophagy plays a vital role in mitochondrial quality control during various physiological and pathological processes. The study aimed to elucidate the role of mitophagy and identify potential mitophagy-related hub genes in SCI pathophysiology. Two datasets (GSE15878 and GSE138637) were analyzed. Firstly, the differentially expressed genes (DEGs) were identified and mitophagy-related genes were obtained from GeneCards, then the intersection between SCI and mitophagy-related genes was determined. Next, we performed gene set enrichment analysis (GSEA), weighted gene co-expression network analysis (WGCNA), protein-protein interaction network (PPI network), least absolute shrinkage and selection operator (LASSO), and cluster analysis to identify and define the hub genes in SCI. Finally, the link between hub genes and infiltrating immune cells was investigated and the potential transcriptional regulation/small molecular compounds to target hub genes were predicted. In total, SKP1 and BAP1 were identified as hub genes of mitophagy-related DEGs during SCI development and regulatory T cells (Tregs)/resting NK cells/activated mast cells may play an essential role in the progression of SCI. LINC00324 and SNHG16 may regulate SKP1 and BAP1, respectively, through miRNAs. Eleven and eight transcriptional factors (TFs) regulate SKP1 and BAP1, respectively, and six small molecular compounds target BAP1. Then, the mRNA expression levels of BAP1 and SKP1 were detected in the injured sites of spinal cord of SD rats at 6 h and 72 h after injury using RT-qPCR, and found that the level were decreased. Therefore, the pathways of mitophagy are downregulated during the pathophysiology of SCI, and SKP1 and BAP1 could be accessible targets for diagnosing and treating SCI.
| S-EPMC10875195 | biostudies-literature
Project description:IntroductionAutophagy, an innate safeguard mechanism for protecting the organism against harmful agents, is implicated in the survival of pancreatic â cells and the development of type 2 diabetes mellitus (T2DM). Potential autophagy-related genes (ARGs) may serve as potential biomarkers for T2DM treatment.MethodsThe GSE25724 dataset was downloaded from the Gene Expression Omnibus (GEO) database, and ARGs were obtained from the Human Autophagy Database. The differentially expressed autophagy-related genes (DEARGs) were screened at the intersection of ARGs and differentially expressed genes (DEGs) between T2DM and non-diabetic islet samples, which were subjected to functional enrichment analyses. A protein-protein interaction (PPI) network was constructed to identify hub DEARGs. Expressions of top 10 DEARGs were validated in human pancreatic â-cell line NES2Y and rat pancreatic INS-1 cells using quantitative reverse transcription polymerase chain reaction (qRT-PCR). Cell viability and insulin secretion were measured after cell transfection with lentiviral vector EIF2AK3 or RB1CC1 into islet cells.ResultsIn total, we discovered 1,270 DEGs (266 upregulated and 1,004 downregulated genes) and 30 DEARGs enriched in autophagy- and mitophagy-related pathways. In addition, we identified GAPDH, ITPR1, EIF2AK3, FOXO3, HSPA5, RB1CC1, LAMP2, GABARAPL2, RAB7A, and WIPI1 genes as the hub ARGs. Next, qRT-PCR analysis revealed that expressions of hub DEARGs were consistent with findings from bioinformatics analysis. EIF2AK3, GABARAPL2, HSPA5, LAMP2, and RB1CC1 were both differentially expressed in the two cell types. Overexpression of EIF2AK3 or RB1CC1 promoted cell viability of islet cells and increased the insulin secretion.DiscussionThis study provides potential biomarkers as therapeutic targets for T2DM.
| S-EPMC10200926 | biostudies-literature
Project description:BackgroundType 2 diabetes mellitus (T2DM) and its related complications contribute to the high morbidity and mortality in worldwide. Skeletal muscle insulin resistance plays a critical role in the onset of T2DM due to the decreasing in the insulin-stimulated glucose uptake. T2DM is associated not only with the inherited factors but also with the noninherited factors. However, the susceptibility genes related with the two factors and the transcription factors (TF) regulating the susceptibility genes in skeletal muscle, which aggravate the development of T2DM were still ill-defined.MethodsIn the present study, the expression profiles by the array of GSE25462 were retrieved from the GEO database. GEO2R was performed to validate the susceptibility differentially expressed genes (SDEG) in skeletal muscle of T2DM. Gene Ontology (GO) analysis and The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted via The Database for Annotation, Visualization, and Integrated Discovery (DAVID). A Protein-Protein Interaction (PPI) network was performed with the STRING.ResultsWith the performance of GEO2R, 229 SDEGs in skeletal muscle of T2DM were identified. The biological processes (BP) of SDEGs was enriched in the cellular response to UV-B most significantly. KEGG pathway analysis revealed that the SDEGs were most significantly enriched in glycosaminoglycan degradation. 5 hub susceptibility genes (GPR84, CALCB, GCG, PTGDR, GNG8) in the skeletal muscle of T2DM were identified. Eventually, the common transcription factors regulating the hub susceptibility genes were identified by means of the online tool PROMO.ConclusionsFive hub susceptibility genes (GPR84, CALCB, GCG, PTGDR, GNG8) in the skeletal muscle of T2DM and the common transcription factors were identified. The outputs would provide new clues on the novel potential targets and the therapeutic strategies for treating T2DM and its related diseases.
| S-EPMC9652898 | biostudies-literature
Project description:BACKGROUND Hepatocellular carcinoma (HCC) is one of the most prevalent cancers in the world. Bioinformatics studies have been widely used for screening genes involved in the initiation and progression of HCC. MATERIAL AND METHODS We obtained liver cancer microarray raw data from the GEO database (GSE54238). Next, weighted gene co-expression network analysis (WGCNA) was used to assess the critical modules. Then, we assessed the gene significance by calculating survival, expression level, and receiver operating characteristic (ROC) in the TCGA database. We also validated the expression of selected genes in the Oncomine database and calculated the relationship between 4 hub genes and immune infiltration. Finally, GSEA enrichment analysis was used to explore the potential mechanism. RESULTS We identified the red and blue modules as the critical modules, and found 176 candidate genes by assessing gene significance. GO and KEEG results suggested that the candidate genes are involved in the cell cycle. Four hub genes - SOX4, STK39, TARBP1, and TDRKH - were eventually screened after validating their expression and power in diagnosing HCC in the TCGA database. Immune infiltration analysis and GSEA enrichment analysis showed that these 4 hub genes were correlated with the immune cell populations infiltration and that multiple mechanisms were involved, such as angiogenesis and epithelial-mesenchymal transition. CONCLUSIONS Our findings revealed that these 4 genes can be regarded as potential prognosticators and therapeutic targets for HCC.
| S-EPMC7111148 | biostudies-literature