Population variation of West Nile virus confers a host-specific fitness benefit in mosquitoes

Loss of West Nile virus genetic diversity during mosquito infection due to species-dependent population bottlenecks

Vector competence (VC) refers to the efficiency of pathogen transmission by vectors. Each step in the infection of a mosquito vector constitutes a barrier to transmission that may impose bottlenecks on virus populations. West Nile virus (WNV) is maintained by multiple mosquito species with varying VC. However, the extent to which bottlenecks and VC are linked is poorly understood. Similarly, quantitative analyses of mosquito-imposed bottlenecks on virus populations are limited. We used molecularly barcoded WNV to quantify tissue-associated population bottlenecks in three variably competent WNV vectors. Our results confirm strong population bottlenecks during mosquito infection that are capable of dramatically reshaping virus population structure in a non-selective manner. In addition, we found that mosquitoes with differing VC uniquely shape WNV population structure: highly competent vectors are more likely to contribute to the maintenance of rare viral genotypes. These findings have important implications for arbovirus emergence and evolution.

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.

An in vitro workflow to create and modify infectious clones using replication cycle reaction

Reverse genetics systems are critical tools in combating emerging viruses which enable a better understanding of the genetic mechanisms by which viruses cause disease. Traditional cloning approaches using bacteria are fraught with difficulties due to the bacterial toxicity of many viral sequences, resulting in unwanted mutations within the viral genome. Here, we describe a novel in vitro workflow that leverages gene synthesis and replication cycle reaction to produce a supercoiled infectious clone plasmid that is easy to distribute and manipulate. We developed two infectious clones as proof of concept: a low passage dengue virus serotype 2 isolate (PUO-218) and the USA-WA1/2020 strain of SARS-CoV-2, which replicated similarly to their respective parental viruses. Furthermore, we generated a medically relevant mutant of SARS-CoV-2, Spike D614G. Results indicate that our workflow is a viable method to generate and manipulate infectious clones for viruses that are notoriously difficult for traditional bacterial-based cloning methods.

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.

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.

Experimental Evolution of West Nile Virus at Higher Temperatures Facilitates Broad Adaptation and Increased Genetic Diversity

West Nile virus (WNV, Flaviviridae, Flavivirus) is a mosquito-borne flavivirus introduced to North America in 1999. Since 1999, the Earth's average temperature has increased by 0.6 °C. Mosquitoes are ectothermic organisms, reliant on environmental heat sources. Temperature impacts vector-virus interactions which directly influence arbovirus transmission. RNA viral replication is highly error-prone and increasing temperature could further increase replication rates, mutation frequencies, and evolutionary rates. The impact of temperature on arbovirus evolutionary trajectories and fitness landscapes has yet to be sufficiently studied. To investigate how temperature impacts the rate and extent of WNV evolution in mosquito cells, WNV was experimentally passaged 12 times in Culex tarsalis cells, at 25 °C and 30 °C. Full-genome deep sequencing was used to compare genetic signatures during passage, and replicative fitness was evaluated before and after passage at each temperature. Our results suggest adaptive potential at both temperatures, with unique temperature-dependent and lineage-specific genetic signatures. Further, higher temperature passage was associated with significantly increased replicative fitness at both temperatures and increases in nonsynonymous mutations. Together, these data indicate that if similar selective pressures exist in natural systems, increases in temperature could accelerate emergence of high-fitness strains with greater phenotypic plasticity.

Pathogenicity and virulence of West Nile virus revisited eight decades after its first isolation

West Nile virus (WNV) is a flavivirus which transmission cycle is maintained between mosquitoes and birds, although it occasionally causes sporadic outbreaks in horses and humans that can result in serious diseases and even death. Since its first isolation in Africa in 1937, WNV had been considered a neglected pathogen until its recent spread throughout Europe and the colonization of America, regions where it continues to cause outbreaks with severe neurological consequences in humans and horses. Although our knowledge about the characteristics and consequences of the virus has increased enormously lately, many questions remain to be resolved. Here, we thoroughly update our knowledge of different aspects of the WNV life cycle: virology and molecular classification, host cell interactions, transmission dynamics, host range, epidemiology and surveillance, immune response, clinical presentations, pathogenesis, diagnosis, prophylaxis (antivirals and vaccines), and prevention, and we highlight those aspects that are still unknown and that undoubtedly require further investigation.

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.

Biological Adaptations Associated with Dehydration in Mosquitoes

Diseases that are transmitted by mosquitoes are a tremendous health and socioeconomic burden with hundreds of millions of people being impacted by mosquito-borne illnesses annually. Many factors have been implicated and extensively studied in disease transmission dynamics, but knowledge regarding how dehydration impacts mosquito physiology, behavior, and resulting mosquito-borne disease transmission remain underdeveloped. The lapse in understanding on how mosquitoes respond to dehydration stress likely obscures our ability to effectively study mosquito physiology, behavior, and vectorial capabilities. The goal of this review is to develop a profile of factors underlying mosquito biology that are altered by dehydration and the implications that are related to disease transmission.

The role of co-infection and swarm dynamics in arbovirus transmission

