Experimental evolution of an RNA virus in wild birds: evidence for host-dependent impacts on population structure and competitive fitness

Experimental evolution of an RNA virus in Caenorhabditis elegans

The discovery of Orsay virus (OrV), the first virus infecting wild populations of Caenorhabditis elegans, has boosted studies of viral immunity pathways in this nematode. Considering the many advantages that C. elegans offers for fundamental research in host-pathogen interactions, this pathosystem has high potential to become a model system for experimental virus evolution studies. However, the evolutionary constraints - i.e, the balance between genetic variation, selection, drift and historical contingency- operating in this pathosystem have barely been explored. Here we describe for the first time an evolution experiment of two different OrV strains in C. elegans. Comparison of the two ancestral strains showed differences in infectivity and sequence, and highlighted the importance of consistently normalize viral inocula for meaningful comparisons among strains. After 10 serial passages of evolution, we report slight changes in infectivity and non-synonymous mutations fixed in the evolved viral populations. In addition, we observed numerous minor variants emerging in the viral population. These minor variants were not randomly distributed along the genome but concentrated in polymorphic genomic regions. Overall, our work established the grounds for future experimental virus evolution studies using Caenorhabditis nematodes.

Maintenance of a host-specific minority mutation in the West Nile virus NS3

West Nile virus (WNV), the most prevalent arthropod-borne virus (arbovirus) in the United States, is maintained in a cycle between Culex spp. mosquitoes and birds. Arboviruses exist within hosts and vectors as a diverse set of closely related genotypes. In theory, this genetic diversity can facilitate adaptation to distinct environments during host cycling, yet host-specific fitness of minority genotypes has not been assessed. Utilizing WNV deep-sequencing data, we previously identified a naturally occurring, mosquito-biased substitution, NS3 P319L. Using both cell culture and experimental infection in natural hosts, we demonstrated that this substitution confers attenuation in vertebrate hosts and increased transmissibility by mosquitoes. Biochemical assays demonstrated temperature-sensitive ATPase activity consistent with host-specific phenotypes. Together these data confirm the maintenance of host-specific minority variants in arbovirus mutant swarms, suggest a unique role for NS3 in viral fitness, and demonstrate that intrahost sequence data can inform mechanisms of host-specific adaptation.

Intracellular Diversity of WNV within Circulating Avian Peripheral Blood Mononuclear Cells Reveals Host-Dependent Patterns of Polyinfection

Arthropod-borne virus (arbovirus) populations exist as mutant swarms that are maintained between arthropods and vertebrates. West Nile virus (WNV) population dynamics are host-dependent. In American crows, purifying selection is weak and population diversity is high compared to American robins, which have 100- to 1000-fold lower viremia. WNV passed in robins leads to fitness gains, whereas that passed in crows does not. Therefore, we tested the hypothesis that high crow viremia allows for higher genetic diversity within individual avian peripheral blood mononuclear cells (PBMCs), reasoning that this could have produced the previously observed host-specific differences in genetic diversity and fitness. Specifically, we infected cells and birds with a molecularly barcoded WNV and sequenced viral RNA from single cells to quantify the number of WNV barcodes in each. Our results demonstrate that the richness of WNV populations within crows far exceeds that in robins. Similarly, rare WNV variants were maintained by crows more frequently than by robins. Our results suggest that increased viremia in crows relative to robins leads to the maintenance of defective genomes and less prevalent variants, presumably through complementation. Our findings further suggest that weaker purifying selection in highly susceptible crows is attributable to this higher viremia, polyinfections and complementation.

Viral intra-host evolutionary dynamics revealed via serial passage of Japanese encephalitis virus in vitro

Analyses of viral inter- and intra-host mutations could better guide the prevention and control of infectious diseases. For a long time, studies on viral evolution have focused on viral inter-host variations. Next-generation sequencing has accelerated the investigations of viral intra-host diversity. However, the theoretical basis and dynamic characteristics of viral intra-host mutations remain unknown. Here, using serial passages of the SA14-14-2 vaccine strain of Japanese encephalitis virus (JEV) as the in vitro model, the distribution characteristics of 1,788 detected intra-host single-nucleotide variations (iSNVs) and their mutated frequencies from 477 deep-sequenced samples were analyzed. Our results revealed that in adaptive (baby hamster kidney (BHK)) cells, JEV is under a nearly neutral selection pressure, and both non-synonymous and synonymous mutations represent an S-shaped growth trend over time. A higher positive selection pressure was observed in the nonadaptive (C6/36) cells, and logarithmic growth in non-synonymous iSNVs and linear growth in synonymous iSNVs were observed over time. Moreover, the mutation rates of the NS4B protein and the untranslated region (UTR) of the JEV are significantly different between BHK and C6/36 cells, suggesting that viral selection pressure is regulated by different cellular environments. In addition, no significant difference was detected in the distribution of mutated frequencies of iSNVs between BHK and C6/36 cells.

Intracellular diversity of WNV within circulating avian peripheral blood mononuclear cells reveals host-dependent patterns of polyinfection

Error-prone replication of RNA viruses generates the genetic diversity required for adaptation within rapidly changing environments. Thus, arthropod-borne virus (arbovirus) populations exist in nature as mutant swarms that are maintained between arthropods and vertebrates. Previous studies have demonstrated that West Nile virus (WNV) population dynamics are host dependent: In American crows, which experience extremely high viremia, purifying selection is weak and population diversity is high compared to American robins, which have 100 to 1000-fold lower viremia. WNV passed in robins experiences fitness gains, whereas that passed in crows does not. Therefore, we tested the hypothesis that high crow viremia allows higher genetic diversity within individual avian peripheral-blood mononuclear cells (PBMCs), reasoning that this could have produced the previously observed host-specific differences in genetic diversity and fitness. Specifically, we infected cells and birds with a novel, barcoded version of WNV and sequenced viral RNA from single cells to quantify the number of WNV barcodes that each contained. Our results demonstrate that the richness of WNV populations within crows far exceeds that in robins. Similarly, rare WNV variants were maintained by crows more frequently than by robins. Our results suggest that increased viremia in crows relative to robins leads to maintenance of defective genomes and less prevalent variants, presumably through complementation. Our findings further suggest that weaker purifying selection in highly susceptible crows is attributable to this higher viremia, polyinfections and complementation. These studies further document the role of particular, ecologically relevant hosts in shaping virus population structure.

Author Summary

WNV mutational diversity in vertebrates is species-dependent. In crows, low frequency variants are common, and viral populations are more diverse. In robins, fewer mutations become permanent fixtures of the overall viral population. We infected crows, robins and a chicken cell line with a genetically marked (barcoded) WNV. Higher levels of virus led to multiple unique WNV genomes infecting individual cells, even when a genotype was present at low levels in the input viral stock. Our findings suggest that higher levels of circulating virus in natural hosts allow less fit viruses to survive in RNA virus populations through complementation by more fit viruses. This is significant as it allows less represented and less fit viruses to be maintained at low levels until they potentially emerge when virus environments change. Overall our data reveal new insights on the relationships between host susceptibility to high viremia and virus evolution.

The Role of the Flavivirus Replicase in Viral Diversity and Adaptation

Flaviviruses include several emerging and re-emerging arboviruses which cause millions of infections each year. Although relatively well-studied, much remains unknown regarding the mechanisms and means by which these viruses readily alternate and adapt to different hosts and environments. Here, we review a subset of the different aspects of flaviviral biology which impact host switching and viral fitness. These include the mechanism of replication and structural biology of the NS3 and NS5 proteins, which reproduce the viral genome; rates of mutation resulting from this replication and the role of mutational frequency in viral fitness; and the theory of quasispecies evolution and how it contributes to our understanding of genetic and phenotypic plasticity.

Translating virus evolution into epidemiology

6 years ago, I published a Cell Host and Microbe paper that described randomization of virus genetic populations during mosquito infection. From within the evolutionary chaos, however, there is an order that can reveal a virus' past. Using these insights, I forged a career harnessing virus evolution to understand epidemiological patterns.

