Intrahost Selection Pressures Drive Rapid Dengue Virus Microevolution in Acute Human Infections.

Dengue, caused by four dengue virus serotypes (DENV-1 to DENV-4), is a highly prevalent mosquito-borne viral disease in humans. Yet, selection pressures driving DENV microevolution within human hosts (intrahost) remain unknown. We employed a whole-genome segmented amplification approach coupled with deep sequencing to profile DENV-3 intrahost diversity in peripheral blood mononuclear cell (PBMC) and plasma samples from 77 dengue patients. DENV-3 intrahost diversity appears to be driven by immune pressures as well as replicative success in PBMCs and potentially other replication sites. Hotspots for intrahost variation were detected in 59%-78% of patients in the viral Envelope and pre-Membrane/Membrane proteins, which together form the virion surface. Dominant variants at the hotspots arose via convergent microevolution, appear to be immune-escape variants, and were evolutionarily constrained at the macro level due to viral replication defects. Dengue is thus an example of an acute infection in which selection pressures within infected individuals drive rapid intrahost virus microevolution.

Chikungunya Virus Sequences Across the First Epidemic in Nicaragua, 2014-2015.

Chikungunya is caused by the mosquito-borne arthrogenic alphavirus, chikungunya virus (CHIKV). Chikungunya was introduced into the Americas in late 2013 and Nicaragua in mid-2014. Here, we sequenced five imported and 30 autochthonous Nicaraguan CHIKV from cases identified in the first epidemic in the country between August 2014 and April 2015. One full-length and two partial genomic sequences were obtained by deep sequencing; Sanger methodology yielded 33 E1 sequences from five imported and 28 autochthonous cases. Phylogenetic analysis indicates that Nicaraguan CHIKV all belonged to the Asian genotype, Caribbean clade. Moreover, E1 gene sequences revealed accumulation of mutations in later months of the epidemic, including four silent mutations in 11 autochthonous cases and three non-synonymous mutations in three autochthonous cases. No mutations contributing to increased transmissibility by Aedes albopictus were identified in the E1 gene. This represents the most comprehensive set of CHIKV sequences available from the Americas to date.

Increased replicative fitness of a dengue virus 2 clade in native mosquitoes: potential contribution to a clade replacement event in Nicaragua.

UNLABELLED: The four dengue virus (DENV) serotypes (DENV serotype 1 [DENV-1] to DENV-4) are transmitted by Aedes aegypti and A. albopictus mosquitoes, causing up to 390 million DENV infections worldwide each year. We previously reported a clade replacement of the DENV-2 Asian-American genotype NI-1 clade by the NI-2B clade in Managua, Nicaragua. Here, we describe our studies of the replicative ability of NI-1 and NI-2B viruses in an A. aegypti cell line (Aag2) and A. aegypti mosquitoes reared from eggs collected in Managua. In coinfection experiments, several different pairs of NI-1 and NI-2B clinical isolates were used to infect Aag2 cells or blood-fed A. aegypti mosquitoes. Results consistently showed a significant replicative advantage of NI-2B over NI-1 viruses early after infection in vitro, and in mosquitoes, NI-2B viruses attained a higher replicative index than NI-1 isolates 3 to 7 days postinfection (dpi). At 7 dpi, NI-2B viruses displayed a significantly higher replicative index in legs and salivary glands; however, this advantage was lost by 14 and 21 dpi. We also found that the percentage of mosquitoes in which NI-2B viruses were dominant was significantly higher than that in which NI-1 viruses were dominant on day 7 but not at later time points. Taken together, these data demonstrate that clade NI-2B holds a replicative advantage over clade NI-1 early in infection that wanes at later time points. This early fitness advantage of NI-2B viruses over NI-1 viruses in the native vector, A. aegypti, suggests a shorter extrinsic incubation period for NI-2B viruses, which could have contributed to the clade replacement event in Managua. IMPORTANCE: Dengue virus (DENV), one of the most medically important arthropod-borne viruses, is transmitted to humans by Aedes aegypti and A. albopictus mosquitoes in tropical and subtropical regions worldwide. Dengue epidemics continue to increase in frequency, geographic range, and severity and are a major public health concern. This is due to globalization, unplanned urbanization, and climate change, as well as host genetics and immune responses and viral genetic changes. DENV consists of four serotypes, in turn composed of genotypes and genetically distinct clades. What drives the frequent replacement of a previously circulating DENV clade by another is unclear. Here, we investigate the replicative fitness of two clades of DENV serotype 2 in Aedes aegypti cells and mosquitoes collected from the region where the viruses circulated and conclude that increased replicative fitness could have contributed to a DENV clade replacement event in Nicaragua. These findings provide insight into vector-driven evolution of DENV epidemics.

Convergent antibody signatures in human dengue.