Arthropod-borne viruses (arboviruses) are transmitted by hematophagous insects, primarily mosquitoes. The geographic range and prevalence of mosquito-borne viruses and their vectors has dramatically increased over the last 50 years. As a result, the most medically important arboviurses now co-exist in many regions, resulting in an increased frequency of co-infections in hosts and vectors. In addition to concurrent infections with human pathogens, mosquito-only viruses and/or enzootic viruses not associated with human disease are ubiquitous in mosquito populations. Moreover, mosquito-borne viruses are largely RNA viruses that exist within individual hosts as a diverse and dynamic swarm of closely related genotypes. Interactions among co-infecting viruses and genotypes can have profound effects on virulence, fitness and evolution. Here, we review our understanding of how these complex interactions influence transmission of mosquito-borne viruses, focusing on the often-neglected virus interactions in the mosquito vector, and identify gaps in our knowledge that should guide future studies.

Growth and adaptation of Zika virus in mammalian and mosquito cells

The recent emergence of Zika virus (ZIKV) in the Americas coincident with increased caseloads of microcephalic infants and Guillain-Barre syndrome has prompted a flurry of research on ZIKV. Much of the research is difficult to compare or repeat because individual laboratories use different virus isolates, growth conditions, and quantitative assays. Here we obtained three readily available contemporary ZIKV isolates and the prototype Ugandan isolate. We generated stocks of each on Vero mammalian cells (ZIKVmam) and C6/36 mosquito cells (ZIKVmos), determined titers by different assays side-by-side, compared growth characteristics using one-step and multi-step growth curves on Vero and C6/36 cells, and examined plaque phenotype. ZIKV titers consistently peaked earlier on Vero cells than on C6/36 cells. Contemporary ZIKV isolates reached peak titer most quickly in a multi-step growth curve when the amplifying cell line was the same as the titering cell line (e.g., ZIKVmam titered on Vero cells). Growth of ZIKVmam on mosquito cells was particularly delayed. These data suggest that the ability to infect and/or replicate in insect cells is limited after growth in mammalian cells. In addition, ZIKVmos typically had smaller, more homogenous plaques than ZIKVmam in a standard plaque assay. We hypothesized that the plaque size difference represented early adaptation to growth in mammalian cells. We plaque purified representative-sized plaques from ZIKVmos and ZIKVmam. ZIKVmos isolates maintained the initial phenotype while plaques from ZIKVmam isolates became larger with passaging. Our results underscore the importance of the cells used to produce viral stocks and the potential for adaptation with minimal cell passages. In addition, these studies provide a foundation to compare current and emerging ZIKV isolates in vitro and in vivo.

The synergistic effect of nsP2-L618, nsP3-R117, and E2-K187 on the large plaque phenotype of chikungunya virus

Chikungunya virus (CHIKV), a mosquito-borne Alphavirus, is the etiological agent of chikungunya fever. CHIKV re-emerged from 2004 onwards, and subsequently caused major outbreaks in many parts of the world including the Indian Ocean islands, Asia, and the Americas. In this study, a large plaque variant of CHIKV isolated from patient in Thailand was subjected to repeated cycles of plaque-purification in Vero cells. The resulting virus produced homogenous large plaques and showed a more pathogenic phenotype than the parental wild-type CHIKV. Whole genome analysis of the large plaque virus in comparison to parental isolate revealed a number of mutations, leading to the following amino acid changes: nsP2 (P⁶¹⁸→L), nsP3 (G¹¹⁷→R), and E2 (N¹⁸⁷→K). Eight recombinant CHIKVs were constructed to determine which amino acids mediated the large plaque phenotype. The results showed the recombinant virus which contains all three mutations, rCHK-L, produced significantly larger plaques than the other recombinant viruses (p < 0.01). Moreover, the plaque size of the other recombinant virus tended to be smaller if they contained only one or two of the large plaque associated mutations in the viral genome. In conclusion, the combination of all three residues (nsP2-L⁶¹⁸, nsP3-R¹¹⁷, and E2-K¹⁸⁷) is required to produce the large plaque phenotype of CHIKV.

Mutagen resistance and mutation restriction of St. Louis encephalitis virus

The error rate of the RNA-dependent RNA polymerase (RdRp) of RNA viruses is important in maintaining genetic diversity for viral adaptation and fitness. Numerous studies have shown that mutagen-resistant RNA virus variants display amino acid mutations in the RdRp and other replicase subunits, which in turn exhibit an altered fidelity phenotype affecting viral fitness, adaptability and pathogenicity. St. Louis encephalitis virus (SLEV), like its close relative West Nile virus, is a mosquito-borne flavivirus that has the ability to cause neuroinvasive disease in humans. Here, we describe the successful generation of multiple ribavirin-resistant populations containing a shared amino acid mutation in the SLEV RdRp (E416K). These E416K mutants also displayed resistance to the antiviral T-1106, an RNA mutagen similar to ribavirin. Structural modelling of the E416K polymerase mutation indicated its location in the pinky finger domain of the RdRp, distant from the active site. Deep sequencing of the E416K mutant revealed lower genetic diversity than wild-type SLEV after growth in both vertebrate and invertebrate cells. Phenotypic characterization showed that E416K mutants displayed similar or increased replication in mammalian cells, as well as modest attenuation in mosquito cells, consistent with previous work with West Nile virus high-fidelity variants. In addition, attenuation was limited to mosquito cells with a functional RNA interference response, suggesting an impaired capacity to escape RNA interference could contribute to attenuation of high-fidelity variants. Our results provide increased evidence that RNA mutagen resistance arises through modulation of the RdRp and give further insight into the consequences of altered fidelity of flaviviruses.

