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

Adaptive evolution of West Nile virus facilitated increased transmissibility and prevalence in New York State

West Nile virus (WNV; Flavivirus, Flaviviridae) was introduced to New York State (NYS) in 1999 and rapidly expanded its range through the continental United States (US). Apart from the displacement of the introductory NY99 genotype with the WN02 genotype, there has been little evidence of adaptive evolution of WNV in the US. WNV NY10, characterized by shared amino acid substitutions R1331K and I2513M, emerged in 2010 coincident with increased WNV cases in humans and prevalence in mosquitoes. Previous studies demonstrated an increase in frequency of NY10 strains in NYS and evidence of positive selection. Here, we present updated surveillance and sequencing data for WNV in NYS and investigate if NY10 genotype strains are associated with phenotypic change consistent with an adaptive advantage. Results confirm a significant increase in prevalence in mosquitoes though 2018, and updated sequencing demonstrates a continued dominance of NY10. We evaluated NY10 strains in Culex pipiens mosquitoes to assess vector competence and found that the NY10 genotype is associated with both increased infectivity and transmissibility. Experimental infection of American robins (Turdus migratorius) was additionally completed to assess viremia kinetics of NY10 relative to WN02. Modelling the increased infectivity and transmissibility of the NY10 strains together with strain-specific viremia demonstrates a mechanistic basis for selection that has likely contributed to the increased prevalence of WNV in NYS.

Deciphering the Tissue Tropism of the RNA Viromes Harbored by Field-Collected Anopheles sinensis and Culex quinquefasciatus

Arboviruses and insect-specific viruses (ISVs) are two major types of viruses harbored by mosquitoes that are distinguished by the involvement of vertebrate hosts in their transmission cycles. While intensive studies have focused on the transmission, tissue tropism, and evolution of arboviruses, these characteristics are poorly investigated in ISVs, which dominate the mosquito virome. Therefore, in this study, we collected two mosquito species, Anopheles sinensis and Culex quinquefasciatus, in the field and used a metatranscriptomics approach to characterize their RNA viromes in different tissues, such as the midgut, legs, salivary gland, eggs, and the remainder of the carcass. Blood-engorged individuals of these species were captured in 3 locations, and 60 mosquitoes were pooled from each species and location. A total of 40 viral species from diverse viral taxa associated with all viral RNA genome types were identified, among which 19 were newly identified in this study. According to the current viral taxonomy, some of these viruses, such as Yancheng Anopheles associated virus 2 (Narnaviridae) and Jiangsu Anopheles-related virus (Ghabrivirales), were novel. The two investigated mosquito species generally harbored distinct viromes. Nevertheless, the viruses were generally shared among different tissue types to various degrees. Specifically, the eggs possessed a viral community with significantly lower diversity and abundance than those in other tissues, whereas the legs and salivary glands exhibited higher viral abundance. The compositions and distributions of the viromes of different mosquito tissues were demonstrated for the first time in our study, providing important insight into the virome dynamics within individual mosquitoes. IMPORTANCE ISVs are considered to be ancestral to arboviruses. Because of their medical importance, arboviruses have been well studied from the aspects of their transmission mode, evolution of dual-host tropism, and genetic dynamics within mosquito vectors. However, the mode of ISV maintenance is poorly understood, even though many novel ISVs have been identified with the emergence of sequencing technology. In our study, in addition to the identification of a diverse virus community, the tissue tropism of RNA viromes harbored by two field-collected mosquito species was demonstrated for the first time. According to the results, the virus communities of different tissues, such as the salivary glands, midguts, legs, and eggs, can help us understand the evolution, transmission routes, and maintenance modes of mosquito-specific viruses in nature.

Fatal Human Infection with Evidence of Intrahost Variation of Eastern Equine Encephalitis Virus, Alabama, USA, 2019

Eastern equine encephalitis virus (EEEV) is an arbovirus in the family Togaviridae, genus Alphavirus, found in North America and associated with freshwater/hardwood swamps in the Atlantic, Gulf Coast, and Great Lakes regions. EEEV disease in humans is rare but causes substantial illness and death. To investigate the molecular epidemiology and microevolution of EEEV from a fatal case in Alabama, USA, in 2019, we used next-generation sequencing of serum and cerebrospinal fluid (CSF). Phylogenetic inference indicated that the infecting strain may be closely related to isolates from Florida detected during 2010-2014, suggesting potential seeding from Florida. EEEV detected in serum displayed a higher degree of variability with more single-nucleotide variants than that detected in the CSF. These data refine our knowledge of EEEV molecular epidemiologic dynamics in the Gulf Coast region and demonstrate potential quasispecies bottlenecking within the central nervous system of a human host.

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.

Specific capture and whole-genome phylogeography of Dolphin morbillivirus

Dolphin morbillivirus (DMV) is considered an emerging threat having caused several epidemics worldwide. Only few DMV genomes are publicly available. Here, we report the use of target enrichment directly from cetacean tissues to obtain novel DMV genome sequences, with sequence comparison and phylodynamic analysis. RNA from 15 tissue samples of cetaceans stranded along the Italian and French coasts (2008-2017) was purified and processed using custom probes (by bait hybridization) for target enrichment and sequenced on Illumina MiSeq. Data were mapped against the reference genome, and the novel sequences were aligned to the available genome sequences. The alignment was then used for phylogenetic and phylogeographic analysis using MrBayes and BEAST. We herein report that target enrichment by specific capture may be a successful strategy for whole-genome sequencing of DMV directly from field samples. By this strategy, 14 complete and one partially complete genomes were obtained, with reads mapping to the virus up to 98% and coverage up to 7800X. The phylogenetic tree well discriminated the Mediterranean and the NE-Atlantic strains, circulating in the Mediterranean Sea and causing two different epidemics (2008-2015 and 2014-2017, respectively), with a limited time overlap of the two strains, sharing a common ancestor approximately in 1998.