Dengue is the most prevalent mosquito-borne viral disease in humans, and the lack of early prognostics, vaccines, and therapeutics contributes to immense disease burden. To identify patterns that could be used for sequence-based monitoring of the antibody response to dengue, we examined antibody heavy-chain gene rearrangements in longitudinal peripheral blood samples from 60 dengue patients. Comparing signatures between acute dengue, postrecovery, and healthy samples, we found increased expansion of B cell clones in acute dengue patients, with higher overall clonality in secondary infection. Additionally, we observed consistent antibody sequence features in acute dengue in the highly variable major antigen-binding determinant, complementarity-determining region 3 (CDR3), with specific CDR3 sequences highly enriched in acute samples compared to postrecovery, healthy, or non-dengue samples. Dengue thus provides a striking example of a human viral infection where convergent immune signatures can be identified in multiple individuals. Such signatures could facilitate surveillance of immunological memory in communities.

Characterization of a model of lethal dengue virus 2 infection in C57BL/6 mice deficient in the alpha/beta interferon receptor.

Dengue virus (DENV) causes dengue fever and dengue haemorrhagic fever/dengue shock syndrome, both considered major public-health problems worldwide. We generated a lethal DENV-2 strain (D220) by 10 additional cycles of subcutaneous inoculation of mice with supernatant from mosquito cells infected with the previously characterized strain D2S10, followed by harvesting of serum. D220 induces mortality at ten-fold lower doses than D2S10 in mice lacking only the alpha/beta interferon (IFN-α/ß) receptor in C57BL/6 or 129 backgrounds under both non-enhanced and antibody-enhanced conditions. Sequence analysis of the complete viral genome revealed five amino acid changes between D220 and D2S10, of which two (K122I in envelope and V115A in NS4B) appear to account for the observed phenotypic differences between the viruses. By causing mortality at lower doses in C57BL/6 mice lacking only the IFN-α/ß receptor, D220 constitutes an improved tool for study of DENV-induced pathogenesis, as well as for testing potential vaccines and antiviral drugs against DENV.

Genome-wide patterns of intrahuman dengue virus diversity reveal associations with viral phylogenetic clade and interhost diversity.

Analogous to observations in RNA viruses such as human immunodeficiency virus, genetic variation associated with intrahost dengue virus (DENV) populations has been postulated to influence viral fitness and disease pathogenesis. Previous attempts to investigate intrahost genetic variation in DENV characterized only a few viral genes or a limited number of full-length genomes. We developed a whole-genome amplification approach coupled with deep sequencing to capture intrahost diversity across the entire coding region of DENV-2. Using this approach, we sequenced DENV-2 genomes from the serum of 22 Nicaraguan individuals with secondary DENV infection and captured ∼75% of the DENV genome in each sample (range, 40 to 98%). We identified and quantified variants using a highly sensitive and specific method and determined that the extent of diversity was considerably lower than previous estimates. Significant differences in intrahost diversity were detected between genes and also between antigenically distinct domains of the Envelope gene. Interestingly, a strong association was discerned between the extent of intrahost diversity in a few genes and viral clade identity. Additionally, the abundance of viral variants within a host, as well as the impact of viral mutations on amino acid encoding and predicted protein function, determined whether intrahost variants were observed at the interhost level in circulating Nicaraguan DENV-2 populations, strongly suggestive of purifying selection across transmission events. Our data illustrate the value of high-coverage genome-wide analysis of intrahost diversity for high-resolution mapping of the relationship between intrahost diversity and clinical, epidemiological, and virological parameters of viral infection.

On the nature of in vivo requirements for rde-4 in RNAi and developmental pathways in C. elegans.

C. elegans RDE-4 is a double-stranded RNA binding protein that has been shown to play a key role in response to foreign double-stranded RNA (dsRNA). We have used diverse tools for analysis of gene function to characterize the domain and organismal foci of RDE-4 action in C. elegans. First, we examined the focus of activity within the RDE-4 protein, by testing a series of RDE-4 deletion constructs for their ability to support dsRNA-triggered gene silencing. These assays indicated a molecular requirement for a linker region and the second dsRNA-binding domain of RDE-4, with ancillary contributions to function from the C and N terminal domains. Second, we used mosaic analysis to explore the cellular focus of action of RDE-4. These experiments indicated an ability of RDE-4 to function non-autonomously in foreign RNA responses. Third, we used growth under stressful conditions to search for evidence of an organismal focus of action for RDE-4 distinct from its role in response to foreign dsRNA. Propagation at high temperatures exposed a conditional requirement for RDE-4 for optimal growth and fertility, indicating at least under these conditions that RDE-4 can serve an essential role in C. elegans.

Distinct phases of siRNA synthesis in an endogenous RNAi pathway in C. elegans soma.