High levels of local inter- and intra-host genetic variation of West Nile virus and evidence of fine-scale evolutionary pressures

West Nile virus (WNV; Flaviviridae, Flavivirus) has been endemic in New York State (NYS) since its 1999 introduction, yet prevalence in Culex mosquitoes varies substantially over small spatial and temporal scales. It is unclear if viral genetics plays a role in this variability, as genetic and phenotypic characterization on local scales has generally been lacking. In addition, intrahost diversity of circulating strains have not been fully characterized despite the documented role of minority variants in viral fitness and virulence. In an effort to characterize WNV variability within epidemiologically relevant scales, we performed phylogenetic analyses on NYS isolates from 1999 to 2012. In addition, we performed full-genome, deep-sequencing and genetic analyses on 15 WNV strains isolated in 2012 from Cx. pipiens in an endemic focus of Suffolk County, NY. Our results indicate continued evolution and seasonal maintenance in NYS, yet also widespread mixing and high levels of genetic diversity within geographic foci and individual seasons. Well supported local clusters with shared amino acid differences were identified and suggest local evolutionary pressures and the potential for phenotypic variability. Intrahost diversity of focal isolates was also high, with polymorphism at levels >1.0% identified in approximately 10% of the WNV genome. Although most minority mutations were unique, mutational hotspots shared among local isolates were identified, particularly in C, NS1 and NS2A genes. The most polymorphic region, positions 3198-3388 of the NS1 gene, was comprised predominately of non-synonymous mutations, suggesting a selective advantage for amino acid diversity in this region.

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.

An Infectious cDNA Clone of Zika Virus to Study Viral Virulence, Mosquito Transmission, and Antiviral Inhibitors

The Asian lineage of Zika virus (ZIKV) has recently caused epidemics and severe disease. Unraveling the mechanisms causing increased viral transmissibility and disease severity requires experimental systems. We report an infectious cDNA clone of ZIKV that was generated using a clinical isolate of the Asian lineage. The cDNA clone-derived RNA is infectious in cells, generating recombinant ZIKV. The recombinant virus is virulent in established ZIKV mouse models, leading to neurological signs relevant to human disease. Additionally, recombinant ZIKV is infectious for Aedes aegypti and thus provides a means to examine virus transmission. The infectious cDNA clone was further used to generate a luciferase ZIKV that exhibited sensitivity to a panflavivirus inhibitor, highlighting its potential utility for antiviral screening. This ZIKV reverse genetic system, together with mouse and mosquito infection models, may help identify viral determinants of human virulence and mosquito transmission as well as inform vaccine and therapeutic strategies.

Analysis of coevolution in nonstructural proteins of chikungunya virus

Background

RNA viruses are characterized by high rate of mutations mainly due to the lack of proofreading repair activities associated with its RNA-dependent RNA-polymerase (RdRp). In case of arboviruses, this phenomenon has lead to the existence of mixed population of genomic variants within the host called quasi-species. The stability of strains within the quasi-species lies on mutations that are positively selected which in turn depend on whether these mutations are beneficial in either or both hosts. Coevolution of amino acids (aa) is one phenomenon that leads to establishment of favorable traits in viruses and leading to their fitness.

Results

Fourteen CHIKV clinical samples collected over three years were subjected to RT-PCR, the four non-structural genes amplified and subjected to various genetic analyses. Coevolution analysis showed 30 aa pairs coevolving in nsP1, 23 aa pairs coevolving in nsP2, 239 in nsP3 and 46 aa coevolving pairs in nsP4 when each non-structural protein was considered independently. Further analysis showed that 705 amino acids pairs of the non-structural polyproteins coevolved together with a correlation coefficient of ≥0.5. Functional relevance of these coevolving amino acids in all the nonstructural proteins of CHIKV were predicted using Eukaryotic Linear Motifs (ELMs) of human.

Conclusions

The present study was undertaken to study co-evolving amino acids in the non-structural proteins of chikungunya virus (CHIKV), an important arbovirus. It was observed that several amino acids residues were coevolving and shared common functions.

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.

Tracking Dengue Virus Intra-host Genetic Diversity during Human-to-Mosquito Transmission

Dengue virus (DENV) infection of an individual human or mosquito host produces a dynamic population of closely-related sequences. This intra-host genetic diversity is thought to offer an advantage for arboviruses to adapt as they cycle between two very different host species, but it remains poorly characterized. To track changes in viral intra-host genetic diversity during horizontal transmission, we infected Aedes aegypti mosquitoes by allowing them to feed on DENV2-infected patients. We then performed whole-genome deep-sequencing of human- and matched mosquito-derived DENV samples on the Illumina platform and used a sensitive variant-caller to detect single nucleotide variants (SNVs) within each sample. >90% of SNVs were lost upon transition from human to mosquito, as well as from mosquito abdomen to salivary glands. Levels of viral diversity were maintained, however, by the regeneration of new SNVs at each stage of transmission. We further show that SNVs maintained across transmission stages were transmitted as a unit of two at maximum, suggesting the presence of numerous variant genomes carrying only one or two SNVs each. We also present evidence for differences in selection pressures between human and mosquito hosts, particularly on the structural and NS1 genes. This analysis provides insights into how population drops during transmission shape RNA virus genetic diversity, has direct implications for virus evolution, and illustrates the value of high-coverage, whole-genome next-generation sequencing for understanding viral intra-host genetic diversity.

