In silico analysis of SARS-CoV-2 genomes: Insights from SARS encoded non-coding RNAs
Neha Periwal
Urvashi Bhardwaj
Sankritya Sarma2 
Pooja Arora
Vikas Sood- 1Department of Biochemistry, Jamia Hamdard, New Delhi, India
- 2Department of Zoology, Hansraj College, University of Delhi, Delhi, India
The recent pandemic caused by Severe Acute Respiratory Syndrome Coronavirus-2 has resulted in enormous deaths around the world. Clues from genomic sequences of parent and their mutants can be obtained to understand the evolving pathogenesis of this virus. Apart from the viral proteins, virus-encoded microRNAs (miRNAs) have been shown to play a vital role in regulating viral pathogenesis. Thus we sought to investigate the miRNAs encoded by SARS-CoV-2, its mutants, and the host. Here, we present the results obtained using a dual approach i.e (i) identifying host-encoded miRNAs that might regulate viral pathogenesis and (ii) identifying viral-encoded miRNAs that might regulate host cell signaling pathways and aid in viral pathogenesis. Analysis utilizing the first approach resulted in the identification of ten host-encoded miRNAs that could target the SARS, SARS-CoV-2, and its mutants. Interestingly our analysis revealed that there is a significantly higher number of host miRNAs that could target the SARS-CoV-2 genome as compared to the SARS reference genome. Results from the second approach resulted in the identification of a set of virus-encoded miRNAs which might regulate host signaling pathways. Our analysis further identified a similar “GA” rich motif in the SARS-CoV-2 and its mutant genomes that was shown to play a vital role in lung pathogenesis during severe SARS infections. In summary, we have identified human and virus-encoded miRNAs that might regulate the pathogenesis of SARS coronaviruses and describe similar non-coding RNA sequences in SARS-CoV-2 that were shown to regulate SARS-induced lung pathology in mice.
Introduction
Coronaviruses are a group of enveloped viruses with positive-sense single-stranded RNA as the genetic material. These viruses cause mild to severe respiratory and intestinal infections in immune-compromised individuals. They were not considered to be pathogenic until the outbreak of Severe Acute Respiratory Syndrome (SARS) in China and the Middle East Respiratory Syndrome Coronavirus (MERS-CoV) in Middle East countries in 2002 and 2012 respectively when the viral infection led to significant deaths among the infected individual’s (Song et al., 2019). As on 24 September 2022, the SARS-CoV-2 pandemic has resulted in more than 611 million infections resulting in 6.5 million deaths globally (https://covid19.who.int/). Data from the ongoing pandemic suggest that mortality caused by this virus has already surpassed that of combined mortality caused by both SARS and MERS coronaviruses. Though SARS-CoV-2 has caused significant mortality and morbidity, its evolution trajectory points towards an increased transmissibility (Volz et al., 2021). SARS-CoV-2 has a slow evolutionary rate i.e. 1 or 2 nucleotide alterations/month/lineage in the viral genome which consists of around 30,000 base pairs (Boehm et al., 2021). Despite this slow mutation rate, SARS-CoV-2 has evolved into around sixteen lineages in a short time since the onset of the pandemic (Tegally et al., 2021). The underlying mutations in these lineages allow the virus to transmit efficiently thereby enhancing the viral fitness. The random mutations in the viral genomes have led to changes in the viral proteins resulting in improved functions (Korber et al., 2020). However, apart from modulating the functions of the viral proteins, the mutations in the viral genomes can lead to changes in the microRNA binding sites thereby resulting in the appearance of the escape mutants (Hosseini Rad and McLellan, 2020).
MicroRNAs (miRNAs) are non-coding RNA molecules that are ~21 nucleotides long and can regulate gene expression primarily via imperfect base pairing with the 3’ untranslated regions of the target RNA molecules (Duchaine and Fabian, 2019; Liu and Wang, 2019). Apart from regulating the biological pathways, miRNAs have been shown to modulate viral infections (Mishra et al., 2020). Following viral infections, innate immune responses are induced leading to the generation of antiviral responses. This arm of antiviral signaling constitutes several potential antiviral factors including antiviral miRNAs. These miRNA molecules can lead to the inhibition of viral replication. For instance, miR-196, miR-351, and miR-448 have been shown to inhibit the Hepatitis C Virus infection (Pedersen et al., 2007). The miR-155 has been shown to regulate type I interferon signaling thereby modulating the replication of several viruses including HIV-1 (Swaminathan et al., 2012), Epstein-Barr virus (Jiang et al., 2006), reticuloendotheliosis virus (Bolisetty et al., 2009), Dengue virus (Su et al., 2020), West-Nile virus (Natekar et al., 2019), Japanese encephalitis virus (Thounaojam et al., 2014; Pareek et al., 2014) and Herpes Simplex Virus 1 (Wang et al., 2019). Conversely, viruses have also been shown to encode small non-coding RNA molecules to facilitate viral replication. It has been shown that several miRNAs encoded by Kaposi’s Sarcoma-associated Herpesvirus help in maintaining the viral latency (Lei et al., 2012). Several RNA viruses including Influenza virus (Li et al., 2018), HIV-1 (Zhang et al., 2014; Li et al., 2016), Dengue virus (Hussain and Asgari, 2014), and Ebola virus (Prasad et al., 2020) are also shown to encode small non-coding RNAs which can modulate the viral infection. Additionally, a complex cross-talk among miRNA and host factors further contributes to the replication of the viruses in their host (Bruscella et al., 2017; Mishra et al., 2020).
The ongoing SARS-CoV-2 pandemic has driven several studies to identify miRNAs as alternative therapeutics against the virus (Panda et al., 2022). Studying miRNAs has become crucial for the development of possible diagnostics, prognostics and therapeutics that could help in treatment strategies in the COVID-19 (Fani et al., 2021; Ying et al., 2021; Ahmed et al., 2022). Several studies aiming to predict miRNAs regulating SARS-CoV-2 pathogenesis were published and various miRNAs are identified that might have an impact on the replication of SARS-CoV-2 (Chen and Zhong, 2020; Marchi et al., 2021; Kucher et al., 2022; Morales et al., 2022; Yang et al., 2022). Apart from directly targeting the viral genomes, the miRNAs can indirectly regulate viral replication via regulating host factors that are critical for viral replication. The expression of proteins implicated in the SARS-CoV-2 life cycle, including ACE2, TMPRSS2, Spike proteins, and Nsp12, can be inhibited by miRNAs (Yang et al., 2022; Paul et al., 2022). The type I Interferon pathway is inhibited by SARS-CoV-2 encoded miRNAs, which also regulate the allelic differential expression of vulnerable genes (Zhu et al., 2021). Additionally, it has been shown that miRNAs can regulate SARS-CoV-2 replication via modulating certain receptors including ACE2, AXL, ADAM17, HAT1, NRP1, VMP1, Cyclophilin A, and Cathepsin B (Zhang et al., 2021; Zeng et al., 2022).
Numerous host- and virus-specific indicators linked to SARS CoV-2 infection have recently emerged (Sevgin and Sevgin, 2021). Host miRNAs were found to target ACE2 and TMPRSS2 genes which regulate SARS-Cov-2 cellular entry (Kaur et al., 2021; Bellae Papannarao et al., 2022) and can also be modified to act in a pro-viral manner by blocking host immune systems (Sevgin and Sevgin, 2021). Recently several host-encoded miRNAs were found to regulate ORF1a/b which is necessary for viral replication and translation (Arghiani et al., 2021). A study on SARS coronavirus successfully characterized the role of viral encoded non-coding RNAs in the virus-induced lung pathology (Morales et al., 2017). Furthermore, phylogenetic and taxonomy studies have confirmed that SARS-CoV-2 forms a sister clade with SARS coronavirus suggesting that rich information on SARS coronavirus can be leveraged to gain insights into SARS-CoV-2 (The species Severe acute respiratory syndrome-related coronavirus, 2020). It was observed that inhibiting viral encoded non-coding RNAs led to a reduction in the lung pathology of virus-infected mice. Taking into consideration the significant similarity between SARS and SARS-CoV-2 and the role of non-coding RNAs in virus pathogenesis and cell signaling pathways, we aimed to study host and viruses encoded miRNAs and compare them among SARS, SARS-CoV-2, and the mutants to understand the pathogenesis of these viruses.
In this study, we utilized a suite of highly popular algorithms to identify host and virus-encoded miRNAs that might have the potential to regulate SARS-CoV-2 pathogenesis. Moreover, we compared our results along with data from other published studies to identify a common set of miRNA candidates. These high-confidence candidates can then be validated in the wet lab as well as provide us with the opportunity to design broadly acting drug molecules against SARS-CoV-2 and its mutants.
Methods
Identification of host miRNAs targeting several SARS-CoV-2 and SARS reference genomes
MicroRNAs are known to regulate diverse biological pathways including various steps in viral replication (Bernier and Sagan, 2018). Hence, we studied the role of host miRNAs in regulating viral replication by predicting the host miRNA binding sites among the genome of these viruses. In order to gain an in-depth understanding of the pathogenesis of SARS-CoV-2, we analyzed the reference genome of SARS-CoV-2 along with its variants namely Alpha (B1.1.7), Delta (B.1.617.2, AY.4.2), Beta (B.1.351), Omicron (B.1.1529, BA.1, BA.3, BA.1.1), Gamma lineage, Epsilon (B.1.27), Eta (B.1.525), Kappa (B.1.525), Lambda (C.37), Mu (B.1.627) and Theta (P.3) that were obtained from NCBI (Brister et al., 2015) and GISAID (Khare et al., 2021) (Supplementary table 1). The SARS reference genome was used for the comparison and the analysis was performed on both the forward and reverse sequences of viral genomes. For the prediction of host miRNAs that could bind to viral genomes, the miRDB algorithm was used (Chen and Wang, 2020). The analysis was performed in the custom mode where the user-defined sequences could be uploaded and processed. All the analyses were done on default settings and only those miRNAs were studied further that had a target score >= 80 so as to follow stringent selection criteria. Thus we obtained a list of miRNAs that could potentially target SARS-CoV-2 and its mutants. Recently, several studies were published to predict host-encoded miRNAs that might regulate SARS-CoV-2. Hence, we present the common miRNAs that are predicted by multiple studies involving diverse set of prediction algorithms.
Identification of SARS-CoV-2 and SARS encoded miRNAs
Several viruses have been shown to regulate host signaling pathways by encoding miRNAs in their genomes. However, there have been debates as to whether RNA viruses could encode miRNAs or not. Recent studies have successfully identified several miRNAs that are encoded by various RNA viruses including HIV-1 (Bennasser et al., 2004), West Nile Virus (Hussain et al., 2012), Dengue (Hussain and Asgari, 2014), Influenza (Li et al., 2018) and Ebola virus (Prasad et al., 2020). It was also observed that genomic sequences of these viruses could fold into stem-loop structures which gave rise to miRNAs (Klaver and Berkhout, 1994). Hence, in order to gain insights into various miRNAs possibly encoded by genomes of SARS coronavirus, SARS-CoV-2, and its variants, a set of computational tools was used to predict miRNAs encoded by these viruses. Since hairpin loops have been shown to produce miRNAs in RNA viruses, we used the miRNAFold algorithm to predict hairpin loops in viral genomes (Tempel and Tahi, 2012). As all hairpin loops do not result in functional miRNAs, hence output file from the miRNAFold algorithm was further processed using miRBoost algorithm to identify whether the predicted hairpin structures could form microRNAs or not (Tran et al., 2015). Once the viral encoded miRNAs were identified, the results were then processed using miRBase to find any similarities among the virus-encoded miRNAs and human miRNAs (Kozomara et al., 2019).
Comparison of SARS encoded non-coding RNA sequences with SARS-CoV-2 genomes
Recently, deep sequencing based approaches identified SARS encoded non-coding RNAs which played a critical role in virus-induced lung pathology (Morales et al., 2017). In order to gain an understanding of the genomic regions responsible for the generation of non-coding RNAs, we aligned SARS-CoV-2 sequences along with the mutants SARS reference sequence using the Clustal omega (Sievers et al., 2011). Aligned SARS genomic sequences that were shown to encode non-coding RNAs were then further studied to identify the percentage homology with that of SARS-CoV-2 and its mutants.
Results
Identification of host miRNAs targeting several SARS, SARS-CoV-2 and its variants
We used miRDB algorithm to gain insights into host encoded miRNAs (Homo sapiens) that could target viral genomes. The analysis was performed on both the forward and reverse genomic sequences of all the genomes that were being studied. Our analysis identified several host miRNAs that could target each of the forward and reverse genomic sequences of SARS, SARS-CoV-2 and its variants (Supplementary Table 2). Figure 1A illustrates the common host miRNAs (n=354) that were predicted in 45% of the genomes used in this study. Similar analysis of common host miRNAs targeting 80% of the genomes identified 54 miRNAs (Figure 1B). A comparison of all the forward and reverse sequences of SARS, SARS-CoV-2 and its variants revealed 10 host miRNAs that could target all the genomes being studied and hence could be potentially used as broad-acting antivirals (Figure 1C). Out of these 10 miRNAs, three of them including hsa-miR-4282, hsa-miR-33a-3p and hsa-miR-4775 were found in several other studies. The hsa-miR-4282 has been shown to target ATF2 gene (Vadivalagan et al., 2022) that arrests the viral infection, hsa-miR-33a-3p targets SMAD4 which is involved in the Wnt signalling pathway (Yousefi et al., 2020), similarly hsa-miR-4775 targets SMAD7 and it plays role in TGF-beta signaling pathway (Yousefi et al., 2020). The hsa-miR-497-3p, inhibits vascular endothelial growth factor A, suppresses angiogenesis, and is found to be involved in the body’s defense mechanisms (Van Laar et al., 2018). Another miRNA, hsa-miR-196a-1-3p was predicted to be able to bind to SARS-CoV, MERS-CoV, and SARS-CoV-2 genomes (Kim et al., 2020). The hsa-miR-548t-3p was found to be down-regulated in the ischemic stroke (Li et al., 2022) whereas hsa-miR-33a-3p was shown to have ~90% sequence identity with VARV and VACV viral genomes. It was observed that hsa-miR-33a-3p targets A46L and A41L genes respectively thereby modulating the host defense mechanisms (Hasan et al., 2016). Recently, other groups have also utilized a diverse set of algorithms to identify the miRNAs that could target SARS-CoV-2 genomes (Mosharaf et al., 2022; Samy et al., 2022; Shirvani et al., 2022; Zeng et al., 2022). In order to further advance our understanding of host miRNAs targeting SARS-CoV-2 genomes, we compared our data with other studies. We observed several host-encoded miRNAs including hsa-miR-3529-3p, has-miR-516b-5p, and hsa-miR-155-5p that were common in many studies pointing toward the potential miRNAs that might be the strong candidates for further wet-lab validation and further understanding how these miRNAs are contributing towards the pathogenesis of SARS-CoV-2.
Figure 1 Shows the host-encoded miRNAs that could target SARS, SARS-CoV-2, and its variants. (A) Common miRNAs that could target 45% of genomes from the study. (B) Host miRNAs that could target 80% of the genomes from the study and (C) Top 10 miRNAs targeting all the forward and reverse genomes used in this study.
Identification of miRNAs encoded by SARS, SARS-CoV-2 and its variants
Viral encoded miRNAs have been shown to regulate cellular signaling pathways. RNA viruses are known to produce miRNAs that are critical for their replication cycle (Li and Zou, 2019). In order to identify miRNAs encoded by SARS, SARS-CoV-2, and its mutants, we utilized a suite of algorithms that could identify stem-loop structures in the viral genomes. This was followed by the prediction of miRNAs among these stem-loop structures and finally comparing the predicted miRNAs with that of human-encoded miRNAs. As mentioned earlier, the miRNAFold algorithm was used to predict the stem loop structures in the viral genomes. The algorithm predicted various hairpin loops in genomes of SARS, SARS-CoV-2, and its variants. We used both the forward and the reverse strand of viral RNA to predict the stem-loop structures. Since all the predicted stem-loop structures do not form miRNAs, we used the miRBoost algorithm to identify potential stem-loop structures that could form miRNAs. The results thus obtained were then processed with the miRBase algorithm to investigate any similarity among virus and host-encoded miRNAs. We used this step with the rationale that any similarity among viral and host-encoded miRNAs will lead to the modulation of similar signaling pathways. The miRbase algorithm identified several viral encoded miRNAs that had similarities with human encoded miRNAs (Supplementary Table 3). Figure 2A represents the miRNAs that were predicted in 45% of the genomes used in this study. Figure 2B presents the top 10 miRNAs that can target various coronavirus genomes used in this study. Comparative analysis of miRNAs predicted from all the genomes used in this study revealed that hsa-miR-4471 was present in 33 out of 34 genomes. This miRNA was predicted to be encoded by SARS-CoV-2 and all its mutants including the forward and the reverse genomic sequences. It was further observed that hsa-miR-4471 was encoded by only the reverse sequence of SARS coronavirus but not by the forward sequence.
Figure 2 Showing the viral encoded miRNAs that had similarity with that of miRNAs encoded by the H sapiens. (A) Heatmap showing miRNAs encoded by 45% of the genomes under study. (B) Top 10 miRNAs that could target various genomes used in the study. (C) Common miRNAs among the forward and the reverse genomic sequences of SARS, SARS-CoV-2 and its variants. (D) Top 10 miRNAs targeting forward and the reverse genomes used in the study.
A literature survey for the top ten miRNAs targeting SARS, SARS-CoV-2, and its mutants revealed that some of the miRNAs have been shown to play a key role in coronavirus pathogenesis. For instance, hsa-miR-3529-5p could modulate SARS-CoV-2 by targeting NF-kb pathway whereas hsa-miR-7-5p was shown to modulate cellular apoptosis via mTOR pathway (Abedi et al., 2021), hsa-miR-145-3p could regulate SARS-CoV-2 via interfering with SMAD3 levels and can also target 5’ end of the Ebola virus genome (Golkar et al., 2016) and hsa-582-5p could target ACE-2 which is an important receptor of SARS-CoV-2 (Kim et al., 2020). Comparative analysis of miRNAs targeting all the forward and reverse genomic sequences led us to the identification of hsa-miR-4502 as a miRNA common in all the forward genomic sequences of SARS, SARS-CoV-2, and its mutants. Since these miRNAs are encoded by viral genomes, the prediction of a single miRNA in all the forward genomes of SARS-SARS-CoV-2 and its mutants points toward the importance of the putative miRNA. Our results are in agreement with the previous study that identified similarities among SARS-CoV-2 encoded pre-miRNAs and human encoded hsa-miR-4502 (Aydemir et al., 2021). Analysis of all the reverse genomic sequences of SARS, SARS-CoV-2 and its mutants led us to the identification of three miRNAs (hsa-miR-8088, hsa-miR-11400, hsa-miR-4471) were found common in all reverse sequences (Figure 2C). The presence of these miRNAs in the reverse sequences of all the viral genomes under study points toward their importance in the viral life cycle.
Then we compared the miRNAs targeting forward and reverse genomes of all the SARS-CoV-2 variants used in the study (Figure 2D). It was observed that some of the predicted miRNAs could target cellular host factors thereby facilitating the replication of SARS-CoV-2. Some of the miRNAs targeting forward genomes of SARS-CoV-2 variants included hsa-miR-3672, hsa-miR-7161-5p, hsa-miR-30b-5p and hsa-miR-3682-5p were have been shown to regulate host factors including ICOSLG, G3BP1/P2, CREB1 respectively (Vastrad et al., 2020; Almutairy et al., 2021; Vadivalagan et al., 2022). Analysis of reverse sequences of SARS-CoV-2 mutants revealed that has-miR-548 could regulate G3BP1/P2 thereby contributing to viral replication (Almutairy et al., 2021).
Comparison of SARS encoded non-coding RNA sequences with SARS-CoV-2 and its mutants
A recent study had identified SARS encoded non-coding RNAs that contributed towards the pathology of the virus (Morales et al., 2017). In an attempt to understand the pathogenesis of SARS-CoV-2 and its variants, we compared its genomic sequence with that of SARS to estimate the sequence similarity among both the viruses in genomic regions that resulted in non-coding RNAs. We observed a very similar “GA” rich region SARS-CoV-2 including its variants that was reported to be enriched in several SARS encoded non-coding RNAs. This region was found to be highly conserved in all the variants except the omicron variant (BA.1, BA.3, BA.1.1, B.1.1.529) which had a deletion in N-9b domain from position 28038-28047 [Figure 3A(ix)]. Furthermore, comparative analysis for other non-coding RNAs revealed that SARS-CoV-2 also had similar stretches albeit with some nucleotide changes (Figure 3A). Figures 3B, C represents the SARS-CoV-2 variants that were used for the comparison. Interestingly genomic regions of SARS-CoV-2 had more than 85% similarity with SARS encoded scRNA-N whose inhibition led to a reduction in SARS-induced lung pathology. Apart from scRNA-N, SARS-CoV-2 showed a high degree of similarity (>80-90%) with several SARS encoded non-coding RNAs from regions Nsp3.3, Nsp13, Membrane protein, 7a, N-9b, N-TDN, and Nucleocapsid protein. This high level of similarity among SARS coronavirus encoded non-coding RNAs with that of SARS-CoV-2 genomes suggest that the virus might also encode similar non-coding RNAs that have been shown to contribute towards lung pathology.
Figure 3 Showing the alignment of genomic regions of SARS-CoV-2 and its mutants with those regions of SARS that were shown to produce non-coding RNAs. (A) (i-xi)) represents the genomic regions where SARS encoded ncRNAs were identified (Morales et al., 2017). (B, C) Represent the sequences of SARS-CoV-2 variants that were used for the comparison.
Discussion
The recent pandemic of COVID-19 has caused significant mortality around the world. Since its origin in Wuhan, China SARS-CoV-2 has been able to spread throughout the world in a very short duration. The efforts of the scientific community to design effective drugs against the virus are hindered by the emergence of SARS-CoV-2 mutants. The virus has evolved into more than sixteen lineages in a very short span of time. Therefore, understanding the evolution of the virus is of utmost importance. Scientists are racing against the time to find an effective cure for the virus. Apart from regulating vital cellular processes (Andergassen and Rinn, 2022), non-coding RNA has been shown to regulate viral pathogenesis (Liu et al., 2019) and they act as a negative regulator and inhibit antiviral responses (Kesheh et al., 2022). It has also been stipulated that the SARS-CoV-2 genome has sequence similarity with the host non-coding RNAs thereby resulting in the possible epigenetic crosstalk among the virus and its host (Talotta et al., 2022). lncRNAs that were found to be functionally linked and implicated in the innate immune response are upregulated as a result of SARS-CoV-2 infection (Enguita et al., 2022). They were found to alter the host’s immune response (Paul et al., 2022). Therefore, the present study was designed to understand the dynamics of non-coding RNA among SARS, SARS-CoV-2, and its mutants. We used a suite of highly popular algorithms to identify host and viral encoded miRNAs that might affect viral replication in host cells. Owing to the significant similarity of SARS-CoV-2 with that of SARS, we sought to leverage the knowledge generated on SARS genomic sequences and apply it to SARS-CoV-2 and its mutant sequences to identify possible overlaps among both the viruses as well as understand ncRNA dynamics as the virus evolves. We compared several SARS-CoV-2 viral genomic sequences and their mutants with that of the SARS reference sequence in order to understand how they might be targeted by host-encoded miRNAs. Interestingly our data suggested that SARS-CoV-2 genomic sequences were targeted by significantly more numbers of host-encoded miRNAs as compared to the SARS reference genome. This might provide a possible explanation as to why the virus causes mild disease in healthy individuals. Moreover, our comparative analysis of our candidates with that of other published candidates led to the identification of ten host-encoded miRNAs that were identified in many studies. Further wet lab validations are required for the proper validation of these potential candidates. These viral encoded miRNAs have been shown to arrest the transition of phases in cell cycle by preventing the replication of specific genes thereby contributing towards the viral replication. For instance, Let-7c-3p which is one of the top ten miRNAs from our study has been shown to be associated with SARS-CoV-2 infection and pathogenesis (Jafarinejad-Farsangi et al., 2020).
Viruses are known to encode miRNAs to dampen the host responses. We utilized a suite of algorithms to identify possible stem-loop structures in viral genomes which could have the potential to be converted into viral encoded miRNAs. Our tools successfully identified several virus-encoded miRNAs that had significant similarities with host miRNAs. Viral encoded miRNAs having very similar sequences to host-encoded miRNAs might be an excellent strategy to regulate cell signaling pathways. Though we identified several miRNAs encoded by SARS and SARS-CoV-2, despite the high degree of genomic similarity among SARS and SARS-CoV-2, we could not identify common miRNAs among both of them.
Recently, SARS-encoded non-coding RNAs were identified in the lungs of infected mice. These non-coding RNAs were distributed along the entire viral genome and antagomir-based inhibition of these viral encoded RNAs resulted in reduced lung pathology. This data indicate an important role of these viral encoded RNAs during viral pathogenesis. Therefore, we aligned all the SARS-CoV-2 and mutant genomic sequences with the SARS reference genome to investigate whether non-coding RNA sequences among both the viruses are conserved or not. We observed that there was a very similar “GA” rich region in all SARS-CoV-2 genomes including the mutants, that was identified as a part of SARS encoded non-coding RNAs. We also found a very high similarity among SARS encoded non-coding RNAs and all SARS-CoV-2 and mutant genomic sequences suggesting a possible role of these non-coding RNAs in SARS-CoV-2 pathogenesis.
To conclude, we have identified potential host encoded miRNAs that could target genomic sequences of SARS, SARS-CoV-2 and its mutants. We have also identified virus-encoded non-coding RNAs that might modulate the pathogenesis of SARS-CoV-2 thereby identifying new molecular targets against the virus.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
VS and PA conceived the idea, designed the experiments, supervised the study wrote and finalized the manuscript. NP, UB, and SS performed the analysis. All authors contributed to the article and approved the submitted version.
Funding
NP is supported by UGC fellowship. VS is supported by UGC under FRP programme. The study was funded by startup grant to VS by UGS and DST PURSE grant.
Acknowledgments
We are thankful to all those scientists who contributed to SARS-CoV-2 genomic sequences.
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2022.966870/full#supplementary-material
References
Abedi, F., Rezaee, R., Hayes, A.W., Nasiripour, S., Karimi, G. (2021). MicroRNAs and SARS-CoV-2 life cycle, pathogenesis, and mutations: biomarkers or therapeutic agents? Cell Cycle 20, 143–153. doi: 10.1080/15384101.2020.1867792
Ahmed, J. Q., Maulud, S. Q., Dhawan, M., Priyanka, , Choudhary, P. O., Jalal, J. P., et al. (2022). MicroRNAs in the development of potential therapeutic targets against COVID-19: A narrative review. J. Infection Public Health 15 (7), 788–799. doi: 10.1016/j.jiph.2022.06.012
Almutairy, B. K., Alshetaili, A., Anwer Khalid, M. N. A. (2021). In silico identification of MicroRNAs targeting the key nucleator of stress granules, G3BP: Promising therapeutics for SARS-CoV-2 infection. Saudi J. Biol. Sci. 28, 7499–7504. doi: 10.1016/j.sjbs.2021.08.056
Andergassen, D., Rinn, J. L. (2022). From genotype to phenotype: genetics of mammalian long non-coding RNAs in vivo. Nat. Rev. Genet. 23, 229–243. doi: 10.1038/s41576-021-00427-8
Arghiani, N., Nissan, T., Matin, M. M. (2021). Role of microRNAs in COVID-19 with implications for therapeutics. Biomed. Pharmacother. 144, 112247. doi: 10.1016/j.biopha.2021.112247
Aydemir, M. N., Aydemir, H. B., Korkmaz, E. M., Budak, M., Cekin, N., Pinarbasi, E, et al. (2021). Computationally predicted SARS-COV-2 encoded microRNAs target NFKB, JAK/STAT and TGFB signaling pathways. Gene Rep. 22, 101012. doi: 10.1016/j.genrep.2020.101012
Bellae Papannarao, J., Schwenke, D. O., Manning, P., Manning, P., Katare, R. (2022). Upregulated miR-200c is associated with downregulation of the functional receptor for severe acute respiratory syndrome coronavirus 2 ACE2 in individuals with obesity. Int. J. Obes. 46, 238–241. doi: 10.1038/s41366-021-00984-2
Bennasser, Y., Le, S.-Y., Yeung, M. L., Jeang, K. T. (2004). HIV-1 encoded candidate micro-RNAs and their cellular targets. Retrovirology 1, 1–5. doi: 10.1186/1742-4690-1-43
Bernier, A., Sagan, S. M. (2018). The diverse roles of microRNAs at the host–virus interface. Viruses 10, 440. doi: 10.3390/v10080440
Boehm, E., Kronig, I., Neher, R. A., Isabella, E., Vetter, P., Kaiser, L., et al. (2021). Novel SARS-CoV-2 variants: The pandemics within the pandemic. Clin. Microbiol. Infection 27, 1109–1117. doi: 10.1016/j.cmi.2021.05.022
Bolisetty, M. T., Dy, G., Tam, W., Beemon, K. L. (2009). Reticuloendotheliosis virus strain T induces miR-155, which targets JARID2 and promotes cell survival. J. Virol. 83, 12009–12017. doi: 10.1128/JVI.01182-09
Brister, J. R., Ako-Adjei, D., Bao, Y., Blinkova, O. (2015). NCBI viral genomes resource. Nucleic Acids Res. 43, D571–D577. doi: 10.1093/nar/gku1207
Bruscella, P., Bottini, S., Baudesson, C., Pawlotsky, J. M., Feray, C., Trabucchi, M. (2017). Viruses and miRNAs: more friends than foes. Front. Microbiol. 8, 824. doi: 10.3389/fmicb.2017.00824
Chen, Y., Wang, X. (2020). miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res. 48, D127–D131. doi: 10.1093/nar/gkz757
Chen, L., Zhong, L. (2020). Genomics functional analysis and drug screening of SARS-CoV-2. Genes Dis. 7, 542–550. doi: 10.1016/j.gendis.2020.04.002
Duchaine, T. F., Fabian, M. R. (2019). Mechanistic insights into microRNA-mediated gene silencing. Cold Spring Harbor Perspect. Biol. 11, a032771. doi: 10.1101/cshperspect.a032771
Enguita, F. J., Leitão, A. L., McDonald, J. T., Zaksas, V., Das, S., Galeano, D., et al. (2022). The interplay between lncRNAs, RNA-binding proteins and viral genome during SARS-CoV-2 infection reveals strong connections with regulatory events involved in RNA metabolism and immune response. Theranostics 12, 3946. doi: 10.7150/thno.73268
Fani, M., Zandi, M., Ebrahimi, S., Soltani, S., Abbasi, S. (2021). The role of miRNAs in COVID-19 disease. Future Virol. 16, 301–306. doi: 10.2217/fvl-2020-0389
Golkar, Z., Battaria, R., Pace, D. G., Bagasra, O. (2016). Inhibition of Ebola virus by anti-Ebola miRNAs in silico. J. Infection Developing Countries 10, 626–634. doi: 10.3855/jidc.7127
Hassan, M., McLean, E., Bagasra, O. (2016). A computational analysis to construct a potential post-exposure therapy against pox epidemic using miRNAs in silico. J. Bioterror Biodef. 7:140. doi: 10.4172/2157-2526.1000140
Hosseini Rad, S. M. A., McLellan, A. D. (2020). Implications of SARS-CoV-2 mutations for genomic RNA structure and host microRNA targeting. Int. J. Mol. Sci. 21, 4807. doi: 10.3390/ijms21134807
Hussain, M., Asgari, S. (2014). MicroRNA-like viral small RNA from dengue virus 2 autoregulates its replication in mosquito cells. Proc. Natl. Acad. Sci. 111, 2746–2751. doi: 10.1073/pnas.1320123111
Hussain, M., Torres, S., Schnettler, E., Funk, A., Grundhoff, A., Pijlman, G. P, et al. (2012). West Nile Virus encodes a microRNA-like small RNA in the 3′ untranslated region which up-regulates GATA4 mRNA and facilitates virus replication in mosquito cells. Nucleic Acids Res. 40, 2210–2223. doi: 10.1093/nar/gkr848
Jafarinejad-Farsangi, S., Jazi, M. M., Rostamzadeh, F., Hadizadeh, M. (2020). High affinity of host human microRNAs to SARS-CoV-2 genome: An in silico analysis. Non-coding RNA Res. 5, 222–231. doi: 10.1016/j.ncrna.2020.11.005
Jiang, J., Lee, E. J., Schmittgen, T. D. (2006). Increased expression of microRNA-155 in Epstein-Barr virus transformed lymphoblastoid cell lines. Genes Chromosomes Cancer 45, 103–106. doi: 10.1002/gcc.20264
Kaur, T., Kapila, S., Kapila, R., Upadhyay, D., Kaur, M., Sharma, C., et al. (2021). Tmprss2 specific miRNAs as promising regulators for SARS-CoV-2 entry checkpoint. Virus Res. 294, 198275. doi: 10.1016/j.virusres.2020.198275
Kesheh, M. M., Mahmoudvand, S., Shokri, S. (2022). Long noncoding RNAs in respiratory viruses: A review. Rev. Med. Virol. 32, e2275. doi: 10.1002/rmv.2275
Khare, S., Gurry, C., Freitas, L., Schultz, M. B., Bach, G., Diallo, A., et al. (2021). GISAID’s role in pandemic response. China CDC Weekly 3, 1049. doi: 10.46234/ccdcw2021.255
Kim, W. R., Park, E. G., Kang, K.-W., Lee, S. M., Kim, B., Kim, H. S., et al. (2020). Expression analyses of microRNAs in hamster lung tissues infected by SARS-CoV-2. Molecules Cells 43, 953. doi: 10.14348/molcells.2020.0177
Klaver, B., Berkhout, B. (1994). Evolution of a disrupted TAR RNA hairpin structure in the HIV-1 virus. EMBO J. 13, 2650–2659. doi: 10.1002/j.1460-2075.1994.tb06555.x
Korber, B., Fischer, W. M., Gnanakaran, S., Yoon, H., Theiler, J., Abfalterer, W., et al. (2020). Tracking changes in SARS-CoV-2 spike: evidence that D614G increases infectivity of the COVID-19 virus. Cell 182, 812–827. doi: 10.1016/j.cell.2020.06.043
Kozomara, A., Birgaoanu, M., Griffiths-Jones, S. (2019). miRBase: From microRNA sequences to function. Nucleic Acids Res. 47, D155–D162. doi: 10.1093/nar/gky1141
Kucher, A., Koroleva, I. A., Zarubin, A., Nazarenko, M. S. (2022). MicroRNAs as the potential regulators of SARS-CoV-2 infection and modifiers of the COVID-19 clinical features. Mol. Biol. 56, 29–45. doi: 10.1134/S0026893322010034
Lei, X., Zhu, Y., Jones, T., Bai, Z., Huang, Y., Gao, S. J., et al. (2012). A kaposi’s sarcoma-associated herpesvirus microRNA and its variants target the transforming growth factor β pathway to promote cell survival. J. Virol. 86, 11698–11711. doi: 10.1128/JVI.06855-11
Li, L., Feng, H., Da, Q., Jiang, H., Chen, L., Xie, L., et al. (2016). Expression of HIV-encoded microRNA-TAR and its inhibitory effect on viral replication in human primary macrophages. Arch. Virol. 161, 1115–1123. doi: 10.1007/s00705-016-2755-5
Li, X., Fu, Z., Liang, H., Wang, Y., Qi, X., Ding, M., et al. (2018). H5N1 influenza virus-specific miRNA-like small RNA increases cytokine production and mouse mortality via targeting poly (rC)-binding protein 2. Cell Res. 28, 157–171. doi: 10.1038/cr.2018.3
Li, K., Shen, L., Zheng, P., Wang, Y., Wang, L., Meng, X., et al. (2022). Identification of MicroRNAs as potential biomarkers for detecting ischemic stroke. Genes Genomics 44, 9–17. doi: 10.1007/s13258-021-01060-9
Liu, S., Liu, X., Li, J., Zhou, H., Carr, M. J., Zhang, Z., et al. (2019). Long noncoding RNAs: novel regulators of virus-host interactions. Rev. Med. Virol. 29, e2046. doi: 10.1002/rmv.2046
Liu, W., Wang, X. (2019). Prediction of functional microRNA targets by integrative modeling of microRNA binding and target expression data. Genome Biol. 20, 1–10. doi: 10.1186/s13059-019-1629-z
Li, X., Zou, X. (2019). An overview of RNA virus-encoded microRNAs. ExRNA 1, 1–5. doi: 10.1186/s41544-019-0037-6
Marchi, R., Sugita, B., Centa, A., Fonseca, A. S., Bortoletto, S., Fiorentin, K., et al. (2021). The role of microRNAs in modulating SARS-CoV-2 infection in human cells: a systematic review. Infection Genet. Evol. 91, 104832. doi: 10.1016/j.meegid.2021.104832
Mishra, R., Kumar, A., Ingle, H., Kumar, H. (2020). The interplay between viral-derived miRNAs and host immunity during infection. Front. Immunol. 10, 3079. doi: 10.3389/fimmu.2019.03079
Morales, L., Oliveros, J. C., Enjuanes, L., Sola, I. (2022). Contribution of host miRNA-223-3p to SARS-CoV-Induced lung inflammatory pathology. Mbio 13, e03135-21. doi: 10.1128/mbio.03135-21
Morales, L., Oliveros, J. C., Fernandez-Delgado, R., Robert tenOever, B., Enjuanes, L., Sola, I., et al. (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
Mosharaf, M., Reza, M., Kibria, M., Ahmed, F.F., Kabir, M.H., Hasan, S., et al. (2022). Computational identification of host genomic biomarkers highlighting their functions, pathways and regulators that influence SARS-CoV-2 infections and drug repurposing. Sci. Rep. 12, 1–22. doi: 10.1038/s41598-022-08073-8
Natekar, J. P., Rothan, H. A., Arora, K., Strate, P.G., Kumar, M. (2019). Cellular microRNA-155 regulates virus-induced inflammatory response and protects against lethal West Nile virus infection. Viruses 12, 9. doi: 10.3390/v12010009
Panda, M., Kalita, E., Singh, S., Kumar, K., Roa, A., Prajapati, V. K., et al. (2022). MiRNA-SARS-CoV-2 dialogue and prospective anti-COVID-19 therapies. Life Sci. 305, 120761. doi: 10.1016/j.lfs.2022.120761
Pareek, S., Roy, S., Kumari, B., Jain, P., Banerjee, A., Vrati, S., et al. (2014). MiR-155 induction in microglial cells suppresses Japanese encephalitis virus replication and negatively modulates innate immune responses. J. Neuroinflamm. 11, 1–13. doi: 10.1186/1742-2094-11-97
Paul, S., Vázquez, L. A. B., Reyes-Pérez, P. R., Estrada-Meza, C., Alburquerque, R. A. A., Pathak, S., et al. (2022). The role of microRNAs in solving COVID-19 puzzle from infection to therapeutics: A mini-review. Virus Res. 308, 198631. doi: 10.1016/j.virusres.2021.198631
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
Prasad, A. N., Ronk, A. J., Widen, S. G., Wood, T.G., Basler, C.F., Bukreyev, A., et al. (2020). Ebola Virus produces discrete small noncoding RNAs independently of the host microRNA pathway which lack RNA interference activity in bat and human cells. J. Virol. 94, e01441–e01419. doi: 10.1128/JVI.01441-19
Samy, A., Maher, M. A., Abdelsalam, N. A., Badr, E. (2022). SARS-CoV-2 potential drugs, drug targets, and biomarkers: a viral-host interaction network-based analysis. Sci. Rep. 12, 1–14. doi: 10.1038/s41598-022-15898-w
Sevgin, O., Sevgin, K. (2021). Systematic review of microRNAs in the SARS-CoV-2 infection: Are microRNAs potential therapy for COVID-19. J. Genet. Genome Res. 8, 053. doi: 10.23937/2378-3648/1410053
Shirvani, H., Jafari, H., Moravveji, S. S., Faranghizadeh, F A., Talebi, M., Ghanavi, J., et al. (2022). Non-coding RNA in SARS-CoV-2: Progress toward therapeutic significance. Int. J. Biol. Macromolecules 222 (Pt A):1538–1550. doi: 10.1016/j.ijbiomac.2022.09.105
Sievers, F., Wilm, A., Dineen, D., Gibson, T.J., Karplus, K., Li, W., et al. (2011). Fast, scalable generation of high-quality protein multiple sequence alignments using clustal omega. Mol. Syst. Biol. 7, 539. doi: 10.1038/msb.2011.75
Song, Z., Xu, Y., Bao, L., Zhang, L., Pin, Y., Yajin, Q., et al. (2019). From SARS to MERS, thrusting coronaviruses into the spotlight. viruses 11, 59. doi: 10.3390/v11010059
Su, Y., Huang, Y., Wu, Y., Chen, H. F., Wu, Y. H., Hsu, C. C., et al. (2020). MicroRNA-155 inhibits dengue virus replication by inducing heme oxygenase-1-mediated antiviral interferon responses. FASEB J. 34, 7283–7294. doi: 10.1096/fj.201902878R
Swaminathan, G., Rossi, F., Sierra, L.-J., Gupta, A., Navas-Martin, S., Martin-Garcia, J., et al. (2012). A role for microRNA-155 modulation in the anti-HIV-1 effects of toll-like receptor 3 stimulation in macrophages. PLoS Pathog. 8 (9), e1002937. doi: 10.1371/journal.ppat.1002937
Talotta, R., Bahrami, S., Laska, M. J. (2022). Sequence complementarity between human noncoding RNAs and SARS-CoV-2 genes: What are the implications for human health? Biochim. Biophys. Acta (BBA)-Molecular Basis Dis. 1868, 166291. doi: 10.1016/j.bbadis.2021.166291
Tegally, H., Wilkinson, E., Lessells, R. J., Giandhari, J., Pillay, S., Msomi, N., et al. (2021). Sixteen novel lineages of SARS-CoV-2 in south Africa. Nat. Med. 27, 440–446. doi: 10.1038/s41591-021-01255-3
Tempel, S., Tahi, F. (2012). A fast ab-initio method for predicting miRNA precursors in genomes. Nucleic Acids Res. 40, e80–e80. doi: 10.1093/nar/gks146
Tran, V. D., Tempel, S., Zerath, B., Zehraoui, F., Tahi, F. (2015). miRBoost: boosting support vector machines for microRNA precursor classification. RNA 21, 775–785. doi: 10.1261/rna.043612.113
Coronaviridae Study Group of the International Committee on Taxonomy of Viruses(2020). The species Severe acute respiratory syndrome-related coronavirus: Classifying 2019-nCoV and naming it SARS-CoV-2. Nat. Microbiol. 5, 536–544. doi: 10.1038/s41564-020-0695-z
Thounaojam, M. C., Kundu, K., Kaushik, D. K., Swaroop, S., Mahadevan, A., Shankar, S.K., et al. (2014). MicroRNA 155 regulates Japanese encephalitis virus-induced inflammatory response by targeting src homology 2-containing inositol phosphatase 1. J. Virol. 88, 4798–4810. doi: 10.1128/JVI.02979-13
Vadivalagan, C., Shitut, A., Kamalakannan, S., Chen, R.M., Serrano-Aroca, A., Mishra, V., et al. (2022). Exosomal mediated signal transduction through artificial microRNA (amiRNA): A potential target for inhibition of SARS-CoV-2. Cell. Signalling 95, 110334. doi: 10.1016/j.cellsig.2022.110334
Van Laar, R., Lincoln, M., Van Laar, B. (2018). Development and validation of a plasma-based melanoma biomarker suitable for clinical use. Br. J. Cancer 118, 857–866. doi: 10.1038/bjc.2017.477
Vastrad, B., Vastrad, C., Tengli, A. (2020). Identification of potential mRNA panels for severe acute respiratory syndrome coronavirus 2 (COVID-19) diagnosis and treatment using microarray dataset and bioinformatics methods. 3 Biotech. 10, 1–65. doi: 10.1007/s13205-020-02406-y
Volz, E., Mishra, S., Chand, M., Johnson, R., Geidelberg, L., Hinsley, W.R., et al. (2021). Assessing transmissibility of SARS-CoV-2 lineage b. 1.1 7 England Nat. 593, 266–269. doi: 10.1038/s41586-021-03470-x
Wang, Z., Li, K., Wang, X., Huang, W., et al. (2019). MiR-155-5p modulates HSV-1 replication via the epigenetic regulation of SRSF2 gene expression. Epigenetics 14, 494–503. doi: 10.1080/15592294.2019.1600388
Yang, C. -Y., Chen, Y. -H., Liu, P. -J., Hu, W. -C., Tsai, K. W. (2022). The emerging role of miRNAs in the pathogenesis of COVID-19: Protective effects of nutraceutical polyphenolic compounds against SARS-CoV-2 infection. Int. J. Med. Sci. 19, 1340. doi: 10.7150/ijms.76168
Ying, H., Ebrahimi, M., Keivan, M., Khoshnam, S. E., Salahi, S., Farzaneh, M., et al. (2021). miRNAs; a novel strategy for the treatment of COVID-19. Cell Biol. Int. 45, 2045–2053. doi: 10.1002/cbin.11653
Yousefi, H., Poursheikhani, A., Bahmanpour, Z., Vatanmakanian, M., Taheri, M., Mashouri, L., et al. (2020). SARS-CoV infection crosstalk with human host cell noncoding-RNA machinery: An in-silico approach. Biomed. Pharmacother. 130, 110548. doi: 10.1016/j.biopha.2020.110548
Zeng, Q., Qi, X., Ma, J., Hu, F., Wang, X., Qin, H., et al. (2022). Distinct miRNAs associated with various clinical presentations of SARS-CoV-2 infection. Iscience 25, 104309. doi: 10.1016/j.isci.2022.104309
Zhang, S., Amahong, K., Sun, X., Lian, X., Liu, J., Sun, H., et al. (2021). The miRNA: a small but powerful RNA for COVID-19. Briefings Bioinf. 22, 1137–1149. doi: 10.1093/bib/bbab062
Zhang, Y., Fan, M., Geng, G., Liu, B., Huang, Z., Luo, H., et al. (2014). A novel HIV-1-encoded microRNA enhances its viral replication by targeting the TATA box region. Retrovirology 11, 1–15. doi: 10.1186/1742-4690-11-23
Keywords: coronavirus, SARS-CoV-2, miRNAs, non-coding RNAs, SARS
Citation: Periwal N, Bhardwaj U, Sarma S, Arora P and Sood V (2022) In silico analysis of SARS-CoV-2 genomes: Insights from SARS encoded non-coding RNAs. Front. Cell. Infect. Microbiol. 12:966870. doi: 10.3389/fcimb.2022.966870
Received: 11 June 2022; Accepted: 05 October 2022;
Published: 28 November 2022.
Edited by:
Rashed Noor, Independent University, BangladeshReviewed by:
Daniela Toro-Ascuy, Autonomous University of Chile, ChileThuy Linh Hoang, Institut Pasteur in Ho Chi Minh City, Vietnam
Copyright © 2022 Periwal, Bhardwaj, Sarma, Arora and Sood. 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: Vikas Sood, vikas1101@gmail.com