key: cord-333080-qytwbsne
authors: Alahari, Suresh K.; Yousefi, Hassan; Poursheikhani, Arash; bahmanpour, Zahra; Vatanmakanian, Mousa; Taheri, Mohammad; Mashouri, Ladan
title: SARS-CoV infection crosstalk with human host cell noncoding-RNA machinery: An in-silico approach
date: 2020-07-28
journal: Biomed Pharmacother
DOI: 10.1016/j.biopha.2020.110548
sha: 
doc_id: 333080
cord_uid: qytwbsne

Although 70% of the genome is transcribed to RNA in humans, only ∼2% of these transcripts are translated into proteins. The rest of the transcripts are defined as noncoding RNAs, including Long noncoding RNAs (LncRNAs) and MicroRNAs (miRNAs) that mostly function post-transcriptionally to regulate the gene expression. The outbreak of a novel coronavirus (SARS-CoV) has caused a major public health concern across the globe. The SARS-CoV is the seventh coronavirus that is known to cause human disease. There are currently no promising antiviral drugs with proven efficacy nor are there vaccines for its prevention. As of July 13, 2020, SARS-CoV has been infected more than 13 million cases in more than 213 countries, with an estimated mortality rate of ∼3%. Thus, it is of utmost important priority to develop novel therapies for COVID-19. It is not fully investigated whether noncoding RNAs regulate signaling pathways that SARS-CoV involved in. Hence, computational analysis of the noncoding RNA interactions and determining importance of key regulatory noncoding RNAs in antiviral defense mechanisms will likely be helpful in developing new drugs to attack SARS-CoV infection.. To elucidate this, we utilized bioinformatic approaches to find the interaction network of SARS-CoV/human proteins, miRNAs, and lncRNAs. We found TGF-beta signaling pathway as one of the potential interactive pathways. Furthermore, potential miRNAs/lncRNAs networks that the virus might engage during infection in human host cells have been shown. Altogether, TGF-beta signaling pathway as well as hub miRNAs, and LncRNAs involve during SARS-CoV pathogenesis can be considered as potential therapeutic targets.

Viral protein complexity is basically limited by their relatively small genomes; therefore, critical functions are mainly accomplished by cooperation between host and viral protein machinery.

This requirement makes the study of host-viral interactions is critical for the understanding of viral pathogenicity [1] . COVID-19 has infected over 13 million people and has killed more than 572,000 people by mid July 2020. This has driven the world into hysteria and panic, causing social and economic disruption. There are currently no promising antiviral drugs or vaccines to combat COVID-19. This is largely due to our limited knowledge about the molecular mechanisms of the virus pathogenesis [2] .

MicroRNAs (miRNAs) are short non-coding RNAs that inhibit the expression of their target genes by directly binding to target mRNAs [3] . Long noncoding RNAs (lncRNAs) are a category of cellular RNAs that are longer than 200 nucleotides in length. Accumulating evidence suggests that miRNAs and lncRNAs play key roles in different disease pathogenesis [4, 5] . More specifically, miRNA have common binding sites for both mRNAs and lncRNAs; thus, lncRNAs can act as competitive endogenous RNA (ceRNA) to suppress the miRNA function or act as precursors and encode miRNAs [4] . These miRNA regulatory functions have been estimated to regulate around 60% of mammalian genes [6] . LncRNAs/ miRNAs play an important regulatory role in the battle between virus and host, especially in regulating the transcription of virus and host genes, stability and translation of mRNAs, and host antiviral response [6] . Developing functional computational models and networks to predict potential SARS-CoV -miRNA/lncRNA association may benefit not only the understanding of COVID-19 mechanism at the noncoding RNA level, but also the detection of disease biomarkers for disease diagnosis, treatment, prognosis, and prevention. Surprisingly, known experimentally confirmed disease-lncRNA associations are still very limited. Since multiple ways of interaction between miRNAs, lncRNAs, and mRNA have been reported to play key roles in determining the cellular functions during viral infection, it is essential to discover these interactions in an integrated fashion to comprehensively decipher the networks and key regulatory noncoding-RNA hubs underpinning the pathology of SARS-CoV.