Analysis of Dengue Virus Genetic Diversity during Human and Mosquito Infection Reveals Genetic Constraints

Dengue viruses (DENV) cause debilitating and potentially life-threatening acute disease throughout the tropical world. While drug development efforts are underway, there are concerns that resistant strains will emerge rapidly. Indeed, antiviral drugs that target even conserved regions in other RNA viruses lose efficacy over time as the virus mutates. Here, we sought to determine if there are regions in the DENV genome that are not only evolutionarily conserved but genetically constrained in their ability to mutate and could hence serve as better antiviral targets. High-throughput sequencing of DENV-1 genome directly from twelve, paired dengue patients’ sera and then passaging these sera into the two primary mosquito vectors showed consistent and distinct sequence changes during infection. In particular, two residues in the NS5 protein coding sequence appear to be specifically acquired during infection in Ae. aegypti but not Ae. albopictus. Importantly, we identified a region within the NS3 protein coding sequence that is refractory to mutation during human and mosquito infection. Collectively, these findings provide fresh insights into antiviral targets and could serve as an approach to defining evolutionarily constrained regions for therapeutic targeting in other RNA viruses.

Dengue viruses cause debilitating and potentially life-threatening acute disease throughout the tropical world. While drug development efforts are underway, there are concerns that drug-resistant strains will emerge rapidly. Indeed, many antiviral drugs for other RNA viruses lose efficacy over time as the virus mutates. Here, we sought to determine if there are regions in the dengue virus genome that are constrained in their ability to mutate and could therefore serve as better targets for antiviral drugs. Deep sequencing of the dengue virus 1 genome directly from the blood of twelve dengue patients and from mosquitoes given this blood showed consistent and distinct mutation patterns during infection. Importantly, we identified regions within the viral genome that are resistant to mutation during human and mosquito infection. Collectively, these findings provide fresh insights into potential antiviral targets and could serve as an approach to defining better regions for therapeutic targeting in other RNA viruses.

Sequence-Specific Fidelity Alterations Associated with West Nile Virus Attenuation in Mosquitoes

High rates of error-prone replication result in the rapid accumulation of genetic diversity of RNA viruses. Recent studies suggest that mutation rates are selected for optimal viral fitness and that modest variations in replicase fidelity may be associated with viral attenuation. Arthropod-borne viruses (arboviruses) are unique in their requirement for host cycling and may necessitate substantial genetic and phenotypic plasticity. In order to more thoroughly investigate the correlates, mechanisms and consequences of arbovirus fidelity, we selected fidelity variants of West Nile virus (WNV; Flaviviridae, Flavivirus) utilizing selection in the presence of a mutagen. We identified two mutations in the WNV RNA-dependent RNA polymerase associated with increased fidelity, V793I and G806R, and a single mutation in the WNV methyltransferase, T248I, associated with decreased fidelity. Both deep-sequencing and in vitro biochemical assays confirmed strain-specific differences in both fidelity and mutational bias. WNV fidelity variants demonstrated host-specific alterations to replicative fitness in vitro, with modest attenuation in mosquito but not vertebrate cell culture. Experimental infections of colonized and field populations of Cx. quinquefaciatus demonstrated that WNV fidelity alterations are associated with a significantly impaired capacity to establish viable infections in mosquitoes. Taken together, these studies (i) demonstrate the importance of allosteric interactions in regulating mutation rates, (ii) establish that mutational spectra can be both sequence and strain-dependent, and (iii) display the profound phenotypic consequences associated with altered replication complex function of flaviviruses.

West Nile virus (WNV) is the most geographically widespread arthropod-borne virus (arbovirus) in the world. Like most arboviruses, WNV is a RNA virus which is highly mutable and exists in nature as genetically diverse mutant swarms. Although many recent studies have investigated the relationship between virus mutation rate and viral fitness, this had not previously been determined for WNV or other flaviviruses. We identified WNV mutations associated with variation in mutation rate using cell culture passage in the presence of a mutagen and engineered these mutations into an infectious WNV clone in order to investigate the causes and consequences of altered fidelity. Our results demonstrate that interactions among proteins which comprise the WNV replication complex can significantly alter both the extent and types of mutations that occur. In addition, we show that both increasing and decreasing WNV fidelity has host-specific effects on replication in cell culture and is associated with nearly complete ablation of WNV infection in mosquito vectors. These results have significant implications for our understanding of arbovirus evolution, replication complex function and arboviral fitness in mosquitoes, and identify important targets to study the determinants and mechanisms of vector competence and arbovirus fidelity.

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

Within hosts, RNA viruses form populations that are genetically and phenotypically complex. Heterogeneity in RNA virus genomes arises due to error-prone replication and is reduced by stochastic and selective mechanisms that are incompletely understood. Defining how natural selection shapes RNA virus populations is critical because it can inform treatment paradigms and enhance control efforts. We allowed West Nile virus (WNV) to replicate in wild-caught American crows, house sparrows and American robins to assess how natural selection shapes RNA virus populations in ecologically relevant hosts that differ in susceptibility to virus-induced mortality. After five sequential passages in each bird species, we examined the phenotype and population diversity of WNV through fitness competition assays and next generation sequencing. We demonstrate that fitness gains occur in a species-specific manner, with the greatest replicative fitness gains in robin-passaged WNV and the least in WNV passaged in crows. Sequencing data revealed that intrahost WNV populations were strongly influenced by purifying selection and the overall complexity of the viral populations was similar among passaged hosts. However, the selective pressures that control WNV populations seem to be bird species-dependent. Specifically, crow-passaged WNV populations contained the most unique mutations (~1.7× more than sparrows, ~3.4× more than robins) and defective genomes (~1.4× greater than sparrows, ~2.7× greater than robins), but the lowest average mutation frequency (about equal to sparrows, ~2.6× lower than robins). Therefore, our data suggest that WNV replication in the most disease-susceptible bird species is positively associated with virus mutational tolerance, likely via complementation, and negatively associated with the strength of selection. These differences in genetic composition most likely have distinct phenotypic consequences for the virus populations. Taken together, these results reveal important insights into how different hosts may contribute to the emergence of RNA viruses.

