Evolutionary dynamics of West Nile virus in the United States, 1999-2011: phylogeny, selection pressure and evolutionary time-scale analysis

Duck Tembusu virus in North Vietnam: epidemiological and genetic analysis reveals novel virus strains

Tembusu virus (TMUV) is an important infectious disease, causing economic losses in duck production. Since the first report of TMUV infection in Vietnam in 2020, the disease has persisted and affected poultry production in the country. This study conducted epidemiological and genetic characterization of the viral strains circulating in north Vietnam based on 130 pooled tissue samples collected in six provinces/cities during 2021. The TMUV genome was examined using conventional PCR. The results indicated that 21 (16.15%) samples and 9 (23.68%) farms were positive for the viral genome. The positive rate was 59.26% for ducks at ages 2-4 weeks, which was significantly higher than for ducks at ages >4 weeks and < 2 weeks. Genetic analysis of the partial envelope gene (891 bp) sequences indicated that the five Vietnamese TMUVs shared 99.55-100% nucleotide identity, while the rates were in the range 99.59-100% based on the pre-membrane gene sequences (498 bp). The five Vietnamese TMUV strains obtained formed a novel single subcluster. These strains were closely related to Chinese strains and differed from the vaccine strain, suggesting that Vietnamese TMUV strains were field viruses. It needs to be further studied on vaccine development to prevent effects of TMUV infection on poultry production across Vietnam.

Many-core algorithms for high-dimensional gradients on phylogenetic trees

MOTIVATION

Advancements in high-throughput genomic sequencing are delivering genomic pathogen data at an unprecedented rate, positioning statistical phylogenetics as a critical tool to monitor infectious diseases globally. This rapid growth spurs the need for efficient inference techniques, such as Hamiltonian Monte Carlo (HMC) in a Bayesian framework, to estimate parameters of these phylogenetic models where the dimensions of the parameters increase with the number of sequences N. HMC requires repeated calculation of the gradient of the data log-likelihood with respect to (wrt) all branch-length-specific (BLS) parameters that traditionally takes O(N2) operations using the standard pruning algorithm. A recent study proposes an approach to calculate this gradient in O(N), enabling researchers to take advantage of gradient-based samplers such as HMC. The CPU implementation of this approach makes the calculation of the gradient computationally tractable for nucleotide-based models but falls short in performance for larger state-space size models, such as Markov-modulated and codon models. Here, we describe novel massively parallel algorithms to calculate the gradient of the log-likelihood wrt all BLS parameters that take advantage of graphics processing units (GPUs) and result in many fold higher speedups over previous CPU implementations.

RESULTS

We benchmark these GPU algorithms on three computing systems using three evolutionary inference examples exploring complete genomes from 997 dengue viruses, 62 carnivore mitochondria and 49 yeasts, and observe a >128-fold speedup over the CPU implementation for codon-based models and >8-fold speedup for nucleotide-based models. As a practical demonstration, we also estimate the timing of the first introduction of West Nile virus into the continental Unites States under a codon model with a relaxed molecular clock from 104 full viral genomes, an inference task previously intractable.

AVAILABILITY AND IMPLEMENTATION

