key: cord-0816510-7t8yv5b8 authors: Bhavya,; Pathak, Ekta; Mishra, Rajeev title: Linking Diabetes mellitus to SARS-CoV-2 infection through differential targeting of the microRNAs in the Pancreas tissue date: 2021-03-31 journal: bioRxiv DOI: 10.1101/2021.03.31.437823 sha: b964524daa84b22f90f78a3dcc2c09eb104ac5a3 doc_id: 816510 cord_uid: 7t8yv5b8 Coronavirus Disease 2019 (COVID-19) severity and Diabetes mellitus affect each other bidirectionally. The plus-sense single-stranded RNA (+ssRNA) genome of the SARS-CoV-2 virus can be targeted and suppressed by the host cell’s microRNAs (miRNAs). Using the differential gene expression analysis between the mock-infected and the SARS-CoV-2-infected pancreatic tissue, we report five Diabetes-associated genes that are upregulated due to SARS-CoV-2 infection in the hESC pancreas tissues. Ten miRNAs regulating these five genes can potentially target the SARS-CoV-2 genome. We hypothesize that the SARS-CoV-2 genome copies in the infected human pancreas cell compete with the host cell’s native genes in being regulated by the native miRNAs. It leads to the reduced miRNA-regulation and, thus, the upregulation of the Diabetes-associated native genes. Thus, the resultant new-onset or elevated Diabetic symptoms may worsen the condition of COVID-19 patients. Coronavirus Disease (COVID-19) is a fast-spreading disease that has caused a global crisis. This pandemic is a highly infectious viral disease caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), also known as novel coronavirus-2019 [1, 2] . According to the World Health Organization's (WHO's) COVID-19 Dashboard [3] , 110,749,023 cases and 2,455,131 deaths have been reported globally due to COVID-19 as of 21 st February 2021. SARS-CoV-2 has a positive-sense single-stranded ribonucleic acid (+ssRNA) genome enclosed in a protein-containing lipid bilayer [4, 5] . Its RNA genome acts as a messenger RNA (mRNA) after entering the host cell and is directly translated by the host cell's ribosomes, resulting in viral proteins. These translated viral proteins are also responsible for viral RNA replication [6] . Patients having either type 1 or type 2 Diabetes mellitus (commonly referred to as Diabetes) show poor prognosis with SARS-CoV-2 infection due to blood glucose level fluctuations and metabolic complications. An increase in blood glucose level helps in the SARS-CoV-2 replication and proliferation in the human monocytes [5, 7] . New-onset Diabetes has also been observed after COVID-19 infection [8] . Thus, COVID-19 severity and Diabetes affect each other in a bidirectional manner. Moreover, the resultant hyperglycaemia can diminish the individual's immune response towards handling the viral infection. COVID-19 mortality is also amplified for diabetic patients due to complications like cardiovascular diseases and kidney-related problems [8] [9] [10] [11] . Thus, regulating the blood glucose levels and preventing diabetic complications is required for diabetic COVID-19 patients to avert severe consequences of the infection on them [5] . Autopsy reports have been known to reveal COVID-19 infection in the human pancreas [12] . It becomes possible due to the angiotensin-converting enzyme 2 (ACE2) receptors on the pancreas cells. ACE2 acts as a very high-affinity receptor on the cell for the spike protein of SARS-CoV-2 [13, 14] . MicroRNAs (miRNAs) play an essential role in the regulation of gene expression in a cell. In a virus-infected cell, the host cell's miRNAs help in cell defense by targeting the viral RNA genome or its transcribed RNAs. In humans, the miRNAs mostly target the 3'UTR of the mRNAs. However, the miRNAs primarily target the 3'UTR and 5'UTR of the viral RNA genomes. During a viral infection, the host cell's miRNAs may target the viral genome rather than the native mRNAs. It leads to a competition for miRNA regulation between the host cell's mRNAs and the viral genome copies in the cell. This differential targeting of the miRNAs may result in the dysregulation of the host cell's genes [15] [16] [17] [18] . In this study, we hypothesize the link between Diabetes and COVID-19 through differential targeting of the miRNAs in the Pancreas tissue. The disease enrichment based on differential gene expression analysis revealed pancreas-and diabetes-associated terms. We identified the miRNAs that can potentially target the SARS-CoV-2 genome as well as the Diabetesassociated genes upregulated in SARS-CoV-2-infected pancreas tissue. Differential gene expression analysis -The RNA-Seq gene expression data for mockinfected and SARS-CoV-2-infected hESC pancreatic tissue were retrieved from NCBI. For the study, three mock-infected and three SARS-CoV-2-infected hESC pancreatic tissue samples were chosen from the Gene Expression Omnibus (GEO) dataset accession, GSE151803 [19] . The Differentially Expressed Genes (DEGs) between the mock-infected and the SARS-CoV-2-infected pancreatic tissue were obtained using the DESeq2 R package (padj < 0.05 and |log2FC| > 1) [20, 21] . Disease enrichment analysis -DEGs-based disease enrichment analysis was done through DAVID -Functional Annotation Tool considering the Gene-Disease Associations Dataset (GAD) [22, 23] . The pancreas-associated diseases and the DEGs associated with them were selected for further studies. A literature search was done for the obtained DEGs to understand their role in the respective diseases. The gene-gene interaction (GGI) network of the DEGs was created using the Cytoscape-GeneMANIA app [24] [25] [26] . Network topological analysis was done using the Cytoscape-NetworkAnalyzer plugin [27] to get an insight about the influence of the genes in the network. The target genes of the CoV-tar-miRNAs were compared to the enriched Diabetes-associated genes. The genes common in both were those Diabetes-associated genes that may get upregulated in pancreas cells due to miRNAs' differential targeting during SARS-CoV-2 infection. Furthermore, the availability of miRNAs in the pancreas tissue was checked from the web-based repository, Tissue-Atlas [33]. pancreatic tissue -Using the DESeq2 R package, 30 differentially expressed genes (DEGs) between the mock-infected and SARS-CoV-2-infected pancreatic tissues of the GEO dataset GSE151803 were revealed [19] [20] [21] . 26 DEGs were upregulated and 4 DEGs were downregulated (Table 1) in the SARS-CoV-2-infected hESC pancreas tissue. Disease enrichment of DEGs -DAVID Tool was used for DEGs-based disease enrichment analysis based on the Gene-Disease Associations Dataset (GAD) [22, 23] . Gene-based disease enrichment of the 30 DEGs resulted in Type I Diabetes as the most significant GAD disease. It was linked to four upregulated DEGs, i.e., CP, SOCS3, AGT and PSMB8 (Table 2) . Two upregulated DEGs, CP and CFB were enriched for the term "insulin." has been linked with the other resultant disease terms, but only the terms associated with the pancreas were selected for the study. network of the DEGs contains 29 genes as HCP5 is a lncRNA, and its data is not present in GeneMANIA. Diabetes-associated (Type 1 Diabetes and insulin-associated) DEGs are connected directly or indirectly through co-expression ( Figure 1 , Table 3 ). Two downregulated DEGs, i.e., AKR1B10 and PKHD1L1, show interaction with the upregulated DEGs ( Figure 1 , Table 4 ). After removing the duplicate edges between the gene nodes, the network's topological parameters were calculated (Table 5 ). Closeness centralities of the SUCNR1 and TACR1 genes were ignored due to their interaction with each other only and not with the significant network. The Diabetes-associated gene CP has the highest closeness and betweenness centrality; thus, it is the most influential gene in the network. The CFB gene is the secondhighest diabetes-associated gene in terms of closeness centrality. 32], 10 and 11 miRNAs, i.e., CoV-tar-miRNAs, were found to be potentially targeting the 3'UTR (Table 6 ) and 5'UTR (Table 7) of the SARS-CoV-2 genome, respectively. CoV-tar-miRNAs targeting Diabetes-associated DEGs -Each of the five Diabetesassociated DEGs was found to be targeted by one or more CoV-tar-miRNAs. Eight 3'UTR and nine 5'UTR CoV-tar-miRNAs were found to be potentially targeting the Diabetesassociated DEGs, respectively. Availability and raw microarray expression of five 3'UTR and six 5'UTR CoV-tar-miRNAs in the human pancreas were obtained from the TissueAtlas webbased repository. Due to the microarray expression data's unavailability in TissueAtlas, hsa-miR-466, hsa-miR-3123, hsa-miR-4691-3p, hsa-miR-3941, hsa-miR-6749-3p and hsa-miR-7851-3p were not considered for further study. hsa-mir-4775 has a negative raw microarray expression value in the human pancreas; thus, it was also not considered for differential targeting analysis. Four 3'UTR and six 5'UTR CoV-tar-miRNAs are expressed in the human pancreas (Table 8 ). The SARS-CoV-2 infection has been detected in many different human organs, including the pancreas. The viral infection may lead to dysregulation of genes in the infected cells. 26 and 4 genes were found to be upregulated and downregulated, respectively, in the SARS-CoV-2- The human miRNAs mostly target the 3'UTR of the mRNAs in the cytoplasm of the cell. However, they primarily target the 3'UTR and 5'UTR of the infecting viral RNA genome. It has been seen that the SARS-CoV-2 virus also infects the pancreas, among other organs of the human body. Since it is a +ss-RNA genome virus, its genome itself acts as an mRNA. Thus, the pancreas cell's miRNAs must target the viral 3'UTR or 5'UTR, causing them to deflect from regulating the pancreas host cell's native genes. This may cause the upregulation of the host cell's native genes. Twenty-one human miRNAs (CoV-tar-miRNAs) were found to be potentially targeting the UTRs of the viral genome. Considering the availability of these miRNAs in the human pancreas tissue, 10 of them also target the Diabetes-associated genes, thus, keeping their expression in check before the infection. After infection, the SARS-CoV-2 genome competes with the mRNAs of these genes in being targeted and regulated by the miRNAs. This differential targeting of the miRNAs explains the upregulation of the Diabetes-associated genes after the viral infection ( Figure 2 ). The SOCS3 (Suppressor of Cytokine Signaling 3) gene codes for a protein that helps regulate cytokine signal transduction [37] . It is the potential target of four 3'UTR CoV-tar-miRNAs, i.e., hsa-miR-1236-3p, hsa-miR-4279, hsa-miR-4717-3p and hsa-miR-5088-5p, and two 5'UTR CoV-tar-miRNAs, i.e., hsa-miR-1303 and hsa-miR-3925-5p. Thus, SOCS3 upregulation after the viral infection may be due to the differential targeting of these miRNAs. Overexpression of SOCS3 has been observed in mice having Type 1 Diabetes. SOCS3-deficiency in pancreas beta cells is associated with increased resistance to apoptosis, thus, preventing Type 1 Diabetes [38] . However, in the SARS-CoV-2-infected Pancreas tissue, the expression of SOCS3 may be increased due to differential targeting of miRNAs; thus, it decreases the beta cells' resistance to the development of Diabetes. The PSMB8 gene is the potential target of three 3'UTR CoV-tar-miRNAs, i.e., hsa-miR-1236-3p, hsa-miR-4279 and hsa-miR-5088-5p, and two 5'UTR CoV-tar-miRNAs, i.e., hsa-miR-3925-5p and hsa-miR-298. PSMB8 codes for Proteasome 20S Subunit Beta 8 and has been found to promote apoptosis [39, 40] . Its upregulation may also contribute to the apoptosis of the pancreas' beta cells. The CP gene is the potential target of two 5'UTR CoV-tar-miRNAs, i.e., hsa-miR-4747-5p and hsa-miR-5196-5p. The CP gene encodes a secretory plasma protein called Ceruloplasmin, the level of which is found to be increased in the Diabetic condition [41] . CP has the highest closeness and betweenness centrality, indicating its strong influence on the network. The CFB gene encodes for the secretory complement factor b, which links obesity to Diabetes. Its level is found to be increased during obesity. It is also linked to insulin resistance [42] . Although in our study, CFB gene expression does not seem to be affected by differential miRNAs targeting, its upregulation is linked with the upregulation of CP, SOCS3, and AGT genes as all these genes are co-expressed according to the gene-gene interaction network ( Figure 1 ). The AGT gene is the potential target of three 5'UTR CoV-tar-miRNAs, i.e., hsa-miR-1303, hsa-miR-298, and hsa-miR-4645-3p. AGT gene encodes for the secretory preangiotensinogen or angiotensinogen precursor, which is essential for the renin-angiotensin system to maintain the blood pressure and fluid and electrolyte homeostasis in the body [43, 44] . Our study suggests that after SARS-CoV-2 infection, these genes associated with Diabetes and cell death get upregulated due to differential miRNA targeting and lead the pancreas cells to death. This hypothesis can be applied to the insulin-producing pancreas beta cells, resulting in low or no insulin secretion. This can worsen the condition of COVID-19 patients due to diabetic complications or new onset of Diabetes, even if the pancreas cell's miRNAs block the viral genome function. The mechanism of differential miRNA targeting must be further validated in vitro. For therapeutic purposes, artificial miRNAs can be designed and inserted into the infected cells to bind with the viral genome, thus, blocking both its function as well as the differential targeting of the host cell's miRNAs. A research study with this concept has been done with four artificial miRNAs for targeting and blocking the Chikungunya virus genome [45] . This study suggests that in a SARS-CoV-2-infected human pancreas cell, the native miRNAs target the viral genome instead of the cell's mRNAs that were being targeted before the infection. This differential miRNA targeting causes the pancreas cell's genes associated with CoV-tar-miRNAs targeting the DEGs CoV-tar-miRNAs' raw microarray expression in human Pancreas WHO Director-General's opening remarks at the media briefing on COVID19 World Health Organization (2020) Naming the coronavirus disease (COVID-19) and the virus that causes it WHO Coronavirus Disease (COVID-19) Dashboard. URL: covid19.who.int Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study, The lancet COVID-19 and diabetes mellitus: from pathophysiology to clinical management Viruses with single-stranded, positive-sense RNA genomes, Molecular virology Elevated glucose levels favor SARS-CoV-2 infection and monocyte response through a HIF-1α/glycolysis-dependent axis New-onset diabetes in Covid-19 Diabetic ketoacidosis precipitated by Covid-19 in a patient with newly diagnosed diabetes mellitus Association of the insulin resistance marker TyG index with the severity and mortality of COVID-19 COVID-19 and Diabetes: a co-conspiracy? A pathological report of three COVID-19 cases by minimally invasive autopsies, Zhonghua bing li xue za zhi= Chinese journal of pathology Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein Binding of SARS coronavirus to its receptor damages islets and causes acute Diabetes On the importance of host microRNAs during viral infection What is the potential function of microRNAs as biomarkers and therapeutic targets in COVID-19? RNA viruses can hijack vertebrate microRNAs to suppress innate immunity Human microRNA hsa-miR-296-5p suppresses enterovirus 71 replication by targeting the viral genome 2020) A human pluripotent stem cell-based platform to study SARS-CoV-2 tropism and model virus infection in human cells and organoids Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 Heavy-tailed prior distributions for sequence count data: removing the noise and preserving large differences The genetic association database Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists GeneMANIA: Fast gene network construction and function prediction for Cytoscape Cytoscape: a software environment for integrated models of biomolecular interaction networks The GeneMANIA prediction server: 1. CP 3'UTR hsa-miR-466* hsa UTR hsa-miR-3123* hsa-miR-4691-3p* hsa UTR hsa-miR-1303 hsa-miR-298 hsa-miR-4645-3p hsa UTR hsa-miR-1236-3p hsa-miR-4279 hsa The authors declare that there are no conflicts of interest with the contents of this article. miRNA Name and its sequence Seed Location on the viral 3'UTR sequence The target sequence of miRNAs