Bioinformatics Analysis of next generation sequencing data for Risk Prediction in Patients with Type 1 diabetes mellitus

Abstract

Type 1 diabetes mellitus (T1DM) comprise the most common forms of autoimmune disease. The aim of this investigation was to apply a bioinformatics approach to reveal related pathways or genes involved in the development of T1DM. The next generation sequencing (NGS) dataset GSE182870 was downloaded from the gene expression omnibus (GEO) database. Differentially expressed gene (DEG) analysis was performed using DESeq2. The g:Profiler was utilized to analyze the functional enrichment, gene ontology (GO) and REACTOME pathway of the differentially expressed genes. Protein-protein interaction (PPI) network, modules, miRNA- hub gene regulatory network, and TF-hub gene regulatory network were conducted via comprehensive target prediction and network analyses. Finally, hub genes were validated by using receiver operating characteristic curve analysis. A total of 860 DEGs were screened out from NGS dataset, among which 477 genes were up regulated and 383 genes were down regulated. GO enrichment analysis indicated that up regulated genes were mainly involved in cellular metabolic process, intracellular anatomical structure and catalytic activity, and the down regulated genes were significantly enriched in cellular nitrogen compound biosynthetic process, protein containing complex and heterocyclic compound binding. REACTOME pathway enrichment analysis showed that the up regulated genes were mainly enriched in metabolism of carbohydrates and the down regulated genes were significantly enriched in metabolism of RNA. A PPI network, modules, miRNA-hub gene regulatory network and TF - hub gene regulatory network were constructed, and hub genes (MAPK14, RHOC, MAD2L1, TAF1, TRAF2, HSP90AA1, TP53, HSP90AB1, UBA52 and RACK1) were identified. This study provides further insights into the underlying pathogenesis of T1DM.