- 1Translational Bioinformatics Group, International Centre for Genetic Engineering and Biotechnology (ICGEB), New Delhi, India
- 2Department of Biochemistry, Jamia Hamdard, New Delhi, India
Understanding the host regulatory mechanisms opposing virus infection and virulence can provide actionable insights to identify novel therapeutics against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We have used a network biology approach to elucidate the crucial factors involved in host responses involving host–microRNA (miRNA) interactions with host and virus genes using recently published experimentally verified protein–protein interaction data. We were able to identify 311 host genes to be potentially targetable by 2,197 human miRNAs. These miRNAs are known to be involved in various biological processes, such as T-cell differentiation and activation, virus replication, and immune system. Among these, the anti-viral activity of 38 miRNAs to target 148 host genes is experimentally validated. Six anti-viral miRNAs, namely, hsa-miR-1-3p, hsa-miR-17-5p, hsa-miR-199a-3p, hsa-miR-429, hsa-miR-15a-5p, and hsa-miR-20a-5p, are previously reported to be anti-viral in respiratory diseases and were found to be downregulated. The interaction network of the 2,197 human miRNAs and interacting transcription factors (TFs) enabled the identification of 51 miRNAs to interact with 77 TFs inducing activation or repression and affecting gene expression of linked genes. Further, from the gene regulatory network analysis, the top five hub genes HMOX1, DNMT1, PLAT, GDF1, and ITGB1 are found to be involved in interferon (IFN)-α2b induction, epigenetic modification, and modulation of anti-viral activity. The comparative miRNAs target identification analysis in other respiratory viruses revealed the presence of 98 unique host miRNAs targeting SARS-CoV-2 genome. Our findings identify prioritized key regulatory interactions that include miRNAs and TFs that provide opportunities for the identification of novel drug targets and development of anti-viral drugs.
Introduction
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a single-stranded positive-sense RNA β-coronavirus of the Coronaviridae family and shares the highest similarity with SARS-CoV, which is the virus responsible for the 2003 SARS outbreak (Wang et al., 2020).
SARS-CoV-2 genome is ∼29–30 kb and is translated into 29 proteins, including structural and non-structural proteins. In SARS-CoV-2 genome, more than two-thirds of the genome constitutes orf1ab, encoding orf1ab polyproteins at 5’, whereas the other one-third is associated with genes encoding structural proteins that include surface (S), envelope (E), membrane (M), and nucleocapsid (N) proteins at the 3′ end genome. Moreover, the SARS-CoV-2 genome codes for six accessory proteins, encoded by ORF3a, ORF6, ORF7a, ORF7b, and ORF8 genes (Khailany et al., 2020).
Viruses are obligate intracellular parasites that use the host cellular machinery to replicate and propagate. RNA viruses tend to evolve rapidly by mutations, enabling its evasion from the host immune response. With evolution, the hosts also develop different ways to fight virus infection, which include innate immunity that provides the first line of defense against viral infections. The host cells have many receptors that recognized virus elements that led to the activation of the interferon system and cytokines (Nakhaei et al., 2009; Girardi et al., 2018). Coronaviruses are known to induce the activation of host pathways linked to stress, apoptosis, autophagy, and innate immunity (Fung and Liu, 2019).
MicroRNAs (miRNAs) are ∼19- to 24-nt non-coding RNAs that regulate gene expression by binding to target messenger RNAs (mRNAs). Several miRNAs play a key role in cell differentiation, development, pathogenesis, cellular growth, apoptosis, and disease progression (Lou et al., 2019). During viral infections, host miRNAs are involved in various signaling pathways modulating host–virus interactions. Host miRNAs regulate viral infectivity and transmission, activate anti-viral immune responses, and play a role in many viral diseases, such as HIV, herpesvirus, and Ebola, by downregulating the host genes (Bernier and Sagan, 2018).
Targeting the host factors has been proven to be useful for restricting viral infections and has a potential for the development of host-directed therapies (HDTs) (Kaufmann et al., 2018). It is predicted that host miRNAs target over two-thirds of all human genes (Friedman et al., 2009). The latest release of miRBase consists of 2,654 mature human miRNAs with the ability to target hundreds of different types of genes and proteins (Griffiths-Jones, 2006).
Transcription factors (TFs) are master regulators that regulate gene transcription by binding to the promoter region of the DNA-binding domains of target genes. Moreover, it is known that there is an interplay between host miRNAs and TFs (Tong et al., 2019) during viral infections as shown in Figure 1, and both are regulated in a coordinated fashion in various biological processes, such as cell proliferation, differentiation, and apoptosis, and many diseases through feed-forward loops (FFLs) and feedback loops (FBLs) (Tsang et al., 2007; Mohamed et al., 2019). Therefore, the identification of key regulatory miRNA–TF interactions can play a crucial role in the development of anti-viral therapies against SARS-CoV-2.
Figure 1. Four most anticipated patterns of host miRNA-mediated virus–host interactions (FFLs). (A) A viral protein and human miRNA have common human target. (B) A viral protein and human miRNA have a protein pair as common target. (C) A viral protein targeting a TF of a human miRNA. (D) A viral protein targeting a host protein forming a feedback loop to a human miRNA.
Recently, using computational tools and databases, few groups of scientists have identified potential host miRNAs that target SARS-CoV-2 genes involved in immune signaling pathways (Ahmed et al., 2020, p. 19; Khan et al., 2020). However, these are limited reports focused on the identification of host miRNA and lack information about their regulatory interactions especially with TFs. Therefore, in the present study, we have investigated miRNAs targeting host genes and its interactions with TFs exploiting the recently published SARS-CoV-2–protein interaction map (Gordon et al., 2020). Out of the 332 host genes interacting with the virus genes, 311 are found to be targeted by 2,197 human miRNAs and involved in various biological processes, such as T-cell differentiation, onco-miRNAs, toxicity, immune system, and viral replications. After an extensive literature survey, we were able to identify 38 miRNAs with anti-viral activity targeting 143 host genes. Six anti-viral miRNAs, namely, hsa-miR-1-3p, hsa-miR-17-5p, hsa-miR-199a-3p, hsa-miR-429, hsa-miR-15a-5p, and hsa-miR-20a-5p, were reported to play a role in respiratory diseases, such as influenza A, adenovirus 2, and respiratory syncytial virus (RSV), and were downregulated that can be used to design anti-viral drugs. From the comparative gene expression analysis, we identified 148 differentially regulated TFs in SARS, whereas in SARS-CoV-2 gene expression analysis, we identified that 48 and 52 TFs are expressed in A549 and normal human bronchial epithelial (NHBE) cell lines. Interestingly, among the differentially expressed TFs, only STAT1 and STAT2 are in the interactome connected with ORF3a and NSP7.
Methodology
To identify miRNA–TF co-regulatory interactions in SARS-CoV-2, we have designed a pipeline as shown in Figure 2. Summarily, a list of host proteins interacting with virus proteins served as a starting dataset for the generation of different regulatory networks for the study.
Figure 2. Pipeline design to identify novel drug targets in SARS-CoV-2 by generation and analysis of gene regulatory networks.
Virus–Human Interactome Dataset
We obtained the protein interaction map of SARS-CoV-2 and human proteins from a recent report of affinity-purification mass spectrometry experiments on the proteins (Gordon et al., 2020). The interactome consists of 26 SARS-CoV-2 (out of 29) proteins interacting with 332 human proteins (Supplementary Table S1). The host–virus interactome is a network where proteins are nodes and interactions between proteins are edges.
Regulatory Interactome Analysis
For the identification of miRNA targeting host genes, its visualization and functional analysis were performed using the miRNet 1.0 server (Fan et al., 2016). The list of genes corresponding to the proteins identified to interact with the virus proteins was analyzed using miRNet to identify host miRNAs regulating the genes. A list of TFs interacting with the miRNAs was generated with the help of the miRNet and TransmiR v2.0 database (Tong et al., 2019). The hub genes were prioritized using network topological property, i.e., degree from the interaction network.
To identify miRNAs targeting SARS-CoV-2 genes, we performed miRNA target predictions using miRanda (version 3.3a) with an energy threshold of -20 kcal/mol (Betel et al., 2007). We downloaded a complete list of all the available mature human miRNAs from miRBase (Release 22.1) (Griffiths-Jones, 2006) and surveyed the literature to identify experimentally validated anti-viral miRNAs among the other miRNAs. To identify and compare miRNAs targeting other respiratory viruses (influenza A, measles, SARS-CoV related), we used miRTarP (Shu et al., 2018).
All human TFs were retrieved from the latest version of AnimalTFdb 3.0, which consists of 1,665 TFs classified into 73 TF families. Regulatory interactions (TF-target) from experimental evidences were extracted from TRRUST (version 2) database that contains 8,444 TF-target regulatory relationships of 800 human TFs (Han et al., 2018).
Gene Expression Analysis
To identify miRNAs expressed in the lung tissues, we used the Tissue Atlas database (Ludwig et al., 2016). For the gene expression analysis of host genes, we obtained two datasets, namely, a microarray dataset for SARS-CoV and an RNA-seq dataset for SARS-CoV-2, from the NCBI GEO database and ViPR host factor data search (Pickett et al., 2012). The first dataset SARS-CoV (GSE17400) (Yoshikawa et al., 2010) is derived from mock-infected Calu-3 subclone 2B4 cells, and SARS-CoV-2 (GSE147507) gene expression datasets are based on the samples derived from the infected primary human lung epithelium (NHBE) and transformed lung alveolar (A549) cells (Blanco-Melo et al., 2020a, b, p. 19). The microarray data of SARS-CoV were normalized by Robust Multi-chip Average (RMA) and SARS-CoV-2 RNA-seq data by DESeq2 R package.
Results
Regulatory Interactome of Host Genes
For the miRNA–host gene interactome analysis, we found that out of 332 host genes, 311 are found to be targeted by 2,197 human miRNAs (Supplementary Table S2). These genes are found to be involved in various biological processes, such as T-cell differentiation, onco-miRNAs, toxicity, regulation of Akt pathway, immune system, and others, as shown in Figure 3A. These miRNAs are distributed across 100 miRNA families, and the top 10 families among them are shown in Figure 3B. Further, from the degree prioritization, we identified the top five hub genes from the miRNA-target network as shown in Figure 4A of the miRNA-target network, and these are DCAF7, G3BP1, SBNO1, PRRC2B, and GGCX, which are reported to interact with SARS-CoV-2 NSP9, N, NSP12, and M proteins.
Figure 3. Gene set enrichment analysis showing (A) functional analysis of miRNAs and (B) miRNA family distribution.
Figure 4. (A) miRNA-target (host genes) interactome. (B) Subnetwork of reported anti-viral miRNA (n = 38) interaction network. (C) miRNA–TF interactome for host genes (green nodes represent TFs, blue nodes represent miRNAs, and yellow nodes represent target host genes). (D) GRN of TF targeting host genes.
Further expression analysis of the miRNA resulted in 49 miRNAs expressed in lung tissues. From the extensive literature survey, we identified 38 out of 2,197 miRNAs targeting 143 number of host genes with anti-viral activity depicted in Figure 4B and (Supplementary Table S3). Interestingly, six, namely, hsa-miR-1-3p, hsa-miR-17-5p, hsa-miR-199a-3p, hsa-miR-429, hsa-miR-15a-5p, and hsa-miR-20a-5p, were reported to play a role in respiratory diseases, such as adenovirus 2, influenza A, and RSV, and were downregulated in lung tissues during viral infection but overexpressed in normal lung tissues (Zhu et al., 2014; Leon-Icaza et al., 2019).
To understand miRNA–TF co-regulation of host genes that interact with SARS-CoV-2 proteins, we generated a miRNA–TF network that led to the identification of 51 miRNAs interacting with 77 TFs by either activation or repression as shown in Figure 4C. The hsa-mir-429 (reported to have anti-viral activity) and TP53 (a TF) were found to have a maximum degree of 28 and 12 in the miRNA–TF network, respectively. From the functional enrichment analysis, we identified the involvement of genes in various important biological pathways, including MAPK signaling, B-cell, Toll-like receptor signaling, and viral diseases, such as hepatitis C or measles.
Gene Regulatory Network Analysis
Out of 332 host genes, 81 genes are involved in the gene regulatory network (GRN) interacting with 38 TFs (Supplementary Table S4). After degree prioritization, we were able to identify the top five hub genes, namely, HMOX1, DNMT1, PLAT, GDF1, and ITGB1, directly targeting ORF3a and ORF8 among SARS-CoV-2 proteins as shown in Table 1. These are reported in modulating several virus activities, such as anti-viral immunity, epigenetic modification, and others (Danesh et al., 2011; Espinoza et al., 2017; Schäfer and Baric, 2017), as shown in Figure 4D.
Gene Expression Analysis
From the gene expression analysis of TF, we identified 148 TFs differentially regulated in SARS, whereas in the SARS-CoV-2 gene expression dataset, 48 and 52 TFs are expressed in A549 and NHBE cell lines. From the differentially regulated TFs in SARS-CoV-2 expression data, we found only two differentially expressed TFs, namely, STAT1 and STAT2, to target HMOX1 and SCRAB1 that interact with the ORF3a and NPS7 virus proteins (Supplementary Table S5).
Comparative Regulatory Interaction Analysis in Other Respiratory Diseases
We identified 1,018 miRNAs to target SARS-CoV-2 genes, among which 98 miRNAs are predicted to target SARS-CoV-2 genes uniquely when compared with other viruses studied here (Supplementary Table S6). 547 miRNAs targeting SARS-CoV-2 genes are found to be common to those in SARS-CoV, influenza A, and measles. Interestingly, we observed that 717 miRNAs targeting SARS-CoV-2 genes are common with those targeting SARS-CoV genes, possibly due to high genome similarity (Xu et al., 2020; Zhang and Holmes, 2020, p. 2) (Supplementary Table S7). The highest common host miRNAs targeting SARS-CoV-2 genes as well as other viruses in the study are those targeting influenza A virus (874 miRNAs), as shown in Table 2. miRNA interacting TFs are preferentially targeted by viral proteins; therefore, we further explored and identified miRNA–TF co-regulatory interactions. From the analysis, we identified 34, 41, 66, and 45 TFs interacting with 21, 30, 43, and 33 in SARS-CoV-2, SARS-CoV, influenza, and measles, respectively (Supplementary Table S8).
Discussion
Host-directed therapies is a recent approach targeting host cell factors that are required by a pathogen for replication (Kaufmann et al., 2018). Studies on SARS-CoV and Middle East respiratory syndrome (MERS) have already explored identifying potential host drug targets to block pathways and genes involved in coronavirus replication (Gassen et al., 2019). From the literature evidence, it is reported that TFs and miRNAs are two important master regulators controlling gene expression at the transcriptional and post-transcriptional levels. Studies also suggested that they play a role in multiple diseases by FFLs or FBLs (Zhang et al., 2015). Viruses tend to increase their genes expression by downregulating the host gene expression either co-transcriptionally in the nucleus or post-transcriptionally in the nucleus or cytoplasm. In herpesvirus, viral factors block the transcription initiation process by inhibiting TAF4. Targeting transcription process by SARS-CoV-2 could help in preventing the assembling of RNA polymerase II on host genes (Saçar Demirci and Adan, 2020). Genome-wide target predictions have previously shown that TFs are susceptible to regulation by miRNAs (Enright et al., 2003). Further, it is believed that miRNAs co-evolve with TFs, and the rapidly evolving TFs preferentially activate miRNAs (Chen and Rajewsky, 2007; Qiu et al., 2010). Thus, a direct relationship can be elucidated between the gene regulation networks controlled by TFs at the transcriptional level and those controlled by miRNAs at the post-transcriptional level. It has been suggested that miRNAs may provide genetic switch mechanisms to essentially inactivate the target genes by regulation of TF functioning and TF-mediated events (Chen et al., 2011). The active concert of these regulators in association to the regulatory FBLs orchestrates various cellular mechanisms. Understanding of miRNA–TF cross-talk during SARS-CoV-2 viral infection can be a more useful information to enhance our knowledge in this direction. The degree of expression of miRNAs is generally correlated in normal and diseased conditions. Many experiments have shown that during viral infection, the host–virus interactions involve overlapping target proteins or protein pairs among viral proteins and human miRNAs (Li et al., 2013). Therefore, the understanding of miRNA–TF interactions can be of relevance to SARS-CoV-2 research aimed at the development of novel therapeutics by modulation of specific TFs. In this study, the miRNA–TF network and the prioritized key interactions, generated on the basis of previous experiments and our predictions, provide useful information about such targetable interactions. From the miRNA–TF interactome, we identified 51 miRNAs interacting with 77 TFs involved in various biological pathways, such as MAPK signaling, B-cell, and Toll-like receptor signaling. It may be speculated that these interactions are also conserved in coronavirus disease 2019 (COVID-19) infection, used by SARS-CoV-2 to evade host immune responses.
During viral infections, anti-viral defense mechanism gets activated by the expression of interferons (IFNs) that are regulated by IFN-regulatory factors (IRF) family and STAT family TFs (Chiang and Liu, 2019). From the gene regulatory analysis, we identified STAT1 and STAT2 to be interacting with HMOX1 and SCRAB1 genes, which further interact with ORF3a and NPS7 in the SARS-CoV-2 PPI network. From the gene expression analysis, these TFs are found to be upregulated during SARS-CoV-2 infection at 24 h, confirming their role in anti-viral host defense by activation of many IFN stimulated genes (ISGs) as observed previously in SARS-CoV infection (Yoshikawa et al., 2010).
Mammalian cells inhibit virus infection by targeting transcripts with cellular miRNA (Song et al., 2010). The first evidence came from the observation that human miR-32 could limit the replication of primate foamy virus type 1 (PFV-1) in cells (Lecellier, 2005). Further, the liver-specific miR-122 was unexpectedly found to enhance the replication of the hepatitis C virus, whereas miR-199a-3p, miR-210, and miR-125a-5p were reported to suppress hepatitis B virus (HBV) replication (Jopling, 2005). In other virus systems, miR-101 was shown to suppress herpes simplex virus type 1 (HSV-1) propagation, whereas miR-3232, miR-491, and miR-654 were shown to inhibit the influenza virus (Lecellier, 2005; Potenza et al., 2011; Zheng et al., 2011).
In the present study, we generated and analyzed regulatory networks involving interacting miRNAs and TFs, targeting host genes based on information from a recently published host–virus interactome data. From the network analysis, we identified 2,197 miRNAs targeting host genes that interact with SARS-CoV-2 proteins. Among them, we identified 38 miRNAs targeting 143 host genes with reported anti-viral activity. Interestingly, six anti-viral miRNAs, namely, hsa-miR-1-3p, hsa-miR-17-5p, hsa-miR-199a-3p, hsa-miR-429, hsa-miR-15a-5p, and hsa-miR-20a-5p, are found to be downregulated post-infection in various viral respiratory diseases infecting the lungs, whereas in normal cases, these miRNAs are overexpressed. Clinical investigations have revealed that patients with cardiac diseases, hypertension, or diabetes who are administered with angiotensin-converting enzyme 2 (ACE2)-enhancing drugs including inhibitors and blockers exhibited elevated expression of ACE2; thus, the risk of getting the SARS-CoV-2 infection increased manifolds in such patients (Fang et al., 2020). Among the identified miRNAs, hsa-mir-9-5p was recently reported to target the 3’ UTR of ACE2 (Chen and Zhong, 2020), and hsa-mir-27b-3p plays a regulatory role in ACE2 signaling (Chen et al., 2015). Thus, a strong correlation between miRNA hsa-mir-27b-3p and ACE2 is established, but further experimental validation in SARS-CoV-2 infection may confirm this finding. Post-infection changes in the cellular miRNA expression profile are well reported (Triboulet et al., 2007; Cameron et al., 2008; Wang et al., 2008). These changes in cellular miRNA expression may be specifically induced by a given virus to create a more favorable intracellular environment for viral replication. On the contrary, these changes could also represent another aspect of the host cell innate immune response, triggered by a viral challenge and may inhibit virus replication (Wang et al., 2008).
To the best of our knowledge, our study is the first-ever attempt to identify key gene regulators and their co-regulation involving immune and MAPK signaling pathways, which may be explored for designing drug targets against SARS-CoV-2.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: GEO: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE17400 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE147507.
Author Contributions
DG and RS conceptualized the study and prepared the manuscript. RS, DG, and DS carried out the computational studies and performed the analysis. All the authors reviewed and approved the final version.
Funding
This work was financially supported by the Department of Biotechnology (DBT), Government of India, grant BT/BI/25/066/2012, awarded to DG. Financial support provided by the Indian Council of Medical Research (ICMR), India to RS as Senior Research Fellowship is duly acknowledged (2019-5850). DS received a fellowship from the Council of Scientific and Industrial Research [CSIR, India, 09/0512(0207)/2016/EMR-1].
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 ICGEB for providing the necessary infrastructure and facilities for the research.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.571274/full#supplementary-material
TABLE S1 | Host–SARS-CoV-2 interaction data.
TABLE S2 | List of miRNA targeting host genes identified using the miRnet server.
TABLE S3 | List of reported anti-viral miRna interactions.
TABLE S4 | GRN interactions.
TABLE S5 | Gene expression data of SARS-CoV and SARS-CoV-2.
TABLE S6 | miRNA identification in SARS-CoV-2 using miRanda.
TABLE S7 | miRNA identification in SARS-CoV, influenza, and measles from miRtarP.
TABLE S8 | miRNA–TF interaction analysis in respiratory viruses.
References
Ahmed, S. S. S. J., Paramasivam, P., Raj, K., Kumar, V., Murugesan, R., and Ramakrishnan, V. (2020). Regulatory cross talk between Sars-CoV-2 receptor binding and replication machinery in the human host. Front. Physiol. 11, 802. doi: 10.3389/fphys.2020.00802
Bernier, A., and Sagan, S. M. (2018). The diverse roles of microRNAs at the host-virus interface. Viruses 10:440. doi: 10.3390/v10080440
Betel, D., Wilson, M., Gabow, A., Marks, D. S., and Sander, C. (2007). The microRNA.org resource: targets and expression. Nucl. Acids Res. 36, D149–D153. doi: 10.1093/nar/gkm995
Blanco-Melo, D., Nilsson-Payant, B. E., Liu, W.-C., Møller, R., Panis, M., Sachs, D., et al. (2020a). SARS-CoV-2 launches a unique transcriptional signature from in vitro, ex vivo, and in vivo systems. bioRxiv[Preprint]. doi: 10.1101/2020.03.24.004655
Blanco-Melo, D., Nilsson-Payant, B. E., Liu, W.-C., Uhl, S., Hoagland, D., Møller, R., et al. (2020b). Imbalanced host response to SARS-CoV-2 drives development of COVID-19. Cell 181, 1036.e9–1045.e9. doi: 10.1016/j.cell.2020.04.026
Cameron, J. E., Fewell, C., Yin, Q., McBride, J., Wang, X., Lin, Z., et al. (2008). Epstein-Barr virus growth/latency III program alters cellular microRNA expression. Virology 382, 257–266. doi: 10.1016/j.virol.2008.09.018
Chen, C.-Y., Chen, S. T., Fuh, C. S., Juan, H.-F., and Huang, H. C. (2011). Coregulation of transcription factors and microRNAs in human transcriptional regulatory network. BMC Bioinformatics 12:S41. doi: 10.1186/1471-2105-12-S1-S41
Chen, K., and Rajewsky, N. (2007). The evolution of gene regulation by transcription factors and microRNAs. Nat. Rev. Genet. 8, 93–103. doi: 10.1038/nrg1990
Chen, L., and Zhong, L. (2020). Lung adenocarcinoma patients own higher risk of SARS-CoV-2 infection. Preprints 2020:2020020386.
Chen, L.-J., Xu, R., Yu, H.-M., Chang, Q., and Zhong, J.-C. (2015). The ACE2/Apelin Signaling, MicroRNAs, and Hypertension. Int. J. Hypertens. 2015, 1–6. doi: 10.1155/2015/896861
Chiang, H.-S., and Liu, H. M. (2019). The molecular basis of viral inhibition of IRF- and STAT-dependent immune responses. Front. Immunol. 9:3086. doi: 10.3389/fimmu.2018.03086
Danesh, A., Cameron, C. M., León, A. J., Ran, L., Xu, L., Fang, Y., et al. (2011). Early gene expression events in ferrets in response to SARS coronavirus infection versus direct interferon-alpha2b stimulation. Virology 409, 102–112. doi: 10.1016/j.virol.2010.10.002
Enright, A. J., John, B., Gaul, U., Tuschl, T., Sander, C., and Marks, D. S. (2003). MicroRNA targets in Drosophila. Genome Biol. 5:R1. doi: 10.1186/gb-2003-5-1-r1
Espinoza, J. A., González, P. A., and Kalergis, A. M. (2017). Modulation of antiviral immunity by heme oxygenase-1. Am. J. Pathol. 187, 487–493. doi: 10.1016/j.ajpath.2016.11.011
Fan, Y., Siklenka, K., Arora, S. K., Ribeiro, P., Kimmins, S., and Xia, J. (2016). miRNet - dissecting miRNA-target interactions and functional associations through network-based visual analysis. Nucl. Acids Res. 44, W135–W141. doi: 10.1093/nar/gkw288
Fang, L., Karakiulakis, G., and Roth, M. (2020). Are patients with hypertension and diabetes mellitus at increased risk for COVID-19 infection? Lancet Respir. Med. 8:e21.
Friedman, R. C., Farh, K. K.-H., Burge, C. B., and Bartel, D. P. (2009). Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 19, 92–105. doi: 10.1101/gr.082701.108
Fung, T. S., and Liu, D. X. (2019). Human coronavirus: host-pathogen interaction. Annu. Rev. Microbiol. 73, 529–557.
Gassen, N. C., Niemeyer, D., Muth, D., Corman, V. M., Martinelli, S., Gassen, A., et al. (2019). SKP2 attenuates autophagy through Beclin1-ubiquitination and its inhibition reduces MERS-coronavirus infection. Nat. Commun. 10: 5770.
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
Gordon, D. E., Jang, G. M., Bouhaddou, M., Xu, J., Obernier, K., White, K. M., et al. (2020). A SARS-CoV-2 protein interaction map reveals targets for drug repurposing. Nature 583, 459-468.
Griffiths-Jones, S. (2006). miRBase: microRNA sequences, targets and gene nomenclature. Nucl. Acids Res. 34, D140–D144. doi: 10.1093/nar/gkj112
Han, H., Cho, J.-W., Lee, S., Yun, A., Kim, H., Bae, D., et al. (2018). TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions. Nucl. Acids Res. 46, D380–D386. doi: 10.1093/nar/gkx1013
Jopling, C. L. (2005). Modulation of hepatitis C Virus RNA abundance by a liver-specific MicroRNA. Science 309, 1577–1581. doi: 10.1126/science.1113329
Kaufmann, S. H. E., Dorhoi, A., Hotchkiss, R. S., and Bartenschlager, R. (2018). Host-directed therapies for bacterial and viral infections. Nat. Rev. Drug Discov. 17, 35–56. doi: 10.1038/nrd.2017.162
Khailany, R. A., Safdar, M., and Ozaslan, M. (2020). Genomic characterization of a novel SARS-CoV-2. Gene Rep. 19:100682. doi: 10.1016/j.genrep.2020.100682
Khan, M. A. A. K., Sany, M. R. U., Islam, M. S., and Islam, A. B. (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
Lecellier, C.-H. (2005). A cellular microRNA mediates antiviral defense in human cells. Science 308, 557–560. doi: 10.1126/science.1108784
Leon-Icaza, S. A., Zeng, M., and Rosas-Taraco, A. G. (2019). microRNAs in viral acute respiratory infections: immune regulation, biomarkers, therapy, and vaccines. ExRNA 1:1.
Li, Z., Cui, X., Li, F., Li, P., Ni, M., Wang, S., et al. (2013). Exploring the role of human miRNAs in virus–host interactions using systematic overlap analysis. Bioinformatics 29, 2375–2379. doi: 10.1093/bioinformatics/btt391
Lou, W., Liu, J., Ding, B., Chen, D., Xu, L., Ding, J., et al. (2019). Identification of potential miRNA-mRNA regulatory network contributing to pathogenesis of HBV-related HCC. J. Transl. Med. 17:7.
Ludwig, N., Leidinger, P., Becker, K., Backes, C., Fehlmann, T., Pallasch, C., et al. (2016). Distribution of miRNA expression across human tissues. Nucleic Acids Res. 44, 3865–3877. doi: 10.1093/nar/gkw116
Mohamed, R. H., Abu-Shahba, N., Mahmoud, M., Abdelfattah, A. M. H., Zakaria, W., and ElHefnawi, M. (2019). Co-regulatory Network of Oncosuppressor miRNAs and Transcription Factors for Pathology of Human Hepatic Cancer Stem Cells (HCSC). Sci. Rep. 9:5564.
Nakhaei, P., Genin, P., Civas, A., and Hiscott, J. (2009). RIG-I-like receptors: sensing and responding to RNA virus infection. Sem. Immunol. 21, 215–222. doi: 10.1016/j.smim.2009.05.001
Pickett, B. E., Sadat, E. L., Zhang, Y., Noronha, J. M., Squires, R. B., Hunt, V., et al. (2012). ViPR: an open bioinformatics database and analysis resource for virology research. Nucleic Acids Res. 40, D593–D598. doi: 10.1093/nar/gkr859
Potenza, N., Papa, U., Mosca, N., Zerbini, F., Nobile, V., and Russo, A. (2011). Human microRNA hsa-miR-125a-5p interferes with expression of hepatitis B virus surface antigen. Nucl. Acids Res. 39, 5157–5163. doi: 10.1093/nar/gkr067
Qiu, C., Wang, J., Yao, P., Wang, E., and Cui, Q. (2010). microRNA evolution in a human transcription factor and microRNA regulatory network. BMC Syst. Biol. 4:90. doi: 10.1186/1752-0509-4-90
Saçar Demirci, M. D., and Adan, A. (2020). Computational analysis of microRNA-mediated interactions in SARS-CoV-2 infection. PeerJ 8:e9369. doi: 10.7717/peerj.9369
Schäfer, A., and Baric, R. S. (2017). Epigenetic landscape during coronavirus infection. Pathogens 6:8. doi: 10.3390/pathogens6010008
Shu, X., Zang, X., Liu, X., Yang, J., and Wang, J. (2018). Predicting MicroRNA mediated gene regulation between human and viruses. Cells 7:100. doi: 10.3390/cells7080100
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. JVI 84, 8849–8860. doi: 10.1128/JVI.00456
Tong, Z., Cui, Q., Wang, J., and Zhou, Y. (2019). TransmiR v2.0: an updated transcription factor-microRNA regulation database. Nucl. Acids Res. 47, D253–D258. doi: 10.1093/nar/gky1023
Triboulet, R., Mari, B., Lin, Y.-L., Chable-Bessia, C., Bennasser, Y., Lebrigand, K., et al. (2007). Suppression of microRNA-silencing pathway by HIV-1 during virus replication. Science 315, 1579–1582. doi: 10.1126/science.1136319
Tsang, J., Zhu, J., and van Oudenaarden, A. (2007). MicroRNA-mediated feedback and feedforward loops are recurrent network motifs in mammals. Mol. Cell 26, 753–767. doi: 10.1016/j.molcel.2007.05.018
Wang, F.-Z., Weber, F., Croce, C., Liu, C.-G., Liao, X., and Pellett, P. E. (2008). Human cytomegalovirus infection alters the expression of cellular microRNA species that affect its replication. J. Virol. 82, 9065–9074.
Wang, H., Li, X., Li, T., Zhang, S., Wang, L., Wu, X., et al. (2020). The genetic sequence, origin, and diagnosis of SARS-CoV-2. Eur. J. Clin. Microbiol. Infect. Dis. 39, 1629–1635. doi: 10.1007/s10096-020-03899-4
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. doi: 10.3390/v12020244
Yoshikawa, T., Hill, T. E., Yoshikawa, N., Popov, V. L., Galindo, C. L., Garner, H. R., et al. (2010). Dynamic innate immune responses of human bronchial epithelial cells to severe acute respiratory syndrome-associated coronavirus infection. PLoS one 5:e8729. doi: 10.1371/journal.pone.0008729
Zhang, H.-M., Kuang, S., Xiong, X., Gao, T., Liu, C., and Guo, A.-Y. (2015). Transcription factor and microRNA co-regulatory loops: important regulatory motifs in biological processes and diseases. Brief. Bioinform. 16, 45–58. doi: 10.1093/bib/bbt085
Zhang, Y.-Z., and Holmes, E. C. (2020). A genomic perspective on the origin and emergence of SARS-CoV-2. Cell 181, 223–227. doi: 10.1016/j.cell.2020.03.035
Zheng, S., Li, Y., Zhang, Y., Li, X., and Tang, H. (2011). MiR-101 regulates HSV-1 replication by targeting ATP5B. Anti. Res. 89, 219–226. doi: 10.1016/j.antiviral.2011.01.008
Keywords: regulatory network, miRNA, SARS-CoV-2, hub genes, TFS
Citation: Sardar R, Satish D and Gupta D (2020) Identification of Novel SARS-CoV-2 Drug Targets by Host MicroRNAs and Transcription Factors Co-regulatory Interaction Network Analysis. Front. Genet. 11:571274. doi: 10.3389/fgene.2020.571274
Received: 10 June 2020; Accepted: 25 August 2020;
Published: 14 October 2020.
Edited by:
Chen Liang, McGill University, CanadaReviewed by:
Ruby C. Y. Lin, Westmead Institute for Medical Research, AustraliaSonia Navas-Martin, Drexel University, United States
Copyright © 2020 Sardar, Satish and Gupta. 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: Dinesh Gupta, dinesh@icgeb.res.in