American alligators are capable of West Nile virus amplification, mosquito infection and transmission

West Nile virus (WNV) overwintering is poorly understood and likely multifactorial. Interest in alligators as a potential amplifying host arose when it was shown that they develop viremias theoretically sufficient to infect mosquitoes. We examined potential ways in which alligators may contribute to the natural ecology of WNV. We experimentally demonstrated that alligators are capable of WNV amplification with subsequent mosquito infection and transmission capability, that WNV-infected mosquitoes readily infect alligators and that water can serve as a source of infection for alligators but does not easily serve as in intermediate means for transmission between birds and alligators. These findings indicate potential mechanisms for maintenance of WNV outside of the primary bird-mosquito transmission cycle.

Impact of extrinsic incubation temperature on natural selection during Zika virus infection of Aedes aegypti and Aedes albopictus

Arthropod-borne viruses (arboviruses) require replication across a wide range of temperatures to perpetuate. While vertebrate hosts tend to maintain temperatures of approximately 37°C—40°C, arthropods are subject to ambient temperatures which can have a daily fluctuation of > 10°C. Temperatures impact vector competence, extrinsic incubation period, and mosquito survival unimodally, with optimal conditions occurring at some intermediate temperature. In addition, the mean and range of daily temperature fluctuations influence arbovirus perpetuation and vector competence. The impact of temperature on arbovirus genetic diversity during systemic mosquito infection, however, is poorly understood. Therefore, we determined how constant extrinsic incubation temperatures of 25°C, 28°C, 32°C, and 35°C control Zika virus (ZIKV) vector competence and population dynamics within Aedes aegypti and Aedes albopictus mosquitoes. We also examined fluctuating temperatures which better mimic field conditions in the tropics. We found that vector competence varied in a unimodal manner for constant temperatures peaking between 28°C and 32°C for both Aedes species. Transmission peaked at 10 days post-infection for Aedes aegypti and 14 days for Aedes albopictus. Conversely, fluctuating temperature decreased vector competence. Using RNA-seq to characterize ZIKV population structure, we identified that temperature alters the selective environment in unexpected ways. During mosquito infection, constant temperatures more often elicited positive selection whereas fluctuating temperatures led to strong purifying selection in both Aedes species. These findings demonstrate that temperature has multiple impacts on ZIKV biology, including major effects on the selective environment within mosquitoes.

Arthropod-borne viruses (arboviruses) have emerged in recent decades due to complex factors that include increases in international travel and trade, the breakdown of public health infrastructure, land use changes, and many others. Climate change also has the potential to shift the geographical ranges of arthropod vectors, consequently increasing the global risk of arbovirus infection. Changing temperatures may alter the virus-host interaction, ultimately resulting in the emergence of new viruses and virus genotypes in new areas. Therefore, we sought to characterize how temperature (both constant and fluctuating) alters the ability of Aedes aegypti and Aedes albopictus to transmit Zika virus, and how it influences virus populations within mosquitoes. We found that intermediate temperatures maximize virus transmission compared to more extreme and fluctuating temperatures. Constant temperatures increased positive selection on virus genomes, while fluctuating temperatures strengthened purifying selection. Our studies provide evidence that in addition to altering vector competence, temperature significantly influences natural selection within mosquitoes.

Infection, dissemination, and transmission efficiencies of Zika virus in Aedes aegypti after serial passage in mosquito or mammalian cell lines or alternating passage in both cell types

Background

Zika virus (ZIKV) is an arthropod-borne virus (arbovirus) with an urban transmission cycle that primarily involves humans and Aedes aegypti. Evidence suggests that the evolution of some arboviruses is constrained by their dependency on alternating between disparate (vertebrate and invertebrate) hosts. The goals of this study are to compare the genetic changes that occur in ZIKV after serial passaging in mosquito or vertebrate cell lines or alternate passaging in both cell types and to compare the replication, dissemination, and transmission efficiencies of the cell culture-derived viruses in Ae. aegypti.

Methods

An isolate of ZIKV originally acquired from a febrile patient in Yucatan, Mexico, was serially passaged six times in African green monkey kidney (Vero) cells or Aedes albopictus (C6/36) cells or both cell types by alternating passage. A colony of Ae. aegypti from Yucatan was established, and mosquitoes were challenged with the cell-adapted viruses. Midguts, Malpighian tubules, ovaries, salivary glands, wings/legs and saliva were collected at various times after challenge and tested for evidence of virus infection.

Results

Genome sequencing revealed the presence of two non-synonymous substitutions in the premembrane and NS1 regions of the mosquito cell-adapted virus and two non-synonymous substitutions in the capsid and NS2A regions of both the vertebrate cell-adapted and alternate-passaged viruses. Additional genetic changes were identified by intrahost variant frequency analysis. Virus maintained by continuous C6/36 cell passage was significantly more infectious in Ae. aegypti than viruses maintained by alternating passage and consecutive Vero cell passage.

Conclusions

Mosquito cell-adapted ZIKV displayed greater in vivo fitness in Ae. aegypti compared to the other viruses, indicating that obligate cycling between disparate hosts carries a fitness cost. These data increase our understanding of the factors that drive ZIKV adaptation and evolution and underscore the important need to consider the in vivo passage histories of flaviviruses to be evaluated in vector competence studies.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04726-1.

Intra-host evolutionary dynamics of the hepatitis C virus among people who inject drugs

Most individuals chronically infected with hepatitis C virus (HCV) are asymptomatic during the initial stages of infection and therefore the precise timing of infection is often unknown. Retrospective estimation of infection duration would improve existing surveillance data and help guide treatment. While intra-host viral diversity quantifications such as Shannon entropy have previously been utilized for estimating duration of infection, these studies characterize the viral population from only a relatively short segment of the HCV genome. In this study intra-host diversities were examined across the HCV genome in order to identify the region most reflective of time and the degree to which these estimates are influenced by high-risk activities including those associated with HCV acquisition. Shannon diversities were calculated for all regions of HCV from 78 longitudinally sampled individuals with known seroconversion timeframes. While the region of the HCV genome most accurately reflecting time resided within the NS3 gene, the gene region with the highest capacity to differentiate acute from chronic infections was identified within the NS5b region. Multivariate models predicting duration of infection from viral diversity significantly improved upon incorporation of variables associated with recent public, unsupervised drug use. These results could assist the development of strategic population treatment guidelines for high-risk individuals infected with HCV and offer insights into variables associated with a likelihood of transmission.

Alternation between taxonomically divergent hosts is not the major determinant of flavivirus evolution

Flaviviruses display diverse epidemiological and ecological features. Tick-borne and mosquito-borne flaviviruses (TBFV and MBFV, respectively) are important human pathogens that alternate replication in invertebrate vectors and vertebrate hosts. The Flavivirus genus also includes insect-specific viruses (ISFVs) and viruses with unknown invertebrate hosts. It is generally accepted that viruses that alternate between taxonomically different hosts evolve slowly and that the evolution of MBFVs and TBFVs is dominated by strong constraints, with limited episodes of positive selection. We exploited the availability of flavivirus genomes to test these hypotheses and to compare their rates and patterns of evolution. We estimated the substitution rates of CFAV and CxFV (two ISFVs) and, by taking into account the time-frame of measurement, compared them with those of other flaviviruses. Results indicated that CFAV and CxFV display relatively different substitution rates. However, these data, together with estimates for single-host members of the Flaviviridae family, indicated that MBFVs do not display relatively slower evolution. Conversely, TBFVs displayed some of lowest substitution rates among flaviviruses. Analysis of selective patterns over longer evolutionary time-frames confirmed that MBFVs evolve under strong purifying selection. Interestingly, TBFVs and ISFVs did not show extremely different levels of constraint, although TBFVs alternate among hosts, whereas ISFVs do not. Additional results showed that episodic positive selection drove the evolution of MBFVs, despite their high constraint. Positive selection was also detected on two branches of the TBFVs phylogeny that define the seabird clade. Thus, positive selection was much more common during the evolution of arthropod-borne flaviviruses than previously thought. Overall, our data indicate that flavivirus evolutionary patterns are complex and most likely determined by multiple factors, not limited to the alternation between taxonomically divergent hosts. The frequency of both positive and purifying selection, especially in MBFVs, suggests that a minority of sites in the viral polyprotein experience weak constraint and can evolve to generate new viral phenotypes and possibly promote adaptation to new hosts.