Modulation of flavivirus population diversity by RNA interference

To test the hypothesis that RNA interference (RNAi) imposes diversifying selection on RNA virus genomes, we quantified West Nile virus (WNV) quasispecies diversity after passage in Drosophila cells in which RNAi was left intact, depleted, or stimulated against WNV. As predicted, WNV diversity was significantly lower in RNAi-depleted cells and significantly greater in RNAi-stimulated cells relative to that in controls. These findings reveal that an innate immune defense can shape viral population structure.

A positively selected mutation in the WNV 2K peptide confers resistance to superinfection exclusion in vivo

Molecular epidemiologic studies of North American (NA) West Nile virus (WNV; Flaviviridae, Flavivirus) have documented the displacement of the introduced NY99 genotype with WN02. In addition, these studies have shown that particular substitutions are under positive selection. One occurs in the C-terminus of the NS4A coding sequence and results in a valine to methionine substitution at position nine of the 2K peptide. 2K-V9M confers the ability to overcome superinfection exclusion in vitro. We hypothesized that WNV strains bearing 2K-V9M have higher fitness than wildtype in Culex quinquefasciatus mosquitoes. Although infection rates and viral titers were not significantly different, virus dissemination rates were significantly higher with WNV 2K-V9M. As a super-infecting virus, WNV 2K-V9M was more successful than wildtype, however, in a mixed infection, 2K-V9M was not. These data support observations that 2K-V9M confers a context-specific selective advantage in mosquitoes and provides an in vivo mechanism for its positive selection.

Chikungunya virus adaptation to a mosquito vector correlates with only few point mutations in the viral envelope glycoprotein

Like most arthropod-borne viruses (arboviruses), chikungunya virus (CHIKV) is a RNA virus maintained in nature in an alternating cycle of replication between invertebrate and vertebrate hosts. It has been assumed that host alternation restricts arbovirus genome evolution and imposes fitness trade-offs. Despite their slower rates of evolution, arboviruses still have the capacity to produce variants capable to exploit new environments. To test whether the evolution of the newly emerged epidemic variant of CHIKV (E1-226V) is constrained by host alternation, the virus was alternately-passaged in hamster-derived BHK-21 cells and Aedes aegypti-derived Aag-2 cells. It was also serially-passaged in BHK-21 or Aag-2 cells to promote adaptation to one cell type and presumably, fitness cost in the bypassed cell type. After 30 passages, obtained CHIKV strains were genetically and phenotypically characterized using in vitro and in vivo systems. Serially- and alternately-passaged strains can be distinguished by amino-acid substitutions in the E2 glycoprotein, responsible for receptor binding. Two substitutions at positions E2-64 and E2-208 only lower the dissemination of the variant E1-226V in Ae. aegypti. These amino-acid changes in the E2 glycoprotein might affect viral infectivity by altering the interaction between CHIKV E1-226V and the cellular receptor on the midgut epithelial cells in Ae. aegypti but not in Aedesalbopictus.

Vector competence of Culex pipiens quinquefasciatus (Diptera: Culicidae) for West Nile virus isolates from Florida

OBJECTIVES

To assess vector competence (infection, dissemination and transmission) of Culex pipiens quinquefasciatus for Florida (FL) West Nile virus (WNV) isolates.

METHODS

West Nile virus isolates (WN-FL-03: NY99 genotype; WN-FL-05-558, WN-FL-05-2186, WN-FL-05-510: WN02 genotype) collected from different regions of FL were used for vector competence experiments in Cx. p. quinquefasciatus from Alachua County and Indian River County in FL. Mosquitoes from both colonies were fed blood containing 7.9 ± 0.2 log10 plaque-forming units WNV/ml ± SE and incubated at 28 °C for 14 days. Vector competence, including rates of infection, dissemination, and transmission, was compared between colonies for WN-FL-03 using chi-squared. Virus titres in bodies, legs and saliva were compared using anova. Daily measurements of in vitro replication of WNV isolates were evaluated in Vero cells so that a standardised virus dose for each isolate could be delivered to mosquitoes.

RESULTS

Infection and dissemination rates were high (≥ 95%) and not affected by isolate or colony (infection, P = 0.679; dissemination, P = 0.799). Transmission rates were low (≤ 20%), detected in one colony and affected by isolate (P = 0.008). Body and leg titres differed between isolates (body titre, P = 0.031; leg titre, P = 0.044) and colonies (body titre, P = 0.001; leg titre, P = 0.013) while saliva titre did not differ between isolates (P = 0.462).

CONCLUSIONS

