Ya Wang
Kai Li3
Shuangyang Mo
Peishan Yao
Shanyu Qin- 1Gastroenterology Department, The First Affiliated Hospital of Guangxi Medical University, Nanning, China
- 2Endocrinology Department, Liuzhou Peoples’ Hospital Affiliated to Guangxi Medical University, Liuzhou, China
- 3Orthopedics Department, The Fourth Affiliated Hospital of Guangxi Medical University, Liuzhou, China
- 4Gastroenterology Department, Liuzhou Peoples’ Hospital Affiliated to Guangxi Medical University, Liuzhou, China
- 5Department of Traumatic Surgery, Microsurgery, and Hand Surgery, Guangxi Zhuang Autonomous Region People’s Hospital, Nanning, Guangxi, China
- 6Department of Pharmacy, Affiliated Tumor Hospital of Guangxi Medical University, Nanning, Guangxi, China
Background: 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.
Methods: 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.
Results: 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.
Conclusion: 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.
1 Introduction
Type 2 diabetes mellitus (T2DM) is a highly prevalent disease worldwide. The number of people with diabetes was 285 million globally in 2010 and 463 million in 2019, and it is estimated to reach 710 million by 2045 (Shaw et al., 2010; Grewal et al., 2020; Sharma et al., 2022). Coronavirus disease 2019 (COVID-19) is an acute respiratory disease caused by SARS-CoV-2 infection. COVID-19 is highly infectious, shows rapid mutation, and has become a global pandemic, resulting in high morbidity and mortality due to complications (To et al., 2021). Previous studies have shown a significant genetic correlation between T2DM and COVID-19 (Wu K. C. H. et al., 2022; Ni et al., 2023). Furthermore, Mendelian randomization has demonstrated a causal relationship between genetic susceptibility to T2DM and SARS-CoV-2 infection and hospitalization due to COVID-19 (Host Genetics Initiative, 2022). Studies have reported that the incidence and severity of COVID-19 are higher in patients with diabetes (Leon-Abarca et al., 2021; Sharma et al., 2022). Oxidative stress and chronic inflammation are pathological mechanisms underlying the occurrence and development of T2DM and its related chronic complications (Wronka et al., 2022). Chronic inflammation leads to insulin resistance, causing the progression of T2DM, which has pro-inflammatory characteristics. Research has shown that hyperglycemia can stimulate immune cells and increase pro-inflammatory cytokines such as tumor necrosis factor α, interleukin 1β, and interleukin −6 (IL-6) (Muniyappa and Gubbi, 2020). Inflammatory markers such as IL-6, C-reactive protein, and ferritin in the blood of diabetic patients are higher than those in non-diabetic patients (Guo et al., 2020). Inflammation and immunity are involved in the progression of diabetic nephropathy (Chen et al., 2022; Rayego-Mateos et al., 2023). Mitochondrial metabolism and immune inflammation are critical factors in the pathogenesis of diabetic cardiomyopathy (Chen et al., 2022). A retrospective cohort study showed that people with COVID-19 combined with diabetes had a more robust inflammatory immune response and increased in-hospital mortality rate due to cardiac damage (Bo et al., 2022). This indicates that inflammation and immunity are closely related to T2DM and its complications. Patients with diabetes combined with COVID-19 may be more susceptible to the cytokine storm, leading to disease exacerbation. Additionally, diabetic patients have higher levels of D-dimer compared to non-diabetic controls, indicating that they are in a hypercoagulable state (Guo et al., 2020). After COVID-19 infection, the coagulation cascade is activated, exacerbating the hypercoagulable state and leading to adverse outcomes of COVID-19.
Research has shown that COVID-19 may affect the pathophysiology of diabetes. Patients with diabetes and COVID-19 may experience elevated blood glucose, possibly due to stress (Wang W. et al., 2020). Additionally, nearly 10% of COVID-19 patients with T2DM experience hypoglycemia (Zhou et al., 2020), and COVID-19 may induce the development of diabetes in people without a history of the disease (Khunti et al., 2021). The progression and prognosis of both diabetes and COVID-19 are interrelated. However, little is known about the molecular pathways underlying the abnormal inflammatory and immune responses in patients with diabetes and COVID-19. Given the vast population of patients with T2DM combined with COVID-19, the complex relationship between T2DM and COVID-19, the lack of dual-purpose treatment strategies for these diseases, and the limited research on common pathways between T2DM and COVID-19, we employed various bioinformatics analysis methods to study the shared genes and pathways between T2DM and COVID-19, providing clues for the development of dual-purpose preventive and therapeutic strategies.
2 Materials and methods
2.1 Microarray chip data
We downloaded datasets from the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo) that met the following inclusion criteria: 1. Studies published between 2012 and 2023; 2. Studies on blood samples from patients with T2DM and healthy controls and those with blood samples from COVID-19 patients and healthy controls; 3. Datasets with a sample size ≥5. We excluded studies with duplicate samples and incomplete data. Two T2DM datasets were included, GSE95849 and GSE15932; two COVID-19 datasets were included, GSE164805 and GSE21331. A flowchart of this study is presented in Figure 1. GSE95849 contains mRNA expression profiles of 18 blood samples, including 12 patients with T2DM and six healthy controls. Detailed information on this dataset is shown in Table 1. GSE15932 contains mRNA expression profiles on 15 peripheral blood mononuclear cell samples, including 10 COVID-19 patients and five healthy controls. GSE213313 contains mRNA expression profiles on 94 blood samples, including 83 COVID-19 patients and 11 healthy controls. GSE15932 contains mRNA expression profiles on 16 blood samples, including eight T2DM patients and 11 healthy controls. We analyzed the GSE95849 and GSE164805 datasets and validated the results in the GSE15932 and GSE213313 datasets.
Figure 1. Flow chart of the study.
Table 1. Details of the GEO datasets.
2.2 Identification of differentially expressed genes shared by T2DM and COVID-19
We performed standardization and preprocessing on the data using the “limma” package in R (v.4.3.0), removed gene probes with missing values and duplicates, took the log2 transformation of the untransformed data, and conducted differential expression analysis to determine the differentially expressed genes (DEGs) between the disease groups and healthy control groups. We used the Benjamini–Hochberg method to correct the p-values for DEGs. We identified DEGs based on the corrected p < 0.05 and |log fold change (FC) |≥ 0.085. We took the intersection of the DEGs from both analyses to obtain the shared DEGs. We visualized the results using a Venn diagram.
2.3 Pathway enrichment analysis
Using the “clusterProfiler” package in R, we performed Gene Ontology (GO) enrichment analysis on the shared DEGs to determine enriched biological processes (BP), molecular functions (MF), and cellular components (CC). We also analyzed Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways to identify shared pathways between T2DM and COVID-19. We set the filtering criteria as p < 0.05 and visualized the results using the “ggpubr” package.
2.4 Protein-protein interaction network construction and identification of common feature genes
We used the online database STRING (http://www.string-db.org/) to construct a protein-protein interaction network (PPI) to predict the functional interactions between proteins. We visualized the results using Cytoscape software and calculated Maximal Clique Centrality scores (Wu J. et al., 2022). We identified hub genes using the CytoHubba plugin in Cytoscape and assessed the most critical modules in the PPI network using the MCODE plugin in Cytoscape.
2.5 ROC curve analysis of common genes
To evaluate the sensitivity and specificity of the shared DEGs in diagnosing T2DM and COVID-19, we used the pROC package in R to perform receiver operating characteristic (ROC) curve analysis. We used the area under the curve (AUC) to evaluate the diagnostic performance of the common genes for T2DM and COVID-19. We also validated the diagnostic performance using ROC curve analysis in the validation dataset.
2.6 Interaction network of common genes
We used the GeneMANIA database (http://genemania.org/) to construct interaction networks of the five common feature genes and their co-expressed genes.
2.7 Correlation analysis of shared DEGs
We used the “corrplot” package in R to perform a correlation analysis of the common feature genes and calculate the correlation coefficients between genes to observe the regulatory relationships among the five common feature genes.
2.8 Regulatory network analysis
We used the miRTarBase (Chou et al., 2018), Starbase (Yang et al., 2011), and Targetscan (Agarwal et al., 2015) databases to predict the target microRNAs (miRNAs) of the shared core genes. We only retained miRNAs predicted by all three databases. We predicted the target transcription factors (TFs) of the shared core genes using the Enrichr website (https://maayanlab.cloud/Enrichr/) by inputting the names of the shared core genes and selecting the “TRANSFAC and JASPAR PWMS” option. We downloaded the data table, retained the TFs with a p < 0.05, and restricted them to genes from humans. Using Cytoscape, we visualized the miRNA-TF-messenger RNA (mRNA) regulatory network.
2.9 Candidate drugs
We used the DGIdb website (DGIdb—Mining the Druggable Genome) to identify candidate drugs based on the integrated score obtained by inputting the common feature genes.
3 Results
3.1 Shared DEGs between T2DM and COVID-19
We performed differential expression analysis and identified 973 DEGs in the T2DM dataset GSE9584 and 143,344 DEGs in the COVID-19 dataset GSE164805. We took the intersection of the DEGs from both datasets and identified 221 shared DEGs (Figure 2), including 81 upregulated genes (Figure 2B) and 140 downregulated genes (Figure 2A).
Figure 2. (A) Common upregulated genes in COVID-19 and T2DM; (B) Common downregulated genes in COVID-19 and T2DM.
3.2 Functional and enrichment analysis of shared DEGs
The GO-BP results showed that these shared DEGs are mainly associated with biological processes such as “negative regulation of phosphate metabolic process,” “negative regulation of phosphorus metabolic process,” “rRNA metabolic process,” “positive regulation of mitotic nuclear division,” “T cell co-stimulation,” “lymphocyte co-stimulation,” “fibrinolysis,” and “alpha-beta T cell activation.” In the GO-CC category, the shared DEGs were enriched in areas such as “secretory granule lumen,” “cytoplasmic vesicle lumen,” “vesicle lumen,” “striated muscle thin filament,” and “myofilament.” In the GO-MF category, the shared DEGs were mainly involved in functions such as “protein serine/threonine phosphatase activity,” “protein tyrosine kinase activity,” and “microtubule binding” (Figure 3A). KEGG enrichment analysis revealed that these shared DEGs were primarily enriched in pathways such as “human T-cell leukemia virus 1 infection,” “type I diabetes mellitus,” “circadian rhythm,” “biosynthesis of nucleotide sugars,” “amino sugar and nucleotide sugar metabolism,” “Th1 and Th2 cell differentiation,” “gastric cancer,” “Cushing’s syndrome,” and “T-cell receptor signaling pathway” (Figure 3B). These results suggest that metabolism and the inflammatory response play essential roles in the crosstalk between T2DM and COVID-19.
Figure 3. (A) GO pathways of shared DEGs between T2DM and COVID-19; (B) KEGG pathways of shared DEGs between T2DM and COVID-19.
3.3 PPI network and extraction of core genes
We imported the shared DEGs of T2DM and COVID-19 into the STRING database, with a minimum required interaction source set to 0.4 and free nodes hidden, resulting in a PPI network with 219 nodes and 133 edges. We visualized the network using Cytoscape software, with red and blue nodes representing upregulated and downregulated genes, respectively (Figure 4). Using the MCODE plugin in Cytoscape, we identified five important modules with parameter settings of Degree Cutoff: 2, Node Score Cutoff: 0.2, K-Core: 2, and Max. Depth: 100. We calculated the Maximal Clique Centrality scores and subsequently used the Cytohubba plugin in Cytoscape to identify high-connectivity hub genes in the network. We identified five common feature genes between T2DM and COVID-19: DHX15, USP14, COPS3, TYK2, and RIOK2 (where DHX15, USP14, COPS3, and TYK2 were upregulated and RIOK2 was downregulated).
Figure 4. PPI network of shared DEGs.
3.4 ROC curves of common feature genes
We evaluated the diagnostic performance of the common feature genes for T2DM and COVID-19 by plotting ROC curves. The results from the ROC curve analysis showed that DHX15, USP14, COPS3, TYK2, and RIOK had AUC values of 0.931, 0.917, 0.986, 0.903, and 0.917, respectively, for T2DM (Figure 5A) and AUC values of 0.960, 0.860, 1.0, 0.9, and 0.90, respectively, for COVID-19 (Figure 5B). In the validation dataset, the AUC values of DHX15, USP14, COPS3, TYK2, and RIOK were 0.922, 0.781, 0.812, 0.844, and 0.984, respectively, for T2DM (Figure 5C), and 0.733, 0.687, 0.737, 0.782, and 0.632, respectively, for COVID-19 (Figure 5D). These results suggest that these five shared core genes can serve as effective diagnostic biomarkers for distinguishing between T2DM and non-diabetes and COVID-19 and healthy controls.
Figure 5. (A) The ROC curve of the diagnostic efficacy in GSE95849; (B) The ROC curve of the diagnostic efficacy in GSE164805; (C) The ROC curve of the diagnostic efficacy in GSE15932; (D) The ROC curve of the diagnostic efficacy in GSE213313.
3.5 Interactome network of common feature genes and their co-expressed genes
GeneMANIA analysis revealed a complex interactome network of common feature genes and their co-expressed genes between T2DM and COVID-19, with “physical interactions” accounting for 77.64%, “predicted” for 5.37%, “co-expression” for 8.01%, “co-localization” for 3.63%, “pathway” for 1.88%, “genetic interactions” for 2.87%, and “shared protein domains” for 0.60% of all interactions. The main biological functions of this network were related to processes such as protein modification by small protein removal, nucleotide-excision repair, and response to type I interferon (Figure 6; Supplementary Material). These results further emphasize the potential involvement of metabolism and inflammatory response and their associated pathways in the development and progression of both T2DM and COVID-19.
Figure 6. Interactome network of common feature genes and their co-expressed genes.
3.6 Correlation analysis of common feature genes
We performed a correlation analysis of the common feature genes using the “corrplot” package in R. Through this analysis, we observed a significant positive correlation between DHX15 and USP14 expression, a significant positive correlation between COPS3 and RIOK2 expression, and a significant negative regulatory relationship between TYK2 and COPS3 as well as RIOK2 expression. This regulatory trend was consistent in both T2DM and COVID-19 (Figure 7).
Figure 7. (A) Regulatory relationships of common feature genes in T2DM; (B) Regulatory relationships of common feature genes in COVID-19 (green represents negative correlation and red represents positive correlation).
3.7 MiRNA-TF-mRNA regulatory network analysis
We retained all predicted miRNAs in miRTarBase, Starbase, and Targetsca, as well as all predicted TFs from the Enrich databases TRANSFAC and JASPAR PWMS for the common feature genes. Using Cytoscape, we visualized the miRNA-TF-mRNA regulatory network, which consisted of five common feature genes, 11 predicted miRNAs, and 12 predicted TFs (Figure 8).
Figure 8. miRNA-TF-mRNA regulatory network. The red oval represents the common core genes, the green triangle represents the predicted miRNA, and the purple rectangle represents the predicted TFs. (Abbreviations: mRNA, messenger RNA; miRNA, microRNA; TF, transcription factor).
We observed that TYK2 was targeted by PITX1, PITX2, CRX, NFYA, SREBF1, RELB, NR1L2, and CEBPB, whereas COPS3 was targeted by NR1L2, RELB, FOS, NFAT2, and NFYA. We also found that miR-124-3p simultaneously regulates THK2, RIOK2, and USP14; miR-221-3p regulates DHX15 and USP14; miR-101-3p regulates RIOK2 and USP14; and miR-155-5p regulates COPS3 and USP14.
3.8 Candidate drugs
We entered the common feature genes into the DGIdb website and obtained results for candidate drugs. The analysis revealed that the non-selective Tyk2 inhibitor, Baricitinib, had the highest overall score (query score of 2.15 and interaction score of 1.55), followed by the Janus kinase inhibitor Tofacitinib (query score of 1.84 and interaction score of 1.32) as potential drugs.
4 Discussion
Evidence suggests that patients with T2DM have a higher incidence, severity, and mortality rate of COVID-19. The mechanisms underlying this relationship are not fully understood but may be related to a dysregulated inflammatory response (Sardu et al., 2020; Nyland et al., 2021). There is limited research exploring the shared molecular mechanisms between T2DM and COVID-19 at the gene level. This study explored the molecular biological functions and pathways common to T2DM and COVID-19, providing clues for developing dual-purpose prevention and treatment strategies.
We analyzed transcriptome data of T2DM and COVID-19 blood samples and identified 221 shared DEGs. GO and KEGG analyses were performed on these DEGs. GO analysis showed that the shared DEGs were closely related to the negative regulation of phosphate and phosphorus metabolic processes, rRNA metabolic processes, positive regulation of mitotic nuclear division, T-cell co-stimulation, and lymphocyte co-stimulation. KEGG enrichment analysis indicated that these shared DEGs were predominantly involved in human T-cell leukemia virus 1 infection, type I diabetes, nucleotide sugar biosynthesis, Th1 and Th2 cell differentiation, and T-cell receptor signaling pathways. Many studies have shown that T-cell function is impaired in patients with T2DM. Our findings suggest that the shared molecular mechanisms between T2DM and COVID-19 may be related to dysregulated T-cell function and nucleotide metabolism. Further investigation of these pathways may provide a basis for developing new prevention and treatment strategies for both diseases (Moura et al., 2017). In obese patients with diabetes, there may be an elevation of activated CD4, CD278 co-stimulatory T-cells, Th17 cells, and cytotoxic T-cells (Richard et al., 2017). Th17 cells are critical in recruiting and invading microorganisms through immune cell recruitment and phagocytosis (Marqués et al., 2012). In addition, Farnsworth et al. reported that patients with T2DM exhibit defects in humoral immunity (Farnsworth et al., 2018). There is evidence to suggest that patients with T2DM exhibit decreased cytokine expression and impaired T-cell responses following viral stimulation (Richard et al., 2017; Tong et al., 2021). CD4 and CD8 T-cells have been detected in patients with COVID-19 (Grifoni et al., 2020; Tong et al., 2021). This suggests the importance of T-cell responses in COVID-19, and we speculate that the dysfunction of T-cells and lymphocytes in T2DM may affect the progression of COVID-19 disease.
We constructed a PPI network and identified five common feature genes (DHX15, USP14, COPS3, TYK2, and RIOK) of T2DM and COVID-19 by extracting the core genes of the network. We used ROC curve analysis to evaluate diagnostic performance, and the AUC values for T2DM diagnosis were 0.931, 0.917, 0.986, 0.903, and 0.917 for DHX15, USP14, COPS3, TYK2, and RIOK2, respectively. The AUC values for COVID-19 diagnosis were 0.960, 0.860, 1.0, 0.9, and 0.90, respectively. We also validated these results using independent T2DM and COVID-19 datasets without overlapping data with the training set, which showed good diagnostic performance. These findings suggest that DHX15, USP14, COPS3, TYK2, and RIOK2 play important roles in developing T2DM and COVID-19.
DHX15 is a member of the RNA helicase DEAH-box family, which is widely expressed in immune cells (Mosallanejad et al., 2014; Detanico et al., 2019). DHX15 has both ATPase and helicase activities (Semlow et al., 2016; Bohnsack et al., 2022; Wang et al., 2022) and plays a role in various cellular processes, including pre-mRNA splicing and ribosomal RNA synthesis (Bohnsack et al., 2021). Most current research on DHX15 is focused on cancer and immunity. Mutations in the SF3B1 splicing factor gene can disrupt the interaction between the splicing factor and SUGP1, leading to splicing errors that can cause cancer. DHX15, as an RNA helicase, is required for SUGP1 and participates in RNA error mediated by SF1B3 in cancer (Zhang et al., 2022). DHX15 forms aggregates with NLRP6 inflammasomes and dsRNA, which participates in host defenses (Shen C. et al., 2021). DHX15 regulates natural killer cell differentiation through the IL-15 signaling pathway, which depends on the ATPase activity of DHX15 rather than its RNA helicase activity. It is a crucial regulatory factor for the homeostasis and function of natural killer cells (Wang et al., 2022). Furthermore, DHX15 is essential in the context of mild/severe COVID-19 and its association with hepatocellular carcinoma and chronic hepatitis B (Sokouti, 2022). However, its ATPase activity and involvement in pre-mRNA splicing have been implicated in pancreatic beta-cell failure (Cnop et al., 2017; Herchuelz and Pachera, 2018). DHX15 may contribute to the progression of diabetes by affecting pancreatic beta cells.
Ubiquitin-specific protease 14 (USP14) is a deubiquitinating enzyme and a member of the ubiquitin-specific protease family. USP14 prevents protein substrate degradation by blocking the ubiquitin chain and promoting protein degradation by activating proteases (Wang F. et al., 2021). USP14 plays a role in the immune response by regulating signaling molecules involved in immunity (Wang F. et al., 2021; Wang D. et al., 2021). In addition to its role in the immune response, USP14 also plays an important role in viral infections, inflammation, tumors, autophagy, and neurodegenerative disorders (Wang D. et al., 2021; Li et al., 2021; Ming et al., 2022; Shi et al., 2022). USP14 is a potential therapeutic target for the treatment of α-herpesvirus-associated diseases. Lowering USP14 levels can inhibit viral replication, whereas agonizing USP14 restores viral replication (Ming et al., 2022). Inhibition of USP14 promotes autophagy in M1 macrophages, which reduces the severity of sepsis induced by cecal ligation and puncture (Xu et al., 2020). Thus, USP14 is closely related to metabolism, immunity, and inflammation. Research has shown that USP14 increases the levels and stability of 3′,5′-cyclic adenosine monophosphate response element-binding protein, enhancing the action of glucagon and hepatic glucose production, thereby promoting hyperglycemia (Liu et al., 2019). Additionally, the long non-coding RNA OGRU competes with miR-320 to regulate the expression of USP14, thereby mediating the progression of diabetic retinopathy (Fu et al., 2021). Furthermore, ubiquitination is crucial in virus infection by modifying viral proteins or host defense factors (Sun et al., 2019; Chen et al., 2024).
COP9 signalosome subunit 3 (COPS3) is an oncogene that promotes lung metastasis in osteosarcoma (Zhang et al., 2018). Autophagy induced by COPS3 promotes cisplatin resistance (Niu et al., 2023), and people with tumors with a low expression of COPS3 and Ras features have a better prognosis (Luo et al., 2009). COPS3 promotes the proliferation, invasion, and epithelial-to-mesenchymal transition of colorectal cancer cells through the MEK/ERK signaling pathway (Xie et al., 2022). It plays a crucial role in linking the Raf-1/MEK/ERK signaling pathway and autophagy regulation, inducing RAB7 to promote autophagosome and lysosome fusion (Zhang et al., 2018; Niu et al., 2023). The role of COPS3 in immunity is not yet clear. Research indicates that capsaicin and quercetin from chili peppers enhance hypoglycemic bioactivity through the RAS/Raf-1/MEK/ERK signaling pathway (Mi et al., 2023). Furthermore, studies suggest that autophagy plays a crucial role in maintaining the normal structure and function of pancreatic islets (Lee, 2014), and increased autophagy may delay the progression of T2DM and protect β-cell function (Yuan et al., 2017; Lambelet et al., 2018; Zhao et al., 2023a). Therefore, COPS3 may regulate the development of T2DM through the Raf-1/MEK/ERK signaling pathway and autophagy regulation. Additionally, SARS-CoV-2 can hinder autophagy by blocking autophagosome-lysosome fusion (Miao et al., 2021), and activation of autophagy can inhibit the replication and spread of SARS-CoV-2 (Gassen et al., 2021; Wei et al., 2021). Thus, COPS3 may regulate the development of COVID-19 by modulating autophagy.
Tyrosine kinase 2 (TYK2) is an intracellular kinase and a member of the Janus kinase family. TYK2 participates in the development of immune-related diseases by mediating cytokine signaling (Mease et al., 2022). The clinical efficacy of TYK2 inhibitors in inflammatory and autoimmune diseases has been demonstrated (Gonciarz et al., 2021), as evidenced by genomic, transcriptomic, and proteomic analyses detecting TYK2 mutations and expression changes in various tumors (Hammarén et al., 2019). TYK2 plays a crucial role in the progression of both diabetes and COVID-19. TYK2 regulates apoptosis in pancreatic islet β-cells and innate immune responses (Marroqui et al., 2015; Chandra et al., 2022). Interestingly, reduced expression of the TYK2 gene in pancreatic β-cells is a contributing factor to virus-induced susceptibility to diabetes (Izumi et al., 2015). The type I interferon response is one of three major pathways that affect susceptibility to and severity of COVID-19, and TYK2 (located at 19p13.2) contains candidate pathogenic genes associated with the type I interferon pathway (Host Genetics Initiative, 2022). The TYK2 locus shows a significant genetic association with severe COVID-19 (Pairo-Castineira et al., 2023). Defects in the TYK2 gene lead to a weakened antiviral response in the body (Shimoda et al., 2000), and downregulation of TYK2 is a molecular mechanism by which SARS-CoV-2 fails to induce a normal interferon response (Akbari et al., 2022).
RIOK2 belongs to the family of RIO kinases, which comprises three members: RIOK1, RIOK2, and RIOK3 (LaRonde-LeBlanc and Wlodawer, 2005). RIO kinases play a role in pre-ribosomal RNA processing and the biogenesis of ribosomes (Cerezo et al., 2021). The role of RIO kinases in cancer has been increasingly investigated, with the RIOK2 ATPase activity essential for cell survival. RIOK2 depletion leads to ribosome instability, decreased protein synthesis in leukemia cells, and apoptosis (Messling et al., 2022). The atypical kinase RIOK2 drives erythroid differentiation, inhibits megakaryocytic and myeloid formation, and plays an essential role in the transcriptional regulation of human hematopoietic differentiation. RIOK2 is closely associated with acute myeloid leukemia, myelodysplastic syndrome, and chronic kidney disease (Ghosh et al., 2022). RIO kinases regulate the cell cycle, AKT signaling, and activation of mutant RAS-driven tumor development (Asquith et al., 2019). Silencing RIOK2 inhibits the migration, invasion, and epithelial-mesenchymal transition of glioma cells, while overexpression of RIOK2 promotes these processes (Song et al., 2020). Currently, there is no research on the role of RIOK2 in T2DM. However, substance metabolism and biosynthesis are closely related to T2DM. Research reports suggest that RIOK is associated with host immune responses. RIOK3 inhibits the antiviral immune response by promoting degradation mediated by TRIM40 of RIG-I and MDA5 (Shen Y. et al., 2021; Zhao et al., 2023b). The immune response plays a crucial role in the progression of COVID-19 and T2DM. The roles of DHX15, USP14, COPS3, TYK2, and RIOK2 in T2DM and COVID-19 remain primarily unknown, emphasizing the importance of future research.
Functional networks of shared genes and their co-expressed genes mainly focus on protein modification by removing small proteins, nucleotide excision repair, and response to type I interferon during viral infection. TFs in the host (including interferon regulatory factors and nuclear factor kappa B) are activated, recruiting specific subpopulations of white blood cells (Hur, 2019; Rajpal et al., 2020). Hyperglycemia, on the other hand, inhibits the production of IFN-1 (Lazear et al., 2019; Rajpal et al., 2020). The impact on the host’s defense against viral infections is likely to involve metabolic pathways and immune responses, which may contribute to the development and progression of both T2DM and COVID-19.
Using the miRTarBase, Starbase, and Targetscan databases, we successfully predicted the targeted miRNAs of the common core genes. We predicted TFs targeted using the Enrichr database to construct a miRNA-TF-mRNA regulatory network. Previous studies have shown that many of the predicted miRNAs are involved in the development of inflammation, immunity, diabetes, and COVID-19. For example, miR-124-3p inhibits high glucose-induced endothelial cell dysfunction (Zhao and He, 2021), whereas miR-1-3p induces endothelial cell dysfunction through targeting SERP1 (Gao et al., 2021). MiR-221 plays a role in diabetes and COVID-19. MiR-221-3p promotes diabetic wound healing by targeting HIPK285 and regulates microvascular dysfunction in diabetic retinopathy by targeting TIMP3 (Wang C. et al., 2020). Dexmedetomidine protects against renal fibrosis in diabetic mice through targeting the miR-101-3p-mediated endothelial-mesenchymal transition (Song et al., 2022). MiR-204-5p regulates endoplasmic reticulum stress in human pancreatic beta cells (Grieco et al., 2019) and reduces renal interstitial cell injury induced by high glucose by blocking the AKT/NF-κB pathway (Qiu et al., 2022). Cyanidin-3-O-glucoside protects endothelial cells by regulating the miR-204-5p/SIRT1-mediated inflammatory response and apoptosis (Wang Z. et al., 2020). MiR-17-5p inhibits the TXNIP/NLRP3 inflammatory pathway and promotes glucose uptake in HTR8/SVneo trophoblast cells (Jiang et al., 2022). Downregulation of miR-17-5p alleviates apoptosis and fibrosis in human renal interstitial cells induced by high glucose (Chen et al., 2021). High glucose also inhibits miR-140-3p, impairing vascular endothelial cell angiogenesis in diabetes (Wang et al., 2019). MiR-215-5p regulates Doxorubicin-induced myocardial cell injury through targeting ZEB2 (Xiong et al., 2021). MiR-25-3p activates autophagy and improves high glucose-induced podocyte injury by inhibiting DUSP1 expression (Huang et al., 2020). Decreased expression of miR-192-5p is associated with diabetic nephropathy (Akpınar et al., 2022). MiR-221 is involved in the inflammatory immune response of severe COVID-19 patients (Gaytán-Pacheco et al., 2022). MiR-101-3p and miR-25-3p exhibit differential expression between severe/critical COVID-19 patients and patients with mild infections (Nicoletti et al., 2022). MiR-140-3p shows differential expression in lung tissues before and after COVID-19 infection (Kim et al., 2020). Significant differences in miR-192-5p expression exist between surviving and deceased patients with severe COVID-19 (de Gonzalo-Calvo et al., 2021). Additionally, miR-320 regulates USP14 expression and mediates the progression of diabetic retinopathy (Fu et al., 2021). MiR-26a-5p regulates the USP14/NF-κB pathway to alleviate inflammation and oxidative stress in diabetic retinopathy (Bian et al., 2024).
In summary, identifying common core genes and pathways in T2DM and COVID-19 may provide new insights into potential therapeutic targets for patients with both conditions. However, common biomarkers, pathways, and therapeutic drugs need to be confirmed through experiments. In this study, we identified common DEGs and common feature genes in T2DM and COVID-19 through bioinformatics analysis methods, predicted miRNAs, and TFs targeting commonly expressed genes.
Our findings suggest that T2DM and COVID-19 may share common pathogenic mechanisms mediated by specific hub genes. This study provides novel avenues for further molecular mechanism research.
Limitations of this study include the need for further experimental validation of the common core genes and pathways identified in T2DM and COVID-19. In addition, the diagnosis of T2DM and COVID-19 cannot be based solely on these common feature genes and pathways but requires consideration of symptoms and laboratory tests.
5 Conclusion
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 molecular mechanism studies.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://www.ncbi.nlm.nih.gov/geo/.
Author contributions
YW designed the study. YW and KL performed statistical analyses. YW wrote the first draft. YW, KL, SM, PY, JZ, SL, and SQ read and approved the final draft. All authors contributed to the article and approved the submitted version.
Acknowledgments
All the authors of the study would like to thank GEO for open access to the database.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The handling editor KW declared a shared parent affiliation with the author SQ at the time of review.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2024.1249501/full#supplementary-material
Abbreviations
T2DM, type 2 diabetes mellitus; COVID-19, coronavirus disease 2019; DEG, differentially expressed genes; GO, gene ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; PPI, protein-protein interaction; ROC, receiver operating characteristic; AUC, area under the curve; BP, biological process; MF, molecular function; CC, cellular component.
References
Agarwal, V., Bell, G. W., Nam, J. W., and Bartel, D. P. (2015). Predicting effective microRNA target sites in mammalian mRNAs. Elife 4, e05005. doi:10.7554/eLife.05005
Akbari, M., Akhavan-Bahabadi, M., Shafigh, N., Taheriazam, A., Hussen, B. M., Sayad, A., et al. (2022). Expression analysis of IFNAR1 and TYK2 transcripts in COVID-19 patients. Cytokine 153, 155849. doi:10.1016/j.cyto.2022.155849
Akpınar, K., Aslan, D., Fenkçi, S. M., and Caner, V. (2022). miR-21-3p and miR-192-5p in patients with type 2 diabetic nephropathy. Diagn. Berl. Ger. 9, 499–507. doi:10.1515/dx-2022-0036
Asquith, C. R. M., East, M. P., and Zuercher, W. J. (2019). RIOK2: straddling the kinase/ATPase line. Nat. Rev. Drug Discov. 18, 574. doi:10.1038/d41573-019-00107-7
Bian, J., Ge, W., and Jiang, Z. (2024). miR-26a-5p attenuates oxidative stress and inflammation in diabetic retinopathy through the USP14/NF-κB signaling pathway. J. Ophthalmol. 2024, 1470898. doi:10.1155/2024/1470898
Bo, Y., Yuli, C., Ye, W., Junfeng, L., Xiaolin, C., Yan, B., et al. (2022). Immune-inflammatory biomarkers and the risk of cardiac injury in COVID-19 patients with diabetes: a retrospective cohort study. Cardiovasc Diabetol. 21, 188. doi:10.1186/s12933-022-01625-2
Bohnsack, K. E., Ficner, R., Bohnsack, M. T., and Jonas, S. (2021). Regulation of DEAH-box RNA helicases by G-patch proteins. Biol. Chem. 402, 561–579. doi:10.1515/hsz-2020-0338
Bohnsack, K. E., Kanwal, N., and Bohnsack, M. T. (2022). Prp43/DHX15 exemplify RNA helicase multifunctionality in the gene expression network. Nucleic Acids Res. 50, 9012–9022. doi:10.1093/nar/gkac687
Cerezo, E. L., Houles, T., Lié, O., Sarthou, M. K., Audoynaud, C., Lavoie, G., et al. (2021). RIOK2 phosphorylation by RSK promotes synthesis of the human small ribosomal subunit. PLoS Genet. 17, e1009583. doi:10.1371/journal.pgen.1009583
Chandra, V., Ibrahim, H., Halliez, C., Prasad, R. B., Vecchio, F., Dwivedi, O. P., et al. (2022). The type 1 diabetes gene TYK2 regulates β-cell development and its responses to interferon-α. Nat. Commun. 13, 6363. doi:10.1038/s41467-022-34069-z
Chen, J., Liu, Q., He, J., and Li, Y. (2022). Immune responses in diabetic nephropathy: pathogenic mechanisms and therapeutic target. Front. Immunol. 13, 958790. doi:10.3389/fimmu.2022.958790
Chen, X., Gu, L., Cheng, X., Xing, J., and Zhang, M. (2021). MiR-17-5p downregulation alleviates apoptosis and fibrosis in high glucose-induced human mesangial cells through inactivation of Wnt/β-catenin signaling by targeting KIF23. Environ. Toxicol. 36, 1702–1712. doi:10.1002/tox.23280
Chen, X., Tian, L., Zhang, L., Gao, W., Yu, M., Li, Z., et al. (2024). Deubiquitinase USP39 promotes SARS-CoV-2 replication by deubiquitinating and stabilizing the envelope protein. Antivir. Res. 221, 105790. doi:10.1016/j.antiviral.2023.105790
Chou, C. H., Shrestha, S., Yang, C. D., Chang, N. W., Lin, Y. L., Liao, K. W., et al. (2018). miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 46, D296-D302–D302. doi:10.1093/nar/gkx1067
Cnop, M., Toivonen, S., Igoillo-Esteve, M., and Salpea, P. (2017). Endoplasmic reticulum stress and eIF2α phosphorylation: the Achilles heel of pancreatic β cells. Mol. Metab. 6, 1024–1039. doi:10.1016/j.molmet.2017.06.001
de Gonzalo-Calvo, D., Benítez, I. D., Pinilla, L., Carratalá, A., Moncusí-Moix, A., Gort-Paniello, C., et al. (2021). Circulating microRNA profiles predict the severity of COVID-19 in hospitalized patients. Transl. Res. J. laboratory Clin. Med. 236, 147–159. doi:10.1016/j.trsl.2021.05.004
Detanico, T., Virgen-Slane, R., Steen-Fuentes, S., Lin, W. W., Rhode-Kurnow, A., Chappell, E., et al. (2019). Co-Expression networks identify DHX15 RNA helicase as a B cell regulatory factor. Front. Immunol. 10, 2903. doi:10.3389/fimmu.2019.02903
Farnsworth, C. W., Schott, E. M., Benvie, A., Kates, S. L., Schwarz, E. M., Gill, S. R., et al. (2018). Exacerbated Staphylococcus aureus foot infections in obese/diabetic mice are associated with impaired germinal center reactions, ig class switching, and humoral immunity. J. Immunol. 201 (201), 560–572. doi:10.4049/jimmunol.1800253
Fu, S., Zheng, Y., Sun, Y., Lai, M., Qiu, J., Gui, F., et al. (2021). Suppressing long noncoding RNA OGRU ameliorates diabetic retinopathy by inhibition of oxidative stress and inflammation via miR-320/USP14 axis. Free Radic. Biol. Med. 169, 361–381. doi:10.1016/j.freeradbiomed.2021.03.016
Gao, M., Yu, T., Liu, D., Shi, Y., Yang, P., Zhang, J., et al. (2021). Sepsis plasma-derived exosomal miR-1-3p induces endothelial cell dysfunction by targeting SERP1. Clin. Sci. 135, 347–365. doi:10.1042/cs20200573
Gassen, N. C., Papies, J., Bajaj, T., Emanuel, J., Dethloff, F., Chua, R. L., et al. (2021). SARS-CoV-2-mediated dysregulation of metabolism and autophagy uncovers host-targeting antivirals. Nat. Commun. 12, 3818. doi:10.1038/s41467-021-24007-w
Gaytán-Pacheco, N., Ibáñez-Salazar, A., Herrera-Van Oostdam, A. S., Oropeza-Valdez, J. J., Magaña-Aquino, M., Adrián López, J., et al. (2022). miR-146a, miR-221, and miR-155 are involved in inflammatory immune response in severe COVID-19 patients. Diagn. Basel. Switz. 13, 133. doi:10.3390/diagnostics13010133
Ghosh, S., Raundhal, M., Myers, S. A., Carr, S. A., Chen, X., Petsko, G. A., et al. (2022). Identification of RIOK2 as a master regulator of human blood cell development. Nat. Immunol. 23, 109–121. doi:10.1038/s41590-021-01079-w
Gonciarz, M., Pawlak-Buś, K., Leszczyński, P., and Owczarek, W. (2021). TYK2 as a therapeutic target in the treatment of autoimmune and inflammatory diseases. Immunotherapy 13, 1135–1150. doi:10.2217/imt-2021-0096
Grewal, A. S., Lather, V., Charaya, N., Sharma, N., Singh, S., and Kairys, V. (2020). Recent developments in medicinal chemistry of allosteric activators of human glucokinase for type 2 diabetes mellitus therapeutics. Curr. Pharm. Des. 26, 2510–2552. doi:10.2174/1381612826666200414163148
Grieco, F. A., Schiavo, A. A., Brozzi, F., Juan-Mateu, J., Bugliani, M., Marchetti, P., et al. (2019). The miRNAs miR-211-5p and miR-204-5p modulate ER stress in human beta cells. J. Mol. Endocrinol. 63, 139–149. doi:10.1530/jme-19-0066
Grifoni, A., Weiskopf, D., Ramirez, S. I., Mateus, J., Dan, J. M., Moderbacher, C. R., et al. (2020). Targets of T Cell responses to SARS-CoV-2 coronavirus in humans with COVID-19 disease and unexposed individuals. Cell 181, 1489–1501. doi:10.1016/j.cell.2020.05.015
Guo, W., Li, M., Dong, Y., Zhou, H., Zhang, Z., Tian, C., et al. (2020). Diabetes is a risk factor for the progression and prognosis of COVID-19. Diabetes Metab. Res. Rev. 36, e3319. doi:10.1002/dmrr.3319
Hammarén, H. M., Virtanen, A. T., Raivola, J., and Silvennoinen, O. (2019). The regulation of JAKs in cytokine signaling and its breakdown in disease. Cytokine 118, 48–63. doi:10.1016/j.cyto.2018.03.041
Herchuelz, A., and Pachera, N. (2018). The Na+/Ca2+ exchanger and the Plasma Membrane Ca2+-ATPase in β-cell function and diabetes. Neurosci. Lett. 663, 72–78. doi:10.1016/j.neulet.2017.08.009
Host Genetics Initiative (2022). A first update on mapping the human genetic architecture of COVID-19. Nature 608, E1–e10. doi:10.1038/s41586-022-04826-7
Huang, H., Liu, H., Tang, J., Xu, W., Gan, H., Fan, Q., et al. (2020). M2 macrophage-derived exosomal miR-25-3p improves high glucose-induced podocytes injury through activation autophagy via inhibiting DUSP1 expression. IUBMB life 72, 2651–2662. doi:10.1002/iub.2393
Hur, S. (2019). Double-stranded RNA sensors and modulators in innate immunity. Annu. Rev. Immunol. 37, 349–375. doi:10.1146/annurev-immunol-042718-041356
Izumi, K., Mine, K., Inoue, Y., Teshima, M., Ogawa, S., Kai, Y., et al. (2015). Reduced Tyk2 gene expression in β-cells due to natural mutation determines susceptibility to virus-induced diabetes. Nat. Commun. 6, 6748. doi:10.1038/ncomms7748
Jiang, Y., Wei, L., Zhang, H., Chen, Y., Gao, P., Zhang, J., et al. (2022). miR-17-5p promotes glucose uptake of HTR8/SVneo trophoblast cells by inhibiting TXNIP/NLRP3 inflammasome pathway. Diabetes, metabolic syndrome Obes. targets Ther. 15, 3361–3374. doi:10.2147/dmso.S385774
Khunti, K., Del Prato, S., Mathieu, C., Kahn, S. E., Gabbay, R. A., and Buse, J. B. (2021). COVID-19, hyperglycemia, and new-onset diabetes. Diabetes care 44, 2645–2655. doi:10.2337/dc21-1318
Kim, W. R., Park, E. G., Kang, K. W., Lee, S. M., Kim, B., and Kim, H. S. (2020). Expression analyses of MicroRNAs in hamster lung tissues infected by SARS-CoV-2. Mol. cells 43, 953–963. doi:10.14348/molcells.2020.0177
Lambelet, M., Terra, L. F., Fukaya, M., Meyerovich, K., Labriola, L., Cardozo, A. K., et al. (2018). Dysfunctional autophagy following exposure to pro-inflammatory cytokines contributes to pancreatic β-cell apoptosis. Cell Death Dis. 9, 96. doi:10.1038/s41419-017-0121-5
LaRonde-LeBlanc, N., and Wlodawer, A. (2005). A family portrait of the RIO kinases. J. Biol. Chem. 280, 37297–37300. doi:10.1074/jbc.R500013200
Lazear, H. M., Schoggins, J. W., and Diamond, M. S. (2019). Shared and distinct functions of type I and type III interferons. Immunity 50, 907–923. doi:10.1016/j.immuni.2019.03.025
Lee, M. S. (2014). Role of islet β cell autophagy in the pathogenesis of diabetes. Trends Endocrinol. metabolism TEM 25, 620–627. doi:10.1016/j.tem.2014.08.005
Leon-Abarca, J. A., Portmann-Baracco, A., Bryce-Alberti, M., Ruiz-Sánchez, C., Accinelli, R. A., Soliz, J., et al. (2021). Diabetes increases the risk of COVID-19 in an altitude dependent manner: an analysis of 1,280,806 Mexican patients. PLoS One 16, e0255144. doi:10.1371/journal.pone.0255144
Li, C., Sun, G., Chen, B., Xu, L., Ye, Y., He, J., et al. (2021). Nuclear receptor coactivator 4-mediated ferritinophagy contributes to cerebral ischemia-induced ferroptosis in ischemic stroke. Pharmacol. Res. 174, 105933. doi:10.1016/j.phrs.2021.105933
Liu, B., Zhang, Z., Hu, Y., Lu, Y., Li, D., Liu, J., et al. (2019). Sustained ER stress promotes hyperglycemia by increasing glucagon action through the deubiquitinating enzyme USP14. Proc. Natl. Acad. Sci. U. S. A. 116, 21732–21738. doi:10.1073/pnas.1907288116
Luo, J., Emanuele, M. J., Li, D., Creighton, C. J., Schlabach, M. R., Westbrook, T. F., et al. (2009). A genome-wide RNAi screen identifies multiple synthetic lethal interactions with the Ras oncogene. Cell 137, 835–848. doi:10.1016/j.cell.2009.05.006
Marqués, J. M., Rial, A., Muñoz, N., Pellay, F. X., Van Maele, L., Léger, H., et al. (2012). Protection against Streptococcus pneumoniae serotype 1 acute infection shows a signature of Th17- and IFN-γ-mediated immunity. Immunobiology 217, 420–429. doi:10.1016/j.imbio.2011.10.012
Marroqui, L., Dos Santos, R. S., Fløyel, T., Grieco, F. A., Santin, I., Op de Beeck, A., et al. (2015). TYK2, a candidate gene for type 1 diabetes, modulates apoptosis and the innate immune response in human pancreatic β-cells. Diabetes 64, 3808–3817. doi:10.2337/db15-0362
Mease, P. J., Deodhar, A. A., van der Heijde, D., Behrens, F., Kivitz, A. J., Neal, J., et al. (2022). Efficacy and safety of selective TYK2 inhibitor, deucravacitinib, in a phase II trial in psoriatic arthritis. Ann. rheumatic Dis. 81, 815–822. doi:10.1136/annrheumdis-2021-221664
Messling, J. E., Agger, K., Andersen, K. L., Kromer, K., Kuepper, H. M., Lund, A. H., et al. (2022). Targeting RIOK2 ATPase activity leads to decreased protein synthesis and cell death in acute myeloid leukemia. Blood 139, 245–255. doi:10.1182/blood.2021012629
Mi, S., Zhu, W., Zhang, X., Wang, Y., Li, T., and Wang, X. (2023). Enhanced hypoglycemic bioactivity via RAS/Raf-1/MEK/ERK signaling pathway by combining capsaicin and QUERCETIN from chili peppers. Mol. Nutr. food Res. 67, e2200577. doi:10.1002/mnfr.202200577
Miao, G., Zhao, H., Li, Y., Ji, M., Chen, Y., Shi, Y., et al. (2021). ORF3a of the COVID-19 virus SARS-CoV-2 blocks HOPS complex-mediated assembly of the SNARE complex required for autolysosome formation. Dev. Cell 56, 427–442.e5. doi:10.1016/j.devcel.2020.12.010
Ming, S. L., Zhang, S., Wang, Q., Zeng, L., Zhou, L. Y., Wang, M. D., et al. (2022). Inhibition of USP14 influences alphaherpesvirus proliferation by degrading viral VP16 protein via ER stress-triggered selective autophagy. Autophagy 18, 1801–1821. doi:10.1080/15548627.2021.2002101
Mosallanejad, K., Sekine, Y., Ishikura-Kinoshita, S., Kumagai, K., Nagano, T., Matsuzawa, A., et al. (2014). The DEAH-box RNA helicase DHX15 activates NF-κB and MAPK signaling downstream of MAVS during antiviral responses. Sci. Signal 7, ra40. doi:10.1126/scisignal.2004841
Moura, J., Rodrigues, J., Gonçalves, M., Amaral, C., Lima, M., and Carvalho, E. (2017). Impaired T-cell differentiation in diabetic foot ulceration. Cell Mol. Immunol. 14, 758–769. doi:10.1038/cmi.2015.116
Muniyappa, R., and Gubbi, S. (2020). COVID-19 pandemic, coronaviruses, and diabetes mellitus. Am. J. Physiol. Endocrinol. Metab. 318, E736-E741–E741. doi:10.1152/ajpendo.00124.2020
Ni, J., Qiu, L. J., Yin, K. J., Chen, G. M., and Pan, H. F. (2023). Shared genetic architecture between type 2 diabetes and COVID-19 severity. J. Endocrinol. investigation 46, 501–507. doi:10.1007/s40618-022-01920-5
Nicoletti, A. S., Visacri, M. B., da Ronda, C. R. d. S. C., Vasconcelos, P. E. d. N. S., Quintanilha, J. C. F., de Souza, R. N., et al. (2022). Differentially expressed plasmatic microRNAs in Brazilian patients with Coronavirus disease 2019 (COVID-19): preliminary results. Mol. Biol. Rep. 49, 6931–6943. doi:10.1007/s11033-022-07338-9
Niu, J., Yan, T., Guo, W., Wang, W., Ren, T., Huang, Y., et al. (2023). The COPS3-FOXO3 positive feedback loop regulates autophagy to promote cisplatin resistance in osteosarcoma. Autophagy 19, 1693–1710. doi:10.1080/15548627.2022.2150003
Nyland, J. E., Raja-Khan, N. T., Bettermann, K., Haouzi, P. A., Leslie, D. L., Kraschnewski, J. L., et al. (2021). Diabetes, drug treatment, and mortality in COVID-19: a multinational retrospective cohort study. Diabetes 70, 2903–2916. doi:10.2337/db21-0385
Pairo-Castineira, E., Rawlik, K., Bretherick, A. D., Qi, T., Wu, Y., Nassiri, I., et al. (2023). GWAS and meta-analysis identifies 49 genetic variants underlying critical COVID-19. Nature 617, 764–768. doi:10.1038/s41586-023-06034-3
Qiu, B., Qi, X., and Wang, J. (2022). CircTLK1 downregulation attenuates high glucose-induced human mesangial cell injury by blocking the AKT/NF-κB pathway through sponging miR-126-5p/miR-204-5p. Biochem. Genet. 60, 1471–1487. doi:10.1007/s10528-021-10146-8
Rajpal, A., Rahimi, L., and Ismail-Beigi, F. (2020). Factors leading to high morbidity and mortality of COVID-19 in patients with type 2 diabetes. J. Diabetes 12, 895–908. doi:10.1111/1753-0407.13085
Rayego-Mateos, S., Rodrigues-Diez, R. R., Fernandez-Fernandez, B., Mora-Fernández, C., Marchant, V., Donate-Correa, J., et al. (2023). Targeting inflammation to treat diabetic kidney disease: the road to 2030. Kidney Int. 103, 282–296. doi:10.1016/j.kint.2022.10.030
Richard, C., Wadowski, M., Goruk, S., Cameron, L., Sharma, A. M., and Field, C. J. (2017). Individuals with obesity and type 2 diabetes have additional immune dysfunction compared with obese individuals who are metabolically healthy. BMJ open diabetes Res. care 5, e000379. doi:10.1136/bmjdrc-2016-000379
Sardu, C., Gargiulo, G., Esposito, G., Paolisso, G., and Marfella, R. (2020). Impact of diabetes mellitus on clinical outcomes in patients affected by Covid-19. Cardiovasc Diabetol. 19, 76. doi:10.1186/s12933-020-01047-y
Semlow, D. R., Blanco, M. R., Walter, N. G., and Staley, J. P. (2016). Spliceosomal DEAH-box ATPases remodel pre-mRNA to activate alternative splice sites. Cell 164, 985–998. doi:10.1016/j.cell.2016.01.025
Sharma, P., Behl, T., Sharma, N., Singh, S., Grewal, A. S., Albarrati, A., et al. (2022). COVID-19 and diabetes: association intensify risk factors for morbidity and mortality. Biomed. Pharmacother. 151, 113089. doi:10.1016/j.biopha.2022.113089
Shaw, J. E., Sicree, R. A., and Zimmet, P. Z. (2010). Global estimates of the prevalence of diabetes for 2010 and 2030. Diabetes Res. Clin. Pract. 87, 4–14. doi:10.1016/j.diabres.2009.10.007
Shen, C., Li, R., Negro, R., Cheng, J., Vora, S. M., Fu, T. M., et al. (2021a). Phase separation drives RNA virus-induced activation of the NLRP6 inflammasome. Cell 184, 5759–5774.e20. doi:10.1016/j.cell.2021.09.032
Shen, Y., Tang, K., Chen, D., Hong, M., Sun, F., Wang, S., et al. (2021b). Riok3 inhibits the antiviral immune response by facilitating TRIM40-mediated RIG-I and MDA5 degradation. Cell Rep. 35, 109272. doi:10.1016/j.celrep.2021.109272
Shi, D., Wu, X., Jian, Y., Wang, J., Huang, C., Mo, S., et al. (2022). USP14 promotes tryptophan metabolism and immune suppression by stabilizing Ido1 in colorectal cancer. Nat. Commun. 13, 5644. doi:10.1038/s41467-022-33285-x
Shimoda, K., Kato, K., Aoki, K., Matsuda, T., Miyamoto, A., Shibamori, M., et al. (2000). Tyk2 plays a restricted role in IFN alpha signaling, although it is required for IL-12-mediated T cell function. Immunity 13, 561–571. doi:10.1016/s1074-7613(00)00055-8
Sokouti, B. (2022). A systems biology approach for investigating significantly expressed genes among COVID-19, hepatocellular carcinoma, and chronic hepatitis B. Egypt. J. Med. Hum. Genet. 23, 146. doi:10.1186/s43042-022-00360-3
Song, L., Feng, S., Yu, H., and Shi, S. (2022). Dexmedetomidine protects against kidney fibrosis in diabetic mice by targeting miR-101-3p-mediated EndMT. Dose-response a Publ. Int. Hormesis Soc. 20, 15593258221083486. doi:10.1177/15593258221083486
Song, Y., Li, C., Jin, L., Xing, J., Sha, Z., Zhang, T., et al. (2020). RIOK2 is negatively regulated by miR-4744 and promotes glioma cell migration/invasion through epithelial-mesenchymal transition. J. Cell Mol. Med. 24, 4494–4509. doi:10.1111/jcmm.15107
Sun, Y., Zheng, H., Yu, S., Ding, Y., Wu, W., Mao, X., et al. (2019). Newcastle disease virus V protein degrades mitochondrial antiviral signaling protein to inhibit host type I interferon production via E3 ubiquitin ligase RNF5. J. virology 93, e00322-19. doi:10.1128/jvi.00322-19
To, K. K., Sridhar, S., Chiu, K. H. Y., Hung, D. L. L., Li, X., Hung, I. F. N., et al. (2021). Lessons learned 1 year after SARS-CoV-2 emergence leading to COVID-19 pandemic. Emerg. microbes Infect. 10, 507–535. doi:10.1080/22221751.2021.1898291
Tong, Z. W. M., Grant, E., Gras, S., Wu, M., Smith, C., Barrett, H. L., et al. (2021). The role of T-cell immunity in COVID-19 severity amongst people living with type II diabetes. FEBS J. 288, 5042–5054. doi:10.1111/febs.16105
Wang, C., Lin, Y., Fu, Y., Zhang, D., and Xin, Y. (2020b). MiR-221-3p regulates the microvascular dysfunction in diabetic retinopathy by targeting TIMP3. Pflugers Archiv Eur. J. physiology 472, 1607–1618. doi:10.1007/s00424-020-02432-y
Wang, D., Ma, H., Zhao, Y., and Zhao, J. (2021b). Ubiquitin-specific protease 14 is a new therapeutic target for the treatment of diseases. J. Cell Physiol. 236, 3396–3405. doi:10.1002/jcp.30124
Wang, D., Wang, H., Liu, C., Mu, X., and Cheng, S. (2019). Hyperglycemia inhibition of endothelial miR-140-3p mediates angiogenic dysfunction in diabetes mellitus. J. diabetes its Complicat. 33, 374–382. doi:10.1016/j.jdiacomp.2019.02.001
Wang, F., Ning, S., Yu, B., and Wang, Y. (2021a). USP14: structure, function, and target inhibition. Front. Pharmacol. 12, 801328. doi:10.3389/fphar.2021.801328
Wang, G., Xiao, X., Wang, Y., Chu, X., Dou, Y., Minze, L. J., et al. (2022). The RNA helicase DHX15 is a critical regulator of natural killer-cell homeostasis and functions. Cell Mol. Immunol. 19, 687–701. doi:10.1038/s41423-022-00852-7
Wang, W., Lu, J., Gu, W., Zhang, Y., Liu, J., and Ning, G. (2020a). Care for diabetes with COVID-19: advice from China. J. Diabetes 12, 417–419. doi:10.1111/1753-0407.13036
Wang, Z., Zhang, M., Guo, Z., and Chen, Q. (2020c). Cyanidin-3-O-glucoside attenuates endothelial cell dysfunction by modulating miR-204-5p/SIRT1-mediated inflammation and apoptosis. BioFactors Oxf. Engl. 46, 803–812. doi:10.1002/biof.1660
Wei, Z., Gao, Y., Meng, F., Chen, X., Gong, Y., Zhu, C., et al. (2021). iDMer: an integrative and mechanism-driven response system for identifying compound interventions for sudden virus outbreak. Briefings Bioinforma. 22, 976–987. doi:10.1093/bib/bbaa341
Wronka, M., Krzeminska, J., Mlynarska, E., Rysz, J., and Franczyk, B. (2022). The influence of lifestyle and treatment on oxidative stress and inflammation in diabetes. Int. J. Mol. Sci. 23, 15743. doi:10.3390/ijms232415743
Wu, J., Li, C., Lei, Z., Cai, H., Hu, Y., Zhu, Y., et al. (2022b). Comprehensive analysis of circRNA-miRNA-mRNA regulatory network and novel potential biomarkers in acute myocardial infarction. Front. Cardiovasc. Med. 9, 850991. doi:10.3389/fcvm.2022.850991
Wu, K. C. H., He, Q., Bennett, A. N., Li, J., and Chan, K. H. K. (2022a). Shared genetic mechanism between type 2 diabetes and COVID-19 using pathway-based association analysis. Front. Genet. 13, 1063519. doi:10.3389/fgene.2022.1063519
Xie, Y., Wei, Z., and Cheng, C. (2022). COPS3 promotes proliferation, invasion, and EMT of colorectal cancer cells by MEK/ERK pathway. Dis. markers 2022, 7594489. doi:10.1155/2022/7594489
Xiong, X., He, Q., Liu, J., Dai, R., Zhang, H., Cao, Z., et al. (2021). MicroRNA miR-215-5p regulates doxorubicin-induced cardiomyocyte injury by targeting ZEB2. J. Cardiovasc. Pharmacol. 78, 622–629. doi:10.1097/fjc.0000000000001110
Xu, F., Ma, Y., Huang, W., Gao, J., Guo, M., Li, J., et al. (2020). Typically inhibiting USP14 promotes autophagy in M1-like macrophages and alleviates CLP-induced sepsis. Cell Death Dis. 11, 666. doi:10.1038/s41419-020-02898-9
Yang, J. H., Li, J. H., Shao, P., Zhou, H., Chen, Y. Q., and Qu, L. H. (2011). starBase: a database for exploring microRNA-mRNA interaction maps from Argonaute CLIP-Seq and Degradome-Seq data. Nucleic Acids Res. 39, D202–D209. doi:10.1093/nar/gkq1056
Yu, Q., Liu, L., Zhang, X., Chang, H., Ma, S., Xie, Z., et al. (2022). MiR-221-3p targets HIPK2 to promote diabetic wound healing. Microvasc. Res. 140, 104306. doi:10.1016/j.mvr.2021.104306
Yuan, T., Rafizadeh, S., Gorrepati, K. D. D., Lupse, B., Oberholzer, J., Maedler, K., et al. (2017). Reciprocal regulation of mTOR complexes in pancreatic islets from humans with type 2 diabetes. Diabetologia 60, 668–678. doi:10.1007/s00125-016-4188-9
Zhang, F., Yan, T., Guo, W., Sun, K., Wang, S., Bao, X., et al. (2018). Novel oncogene COPS3 interacts with Beclin1 and Raf-1 to regulate metastasis of osteosarcoma through autophagy. J. Exp. Clin. cancer Res. CR 37, 135. doi:10.1186/s13046-018-0791-6
Zhang, J., Huang, J., Xu, K., Xing, P., Huang, Y., Liu, Z., et al. (2022). DHX15 is involved in SUGP1-mediated RNA missplicing by mutant SF3B1 in cancer. Proc. Natl. Acad. Sci. U. S. A. 119, e2216712119. doi:10.1073/pnas.2216712119
Zhao, H., and He, Y. (2021). MiR-124-3p suppresses the dysfunction of high glucose-stimulated endothelial cells by targeting G3BP2. Front. Genet. 12, 723625. doi:10.3389/fgene.2021.723625
Zhao, X., Bie, L. Y., Pang, D. R., Li, X., Yang, L. F., Chen, D. D., et al. (2023a). The role of autophagy in the treatment of type II diabetes and its complications: a review. Front. Endocrinol. (Lausanne) 14, 1228045. doi:10.3389/fendo.2023.1228045
Zhao, X., Dan, C., Gong, X. Y., Li, Y. L., Qu, Z. L., Sun, H. Y., et al. (2023b). Yellow catfish RIO kinases (RIOKs) negatively regulate fish interferon-mediated antiviral response. Dev. Comp. Immunol. 142, 104656. doi:10.1016/j.dci.2023.104656
Keywords: T2DM, COVID-19, common feature genes, bioinformatics, pathways
Citation: Wang Y, Li K, Mo S, Yao P, Zeng J, Lu S and Qin S (2024) Identification of common genes and pathways between type 2 diabetes and COVID-19. Front. Genet. 15:1249501. doi: 10.3389/fgene.2024.1249501
Received: 28 June 2023; Accepted: 21 March 2024;
Published: 18 April 2024.
Edited by:
Ka Wu, Nanning Second People’s Hospital, ChinaReviewed by:
Alex Cleber Improta Caria, University of São Paulo, BrazilZhongshan Cheng, St. Jude Children’s Research Hospital, United States
Copyright © 2024 Wang, Li, Mo, Yao, Zeng, Lu and Qin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shunyu Lu, nnlsy2010@163.com; Shanyu Qin, qinshanyu@gxmu.edu.cn