article: We have used multiplexed high-throughput sequencing to characterize changes in small RNA populations that occur during viral infection in animal cells. Small RNA-based mechanisms such as RNA interference (RNAi) have been shown in plant and invertebrate systems to play a key role in host responses to viral infection. Although homologs of the key RNAi effector pathways are present in mammalian cells, and can launch an RNAi-mediated degradation of experimentally targeted mRNAs, any role for such responses in mammalian host-virus interactions remains to be characterized. Six different viruses were examined in 41 experimentally susceptible and resistant host systems. We identified virus-derived small RNAs (vsRNAs) from all six viruses, with total abundance varying from “vanishingly rare” (less than 0. 1% of cellular small RNA) to highly abundant (comparable to abundant micro-RNAs “miRNAs”). In addition to the appearance of vsRNAs during infection, we saw a number of specific changes in host miRNA profiles. For several infection models investigated in more detail, the RNAi and Interferon pathways modulated the abundance of vsRNAs. We also found evidence for populations of vsRNAs that exist as duplexed siRNAs with zero to three nucleotide 3′ overhangs. Using populations of cells carrying a Hepatitis C replicon, we observed strand-selective loading of siRNAs onto Argonaute complexes. These experiments define vsRNAs as one possible component of the interplay between animal viruses and their hosts.
  Biological systems are protected by innate immune mechanisms initiated by host sensors called pattern recognition receptors (‘PRRs’) that recognize specific “foreign” features of invading pathogens to initiate multiple downstream anti-pathogen cascades. PRRs that detect nucleic acid structures characteristic of viral infection (such as single- or double-stranded RNA or DNA) are among the innate responders that protect diverse cell types from viral pathogenesis (for review, see [1], [2]). How the cell handles viral double-stranded RNA (dsRNA) is of special interest because dsRNA is a necessary intermediate in the replication of RNA viruses. In addition to dsRNA that forms during replication of the virus genome, RNA duplexes can form due to self-complementarity in the virus genome, and in some instances, from sense-antisense transcription of overlapping genes. Four of the most studied families of PRRs for dsRNA are: (a) cytoplasmic RNA helicases like Retinoic acid-inducible gene I & Melanoma differentiation-associated gene-5 (“RIG-I” & “Mda-5, ” which trigger mitochondrial-localized antiviral pathways); (b) Protein Kinase R (“PKR, ” which induces a translational arrest state in cells after sensing dsRNA); (c) 2′–5′ oligoadenylate synthetase (“OAS, ” which stimulates the ssRNase activity of RNase L in response to dsRNA); and (d) Toll-like receptors (“TLRs, ” which bind various forms of RNA or DNA). All of these PRRs trigger the Interferon (IFN) responses, and activate IFN-stimulated genes (ISGs) that establish an antiviral state in the infected cell (for review, see [3]). The IFN signaling pathway is central to the detection of, and response to, virus infections in cells. Type I IFNs (IFN-α and IFN-β) make up one of the first lines of defense in the innate immune response to viruses by inducing antiviral ISGs, modulating the levels of specific host-encoded miRNAs [4], and in a feedback loop, that of PKR and OAS. Many viruses are also susceptible to treatment with Type I IFNs, and conversely, cells that have higher basal activity of ISGs seem to mount a more successful antiviral response, and are not targeted by viruses [5]. Dicer is another PRR that recognizes dsRNA, chopping it into smaller duplexes called siRNAs that are 19–27 nucleotides (nt) long [6], [7]. These siRNAs have a terminal 5′ mono-phosphate and a terminal 3′ hydroxyl on both strands, generally have 2 nt 3′ overhangs, and are fed into an RNA-induced silencing complex “RISC” (for review on Dicer and Argonautes, see [8], [9]). siRNA duplexes are unwound, and only one strand remains associated with RISC (the mechanism of unwinding and choice of strand is poorly understood; for review, see [10]). One of the key components of RISC is a protein called Argonaute-2 (Ago-2), which belongs to the Argonaute family of proteins. Ago-2 is the only member of the family that has cleavage activity, and is the designated ‘slicer’ protein in RISC that mediates cleavage of mRNA in a sequence-directed manner by a process termed RNA interference, or ‘RNAi’ [11], [12], [13], [14]. There is strong evidence for an antiviral role for RNAi in plant and invertebrate systems (for review, see [15], [16], [17]). Viruses replicate most effectively in these systems in the absence of key elements of the RNAi pathway: either in cells lacking components of the RNAi machinery, or in the presence of virus-encoded suppressors of the silencing pathway (for review, see [18], [19]). As expected, virus-derived siRNAs (vsRNAs) can be detected in some plant and invertebrate systems that are capable of mounting a successful/partially successful RNAi response [15], [16], [17]. A population of vsRNAs would be an expected component of any viral defense pathway that acted through an RNAi mechanism. In mammalian cells, short duplex RNAs can effectively enter the RNAi pathway and function in sequence-specific silencing, while duplexes longer than 30 nt generally produce a more complex response including the induction of multiple non-specific pathways including the IFN response (for review, see [20], [21]). Indeed, RNA and DNA viruses have evolved a host of defense mechanisms to counteract the nonspecific signaling effects of dsRNA. For example, Adenovirus VA RNA sequesters PKR [22], while proteins from Vaccinia virus (E3L), Porcine Rotaviruses (NSP3), and Influenza A virus (NS1) sequester dsRNA and prevent stimulation of the IFN response [23], [24], [25], [26]. Viral proteins can also inhibit signaling downstream of dsRNA binding, as in the case of the HCV protease NS3/4A, which cleaves IPS-1 (the RIG-I/MDA-5 signaling partner) to consequently disrupt induction of IFN responses [27]. Several of these dsRNA-binding proteins may also facilitate viral evasion of host immune responses by inhibiting RNAi [28]. Additionally, some viruses make their genomes inaccessible to PRRs of various types including IFN effectors and the siRNA-programmed RISC complex (e. g. [29]). Viruses may also perturb another class of effectors involved in RNAi called micro-RNAs (miRNAs), which are a class of cellular small RNAs generated by Dicer from hairpin structures. Cellular miRNA profiles are frequently modulated upon infection by viruses, and this may contribute in some cases to infectivity and pathogenesis [30]. Conversely, some viruses usurp the host miRNA machinery for processing miRNA-like structures encoded in the viral genome, potentially using these molecules for regulation of virus/host gene expression [31]. With so much potential for RNA-mediated cross talk between the IFN response, the RNAi pathway, and the virus itself, it has been difficult to demonstrate a precise role for the RNAi pathway in vertebrate antiviral defense. The difficulties in segregating IFN and RNAi functions have given rise to speculations that the antiviral role of RNAi may have been lost during evolution, or alternatively, that RNAi-based defense may only be harnessed by triggers such as short hairpins and siRNAs that do not stimulate the IFN pathway. There has been some attempt at demonstrating recognition of viral RNA by the RNAi machinery. For instance, in Vero cells (which lack IFNα/β), inhibition of RNAi by Dicer knockdown increases replication of an RNA virus, the Influenza A Virus [32]. Additionally, there are cases where short virus-derived RNAs can be detected in vertebrate systems (e. g. from HDV [33] by high-throughput sequencing, and the HCV replicon [34], by bulk analysis methods). However, it is still not clear how general the presence of such RNAs is, and whether these RNAs can participate in host defense mechanisms. To complicate this issue, many of the classically-studied virus-host systems have been chosen based on the ability of the virus to rapidly replicate and kill host cells; these experimental infection systems may artificially under-represent the capacity of vertebrate cells to protect themselves, hence biasing against systems where RNAi might have a significant role in host-virus interactions. Here, we sought a broader survey of potential RNA-derived defenses in viral infection systems. Given no knowledge of which virus type might engage the RNAi machinery, and which cell types might efficiently use this machinery in defense, we cast a wide net in terms of both virus families and host cells. In this study, we describe small RNA populations from six different RNA viral pathogens, each in a variety of animal cell infection systems (including both immune-competent and immune-compromised hosts). Upon examining small RNA populations from ∼150 samples with sample-specific DNA barcodes, we found viral-derived small RNAs (vsRNAs) from each virus, with vsRNA populations sensitive to both viral and host characteristics. A more detailed analysis of vsRNAs in two viral infection models (Hepatitis C Virus and Poliovirus) in various host types revealed that multiple distinct pools of vsRNAs may co-exist during infection: as single strands, as part of duplexes, and in complexes that may contain Argonautes. We also observed specific changes in cell-derived miRNA populations, providing a clear indication of host perturbation by the virus. The characterization of small RNA populations during RNA virus infections provides both an experimental entry point, and an indication of the complexity that will need to be addressed in understanding roles for small RNAs in host and viral processes.
  In the following sections, we will describe small RNA populations present during infection of animal cells with six different viruses. In each case, we have taken infected cells, extracted small RNA populations, and characterized these populations using high-throughput sequencing methods. Two high-throughput sequencing platforms were used: Roche/454 pyrosequencing (http: //www. 454. com/), to obtain several hundred thousand sequences from pools of appropriately linkered amplicon templates; and Solexa/Illumina technology, which yields larger datasets of shorter reads (http: //www. illumina. com/). Due to the large number of samples to be analyzed, we used DNA barcodes to ‘tag’ RNA samples from individual experiments, which facilitated sequencing in parallel from multiple samples. This allowed us to work with samples from different viral systems and diverse experimental conditions in a cost-effective manner, with a small number of instrument runs. Viral-derived sequences were identified in sequence datasets through pattern matching using standard software (BLAT [35] and BLAST [36]). We use the term ‘vsRNA’ to refer to small RNA segments whose sequences show perfect complementarity to the infecting viral genome at every base position (reference genomes listed in Table S1). vsRNAs are distinct from host-derived miRNAs that may show partial complementarity to sites in the viral genome (e. g. [37], [38], [39], Fig. S1). We detected 77,609 vsRNAs out of 19,425,777 sequences from 151 datasets (Fig. 1, Figure S2: length distributions). The most abundant vsRNAs from each virus are listed in Table S2. For a small number of vsRNAs (0. 033%) we observed a perfect match to both the host and viral genomes (Table S3, Table S4). The fractions of vsRNAs that matched host genomes were approximately as expected by random sequence coincidence (for example: the human genome, with a unique genome complexity of 2×109 bp, would match approximatly 1 in 4000 arbitrary 22-mer sequences). The perfect nature of the homology makes it difficult to determine whether this minor class of sRNAs was derived from the host or from the virus. Furthermore, to validate the specificity of the barcoding and sequencing assays, we carried out sequence comparisons to the full set of viruses for each experimental sample. We identified 13 vsRNAs that were ‘rogue’ hits i. e. mapped to one of the other 5 viruses not used in that particular experiment. In no sample were the ‘rogue’ matches present at more than 0. 008% of all parsed sequences (Table S3). For certain purposes, it will be of interest to compare vsRNA incidences in different samples. Such comparisons require some normalization for total depth of RNA sequencing. In Table S3, we provide two distinct normalizations for each sample: normalization to total small RNAs recovered and sequenced (v/sRNA), and normalization to the population of cellular miRNAs that are expected to represent a large proportion of bona-fide small RNA effectors (v/miR; miRNAs are defined as documented in miRBase ver9. 2 [40], [41], [42]). There is a substantial challenge in choosing and interpreting appropriate normalization schemes: any change in sample character that results in increased levels of non-specific degradation of RNA will increase the levels of non-specific decay products (which may include decay products of both cellular and viral long RNAs), and impact both v/miR and v/sRNA ratios. Another important consideration is whether the difference in v/miR (or v/sRNA) ratios between two samples is above the variance in ratios observed between technical replicates. For all relevant technical replicates in our analysis, the variances in v/miR and v/sRNA ratios were <3-fold and <1. 5-fold, respectively. These notes provide caution in interpreting small differences in normalization values between samples. In describing the results of this work, we have taken care to avoid any a-priori assumption that small RNAs identified by sequencing play a functional role in gene silencing, viral pathogenesis, or host response. In the Discussion section, we will summarize arguments pertaining to this question. Components of the worm RNAi machinery such as the argonaute, rde-1 [43], the dsRNA binding protein, rde-4 [44], [45], and the RNA-dependent RNA Polymerase or RdRP, rrf-1 [46] are essential for protection against Vesicular Stomatitis Virus ‘VSV’ [47], and Flock House Virus ‘FHV’ replication [48]. To characterize small RNA populations in an animal system known to utilize the RNAi machinery in antiviral defense, we used C. elegans experimentally infected with FHV RNA1ΔB2 (FHV RNA1 that expresses a mutant version of the RNAi suppressor protein, B2; [48]). Two different vsRNA capture and library production schemes were used to enrich for Dicer products or for RdRP products, both of which have structures distinct from those of RNA fragments generated by alkali-induced degradation. The first (5′-phosphate-dependent cloning) requires a single phosphate at the 5′ end of the RNA, and allows for the capture of Dicer products (which have a mono-Phosphate and a hydroxyl moiety at their 5′ and 3′ termini). The second (5′-phosphate-independent cloning; [49]) is designed to capture RNA populations with any number of 5′ phosphates (zero, mono, di, tri), including both RdRP products (which have a tri-Phosphate and a hydroxyl moiety at their 5′ and 3′ termini) and Dicer products. Both procedures require a 3′ end that can ligate to a pre-adenylated linker, and allow for the capture of 3′-OH and 2′-O-Methyl structures but not 3′ phosphate termini, thus minimizing the extent of capture of degradation products (many of which have 3′ mono-phosphate termini). 5′ mono-phosphorylated (5′-P) vsRNAs were present during abortive FHV RNA1ΔB2 replication in wild-type animals (v/miR = 0. 007; Fig. 2B). vsRNAs were absent in two RNAi-defective mutants, rrf-1 (pk1417) I and rde-4 (ne299) III (Table S3), while as predicted, genomic viral RNA replicated to high levels in these mutants (Parameswaran P, unpublished). Similarly, vsRNAs were much reduced (19-fold; P-value = 2. 3E-227) in rde-1 (ne300) V mutants (Fig. 2C). We also observed a difference in strand ratios of vsRNAs (Positive∶Negative) between strains: 1∶2. 4 in wild-type, versus 1∶1. 1 in the mutant, rde-1 (P-value = 0. 0016). The population of RNAs captured with no requirement for a 5′-P terminus (i. e. 5′-xP RNAs) yielded a stronger signature for vsRNAs in wild-type worms with replicating RNA1ΔB2 (v/miR = 0. 019; Fig. 2F). Fewer vsRNAs mapped to the positive strand of FHV than to the negative strand, with a Positive∶Negative vsRNA strand ratio of 1∶3. 5 (P-value = 1. 1E-48). rde-4−/− was the only RNAi-defective mutant that yielded a detectable signature for 5′-xP vsRNAs (v/miR = 0. 0014), with a strand ratio (Positive∶Negative) of 1. 3∶1 (Fig. 2G). Interestingly, in wild-type worms, both 5′-P and 5′-xP vsRNAs were distributed throughout the length of the genome, with increased frequencies of positive-strand vsRNAs detected in the 3′ region that also encodes the subgenomic RNA species RNA3 (Fig. 2B, 2F). To identify virus-host systems in which RNAi might participate as an antiviral defense mechanism, we sequenced small RNAs from diverse populations of cells (of human or mouse origin) infected with one of five viruses: Vesicular Stomatitis Virus (VSV), Poliovirus, West Nile Virus (WNV), Dengue Virus, or Hepatitis C Virus (HCV). These viruses were purposefully chosen as token members of diverse families (Fig. 1), and are mostly positive-stranded (except for Vesicular Stomatitis Virus, which is negative-stranded). We identified vsRNAs from all six surveyed viruses (Fig. 1; Table S3), albeit in only a fraction of all infected samples investigated. From this initial survey, we made a choice of a single host-virus system in which to further investigate vsRNA biogenesis. The viral system chosen for this purpose was HCV infection of human Hepatoma cells. While the remainder of the Results section will focus primarily on HCV, we will briefly summarize our observations in the four other virus systems. For the Polio, VSV, West Nile and Dengue (Fig. S3) systems (Table S3), the abundance and molecular features of vsRNAs were dependent on the nature of the host and/or the virus, with some notable trends: A more detailed description of small RNA profiles from West Nile Virus, Dengue, Vesicular Stomatitis Virus, and Poliovirus is provided in the supporting document (Text S1), and in Supplementary Tables & Figures. HCV is an enveloped, positive-stranded RNA virus that is a member of the Flaviviridae family. Its genome is flanked by short stretches of structured RNA in the 5′ and 3′ UTRs, is uncapped, and lacks a 3′ poly-A tail. A previous study with HCV-1b-infected Huh7. 5 cells failed to identify HCV-derived vsRNAs using standard sequencing protocols [57]. We expanded on this work by choosing two cell-culture-based systems that are used for studying HCV replication: an HCC cell line (Huh7) harboring a subgenomic replicon of genotype 1b [55], and an infectious virion system (Huh7. 5 cells infected with tissue-culture-produced virions of genotype 2a [58]). v/miR levels of 5′-P vsRNAs from replicon cells varied between 0. 03 and 0. 14 (Fig. S8, Table S3), with an estimated 7300 +/− 2200 vsRNA molecules per cell (based on the approximation that the most abundant miRNA, miR-122a, is present at 15,000 copies per hepatoma cell in culture [59]). In virus-infected Huh7. 5 cells, we detected a very low incidence of vsRNAs at early time points, with an increase over time (Fig. S9, S10). vsRNAs were found starting at 1 day post-infection (dpi) in Huh7. 5 cells (v/miR = 0. 000025), and steadily increased (excluding a possible dip at 11dpi), reaching a v/miR value of 0. 056 at 15 dpi. vsRNAs from the sense (positive) and antisense (negative) viral strands were roughly equally abundant in both the replicon and the infectious virion systems (sense-to-antisense ratios of 1. 07 to 1. 9; Fig. 5, Table S3). This contrasts with the observed ratios of genome-length viral RNAs, where the sense strand is 5- to 10-fold more abundant in replicon-harboring cells and in infected hepatocytes [53], [54], [55]. The near-equivalent abundance of vsRNA strands is consistent with vsRNAs deriving from cleavage of a double-stranded replication intermediate, which has an equimolar ratio of positive and negative strands. HCV vsRNAs from both the 1b and 2a genotypes were distributed throughout the length of the genome, with several ‘hotspots, ’ where many vsRNAs were found clustered in specific regions of the genome (Fig. 5). Direct comparison of vsRNA distributions, and of individual vsRNA hotspot species between the replicates demonstrated that both were reproducible properties of HCV infection (Fig. S11A–D). Additionally, the ability of structured sequences such as those found in the HCV IRES and EMCV IRES to produce specific small RNA populations is of considerable interest. A comparison of vsRNA localization and published secondary structures of the HCV IRES and EMCV IRES is shown in Fig. S12 & S13. We also found some evidence for nucleotide bias among HCV replicon (HCVrep) -derived positive-strand and negative-strand vsRNA populations, including a bias toward strings of Cs and Gs at the 5′ and 3′ termini respectively (Fig. S14). This potentially creates favorable conditions either for intramolecular base pairing within a vsRNA, or for base pairing between overlapping sense-antisense vsRNA pairs at their termini. We further investigated the potential for duplexed structures of sense and antisense vsRNAs from HCVrep, by comparing sequence placement for sense and antisense vsRNAs within the viral genome. There are examples of independently captured sense and antisense HCVrep-derived vsRNAs that could derive from a dsRNA duplex with a 0–3 base 3′ overhang (Fig. 6A–C). These are similar to canonical overhangs in Dicer-generated siRNAs. If we first separated HCVrep-derived vsRNAs into different size ranges and then calculated the distribution of overhangs, duplexes formed by overlapping sets of 20–21 nt sense and antisense vsRNAs had a strong bias for one or two nt 3′ overhangs (Fig. 6A). On the other hand, duplexes formed by vsRNAs that are 24–26 nt long have a wider overhang range of zero to three nucleotides (Fig. 6B; False Discovery Rate is less than 0. 01%). To explore the possibility that vsRNAs may associate with core components of the RISC machinery (the Argonaute, or “Ago” proteins) despite our inability to detect a role for vsRNAs in silencing pathways, we used transient transfection to express FLAG/HA-tagged Ago-1, Ago-2, Ago-3 or Ago-4 [12] in HCVrep cell lines. We note a limitation of the Argonaute immunoprecipitation (IP) assays in that a large fraction of small RNAs from the cell may be capable of associating with Argonautes in a specific or non-specific manner; nonetheless, the expectation of such experiments is that immunoprecipitation will lead to enrichment for small RNAs that specifically associate with the tagged Argonaute. We compared RNA populations from each of the four Argonaute IPs to RNAs from the Mock-IP (i. e. IP with FLAG Ab, using lysates from mock-transfected cells), to give us an indication of the specificity of the IPs, and conversely, of the degree of non-specificity due to “stickiness” of the α-FLAG-M2 Antibody (Fig. 7A). Specifically, we compared the enrichment for vsRNAs in the Ago IPs (relative to Mock IPs), first to the enrichment for RNAs previously known to be Ago-associated (miRNAs, some miRNA*s) [12], [60], [61], [62], and second to the de-enrichment for RNAs which have less (or no) specific association with Argonaute (ribosomal RNAs). In the Ago IPs, we observed a several-fold enrichment for vsRNAs (similar to that observed for miRNAs), accompanied by a marked de-enrichment for rRNA fragments (Fig. 7A; for raw data, see Fig. S15B–C). This indicates that at least a subpopulation of vsRNAs associates specifically with all four Argonautes. A striking feature of Ago association in general is the rapid reduction of the initial dsRNA duplex to a single-stranded guide RNA [63]. We compared the duplex properties of vsRNA populations in total cell lysates to those of vsRNAs that are specifically associated with the Argonautes (Fig. 7B). The IP datasets for Ago-2 and Ago-4 showed a notable feature: a striking de-enrichment for duplexes with 0–3 nt 3′ overhangs, compared to their respective total RNA samples (P-values of 0 and 2. 1E-49 respectively). Ago-1 IP and Ago-3 IP showed a de-enrichment for such duplexes, but total RNA samples for Ago-1 and Ago-3 did not have sufficient sequence coverage to allow for a comparison. These data suggest that HCVrep cell lysates have populations of duplexed (and some single-stranded) vsRNAs, with only a single strand of each duplex reproducibly incorporated into an Ago complex.
  We consider two possible sources for the vsRNA populations that were observed during infection: (i) the vsRNAs could be participants in a specific pathway (or pathways) in which small RNAs are generated from the viral genome for host or viral functions; and (ii) the vsRNAs could be products of non-specific degradation of longer (e. g. full-length or subgenomic) viral RNAs mediated by ssRNA nucleases, chemicals, pH, mechanical shear etc. We note that the small RNA populations characterized by sequencing may be a mixture of (i) biologically relevant small RNAs, and (ii) degradation products with limited significance. In particular, any population of larger RNAs, on extraction and experimental manipulation, can yield a sub-population of RNAs in every size range, including the miRNA and siRNA size range of 19–30 nt. Since viral genomic RNAs and mRNAs are abundant in infected cells, we would certainly expect degraded derivatives to contribute to sequenced pools. Despite the likely capture of some degradation products, there are several strong indications of vsRNA populations that are not simply the result of degradative mechanisms. The above characteristics of vsRNA populations strongly argue that mammalian cells retain the ability (present in lower organisms; for review, see [16], [17]) to utilize vsRNAs as RNAi effectors. As for any host-virus interaction, the expectation would be that the utilization of small RNA-based mechanisms would be highly dependent on the biology of the virus and the host. This is evident from observed differences in the strandedness of vsRNAs in different systems infected with the same virus, and may be accounted for by several factors such as (a) substantial contributions from potential degradation mechanisms; (b) nuclease activity on ssRNA templates, especially if the secondary structures in one strand are targeted preferentially; (c) strand accessibility; (d) Dicer processivity and activity, and (e) selective RISC loading. Some infectious systems (e. g. Poliovirus or VSV in HeLa cells; Fig. S18B, S17C) yield a vsRNA profile with a skew towards positive strand vsRNAs. By contrast, productive infections with the same viruses in other host environments (e. g. Poliovirus in mouse muscle (Fig. 4C, S18H) or VSV in MEF/BHK cells (Fig. 3D–E, S17B) can show strong signatures from both strands, consistent with potential generation by dsRNA-related mechanisms. Host-virus interactions may also determine bulk abundance of vsRNA populations. In the systems we surveyed, vsRNA abundance varied between 0 and 12. 8 (v/miR; relative to miRNAs), or between 0 and 0. 02 (v/sRNA; relative to all small RNAs; Table S3). Even within a single host, tissue specific factors seem to govern the efficacies of small RNA-related mechanisms. This was observed in IFNαβR−/− mice infected with West Nile Virus: even though levels of full-length viral RNA in the brain were comparable to levels in spleen & lymph nodes [5], [50], vsRNAs in the brain were undetectable (454-141,454-144; Table S3). Despite the complex nature of the factors that govern vsRNA abundance, a number of trends are suggested by comparison of vsRNA levels in paired “wild-type” and “mutant” hosts infected with various viruses. In each of the six sets of experiments where we compared vsRNA levels in parallel infections in IFNαβR (+/−), or ago-2 (+/−) hosts, we observed an increase (ranging from 1. 7-fold to >30-fold) in vsRNA levels in the infected “mutant” host. For five out of these six comparisons, the variation in vsRNA abundance between “wild-type” and “mutant” hosts is higher than the (maximum) 3-fold difference we observed in technical replicates. Viruses can hijack the RNAi machinery at various levels: at the level of Dicer, Argonaute, RISC-mediated silencing, or a combination of the above (for review, see [19]). Downregulation of Dicer has an attenuating effect on Hepatitis C Virus replication [88]. Several studies have shown that HCV proteins inhibit Dicer and Argonaute [34], [89], [90]. In these experiments, inhibition is not complete (∼60–70%; [89]), concordant with our ability to detect vsRNA populations with structures consistent with synthesis by Dicer. We also see evidence for loading of vsRNAs onto Ago complexes, indicating that some downstream steps are not entirely impaired. Key players in the piRNA pathway, Piwi and Aubergine, are required for protection against Drosophila X Virus infections [84]. The PIWI-piRNA pathway produces 26–31 nt RNAs in a Dicer-independent manner, and is mostly active in the germline [91], and in adjacent somatic tissues [92], [93]. We note that animal systems with PIWI proteins (vertebrates), we observe a wide size range for vsRNAs (Fig. S2). These diverse size classes, together with our observation that vsRNAs are still present in systems that lack Dicer (Poliovirus infections in dcr−/− MEFs), suggest that the Dicer pathway (which is thought to primarily produce RNAs shorter than 27 nt) may not be the only source for these vsRNAs, and that there may some contribution from baseline levels of PIWI proteins (or other novel proteins) in the various systems. We also identified a population of Polio-derived vsRNAs in dcr-1−/− and in eri-1+/+ MEFs that can form duplexes with 9 or 10 bp overhangs (Fig. S16I, S16L), which are hallmarks of piRNAs that are formed by a ping-pong mechanism [94], [95]. This observation brings up the question of how the piRNA pathway interfaces with the RNAi pathway during viral infections, and whether in the absence of the RNAi (and perhaps the IFN) machineries, we could uncover a role for the piRNA pathway in antiviral immunity. Long dsRNAs, such as those found in the viral replication intermediate, primarily induce the non-specific interferon response in mammalian cells. In the absence of the robust IFN pathway, long dsRNA becomes a trigger for the sequence-directed RNAi pathway [96], [97], [98], [99]. Accordingly, we detected a more abundant signature for dsRNA-derived vsRNAs in IFNαβR−/− mice (Poliovirus and West Nile Virus; Table S3). Conversely, when we stimulated the Interferon pathway by providing an exogenous supply of IFN-α to a culture of HCVrep-harboring cells, we found that concurrent with reduction in full-length genomic RNA (as reported in: [55], [100]), HCV-derived vsRNAs also dropped several fold (Fig. S20). Thus vsRNA abundance seems to associate with the strength of the IFN response: the more robust the IFN response, the fewer the number of vsRNAs. The observed boost in vsRNA abundance in IFN knockout conditions could conceivably reflect a number of distinct effects including augmented viral replication in the absence of IFN, and/or specific interactions between IFN stimulation and the RNAi machineries. We observed consistent effects of viral infection on miRNA profiles in distinct experimental systems (Fig. S7). For example: Downregulation of the miR-17-92 cluster (i) in virus-infected cells may be a pro-apoptotic indicator, since an increase in expression of the miR-17-92 cluster is associated with an inhibition of apoptosis [101], [102], [103]. A decrease in miR-125b (ii) is observed post-LPS-stimulation of macrophages, and causes de-repression of TNF-α in a sequence-specific manner [104]. We speculate that miR-125b may be one of the regulators of the TNF-α response during viral infection, and that regulation of miR-125b may require an intact IFN response. Increased levels of the anti-apoptotic miR-21 (iii) were found in memory and effector T cells, compared to naïve T cells [105]. It is tempting to speculate that this upregulation of miR-21 may be indicative of proliferating immune cells post-recognition of viral antigens. The identification of multiple vsRNAs, some of which are derived from ‘hotspot’ locations in diverse viral genomes may be useful for designing cocktails of siRNAs for therapeutic purposes, and for mapping areas of the viral genome that are more susceptible to RNAi machineries. Much work still remains to be done in designing siRNA duplexes such that the most accessible strand of the virus may be successfully targeted by siRNAs. A major concern is whether RNAi would be effective in combating viruses with fast replication kinetics. Poliovirus [106] and Semliki Forest Virus [107] replication in cultured cells have been effectively attenuated by an exogenous supply of siRNA triggers against the virus. Whether this holds true for clearing infections in whole organisms remains to be tested.
  Mice used for experimental infections with West Nile Virus were genotyped and bred in the animal facilities of the Washington University School of Medicine, and experiments were performed with approval from, and according to the guidelines of, the Washington University Animal Studies Committee (which is IACUC approved). For infections with Poliovirus, mice that express the human Poliovirus Receptor gene were maintained in BSL-2 animal facilities at Stanford University. The methods for mouse use and care were approved by the Stanford University Administrative Panel on Laboratory Animal Care (APLAC), and are in accordance with the USDA Animal Welfare Act and the Public Health Service Policy on Human Care and Use of Laboratory Animals. Worms and mouse tissues were flash frozen in liquid nitrogen, powdered and lysed. Cell lines were trypsinized and washed in PBS pre-lysis. The mirVana kit (Ambion) was used for isolation of RNAs shorter than 200 nt from all samples. Different protocols were used to prepare libraries of RNAs with mono-phosphorylated, or with modified 5′ termini. Briefly, in the 5′-P-Dep protocol, RNA was linkered at the 3′ terminus, size-selected, linkered at the 5′ terminus, reverse-transcribed, amplified using ten-nucleotide barcoded PCR primers [115], and sequenced on the Roche/454 GS-20 or the GS-FLX platforms. In one version of the 5′-P-IND protocol (used for all samples sequenced on the GS-20/FLX), the RNA was linkered on the 3′ end, size-selected and reverse-transcribed before addition of the second linker [49]. In the 5′-P-IND protocol for preparing Solexa libraries, RNA was linkered on the 3′ end, dephosphorylated and re-phosphorylated (to replace any multi-phosphate moieties with a monophosphate), linkered on the 5′ end, reverse-transcribed and amplified (method courtesy of Guoping Gu). All 5′-P-Dep and 5′-P-IND libraries for Solexa were prepared by introducing four-nucleotide barcodes as part of the 5′ linker, rather than during PCR amplification (Lui WO, Parameswaran P; unpublished). Shorter barcodes allowed for the use of non-barcoded primers for amplification, and most importantly, for greater allocation of sequence space to the sequence of interest. After preparation of barcoded libraries, the libraries were pooled together in molar ratios that were proportional to the sequencing depth required from each sample. The presence of a large number of libraries required multiple sequencing runs on the Roche/454 and the Illumina sequencers. Sequences obtained from the Roche/454 platform were handled differently from those obtained on the Solexa platform due to different amplicon structures. Sequences from the Roche/454 platform were binned based on their barcodes using Barsort [115], trimmed using perl scripts (PP) and aligned using a local copy of the multiple alignment program Blast (word size = 11). Sequences from Illumina' s platform were segregated into individual datasets based on perfect barcode match, trimmed to ensure removal of the flanking adapter sequences before analysis, and aligned using Blat (tile size = 11; step size = 5, run on Mac OS X). Alignments were performed to databases of species-specific miRNAs, and of the various viral genomes. The alignments were subsequently parsed to yield unique hits with the highest homology for each matching read. We filtered for matches of >16 nt for Roche/454 sequencing, and of >19 nt for Solexa sequencing. For comparing incidence of miRNAs across samples, the frequency of miRNAs and standard error were computed as per the following formulae [116]: The distribution, length and strandedness of vsRNAs (i. e. RNAs that mapped to the genome of the infecting virus) were plotted as a function of the length of the viral genome. Starts and ends of positive-strand and negative-strand vsRNAs were then compared to generate a matrix of the percent incidence of different types of overhangs. For generating the nucleotide bias, the vsRNAs were aligned at the 5′ termini (or at the 3′ termini), and the nucleotide frequencies were counted at each position ten nucleotides upstream and ten nucleotides from the 5′ end of the vsRNA population (or ten nucleotides downstream and ten nucleotides from the 3′ end of the vsRNA population). These numbers were then fed into Pictogram (Chris Burge, MIT; http: //genes. mit. edu/pictogram. html), and were normalized to the total numbers of A, C, G and T in the population of cloned vsRNAs, thus circumventing any bias that arose from a skewed ratio of nucleotides inherent to either the cloning process, or to the genome of the virus. Calculation of P-values: For establishing statistical significance, two-tailed P-values were calculated using either the Z-test for two proportions, or Fisher' s Exact test. Fisher' s Exact test was used in test cases with small sample or population sizes, while the Z-test was used if sample or population sizes were large. Calculation of False Discovery Rate (FDR): Populations of 20,21-mers and 24,25,26-mer HCV vsRNAs were pooled, and positive strand and negative strand vsRNAs were randomly selected from the pool. 10,000 datasets that mirrored positive strand and negative strand vsRNA abundances from the original 20,21-mer population, and another 10,000 datasets that mirrored vsRNA abundances (of both polarities) from the original 24,25,26-mer population were generated. The percent of randomly generated datasets wherein 1–2 nt 3′overhang duplexes were as abundant as in the original 20,21-mer population, or as depleted as in the 24,25,26-mer population was calculated (for example, FDR of 0. 01% indicates that we were unable to detect percent incidence of duplexes similar to what was observed, in at least 10,000 randomly-generated datasets). Cultures of HCVrep cells were grown to ∼80% confluency in 10-cm. plates, and transfected with 20ug each of FLAG/HA-tagged Ago-1, Ago-2, Ago-3 or Ago-4 (codon-optimized), and harvested 24 hours post-transfection. Cells were washed twice with 4mL PBS, and 0. 75mL of ISOB/NP40 (10mM Tris pH 7. 9,0. 15M NaCl, 1. 5mM MgCl2,0. 8% NP40, proteinase inhibitor at 1 tablet per 13mL solution) was added to each plate. Cell lysates were vortexed, incubated on ice for 20 minutes, spun at 4°C for 10 min at 13,000 rpm, and the supernatants were transferred to new tubes. 10 uL ‘Fake’ (uncoated) Sepharose 4B beads were washed twice with 1. 5ml NET-1 buffer (1×TBS-0. 2% Tween), and incubated with supernatants at 4°C for 2 hrs. The supernatants post-‘fake-bead’ IP were subsequently incubated with 10uL of washed Anti-Flag M2 beads [Sigma] at 4°C for 4 hrs. Finally, beads were washed 3× with 500uL NET-1 buffer, and RNA was extracted with lysis buffer (mirVana kit, Ambion). Transfections were performed in duplicate: 1 plate was used as input for IPs, while the other was used for cloning of total small RNA populations. For Mock-IPs, lysates from mock-transfected cells were used as described above. The loading control was a 5. 7 kb region of the pGEM plasmid carrying the Poliovirus 1 genome (pGEM-PV1), obtained by digestion of the plasmid with HindIII and AgeI. Equivalent amounts of RNA (as measured using a NanoDrop) from each sample were run on a 1% agarose-formaldehyde gel, and transferred under basic conditions to a Hybond+ membrane. The loading control was visualized by staining the blot with methylene blue. 60-mer DNA probes to specifically detect either the positive strand (AF-PP-350: GTCACCGCTTGTAGAATTGTCATTGCCCTGTTGATGTTCCTTTCTGTTTGAACCTGGCTG & AF-PP-351: TCATCTATGGTTTGCCGATACGTGGTGTTGCTAATCCATGGCACT ACCATAGTACATGAG), or the negative strand (AF-PP-84: TTCACGGGTACGTTC ACTCCTGACAACAACCAGACATCACCTGCCCGCAGGTTCTGCCCG, AF-PP-86: ATTC GGACACCAAAACAAAGCGGTGTACACTGCAGGTTACAAAATTTGCAACTACCACTT) were end-labeled with 32P, and were successively used on the same blot (the blot was stripped in between the two hybridizations). The blots were hybridized and washed at 50°C, and exposed using a phosphorimager screen.
summary: Short RNAs derived from invading viruses with RNA genomes are important components of antiviral immunity in plants, worms and flies. The regulated generation of these short RNAs, and their engagement by the immune apparatus, is essential for inhibiting viral growth in these organisms. Mammals have the necessary protein components to generate these viral-derived short RNAs (“vsRNAs”), raising the question of whether vsRNAs in mammals are a general feature of infections with RNA viruses. Our work with Hepatitis C, Polio, Dengue, Vesicular Stomatitis, and West Nile viruses in a broad host repertoire demonstrates the generality of RNA virus-derived vsRNA production, and the ability of the cellular short RNA apparatus to engage these vsRNAs in mammalian cells. Detailed analyses of vsRNA and host-derived short RNA populations demonstrate both common and virus-specific features of the interplay between viral infection and short RNA populations. The vsRNA populations described in this work represent a novel dimension in both viral pathogenesis and host response.
keywords: immunology/immune response
  microbiology/innate immunity
  immunology/innate immunity
  molecular biology/bioinformatics
  infectious diseases/viral infections
  virology/effects of virus infection on host gene expression
  genetics and genomics/bioinformatics
  virology/host antiviral responses
title: Six RNA Viruses and Forty-One Hosts: Viral Small RNAs and Modulation of Small RNA Repertoires in Vertebrate and Invertebrate Systems