Genome Number and Size Polymorphism in Zika Virus Infectious Units

The arthropod-borne Zika virus (ZIKV) infects humans and can cause severe neurological sequelae, particularly in fetuses infected in utero. How this virus has been able to spread across vast geological ranges and evolve in new host populations is not yet understood.

Zika virus (ZIKV; Flaviviridae, Flavivirus) is an arthropod-borne infection that can result in severe outcomes, particularly in fetuses infected in utero. It has been assumed that infection by ZIKV, as well as other viruses, is largely initiated by individual virus particles binding to and entering a cell. However, recent studies have demonstrated that multiple virus particles are frequently delivered to a cell simultaneously and that this collective particle delivery enhances infection. ZIKV is maintained in nature between Aedes aegypti mosquitos and vertebrate hosts, including humans. Human infection is initiated through the injection of a relatively small initial inoculum comprised of a genetically complex virus population. Since most mutations decrease virus fitness, collective particle transmission could benefit ZIKV and other arthropod-borne diseases by facilitating the maintenance of genetic complexity and adaptability during infection or through other mechanisms. Therefore, we utilized a barcoded ZIKV to quantify the number of virus genomes that initiate a plaque. We found that individual plaques contain a mean of 10 infecting viral genomes (range, 1 to 212). Few plaques contained more than two dominant genomes. To determine whether multigenome infectious units consist of collectively transmitting virions, infectious units of ZIKV were then separated mechanically by centrifugation, and heavier fractions were found to contain more genomes per plaque-forming unit, with larger diameters. Finally, larger/heavier infectious units reformed after removal. These data suggest that ZIKV populations consist of a variety of infectious unit sizes, likely mostly made up of aggregates, and only rarely begin with a single virus genome.

IMPORTANCE The arthropod-borne Zika virus (ZIKV) infects humans and can cause severe neurological sequelae, particularly in fetuses infected in utero. How this virus has been able to spread across vast geological ranges and evolve in new host populations is not yet understood. This research demonstrates a novel mechanism of ZIKV transmission through multigenome aggregates, providing insight into ZIKV evolution, immunologic evasion, and better future therapeutic design. This study shows that ZIKV plaques result from collections of genomes rather than individual genomes, increasing the potential for interactions between ZIKV genotypes.

Recurrent evolution of high virulence in isolated populations of a DNA virus

Hosts and viruses are constantly evolving in response to each other: as a host attempts to suppress a virus, the virus attempts to evade and suppress the host's immune system. Here, we describe the recurrent evolution of a virulent strain of a DNA virus, which infects multiple Drosophila species. Specifically, we identified two distinct viral types that differ 100-fold in viral titer in infected individuals, with similar differences observed in multiple species. Our analysis suggests that one of the viral types recurrently evolved at least four times in the past ~30,000 years, three times in Arizona and once in another geographically distinct species. This recurrent evolution may be facilitated by an effective mutation rate which increases as each prior mutation increases viral titer and effective population size. The higher titer viral type suppresses the host-immune system and an increased virulence compared to the low viral titer type.

Evolution of resistance to fluoroquinolones by dengue virus serotype 4 provides insight into mechanism of action and consequences for viral fitness

Drugs against flaviviruses such as dengue (DENV) and Zika (ZIKV) virus are urgently needed. We previously demonstrated that three fluoroquinolones, ciprofloxacin, enoxacin, and difloxacin, suppress replication of six flaviviruses. To investigate the barrier to resistance and mechanism(s) of action of these drugs, DENV-4 was passaged in triplicate in HEK-293 cells in the presence or absence of each drug. Resistance to ciprofloxacin was detected by the seventh passage and to difloxacin by the tenth, whereas resistance to enoxacin did not occur within ten passages. Two putative resistance-conferring mutations were detected in the envelope gene of ciprofloxacin and difloxacin-resistant DENV-4. In the absence of ciprofloxacin, ciprofloxacin-resistant viruses sustained a significantly higher viral titer than control viruses in HEK-293 and HuH-7 cells and resistant viruses were more stable than control viruses at 37 °C. These results suggest that the mechanism of action of ciprofloxacin and difloxacin involves interference with virus binding or entry.

Species-Specific Evolution of Ebola Virus during Replication in Human and Bat Cells

Ebola virus (EBOV) causes a severe, often fatal disease in humans and nonhuman primates. Within the past decade, EBOV has caused two large and difficult-to-control outbreaks, one of which recently ended in the Democratic Republic of the Congo. Bats are the likely reservoir of EBOV, but little is known of their relationship with the virus. We perform serial passages of EBOV in human and bat cells and use circular sequencing to compare the short-term evolution of the virus. Virus populations passaged in bat cells have sequence markers indicative of host RNA editing enzyme activity, including evidence for ADAR editing of the EBOV glycoprotein. Multiple regions in the EBOV genome appear to have undergone adaptive evolution when passaged in bat and human cells. Individual mutated viruses are rescued and characterized. Our results provide insight into the host species-specific evolution of EBOV and highlight the adaptive flexibility of the virus.

West Nile Virus fidelity modulates the capacity for host cycling and adaptation

The fidelity of flaviviruses is thought to be tightly regulated for optimal fitness within and between hosts. West Nile virus (WNV) high-fidelity (HiFi) mutations V793I and G806R within the RNA-dependent RNA polymerase, and low-fidelity (LoFi) mutation T248I within the methyltransferase, were previously shown to attenuate infectivity and replicative fitness in Culex mosquitoes and Culex tarsalis (CXT) cells but not in mammalian cells. We hypothesized that fidelity alterations would modify adaptation and maintenance in a host-specific manner. To test this hypothesis, wild-type (WT), HiFi (V793I/G806R) and LoFi (T248I) variants were sequentially passaged eight times in avian (PDE) or mosquito cells, or alternately between the two. Initial characterization confirmed that fidelity mutants are attenuated in mosquito, but not avian, cells. Deep sequencing revealed mutations unique to both cell lines and fidelity mutants, including ENV G1378A, a mutation associated with avian cell adaptation. To characterize maintenance and adaptation, viral outputs were monitored throughout passaging and viral fitness was assessed. The results indicate that fidelity mutants can at times recover fitness during mosquito cell passage, but remain attenuated relative to WT. Despite similar initial fitness, LoFi mutants were impaired during sequential passage in avian cells. Conversely, HiFi mutants passaged in avian cells showed increased adaptation, suggesting that increased fidelity may be advantageous in avian hosts. Although some adaptation occurred with individual mutants, the output titres of fidelity mutants were on average lower and were often lost during host switching. These data confirm that arbovirus fidelity is likely fine-tuned to maximize survival in disparate hosts.

An RNA Thermometer Activity of the West Nile Virus Genomic 3'-Terminal Stem-Loop Element Modulates Viral Replication Efficiency during Host Switching

The 3′-terminal stem-loop (3′SL) of the RNA genome of the flavivirus West Nile (WNV) harbors, in its stem, one of the sequence elements that are required for genome cyclization. As cyclization is a prerequisite for the initiation of viral replication, the 3′SL was proposed to act as a replication silencer. The lower part of the 3′SL is metastable and confers a structural flexibility that may regulate the switch from the linear to the circular conformation of the viral RNA. In the human system, we previously demonstrated that a cellular RNA-binding protein, AUF1 p45, destabilizes the 3′SL, exposes the cyclization sequence, and thus promotes flaviviral genome cyclization and RNA replication. By investigating mutant RNAs with increased 3′SL stabilities, we showed the specific conformation of the metastable element to be a critical determinant of the helix-destabilizing RNA chaperone activity of AUF1 p45 and of the precision and efficiency of the AUF1 p45-supported initiation of RNA replication. Studies of stability-increasing mutant WNV replicons in human and mosquito cells revealed that the cultivation temperature considerably affected the replication efficiencies of the viral RNA variants and demonstrated the silencing effect of the 3′SL to be temperature dependent. Furthermore, we identified and characterized mosquito proteins displaying similar activities as AUF1 p45. However, as the RNA remodeling activities of the mosquito proteins were found to be considerably lower than those of the human protein, a potential cell protein-mediated destabilization of the 3′SL was suggested to be less efficient in mosquito cells. In summary, our data support a model in which the 3′SL acts as an RNA thermometer that modulates flavivirus replication during host switching.