Although we learned a great deal of coronavirus infection, the effective approaches for treatment are still scarce. Unsurprisingly, there is little knowledge of the detailed molecular interactions for J o u r n a l P r e -p r o o f the SARS-CoV infection crosstalk with noncoding RNAs network in the cell. Several signaling pathways regulate miRNAs and lncRNAs that affect the viral function and the infection progress, including replication, translation, and even modulating the host expression. Since it is also important to consider the virus-host protein-protein interaction in order to understand its pathogenesis during infection [7, 8] ; we used Viruses-STRING, a protein-protein interaction database specifically designed to analyze the virus-host interactions [9] . Subsequently, the Gene Set Enrichment Analysis (GSEA) was used to assess and interpret the coordinate pathways that might be potentially involved in these interactions [10] . Furthermore, we used the miRWalk2.0 database to generate a network for a comprehensive overview of miRNA-lncRNA interactions using experimentally verified data [11] . Our findings revealed a network of miRNAs/lncRNAs which regulate these pathways. Based on this, we designed a network-based model to propose the most important regulatory factors and hubs with a functional role in regulating the viral infection process using a noncoding RNA perspective.

To generate an interactome network of Human SARS coronavirus (SARS-CoV) with human proteins, we used the STRING database http://viruses.string-db.org/ (version 10.5) [9] . Next, the Gene Set Enrichment Analysis (GSEA) (https://www.gsea-msigdb.org/) was performed for the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotations [10, 12] . Then, the most important canonical pathways with the highest number of genes were selected. In order to disclose the most possible miRNAs regulating the selected mRNAs, the MiRTarBase (http://miRTarBase.cuhk.edu.cn/), one of the most important experimentally validated miRNA-target interactions (MTIs) databases, was employed. For the sake of accuracy, only functional MTIs were chosen for further study [12] . Further, miRNA-lncRNA interactions were determined using the miRWalk2.0 database (http://zmf.umm.uniheidelberg.de/apps/zmf/mirwalk2/), a database on predicted and published noncoding RNAs; and seed length (SL) ≥ 13 were selected as the most significant interactions [11] . Finally, the J o u r n a l P r e -p r o o f lncRNA /miRNA/mRNA network was virtualized by Cytoscape 3.8.0. Hub targets were analyzed by cytoHubba package, a Cytoscape plug-in for finding hub objects in a network, and 5 top hub nodes with higher degrees were recognized.

Human SARS coronavirus (SARS-CoV) and the Homo sapiens interaction network was downloaded as a vector image from the viruses STRING dataset. Protein-protein interactions were obtained and shaped a network with median confidence (scores greater than 0.4) between SARS-CoV and human host cells. The obtained protein-protein interactions are summarized in 

As many as 100 terms were examined in the canonical pathway analysis. It was found that many of these pathways are central pathways, meaning that most of the genes involved in this network are associated with canonical pathways. 20 top canonical pathway terms such as TGF-beta signaling pathway (P value=6.99E-21; gene count= 11), Wnt signaling pathway (P value=4.09E-14; gene count=9), MAPK signaling pathway (P value=7.51E-12; gene count=9) are shown in Figure 2 and the genes involved in these pathways are shown in Table 1 . Among these 20 pathways, the TGF-beta signaling pathway with the highest number of genes was selected for further study. In addition, based on gene ontology, 100 terms were selected such as those associated with response to transforming growth factor-beta (P value=3.22E-23; gene count=15), transforming growth factor-beta receptor cytoplasmic mediator activity (P value=2.95E-23; gene count=8), cellular response to endogenous stimulus (P value=8.03E-23; gene count=23), transforming growth factor-beta receptor signaling pathway (P value=7.79E-23; gene count=14).

The 20 top enriched GO functions are shown in Figure 3 . Table 2 .

Over the last few decades, an overgrowing line of evidence indicated a link of noncoding RNA molecules with the gene regulatory networks including those controlling the disease pathogenesis [6] . Having the tight interaction between viral proteins with their host proteome, it is expected that the viral infection could influence the noncoding RNA network in humans.

