Abstract
Viral proteins encode numerous antiviral activities to modify the host immunity. In this article, we hypothesize that viral genomes and gene transcripts interfere with host gene expression using passive mechanisms to deregulate host microRNA (miRNA) activity. We postulate that various RNA viruses mimic or block binding between a host miRNA and its target transcript, a phenomenon mediated by the miRNA seed site at the 5′ end of miRNA. Virus-encoded miRNA seed sponges (vSSs) can potentially bind to host miRNA seed sites and prevent interaction with their native targets thereby relieving native miRNA suppression. In contrast, virus-encoded miRNA seed mimics (vSMs) may mediate considerable downregulation of host miRNA activity. We analyzed genomes from diverse RNA viruses for vSS and vSM signatures and found an abundance of these motifs indicating that this may be a mechanism of deceiving host immunity. Employing respiratory syncytial virus and measles virus as models, we reveal that regions surrounding vSS or vSM motifs have features characteristics of pre-miRNA templates and show that RSV viral transcripts are processed into small RNAs that may behave as vSS or vSM effectors. These data suggest that complex molecular interactions likely occur at the host-virus interface. Identifying the mechanisms in the network of interactions between the host and viral transcripts can help uncover ways to improve vaccine efficacy, therapeutics, and potentially mitigate the adverse events that may be associated with some vaccines.
Introduction
Regulation of gene expression is a complex occurrence involving transcriptional and post-transcriptional mechanisms (Burgess, 2017). Translation of host mRNAs is regulated by small evolutionarily conserved small non-coding RNAs or miRNAs (Chen et al., 2017; Morales et al., 2017; Murashov, 2017) principally by sequence complementarity. Sequence complementarity is also fundamental to pathways such as clustered regularly interspaced short palindromic repeats (CRISPR) and RNA interference (RNAi) pathways (Panek et al., 2016; Kaikkonen and Adelman, 2018; Damas et al., 2019; Hussain et al., 2019) of which miRNAs have an integral function. Host miRNA genes can be intronic, intergenic, or independent transcription units typically processed by RNA polymerase II. Initial transcription produces primary miRNA transcripts (pri-miRNA) which are processed by a nuclear microprocessor complex to produce pre-miRNAs (Bartel, 2004; Finnegan and Pasquinelli, 2013; Ha and Kim, 2014). Pre-miRNAs are exported out into the cytoplasm and further processed into a 18–25nt long dsRNA by class III RNAse enzyme Dicer to form the mature mRNA duplex. One strand of the duplex is thermodynamically stable and referred to as the guide strand while the other strand is the passenger strand. The mature guide strand at the 5' end encodes a 6 nt ‘miRNA seed’ that is complementary to a miRNA recognition element (MRE) in target transcripts (Bartel, 2004; Finnegan and Pasquinelli, 2013; Ha and Kim, 2014). Sequence dependent pairing of the seed site with MRE can lead to either mRNA decay or blocked translation (Fabian et al., 2010). Each miRNA can target several genes although the stoichiometry of this interaction is inadequately defined (Weill et al., 2015). Non-seed-mediated miRNA regulation has also been recognized (Fabian et al., 2010; Ghosal et al., 2014; Li et al., 2014; Kumari et al., 2016), is an alternate mechanism of miRNA function (Cloonan, 2015), and multiple miRNAs may regulate a single gene (Bartel, 2009; Friedman et al., 2009; Fabian et al., 2010). There are several features that may contribute to function as different motifs may be able to interact among viral and host transcripts given each position can have one of the four nucleotides (A/U/G/C) and the six nt in the seed site can be promiscuous and bind to many targets (Friedman et al., 2009). Host gene miRNA regulation is well established and has been demonstrated to modulate ~60% of the human transcriptome (Friedman et al., 2009; Dong and Lou, 2012; Hashimoto et al., 2013; Jia et al., 2014). It is well known that miRNAs regulate cell physiology during normal homeostasis as well as during disease states (Karnati et al., 2015; Alipoor et al., 2016; Vishnoi and Rani, 2017; Olejniczak et al., 2018). Importantly, differential expression of miRNAs has been used as biomarkers for diagnosis, treatment and prognosis, and miRNA expression is modified in response to viral infection (Sullivan and Ganem, 2005; Piedade and Azevedo-Pereira, 2016; Auvinen, 2017; Trobaugh and Klimstra, 2017).
DNA viruses typically replicate with a high fidelity and encode their own miRNAs (Klinke et al., 2014; Flor and Blom, 2016; Albanese et al., 2017; Qin et al., 2017) to regulate virus and host gene expression. In contrast, RNA viruses lack replicative fidelity and arise in the host following infection as swarms of quasispecies (Steinhauer et al., 1992). The quasispecies generally have poor replicative fitness. RNA virus replication and gene transcription are catalyzed by RNA-dependent RNA polymerase (RdRP) which is error-prone (Perez-Rodriguez et al., 2016). It is not understood what proportion of quasispecies are translated during infection, though non-canonical start codons can be used during RNA viral protein synthesis (Firth and Brierley, 2012). Thus, quasispecies can contribute to antigenic diversity via non-canonical translation of alternative viral proteins, and inhibition of these processes can improve viral yield during vaccine production via molecular breeding, increase vaccine safety and stability and reduce potential adverse events (Perez-Rodriguez et al., 2016). While a lack of exonuclease proof reading activity occurs in the RdRP, a feature explaining how low fidelity arises, it remains unclear why RNA viruses have and maintain low fidelity RdRP. It is possible that low RdRP fidelity facilitates virus replication (Hopfield, 1974; Back et al., 1996) by allowing RNA viruses to escape situations where unfit mutations predominate fitness leading to species collapse and attenuation (Lauring and Andino, 2010).
In this article, we speculate that RNA viruses interfere with the host RNAi machinery that regulates both foreign and endogenous gene expression by miRNAs. We postulate that RNA viral genomes and gene transcripts encode motifs that can either mimic or block native miRNA activity through sequence homology or complementarity. We analyzed the genomes of several RNA virus families (Supplementary Tables 1, 2) and identified vSMs that potentially may mimic miRNA seed activity and thus increase native miRNA repression of host antiviral pathways. We also identified several viral seed sponges that can potentially block host miRNAs and relieve native miRNA suppression of pro-viral host genes. We focused our analysis on the Paramyxoviridae family of RNA viruses owing to their impact on human and animal health. Paramyxoviruses have negative sense, non-segmented, single-stranded RNA genomes that are transcribed in a gradient leading to a differential abundance of viral transcripts with all steps in the viral life cycle occurring in the cytosol where host miRNAs also regulate gene expression. Paramyxoviruses are classified into two subfamilies Paramyxovirinae (e.g., Avulavirus, Henipavirus, Morbilivirus, Respirovirus, and Rubulavirus genera) and Pneumovirinae (e.g., Pneumovirus and Metapneumovirus genera) (Aguilar and Lee, 2011; Amarasinghe et al., 2017; Rima et al., 2017). Paramyxovirinae members causing morbidity and mortality include measles (MV), Mumps (MuV), Hendravirus (HV), Nipah virus (NiV), and the Pneumoviruses, i.e., respiratory syncytial virus (RSV) and human metapneumovirus (hMPV).
In this article, we propose that quasispecies enable RNA viruses to modulate host gene expression by regulating miRNA function via sequence complementarity or identity with the miRNA seed sites. We also suppose that vSM function to increase native miRNA-based suppression, while vSS inhibit native miRNA activity and increase host gene expression to the advantage of the virus. Preliminary analysis has identified a number of vSM or vSS in several Paramyxovirus genomes. For example, for RSV the regions that neighbor potential vSS or vSM are predicted to form stable stem loop structures that are typically substrates for nuclear and cytosolic RNAses of the RNAi pathway (Cai et al., 2004; Ritchie et al., 2007; Shu et al., 2007; Kurihara and Watanabe, 2010). These findings suggest that these regions in the viral genome can be templates for cytosolic Dicer activity. We have confirmed that during RSV infection gene transcripts are processed into sncRNAs and have identified viral transcripts that harbor vSS or vSM using next generation sequencing (NGS). Genomic analyses show that these motifs are more abundant in genes that have known or predicted immunomodulatory function or are involved in viral replication/transcription. These data suggest that interactions with host miRNAs may be part of a mechanism to modulate miRNA-mediated host regulatory pathways and regulate viral gene expression. These findings have important implications for better understanding of host-virus interaction as well as rational vaccine design strategies.
Materials and Methods
Viruses and Cell Culture
Mycoplasma-free virus stocks of wild type RSV strain A2 were expanded in Vero cells (ATCC CCL81) and maintained in DMEM (Hyclone, Salt Lake City, Utah USA) supplemented with 5% heat-inactivated fetal bovine serum (Hyclone, USA) as previously described (Oshansky et al., 2009). A549 cells (ATCC-CCL185) grown in DMEM supplemented with 5% serum as above were used for all infections. A549 cells were infected at a multiplicity of infection (MOI) of 1.0 as previously described (Oshansky et al., 2009). SHSY5Y cells were maintained in DMEM with 10% heat inactivated FBS.
Nucleotide Sequence Analysis
Complete genome sequences for RSV, MV, HMPV, MuPV, NDV, HV, NiV, MuV were from the National Center for Biotechnology Information (NCBI). Accession numbers for all sequences analyzed are given in Supplementary Table 1. A local database of human mature miRNA sequences version 21.0 was constructed locally in BioEdit. Viral sequences were analyzed using BLASTN (Altschul et al., 1990; Altschul and Gish, 1996; Altschul and Pop, 2017) against this local database using the parameters Expect value (E) = 10, matrix (M) = BLOSUM62, Low complexity Repeat masking = OFF and output = Tabular. CSV files were imported into Microsoft Excel 2010 and filtered to identify hits where miRNA start or miRNA end was ≤3. Hits in the same orientation as the miRNA (5′-3′) were designated vSMs while vSSs were in anti-sense orientation.
Structure Prediction of RSV vSMs (mimics) and vSSs (sponges)
Nucleotide sequences (100 nt) flanking each predicted vSM or vSS for RSV were extracted from parental genomic sequence and analyzed by miRNAfold (Sullivan and Ganem, 2005) or RNA structure (Xu and Mathews, 2016). Structures were visualized with VARNA GUI (Hashimoto et al., 2013). Pre-miRNA sequences were used as controls in prediction. Hybridization stability was calculated using RNA hybrid (Jia et al., 2014).
Analysis of Small RNA Processing
Total RNA from mock-treated or RSV A2-infected (MOI = 1.0) Vero cells was isolated and fractionated using RNAzol RT according to the manufacturer protocol (MRCgene, Ohio). Fractionated small RNA was polyadenylated and then reverse transcribed with Protoscript II (New England Biolabs, MA) RSV G- and L-specific reverse oligomers followed by PCR using gene specific forward and reverse oligomers in a reaction containing 10 uM final primer concentration as per manufacturer's recommendations. PCR amplicons were resolved on a denaturing 12% PAGE gel in 1X TBE buffer and stained with SYBR Gold. Sequences for RSV G and L oligomers are provided in Supplemental Methods.
Next Generation Sequencing (NGS)
Type II respiratory epithelial (A549) cells or neuronal (SHSY5Y) cells were mock-treated or infected with RSV A2 (MOI = 1.0) for 24 h. Total RNA was isolated using RNAzol RT (MRCgene, Ohio) per the manufacturer protocol and size fractionated. Small RNA was quantified using Qubit broad range RNA kit (Invitrogen, USA). Size fractionated small RNA was ligated to proprietary Illumina adaptors using T4 RNA Ligase deletion mutant (Epicenter, USA). Adaptor ligated RNA was reverse transcribed and amplified for limited number of cycles during which index barcodes were attached to each cDNA pool. cDNA libraries for each sample at the end of the incubation temperature were again analyzed on a Tapestation 2200, quantified and then pooled at equimolar ratios as recommended by manufacturer. Pooled cDNAs were denatured using freshly prepared 0.2 N NaOH at room temperature for 3 min and then mixed with hybridization buffer. Diluted libraries at 20 pM concentration were loaded onto Illumina MiSeq and run to generate fastq files as per Illumina MiSeq protocol. Sequencing reads obtained had quality scores >30, were trimmed of adaptors, and were then analyzed by BLAST against RSV A2 genome to determine virus specific transcripts and determine gene location for these transcripts.
Results
Paramyxovirus Transcripts Mimic or Bind to Host miRNAs by Seed Sequence
Paramyxovirus genome replication, gene transcription, and translation occur in the cytoplasm, which is the primary site of host post-transcriptional gene regulation by miRNAs. To identify the sequence motifs in Paramyxovirus genomes that could mimic or inhibit host miRNA function based on seed sequence homology/complementarity, RSV, HMPV, MuPV, NDV, HV, NiV, MV, and MuV genomes were analyzed by basic local alignment search tool (BLAST) (Altschul et al., 1990) against a database of either human mature miRNAs or cognate species using criteria optimized for finding short matches (matrix = BLOSUM 62, Expect value = 10). Hits were filtered to identify matches in the miRNA seed site that were either in the sense (+/+) or antisense (±) orientation relative to the miRNA. Sense orientation hits were designated as vSMs (Table 2), while antisense hits were designated vSSs (Table 3). Across genera, we identified 85 vSSs compared to 60 vSMs for the type strains (Figure 1A). The ratio of the number of vSSs to vSMs was typically in the 1.5–2.0 + range for most genera except HMPV where they were in approximately equal numbers (Figure 1B). The ratio of vSS/vSM was determined to determine if vSMs vs. vSSs were selected for or against during evolution. As noted in the Introduction, vSMs are expected to enhance native miRNA mediated suppression while vSSs are expected to bind to miRNAs and relieve their repressor effect. A high vSM/vSS ratio would indicate that mimicking host miRNA activity could be a predominant mechanisms of modulating host miRNA activity. Alternatively, a low vSM/vSS ratio might indicate that inhibiting miRNA activity would be a more conserved mode of regulating miRNA function. Most vSMs or vSSs were located in the L gene across the Paramyxovirus genera followed by F protein and G protein genes. Details of vSMs and vSSs in are discussed below.
Figure 1
Respiratory Syncytial Virus
RSV is grouped as A and B strains based on the diversity in the G protein (Johnson and Graham, 2004; Papenburg and Boivin, 2010; Pangesti et al., 2018). For strain A2 (GenBank accession number M74568), we did not identify any vSM or vSS while for strain B1 (GenBank accession number AF013254), one vSM was identified in the L gene. Analysis of clade A (326 sequences identified) and clade B (60 sequence strains) (accession numbers in Supplementary Table 1) identified two vSMs, one mapping to miR-4719 in the G gene and the other mapping to miR-556-3p in L gene (Figure 2A). Strains that encoded vSM-4719 did not encode vSM-556-3p and vice-versa. For B1, the single vSM corresponding to miR-4311 mapped to the L gene. The viral gene segments that encoded vSSs included miR-2278 in the NS1/1c gene and miR-4280 and miR-1273f in the G gene (Figure 3A).
Figure 2
Figure 3
Human Metapneumovirus
HMPV, like RSV, is grouped as A and B strains based on G gene diversity (Kahn, 2006; Papenburg and Boivin, 2010; Schuster and Williams, 2014). A single vSM for miRNA let-7f was located in the L gene (Figure 2B), and a single vSS complementary to multiple members of the miR-548 family was found in the M2 gene (Figure 3B).
Hendravirus
Analysis of the HV genome identified three major vSMs for miR-4661-3p, miR-30e-5p, and miR-580-3p in F, G, and P genes (Figure 2C). vSSs for miR-6815-3p and miR-29c-3p in the N gene, miR-374a-5p, and miR-26a-5p in the L gene, and miR-1250-5p in the G gene were identified (Figure 3C).
Murine Pneumovirus
Analysis of MuPV genome identified several conserved vSMs containing seed sequence of miR-4311 in the N gene, miR-539-5p in the SH gene, and miR-2681-3p in the G gene (Figure 2D), as well as vSSs complementary to miR-6825 and miR-138-5p in the P and G genes, respectively (Figure 3D).
New Castle Disease Virus
Avian Paramyxovirus-1 (APMV-1) is the etiological agent of NDV. Analysis of the APMV-1 genome identified four highly conserved vSMs for miR-1681, miR-6703-3p, and miR-6648-3p in the L gene and miR-146b-3p in the M gene (Figure 2E). The genome also encoded two vSSs each for miR-1465, miR-7456-5p in the F gene, and miR-6516-3p and miR-6662-3p in the L gene, and miR-1636 in the intergenic region (Figure 3E). APMC-1 encoded the highest number of vSMs and vSSs among all paramyxoviruses.
Nipah Virus
Analysis of NiV genomes identified two vSMs for miR-3925-5p and miR-4698 in the L gene (Figure 2F). The genomes also encoded vSSs against miR6846-3p in the G gene, miR-6817-5p in the P gene and miR-3173-5p in the F gene (Figure 3F).
Measles Virus
The MV genome had a vSM for miR-569 in the L gene (Figure 2G), and a vSS for miR-6770-3p in the P gene (Figure 3G).
Mumps Virus
MuV genome has a vSM in the L gene for mir-548i (Figure 2H). The genome also encodes five vSSs against miR-4781-3p in the NP gene, miR-6877-3p in the P gene, miR-376c-3p in the M gene, miR-143-5p in the SH gene, and miR-6861-5p in the HN gene (Figure 3H).
These comparative analyses of Paramyxoviruses revealed an enrichment of vSSs relative to vSMs and localization of vSM sequences in L or G genes. It is not clear if these miRNA sequences evolved as an adaptation to a host response or are an outcome of low polymerase (L gene) fidelity. The conservation of these sequences in the circulating viral strains of multiple paramyxoviruses suggests that they may facilitate virus replication possible by contributing to immune evasion.
Viral RNA Is Processed to Smaller Transcripts During Infection
Stem-loop structures are a feature of miRNA templates and are recognized by both nuclear and cytosolic class III RNAses (Fabian et al., 2010; Weill et al., 2015). Several studies suggested RSV genome/antigenome might not fold to form secondary structures due to N protein encapsidation (Ghosal et al., 2014; Kumari et al., 2016). However, recent data has shown that RNA probes can bind to RSV genomic RNA (Bartel, 2009; Li et al., 2014; Cloonan, 2015) support the hypothesis that the RSV genome/antigenome is exposed during infection and may be able to form secondary structural motifs that can act as templates for Dicer, a class III cytosolic RNAse (Friedman et al., 2009). There is no published data to indicate that viral transcripts are protein encapsidated in the cytosol during infection. We analyzed genomic sequences near predicted vSMs and vSSs, and RNA secondary structure predictions indicated the potential of these regions to fold into stable stem loop structures similar to pre-miRNA transcripts (Figures 4A,B). Randomization of the sequence disrupted the structure. Stability of interaction between miRNAs and vSSs were computationally determined with RNAhybrid, a tool for finding the minimum free energy hybridization of a long and a short RNA to predict miRNA targets (Jia et al., 2014). Interactions between the predicted vSS and complementary miRNA are very stable as evident from the low mean free energy of the hybrid (Table 1). To establish if the RSV transcripts are processed into small RNAs during infection, Vero cells were infected with RSV A2 (MOI = 5.0) for 24 or 48 h, and size-fractionated small RNA was reverse transcribed using G- or L- specific oligonucleotides and amplified using G- or L-specific forward and reverse oligomers. MOI of 1.0 was used to infect A549 or SHY5Y cells with RSV or MV followed by small RNA isolation at 24hpi and NGS sequencing. MOI of 5.0 was used to infect Vero cells for RT-PCR analysis of small RNA products. The different MOIs were used because of the different sensitivities of the assay methods. With NGS, we expected to observe even small differences owing to the extremely high sensitivity of the method. In contrast, for Vero cells, RSV infection produces little to no cytopathic effect at 24–48 hpi, thus necessitating higher MOI infection. Electropherograms clearly showed bands of ~40, 50, 70, and 80 bp for the G gene, and 40 and 80 bp for the L gene along with several bands in the ~120 bp size range that are similar to host pre-miRNA sizes (Figure 4C). These data support the hypothesis that vSM- and vSS-encoding regions of the genome are folded and processed into smaller transcripts potentially via cellular RNAi machinery. Experimental validation of predicted structures is out of the scope of this hypothesis article.
Figure 4
Table 1
| Virus | vSS | miRNA | MFE (Kcal/mol) |
|---|---|---|---|
| RSV A2 | NS1-2278 | miR-2278 | −32.4 |
| G-4280 | miR-4280 | −25.2 | |
| G-1273f | miR-1273f | −29.7 | |
| RSV B1 | NS1-2278 | miR-2278 | −31.9 |
| M/IGR4-7161-3p | miR-7161-3p | −34.0 | |
| G-2113 | miR-2113 | −29.0 |
Binding energy calculations for vSSs per RNA duplex.
Stability of vSS miRNA was determined using RNAhybrid. Mean free energy (MFE) denotes stability of interaction. NS1, Non-structural gene-1; G, RSV G gene; M, Matrix gene. Number indicate the potential miRNA target.
Table 2
| Genera | vSMs |
|---|---|
| RSV A | miR-4719 > miR-556-3p > miR-4801 > miR-8074 > miR-3613-3p > miR-1253 > miR-3618 > miR-182-5p > miR-8485 > miR-4770 |
| RSV B | miR-4311 > miR-6780a-5p > let-7b-3p > miR-499a-3p |
| HMPV A | miR-4662a-5p>miR-20b-5p> let-7f-5p > miR-5194 > miR-548u > miR-1267 > miR-3925-5p > miR-3148 > miR-4799-5p > miR-6884-5p > miR-335-5p |
| HMPV B | miR-6884-5p > miR-592 > miR-1267 > miR-20b-5p |
| NDV | miR-6876-5p > miR-1306-5p > miR-1264 > miR-4280 > miR-6888-3p > miR-6855-3p > miR-338-3p > miR-548aq-5p > miR-6875-5p > miR-8069 > miR-3925-3p > miR-1249-5p > miR-892b > miR-4656 > miR-4679 |
| HV | miR-4661-3p > miR-4651 > miR-30e-5p > miR-6734-5p |
| NiV | miR-3925-5p> miR-4698 > miR-4693-3p |
| MV | miR-374a-3p > miR-5571-5p > miR-3942-3p > miR-4662a-5p> miR-5688 > miR-6861-3p > miR-4275 > miR-3617-3p > miR-4273 > miR-569 > miR-3202 > miR-1297 > miR-26b-5p > miR-3934-3p > miR-1255b-5p > miR-6755-3p > miR-1469 > miR-324-3p > miR-548ah-5p > miR-593-3p > miR-4518 > miR-6803-3p > miR-4690-5p > miR-4466 > miR-210-3p > miR-2909 > miR-3909 > miR-4505 > let-7b-3p > miR-2276-5p > miR-4420 > miR-4498 > miR-6075 > miR-2276-3p > miR-4486 > miR-6793-3p |
| MuV | miR-6861-5p > miR-6877-3p > miR-143-5p > miR-4433a-5p > miR-1236-3p > miR-4781-3p > miR-219a-2-3p > miR-1233-5p > miR-4645-5p > miR-3131 > miR-6881-5p > miR-4269 > miR-7152-3p > miR-4286 > miR-4323 |
vSMs in Paramyxovirus genomes.
vSMs in bold are conserved across two or more viruses. vSM hits are listed in order of their decreasing abundance among viral genomes. RSV, respiratory syncytial virus; hMPV, human metapneumovirus; NDV, New Castle disease virus; HV, Hendravirus; NiV, Nipah virus; MV, Measles virus; MuV, Mumps virus.
Table 3
| Genera | vSSs |
|---|---|
| RSV A | miR-2278 > miR-4280 > miR-1273f > miR-4753-3p > miR-592 > miR-4639-3p |
| RSV B | miR-2278>miR-7161-3p > miR-2113 > miR-7151-5p > miR-8068 > miR-369-3p |
| HMPV A | miR-8063>miR-323b-5p>miR-6831-5p>miR-337-5p> miR-5584-3p > miR-494-5p>miR-3678-3p |
| HMPV B | miR-891a-3p>miR-6831-5p>miR-323b-5p>miR-337-5p>miR-8063>miR-3678-3p > miR-6730-3p |
| NDV | miR-219a-1-3p > miR-3165 > miR-1257 > miR-758-5p > miR-7161-3p > miR-7844-5p > miR-6868-3p |
| HV | miR-374a-3p > miR-1250-5p > miR-6815-3p > miR-29b-3p > miR-29c-3p > miR-26a-5p > miR-6740-3p |
| NiV | miR-6817-5p > miR-6846-3p > miR-3173-5p > miR-1205 > miR-2054 |
| MV | miR-6770-3p > miR-4455 > miR-514b-3p > miR-6506-5p > miR-1183 > miR-6841-3p > miR-4714-5p |
| MuV | miR-6861-5p > miR-6877-3p > miR-143-5p > miR-4433a-5p > miR-1236-3p > miR-4781-3p > miR-219a-2-3p |
vSSs in Paramyxovirus genomes.
vSSs in bold are conserved across two or more viruses. vSSs identified Paramyxovirus genomes. vSMs in bold are conserved across two or more viruses. vSS hits are listed in order of their decreasing abundance amongst viral genomes. RSV, respiratory syncytial virus; hMPV, human metapneumovirus; NDV, New Castle disease virus; HV, hendravirus; NiV, nipah virus; MV, measles virus; MuV, mumps virus.
Small Viral RNA Transcripts Encode vSMs and vSSs
The results show that RSV infection of Vero cells, which lack a functional IFN α/β locus (Karnati et al., 2015; Alipoor et al., 2016), produces small RNAs from the G and L genes (Figure 4C). We also analyzed small RNAs transcriptomes of mock-treated, RSV-infected A549 cells, or RSV-infected neuronal SHSY5Y cells by NGS to identify if these transcripts encoded vSMs and vSSs. Transcripts mapping to viral genome over their entire length (36 nts) were also compared to mature human miRNAs to identify vSMs and vSSs. In RSV-infected A549 cells, we identified viral transcripts that encoded three overlapping full-length vSMs (miR-384, miR-6508-5p, and miR-642a-3p) and others that were perfectly complementary to miR-4483 (Figure 5A). Viral transcripts from RSV infected SHSY5Y cells encoded a near full length miR-3201 or miR-642a-3p transcript and vSSs against miR-642b-5p and miR-182-3p (Figure 5C). Many of transcripts encoded only vSMs or vSSs in contrast to the transcripts noted above which encoded nearly full-length host miRNAs (Figures 5B,D). These findings support our hypothesis that Paramyxovirus gene transcripts may modulate host miRNA activity via seed match or complementarity. While host encoded competing endogenous RNA are known, this is the first such analysis done for Paramyxovirus genomes thereby opening new areas of investigation.
Figure 5
Discussion
The low replicative fidelity of RdRp of Paramyxoviruses appears to be an evolutionary strategy since it allows for quasispecies generation, favors emergence of variants that can escape immune pressure, and may allow the virus to modulate host responses and gene expression via regulation of miRNA function. Tempering host gene expression is important for viral replication, and viruses have evolved several means to counter host responses to virus infection using mechanisms including shutdown of host transcription, translation, and modification of sncRNA and miRNA expression (Vishnoi and Rani, 2017; Olejniczak et al., 2018).
RIG-I and MDA-5 are prototypical pattern recognition receptors that detect and respond to presence of 5′ triphosphorylated and double-stranded RNA, respectively, during RNA virus infection (Piedade and Azevedo-Pereira, 2016; Auvinen, 2017). These receptors trigger a signaling cascade that culminates in establishment of an antiviral state in the infected and neighboring cells (Piedade and Azevedo-Pereira, 2016). Viruses avoid host responses to replicate and moderate these responses to facilitate replication. In addition to modifying the antiviral response, we propose that quasispecies generation during RNA virus replication helps regulate host gene expression by modulation of host miRNA function and activity. We speculate that RNA viral quasispecies produce a cloud of molecular signatures that mimic or inhibit host miRNA activity via sequence complementarity and alter the expression and function of several host pathways during infection. Our data suggest that vSMs and vSSs encoded by Paramyxovirus genomes and transcripts aid quasispecies generation by modulating various host pathways (Table 4) (Figure 6). It is likely that the generation of quasispecies is dynamic and adapts to host cell pressure. Thus, viral transcripts can mutate to produce vSMs or vSSs that modulate different miRNAs in different cell types as well as carry one or more frameshift mutations to produce alternative proteins/peptides that provoke modified or altered responses as these data suggest where viral transcripts can potentially bind to or mimic host miRNAs perturbing their expression and diverting host resources.
Table 4
| # of vSMs | # of vSSs | |||||
|---|---|---|---|---|---|---|
| Pathways deregulated by virus infection | RSV A | RSV B | Total | RSV A | RSV B | Total |
| Prostate cancer | 1 | 1 | 2 | 1 | 1 | 2 |
| Axon guidance | 1 | 1 | 2 | 1 | 1 | |
| PI3K-Akt signaling pathway | 1 | 1 | 2 | 1 | 1 | 2 |
| Chronic myeloid leukemia | 1 | 1 | 2 | 1 | 1 | |
| Non-small cell lung cancer | 1 | 1 | 1 | 1 | 2 | |
| Transcriptional mis-regulation in cancer | 1 | 1 | 2 | 1 | 1 | |
| Endometrial cancer | 1 | 1 | 1 | 1 | 2 | |
| T cell receptor signaling pathway | 1 | 1 | 2 | 1 | 1 | |
| Focal adhesion | 1 | 1 | 2 | 0 | ||
| Neurotrophin signaling pathway | 1 | 1 | 1 | 1 | ||
| Vasopressin-regulated water reabsorption | 1 | 1 | 2 | 2 | ||
| ErbB signaling pathway | 1 | 1 | 1 | 1 | ||
| Regulation of actin cytoskeleton | 1 | 1 | 0 | |||
| Prion diseases | 1 | 1 | 2 | 0 | ||
| Glioma | 1 | 1 | 1 | 1 | ||
| Endocytosis | 1 | 1 | 2 | 0 | ||
| Wnt signaling pathway | 1 | 1 | 2 | 0 | ||
| Colorectal cancer | 1 | 1 | 2 | 0 | ||
| Ubiquitin mediated proteolysis | 1 | 1 | 2 | 0 | ||
| Dopaminergic synapse | 1 | 1 | 2 | 0 | ||
| MAPK signaling pathway | 1 | 1 | 2 | 0 | ||
| Renal cell carcinoma | 1 | 1 | 1 | 1 | ||
| Melanogenesis | 1 | 1 | 1 | 1 | ||
| TGF-beta signaling pathway | 1 | 1 | 1 | 1 | ||
| B cell receptor signaling pathway | 1 | 1 | 1 | 1 | ||
| mRNA surveillance pathway | 1 | 1 | 2 | 0 | ||
miRs targeted by vSMs and vSSs for HMPV or RSV strain A or B were analyzed by DIANA miRPath server.
Significance of association based on p value ≤ 0.05 and False Discovery Rate (FDR) correction. Pathways with p-values < 0.05 are shown which are significant. Numbers indicate number of vSMs or vSSs.
Figure 6
Statements
Author contributions
AB and AM collected and analyzed the data and wrote the manuscript with RT.
Funding
Studies in the Tripp laboratory are supported by funds from the Georgia Research Alliance (GRA). RT is a GRA scholar. The funders had no role in design of the analyses in the manuscript. All views and opinions are of the authors only.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2019.00152/full#supplementary-material
Supplementary Table 1Lists of all raw data.
Supplementary Table 2Lists analysis of other RNA viruses that also show vSMs and vSSs.
References
1
AguilarH. C.LeeB. (2011). Emerging paramyxoviruses: molecular mechanisms and antiviral strategies. Expert Rev. Mol. Med.13:e6. 10.1017/S1462399410001754
2
AlbaneseM.TagawaT.BuschleA.HammerschmidtW. (2017). MicroRNAs of epstein-barr virus control innate and adaptive antiviral immunity. J. Virol.91:01667—16. 10.1128/JVI.01667-16
3
AlipoorS. D.AdcockI. M.GarssenJ.MortazE.VarahramM.MirsaeidiM.et al. (2016). The roles of miRNAs as potential biomarkers in lung diseases. Eur. J. Pharmacol.791, 395–404. 10.1016/j.ejphar.2016.09.015
4
AltschulS. F.GishW. (1996). Local alignment statistics. Meth. Enzymol.266, 460–480. 10.1016/S0076-6879(96)66029-7
5
AltschulS. F.GishW.MillerW.MyersE. W.LipmanD. J. (1990). Basic local alignment search tool. J. Mol. Biol.215, 403–410. 10.1016/S0022-2836(05)80360-2
6
AltschulS. F.PopM. (2017). Chapter 20.1. Sequence alignment, in Handbook of Discrete and Combinatorial Mathematics, eds RosenK. H.ShierD. R.GoddardW. (Boca Raton, FL: CRC Press/Taylor & Francis).
7
AmarasingheG. K.BaoY.BaslerC. F.BavariS.BeerM.BejermanN.et al. (2017). Taxonomy of the order mononegavirales: update 2017. Arch. Virol.162, 2493–2504. 10.1007/s00705-017-3311-7
8
AuvinenE. (2017). Diagnostic and prognostic value of microRNA in viral diseases. Mol. Diagn. Ther.21, 45–57. 10.1007/s40291-016-0236-x
9
BackN. K.NijhuisM.KeulenW.BoucherC. A.Oude EssinkB. O.van KuilenburgA. B.et al. (1996). Reduced replication of 3TC-resistant HIV-1 variants in primary cells due to a processivity defect of the reverse transcriptase enzyme. EMBO J.15, 4040–4049. 10.1002/j.1460-2075.1996.tb00777.x
10
BartelD. P. (2004). MicroRNAs: genomics, biogenesis, mechanism, and function. Cell116, 281–297. 10.1016/S0092-8674(04)00045-5
11
BartelD. P. (2009). MicroRNAs: target recognition and regulatory functions. Cell136, 215–233. 10.1016/j.cell.2009.01.002
12
BurgessD. J. (2017). Gene expression: principles of gene regulation across tissues. Nat. Rev. Genet.18:701. 10.1038/nrg.2017.94
13
CaiX.HagedornC. H.CullenB. R. (2004). Human microRNAs are processed from capped, polyadenylated transcripts that can also function as mRNAs. RNA10, 1957–1966. 10.1261/rna.7135204
14
ChenZ.LiS.SubramaniamS.ShyyJ. Y.ChienS. (2017). Epigenetic regulation: a new frontier for biomedical engineers. Annu. Rev. Biomed. Eng.19, 195–219. 10.1146/annurev-bioeng-071516-044720
15
CloonanN. (2015). Re-thinking miRNA-mRNA interactions: intertwining issues confound target discovery. Bioessays37, 379–388. 10.1002/bies.201400191
16
DamasN. D.FossatN.ScheelT. K. H. (2019). Functional interplay between RNA viruses and non-coding RNA in mammals. Non-Coding RNA5:7. 10.3390/ncrna5010007
17
DongF.LouD. (2012). MicroRNA-34b/c suppresses uveal melanoma cell proliferation and migration through multiple targets. Mol. Vis.18, 537–546.
18
FabianM. R.SonenbergN.FilipowiczW. (2010). Regulation of mRNA translation and stability by microRNAs. Annu. Rev. Biochem.79, 351–379. 10.1146/annurev-biochem-060308-103103
19
FinneganE. F.PasquinelliA. E. (2013). MicroRNA biogenesis: regulating the regulators. Crit. Rev. Biochem. Mol. Biol.48, 51–68. 10.3109/10409238.2012.738643
20
FirthA. E.BrierleyI. (2012). Non-canonical translation in RNA viruses. J. Gen. Virol.93, 1385–1409. 10.1099/vir.0.042499-0
21
FlorT. B.BlomB. (2016). Pathogens use and abuse microRNAs to deceive the immune system. Int. J. Mol. Sci.17:538. 10.3390/ijms17040538
22
FriedmanR. C.FarhK. K.BurgeC. B.BartelD. P. (2009). Most mammalian mRNAs are conserved targets of microRNAs. Genome Res.19, 92–105. 10.1101/gr.082701.108
23
GhosalS.DasS.SenR.ChakrabartiJ. (2014). HumanViCe: host ceRNA network in virus infected cells in human. Front. Genet.5:249. 10.3389/fgene.2014.00249
24
HaM.KimV. N. (2014). Regulation of microRNA biogenesis. Nat. Rev. Mol. Cell Biol.15, 509–524. 10.1038/nrm3838
25
HashimotoY.AkiyamaY.YuasaY. (2013). Multiple-to-multiple relationships between microRNAs and target genes in gastric cancer. PLoS ONE8:e62589. 10.1371/journal.pone.0062589
26
HopfieldJ. J. (1974). Kinetic proofreading: a new mechanism for reducing errors in biosynthetic processes requiring high specificity. Proc. Natl. Acad. Sci. U.S.A.71, 4135–4139. 10.1073/pnas.71.10.4135
27
HussainW.MahmoodT.HussainJ.AliN.ShahT.QayyumS.et al. (2019). CRISPR/Cas system: a game changing genome editing technology, to treat human genetic diseases. Gene685, 70–75. 10.1016/j.gene.2018.10.072
28
JiaS.ZhaiH.ZhaoM. (2014). MicroRNAs regulate immune system via multiple targets. Discov. Med.18, 237–247.
29
JohnsonT. R.GrahamB. S. (2004). Contribution of respiratory syncytial virus G antigenicity to vaccine-enhanced illness and the implications for severe disease during primary respiratory syncytial virus infection. Pediatr. Infect. Dis. J.23, S46–S57. 10.1097/01.inf.0000108192.94692.d2
30
KahnJ. S. (2006). Epidemiology of human metapneumovirus. Clin. Microbiol. Rev.19, 546–557. 10.1128/CMR.00014-06
31
KaikkonenM. U.AdelmanK. (2018). Emerging roles of non-coding RNA transcription. Trends Biochem. Sci.43, 654–667. 10.1016/j.tibs.2018.06.002
32
KarnatiH. K.PanigrahiM. K.GuttiR. K.GreigN. H.TamargoI. A. (2015). miRNAs: key players in neurodegenerative disorders and epilepsy. J. Alzheimer's Dis.48, 563–580. 10.3233/JAD-150395
33
KlinkeO.FeederleR.DelecluseH. J. (2014). Genetics of Epstein-Barr virus microRNAs. Semin. Cancer Biol.26, 52–59. 10.1016/j.semcancer.2014.02.002
34
KumariB.JainP.DasS.GhosalS.HazraB.TrivediA. C.et al. (2016). Dynamic changes in global microRNAome and transcriptome reveal complex miRNA-mRNA regulated host response to Japanese encephalitis virus in microglial cells. Sci. Rep.6:20263. 10.1038/srep20263
35
KuriharaY.WatanabeY. (2010). Processing of miRNA precursors. Methods Mol. Biol.592, 231–241. 10.1007/978-1-60327-005-2_15
36
LauringA. S.AndinoR. (2010). Quasispecies theory and the behavior of RNA viruses. PLoS Pathog.6:e1001005. 10.1371/journal.ppat.1001005
37
LiC.HuJ.HaoJ.ZhaoB.WuB.SunL.et al. (2014). Competitive virus and host RNAs: the interplay of a hidden virus and host interaction. Protein Cell5, 348–356. 10.1007/s13238-014-0039-y
38
MoralesS.MonzoM.NavarroA. (2017). Epigenetic regulation mechanisms of microRNA expression. Biomol. Concepts8, 203–212. 10.1515/bmc-2017-0024
39
MurashovA. K. (2017). RNAi and MicroRNA-mediated gene regulation in stem cells. Methods Mol. Biol.1622, 15–25. 10.1007/978-1-4939-7108-4_2
40
OlejniczakM.Kotowska-ZimmerA.KrzyzosiakW. (2018). Stress-induced changes in miRNA biogenesis and functioning. Cell. Mol. Life Sci.75, 177–191. 10.1007/s00018-017-2591-0
41
OshanskyC. M.KrunkoskyT. M.BarberJ.JonesL. P.TrippR. A. (2009). Respiratory syncytial virus proteins modulate suppressors of cytokine signaling 1 and 3 and the type I interferon response to infection by a toll-like receptor pathway. Viral Immunol.22, 147–161. 10.1089/vim.2008.0098
42
PanekJ.KolarM.HerrmannovaA.ValasekL. S. (2016). A systematic computational analysis of the rRNA-3' UTR sequence complementarity suggests a regulatory mechanism influencing post-termination events in metazoan translation. RNA22, 957–967. 10.1261/rna.056119.116
43
PangestiK.Abd El GhanyN. A. M.WalshM. G.KessonA. M.Hill-CawthorneG. A. (2018). Molecular epidemiology of respiratory syncytial virus. Rev. Med. Virol.28:e1968. 10.1002/rmv.1968
44
PapenburgJ.BoivinG. (2010). The distinguishing features of human metapneumovirus and respiratory syncytial virus. Rev. Med. Virol.20, 245–260. 10.1002/rmv.651
45
Perez-RodriguezF. J.D'AndreaL.de CastellarnauM.CostafredaM. I.GuixS.RibesE.et al. (2016). Improving virus production through quasispecies genomic selection and molecular breeding. Sci. Rep.6:35962. 10.1038/srep35962
46
PiedadeD.Azevedo-PereiraJ. M. (2016). The role of microRNAs in the pathogenesis of herpesvirus infection. Viruses8:156. 10.3390/v8060156
47
QinJ.LiW.GaoS. J.LuC. (2017). KSHV microRNAs: tricks of the devil. Trends Microbiol.25, 648–661. 10.1016/j.tim.2017.02.002
48
RimaB.CollinsP.EastonA.FouchierR.KurathG.LambR. A.et al. (2017). ICTV virus taxonomy profile: pneumoviridae. J. Gen. Virol.98, 2912–2913. 10.1099/jgv.0.000959
49
RitchieW.LegendreM.GautheretD. (2007). RNA stem-loops: to be or not to be cleaved by RNAse III. RNA13, 457–462. 10.1261/rna.366507
50
SchusterJ. E.WilliamsJ. V. (2014). Human metapneumovirus. Microbiol. Spectr. 2. 10.1201/b16778-33
51
ShuW.BoX.NiM.ZhengZ.WangS. (2007). In silico genetic robustness analysis of microRNA secondary structures: potential evidence of congruent evolution in microRNA. BMC Evol. Biol.7:223. 10.1186/1471-2148-7-223
52
SteinhauerD. A.DomingoE.HollandJ. J. (1992). Lack of evidence for proofreading mechanisms associated with an RNA virus polymerase. Gene122, 281–288. 10.1016/0378-1119(92)90216-C
53
SullivanC. S.GanemD. (2005). MicroRNAs and viral infection. Mol. Cell20, 3–7. 10.1016/j.molcel.2005.09.012
54
TrobaughD. W.KlimstraW. B. (2017). MicroRNA regulation of RNA virus replication and pathogenesis. Trends Mol. Med.23, 80–93. 10.1016/j.molmed.2016.11.003
55
VishnoiA.RaniS. (2017). MiRNA biogenesis and regulation of diseases: an overview. Methods Mol. Biol.1509, 1–10. 10.1007/978-1-4939-6524-3_1
56
WeillN.LisiV.ScottN.DallaireP.PellouxJ.MajorF. (2015). MiRBooking simulates the stoichiometric mode of action of microRNAs. Nucleic Acids Res.43, 6730–6738. 10.1093/nar/gkv619
57
XuZ. Z.MathewsD. H. (2016). Secondary structure prediction of single sequences using rnastructure. Methods Mol. Biol.1490, 15–34. 10.1007/978-1-4939-6433-8_2
Summary
Keywords
RNA viruses, microRNA, miRNA, miR, seed sequence, mimics, sponges
Citation
Bakre AA, Maleki A and Tripp RA (2019) MicroRNA and Nonsense Transcripts as Putative Viral Evasion Mechanisms. Front. Cell. Infect. Microbiol. 9:152. doi: 10.3389/fcimb.2019.00152
Received
07 February 2019
Accepted
23 April 2019
Published
08 May 2019
Volume
9 - 2019
Edited by
Jonatas Abrahao, Federal University of Minas Gerais, Brazil
Reviewed by
Eric Roberto Guimarães Rocha Aguiar, Federal University of Bahia, Brazil; Kai Huang, The University of Texas Medical Branch at Galveston, United States; Rafael Campos, State University of Ceará, Brazil
Updates
Copyright
© 2019 Bakre, Maleki and Tripp.
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: Ralph A. Tripp ratripp@uga.edu
This article was submitted to Virus and Host, a section of the journal Frontiers in Cellular and Infection Microbiology
Disclaimer
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.