Comparative Pathology of West Nile Virus in Humans and Non-Human Animals

West Nile virus (WNV) continues to be a major cause of human arboviral neuroinvasive disease. Susceptible non-human vertebrates are particularly diverse, ranging from commonly affected birds and horses to less commonly affected species such as alligators. This review summarizes the pathology caused by West Nile virus during natural infections of humans and non-human animals. While the most well-known findings in human infection involve the central nervous system, WNV can also cause significant lesions in the heart, kidneys and eyes. Time has also revealed chronic neurologic sequelae related to prior human WNV infection. Similarly, neurologic disease is a prominent manifestation of WNV infection in most non-human non-host animals. However, in some avian species, which serve as the vertebrate host for WNV maintenance in nature, severe systemic disease can occur, with neurologic, cardiac, intestinal and renal injury leading to death. The pathology seen in experimental animal models of West Nile virus infection and knowledge gains on viral pathogenesis derived from these animal models are also briefly discussed. A gap in the current literature exists regarding the relationship between the neurotropic nature of WNV in vertebrates, virus propagation and transmission in nature. This and other knowledge gaps, and future directions for research into WNV pathology, are addressed.

In vitro evolution of herpes simplex virus 1 (HSV-1) reveals selection for syncytia and other minor variants in cell culture

The large dsDNA virus herpes simplex virus 1 (HSV-1) is considered to be genetically stable, yet it can rapidly evolve in response to strong selective pressures such as antiviral treatment. Deep sequencing has revealed that clinical and laboratory isolates of this virus exist as populations that contain a mixture of minor alleles or variants, similar to many RNA viruses. The classic virology approach of plaque purifying virus creates a genetically homogenous population, but it is not clear how closely this represents the mixed virus populations found in nature. We sought to study the evolution of mixed versus highly purified HSV-1 populations in controlled cell culture conditions, to examine the impact of this genetic diversity on evolution. We found that a mixed population of HSV-1 acquired more genetic diversity and underwent a more dramatic phenotypic shift than a plaque-purified population, producing a viral population that was almost entirely syncytial after just ten passages. At the genomic level, adaptation and genetic diversification occurred at the level of minor alleles or variants in the viral population. Certain genetic variants in the mixed viral population appeared to be positively selected in cell culture, and this shift was also observed in clinical samples during their first passages in vitro. In contrast, the plaque-purified viral population did not appear to change substantially in phenotype or overall quantity of minor allele diversity. These data indicate that HSV-1 is capable of evolving rapidly in a given environment, and that this evolution is facilitated by diversity in the viral population.

On the Fly: Interactions Between Birds, Mosquitoes, and Environment That Have Molded West Nile Virus Genomic Structure Over Two Decades

West Nile virus (WNV) was first identified in North America almost 20 yr ago. In that time, WNV has crossed the continent and established enzootic transmission cycles, resulting in intermittent outbreaks of human disease that have largely been linked with climatic variables and waning avian seroprevalence. During the transcontinental dissemination of WNV, the original genotype has been displaced by two principal extant genotypes which contain an envelope mutation that has been associated with enhanced vector competence by Culex pipiens L. (Diptera: Culicidae) and Culex tarsalis Coquillett vectors. Analyses of retrospective avian host competence data generated using the founding NY99 genotype strain have demonstrated a steady reduction in viremias of house sparrows over time. Reciprocally, the current genotype strains WN02 and SW03 have demonstrated an inverse correlation between house sparrow viremia magnitude and the time since isolation. These data collectively indicate that WNV has evolved for increased avian viremia while house sparrows have evolved resistance to the virus such that the relative host competence has remained constant. Intrahost analyses of WNV evolution demonstrate that selection pressures are avian species-specific and purifying selection is greater in individual birds compared with individual mosquitoes, suggesting that the avian adaptive and/or innate immune response may impose a selection pressure on WNV. Phylogenomic, experimental evolutionary systems, and models that link viral evolution with climate, host, and vector competence studies will be needed to identify the relative effect of different selective and stochastic mechanisms on viral phenotypes and the capacity of newly evolved WNV genotypes for transmission in continuously changing landscapes.

Twenty years of West Nile virus spread and evolution in the Americas visualized by Nextstrain

It has been 20 years since West Nile virus first emerged in the Americas, and since then, little progress has been made to control outbreaks caused by this virus. After its first detection in New York in 1999, West Nile virus quickly spread across the continent, causing an epidemic of human disease and massive bird die-offs. Now the virus has become endemic to the United States, where an estimated 7 million human infections have occurred, making it the leading mosquito-borne virus infection and the most common cause of viral encephalitis in the country. To bring new attention to one of the most important mosquito-borne viruses in the Americas, we provide an interactive review using Nextstrain: a visualization tool for real-time tracking of pathogen evolution (nextstrain.org/WNV/NA). Nextstrain utilizes a growing database of more than 2,000 West Nile virus genomes and harnesses the power of phylogenetics for students, educators, public health workers, and researchers to visualize key aspects of virus spread and evolution. Using Nextstrain, we use virus genomics to investigate the emergence of West Nile virus in the U S, followed by its rapid spread, evolution in a new environment, establishment of endemic transmission, and subsequent international spread. For each figure, we include a link to Nextstrain to allow the readers to directly interact with and explore the underlying data in new ways. We also provide a brief online narrative that parallels this review to further explain the data and highlight key epidemiological and evolutionary features (nextstrain.org/narratives/twenty-years-of-WNV). Mirroring the dynamic nature of outbreaks, the Nextstrain links provided within this paper are constantly updated as new West Nile virus genomes are shared publicly, helping to stay current with the research. Overall, our review showcases how genomics can track West Nile virus spread and evolution, as well as potentially uncover novel targeted control measures to help alleviate its public health burden.

Effects of Arbovirus Multi-Host Life Cycles on Dinucleotide and Codon Usage Patterns

Arthropod-borne viruses (arboviruses) of vertebrates including dengue, zika, chikungunya, Rift Valley fever, and blue tongue viruses cause extensive morbidity and mortality in humans, agricultural animals, and wildlife across the globe. As obligate intercellular pathogens, arboviruses must be well adapted to the cellular and molecular environment of both their arthropod (invertebrate) and vertebrate hosts, which are vastly different due to hundreds of millions of years of separate evolution. Here we discuss the comparative pressures on arbovirus RNA genomes as a result of a dual host life cycle, focusing on pressures that do not alter amino acids. We summarize what is currently known about arboviral genetic composition, such as dinucleotide and codon usage, and how cyclical infection of vertebrate and invertebrate hosts results in different genetic profiles compared with single-host viruses. To serve as a comparison, we compile what is known about arthropod tRNA, dinucleotide, and codon usages and compare this with vertebrates. Additionally, we discuss the potential roles of genetic robustness in arboviral evolution and how it may vary from other viruses. Overall, both arthropod and vertebrate hosts influence the resulting genetic composition of arboviruses, but a great deal remains to be investigated.

Endless Forms: Within-Host Variation in the Structure of the West Nile Virus RNA Genome during Serial Passage in Bird Hosts