We provide an implementation of our GPU algorithms in BEAGLE v4.0.0 (https://github.com/beagle-dev/beagle-lib), an open-source library for statistical phylogenetics that enables parallel calculations on multi-core CPUs and GPUs. We employ a BEAGLE-implementation using the Bayesian phylogenetics framework BEAST (https://github.com/beast-dev/beast-mcmc).

Recent Population Dynamics of Japanese Encephalitis Virus

Japanese encephalitis virus (JEV) causes acute viral encephalitis in humans and reproductive disorders in pigs. JEV emerged during the 1870s in Japan, and since that time, JEV has been transmitted exclusively throughout Asia, according to known reporting and sequencing records. A recent JEV outbreak occurred in Australia, affecting commercial piggeries across different temperate southern Australian states, and causing confirmed infections in humans. A total of 47 human cases and 7 deaths were reported. The recent evolving situation of JEV needs to be reported due to its continuous circulation in endemic regions and spread to non-endemics areas. Here, we reconstructed the phylogeny and population dynamics of JEV using recent JEV isolates for the future perception of disease spread. Phylogenetic analysis shows the most recent common ancestor occurred about 2993 years ago (YA) (95% Highest posterior density (HPD), 2433 to 3569). Our results of the Bayesian skyline plot (BSP) demonstrates that JEV demography lacks fluctuations for the last two decades, but it shows that JEV genetic diversity has increased during the last ten years. This indicates the potential JEV replication in the reservoir host, which is helping it to maintain its genetic diversity and to continue its dispersal into non-endemic areas. The continuous spread in Asia and recent detection from Australia further support these findings. Therefore, an enhanced surveillance system is needed along with precautionary measures such as regular vaccination and mosquito control to avoid future JEV outbreaks.

Re-Emerged Genotype IV of Japanese Encephalitis Virus Is the Youngest Virus in Evolution

An outbreak of viral encephalitis caused by a Japanese encephalitis virus (JEV) genotype IV infection occurred in Australia between 2021 and 2022. A total of 47 cases and seven deaths were reported as of November 2022. This is the first outbreak of human viral encephalitis caused by JEV GIV since it was first isolated in Indonesia in the late 1970s. Here, a comprehensive phylogenetic analysis based on the whole genome sequences of JEVs revealed it emerged 1037 years ago (95% HPD: 463 to 2100 years). The evolutionary order of JEV genotypes is as follows: GV, GIII, GII, GI, and GIV. The JEV GIV emerged 122 years ago (95% HPD: 57-233) and is the youngest viral lineage. The mean substitution rate of the JEV GIV lineage was 1.145 × 10-3 (95% HPD values, 9.55 × 10-4, 1.35 × 10-3), belonging to rapidly evolving viruses. A series of amino acid mutations with the changes of physico-chemical properties located in the functional important domains within the core and E proteins distinguished emerging GIV isolates from old ones. These results demonstrate the JEV GIV is the youngest JEV genotype at a rapid evolution stage and has good host/vector adaptability for introduction to non-endemic areas. Thus, surveillance of JEVs is highly recommended.

Development of Toehold Switches as a Novel Ribodiagnostic Method for West Nile Virus

West Nile virus (WNV) is an emerging neurotropic RNA virus and a member of the genus Flavivirus. Naturally, the virus is maintained in an enzootic cycle involving mosquitoes as vectors and birds that are the principal amplifying virus hosts. In humans, the incubation period for WNV disease ranges from 3 to 14 days, with an estimated 80% of infected persons being asymptomatic, around 19% developing a mild febrile infection and less than 1% developing neuroinvasive disease. Laboratory diagnosis of WNV infection is generally accomplished by cross-reacting serological methods or highly sensitive yet expensive molecular approaches. Therefore, current diagnostic tools hinder widespread surveillance of WNV in birds and mosquitoes that serve as viral reservoirs for infecting secondary hosts, such as humans and equines. We have developed a synthetic biology-based method for sensitive and low-cost detection of WNV. This method relies on toehold riboswitches designed to detect WNV genomic RNA as transcriptional input and process it to GFP fluorescence as translational output. Our methodology offers a non-invasive tool with reduced operating cost and high diagnostic value that can be used for field surveillance of WNV in humans as well as in bird and mosquito populations.

An advanced sequence clustering and designation workflow reveals the enzootic maintenance of a dominant West Nile virus subclade in Germany

West Nile virus (WNV) is the most widespread arthropod-borne (arbo) virus and the primary cause of arboviral encephalitis globally. Members of WNV species genetically diverged and are classified into different hierarchical groups below species rank. However, the demarcation criteria for allocating WNV sequences into these groups remain individual and inconsistent, and the use of names for different levels of the hierarchical levels is unstructured. In order to have an objective and comprehensible grouping of WNV sequences, we developed an advanced grouping workflow using the 'affinity propagation clustering' algorithm and newly included the 'agglomerative hierarchical clustering' algorithm for the allocation of WNV sequences into different groups below species rank. In addition, we propose to use a fixed set of terms for the hierarchical naming of WNV below species level and a clear decimal numbering system to label the determined groups. For validation, we applied the refined workflow to WNV sequences that have been previously grouped into various lineages, clades, and clusters in other studies. Although our workflow regrouped some WNV sequences, overall, it generally corresponds with previous groupings. We employed our novel approach to the sequences from the WNV circulation in Germany 2020, primarily from WNV-infected birds and horses. Besides two newly defined minor (sub)clusters comprising only three sequences each, Subcluster 2.5.3.4.3c was the predominant WNV sequence group detected in Germany from 2018 to 2020. This predominant subcluster was also associated with at least five human WNV infections in 2019-20. In summary, our analyses imply that the genetic diversity of the WNV population in Germany is shaped by enzootic maintenance of the dominant WNV subcluster accompanied by sporadic incursions of other rare clusters and subclusters. Moreover, we show that our refined approach for sequence grouping yields meaningful results. Although we primarily aimed at a more detailed WNV classification, the presented workflow can also be applied to the objective genotyping of other virus species.

The first three waves of the Covid-19 pandemic hint at a limited genetic repertoire for SARS-CoV-2

The genomic diversity of SARS-CoV-2 is the result of a relatively low level of spontaneous mutations introduced during viral replication. With millions of SARS-CoV-2 genome sequences now available, we can begin to assess the overall genetic repertoire of this virus. We find that during 2020, there was a global wave of one variant that went largely unnoticed, possibly because its members were divided over several sublineages (B.1.177 and sublineages B.1.177.XX). We collectively call this Janus, and it was eventually replaced by the Alpha (B.1.1.7) variant of concern (VoC), next replaced by Delta (B.1.617.2), which itself might soon be replaced by a fourth pandemic wave consisting of Omicron (B.1.1.529). We observe that splitting up and redefining variant lineages over time, as was the case with Janus and is now happening with Alpha, Delta and Omicron, is not helpful to describe the epidemic waves spreading globally. Only ∼5% of the 30 000 nucleotides of the SARS-CoV-2 genome are found to be variable. We conclude that a fourth wave of the pandemic with the Omicron variant might not be that different from other VoCs, and that we may already have the tools in hand to effectively deal with this new VoC.

The central role of Italy in the spatial spread of USUTU virus in Europe

USUTU virus (USUV) is an arbovirus maintained in the environment through a bird-mosquito enzootic cycle. Previous surveillance plans highlighted the endemicity of USUV in North-eastern Italy. In this work, we sequenced 138 new USUV full genomes from mosquito pools (Culex pipiens) and wild birds collected in North-eastern Italy and we investigated the evolutionary processes (phylogenetic analysis, selection pressure and evolutionary time-scale analysis) and spatial spread of USUV strains circulating in the European context and in Italy, with a particular focus on North-eastern Italy. Our results confirmed the circulation of viruses belonging to four different lineages in Italy (EU1, EU2, EU3 and EU4), with the newly sequenced viruses from the North-eastern regions, Veneto and Friuli Venezia Giulia, belonging to the EU2 lineage and clustering into two different sub-lineages, EU2-A and EU2-B. Specific mutations characterize each European lineage and geographic location seem to have shaped their phylogenetic structure. By investigating the spatial spread in Europe, we were able to show that Italy acted mainly as donor of USUV to neighbouring countries. At a national level, we identified two geographical clusters mainly circulating in Northern and North-western Italy, spreading both northward and southward. Our analyses provide important information on the spatial and evolutionary dynamics of USUTU virus that can help to improve surveillance plans and control strategies for this virus of increasing concern for human health.

High genetic variability of Schmallenberg virus M-segment leads to efficient immune escape from neutralizing antibodies

Schmallenberg virus (SBV) is the cause of severe fetal malformations when immunologically naïve pregnant ruminants are infected. In those malformed fetuses, a "hot-spot"-region of high genetic variability within the N-terminal region of the viral envelope protein Gc has been observed previously, and this region co-localizes with a known key immunogenic domain. We studied a series of M-segments of those SBV variants from malformed fetuses with point mutations, insertions or large in-frame deletions of up to 612 nucleotides. Furthermore, a unique cell-culture isolate from a malformed fetus with large in-frame deletions within the M-segment was analyzed. Each Gc-protein with amino acid deletions within the "hot spot" of mutations failed to react with any neutralizing anti-SBV monoclonal antibodies or a domain specific antiserum. In addition, in vitro virus replication of the natural deletion variant could not be markedly reduced by neutralizing monoclonal antibodies or antisera from the field. The large-deletion variant of SBV that could be isolated in cell culture was highly attenuated with an impaired in vivo replication following the inoculation of sheep. In conclusion, the observed amino acid sequence mutations within the N-terminal main immunogenic domain of glycoprotein Gc result in an efficient immune evasion from neutralizing antibodies in the special environment of a developing fetus. These SBV-variants were never detected as circulating viruses, and therefore should be considered to be dead-end virus variants, which are not able to spread further. The observations described here may be transferred to other orthobunyaviruses, particularly those of the Simbu serogroup that have been shown to infect fetuses. Importantly, such mutant strains should not be included in attempts to trace the spatial-temporal evolution of orthobunyaviruses in molecular-epidemiolocal approaches during outbreak investigations.

Epidemiological hypothesis testing using a phylogeographic and phylodynamic framework

Computational analyses of pathogen genomes are increasingly used to unravel the dispersal history and transmission dynamics of epidemics. Here, we show how to go beyond historical reconstructions and use spatially-explicit phylogeographic and phylodynamic approaches to formally test epidemiological hypotheses. We illustrate our approach by focusing on the West Nile virus (WNV) spread in North America that has substantially impacted public, veterinary, and wildlife health. We apply an analytical workflow to a comprehensive WNV genome collection to test the impact of environmental factors on the dispersal of viral lineages and on viral population genetic diversity through time. We find that WNV lineages tend to disperse faster in areas with higher temperatures and we identify temporal variation in temperature as a main predictor of viral genetic diversity through time. By contrasting inference with simulation, we find no evidence for viral lineages to preferentially circulate within the same migratory bird flyway, suggesting a substantial role for non-migratory birds or mosquito dispersal along the longitudinal gradient.

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 Epidemic in Germany Triggered by Epizootic Emergence, 2019

One year after the first autochthonous transmission of West Nile virus (WNV) to birds and horses in Germany, an epizootic emergence of WNV was again observed in 2019. The number of infected birds and horses was considerably higher compared to 2018 (12 birds, two horses), resulting in the observation of the first WNV epidemy in Germany: 76 cases in birds, 36 in horses and five confirmed mosquito-borne, autochthonous human cases. We demonstrated that Germany experienced several WNV introduction events and that strains of a distinct group (Eastern German WNV clade), which was introduced to Germany as a single introduction event, dominated mosquito, birds, horse and human-related virus variants in 2018 and 2019. Virus strains in this clade are characterized by a specific-Lys2114Arg mutation, which might lead to an increase in viral fitness. Extraordinary high temperatures in 2018/2019 allowed a low extrinsic incubation period (EIP), which drove the epizootic emergence and, in the end, most likely triggered the 2019 epidemic. Spatiotemporal EIP values correlated with the geographical WNV incidence. This study highlights the risk of a further spread in Germany in the next years with additional human WNV infections. Thus, surveillance of birds is essential to provide an early epidemic warning and thus, initiate targeted control measures.

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.

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.

Phylogenetic analysis of West Nile virus in Quebec, Canada, 2004-2016: Co-circulation of distinct variants harbouring conserved amino acid motifs in North America

West Nile virus (WNV) was introduced for the first time in the western hemisphere in 1999 in New York City. In 2002, a phenotype-modifying mutation (Env-V159A) defined the first North American genotype WN02. So far, three genotypes has been described in North America but little is known about WNV evolution in Canada. We report the phylogenetic characterization of twenty-six WNV genomes isolated from mosquitoes in the province of Quebec. WNV strains found in Quebec are phylogenetically related to American strains collected in northern and southern regions. We also noted the presence of two robust monophyletic groups of isolates characterized by distinct conserved amino acid motifs. These emerging genotypes were detected for several years in different ecosystems. These results highlight the need for the maintenance of a nationwide surveillance to follow the dispersion of emergent WNV genotypes.

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.

Terrestrial Bird Migration and West Nile Virus Circulation, United States

Host migration and emerging pathogens are strongly associated, especially with regard to zoonotic diseases. West Nile virus (WNV), a mosquitoborne pathogen capable of causing severe, sometimes fatal, neuroinvasive disease in humans, is maintained in highly mobile avian hosts. Using phylogeographic approaches, we investigated the relationship between WNV circulation in the United States and the flight paths of terrestrial birds. We demonstrated southward migration of WNV in the eastern flyway and northward migration in the central flyway, which is consistent with the looped flight paths of many terrestrial birds. We also identified 3 optimal locations for targeted WNV surveillance campaigns in the United States—Illinois, New York, and Texas. These results illustrate the value of multidisciplinary approaches to surveillance of infectious diseases, especially zoonotic diseases.

West Nile Virus Prevalence in Wild Birds from Mexico

West Nile virus (WNV) emerged in the Americas with its introduction in 1999 and now is considered endemic across the continent. In 2002, WNV was detected in Mexico, where its occurrence and mortality are considerably lower compared with the US. However, continuous national surveillance programs in Mexico are nonexistent. Birds are considered the primary hosts and primary geographic dispersers of this pathogen. A total of 200 cloacal and tracheal samples from wild migratory or resident birds were retrospectively analyzed using reverse transcription PCR to detect WNV from birds collected in Mexico from 2008 to 2009. The overall prevalence was 8% (16/200), and positive samples were from Oaxaca, Chiapas, and Tamaulipas in Ruby-throated Hummingbird ( Archilochus colubris), Double-crested Cormorant ( Phalacrocorax auritus), Ring-billed Gull ( Larus delawarensis), and Mourning Dove ( Zenaida macroura). Analysis of the partial sequence of the envelope gene from one of the samples from Oaxaca provided evidence that the virus belonged to the WN99 genotype. Taken together, these results demonstrated that WNV circulated in wild birds from northern and southern Mexico during the 2008-09 season, providing further information about the presence of WNV in Mexico.

Differing epidemiological dynamics of Chikungunya virus in the Americas during the 2014-2015 epidemic

Chikungunya virus (CHIKV) has been detected sporadically since the 1950s and includes three distinct co-circulating genotypes. In late 2013, the Asian genotype of CHIKV was responsible for the Caribbean outbreak (CO) that rapidly became an epidemic throughout the Americas. There is a limited understanding of the molecular evolution of CHIKV in the Americas during this epidemic. We sequenced 185 complete CHIKV genomes collected mainly from Nicaragua in Central America and Florida in the United States during the 2014-2015 Caribbean/Americas epidemic. Our comprehensive phylogenetic analyses estimated the epidemic history of the Asian genotype and the recent Caribbean outbreak (CO) clade, revealed considerable genetic diversity within the CO clade, and described different epidemiological dynamics of CHIKV in the Americas. Specifically, we identified multiple introductions in both Nicaragua and Florida, with rapid local spread of viruses in Nicaragua but limited autochthonous transmission in Florida in the US. Our phylogenetic analysis also showed phylogeographic clustering of the CO clade. In addition, we identified the significant amino acid substitutions that were observed across the entire Asian genotype during its evolution and examined amino acid changes that were specific to the CO clade. Deep sequencing analysis identified specific minor variants present in clinical specimens below-consensus levels. Finally, we investigated the association between viral phylogeny and geographic/clinical metadata in Nicaragua. To date, this study represents the largest single collection of CHIKV complete genomes during the Caribbean/Americas epidemic and significantly expands our understanding of the emergence and evolution of CHIKV CO clade in the Americas.

Large-Scale Complete-Genome Sequencing and Phylodynamic Analysis of Eastern Equine Encephalitis Virus Reveals Source-Sink Transmission Dynamics in the United States

Eastern equine encephalitis virus (EEEV) has a high case-fatality rate in horses and humans, and Florida has been hypothesized to be the source of EEEV epidemics for the northeastern United States. To test this hypothesis, we sequenced complete genomes of 433 EEEV strains collected within the United States from 1934 to 2014. Phylogenetic analysis suggested EEEV evolves relatively slowly and that transmission is enzootic in Florida, characterized by higher genetic diversity and long-term local persistence. In contrast, EEEV strains in New York and Massachusetts were characterized by lower genetic diversity, multiple introductions, and shorter local persistence. Our phylogeographic analysis supported a source-sink model in which Florida is the major source of EEEV compared to the other localities sampled. In sum, this study revealed the complex epidemiological dynamics of EEEV in different geographic regions in the United States and provided general insights into the evolution and transmission of other avian mosquito-borne viruses in this region.IMPORTANCE Eastern equine encephalitis virus (EEEV) infections are severe in horses and humans on the east coast of the United States with a >90% mortality rate in horses, an ∼33% mortality rate in humans, and significant brain damage in most human survivors. However, little is known about the evolutionary characteristics of EEEV due to the lack of genome sequences. By generating large collection of publicly available complete genome sequences, this study comprehensively determined the evolution of the virus, described the epidemiological dynamics of EEEV in different states in the United States, and identified Florida as one of the major sources. These results may have important implications for the control and prevention of other mosquito-borne viruses in the Americas.

Whole-genome sequencing and analyses identify high genetic heterogeneity, diversity and endemicity of rotavirus genotype P[6] strains circulating in Africa

Rotavirus A (RVA) exhibits a wide genotype diversity globally. Little is known about the genetic composition of genotype P[6] from Africa. This study investigated possible evolutionary mechanisms leading to genetic diversity of genotype P[6] VP4 sequences. Phylogenetic analyses on 167 P[6] VP4 full-length sequences were conducted, which included six porcine-origin sequences. Of the 167 sequences, 57 were newly acquired through whole genome sequencing as part of this study. The other 110 sequences were all publicly-available global P[6] VP4 full-length sequences downloaded from GenBank. The strength of association between the phenotypic features and the phylogeny was also determined. A number of reassortment and mixed infections of RVA genotype P[6] strains were observed in this study. Phylogenetic analyses demostrated the extensive genetic diversity that exists among human P[6] strains, porcine-like strains, their concomitant clades/subclades and estimated that P[6] VP4 gene has a higher substitution rate with the mean of 1.05E-3 substitutions/site/year. Further, the phylogenetic analyses indicated that genotype P[6] strains were endemic in Africa, characterised by an extensive genetic diversity and long-time local evolution of the viruses. This was also supported by phylogeographic clustering and G-genotype clustering of the P[6] strains when Bayesian Tip-association Significance testing (BaTS) was applied, clearly supporting that the viruses evolved locally in Africa instead of spatial mixing among different regions. Overall, the results demonstrated that multiple mechanisms such as reassortment events, various mutations and possibly interspecies transmission account for the enormous diversity of genotype P[6] strains in Africa. These findings highlight the need for continued global surveillance of rotavirus diversity.

Full-Genome Characterization and Genetic Evolution of West African Isolates of Bagaza Virus

Bagaza virus is a mosquito-borne flavivirus, first isolated in 1966 in Central African Republic. It has currently been identified in mosquito pools collected in the field in West and Central Africa. Emergence in wild birds in Europe and serological evidence in encephalitis patients in India raise questions on its genetic evolution and the diversity of isolates circulating in Africa. To better understand genetic diversity and evolution of Bagaza virus, we describe the full-genome characterization of 11 West African isolates, sampled from 1988 to 2014. Parameters such as genetic distances, N-glycosylation patterns, recombination events, selective pressures, and its codon adaptation to human genes are assessed. Our study is noteworthy for the observation of N-glycosylation and recombination in Bagaza virus and provides insight into its Indian origin from the 13th century. Interestingly, evidence of Bagaza virus codon adaptation to human house-keeping genes is also observed to be higher than those of other flaviviruses well known in human infections. Genetic variations on genome of West African Bagaza virus could play an important role in generating diversity and may promote Bagaza virus adaptation to other vertebrates and become an important threat in human health.

The Effects of Sampling Location and Predictor Point Estimate Certainty on Posterior Support in Bayesian Phylogeographic Generalized Linear Models

The use of generalized linear models in Bayesian phylogeography has enabled researchers to simultaneously reconstruct the spatiotemporal history of a virus and quantify the contribution of predictor variables to that process. However, little is known about the sensitivity of this method to the choice of the discrete state partition. Here we investigate this question by analyzing a data set containing 299 sequences of the West Nile virus envelope gene sampled in the United States and fifteen predictors aggregated at four spatial levels. We demonstrate that although the topology of the viral phylogenies was consistent across analyses, support for the predictors depended on the level of aggregation. In particular, we found that the variance of the predictor support metrics was minimized at the most precise level for several predictors and maximized at more sparse levels of aggregation. These results suggest that caution should be taken when partitioning a region into discrete locations to ensure that interpretable, reproducible posterior estimates are obtained. These results also demonstrate why researchers should use the most precise discrete states possible to minimize the posterior variance in such estimates and reveal what truly drives the diffusion of viruses.

The presence of duck Tembusu virus in Thailand since 2007: A retrospective study

Duck Tembusu virus (DTMUV), a newly emerging virus in ducks, was first reported in China in 2010. However, an unknown severe contagious disease associated with severe neurological signs and egg production losses in ducks, resembling to DTMUV infection, was observed in Thailand since 2007. To determine the presence of DTMUV in 2007, the clinical samples from affected ducks collected in 2007 were tested for DTMUV using pathological and virological analyses. Gross and histopathological lesions of affected ducks were mostly restricted to the ovary, brain and spinal cord, and correlated with the presence of flavivirus antigen in the brain and spinal cord samples. Subsequently, DTMUV was identified by RT-PCR and nucleotide sequencing of the polyprotein gene. Phylogenetic analysis of the polyprotein gene sequence revealed that the 2007 Thai DTMUV was a unique virus, belonged within DTMUV cluster 1, but distinctively separated from the Malaysian DTMUV, which was the most closely related DTMUV. It is interesting to note that the 2007 Thai DTMUV was genetically different from the currently circulating Thai and Chinese DTMUVs, which belonged to cluster 2. Our findings indicated that the 2007 Thai DTMUV emerged earlier from a common ancestor with the recently reported DTMUVs; however, it was genetically distinctive to any of the currently circulating DTMUVs. In conclusion, our data demonstrated the presence of DTMUV in the Thai ducks since 2007, prior to the first report of DTMUV in China in 2010. This study indicates that DTMUV may have circulated in the region long before 2010 and highlights high genetic diversity of DTMUVs in Asia.

Are RNA Viruses Candidate Agents for the Next Global Pandemic? A Review

Pathogenic RNA viruses are potentially the most important group involved in zoonotic disease transmission, and they represent a challenge for global disease control. Their biological diversity and rapid adaptive rates have proved to be difficult to overcome and to anticipate by modern medical technology. Also, the anthropogenic change of natural ecosystems and the continuous population growth are driving increased rates of interspecies contacts and the interchange of pathogens that can develop into global pandemics. The combination of molecular, epidemiological, and ecological knowledge of RNA viruses is therefore essential towards the proper control of these emergent pathogens. This review outlines, throughout different levels of complexity, the problems posed by RNA viral diseases, covering some of the molecular mechanisms allowing them to adapt to new host species-and to novel pharmaceutical developments-up to the known ecological processes involved in zoonotic transmission.

Molecular Evolution of the VP1 Gene in Human Norovirus GII.4 Variants in 1974-2015

Human norovirus (HuNoV) is a leading cause of viral gastroenteritis worldwide, of which GII.4 is the most predominant genotype. Unlike other genotypes, GII.4 has created various variants that escaped from previously acquired immunity of the host and caused repeated epidemics. However, the molecular evolutionary differences among all GII.4 variants, including recently discovered strains, have not been elucidated. Thus, we conducted a series of bioinformatic analyses using numerous, globally collected, full-length GII.4 major capsid (VP1) gene sequences (466 strains) to compare the evolutionary patterns among GII.4 variants. The time-scaled phylogenetic tree constructed using the Bayesian Markov chain Monte Carlo (MCMC) method showed that the common ancestor of the GII.4 VP1 gene diverged from GII.20 in 1840. The GII.4 genotype emerged in 1932, and then formed seven clusters including 14 known variants after 1980. The evolutionary rate of GII.4 strains was estimated to be 7.68 × 10⁻³ substitutions/site/year. The evolutionary rates probably differed among variants as well as domains [protruding 1 (P1), shell, and P2 domains]. The Osaka 2007 variant strains probably contained more nucleotide substitutions than any other variant. Few conformational epitopes were located in the shell and P1 domains, although most were contained in the P2 domain, which, as previously established, is associated with attachment to host factors and antigenicity. We found that positive selection sites for the whole GII.4 genotype existed in the shell and P1 domains, while Den Haag 2006b, New Orleans 2009, and Sydney 2012 variants were under positive selection in the P2 domain. Amino acid substitutions overlapped with putative epitopes or were located around the epitopes in the P2 domain. The effective population sizes of the present strains increased stepwise for Den Haag 2006b, New Orleans 2009, and Sydney 2012 variants. These results suggest that HuNoV GII.4 rapidly evolved in a few decades, created various variants, and altered its evolutionary rate and antigenicity.

Evidence of increasing diversification of Zika virus strains isolated in the American continent

Zika virus (ZIKV) is a member of the family Flaviviridae. ZIKV emerged in Brazil in 2015, causing an unprecedented epidemic and since then the virus has rapidly spread throughout the Americas. These facts highlight the need of detailed phylogenetic studies to understand the emergence, spread, and evolution of ZIKV populations. For these reasons, a Bayesian coalescent Markov Chain Monte Carlo analysis of complete genome sequences of ZIKV strains recently isolated in the American continent was performed. The results of these studies revealed an increasing diversification of ZIKV strains in different genetic lineages and co-circulation of distinct genetic lineages in several countries in the region. The time of the most recent common ancestor (tMRCA) was established to be around February 20, 2014 for ZIKV strains circulating in the American region. A mean rate of evolution of 1.55 × 10^-3 substitutions/site/year was obtained for ZIKV strains included in this study. A Bayesian skyline plot indicate a sharp increase in population size from February 2014 to July 2015 and a decline during 2016. These results are discussed in terms of the emergence and evolution of ZIKV populations in the American continent.

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.

Correlation of West Nile Virus Incidence in Donated Blood with West Nile Neuroinvasive Disease Rates, United States, 2010-2012

Over the past decade, West Nile virus (WNV) has spread across the United States. We aggregated blood donor data from 2010–2012 and then calculated the incidence of WNV RNA–positive donations and compared the incidence with neuroinvasive disease (NID) case data from the ArboNET surveillance system. Of 10,107,853 donations, 640 were confirmed positive. The seasonal WNV incidence rate per 100,000 persons was 33.4 (95% CI 22–45) in 2010, 25.7 (95% CI 15–34) in 2011, and 119.9 (95% CI 98–141) in 2012. NID to blood donor ratios were 1 in 164 (95% CI 152–178) in 2010, 1 in 158 (95% CI 145–174) in 2011, and 1 in 131 (95% CI 127–136) in 2012. We updated estimates of the ratio of NID to WNV infection rates, demonstrating stable disease penetrance over the study period. Blood donor WNV RNA screening is a valuable public health tool for WNV surveillance.

Molecular evolution of emerging Banna virus

Banna virus (BAV) is an emerging pathogen that causes human viral encephalitis and has been isolated from types of blood-sucking insects and mammals in Asia. However, there are no reported systematic studies that describe the origin and evolution of BAV. Here, a phylogenetic analysis of BAVs isolated from a variety of potential vectors and vertebrate hosts worldwide revealed that BAVs emerged in the beginning of the 20th century and do not exhibit a species barrier. The mean substitution rate of BAVs was 2.467×10^-2substitution/site/year (95% HPD, 1.093×10^-3 to 5.628×10^-2). The lineage is mainly composed of BAVs from high-latitude regions, which are the most recently emerged viruses with significantly higher substitution rates compared with the lineage comprised of the isolates from middle or low-latitude regions. The genetic differences between BAV strains are positively correlated with the geographic distribution. Strains from the same latitude regions are almost 100% identical, whereas the differences between strains from long distance regions with different latitudes could be >60%. Our results demonstrate that BAV is an emerging virus at a stage that involves rapid evolution and has great potential for introduction into non-endemic areas. Thus, enhanced surveillance of BAV is highly recommended worldwide.

Genetic Variability of West Nile Virus in U.S. Blood Donors from the 2012 Epidemic Season

West Nile virus (WNV) is an arbovirus maintained in nature in a bird-mosquito enzootic cycle which can also infect other vertebrates including humans. WNV is now endemic in the United States (U.S.), causing yearly outbreaks that have resulted in an estimated total of 4-5 million human infections. Over 41,700 cases of West Nile disease, including 18,810 neuroinvasive cases and 1,765 deaths, were reported to the CDC between 1999 and 2014. In 2012, the second largest West Nile outbreak in the U.S. was reported, which caused 5,674 cases and 286 deaths. WNV continues to evolve, and three major WNV lineage I genotypes (NY99, WN02, and SW/WN03) have been described in the U.S. since introduction of the virus in 1999. We report here the WNV sequences obtained from 19 human samples acquired during the 2012 U.S. outbreak and our examination of the evolutionary dynamics in WNV isolates sequenced from 1999-2012. Maximum-likelihood and Bayesian methods were used to perform the phylogenetic analyses. Selection pressure analyses were performed with the HyPhy package using the Datamonkey web-server. Using different codon-based and branch-site selection models, we detected a number of codons subjected to positive pressure in WNV genes. Thirteen of the 19 completely sequenced isolates from 10 U.S. states were genetically similar, sharing up to 55 nucleotide mutations and 4 amino acid substitutions when compared with the prototype isolate WN-NY99. Overall, these analyses showed that following a brief contraction in 2008-2009, WNV genetic divergence in the U.S. continued to increase in 2012, and that closely related variants were found across a broad geographic range of the U.S., coincident with the second-largest WNV outbreak in U.S.

Reconstruction of the Evolutionary History and Dispersal of Usutu Virus, a Neglected Emerging Arbovirus in Europe and Africa

Usutu virus (USUV), one of the most neglected Old World encephalitic flaviviruses, causes epizootics among wild and captive birds and sporadic infection in humans. The dynamics of USUV spread and evolution in its natural hosts are unknown. Here, we present the phylogeny and evolutionary history of all available USUV strains, including 77 newly sequenced complete genomes from a variety of host species at a temporal and spatial scaled resolution. The results showed that USUV can be classified into six distinct lineages and that the most recent common ancestor of the recent European epizootics emerged in Africa at least 500 years ago. We demonstrated that USUV was introduced regularly from Africa into Europe in the last 50 years, and the genetic diversity of European lineages is shaped primarily by in situ evolution, while the African lineages have been driven by extensive gene flow. Most of the amino acid changes are deleterious polymorphisms removed by purifying selection, with adaptive evolution restricted to the NS5 gene and several others evolving under episodic directional selection, indicating that the ecological or immunological factors were mostly the key determinants of USUV dispersal and outbreaks. Host-specific mutations have been detected, while the host transition analysis identified mosquitoes as the most likely origin of the common ancestor and birds as the source of the recent European USUV lineages. Our results suggest that the major migratory bird flyways could predict the continental and intercontinental dispersal patterns of USUV and that migratory birds might act as potential long-distance dispersal vehicles.

Usutu virus (USUV), a mosquito-borne flavivirus of the Japanese encephalitis virus antigenic group, caused massive bird die-offs, mostly in Europe. There is increasing evidence that USUV appears to be pathogenic for humans, becoming a potential public health problem. The emergence of USUV in Europe allows us to understand how an arbovirus spreads, adapts, and evolves in a naive environment. Thus, understanding the epidemiological and evolutionary processes that contribute to the emergence, maintenance, and further spread of viral diseases is the sine qua non to develop and implement surveillance strategies for their control. In this work, we performed an expansive phylogeographic and evolutionary analysis of USUV using all published sequences and those generated during this study. Subsequently, we described the genetic traits, reconstructed the potential pattern of geographic spread between continents/countries of the identified viral lineages and the drivers of viral migration, and traced the origin of outbreaks and transition events between different hosts.

Incidence of West Nile virus infection in the Dallas–Fort Worth metropolitan area during the 2012 epidemic

The 2012 West Nile virus (WNV) epidemic was the largest since 2003 and the North Texas region was the most heavily impacted. We conducted a serosurvey of blood donors from four counties in the Dallas–Fort Worth area to characterize the epidemic. Blood donor specimens collected in November 2012 were tested for WNV-specific antibodies. Donors positive for WNV-specific IgG, IgM, and neutralizing antibodies were considered to have been infected in 2012. This number was adjusted using a multi-step process that accounted for timing of IgM seroreversion determined from previous longitudinal studies of WNV-infected donors. Of 4971 donations screened, 139 (2·8%) were confirmed WNV IgG positive, and 69 (1·4%) had IgM indicating infection in 2012. After adjusting for timing of sampling and potential seroreversion, we estimated that 1·8% [95% confidence interval (CI) 1·5–2·2] of the adult population in the Dallas–Fort Worth area were infected during 2012. The resulting overall estimate for the ratio of infections to reported WNV neuroinvasive disease (WNND) cases was 238:1 (95% CI 192–290), with significantly increased risk of WNND in older age groups. These findings were very similar to previous estimates of infections per WNND case, indicating no change in virulence as WNV evolved into an endemic infection in the United States.

Purpura fulminans associated with acute West Nile virus encephalitis

Purpura fulminans is a progressive thrombotic disorder that presents with widespread purpura due to deficiency or dysfunction of protein C or protein S. Lesions present as well-demarcated erythematous macules that progress to irregular areas of hemorrhagic necrosis.West Nile virus is a member of the Flaviviridae family transmitted to humans through the bite of various mosquito species. It manifests as West Nile fever in 25% of those infected and less commonly as neuroinvasive disease. An African American man in his fortiespresented with altered mental status and was noted to have evidence of disseminated intravascular coagulation according to his lab data. He then developed dusky skin discoloration and systemic flaccid bullae with desquamation. Biopsy was consistent with purpura fulminans and the patient eventually developed symmetric peripheral gangrene, requiring amputations of all four extremities. Infectious work up revealed positive testing for IgM and IgG antibodies in serum and cerebrospinal fluid leading to the diagnosis of acute West Nile Virus encephalitis. We present this case to describe the rarely reported association of purpura fulminans with West Nile Virus infection.

Mosquito Surveillance for 15 Years Reveals High Genetic Diversity Among West Nile Viruses in Israel

West Nile Virus (WNV) is endemic in Israel and has been the cause of several outbreaks in recent years. In 2000, a countrywide mosquito survey was established to monitor WNV activity and characterize viral genotypes in Israel. We analyzed data from 7135 pools containing 277 186 mosquitoes collected over the past 15 years and, here, report partial sequences of WNV genomes obtained from 102 of the 336 positive mosquito pools. Phylogenetic analysis demonstrated that cluster 4 and the Mediterranean and Eastern European subtypes of cluster 2 within WNV lineage 1 circulated in Israel, as did WNV lineage 2, highlighting a high genetic diversity of WNV genotypes in our region. As a major crossroads for bird migration between Africa and Eurasia and with a long history of human infection, Israel serves as a resource hub for WNV in Africa and Eurasia and provides valuable information on WNV circulation in these regions.

Fluid Spatial Dynamics of West Nile Virus in the United States: Rapid Spread in a Permissive Host Environment

The introduction of West Nile virus (WNV) into North America in 1999 is a classic example of viral emergence in a new environment, with its subsequent dispersion across the continent having a major impact on local bird populations. Despite the importance of this epizootic, the pattern, dynamics, and determinants of WNV spread in its natural hosts remain uncertain. In particular, it is unclear whether the virus encountered major barriers to transmission, or spread in an unconstrained manner, and if specific viral lineages were favored over others indicative of intrinsic differences in fitness. To address these key questions in WNV evolution and ecology, we sequenced the complete genomes of approximately 300 avian isolates sampled across the United States between 2001 and 2012. Phylogenetic analysis revealed a relatively star-like tree structure, indicative of explosive viral spread in the United States, although with some replacement of viral genotypes through time. These data are striking in that viral sequences exhibit relatively limited clustering according to geographic region, particularly for those viruses sampled from birds, and no strong phylogenetic association with well-sampled avian species. The genome sequence data analyzed here also contain relatively little evidence for adaptive evolution, particularly of structural proteins, suggesting that most viral lineages are of similar fitness and that WNV is well adapted to the ecology of mosquito vectors and diverse avian hosts in the United States. In sum, the molecular evolution of WNV in North America depicts a largely unfettered expansion within a permissive host and geographic population with little evidence of major adaptive barriers.

IMPORTANCE

How viruses spread in new host and geographic environments is central to understanding the emergence and evolution of novel infectious diseases and for predicting their likely impact. The emergence of the vector-borne West Nile virus (WNV) in North America in 1999 represents a classic example of this process. Using approximately 300 new viral genomes sampled from wild birds, we show that WNV experienced an explosive spread with little geographical or host constraints within birds and relatively low levels of adaptive evolution. From its introduction into the state of New York, WNV spread across the United States, reaching California and Florida within 4 years, a migration that is clearly reflected in our genomic sequence data, and with a general absence of distinct geographical clusters of bird viruses. However, some geographically distinct viral lineages were found to circulate in mosquitoes, likely reflecting their limited long-distance movement compared to avian species.

Phylogenetic characterization of Central/Southern European lineage 2 West Nile virus: analysis of human outbreaks in Italy and Greece, 2013-2014

In recent years, West Nile virus (WNV) lineage 2 has been spreading and causing disease outbreaks in humans and animals in Europe. In order to characterize viral diversity, we performed full-length genome sequencing of WNV lineage 2 from human samples collected during outbreaks in Italy and Greece in 2013 and 2014. Phylogenetic analysis showed that these WNV lineage 2 genomes belonged to a monophyletic clade derived from a single introduction into Europe of the prototype Hungarian strain. Correlation of phylogenetic data with geospatial information showed geographical clustering of WNV genome sequences both in Italy and in Greece, indicating that the virus had evolved and diverged during its dispersal in Europe, leading to the emergence of novel genotypes, as it adapted to local ecological niches. These genotypes carried divergent conserved amino acid substitutions, which might have been relevant for viral adaptation, as suggested by selection pressure analysis and in silico and experimental modelling of sequence changes. In conclusion, the results of this study provide further information on WNV lineage 2 transmission dynamics in Europe, and emphasize the need for WNV surveillance activities to monitor viral evolution and diversity.

Evolutionary History of Cathepsin L (L-like) Family Genes in Vertebrates

Cathepsin L family, an important cysteine protease found in lysosomes, is categorized into cathepsins B, F, H, K, L, S, and W in vertebrates. This categorization is based on their sequence alignment and traditional functional classification, but the evolutionary relationship of family members is unclear. This study determined the evolutionary relationship of cathepsin L family genes in vertebrates through phylogenetic construction. Results showed that cathepsins F, H, S and K, and L and V were chronologically diverged. Tandem-repeat duplication was found to occur in the evolutionary history of cathepsin L family. Cathepsin L in zebrafish, cathepsins S and K in xenopus, and cathepsin L in mice and rats underwent evident tandem-repeat events. Positive selection was detected in cathepsin L-like members in mice and rats, and amino acid sites under positive selection pressure were calculated. Most of these sites appeared at the connection of secondary structures, suggesting that the sites may slightly change spatial structure. Severe positive selection was also observed in cathepsin V (L2) of primates, indicating that this enzyme had some special functions. Our work provided a brief evolutionary history of cathepsin L family and differentiated cathepsins S and K from cathepsin L based on vertebrate appearance. Positive selection was the specific cause of differentiation of cathepsin L family genes, confirming that gene function variation after expansion events was related to interactions with the environment and adaptability.

Evolutionary analysis of Tembusu virus: evidence for the emergence of a dominant genotype

We recently identified Tembusu virus (TMUV) as a causative agent of duck infectious disease, which has spread in China since 2010. A recent study has indicated a potential case of human infection by TMUV, highlighting the need for further study of TMUV, especially its evolution. Because no information exists regarding the evolution of TMUV, we conducted comprehensive phylogenetic analyses using the largest collection of complete open reading frame (ORF) sequences of TMUV. Our results indicated that two lineages of viruses were associated with the 2010 outbreak in China, and lineage II, in particular sublineage II-c, has arisen as the dominant lineage currently circulating. We inferred that the most recent common ancestor (MRCA) of this TMUV was emerged around 1996. Evidence of natural recombination was also detected in TMUV. Molecular adaptation analyses revealed that strong negative selection shaped the evolution of TMUV, while a number of codons subjected to positive pressure were also identified. Our study, for the first time, illustrated the evolutionary history and character of TMUV and will be helpful for vaccine and diagnostic development.

West nile virus and equine encephalitis viruses: new perspectives

Mosquito-borne diseases affect horses worldwide. Mosquito-borne diseases generally cause encephalomyelitis in the horse and can be difficult to diagnose antemortem. In addition to general disease, and diagnostic and treatment aspects, this review article summarizes the latest information on these diseases, covering approximately the past 5 years, with a focus on new equine disease encroachments, diagnostic and vaccination aspects, and possible therapeutics on the horizon.

Complete Genome Sequence of West Nile Virus Strains Used for the Formulation of CBER/FDA RNA Reference Reagents and Lot Release Panels for Nucleic Acid Testing

We report the complete sequences of two West Nile virus strains (FDA-Hu02 and NY99) used for the formulation of CBER/FDA RNA reference reagents and lot release panels for use with nucleic acid technology testing.

Co-circulation of West Nile virus variants, Arizona, USA, 2010

Molecular analysis of West Nile virus (WNV) isolates obtained during a 2010 outbreak in Maricopa County, Arizona, USA, demonstrated co-circulation of 3 distinct genetic variants, including strains with novel envelope protein mutations. These results highlight the continuing evolution of WNV in North America and the current complexity of WNV dispersal and transmission.

Characterization of rubella virus genotypes among pregnant women in northern Vietnam, 2011-2013

Rubella virus (RV) infection is an unresolved clinical complication that affects children in developing countries including Vietnam. RV infection during the first trimester of pregnancy causes severe birth defects known as congenital rubella syndrome. This study reports on the genomic characterization of RV strains circulating in northern Vietnam during 2011-2013. RV-IgM positive amniotic fluid specimens were collected from 38 women from northern Vietnam who presented with clinical rubella at the National Hospital of Obstetrics and Gynecology in Hanoi, Vietnam. The RV genes were determined by nested PCR with primers amplifying the 739-nucleotide coding region of the E1 gene. The sequences from the amplified DNA fragments were phylogenetically analyzed and compared to reference RV strains. Seventeen out of 38 samples are positive for RV detecting. All new RV isolates are clustered to genotype 2B. Eighteen amino acid mutations were found in the T and B cell epitopes. These results suggest that genotype 2B RV strains frequently circulate in northern Vietnam. These data describe the RV genotype in Vietnam with the aim of improving maternal and child health in this country.

Phylogeography and molecular epidemiology of an epidemic strain of dengue virus type 1 in Sri Lanka

In 2009, a severe epidemic of dengue disease occurred in Sri Lanka, with higher mortality and morbidity than any previously recorded epidemic in the country. It corresponded to a shift to dengue virus 1 as the major disease-causing serotype in Sri Lanka. Dengue disease reached epidemic levels in the next 3 years. We report phylogenetic evidence that the 2009 epidemic DENV-1 strain continued to circulate within the population and caused severe disease in the epidemic of 2012. Bayesian phylogeographic analyses suggest that the 2009 Sri Lankan epidemic DENV-1 strain may have traveled directly or indirectly from Thailand through China to Sri Lanka, and after spreading within the Sri Lankan population, it traveled to Pakistan and Singapore. Our findings delineate the dissemination route of a virulent DENV-1 strain in Asia. Understanding such routes will be of particular importance to global control efforts.

Evolutionary dynamics of West Nile virus in Georgia, 2001-2011

From 1999-2001, West Nile virus (WNV) spread throughout the eastern United States (US) and was first detected in Georgia in 2001. To date, the virus has been detected in over 2,500 dead wild bird and mosquito samples from across Georgia. We sequenced the premembrane (preM) and envelope gene (E) (2004 bp) from 111 isolates collected from 2001 to 2011. To assess viral gene flow from other geographic regions in the US, we combined our data with WNV sequences available at the National Center for Biotechnology Information (NCBI) and performed phylogenetic analysis. We found evidence that WNV isolates detected in Chatham County Georgia most likely originated from the Northeastern United States. These results highlight the growing importance of adequate genetic surveillance for monitoring and controlling viruses of public health concern.