Variation in vector competence of mosquito populations may be attributed, in part, to exposures to WNV with genetic differences leading to different rates of replication in mosquitoes. Evaluation of vector competence for different WNV isolates may help us understand vector-virus interactions and, hence, the role of vectors in complex virus transmission cycles in nature.

Vector-virus interactions and transmission dynamics of West Nile virus

West Nile virus (WNV; Flavivirus; Flaviviridae) is the cause of the most widespread arthropod-borne viral disease in the world and the largest outbreak of neuroinvasive disease ever observed. Mosquito-borne outbreaks are influenced by intrinsic (e.g., vector and viral genetics, vector and host competence, vector life-history traits) and extrinsic (e.g., temperature, rainfall, human land use) factors that affect virus activity and mosquito biology in complex ways. The concept of vectorial capacity integrates these factors to address interactions of the virus with the arthropod host, leading to a clearer understanding of their complex interrelationships, how they affect transmission of vector-borne disease, and how they impact human health. Vertebrate factors including host competence, population dynamics, and immune status also affect transmission dynamics. The complexity of these interactions are further exacerbated by the fact that not only can divergent hosts differentially alter the virus, but the virus also can affect both vertebrate and invertebrate hosts in ways that significantly alter patterns of virus transmission. This chapter concentrates on selected components of the virus-vector-vertebrate interrelationship, focusing specifically on how interactions between vector, virus, and environment shape the patterns and intensity of WNV transmission.

Characterisation of divergent flavivirus NS3 and NS5 protein sequences detected in Rhipicephalus microplus ticks from Brazil

Transcripts similar to those that encode the nonstructural (NS) proteins NS3 and NS5 from flaviviruses were found in a salivary gland (SG) complementary DNA (cDNA) library from the cattle tick Rhipicephalus microplus. Tick extracts were cultured with cells to enable the isolation of viruses capable of replicating in cultured invertebrate and vertebrate cells. Deep sequencing of the viral RNA isolated from culture supernatants provided the complete coding sequences for the NS3 and NS5 proteins and their molecular characterisation confirmed similarity with the NS3 and NS5 sequences from other flaviviruses. Despite this similarity, phylogenetic analyses revealed that this potentially novel virus may be a highly divergent member of the genus Flavivirus. Interestingly, we detected the divergent NS3 and NS5 sequences in ticks collected from several dairy farms widely distributed throughout three regions of Brazil. This is the first report of flavivirus-like transcripts in R. microplus ticks. This novel virus is a potential arbovirus because it replicated in arthropod and mammalian cells; furthermore, it was detected in a cDNA library from tick SGs and therefore may be present in tick saliva. It is important to determine whether and by what means this potential virus is transmissible and to monitor the virus as a potential emerging tick-borne zoonotic pathogen.

Allele-specific qRT-PCR demonstrates superior detection of single nucleotide polymorphisms as genetic markers for West Nile virus compared to Luminex® and quantitative sequencing

To enable in vivo and in vitro competitive fitness comparisons among West Nile viruses (WNV), three reference viruses were marked genetically by site-directed mutagenesis with five synonymous nucleotide substitutions in the envelope gene region of the genome. Phenotypic neutrality of the mutants was assessed experimentally by competitive replication in cell culture and genetic stability of the substituted nucleotides was confirmed by direct sequencing. Luminex(®) technology, quantitative sequencing and quantitative RT-PCR (qRT-PCR) were compared in regard to specificity, sensitivity and accuracy for quantitation of wildtype and genetically marked viruses in mixed samples based on RNA obtained from samples of known viral titers. Although Luminex(®) technology and quantitative sequencing provided semi-quantitative or qualitative measurements, a sequence-specific primer extension approach using a specific reverse primer set in singleplex qRT-PCR demonstrated the best quantitation and specificity in the detection of RNA from wildtype and mutant viruses.

Factors shaping the adaptive landscape for arboviruses: implications for the emergence of disease

Many examples of the emergence or re-emergence of infectious diseases involve the adaptation of zoonotic viruses to new amplification hosts or to humans themselves. These include several instances of simple mutational adaptations, often to hosts closely related to the natural reservoirs. However, based on theoretical grounds, arthropod-borne viruses, or arboviruses, may face several challenges for adaptation to new hosts. Here, we review recent findings regarding adaptive evolution of arboviruses and its impact on disease emergence. We focus on the zoonotic alphaviruses Venezuelan equine encephalitis and chikungunya viruses, which have undergone adaptive evolution that mediated recent outbreaks of disease, as well as the flaviviruses dengue and West Nile viruses, which have emerged via less dramatic adaptive mechanisms.

The role of virulence in in vivo superinfection fitness of the vertebrate RNA virus infectious hematopoietic necrosis virus

We have developed a novel in vivo superinfection fitness assay to examine superinfection dynamics and the role of virulence in superinfection fitness. This assay involves controlled, sequential infections of a natural vertebrate host, Oncorhynchus mykiss (rainbow trout), with variants of a coevolved viral pathogen, infectious hematopoietic necrosis virus (IHNV). Intervals between infections ranged from 12 h to 7 days, and both frequency of superinfection and viral replication levels were examined. Using virus genotype pairs of equal and unequal virulence, we observed that superinfection generally occurred with decreasing frequency as the interval between exposures to each genotype increased. For both the equal-virulence and unequal-virulence genotype pairs, the frequency of superinfection in most cases was the same regardless of which genotype was used in the primary exposure. The ability to replicate in the context of superinfection also did not differ between the genotypes of equal or unequal virulence tested here. For both genotype pairs, the mean viral load of the secondary virus was significantly reduced in superinfection while primary virus replication was unaffected. Our results demonstrate, for the two genotype pairs examined, that superinfection restriction does occur for IHNV and that higher virulence did not correlate with a significant difference in superinfection fitness. To our knowledge, this is the first assay to examine the role of virulence of an RNA virus in determining superinfection fitness dynamics within a natural vertebrate host.

