Identification of common genes and pathways between type 2 diabetes and COVID-19.

Abstract

Numerous studies have reported a high incidence and risk of severe illness due to coronavirus disease 2019 (COVID-19) in patients with type 2 diabetes (T2DM). COVID-19 patients may experience elevated or decreased blood sugar levels and may even develop diabetes. However, the molecular mechanisms linking these two diseases remain unclear. This study aimed to identify the common genes and pathways between T2DM and COVID-19. Two public datasets from the Gene Expression Omnibus (GEO) database (GSE95849 and GSE164805) were analyzed to identify differentially expressed genes (DEGs) in blood between people with and without T2DM and COVID-19. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed on the common DEGs. A protein-protein interaction (PPI) network was constructed to identify common genes, and their diagnostic performance was evaluated by receiver operating characteristic (ROC) curve analysis. Validation was performed on the GSE213313 and GSE15932 datasets. A gene co-expression network was constructed using the GeneMANIA database to explore interactions among core DEGs and their co-expressed genes. Finally, a microRNA (miRNA)-transcription factor (TF)-messenger RNA (mRNA) regulatory network was constructed based on the common feature genes. In the GSE95849 and GSE164805 datasets, 81 upregulated genes and 140 downregulated genes were identified. GO and KEGG enrichment analyses revealed that these DEGs were closely related to the negative regulation of phosphate metabolic processes, the positive regulation of mitotic nuclear division, T-cell co-stimulation, and lymphocyte co-stimulation. Four upregulated common genes (DHX15, USP14, COPS3, TYK2) and one downregulated common feature gene (RIOK2) were identified and showed good diagnostic accuracy for T2DM and COVID-19. The AUC values of DHX15, USP14, COPS3, TYK2, and RIOK2 in T2DM diagnosis were 0.931, 0.917, 0.986, 0.903, and 0.917, respectively. In COVID-19 diagnosis, the AUC values were 0.960, 0.860, 1.0, 0.9, and 0.90, respectively. Validation in the GSE213313 and GSE15932 datasets confirmed these results. The miRNA-TF-mRNA regulatory network showed that TYH2 was targeted by PITX1, PITX2, CRX, NFYA, SREBF1, RELB, NR1L2, and CEBP, whereas miR-124-3p regulates THK2, RIOK2, and USP14. We identified five common feature genes (DHX15, USP14, COPS3, TYK2, and RIOK2) and their co-regulatory pathways between T2DM and COVID-19, which may provide new insights for further molecular mechanism studies.