The enzymes that copy RNA genomes lack proofreading, and viruses that possess RNA genomes, such as West Nile virus, rapidly diversify into swarms of mutant lineages within a host. Intrahost variation of the primary genomic sequence of RNA viruses has been studied extensively because the extent of this variation shapes key virus phenotypes. However, RNA genomes also form complex secondary structures based on within-genome nucleotide complementarity, which are critical regulators of the cyclization of the virus genome that is necessary for efficient replication and translation. We sought to characterize variation in these secondary structures within populations of West Nile virus during serial passage in three bird species. Our study indicates that the intrahost population of West Nile virus is a diverse assortment of RNA secondary structures that should be considered in future analyses of intrahost viral diversity, but some regions that are critical for genome cyclization are conserved within hosts. Besides potential impacts on viral replication, structural diversity can influence the efficacy of small RNA antiviral therapies.

RNA viruses are infamous for their high rates of mutation, which produce swarms of genetic variants within individual hosts. To date, analyses of intrahost genetic diversity have focused on the primary genome sequence. However, virus phenotypes are shaped not only by primary sequence but also by the secondary structures into which this sequence folds. Such structures enable viral replication, translation, and binding of small RNAs, yet within-host variation at the structural level has not been adequately explored. We characterized the structural diversity of the 5′ untranslated region (UTR) of populations of West Nile virus (WNV) that had been subject to five serial passages in triplicate in each of three bird species. Viral genomes were sampled from host serum samples at each passage (n = 45 populations) and subjected to next-generation sequencing. For populations derived from passages 1, 3, and 5 (n = 9 populations), we predicted the impact of each mutation occurring at a frequency of ≥1% on the secondary structure of the 5′ UTR. As expected, mutations in double-stranded (DS) regions of the 5′ UTR stem structures caused structural changes of significantly greater magnitude than did mutations in single-stranded (SS) regions. Despite the greater impact of mutations in DS regions, mutations in DS and SS regions occurred at similar frequencies, with no evidence of enhanced selection against mutation in DS regions. In contrast, mutations in two regions that mediate genome cyclization and thereby regulate replication and translation, the 5′ cyclization sequence and the UAR flanking stem (UFS), were suppressed in all three hosts.

IMPORTANCE The enzymes that copy RNA genomes lack proofreading, and viruses that possess RNA genomes, such as West Nile virus, rapidly diversify into swarms of mutant lineages within a host. Intrahost variation of the primary genomic sequence of RNA viruses has been studied extensively because the extent of this variation shapes key virus phenotypes. However, RNA genomes also form complex secondary structures based on within-genome nucleotide complementarity, which are critical regulators of the cyclization of the virus genome that is necessary for efficient replication and translation. We sought to characterize variation in these secondary structures within populations of West Nile virus during serial passage in three bird species. Our study indicates that the intrahost population of West Nile virus is a diverse assortment of RNA secondary structures that should be considered in future analyses of intrahost viral diversity, but some regions that are critical for genome cyclization are conserved within hosts. Besides potential impacts on viral replication, structural diversity can influence the efficacy of small RNA antiviral therapies.

Sustained Wolbachia-mediated blocking of dengue virus isolates following serial passage in Aedes aegypti cell culture

Wolbachia is an intracellular endosymbiont of insects that inhibits the replication of a range of pathogens in its arthropod hosts. The release of Wolbachia into wild populations of mosquitoes is an innovative biocontrol effort to suppress the transmission of arthropod-borne viruses (arboviruses) to humans, most notably dengue virus. The success of the Wolbachia-based approach hinges upon the stable persistence of the ‘pathogen blocking’ effect, whose mechanistic basis is poorly understood. Evidence suggests that Wolbachia may affect viral replication via a combination of competition for host resources and activation of host immunity. The evolution of resistance against Wolbachia and pathogen blocking in the mosquito or the virus could reduce the public health impact of the symbiont releases. Here, we investigate if dengue 3 virus (DENV-3) is capable of accumulating adaptive mutations that improve its replicative capacity during serial passage in Wolbachia wMel-infected cells. During the passaging regime, viral isolates in Wolbachia-infected cells exhibited greater variation in viral loads compared to controls. The viral loads of these isolates declined rapidly during passaging due to the blocking effects of Wolbachia carriage, with several being lost all together and the remainder recovering to low but stable levels. We attempted to sequence the genomes of the surviving passaged isolates but, given their low abundance, were unable to obtain sufficient depth of coverage for evolutionary analysis. In contrast, viral loads in Wolbachia-free control cells were consistently high during passaging. The surviving isolates passaged in the presence of Wolbachia exhibited a reduced ability to replicate even in Wolbachia-free cells. These experiments demonstrate the challenge for dengue in evolving resistance to Wolbachia-mediated blocking.

Experimental evolution: its principles and applications in developing stress-tolerant yeasts

Stress tolerance and resistance in industrial yeast strains are important attributes for cost-effective bioprocessing. The source of stress-tolerant yeasts ranges from extremophilic environments to laboratory engineered strains. However, industrial stress-tolerant yeasts are very rare in nature as the natural environment forces them to evolve traits that optimize survival and reproduction and not the ability to withstand harsh habitat-irrelevant industrial conditions. Experimental evolution is a frequent method used to uncover the mechanisms of evolution and microbial adaption towards environmental stresses. It optimizes biological systems by means of adaptation to environmental stresses and thus has immense power of development of robust stress-tolerant yeasts. This mini-review briefly outlines the basics and implications of evolution experiments and their applications to industrial biotechnology. This work is meant to serve as an introduction to those new to the field of experimental evolution, and as a guide to biologists working in the field of yeast stress response. Future perspectives of experimental evolution for potential biotechnological applications have also been elucidated.

An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar

How viruses evolve within hosts can dictate infection outcomes; however, reconstructing this process is challenging. We evaluate our multiplexed amplicon approach, PrimalSeq, to demonstrate how virus concentration, sequencing coverage, primer mismatches, and replicates influence the accuracy of measuring intrahost virus diversity. We develop an experimental protocol and computational tool, iVar, for using PrimalSeq to measure virus diversity using Illumina and compare the results to Oxford Nanopore sequencing. We demonstrate the utility of PrimalSeq by measuring Zika and West Nile virus diversity from varied sample types and show that the accumulation of genetic diversity is influenced by experimental and biological systems.

Mechanisms and Concepts in RNA Virus Population Dynamics and Evolution

RNA viruses are unique in their evolutionary capacity, exhibiting high mutation rates and frequent recombination. They rapidly adapt to environmental changes, such as shifts in immune pressure or pharmacological challenge. The evolution of RNA viruses has been brought into new focus with the recent developments of genetic and experimental tools to explore and manipulate the evolutionary dynamics of viral populations. These studies have uncovered new mechanisms that enable viruses to overcome evolutionary challenges in the environment and have emphasized the intimate relationship of viral populations with evolution. Here, we review some of the emerging viral and host mechanisms that underlie the evolution of RNA viruses. We also discuss new studies that demonstrate that the relationship between evolutionary dynamics and virus biology spans many spatial and temporal scales, affecting transmission dynamics within and between hosts as well as pathogenesis.

Mapping the Evolutionary Potential of RNA Viruses

The deterministic force of natural selection and stochastic influence of drift shape RNA virus evolution. New deep-sequencing and microfluidics technologies allow us to quantify the effect of mutations and trace the evolution of viral populations with single-genome and single-nucleotide resolution. Such experiments can reveal the topography of the genotype-fitness landscapes that shape the path of viral evolution. By combining historical analyses, like phylogenetic approaches, with high-throughput and high-resolution evolutionary experiments, we can observe parallel patterns of evolution that drive important phenotypic transitions. These developments provide a framework for quantifying and anticipating potential evolutionary events. Here, we examine emerging technologies that can map the selective landscapes of viruses, focusing on their application to pathogenic viruses. We identify areas where these technologies can bolster our ability to study the evolution of viruses and to anticipate and possibly intervene in evolutionary events and prevent viral disease.

Low-fidelity Venezuelan equine encephalitis virus polymerase mutants to improve live-attenuated vaccine safety and efficacy