Viral fitness: definitions, measurement, and current insights

Viral fitness is an active area of research, with recent work involving an expanded number of human, non-human vertebrate, invertebrate, plant, and bacterial viruses. Many publications deal with RNA viruses associated with major disease emergence events, such as HIV-1, influenza virus, and Dengue virus. Study topics include drug resistance, immune escape, viral emergence, host jumps, mutation effects, quasispecies diversity, and mathematical models of viral fitness. Important recent trends include increasing use of in vivo systems to assess vertebrate virus fitness, and a broadening of research beyond replicative fitness to also investigate transmission fitness and epidemiologic fitness. This is essential for a more integrated understanding of overall viral fitness, with implications for disease management in the future.

Point mutations in the West Nile virus (Flaviviridae; Flavivirus) RNA-dependent RNA polymerase alter viral fitness in a host-dependent manner in vitro and in vivo

The West Nile virus (WNV) genome contains a single RNA-dependent RNA polymerase (RdRp) gene, which is responsible for replication of the viral genome and, as such, is an important target for antiviral therapy. Viral RdRps are known to lack proofreading capabilities and as a result viruses such as WNV exist as a mixture of viral genotypes within an infection, enabling the virus to readily emerge and adapt to new host environments. To test the consequences of subtle structural alterations remote from the RdRp active-site, the following single point mutations were engineered in the WNV NS5 RdRp coding region: T363N, A365N, and T537I; these mutations were selected in an effort to stabilize the secondary structural elements near the rNTP binding pocket of the RdRp. Mutant viruses were tested in vitro on Vero, C6/36, Culex tarsalis and DF-1 cell types and in vivo in one day old chickens and Culex pipiens mosquitoes. Plaque morphology was affected by each mutation and growth and RNA replication kinetics were altered as well. Our results demonstrate that subtle alteration of the RdRp protein away from the active site can have a significant overall biological effect on WNV fitness, and that this effect can be host-dependent.

Quantification of intrahost bottlenecks of West Nile virus in Culex pipiens mosquitoes using an artificial mutant swarm

Mosquito-borne viruses are predominantly RNA viruses which exist within hosts as diverse mutant swarms. Defining the way in which stochastic forces within mosquito vectors shape these swarms is critical to advancing our understanding of the evolutionary and adaptive potential of these pathogens. There are multiple barriers within a mosquito which a viral swarm must traverse in order to ultimately be transmitted. Here, using artificial mutant swarms composed of neutral variants of West Nile virus (WNV), we tracked changes to swarm breadth over time and space in Culex pipiens mosquitoes. Our results demonstrate that all variants have the potential to survive intrahost bottlenecks, yet mean swarm breadth decreases during both midgut infection and transmission when starting populations contain higher levels of minority variants. In addition, WNV swarms are subject to temporal sweeps which act to significantly decrease intrahost diversity over time. Taken together, these data demonstrate the profound effects that stochastic forces can have in shaping arboviral mutant swarms.

Specific and nonspecific host adaptation during arboviral experimental evolution

During the past decade or so, there has been a substantial body of work to dissect arboviral evolution and to develop models of adaptation during host switching. Regardless of what species serve as host or vectors, and of the geographic distribution and the mechanisms of replication, arboviruses tend to have slow evolutionary rates in nature. The hypothesis that this is the result of replication in the disparate environments provided by host and vector did not receive solid experimental support in any of the many viral species tested. Instead, it seems that from the virus's point of view, either the two environments are sufficiently similar or one of the environments so dominates viral evolution that there is tolerance for suboptimal adaptation to the other environment. Replication in alternating environments has an unexpected cost in that there is decreased genetic variance that translates into a compromised adaptability for bypassed environments. Arboviruses under strong and continuous positive selection may have unusual patterns of genomic changes, with few or no mutations accumulated in the consensus sequence or with dN/dS values typically consistent with random drift in DNA-based organisms.

West Nile virus experimental evolution in vivo and the trade-off hypothesis

In nature, arthropod-borne viruses (arboviruses) perpetuate through alternating replication in vertebrate and invertebrate hosts. The trade-off hypothesis proposes that these viruses maintain adequate replicative fitness in two disparate hosts in exchange for superior fitness in one host. Releasing the virus from the constraints of a two-host cycle should thus facilitate adaptation to a single host. This theory has been addressed in a variety of systems, but remains poorly understood. We sought to determine the fitness implications of alternating host replication for West Nile virus (WNV) using an in vivo model system. Previously, WNV was serially or alternately passed 20 times in vivo in chicks or mosquitoes and resulting viruses were characterized genetically. In this study, these test viruses were competed in vivo in fitness assays against an unpassed marked reference virus. Fitness was assayed in chicks and in two important WNV vectors, Culex pipiens and Culex quinquefasciatus. Chick-specialized virus displayed clear fitness gains in chicks and in Cx. pipiens but not in Cx. quinquefasciatus. Cx. pipiens-specialized virus experienced reduced fitness in chicks and little change in either mosquito species. These data suggest that when fitness is measured in birds the trade-off hypothesis is supported; but in mosquitoes it is not. Overall, these results suggest that WNV evolution is driven by alternate cycles of genetic expansion in mosquitoes, where purifying selection is weak and genetic diversity generated, and restriction in birds, where purifying selection is strong.