There is a current urgent need for studying the pathogenesis of human SARS-CoV which has hit the world in a tremendously ascending trend. Thus, a huge number of research groups are trying very hard to find rapid diagnostic methods, studying current treatment options, developing novel therapeutic options, and hopefully developing an effective vaccine. However, fine-tuning crosstalk between viral proteins and host noncoding RNA machinery is almost an underappreciated layer in this investigation race. On the other hand, studying the noncoding RNAs in response to SARS-CoV infection using wet-lab strategies or low-scale methods might not fit the current emergency call. Therefore, we aimed to follow a comprehensive J o u r n a l P r e -p r o o f bioinformatics approach to dissect the strongest possible changes in the noncoding RNA network in response to the viral infection in humans using the available online databases. Our in-silico analysis has built a network of protein-protein interaction between the Human SARS coronavirus (SARS-CoV) and host proteome (Figure 1) , as well as strong miRNA-mRNA-lncRNA crosstalk ( Figure 4 and Table 2 ) possibly modulating the human response to the viral infection.

In this study, we have performed multiple layers of bioinformatics analyses in order to obtain a big picture of the most significant events that occurred following the virus-host interaction in terms of changes in cell signaling pathways and gene expression via noncoding RNAs. We applied the viruses. STRING online database as a powerful platform to explore the most significant protein-protein interactions between SARS-CoV viral particles and host proteome.

Considering the median confidence scores greater than 0.4, our analysis predicted that 15 viral proteins might be in direct interactions with 40 host proteins, as depicted in the network ( Figure   1 ). In order to gain a functional prediction of the interacting proteins, the list of obtained strong interactions was then reanalyzed using gene ontology and gene set enrichment analysis (GSEA).

A more concentrated list of the major pathways involved were provided using these in silico techniques, among which the interacting genes showed a great association with the TGF-β signaling and its receptor cytoplasmic activity. Also, significant GO categories for overlapped genes using the threshold of P < 0.05 and FDR < 0.05 showed meaningful predicted changes in enriched GO functions of the genes involved in interaction network, including those controlling the transforming growth factor-beta, transforming growth factor-beta receptor cytoplasmic mediator activity, and cellular response to endogenous stimulus. Based on our data, it is assumed that SARS-CoV infection in mammals raises a high level of inflammatory and immune responses by direct viral-host protein interactions. There are also studies that have been shown the role of TGF-beta signaling pathway in SARS-CoV infection [13] [14] [15] .

In an effort to decipher the probable non-coding RNAs molecules involved, we further analyzed the most significant abovementioned gene sets using the miRTarBase and miRWalk algorithm.

This analysis resulted in a network of miRNA-mRNA and miRNA-lncRNA interaction illustrated in Figure 4 . Various mechanisms based on miRNA-mediated regulation of viruses have been reported [16] . This highly coordinated interaction may eventually result in positive outcomes with regard to the host immune response or in favor of viral proliferation [17, 18] . This J o u r n a l P r e -p r o o f is highly depends on the miRNA target [19] . More precisely, the small noncoding RNA might act directly by guiding an Argonaute protein on the viral mRNA or might target a cellular transcript that is important in the host immune response against viral infection. Consistently, our analysis confirmed that a change in miRNA network targeting the main immune hub signaling pathways responding to SARS-CoV infection potentially act as important immune-modulatory mechanisms in mammals. These findings predict a highly immunogenic viral-host interaction which may be very promising with regard to vaccine development. J o u r n a l P r e -p r o o f

Experimental Analysis of Viral-Host Interactions

The COVID-19 pandemic in the US: a clinical update

Regulation of microRNA biogenesis and its crosstalk with other cellular pathways

Long noncoding RNAs and exosomal lncRNAs: classification, and mechanisms in breast cancer metastasis and drug resistance

Non-coding RNAs in human disease

Functional Interplay between RNA Viruses and Non-Coding RNA in Mammals

Viral perturbations of host networks reflect disease etiology

Exploring the SARS-CoV-2 virus-host-drug interactome for drug repurposing

STRING: A virus-host protein-protein interaction database. Viruses

Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles

miRWalk2. 0: a comprehensive atlas of microRNA-target interactions

PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes

SARS coronavirus papain-like protease induces Egr-1-dependent up-regulation of TGF-β1 via ROS/p38 MAPK/STAT3 pathway

A potential treatment of COVID-19 with TGF-β blockade

Cellular microRNA networks regulate host dependency of hepatitis C virus infection

Modulation of hepatitis C virus RNA abundance by a liver-specific MicroRNA

microRNA-122 stimulates translation of hepatitis C virus RNA. Embo j

On the Importance of Host MicroRNAs During Viral Infection

The authors declare that there is no conflict of interest.J o u r n a l P r e -p r o o f