During RNA virus replication, there is the potential to incorporate mutations that affect virulence or pathogenesis. For live-attenuated vaccines, this has implications for stability, as replication may result in mutations that either restore the wild-type phenotype via reversion or compensate for the attenuating mutations by increasing virulence (pseudoreversion). Recent studies have demonstrated that altering the mutation rate of an RNA virus is an effective attenuation tool. To validate the safety of low-fidelity mutations to increase vaccine attenuation, several mutations in the RNA-dependent RNA-polymerase (RdRp) were tested in the live-attenuated Venezuelan equine encephalitis virus vaccine strain, TC-83. Next generation sequencing after passage in the presence of mutagens revealed a mutant containing three mutations in the RdRp, TC-83 3x, to have decreased replication fidelity, while a second mutant, TC-83 4x displayed no change in fidelity, but shared many phenotypic characteristics with TC-83 3x. Both mutants exhibited increased, albeit inconsistent attenuation in an infant mouse model, as well as increased immunogenicity and complete protection against lethal challenge of an adult murine model compared with the parent TC-83. During serial passaging in a highly permissive model, the mutants increased in virulence but remained less virulent than the parent TC-83. These results suggest that the incorporation of low-fidelity mutations into the RdRp of live-attenuated vaccines for RNA viruses can confer increased immunogenicity whilst showing some evidence of increased attenuation. However, while in theory such constructs may result in more effective vaccines, the instability of the vaccine phenotype decreases the likelihood of this being an effective vaccine strategy.

Mosquitoes Transmit Unique West Nile Virus Populations during Each Feeding Episode

Arthropod-borne viruses (arboviruses), such as Zika virus, chikungunya virus, and West Nile virus (WNV), pose continuous threats to emerge and cause large epidemics. Often, these events are associated with novel virus variants optimized for local transmission that first arise as minorities within a host. Thus, the conditions that regulate the frequency of intrahost variants are important determinants of emergence. Here, we describe the dynamics of WNV genetic diversity during its transmission cycle. By temporally sampling saliva from individual mosquitoes, we demonstrate that virus populations expectorated by mosquitoes are highly diverse and unique to each feeding episode. After transmission to birds, however, most genetic diversity is removed by strong purifying selection. Further, transmission of potentially mosquito-adaptive WNV variants is strongly influenced by genetic drift in mosquitoes. These results highlight the complex evolutionary forces a novel virus variant must overcome to alter infection phenotypes at the population level.

Selective constraint and adaptive potential of West Nile virus within and among naturally infected avian hosts and mosquito vectors