Divergent Mutational Landscapes of Consensus and Minority Genotypes of West Nile Virus Demonstrate Host and Gene-Specific Evolutionary Pressures

Our current understanding of the natural evolution of RNA viruses comes largely from consensus level genetic analyses which ignore the diverse mutant swarms that comprise within-host viral populations. The breadth and composition of viral mutant swarms impact viral fitness and adaptation, and the capacity for swarm plasticity is likely to be particularly important for arthropod-borne viruses (arboviruses) that cycle between taxonomically divergent hosts. Despite this, characterization of the relationship between the selective pressures and genetic signatures of the mutant swarm and consensus sequences is lacking. To clarify this, we analyzed previously generated whole genome, deep-sequencing data from 548 West Nile virus samples isolated from avian tissues or mosquitoes in New York State from 1999-2018. Both consensus level (interhost) and minority level (intrahost) nucleotide and amino acid sequences were analyzed, and diversity at each position was calculated across the genome in order to assess the relationship between minority and consensus sequences for individual genes and hosts. Our results indicate that consensus sequences are an inept representation of the overall genetic diversity. Unique host and gene-specific signatures and selective pressures were identified. These data demonstrate that an accurate and comprehensive understanding of arbovirus evolution and adaptation within and between hosts requires consideration of minority genotypes.

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.

Introduction, Spread, and Establishment of West Nile Virus in the Americas

The introduction of West Nile virus (WNV) to North America in 1999 and its subsequent rapid spread across the Americas demonstrated the potential impact of arboviral introductions to new regions, and this was reinforced by the subsequent introductions of chikungunya and Zika viruses. Extensive studies of host-pathogen-vector-environment interactions over the past two decades have illuminated many aspects of the ecology and evolution of WNV and other arboviruses, including the potential for pathogen adaptation to hosts and vectors, the influence of climate, land use and host immunity on transmission ecology, and the difficulty in preventing the establishment of a zoonotic pathogen with abundant wildlife reservoirs. Here, we focus on outstanding questions concerning the introduction, spread, and establishment of WNV in the Americas, and what it can teach us about the future of arboviral introductions. Key gaps in our knowledge include the following: viral adaptation and coevolution of hosts, vectors and the virus; the mechanisms and species involved in the large-scale spatial spread of WNV; how weather modulates WNV transmission; the drivers of large-scale variation in enzootic transmission; the ecology of WNV transmission in Latin America; and the relative roles of each component of host-virus-vector interactions in spatial and temporal variation in WNV transmission. Integrative studies that examine multiple factors and mechanisms simultaneously are needed to advance our knowledge of mechanisms driving transmission.

Many human RNA viruses show extraordinarily stringent selective constraints on protein evolution

How negative selection, positive selection, and population size contribute to the large variation in nucleotide substitution rates among RNA viruses remains unclear. Here, we studied the ratios of nonsynonymous-to-synonymous substitution rates (d _N/d _S) in protein-coding genes of human RNA and DNA viruses and mammals. Among the 21 RNA viruses studied, 18 showed a genome-average d _N/d _S from 0.01 to 0.10, indicating that over 90% of nonsynonymous mutations are eliminated by negative selection. Only HIV-1 showed a d _N/d _S (0.31) higher than that (0.22) in mammalian genes. By comparing the d _N/d _S values among genes in the same genome and among species or strains, we found that both positive selection and population size play significant roles in the d _N/d _S variation among genes and species. Indeed, even in flaviviruses and picornaviruses, which showed the lowest ratios among the 21 species studied, positive selection appears to have contributed significantly to d _N/d _S We found the view that positive selection occurs much more frequently in influenza A subtype H3N2 than subtype H1N1 holds only for the hemagglutinin and neuraminidase genes, but not for other genes. Moreover, we found no support for the view that vector-borne RNA viruses have lower d _N/d _S ratios than non-vector-borne viruses. In addition, we found a correlation between d _N and d _S, implying a correlation between d _N and the mutation rate. Interestingly, only 2 of the 8 DNA viruses studied showed a d _N/d _S < 0.10, while 4 showed a d _N/d _S > 0.22. These observations increase our understanding of the mechanisms of RNA virus evolution.

Evolutionary dynamics and molecular epidemiology of West Nile virus in New York State: 1999-2015

Following its introduction into New York State (NYS) in 1999, West Nile virus (WNV; Flavivirus, Flaviviridae) underwent a rapid expansion throughout the USA and into Canada and Latin America. WNV has been characterized as being evolutionarily stable, with weak geographic structure, a dominance of purifying selection and limited adaptive change. We analyzed all available full-genome WNV sequences, focusing on the 543 available sequences from NYS, which included 495 newly sequenced 2000–15 isolates. In addition, we analyzed deep-sequencing data from 317 of these isolates. While our data are generally in agreement with the limited pace of evolutionary change and broad geographic and temporal mixing identified in other studies, we have identified some important exceptions. Most notably, there are 14 codons which demonstrated evidence of positive selection as determined by multiple models, including some positions with evidence of selection in NYS exclusively. Coincident with increased WNV activity, genotypes possessing one or more of these mutations, designated NY01, NY07, and NY10, have increased in prevalence in recent years and displaced historic strains. In addition, we have found a geographical bias with many of these mutations, which suggests selective pressures and adaptations could be regional. Lastly, our deep-sequencing data suggest both increased overall diversity in avian tissue isolates relative to mosquito isolates and multiple non-synonymous minority variants that are both host-specific and retained over time and space. Together, these data provide novel insight into the evolutionary pressures on WNV and the need for continued genetic surveillance and characterization of emergent strains.

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.

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.