Nonconsensus West Nile virus genomes arising during mosquito infection suppress pathogenesis and modulate virus fitness in vivo

West Nile virus (WNV) is similar to other RNA viruses in that it forms genetically complex populations within hosts. The virus is maintained in nature in mosquitoes and birds, with each host type exerting distinct influences on virus populations. We previously observed that prolonged replication in mosquitoes led to increases in WNV genetic diversity and diminished pathogenesis in mice without remarkable changes to the consensus genome sequence. We therefore sought to evaluate the relationships between individual and group phenotypes in WNV and to discover novel viral determinants of pathogenesis in mice and fitness in mosquitoes and birds. Individual plaque size variants were isolated from a genetically complex population, and mutations conferring a small-plaque and mouse-attenuated phenotype were localized to the RNA helicase domain of the NS3 protein by reverse genetics. The mutation, an Asp deletion, did not alter type I interferon production in the host but rendered mutant viruses more susceptible to interferon compared to wild type (WT) WNV. Finally, we used an in vivo fitness assay in Culex quinquefasciatus mosquitoes and chickens to determine whether the mutation in NS3 influenced fitness. The fitness of the NS3 mutant was dramatically lower in chickens and moderately lower in mosquitoes, indicating that RNA helicase is a major fitness determinant of WNV and that the effect on fitness is host specific. Overall, this work highlights the complex relationships that exist between individual and group phenotypes in RNA viruses and identifies RNA helicase as an attenuation and fitness determinant in WNV.

West Nile virus genetic diversity is maintained during transmission by Culex pipiens quinquefasciatus mosquitoes

Due to error-prone replication, RNA viruses exist within hosts as a heterogeneous population of non-identical, but related viral variants. These populations may undergo bottlenecks during transmission that stochastically reduce variability leading to fitness declines. Such bottlenecks have been documented for several single-host RNA viruses, but their role in the population biology of obligate two-host viruses such as arthropod-borne viruses (arboviruses) in vivo is unclear, but of central importance in understanding arbovirus persistence and emergence. Therefore, we tracked the composition of West Nile virus (WNV; Flaviviridae, Flavivirus) populations during infection of the vector mosquito, Culex pipiens quinquefasciatus to determine whether WNV populations undergo bottlenecks during transmission by this host. Quantitative, qualitative and phylogenetic analyses of WNV sequences in mosquito midguts, hemolymph and saliva failed to document reductions in genetic diversity during mosquito infection. Further, migration analysis of individual viral variants revealed that while there was some evidence of compartmentalization, anatomical barriers do not impose genetic bottlenecks on WNV populations. Together, these data suggest that the complexity of WNV populations are not significantly diminished during the extrinsic incubation period of mosquitoes.

Insights into arbovirus evolution and adaptation from experimental studies

Arthropod-borne viruses (arboviruses) are maintained in nature by cycling between vertebrate hosts and haematophagous invertebrate vectors. These viruses are responsible for causing a significant public health burden throughout the world, with over 100 species having the capacity to cause human disease. Arbovirus outbreaks in previously naïve environments demonstrate the potential of these pathogens for expansion and emergence, possibly exacerbated more recently by changing climates. These recent outbreaks, together with the continued devastation caused by endemic viruses, such as Dengue virus which persists in many areas, demonstrate the need to better understand the selective pressures that shape arbovirus evolution. Specifically, a comprehensive understanding of host-virus interactions and how they shape both host-specific and virus-specific evolutionary pressures is needed to fully evaluate the factors that govern the potential for host shifts and geographic expansions. One approach to advance our understanding of the factors influencing arbovirus evolution in nature is the use of experimental studies in the laboratory. Here, we review the contributions that laboratory passage and experimental infection studies have made to the field of arbovirus adaptation and evolution, and how these studies contribute to the overall field of arbovirus evolution. In particular, this review focuses on the areas of evolutionary constraints and mutant swarm dynamics; how experimental results compare to theoretical predictions; the importance of arbovirus ecology in shaping viral swarms; and how current knowledge should guide future questions relevant to understanding arbovirus evolution.

Host alternation of chikungunya virus increases fitness while restricting population diversity and adaptability to novel selective pressures

The mechanisms by which RNA arboviruses, including chikungunya virus (CHIKV), evolve and maintain the ability to infect vertebrate and invertebrate hosts are poorly understood. To understand how host specificity shapes arbovirus populations, we studied CHIKV populations passaged alternately between invertebrate and vertebrate cells (invertebrate ↔ vertebrate) to simulate natural alternation and contrasted the results with those for populations that were artificially released from cycling by passage in single cell types. These CHIKV populations were characterized by measuring genetic diversity, changes in fitness, and adaptability to novel selective pressures. The greatest fitness increases were observed in alternately passaged CHIKV, without drastic changes in population diversity. The greatest increases in genetic diversity were observed after serial passage and correlated with greater adaptability. These results suggest an evolutionary trade-off between maintaining fitness for invertebrate ↔ vertebrate cell cycling, where maximum adaptability is possible only via enhanced population diversity and extensive exploration of sequence space.