Arthropod-borne viruses are among the most genetically constrained RNA viruses, yet they have a remarkable propensity to adapt and emerge. We studied wild birds and mosquitoes naturally infected with West Nile virus (WNV) in a 'hot spot' of virus transmission in Chicago, IL, USA. We generated full coding WNV genome sequences from spatiotemporally matched bird and mosquito samples using high-throughput sequencing, allowing a molecular evolutionary assessment with deep coverage. Mean FST among samples was 0.66 (±0.02 SE) and was bimodal, with mean nucleotide diversity being higher between samples (interhost πN = 0.001; πS = 0.024) than within them (intrahost πN < 0.0001; πS < 0.001). Eight genomic sites with FST > 1.01 (in the PrM, NS2a, NS3, NS4b, and 5'-noncoding genomic regions) showed bird versus mosquito variant frequency differences of >30 per cent and/or polymorphisms fixed in ≥5 host or vector individuals, suggesting host tropism for these variants. However, phylogenetic analyses demonstrated a lack of grouping by bird or mosquito, most inter-sample differences were synonymous (mean interhost πN/πS = 0.04), and there was no significant difference between hosts and vectors in either their nucleotide diversities or levels of purifying selection (mean intrahost πN/πS = 0.28 in birds and πN/πS = 0.21 in mosquitoes). This finding contrasts with the 'trade-off' and 'selective sieve' hypotheses that have been proposed and tested in the laboratory, which predict strong host versus vector effects on WNV genetic variation, with heightened selective constraint in birds alternating with heightened viral diversity in mosquitoes. Overall, our data show WNV to be highly selectively constrained within and between both hosts and vectors but still able to vary at a limited number of sites across the genome. Such site-specific plasticity in the face of overall selective constraint may offer a mechanism whereby highly constrained viruses such as WNV and its relatives can still adapt and emerge.

African and Asian Zika Virus Isolates Display Phenotypic Differences Both In Vitro and In Vivo

Zika virus (ZIKV) is a mosquito-borne member of the genus Flavivirus that has emerged since 2007 to cause outbreaks in Africa, Asia, Oceania, and most recently, in the Americas. Here, we used an isolate history as well as genetic and phylogenetic analyses to characterize three low-passage isolates representing African (ArD 41525) and Asian (CPC-0740, SV0127-14) lineages to investigate the potential phenotypic differences in vitro and in vivo. The African isolate displayed a large plaque phenotype (∼3-4 mm) on Vero and HEK-293 cells, whereas the Asian isolates either exhibited a small plaque phenotype (∼1-2 mm) or did not produce any plaques. In multistep replication kinetics in nine different vertebrate and insect cell lines, the African isolate consistently displayed faster replication kinetics and yielded ∼10- to 10,000-fold higher peak virus titers (infectious or RNA copies) compared with the Asian isolates. Oral exposure of Aedes aegypti mosquitoes with the African isolate yielded higher infection and dissemination rates compared with the Asian isolates. Infection of Ifnar1^-/- mice with the African isolate produced a uniformly fatal disease, whereas infection with the Asian isolates produced either a delay in time-to-death or a significantly lower mortality rate. Last, the African isolate was > 10,000-fold more virulent than the Asian isolates in an interferon type I antibody blockade mouse model. These data demonstrate substantial phenotypic differences between low-passage African and Asian isolates both in vitro and in vivo and warrant further investigation. They also highlight the need for basic characterization of ZIKV isolates, as the utilization of the uncharacterized isolates could have consequences for animal model and therapeutic/vaccine development.

Metagenomic Investigation of Idiopathic Meningoencephalomyelitis in Dogs

BACKGROUND

Meningoencephalomyelitis of unknown origin (MUO) is a common and life-threatening neuroinflammatory disease in dogs. Features of the disease are suggestive of an underlying immune-mediated process, but the association of this disease with a pathogen is still unknown.

HYPOTHESIS/OBJECTIVES

To search for candidate etiologic agent associated with cases if MUO using next generation metagenomic sequencing.

ANIMALS

Twenty-two dogs diagnosed with either MUO (11/22; 10 CSF and 3 brain), or noninflammatory CNS diseases inconsistent with MUO (11/22; 11 CSF and 2 brain) that served as negative controls.

METHODS

A case control study was performed by identifying MUO and non-MUO cases. Samples were blindly processed and then unblinded for comparative analyses. Inclusion criteria for MUO cases included consistent MRI lesions and inflammatory CSF with a negative PCR panel for infectious agents or histopathologic diagnosis. Dogs with glucocorticoid therapy within 2 weeks of sample collection were excluded. Fresh-frozen cerebrospinal fluid (CSF; 21) and brain (5) samples were collected and RNA and DNA were extracted separately for shotgun metagenomic sequencing. Known positive samples were used as controls to validate our sequencing and analysis pipelines and to establish limits of detection. Sequencing results were analyzed at a nucleotide and protein level for broad comparison to known infectious organisms.

RESULTS

No candidate etiologic agents were identified in dogs with MUO.

CONCLUSIONS AND CLINICAL IMPORTANCE

These results support but do not prove the hypothesis that MUO is not associated with infectious agents and might be an autoimmune disease.

Full-genome dengue virus sequencing in mosquito saliva shows lack of convergent positive selection during transmission by Aedes aegypti

Like other pathogens with high mutation and replication rates, within-host dengue virus (DENV) populations evolve during infection of their main mosquito vector, Aedes aegypti. Within-host DENV evolution during transmission provides opportunities for adaptation and emergence of novel virus variants. Recent studies of DENV genetic diversity failed to detect convergent evolution of adaptive mutations in mosquito tissues such as midgut and salivary glands, suggesting that convergent positive selection is not a major driver of within-host DENV evolution in the vector. However, it is unknown whether this conclusion extends to the transmitted viral subpopulation because it is technically difficult to sequence DENV genomes in mosquito saliva. Here, we achieved DENV full-genome sequencing by pooling saliva samples collected non-sacrificially from 49 to 163 individual Ae. aegypti mosquitoes previously infected with one of two DENV-1 genotypes. We compared the transmitted viral subpopulations found in the pooled saliva samples collected in time series with the input viral population present in the infectious blood meal. In all pooled saliva samples examined, the full-genome consensus sequence of the input viral population was unchanged. Although the pooling strategy prevents analysis of individual saliva samples, our results demonstrate the lack of strong convergent positive selection during a single round of DENV transmission by Ae. aegypti. This finding reinforces the idea that genetic drift and purifying selection are the dominant evolutionary forces shaping within-host DENV genetic diversity during transmission by mosquitoes.

Promiscuous viruses-how do viruses survive multiple unrelated hosts?

Arthropod-borne viruses (arboviruses) require efficient replication in taxonomically divergent hosts in order to perpetuate in nature. This review discusses recent advances in our understanding of the phylogenetic position of arthropod-borne viruses relative to insect-specific viruses, which appear to be more common and ecological requirements for successful adoption of the 'arbovirus phenotype.' Several molecular and other mechanisms that permit replication in divergent hosts are also discussed.

Transmission bottlenecks and RNAi collectively influence tick-borne flavivirus evolution

Arthropod-borne RNA viruses exist within hosts as heterogeneous populations of viral variants and, as a result, possess great genetic plasticity. Understanding the micro-evolutionary forces shaping these viruses can provide insights into how they emerge, adapt, and persist in new and changing ecological niches. While considerable attention has been directed toward studying the population dynamics of mosquito-borne viruses, little is known about tick-borne virus populations. Therefore, using a mouse and Ixodes scapularis tick transmission model, we examined Powassan virus (POWV; Flaviviridae, Flavivirus) populations in and between both the vertebrate host and arthropod vector. We found that genetic bottlenecks, RNAi-mediated diversification, and selective constraints collectively influence POWV evolution. Together, our data provide a mechanistic explanation for the slow, long-term evolutionary trends of POWV, and suggest that all arthropod-borne viruses encounter similar selective pressures at the molecular level (i.e. RNAi), yet evolve much differently due to their unique rates and modes of transmission.

Vector Competence of American Mosquitoes for Three Strains of Zika Virus

In 2015, Zika virus (ZIKV; Flaviviridae; Flavivirus) emerged in the Americas, causing millions of infections in dozens of countries. The rapid spread of the virus and the association with disease outcomes such as Guillain-Barré syndrome and microcephaly make understanding transmission dynamics essential. Currently, there are no reports of vector competence (VC) of American mosquitoes for ZIKV isolates from the Americas. Further, it is not clear whether ZIKV strains from other genetic lineages can be transmitted by American Aedes aegypti populations, and whether the scope of the current epidemic is in part facilitated by viral factors such as enhanced replicative fitness or increased vector competence. Therefore, we characterized replication of three ZIKV strains, one from each of the three phylogenetic clades in several cell lines and assessed their abilities to be transmitted by Ae. aegypti mosquitoes. Additionally, laboratory colonies of different Culex spp. were infected with an American outbreak strain of ZIKV to assess VC. Replication rates were variable and depended on virus strain, cell line and MOI. African strains used in this study outcompeted the American strain in vitro in both mammalian and mosquito cell culture. West and East African strains of ZIKV tested here were more efficiently transmitted by Ae. aegypti from Mexico than was the currently circulating American strain of the Asian lineage. Long-established laboratory colonies of Culex mosquitoes were not efficient ZIKV vectors. These data demonstrate the capacity for additional ZIKV strains to infect and replicate in American Aedes mosquitoes and suggest that neither enhanced virus replicative fitness nor virus adaptation to local vector mosquitoes seems likely to explain the extent and intensity of ZIKV transmission in the Americas.

The mechanisms contributing to the explosive nature of the current ZIKV outbreak in the Americas are poorly understood. Therefore, we characterized the replication of three strains, one from each phylogenetic clade of ZIKV and evaluated virus strain differences in transmission efficiency by American mosquitoes. Our results suggest that the strain currently circulating in the Americas does not have unusually high infectivity for American Ae. aegypti as compared to the African strains used in this study. Colonized Culex mosquitoes were inefficient vectors. In vitro data suggested slower replication and decreased fitness for the currently circulating American strain compared to African strains isolated decades ago. Therefore, viral adaptation to local mosquitoes does not appear, at present, to be driving the current ZIKV outbreak in the Americas.

Dynamics of West Nile virus evolution in mosquito vectors

West Nile virus remains the most common cause of arboviral encephalitis in North America. Since it was introduced, it has undergone adaptive genetic change as it spread throughout the continent. The WNV transmission cycle is relatively tractable in the laboratory. Thus the virus serves as a convenient model system for studying the population biology of mosquito-borne flaviviruses as they undergo transmission to and from mosquitoes and vertebrates. This review summarizes the current knowledge regarding the population dynamics of this virus within mosquito vectors.

Transmission and evolution of tick-borne viruses

Ticks transmit a diverse array of viruses such as tick-borne encephalitis virus, Powassan virus, and Crimean-Congo hemorrhagic fever virus that are reemerging in many parts of the world. Most tick-borne viruses (TBVs) are RNA viruses that replicate using error-prone polymerases and produce genetically diverse viral populations that facilitate their rapid evolution and adaptation to novel environments. This article reviews the mechanisms of virus transmission by tick vectors, the molecular evolution of TBVs circulating in nature, and the processes shaping viral diversity within hosts to better understand how these viruses may become public health threats. In addition, remaining questions and future directions for research are discussed.

Measurements of Intrahost Viral Diversity Are Extremely Sensitive to Systematic Errors in Variant Calling

With next-generation sequencing technologies, it is now feasible to efficiently sequence patient-derived virus populations at a depth of coverage sufficient to detect rare variants. However, each sequencing platform has characteristic error profiles, and sample collection, target amplification, and library preparation are additional processes whereby errors are introduced and propagated. Many studies account for these errors by using ad hoc quality thresholds and/or previously published statistical algorithms. Despite common usage, the majority of these approaches have not been validated under conditions that characterize many studies of intrahost diversity. Here, we use defined populations of influenza virus to mimic the diversity and titer typically found in patient-derived samples. We identified single-nucleotide variants using two commonly employed variant callers, DeepSNV and LoFreq. We found that the accuracy of these variant callers was lower than expected and exquisitely sensitive to the input titer. Small reductions in specificity had a significant impact on the number of minority variants identified and subsequent measures of diversity. We were able to increase the specificity of DeepSNV to >99.95% by applying an empirically validated set of quality thresholds. When applied to a set of influenza virus samples from a household-based cohort study, these changes resulted in a 10-fold reduction in measurements of viral diversity. We have made our sequence data and analysis code available so that others may improve on our work and use our data set to benchmark their own bioinformatics pipelines. Our work demonstrates that inadequate quality control and validation can lead to significant overestimation of intrahost diversity.

IMPORTANCE

Advances in sequencing technology have made it feasible to sequence patient-derived viral samples at a level sufficient for detection of rare mutations. These high-throughput, cost-effective methods are revolutionizing the study of within-host viral diversity. However, the techniques are error prone, and the methods commonly used to control for these errors have not been validated under the conditions that characterize patient-derived samples. Here, we show that these conditions affect measurements of viral diversity. We found that the accuracy of previously benchmarked analysis pipelines was greatly reduced under patient-derived conditions. By carefully validating our sequencing analysis using known control samples, we were able to identify biases in our method and to improve our accuracy to acceptable levels. Application of our modified pipeline to a set of influenza virus samples from a cohort study provided a realistic picture of intrahost diversity and suggested the need for rigorous quality control in such studies.

Genetic Drift, Purifying Selection and Vector Genotype Shape Dengue Virus Intra-host Genetic Diversity in Mosquitoes

Due to their error-prone replication, RNA viruses typically exist as a diverse population of closely related genomes, which is considered critical for their fitness and adaptive potential. Intra-host demographic fluctuations that stochastically reduce the effective size of viral populations are a challenge to maintaining genetic diversity during systemic host infection. Arthropod-borne viruses (arboviruses) traverse several anatomical barriers during infection of their arthropod vectors that are believed to impose population bottlenecks. These anatomical barriers have been associated with both maintenance of arboviral genetic diversity and alteration of the variant repertoire. Whether these patterns result from stochastic sampling (genetic drift) rather than natural selection, and/or from the influence of vector genetic heterogeneity has not been elucidated. Here, we used deep sequencing of full-length viral genomes to monitor the intra-host evolution of a wild-type dengue virus isolate during infection of several mosquito genetic backgrounds. We estimated a bottleneck size ranging from 5 to 42 founding viral genomes at initial midgut infection, irrespective of mosquito genotype, resulting in stochastic reshuffling of the variant repertoire. The observed level of genetic diversity increased following initial midgut infection but significantly differed between mosquito genetic backgrounds despite a similar initial bottleneck size. Natural selection was predominantly negative (purifying) during viral population expansion. Taken together, our results indicate that dengue virus intra-host genetic diversity in the mosquito vector is shaped by genetic drift and purifying selection, and point to a novel role for vector genetic factors in the genetic breadth of virus populations during infection. Identifying the evolutionary forces acting on arboviral populations within their arthropod vector provides novel insights into arbovirus evolution.

During infection of their arthropod vectors, arthropod-borne viruses (arboviruses) such as dengue viruses traverse several anatomical barriers that are believed to cause dramatic reductions in population size. Such population bottlenecks challenge the maintenance of viral genetic diversity, which is considered critical for fitness and adaptability of arboviruses. Anatomical barriers in the vector were previously associated with both maintenance of arboviral genetic diversity and alteration of the variant repertoire. However, the relative role of random processes and natural selection, and the influence of vector genetic heterogeneity have not been elucidated. In this study, we used high-throughput sequencing to monitor dengue virus genetic diversity during infection of several genetic backgrounds of their mosquito vector. Our results show that initial infection of the vector is randomly founded by only a few tens of individual virus genomes. The overall level of viral genetic diversity generated during infection was predominantly under purifying selection but differed significantly between mosquito genetic backgrounds. Thus, in addition to random evolutionary forces and the purging of deleterious mutations that shape dengue virus genetic diversity during vector infection, our results also point to a novel role for vector genetic factors in the genetic breadth of virus populations.

West Nile Virus Fitness Costs in Different Mosquito Species

West Nile virus (WNV) remains an important public health problem causing annual epidemics in the United States. Grubaugh et al. observed that WNV genetic divergence is dependent on the vector mosquito species. This suggests that specific WNV vector-bird species pairings may generate novel genotypes that could promote outbreaks.

Genetic Drift during Systemic Arbovirus Infection of Mosquito Vectors Leads to Decreased Relative Fitness during Host Switching

The emergence of mosquito-borne RNA viruses, such as West Nile virus (WNV), is facilitated by genetically complex virus populations within hosts. Here, we determine whether WNV enzootic (Culex tarsalis, Cx. quinquefasciatus, and Cx. pipiens) and bridge vectors (Aedes aegypti) have differential impacts on viral mutational diversity and fitness. During systemic mosquito infection, WNV faced stochastic reductions in genetic diversity that rapidly was recovered during intra-tissue population expansions. Interestingly, this intrahost selection and diversification was mosquito species dependent with Cx. tarsalis and Cx. quinquefasciatus exhibiting greater WNV divergence. However, recovered viral populations contained a preponderance of potentially deleterious mutations (i.e., high mutational load) and had lower relative fitness in avian cells compared to input virus. These findings demonstrate that the adaptive potential associated with mosquito transmission varies depending on the mosquito species and carries a significant fitness cost in vertebrates.

Isolation of a Novel Insect-Specific Flavivirus from Culiseta melanura in the Northeastern United States

The genus Flavivirus includes a number of newly recognized viruses that infect and replicate only within mosquitoes. To determine whether insect-specific flaviviruses (ISFs) may infect Culiseta (Cs.) melanura mosquitoes, we screened pools of field-collected mosquitoes for virus infection by RT-PCR targeting conserved regions of the NS5 gene. NS5 nucleotide sequences amplified from Cs. melanura pools were genetically similar to other ISFs and most closely matched Calbertado virus from Culex tarsalis, sharing 68.7% nucleotide and 76.1% amino acid sequence identity. The complete genome of one virus isolate was sequenced to reveal a primary open reading frame (ORF) encoding a viral polyprotein characteristic of the genus Flavivirus. Phylogenetic analysis showed that this virus represents a distinct evolutionary lineage that belongs to the classical ISF group. The virus was detected solely in Cs. melanura pools, occurred in sampled populations from Connecticut, New York, New Hampshire, and Maine, and infected both adult and larval stages of the mosquito. Maximum likelihood estimate infection rates (MLE-IR) were relatively stable in overwintering Cs. melanura larvae collected monthly from November of 2012 through May of 2013 (MLE-IR = 0.7-2.1/100 mosquitoes) and in host-seeking females collected weekly from June through October of 2013 (MLE-IR = 3.8-11.5/100 mosquitoes). Phylogenetic analysis of viral sequences revealed limited genetic variation that lacked obvious geographic structure among strains in the northeastern United States. This new virus is provisionally named Culiseta flavivirus on the basis of its host association with Cs. melanura.

West Nile Virus Population Structure, Injury, and Interferon-Stimulated Gene Expression in the Brain From a Fatal Case of Encephalitis

Background. West Nile virus (WNV) infection in humans can result in severe, acute encephalitis typically involving subcortical gray matter brain regions. West Nile virus replication within specific human brain regions from a human case of acute encephalitis has not been studied.

Methods. We describe a fatal case of WNV encephalitis in which we obtained tissue from specific brain regions at autopsy to evaluate viral-host interactions using next-generation sequencing and immunohistochemistry analysis.

Results. We found that WNV populations in the injured subcortical brain regions exhibited increased amino acid variation and increased expression of specific interferon genes compared with cortical tissues despite similar viral burden.

Conclusions. These observational, patient-based data suggest that neuronal injury and the strength of viral selection pressure may be associated with the level of the innate immune responses. Further studies in human and animal models evaluating the role of innate immune responses on injury patterns and viral selection pressure are needed.

Next-generation sequencing shows West Nile virus quasispecies diversification after a single passage in a carrion crow (Corvus corone) in vivo infection model

West Nile virus (WNV) occurs as a population of genetic variants (quasispecies) infecting a single animal. Previous low-resolution viral genetic diversity estimates in sampled wild birds and mosquitoes, and in multiple-passage adaptation studies in vivo or in cell culture, suggest that WNV genetic diversification is mostly limited to the mosquito vector. This study investigated genetic diversification of WNV in avian hosts during a single passage using next-generation sequencing. Wild-captured carrion crows were subcutaneously infected using a clonal Middle-East WNV. Blood samples were collected 2 and 4 days post-infection. A reverse-transcription (RT)-PCR approach was used to amplify the WNV genome directly from serum samples prior to next-generation sequencing resulting in an average depth of at least 700 ×  in each sample. Appropriate controls were sequenced to discriminate biologically relevant low-frequency variants from experimentally introduced errors. The WNV populations in the wild crows showed significant diversification away from the inoculum virus quasispecies structure. By contrast, WNV populations in intracerebrally infected day-old chickens did not diversify from that of the inoculum. Where previous studies concluded that WNV genetic diversification is only experimentally demonstrated in its permissive insect vector species, we have experimentally shown significant diversification of WNV populations in a wild bird reservoir species.