Md. Abdullah-Al-Kamran Khan1
Md. Rabi Us Sany2
Md. Shafiqul Islam2
Abul Bashar Mir Md. Khademul Islam2*- 1Department of Mathematics and Natural Sciences, BRAC University, Dhaka, Bangladesh
- 2Department of Genetic Engineering & Biotechnology, University of Dhaka, Dhaka, Bangladesh
A detailed understanding of the molecular mechanism of SARS-CoV-2 pathogenesis is still elusive, and there is a need to address its deadly nature and to design effective therapeutics. Here, we present a study that elucidates the interplay between the SARS-CoV and SARS-CoV-2 viruses' and host's miRNAs, an epigenetic regulator, as a mode of pathogenesis; and we explored how the SARS-CoV and SARS-CoV-2 infections differ in terms of their miRNA-mediated interactions with the host and the implications this has in terms of disease complexity. We have utilized computational approaches to predict potential host and viral miRNAs and their possible roles in different important functional pathways. We have identified several putative host antiviral miRNAs that can target the SARS viruses and also predicted SARS viruses-encoded miRNAs targeting host genes. In silico predicted targets were also integrated with SARS-infected human cell microarray and RNA-seq gene expression data. A comparison between the host miRNA binding profiles on 67 different SARS-CoV-2 genomes from 24 different countries with respective country's normalized death count surprisingly uncovered some miRNA clusters, which are associated with increased death rates. We have found that induced cellular miRNAs can be both a boon and a bane to the host immunity, as they have possible roles in neutralizing the viral threat; conversely, they can also function as proviral factors. On the other hand, from over representation analysis, our study revealed that although both SARS-CoV and SARS-CoV-2 viral miRNAs could target broad immune-signaling pathways; only some of the SARS-CoV-2 miRNAs are found to uniquely target some immune-signaling pathways, such as autophagy, IFN-I signaling, etc., which might suggest their immune-escape mechanisms for prolonged latency inside some hosts without any symptoms of COVID-19. Furthermore, SARS-CoV-2 can modulate several important cellular pathways that might lead to the increased anomalies in patients with comorbidities like cardiovascular diseases, diabetes, breathing complications, etc. This might suggest that miRNAs can be a key epigenetic modulator behind the overcomplications amongst the COVID-19 patients. Our results support that miRNAs of host and SARS-CoV-2 can indeed play a role in the pathogenesis which can be further concluded with more experiments. These results will also be useful in designing RNA therapeutics to alleviate the complications from COVID-19.
Introduction
Coronavirus outbreaks have been reported over the past three decades, but the recent SARS-CoV-2 pandemic has reached more than 200 countries, has been the causative agent for the death of 58,392 people around the globe, and 1,087,374 coronavirus cases have been filed as of the date of writing this article (Worldometer, 2020). Among closed cases of SARS-CoV-2, 20% of the patients have died, and 5% of patients within active cases are in critical situations (Worldometer, 2020). The initial estimation of the SARS-CoV-2 death rate of 3.4%, as declared by the WHO (2020), requires refreshing, as the global casualty is on the rise. This novel virus requires novel and in-depth studies to promote new strategies for the management of this pandemic.
The coronavirus subfamily is a single-stranded positive-sense (+ssRNA) virus with a genome size of around 30 kb (Lu et al., 2020). The family is categorized into four subgenera as alpha, beta, gamma, and delta coronavirus (Cheng and Shan, 2020). SARS-CoV-2 is a beta coronavirus with a genome size of 29.9 kb (Accession no. NC_045512.2), 11 genes being reported in NCBI-Gene (2020). Phylogenetic analysis between SARS-CoV-2 and SARS-CoV showed ~79% similarity. Meanwhile, the distance is much greater for MERS-CoV (~50% similarity), but the closest relative to the SARS-CoV-2 is bat-derived SARS-like coronavirus (~90% similarity) (Jiang et al., 2020; Lu et al., 2020; Ren et al., 2020). Genomic analysis of SARS-CoV and SARS-CoV-2 has shown substitution of 380 amino acids and deletion of ORF8a, elongation of ORF8b (84 vs. 121 amino acid residues) and truncation of ORF3b (154aa in SARS-CoV whereas 22aa in SARS-CoV-2) (Lu et al., 2020).
MicroRNAs are small ncRNA molecules that regulate post-transcriptional-level gene expression; it has already been established that viruses use host machinery to produce miRNAs (Ambros, 2001). Although miRNA can be an important antiviral tool (Trobaugh and Klimstra, 2017) that can stimulate the innate and adaptive immune system (Ambros, 2001; Trobaugh and Klimstra, 2017) but that can also be a back door for viral propagation; it is non-antigenic, thereby modulating cellular pathways without triggering the host immune response (Cullen, 2013; Głobinska et al., 2014). For example, nucleocapsid protein of coronavirus OC43 binds miR-9 and activates NF-κB (Lai et al., 2014). Although host microRNAs are either utilized or regulated by viruses, viral miRNAs are another side of the coin; they regulate host gene expression, cellular proliferation, stress-related genes, and even viral gene expression (Cullen, 2010; Haasnoot and Berkhout, 2011; Lai et al., 2014). A summary explored how a number of DNA and RNA viruses produce miRNAs known as viral miRNAs (v-miRNAs) to evade the host immune response (Mishra et al., 2020). Novel viral miRNAs have been predicted to play an important role in neurological disorders as well (Islam et al., 2019). Among RNA viruses, HIV-1-encoded miR-H1, for example, can cause mononuclear cells apoptosis; H5N1 influenza virus-encoded miR-HA-3p targets host PCBP2 and contributes to the “cytokine storm” and mortality; and KUN-miR-1 of the West Nile virus targets the host's GATA4, which facilitates virus replication (Li and Zou, 2019). The interactions of host miRNAs with the SARS-CoV genome and viral proteins have been elucidated to suppress viral growth and immune evasion (Mallick et al., 2009). Novel classes of ncRNAs have been also observed by studies as possibly playing a definitive role in pathogenesis and survival (Liu et al., 2018). Respiratory viral infections caused by influenza, rhinovirus, adenovirus, RSV, and coronaviruses can be related to aberrant host miRNA expression, and their effect on the host can results in cell apoptosis, inhibition of immunologic pathways, downregulation of host antiviral responses, etc. (Mallick et al., 2009; Bondanese et al., 2014; Islam et al., 2019; Li and Zou, 2019; Mishra et al., 2020). Transmissible gastroenteritis virus (TGEV), though it induces significant IFN-I production after infection by inducing endoplasmic reticulum (ER), can evade antiviral effect of IFN-I by downregulating miR-30a-5p, which normally enhances IFN-I antiviral activity (Ma et al., 2018).
On the other hand, host miRNA expression plays a major role in controlling viral pathogenesis by mediating T cells and antiviral effector functions (Dickey et al., 2016). The first-reported example of a cellular miRNA that targets a viral RNA genome is miR-32, which targets the retrovirus PFV-1 transcripts and results in reduced virus replication (Lecellier et al., 2005). Similarly, miR-24, miR-93 can target VSV virus L and P proteins (Otsuka et al., 2007); miR-29a targets HIV Nef proteins (Ahluwalia et al., 2008) to inhibit replication; and miR-1, miR-30, miR-128, miR-196, miR-296, miR-351, miR-431, and miR-448 targets HCV C and NS5A proteins to inhibit translation/replication by inducing IFN signaling (Pedersen et al., 2007). miRNAs can therefore provide a different perspective in explaining the pathogenesis and infectivity of the novel SARS-CoV-2. Although SARS-CoV is distantly related to SARS-CoV-2, there are some similarities in their signs and symptoms, and they might be similar in pathogenesis, but there are crucial differences between the two diseases too (Xu et al., 2020). SARS-CoV-2 has infected many countries, and this has resulted in a stable mutation rate and some degree of variation (Xu et al., 2020). There is evidence that viral pathogens can have a novel immune evasion role by utilizing host miRNAs (Islam et al., 2019; Mishra et al., 2020).
The detailed miRNA-mediated epigenetic interplay between SARS-CoV-2 and its host is yet to be elucidated. It is not known what probable miRNAs produced by SARS-CoV-2 are affecting which human processes. Additionally, we do not know which antiviral miRNAs are taking part in host immunity. The genomic difference that, as a result, controls the host miRNA target sites and viral miRNAs might explain the difference between SARS-CoV and various isolated of SARS-CoV-2 in terms of pathogenesis and infectivity. Here in this study, we hypothesize on three potential effects of host and viral miRNA: (1) genomic differences between SARS-CoV and SARS-CoV-2 can lead to variations in host miRNA binding and differences in hence pathogenicity, signs and symptoms of these diseases and might explain the relatively longer incubation period of SARS-CoV-2; (2) on the other hand, there might be differences in viral miRNAs that can regulate expressions of different sets of host genes, which can in turn be advantageous to the virus or the host; and (3), due to a rapid mutation rate, observed variations among SARS-CoV-2 isolates in different regions of the world might result in variation in host capacities to target the virus with its miRNAs. This, in turn, might play a significant role in varying degrees of disease severity, symptoms, and mortality rate in different regions. In this study, we have carried out comparative analyses between SARS-CoV and SARS-CoV-2 with respect to host miRNA–viral genome interaction as well their differences based on region-specific isolates of SARS-CoV-2 and viral miRNA–host mRNA interactions to delineate the exclusive features of COVID-19 and their roles in viral survival and pathogenicity with respect to SARS-CoV (Figure 1).
Figure 1. Overall workflow of the whole study.
Materials and Methods
Obtaining SARS-CoV and SARS-CoV-2 Genome Sequences
The reference genomes of SARS-CoV (RefSeq Accession no. NC_004718.3) and SARS-CoV-2 (RefSeq Accession no. NC_045512.2) were retrieved from the NCBI RefSeq database (NCBI-RefSeq, 2020). A total of 67 whole-genome sequences of SARS-CoV-2 isolates covering 24 different countries (Supplementary File 1) were retrieved from NCBI-Virus (2020) and Shu and McCauley (2017).
Obtaining Human 3′UTR and Mature miRNA Sequences
Human miRNAs were accessed from microRNA database miRBase (Kozomara et al., 2018), and 3′UTR sequences of human protein-coding genes were obtained from Ensembl-Biomart (Hunt et al., 2018) (release 99).
Prediction of Viral Pre-miRNA and Validation
We used miRNAFold (Tav et al., 2016) for de novo prediction of all possible precursor-miRNAs from the obtained reference sequences of SARS-CoV and SARS-CoV-2 with a sliding window size of 150 and minimum hairpin size as 0. The results were validated using three different tools. First, RNAfold (Gruber et al., 2008) was used with minimum free energy (MFE) and partition function fold algorithm to find stable secondary structures. Second, a fixed-order Markov model-based algorithm namely FOMmiR (Shen et al., 2012) was used. Finally, a SVM-based tool iMiRNA-SSF (Chen et al., 2016) was used that calculates minimum free energy (MFE), p-value of randomization test (P-value), and the local triplet sequence-structure features. The common predictions from these three tools were utilized for further analysis.
Prediction of Mature miRNA
A Naive Bayes classifier algorithm implemented in tool MatureBayes (Gkirtzou et al., 2010) was used to identify mature miRNA candidates within the miRNA precursor sequences.
RNA-RNA Interaction Analysis
Three different tools were used to analyze RNA–RNA interactions for the host miRNA-viral genome and viral miRNA-host 3′UTR of coding sequences. IntaRNA 2.0 (Mann et al., 2017) was used considering sites with parameters –mode=H, –model=X, –outMode=C, and ΔΔG ≤ −10 kcal/mol, with seed 2–8, allowing for G:U base pairs. microRNA.org (Betel et al., 2008) was used with a score cutoff ≥ 140, energy cutoff ≤ −20 kcal/mol, gap opening = −9.0, and gap extension = −4.0; psRNATarget (Dai and Zhao, 2011) with the default parameter was also used to determine RNA-RNA interactions. Finally, the common predictions from these three tools were considered for downstream analysis.
Extraction of Targets of Host miRNAs
Targets of host miRNAs were mainly obtained using Funrich software (Pathan et al., 2017), which curates the experimentally validated targets of host miRNAs from different databases. The targets were also cross checked with the experimentally validated targets from miRTarBase database (Huang et al., 2019).
Target Genes Functional Enrichment Analysis
Enrichment Analysis in Gitools
The functional annotation of target genes is based on Gene Ontology (GO) (Ashburner et al., 2000), as extracted from the EnsEMBL (Hubbard et al., 2007) and KEGG pathway database (Kanehisa and Goto, 2000). Accordingly, all genes are classified into the ontology categories' biological process (GOBP) and pathways when possible. We have taken only the GO/pathway categories that have at least 10 genes annotated. We used Gitools for enrichment analysis and heatmap generation (Perez-Llamas and Lopez-Bigas, 2011). Resulting p-values were adjusted for multiple testing using the Benjamini and Hochberg's method of False Discovery Rate (FDR) (Benjamini and Hochberg, 1995).
Enrichment Analysis Using Web-Based Tools
The host miRNAs targeting SARS-CoV and SARS-CoV-2 were used for functional over-representation analysis to visualize and predict the roles of these miRNAs in human diseases and find enriched pathways. Besides Gitools, functional enrichment analyses of the target human genes were conducted using EnrichR (Kuleshov et al., 2016), DAVID 6.8 (Huang et al., 2009; Sherman and Lempicki, 2009), WebGestalt 2019 (Liao et al., 2019), and FunRich 3.1.3 (Pathan et al., 2017). The targeted genes are analyzed to determine their role in viral pathogenesis, infectivity, and immune evasion.
Microarray Expression Data Analysis
Microarray data for changes in gene expression induced by SARS-CoV on 2B4 cells, comparing the infected cells with SARS-CoV with the uninfected cells (for 12, 24, and 48 h), obtained from the Gene Expression Omnibus (GEO), ID GSE17400 (https://www.ncbi.nlm.nih.gov/geo) (Barrett et al., 2012). Raw Affymatrix CEL files were background corrected and normalized using the “rma” algorithm of Bioconductor package “affy” (version 1.28.1). The quality of microarray experiment (data not shown) was verified by Bioconductor package “arrayQualityMetrics” (Kauffmann et al., 2009) (version 3.2.4 under Bioconductor version 3.10; R version 3.6.0). To determine genes that are differentially expressed (DE) between two experimental conditions, Bioconductor package Limma (Smyth, 2005) was utilized to generate contrast matrices and fit the corresponding linear model. Probe annotations of genes were done using the Ensembl gene model (Ensembl version 99) as extracted from Biomart (Flicek et al., 2007) and by using in-house python script. When more than one probe was annotated to the same gene, the highest absolute expression value was considered (maximizing). To consider a gene is differentially expressed, multiple tests corrected FDR (Benjamini and Hochberg, 1995) p-value ≤ 0.05 was used as a cut-off.
RNA-seq Expression Data Analysis
RNA-seq raw read-count data on SARS-CoV-2 mediated expression changes in primary human lung epithelium (NHBE), and transformed lung alveolar (A549) cells were obtained from the GEO database (GSE147507) (Barrett et al., 2012). For the differential expression (DE) analysis, we used the Bioconductor package DESeq2 (version 1.38.0) (Anders and Huber, 2010) with R version 3.6.0 (R Core Team, 2016) with a model based on the negative binomial distribution. To avoid false positives, we considered only those transcripts where at least 10 reads were annotated and a p-value of 0.01.
MicroRNA Clustering
The hierarchal clustering of human miRNAs that could target SARS-CoV-2 genomes (binary mode) obtained from various countries was done using Manhattan distance and complete linkage analysis with the Genesis tool (Sturn et al., 2002). The human death count (per million population) due to SARS-CoV-2 infection was obtained on the 2nd of April, 2020, from the “Worldometer” website (Worldometer, 2020).
Overlap Analysis
Two- or three-way overlap analysis was done using online venn-diagram program Venny 2.1.0 (Oliveros, 2018). Multiple pairwise overlaps, correlation analyses, and heatmap generations were carried out using Gitools (Perez-Llamas and Lopez-Bigas, 2011).
Data Visualization
We have visualized human miRNAs that bind to the virus genome in web-genome browser NCBI genome data viewer (NCBI's-Genome-Browser, 2020).
Results
Several Human miRNAs Are Found to Target SARS-CoV and SARS-CoV-2
It is possible that, during viral infections, host-encoded miRNAs can modulate viral infections as a means of a host immune response (Girardi et al., 2018). To identify possible host miRNAs that can get induced during the SARS-CoV (R) and SARS-CoV-2 (R) infections, we have utilized a bioinformatics approach. From our rigorous analysis pipeline, which covers three different well-established algorithms (IntaRNA, miRanda, and psRNATarget) to predict RNA-RNA interactions, we have identified 122 and 106 host antiviral miRNAs against SARS-CoV (R) and SARS-CoV-2 (R), respectively (Figures 2A,B, Supplementary File 2). Amongst these, 27 miRNAs were found to be targeting both viruses (Figure 2A). Whilst comparing these miRNAs with the antiviral miRNAs from VIRmiRNA (Qureshi et al., 2014), we have found four (hsa-miR-654-5p, hsa-miR-198, hsa-miR-622, and hsa-miR-323a-5p) and three (hsa-miR-17-5p, hsa-miR-20b-5p, and hsa-miR-323a-5p) host miRNAs against SARS-CoV (R) and SARS-CoV-2 (R), respectively, to exhibit experimental evidence of having antiviral roles during infections (Figures 2A–C).
Figure 2. The host miRNAs targeting SARS-CoV and SARS-CoV-2. (A) Common host miRNAs and their target genes in SARS-CoV and SARS-CoV-2. (B) Host miRNAs and their target genes which uniquely target either SARS-CoV or SARS-CoV-2. (C) Venn diagram showing the common and unique host miRNAs targeting SARS-CoV and SARS-CoV-2 and host miRNAs that have exhibited experimental evidence as antiviral miRNAs.
Moreover, we compared the miRNAs targeting the two reference genomes of SARS-CoV (R) and SARS-CoV-2 (R), and we found most of the host miRNAs can target the ORF1ab region, followed by the S region as the second-most targeted (Figures 3A,B). Additionally, the M, N, ORF3a, ORF7a, ORF8 (ORF8a, ORF8b for SARS-CoV), 5′ UTR, and 3′ UTR regions of both viruses were targeted by host miRNAs. The significant variance was observed in the targeting positions of the host miRNAs between these two viruses (Figures 3A,B).
Figure 3. Genome browser view of host miRNAs targeting the regions of (A) SARS-CoV (Reference) and (B) SARS-CoV-2 (Reference) genomes.
Since the RNA virus mutates fast, it is conceivable that mutations in crucial genomic locations would lead to differences in host miRNA binding patterns. Therefore, the ability of the host miRNAs in targeting genomes of 67 SARS-CoV-2 isolates covering 24 different countries was also seen. Although, as expected, most of the identified host miRNAs' binding profiles across these isolates remained somewhat similar to that of SARS-CoV-2 reference sequence; interestingly, we have identified 24 host miRNAs that bind differentially across the isolates (Figure 4A), which might have occurred due to the genomic variations between these isolates. Complete linkage agglomerative hierarchal cluster (HCL) analysis with Manhattan distance of these miRNAs (binary mode, bind or not bind) revealed two major clusters with a side cluster for one South Korean and two Singaporean isolates (Figure 4B). As miRNA is crucial in both host defense and viral pathogenesis, to understand the significance of this cluster, we have also compared the host miRNA clusters with the death rate (normalized by per million population) from different countries. Surprisingly, relatively higher deaths are found to be more prominent in the European major clusters (right side cluster) compared to the other major cluster (left side), and we also found much lower deaths in side clusters (Figure 4B). However, many more factors are there which can also play pivotal roles in the susceptibility of the patients, this suggests only one interesting aspect for this correlation between the miRNA binding pattern and host susceptibility.
Figure 4. Differences of host miRNA-binding profiles, (A) representing only uncommon miRNAs binding pattern in 67 different SARS-CoV-2 genomes from 24 different countries, and (B) hierarchal clustering of all miRNAs binding in 67 genomes (upper panel, same country with same color code) and association of country-specific death rates (in color-coded scale) per million population (lower panel).
Host miRNAs Targeting SARS-CoV and SARS-CoV-2 Play Crucial Roles in Neutralizing the Virus
Though the primary action elicited by host miRNAs is to silence the viral RNA, they might also modulate some host factors which provide an edge to the viral pathogenesis. To find out if these particular pathways are also targeted by the host miRNAs induced by SARS-CoV and SARS-CoV-2 infections, we have performed miRNA pathway enrichment analysis. We have found that several such pathways might be deregulated by the host miRNAs to suppress the entry of the virus, prevent the spread of the virions, and to minimize the systemic symptoms resulting from the infection (Figure 5A).
Figure 5. Enrichment analysis and comparison between host miRNA targets induced by SARS-CoV and SARS-CoV-2 infections. (A) Heatmap representation of enriched pathways involved in host defense obtained using Funrich software. (B) Enriched pathways that might act as proviral mechanisms obtained using Funrich software. Enrichment of downregulated host miRNA target genes in SARS-CoV and SARS-CoV-2 using gitools. (C) GO Biological Process module. (D) GO Molecular Function module. (E) KEGG pathway modules. The significance of enrichment in terms of adjusted p-value (<0.05) is represented in a color-coded P-value scale for all heatmaps. Color toward red indicates higher significance and color toward yellow indicates less significance, while gray means non-significant. Only selected significant enriched terms are shown.
Host miRNAs might have a probable role in blocking the entry of the virus, as they are found to be targeting the pathways needed for viral entry- PDGF receptor-like signaling (Soroceanu et al., 2008), Arf-6 signaling (García-Expósito et al., 2011), PI3K-Akt signaling (Diehl and Schaal, 2013), EGFR signaling (Zheng et al., 2014), signaling events mediated by focal adhesion kinase (Elbahesh et al., 2014), CDC42 signaling (Swaine and Dittmar, 2015), the EphrinB-EPHB pathway (Wang et al., 2019), Cadherin signaling (Mateo et al., 2015), RTK signaling (Haqshenas and Doerig, 2019), etc. (Figure 5A).
They can also block some machinery like-p38 MAPK signaling (Hirasawa et al., 2003), FAK signaling (Elbahesh et al., 2014), PI3K-Akt signaling (Diehl and Schaal, 2013), etc., which can be hijacked by viruses for their efficient replication, pre-mRNA processing, and translation (Figure 5A). These host miRNAs might also try to reduce some host-induced inflammatory responses to prevent acute lung damage by targeting IGF1 signaling (Li et al., 2019), VEGF signaling (Alkharsah, 2018), PAR1 signaling (Heuberger and Schuepbach, 2019), integrin signaling (Teoh et al., 2015), TGF-beta signaling (Denney et al., 2018), TRAIL signaling (Cummins and Badley, 2009), etc. (Figure 5A). Some signaling pathways, such as CXCR4 signaling (Arnolds and Spencer, 2014), TGF-beta signaling (Denney et al., 2018), mTOR signaling (Le Sage et al., 2016), PI3K-Akt signaling (Diehl and Schaal, 2013), etc., can facilitate viral survival in infected cells by inhibiting apoptosis, autophagy, early immune responses, etc. Host miRNAs may function to downregulate these to invoke a proper immune response against the viruses (Figure 5A).
Infection Induced Host miRNAs Can Function as a Proviral Factor by Inhibiting Host Immune Surveillance Pathways
Host miRNAs can be like double-edged swords, as sometimes, they can facilitate viral immune evasion by targeting some important host immune responses (Bruscella et al., 2017). Our host miRNA enrichment analysis showed several significant pathways, such as IFN-gamma signaling (Kang et al., 2018), TGF-beta signaling (Mogensen and Paludan, 2001), Interleukin signaling (Kimura et al., 2013), IGF1 signaling (Li et al., 2019), TRAIL signaling (Cummins and Badley, 2009), etc. These are involved in important proinflammatory cytokine signaling during viral infections (Figure 5B). Interestingly, we have found out that host miRNAs induced during SARS-CoV-2 infection may particularly downregulate the signaling of different Toll-Like Receptors (TLRs) (Kimura et al., 2013), which are considered as the primary stimulatory molecules for producing host antiviral responses (i.e., production of interferons and other inflammatory cytokines) (Figure 5B). Also, other receptor signaling involved in antiviral responses like- uPA-UPAR signaling (Alfano et al., 2003), TRAF6 signaling (Konno et al., 2009), S1P1 signaling (Oldstone et al., 2013), Estrogen receptor signaling (Kovats, 2015), Protease-activated Receptor (PAR) signaling (Antoniak et al., 2013), Bone morphogenetic protein (BMP) signaling (Eddowes et al., 2019), etc. can also be deregulated by the host miRNAs, leading to the host's immune suppression (Figure 5B).
Host miRNAs' Targeted Downregulated Pathways Are Related to the Comorbidities of COVID-19
SARS-CoV-2-infected patients with comorbidities (i.e., cardiovascular diseases, diabetes, and renal problems) are found to be more susceptible to COVID-19. To find out whether host miRNAs play a role in these, we performed enrichment analyses using the downregulated targets genes of the host miRNAs using the expression data obtained from GEO dataset (GSE17400 for SARS-CoV and GSE147507 for SARS-CoV-2). These revealed that the downregulated targets of host miRNAs are involved in functions and pathways, such as heart development, kidney development, several neuronal processes, metabolic process, regulation of cellular ketone metabolism, insulin resistance, glucagon signaling pathway, fatty acid metabolism, and PPAR signaling (Figures 5C–E). Aberrant regulation of these processes can overcomplicate the disease conditions of patients having existing disorders.
Viral miRNAs Encoded by SARS-CoV and SARS-CoV-2 Can Target Several Host Genes
Many human viruses were found to produce miRNAs to assist in their overall pathogenesis by modulating host factors (Bruscella et al., 2017). Previous study on SARS-CoV also suggests that viral small non-coding RNAs can help its efficient pathogenesis (Morales et al., 2017). Our bioinformatics approach suggests that SARS-CoV and SARS-CoV-2 can also encode some viral miRNAs. The miRNAfold tool (Tav et al., 2016) yielded 529 and 519 putative pre-miRNAs from the genome of SARS-CoV and SARS-CoV-2, respectively. The RNAfold tool (Gruber et al., 2008) predicted 303 and 308 of these precursors of SARS-CoV and SARS-CoV-2, respectively are highly stable for forming hairpin structure which is a prerequisite of mature miRNA formation. Using FomMiR (Shen et al., 2012) and IMiRNA-SSF (Chen et al., 2016), we then predicted which of these highly stable precursors can truly produce mature miRNAs. We have found 63 and 85 such precursors, respectively for SARS-CoV and SARS-CoV-2. Using the Maturebayes tool from these precursors, we identified 126 and 170 mature miRNAs from SARS-CoV and SARS-CoV-2, respectively (Supplementary File 3). We predicted the human target genes by utilizing three different target prediction tools, and, reducing false positives, we have taken only the common set. This returned 5,292 and 6,369 human target genes for SARS-CoV and SARS-CoV-2, respectively (Supplementary File 4). Out of these, 2,992 genes are found to be common in both, while 2,300 and 3,377 genes were found to be unique targets of SARS-CoV and SARS-CoV, respectively. An apparent difference of the coding regions of miRNAs between SARS-CoV and SARS-CoV-2 was observed (Figures 6A,B).
Figure 6. Genome browser view of viral miRNAs transcribed from the regions of (A) SARS-CoV (Reference) and (B) SARS-CoV-2 (Reference) genomes.
SARS-CoV and SARS-CoV-2 Can Evade Host's Immune Surveillance Pathway by Utilizing Its miRNAs
Many viruses use their miRNAs to suppress or escape host's immune responses (Mishra et al., 2020). To identify which pathways are associated with SARS-CoV and SARS-CoV infection, we have performed the gene ontology (GO) and pathway functional enrichment of the targeted genes using different tools. This revealed a myriad of significant functions and pathways involved in host immune responses, such as Wnt signaling (Ljungberg et al., 2019), MAPK signaling (Kimura et al., 2013), T-cell-mediated immunity (Channappanavar et al., 2014), autophagy (Yordy and Iwasaki, 2011), FGF receptor binding (van Asten et al., 2018), TGF-beta signaling (Denney et al., 2018), VEGF signaling (Alkharsah, 2018), ErbB signaling (Zheng et al., 2014), mTOR signaling (Le Sage et al., 2016), and TNF-alpha signaling (Kimura et al., 2013) are particularly targeted by SARS-CoV-2 (Figures 7A–E).
Figure 7. Enrichment analysis and comparison between the SARS-CoV and SARS-CoV-2 encoded viral miRNAs' target human genes. Functional enrichment using gitools: the (A) GO Biological Process module and (B) GO Molecular Funtion module. Enriched pathways obtained from (C) the Webgestalt (KEGG and Wikipathways) tool, (D) DAVID (KEGG pathways) tool, and (E) EnrichR (KEGG, Wikipathways, BioPlanet pathways) tool. Color codes are as in Figure 5. Only selected significant enriched terms are shown.
Functions and pathways, such as heart development, brain development, and the insulin signaling pathway, etc. (Figures 7A–E), were also enriched for SARS-CoV-2 only, which can be targeted by the viral miRNAs, making the patients with previous complications more susceptible to COVID-19, and it can also lead to several signs uniquely found in SARS-CoV-2-infected patients.
We have also identified the downregulated target genes by curating the GEO expression datasets (GSE17400 for SARS-CoV and GSE147507 for SARS-CoV-2) and found 120 and 35 downregulated target genes in SARS-CoV and SARS-CoV-2, respectively (Supplementary File 5). These downregulated target genes are found to be involved in different immune-signaling pathways as well as different pathways related to organ-specific function (Figure 8).
Figure 8. Enrichment analysis and comparison between the enriched pathways of (A) SARS-CoV and (B) SARS-CoV-2 encoded viral miRNAs' downregulated target genes, obtained using the EnrichR (KEGG, Wikipathways, and BioPlanet pathways) tool. –Log10 (P-value) scale is utilized for all in a bar graph. The higher the bar height, the more significant an enriched term is. Only selected significant enriched terms are shown.
Discussion
Cellular miRNAs play a crucial role during the viral infection to strengthen host immunity by targeting virus genes as well as pathways that viruses utilize for their survival and immune evasion (Girardi et al., 2018). Viruses themselves can encode their miRNAs to target these immune-signaling pathways (Bruscella et al., 2017). COVID-19 has become a serious public health issue, though the complete molecular mechanism of pathogenesis is not fully understood yet. In this context, we have carried out this whole study to investigate the miRNA-mediated interactions between the host and the SARS-CoV-2 virus, which might shed some light on the tug-of-war between host's immune responses and virus's circumvention strategies. Though the disease conditions caused by SARS-CoV and SARS-CoV-2 are more or less similar, several unique features [i.e., long incubation, enhanced latency, asymptomatic infection, intense pain, severe lung damage, etc. (Ceccarelli et al., 2020)] of SARS-CoV-2 make it more challenging to manage compared to SARS-CoV. We also sought to find out if there are any existing differences between SARS-CoV and SARS-CoV-2 in the context of miRNA-mediated regulation of host responses.
As host miRNAs are one of the key sources of immune protection against viral infections, we have tried to find out which cellular miRNAs can target SARS-CoV and SARS-CoV-2 genes. Due to differences in the genome sequences between these two viruses, there was a significant difference between cellular miRNAs and their targeting viral genes. Likewise, some of the commonly found cellular miRNAs were showing differential binding preferences for these viral genes (Figure 2A). Previous study by Mallick et al. showed that cellular miRNAs can boost up host's immune response as well as they can assist in viral immune evasion mechanisms (Mallick et al., 2009). Another study by Morales et al. suggested that SARS-CoV can encode small non-coding RNAs that can play a role in inflammatory lung pathology (Morales et al., 2017). We identified that some of our predicted miRNAs have partial sequence similarities with the SARS-CoV svRNAs reported by Morales et al. (2017) (Supplementary File 6). We also compared the induced host miRNAs' profiles of 67 SARS-CoV-2 isolates from 24 different countries across the globe. From this analysis, we have identified several clusters and associated miRNAs, and our correlation study between these clusters with the death counts all over the world shed some light on the burning question and suggests why Europeans are more prone to COVID-19 (Figure 3B).
We found several miRNAs with experimentally validated antiviral roles; among those, hsa-miR-323a-5p and hsa-miR-654-5p (predicted for SARS-CoV) were found to inhibit viral replication in H1N1 Influenza virus infection (Song et al., 2010), while hsa-miR-17-5p and hsa-miR-20b-5p (predicted for SARS-CoV-2) were found to be upregulated in H7N9 Influenza virus infection (Zhu et al., 2014).
Apart from the basic role of cellular miRNAs in eliminating the transcripts of viruses, they can also modulate some host pathways which supposedly can be utilized by the infecting virus to avoid host's immune response. We also identified several such pathways involved in viral entry, replication, translation mechanisms, etc. These can be targeted by the cellular miRNAs induced by SARS-CoV and SARS-CoV-2 infection. Moreover, several immune-response pathways, such as TLR signaling, interleukin signaling, TRAF6 signaling, etc., were exclusively found to be targeted by SARS-CoV-2-induced host miRNAs (Figure 5B), and SARS-CoV-2-encoded miRNAs can target pathways, such as autophagy, IFN-I signaling, wnt signaling, mTOR signaling, etc. SARS-CoV-encoded miRNAs targets were, however, not found to be enriched in these pathways (Figures 7A–E). Target genes downregulated by SARS-CoV-2 miRNAs are found to be involved in the Ca2+ signaling pathway, and these are considered important activators of many signaling pathways (Zhou et al., 2009) (Figure 8B). All of these suggest why SARS-CoV-2 infections might be fatal for those who are immunosuppressed (D'Antiga, 2020).
Interestingly, our findings have enlightened several poorly understood mechanisms behind many of the unique clinical and pathological features of SARS-CoV-2, which has made it significantly different from SARS-CoV. We predicted both cellular miRNAs and viral encoded miRNAs, induced during SARS-CoV and SARS-CoV-2 infection, were found to target cytokine-signaling pathways involved in immune responses leading to the improved viral pathogenesis. Also, we found that SARS-CoV-2 miRNAs can target different important organ-specific cellular functions and pathways. We showed that SARS-CoV-2-encoded miRNAs can target the insulin-signaling pathway (Figure 7A, Supplementary Figure 1), and aberration of this pathway might overcomplicate the whole disease condition for COVID-19 patients with existing diabetic problems (Shimizu et al., 1980; del Campo et al., 2012). Our data also suggests that the SARS-CoV-2 miRNAs can target heart development-related pathways (Figure 7A, Supplementary Figure 1), which might lead to similar consequences like viral myocarditis (Dennert et al., 2008), making the disease more fatal for the patients with existing cardiovascular complications. These SARS-CoV-2-encoded miRNAs might also target genes associated with brain development (Figure 7A, Supplementary Figure 1), which might provide a clue about the neurological signs like headaches, vomiting, and nausea. SARS-CoV-2-induced host miRNAs can also downregulate kidney development and regulation of cellular ketone metabolic processes, etc. (Figure 5C), increasing the burden upon the kidneys (Kanikarla-Marie and Jain, 2016), which might be fatal for patients who have diabetes and kidney complications. HIF-1 signaling was also found to be targeted by SARS-CoV-2 miRNAs (Figure 7E, Supplementary Figure 1). This pathway is found to be associated with many viral infections, as HIF-1 plays an important role in cellular survival during hypoxic conditions (Santos and Andrade, 2017); COVID-19 patients suffer from the lack of oxygens due to breathing complications. This pathway might therefore play a crucial role in mitigating the condition, but viral miRNA-mediated deregulation of this pathway might result in severe consequences.
Our findings can explain that the interplay between host miRNAs and SARS-CoV-2 can promote viral pathogenesis by deregulating major antiviral immune-signaling pathways; furthermore, the resulting abnormal regulation of several host pathways might lead to an increased complications in the infected patients. Our study, which was conducted using machine learning and knowledgebase approaches, alongside further experiments, has the full potential to provide a more detailed understanding of the disease progression, and, based on these results, novel therapeutic interventions using RNA interference (RNAi) can be designed.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: Gene Expression Omnibus (GEO), https://www.ncbi.nlm.nih.gov/geo/, Accession: GSE147507 and GSE17400.
Author Contributions
AI conceived the project. AI and MK designed the workflow. MK, MS, and MI collected the data. All authors performed the analyses and wrote the manuscript. All authors contributed to the article and approved the submitted version.
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.
Acknowledgments
We acknowledge Rafeed Rahman Turjya for valuable suggestions. We would also like to thank all the authors who have kindly deposited and shared genome data on GISAID (https://www.gisaid.org/).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00765/full#supplementary-material
References
Ahluwalia, J. K., Khan, S. Z., Soni, K., Rawat, P., Gupta, A., Hariharan, M., et al. (2008). Human cellular microRNA hsa-miR-29a interferes with viral nef protein expression and HIV-1 replication. Retrovirology 5, 117–117. doi: 10.1186/1742-4690-5-117
Alfano, M., Sidenius, N., Blasi, F., and Poli, G. (2003). The role of urokinase-type plasminogen activator (uPA)/uPA receptor in HIV-1 infection. J. Leuk. Biol. 74, 750–756. doi: 10.1189/jlb.0403176
Alkharsah, K. R. (2018). VEGF upregulation in viral infections and its possible therapeutic implications. Int. J. Mol. Sci. 19:1642. doi: 10.3390/ijms19061642
Ambros, V. (2001). microRNAs: tiny regulators with great potential. Cell 107, 823–826. doi: 10.1016/S0092-8674(01)00616-X
Anders, S., and Huber, W. (2010). Differential expression analysis for sequence count data. Genome Biol. 11:R106. doi: 10.1186/gb-2010-11-10-r106
Antoniak, S., Owens, A. P. III., Baunacke, M., Williams, J. C., Lee, R. D., Weithäuser, A., et al. (2013). PAR-1 contributes to the innate immune response during viral infection. J. Clin. Invest. 123, 1310–1322. doi: 10.1172/JCI66125
Arnolds, K. L., and Spencer, J. V. (2014). CXCR4: a virus's best friend? Infect. Genet. Evol. 25, 146–156. doi: 10.1016/j.meegid.2014.04.018
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. Nat. Genet. 25, 25–29. doi: 10.1038/75556
Barrett, T., Wilhite, S. E., Ledoux, P., Evangelista, C., Kim, I. F., Tomashevsky, M., et al. (2012). NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 41, D991–D995. doi: 10.1093/nar/gks1193
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. B Methodol. 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
Betel, D., Wilson, M., Gabow, A., Marks, D. S., and Sander, C. (2008). The microRNA.org resource: targets and expression. Nucleic Acids Res. 36, D149–D153. doi: 10.1093/nar/gkm995
Bondanese, V. P., Francisco-Garcia, A., Bedke, N., Davies, D. E., and Sanchez-Elsner, T. (2014). Identification of host miRNAs that may limit human rhinovirus replication. World J. Biol. Chem. 5:437. doi: 10.4331/wjbc.v5.i4.437
Bruscella, P., Bottini, S., Baudesson, C., Pawlotsky, J.-M., Feray, C., and Trabucchi, M. (2017). Viruses and miRNAs: more friends than foes. Front. Microbiol. 8:824. doi: 10.3389/fmicb.2017.00824
Ceccarelli, M., Berretta, M., Rullo, E. V., Nunnari, G., and Cacopardo, B. (2020). Editorial–differences and similarities between severe acute respiratory syndrome (SARS)-CoronaVirus (CoV) and SARS-CoV-2. Would a rose by another name smell as sweet? Eur. Rev. Med. Pharmacol. Sci. 24, 2781–2783. doi: 10.26355/eurrev_202003_20551
Channappanavar, R., Zhao, J., and Perlman, S. (2014). T cell-mediated immune response to respiratory coronaviruses. Immunol. Res. 59, 118–128. doi: 10.1007/s12026-014-8534-z
Chen, J., Wang, X., and Liu, B. (2016). IMiRNA-SSF: improving the identification of MicroRNA precursors by combining negative sets with different distributions. Sci. Rep. 6, 19062–19062. doi: 10.1038/srep19062
Cheng, Z. J., and Shan, J. (2020). 2019 Novel coronavirus: where we are and what we know. Infection 48, 155–163. doi: 10.1007/s15010-020-01401-y
Cullen, B. R. (2010). Five questions about viruses and microRNAs. PLoS Pathog. 6:e1000787. doi: 10.1371/journal.ppat.1000787
Cullen, B. R. (2013). MicroRNAs as mediators of viral evasion of the immune system. Nat. Immunol. 14:205. doi: 10.1038/ni.2537
Cummins, N., and Badley, A. (2009). The TRAIL to viral pathogenesis: the good, the bad and the ugly. Curr. Mol. Med. 9, 495–505. doi: 10.2174/156652409788167078
Dai, X., and Zhao, P. X. (2011). psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res. 39, W155–W159. doi: 10.1093/nar/gkr319
D'Antiga, L. (2020). Coronaviruses and immunosuppressed patients. The facts during the third epidemic. Liver Transpl. 26, 832–834. doi: 10.1002/lt.25756
del Campo, J. A., García-Valdecasas, M., Rojas, L., Rojas, Á., and Romero-Gómez, M. (2012). The hepatitis C virus modulates insulin signaling pathway in vitro promoting insulin resistance. PLoS ONE 7:e47904. doi: 10.1371/journal.pone.0047904
Dennert, R., Crijns, H. J., and Heymans, S. (2008). Acute viral myocarditis. Eur. Heart J. 29, 2073–2082. doi: 10.1093/eurheartj/ehn296
Denney, L., Branchett, W., Gregory, L. G., Oliver, R. A., and Lloyd, C. M. (2018). Epithelial-derived TGF-β1 acts as a pro-viral factor in the lung during influenza A infection. Mucos. Immunol. 11, 523–535. doi: 10.1038/mi.2017.77
Dickey, L. L., Worne, C. L., Glover, J. L., Lane, T. E., and O'Connell, R. M. (2016). MicroRNA-155 enhances T cell trafficking and antiviral effector function in a model of coronavirus-induced neurologic disease. J. Neuroinflamm. 13, 240–240. doi: 10.1186/s12974-016-0699-z
Diehl, N., and Schaal, H. (2013). Make yourself at home: viral hijacking of the PI3K/Akt signaling pathway. Viruses 5, 3192–3212. doi: 10.3390/v5123192
Eddowes, L. A., Al-Hourani, K., Ramamurthy, N., Frankish, J., Baddock, H. T., Sandor, C., et al. (2019). Antiviral activity of bone morphogenetic proteins and activins. Nat. Microbiol. 4, 339–351. doi: 10.1038/s41564-018-0301-9
Elbahesh, H., Cline, T., Baranovich, T., Govorkova, E. A., Schultz-Cherry, S., and Russell, C. J. (2014). Novel roles of focal adhesion kinase in cytoplasmic entry and replication of influenza A viruses. J. Virol. 88, 6714–6728. doi: 10.1128/JVI.00530-14
Flicek, P., Aken, B. L., Beal, K., Ballester, B., Cáccamo, M., Chen, Y., et al. (2007). Ensembl 2008. Nucleic Acids Res. 36, D707–D714. doi: 10.1093/nar/gkm988
García-Expósito, L., Barroso-González, J., Puigdomènech, I., Machado, J.-D., Blanco, J., and Valenzuela-Fernández, A. (2011). HIV-1 requires Arf6-mediated membrane dynamics to efficiently enter and infect T lymphocytes. Mol. Biol. Cell 22, 1148–1166. doi: 10.1091/mbc.e10-08-0722
Girardi, E., López, P., and Pfeffer, S. (2018). On the importance of host microRNAs during viral infection. Front. Genet. 9:439. doi: 10.3389/fgene.2018.00439
Gkirtzou, K., Tsamardinos, I., Tsakalides, P., and Poirazi, P. (2010). MatureBayes: a probabilistic algorithm for identifying the mature miRNA within novel precursors. PLoS ONE 5:e11843. doi: 10.1371/journal.pone.0011843
Głobinska, A., Pawełczyk, M., and Kowalski, M. L. (2014). MicroRNAs and the immune response to respiratory virus infections. Expert Rev. Clin. Immunol. 10, 963–971. doi: 10.1586/1744666X.2014.913482
Gruber, A. R., Lorenz, R., Bernhart, S. H., Neuböck, R., and Hofacker, I. L. (2008). The Vienna RNA websuite. Nucleic Acids Res. 36, W70–W74. doi: 10.1093/nar/gkn188
Haasnoot, J., and Berkhout, B. (2011). RNAi and cellular miRNAs in infections by mammalian viruses. Methods Mol. Biol. 721, 23–41. doi: 10.1007/978-1-61779-037-9_2
Haqshenas, G., and Doerig, C. (2019). Targeting of host cell receptor tyrosine kinases by intracellular pathogens. Sci. Signal. 12:eaau9894. doi: 10.1126/scisignal.aau9894
Heuberger, D. M., and Schuepbach, R. A. (2019). Protease-activated receptors (PARs): mechanisms of action and potential therapeutic modulators in PAR-driven inflammatory diseases. Thromb. J. 17:4. doi: 10.1186/s12959-019-0194-8
Hirasawa, K., Kim, A., Han, H.-S., Han, J., Jun, H.-S., and Yoon, J.-W. (2003). Effect of p38 mitogen-activated protein kinase on the replication of encephalomyocarditis virus. J. Virol. 77, 5649–5656. doi: 10.1128/JVI.77.10.5649-5656.2003
Huang, D. W., Sherman, B. T., and Lempicki, R. A. (2009). Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 37, 1–13. doi: 10.1093/nar/gkn923
Huang, H.-Y., Lin, Y.-C.-D., Li, J., Huang, K.-Y., Shrestha, S., Hong, H.-C., et al. (2019). miRTarBase 2020: updates to the experimentally validated microRNA–target interaction database. Nucleic Acids Res. 48, D148–D154. doi: 10.1093/nar/gkz896
Hubbard, T. J. P., Aken, B. L., Beal, K., Ballester, B., Cáccamo, M., Chen, Y., et al. (2007). Ensembl 2007. Nucleic Acids Res. 35, D610–D617. doi: 10.1093/nar/gkl996
Hunt, S. E., McLaren, W., Gil, L., Thormann, A., Schuilenburg, H., Sheppard, D., et al. (2018). Ensembl variation resources. Database 2018:bay119. doi: 10.1093/database/bay119
Islam, M. S., Khan, M. A.-A.-K., Murad, M. W., Karim, M., and Islam, A. B. M. M.K. (2019). In silico analysis revealed Zika virus miRNAs associated with viral pathogenesis through alteration of host genes involved in immune response and neurological functions. J. Med. Virol. 91, 1584–1594. doi: 10.1002/jmv.25505
Jiang, S., Du, L., and Shi, Z. (2020). An emerging coronavirus causing pneumonia outbreak in Wuhan, China: calling for developing therapeutic and prophylactic strategies. Emerg. Microbes Infect. 9, 275–277. doi: 10.1080/22221751.2020.1723441
Kanehisa, M., and Goto, S. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30. doi: 10.1093/nar/28.1.27
Kang, S., Brown, H. M., and Hwang, S. (2018). Direct antiviral mechanisms of interferon-gamma. Immune Netw. 18:e33. doi: 10.4110/in.2018.18.e33
Kanikarla-Marie, P., and Jain, S. K. (2016). Hyperketonemia and ketosis increase the risk of complications in type 1 diabetes. Free Radic. Biol. Med. 95, 268–277. doi: 10.1016/j.freeradbiomed.2016.03.020
Kauffmann, A., Gentleman, R., and Huber, W. (2009). arrayQualityMetrics—a bioconductor package for quality assessment of microarray data. Bioinformatics 25, 415–416. doi: 10.1093/bioinformatics/btn647
Kimura, H., Yoshizumi, M., Ishii, H., Oishi, K., and Ryo, A. (2013). Cytokine production and signaling pathways in respiratory virus infection. Front. Microbiol. 4:276. doi: 10.3389/fmicb.2013.00276
Konno, H., Yamamoto, T., Yamazaki, K., Gohda, J., Akiyama, T., Semba, K., et al. (2009). TRAF6 establishes innate immune responses by activating NF-κB and IRF7 upon sensing cytosolic viral RNA and DNA. PLoS ONE 4:e5674. doi: 10.1371/journal.pone.0005674
Kovats, S. (2015). Estrogen receptors regulate innate immune cells and signaling pathways. Cell. Immunol. 294, 63–69. doi: 10.1016/j.cellimm.2015.01.018
Kozomara, A., Birgaoanu, M., and Griffiths-Jones, S. (2018). miRBase: from microRNA sequences to function. Nucleic Acids Res. 47, D155–D162. doi: 10.1093/nar/gky1141
Kuleshov, M. V., Jones, M. R., Rouillard, A. D., Fernandez, N. F., Duan, Q., Wang, Z., et al. (2016). Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 44, W90–W97. doi: 10.1093/nar/gkw377
Lai, F. W., Stephenson, K. B., Mahony, J., and Lichty, B. D. (2014). Human coronavirus OC43 nucleocapsid protein binds microRNA 9 and potentiates NF-κB activation. J. Virol. 88, 54–65. doi: 10.1128/JVI.02678-13
Le Sage, V., Cinti, A., Amorim, R., and Mouland, A. J. (2016). Adapting the stress response: viral subversion of the mTOR signaling pathway. Viruses 8:152. doi: 10.3390/v8060152
Lecellier, C.-H., Dunoyer, P., Arar, K., Lehmann-Che, J., Eyquem, S., Himber, C., et al. (2005). A cellular microRNA mediates antiviral defense in human cells. Science 308, 557–560. doi: 10.1126/science.1108784
Li, G., Zhou, L., Zhang, C., Shi, Y., Dong, D., Bai, M., et al. (2019). Insulin-like growth factor 1 regulates acute inflammatory lung injury mediated by influenza virus infection. Front. Microbiol. 10:2541. doi: 10.3389/fmicb.2019.02541
Li, X., and Zou, X. (2019). An overview of RNA virus-encoded microRNAs. ExRNA 1:37. doi: 10.1186/s41544-019-0037-6
Liao, Y., Wang, J., Jaehnig, E. J., Shi, Z., and Zhang, B. (2019). WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res. 47, W199–W205. doi: 10.1093/nar/gkz401
Liu, C., Chen, Z., Hu, Y., Ji, H., Yu, D., Shen, W., et al. (2018). Complemented palindromic small RNAs first discovered from SARS coronavirus. Genes 9, 442–442. doi: 10.3390/genes9090442
Ljungberg, J. K., Kling, J. C., Tran, T. T., and Blumenthal, A. (2019). Functions of the WNT signaling network in shaping host responses to infection. Front. Immunol. 10:2521. doi: 10.3389/fimmu.2019.02521
Lu, R., Zhao, X., Li, J., Niu, P., Yang, B., Wu, H., et al. (2020). Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding. Lancet 395, 565–574. doi: 10.1016/S0140-6736(20)30251-8
Ma, Y., Wang, C., Xue, M., Fu, F., Zhang, X., Li, L., et al. (2018). The coronavirus transmissible gastroenteritis virus evades the type I interferon response through IRE1α-mediated manipulation of the microRNA miR-30a-5p/SOCS1/3 axis. J. Virol. 92:e00728-18. doi: 10.1128/JVI.00728-18
Mallick, B., Ghosh, Z., and Chakrabarti, J. (2009). MicroRNome analysis unravels the molecular basis of SARS infection in bronchoalveolar stem cells. PLoS ONE 4:e7837. doi: 10.1371/journal.pone.0007837
Mann, M., Wright, P. R., and Backofen, R. (2017). IntaRNA 2.0: enhanced and customizable prediction of RNA-RNA interactions. Nucleic Acids Res. 45, W435–W439. doi: 10.1093/nar/gkx279
Mateo, M., Generous, A., Sinn, P. L., and Cattaneo, R. (2015). Connections matter—how viruses use cell–cell adhesion components. J. Cell Sci. 128, 431–439. doi: 10.1242/jcs.159400
Mishra, R., Kumar, A., Ingle, H., and Kumar, H. (2020). The interplay between viral-derived miRNAs and host immunity during infection. Front. Immunol. 10:3079. doi: 10.3389/fimmu.2019.03079
Mogensen, T. H., and Paludan, S. R. (2001). Molecular pathways in virus-induced cytokine production. Microbiol. Mol. Biol. Rev. 65, 131–150. doi: 10.1128/MMBR.65.1.131-150.2001
Morales, L., Oliveros, J. C., Fernandez-Delgado, R., tenOever, B. R., Enjuanes, L., and Sola, I. (2017). SARS-CoV-encoded small RNAs contribute to infection-associated lung pathology. Cell Host Microbe 21, 344–355. doi: 10.1016/j.chom.2017.01.015
NCBI-Gene (2020). RE: Gene Links for Nucleotide (Select 1798174254)—Gene—NCBI. Bethesda, MD: National Center for Biotechnology Information (NCBI).
NCBI-RefSeq (2020). RE: RefSeq: NCBI Reference Sequence Database. Bethesda, MD: National Center for Biotechnology Information (NCBI).
NCBI's-Genome-Browser (2020). RE: NCBI's Genome Browser for Human (Homo sapiens)-Genome Data Viewer. Bethesda, MD: National Center for Biotechnology Information (NCBI).
NCBI-Virus (2020). RE: NCBI Virus. Bethesda, MD: National Center for Biotechnology Information (NCBI).
Oldstone, M. B. A., Teijaro, J. R., Walsh, K. B., and Rosen, H. (2013). Dissecting influenza virus pathogenesis uncovers a novel chemical approach to combat the infection. Virology 435, 92–101. doi: 10.1016/j.virol.2012.09.039
Otsuka, M., Jing, Q., Georgel, P., New, L., Chen, J., Mols, J., et al. (2007). Hypersusceptibility to vesicular stomatitis virus infection in Dicer1-deficient mice is due to impaired miR24 and miR93 expression. Immunity 27, 123–134. doi: 10.1016/j.immuni.2007.05.014
Pathan, M., Keerthikumar, S., Chisanga, D., Alessandro, R., Ang, C.-S., Askenase, P., et al. (2017). A novel community driven software for functional enrichment analysis of extracellular vesicles data. J. Extracell. Vesicles 6, 1321455–1321455. doi: 10.1080/20013078.2017.1321455
Pedersen, I. M., Cheng, G., Wieland, S., Volinia, S., Croce, C. M., Chisari, F. V., et al. (2007). Interferon modulation of cellular microRNAs as an antiviral mechanism. Nature 449, 919–922. doi: 10.1038/nature06205
Perez-Llamas, C., and Lopez-Bigas, N. (2011). Gitools: analysis and visualisation of genomic data using interactive heat-maps. PLoS ONE 6:e19541. doi: 10.1371/journal.pone.0019541
Qureshi, A., Thakur, N., Monga, I., Thakur, A., and Kumar, M. (2014). VIRmiRNA: a comprehensive resource for experimentally validated viral miRNAs and their targets. Database (Oxford) 2014:bau103. doi: 10.1093/database/bau103
R Core Team (2016). R: A Language and Environment for Statistical Computing. Vienna. Available online at: https://www.R-project.org/
Ren, L.-L., Wang, Y.-M., Wu, Z.-Q., Xiang, Z.-C., Guo, L., Xu, T., et al. (2020). Identification of a novel coronavirus causing severe pneumonia in human: a descriptive study. Chin. Med. J. 133, 1015–1024. doi: 10.1097/CM9.0000000000000722
Santos, S. A. D., and Andrade, D. R. J. (2017). HIF-1alpha and infectious diseases: a new frontier for the development of new therapies. Rev. Inst. Med. Trop. Sao Paulo 59:e92. doi: 10.1590/s1678-9946201759092
Shen, W., Chen, M., Wei, G., and Li, Y. (2012). MicroRNA prediction using a fixed-order Markov model based on the secondary structure pattern. PLoS ONE 7:e48236. doi: 10.1371/journal.pone.0048236
Sherman, B. T., and Lempicki, R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4:44. doi: 10.1038/nprot.2008.211
Shimizu, F., Hooks, J. J., Kahn, C. R., and Notkins, A. L. (1980). Virus-induced decrease of insulin receptors in cultured human cells. J. Clin. Invest. 66, 1144–1151. doi: 10.1172/JCI109944
Shu, Y., and McCauley, J. (2017). GISAID: Global initiative on sharing all influenza data - from vision to reality. Euro Surveill. 22:30494. doi: 10.2807/1560-7917.ES.2017.22.13.30494
Smyth, G. K. (2005). “limma: linear models for microarray data,” in Bioinformatics and Computational Biology Solutions Using R and Bioconductor. Statistics for Biology and Health, eds R. Gentleman, V. J. Carey, W. Huber, R. A. Irizarry, and S. Dudoit (New York, NY: Springer), 397–420.
Song, L., Liu, H., Gao, S., Jiang, W., and Huang, W. (2010). Cellular microRNAs inhibit replication of the H1N1 influenza A virus in infected cells. J. Virol. 84, 8849–8860. doi: 10.1128/JVI.00456-10
Soroceanu, L., Akhavan, A., and Cobbs, C. S. (2008). Platelet-derived growth factor-alpha receptor activation is required for human cytomegalovirus infection. Nature 455, 391–395. doi: 10.1038/nature07209
Sturn, A., Quackenbush, J., and Trajanoski, Z. (2002). Genesis: cluster analysis of microarray data. Bioinformatics 18, 207–208. doi: 10.1093/bioinformatics/18.1.207
Swaine, T., and Dittmar, M. T. (2015). CDC42 use in viral cell entry processes by RNA viruses. Viruses 7, 6526–6536. doi: 10.3390/v7122955
Tav, C., Tempel, S., Poligny, L., and Tahi, F. (2016). miRNAFold: a web server for fast miRNA precursor prediction in genomes. Nucleic Acids Res. 44, W181–W184. doi: 10.1093/nar/gkw459
Teoh, C. M., Tan, S. S. L., and Tran, T. (2015). Integrins as therapeutic targets for respiratory diseases. Curr. Mol. Med. 15, 714–734. doi: 10.2174/1566524015666150921105339
Trobaugh, D. W., and Klimstra, W. B. (2017). MicroRNA regulation of RNA virus replication and pathogenesis. Trends Mol. Med. 23, 80–93. doi: 10.1016/j.molmed.2016.11.003
van Asten, S. D., Raaben, M., Nota, B., and Spaapen, R. M. (2018). Secretome screening reveals fibroblast growth factors as novel inhibitors of viral replication. J. Virol. 92:e00260-18. doi: 10.1128/JVI.00260-18
Wang, J., Zheng, X., Peng, Q., Zhang, X., and Qin, Z. (2019). Eph receptors: the bridge linking host and virus. Cell. Mol. Life Sci. 77, 2355–2365. doi: 10.1007/s00018-019-03409-6
WHO (2020). WHODirector-General's Opening Remarks at the Media Briefing on COVID-19. Geneva: WHO. Available online at: https://www.who.int/dg/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19—6-april-2020 WHO (2020). RE: WHODirector-General's Opening Remarks at the Media Briefing on COVID-19.
Worldometer (2020). Coronavirus Cases. Worldometer. Available online at: https://www.worldometers.info/coronavirus/ (accessed April 4, 2020).
Xu, J., Zhao, S., Teng, T., Abdalla, A. E., Zhu, W., Xie, L., et al. (2020). Systematic comparison of two animal-to-human transmitted human coronaviruses: SARS-CoV-2 and SARS-CoV. Viruses 12, 244–244. doi: 10.3390/v12020244
Yordy, B., and Iwasaki, A. (2011). Autophagy in the control and pathogenesis of viral infection. Curr. Opin. Virol. 1, 196–203. doi: 10.1016/j.coviro.2011.05.016
Zheng, K., Kitazato, K., and Wang, Y. (2014). Viruses exploit the function of epidermal growth factor receptor. Rev. Med. Virol. 24, 274–286. doi: 10.1002/rmv.1796
Zhou, Y., Frey, T. K., and Yang, J. J. (2009). Viral calciomics: interplays between Ca2+ and virus. Cell. Calcium 46, 1–17. doi: 10.1016/j.ceca.2009.05.005
Keywords: SARS-CoV-2, COVID-19, miRNA-microRNA, viral pathogenesis, immune regulation
Citation: Khan MA-A-K, Sany MRU, Islam MS and Islam ABMMK (2020) Epigenetic Regulator miRNA Pattern Differences Among SARS-CoV, SARS-CoV-2, and SARS-CoV-2 World-Wide Isolates Delineated the Mystery Behind the Epic Pathogenicity and Distinct Clinical Characteristics of Pandemic COVID-19. Front. Genet. 11:765. doi: 10.3389/fgene.2020.00765
Received: 06 April 2020; Accepted: 29 June 2020;
Published: 10 July 2020.
Edited by:
William Cho, Queen Elizabeth Hospital (QEH), Hong KongReviewed by:
Suman Ghosal, National Institutes of Health (NIH), United StatesShaoli Das, National Institutes of Health (NIH), United States
Dong Xu, University of Missouri, United States
Copyright © 2020 Khan, Sany, Islam and Islam. 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: Abul Bashar Mir Md. Khademul Islam, khademul@du.ac.bd