Endogenous RNA-directed RNA polymerases (RdRPs) are cellular components capable of synthesizing new complementary RNAs from existing RNA templates. We present evidence for successive engagement of two different RdRPs in an endogenous siRNA-based mechanism targeting specific mRNAs in C. elegans soma. In the initiation stage of this process, a group of mRNA species are chosen as targets for downregulation, leading to accumulation of rare 26 nt 5'-phosphorylated antisense RNAs that depend on the RdRP homolog RRF-3, the Argonaute ERGO-1, DICER, and a series of associated ("ERI") factors. This primary process leads to production of a much more abundant class of 22 nt antisense RNAs, dependent on a secondary RdRP (RRF-1) and associating with at least one distinct Argonaute (NRDE-3). The requirement for two RdRP/Argonaute combinations and initiation by a rare class of uniquely structured siRNAs in this pathway illustrate the caution and flexibility used as biological systems exploit the physiological copying of RNA.

Six RNA viruses and forty-one hosts: viral small RNAs and modulation of small RNA repertoires in vertebrate and invertebrate systems.

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.

Human tRNA-derived small RNAs in the global regulation of RNA silencing.

Competition between mammalian RNAi-related gene silencing pathways is well documented. It is therefore important to identify all classes of small RNAs to determine their relationship with RNAi and how they affect each other functionally. Here, we identify two types of 5'-phosphate, 3'-hydroxylated human tRNA-derived small RNAs (tsRNAs). tsRNAs differ from microRNAs in being essentially restricted to the cytoplasm and in associating with Argonaute proteins, but not MOV10. The first type belongs to a previously predicted Dicer-dependent class of small RNAs that we find can modestly down-regulate target genes in trans. The 5' end of type II tsRNA was generated by RNaseZ cleavage downstream from a tRNA gene, while the 3' end resulted from transcription termination by RNA polymerase III. Consistent with their preferential association with the nonslicing Argonautes 3 and 4, canonical gene silencing activity was not observed for type II tsRNAs. The addition, however, of an oligonucleotide that was sense to the reporter gene, but antisense to an overexpressed version of the type II tsRNA, triggered robust, >80% gene silencing. This correlated with the redirection of the thus reconstituted fully duplexed double-stranded RNA into Argonaute 2, whereas Argonautes 3 and 4 were skewed toward less structured small RNAs, particularly single-strand RNAs. We observed that the modulation of tsRNA levels had minor effects on the abundance of microRNAs, but more pronounced changes in the silencing activities of both microRNAs and siRNAs. These findings support that tsRNAs are involved in the global control of small RNA silencing through differential Argonaute association, suggesting that small RNA-mediated gene regulation may be even more finely regulated than previously realized.

Capped small RNAs and MOV10 in human hepatitis delta virus replication.

The evolutionary origin of human hepatitis delta virus (HDV) replication by RNA-directed transcription is unclear. Here we identify two species of 5'-capped, approximately 18-25-nucleotide small RNAs. One was of antigenomic polarity, corresponding to the 5' end of hepatitis delta antigen (HDAg) mRNA, and interacted with HDAg and RNA polymerase II (Pol II), whereas the other mapped to a structurally analogous region on the genomic RNA hairpin. An HDAg-interaction screen indicated that HDAg interacts with MOV10, the human homolog of the Arabidopsis thaliana RNA amplification factor gene SDE3 and Drosophila melanogaster RISC-maturation factor gene Armitage (armi). MOV10 knockdown inhibited HDV replication, but not HDAg mRNA translation, further supporting a role for MOV10 in RNA-directed transcription. Together, our studies define RNA hairpins as critical elements for the initiation of HDV-related, RNA-directed transcription. The identification of capped small RNAs and the involvement of MOV10 in HDV replication further suggest a conserved mechanism related to RNA-directed transcription in lower eukaryotes.

A pyrosequencing-tailored nucleotide barcode design unveils opportunities for large-scale sample multiplexing.

Multiplexed high-throughput pyrosequencing is currently limited in complexity (number of samples sequenced in parallel), and in capacity (number of sequences obtained per sample). Physical-space segregation of the sequencing platform into a fixed number of channels allows limited multiplexing, but obscures available sequencing space. To overcome these limitations, we have devised a novel barcoding approach to allow for pooling and sequencing of DNA from independent samples, and to facilitate subsequent segregation of sequencing capacity. Forty-eight forward-reverse barcode pairs are described: each forward and each reverse barcode unique with respect to at least 4 nt positions. With improved read lengths of pyrosequencers, combinations of forward and reverse barcodes may be used to sequence from as many as n(2) independent libraries for each set of 'n' forward and 'n' reverse barcodes, for each defined set of cloning-linkers. In two pilot series of barcoded sequencing using the GS20 Sequencer (454/Roche), we found that over 99.8% of obtained sequences could be assigned to 25 independent, uniquely barcoded libraries based on the presence of either a perfect forward or a perfect reverse barcode. The false-discovery rate, as measured by the percentage of sequences with unexpected perfect pairings of unmatched forward and reverse barcodes, was estimated to be <0.005%.