Patterns of intra- and interhost nonsynonymous variation reveal strong purifying selection in dengue virus

Binding Evolution of the Dengue Virus Envelope Against DC-SIGN: A Combined Approach of Phylogenetics and Molecular Dynamics Analyses Over 30 Years of Dengue Virus in Brazil

The Red Queen Hypothesis (RQH), derived from Lewis Carroll's "Through the Looking-Glass", postulates that organisms must continually adapt in response to each other to maintain relative fitness. Within the context of host-pathogen interactions, the RQH implies an evolutionary arms race, wherein viruses evolve to exploit hosts and hosts evolve to resist viral invasion. This study delves into the dynamics of the RQH in the context of virus-cell interactions, specifically focusing on virus receptors and cell receptors. We observed multiple virus-host systems and noted patterns of co-evolution. As viruses evolved receptor-binding proteins to effectively engage with cell receptors, cells countered by altering their receptor genes. This ongoing mutual adaptation cycle has influenced the molecular intricacies of receptor-ligand interactions. Our data supports the RQH as a driving force behind the diversification and specialization of both viral and host cell receptors. Understanding this co-evolutionary dance offers insights into the unpredictability of emerging viral diseases and potential therapeutic interventions. Future research is crucial to dissect the nuanced molecular changes and the broader ecological consequences of this ever-evolving battle. Here, we combine phylogenetic inferences, structural modeling, and molecular dynamics analyses to describe the epidemiological characteristics of major Brazilian DENV strains that circulated from 1990 to 2022 from a combined perspective, thus providing us with a more detailed picture on the dynamics of such interactions over time.

Molecular epidemiology of dengue in Malaysia: 2015-2021

Dengue has been one of the major public health problems in Malaysia for decades. Over 600,000 dengue cases and 1,200 associated fatalities have been reported in Malaysia from 2015 to 2021, which was 100% increase from the cumulative total of dengue cases reported during the preceding 07-year period from 2008 to 2014. However, studies that describe the molecular epidemiology of dengue in Malaysia in recent years are limited. In the present study, we describe the genetic composition and dispersal patterns of Dengue virus (DENV) by using 4,004 complete envelope gene sequences of all four serotypes (DENV-1 = 1,567, DENV-2 = 1,417, DENV-3 = 762 and DENV-4 = 258) collected across Malaysia from 2015 to 2021. The findings revealed that DENV populations in Malaysia were highly diverse, and the overall heterogeneity was maintained through repetitive turnover of genotypes. Phylogeography analyses suggested that DENV dispersal occurred through an extensive network, mainly among countries in South and East Asia and Malaysian states, as well as among different states, especially within Peninsular Malaysia. The results further suggested Selangor and Johor as major hubs of DENV emergence and spread in Malaysia.

Evidence of Zika Virus Reinfection by Genome Diversity and Antibody Response Analysis, Brazil

We generated 238 Zika virus (ZIKV) genomes from 135 persons in Brazil who had samples collected over 1 year to evaluate virus persistence. Phylogenetic inference clustered the genomes together with previously reported ZIKV strains from northern Brazil, showing that ZIKV has been remained relatively stable over time. Temporal phylogenetic analysis revealed limited within-host diversity among most ZIKV-persistent infected associated samples. However, we detected unusual virus temporal diversity from >5 persons, uncovering the existence of divergent genomes within the same patient. All those patients showed an increase in neutralizing antibody levels, followed by a decline at the convalescent phase of ZIKV infection. Of interest, in 3 of those patients, titers of neutralizing antibodies increased again after 6 months of ZIKV infection, concomitantly with real-time reverse transcription PCR re-positivity, supporting ZIKV reinfection events. Altogether, our findings provide evidence for the existence of ZIKV reinfection events.

Genetic Diversity and Dispersal of DENGUE Virus among Three Main Island Groups of the Philippines during 2015-2017

Dengue has been one of the major public health concerns in the Philippines for more than a century. The annual dengue case burden has been increasing in recent years, exceeding 200,000 in 2015 and 2019. However, there is limited information on the molecular epidemiology of dengue in the Philippines. We, therefore, conducted a study to understand the genetic composition and dispersal of DENV in the Philippines from 2015 to 2017 under UNITEDengue. Our analyses included 377 envelope (E) gene sequences of all 4 serotypes obtained from infections in 3 main island groups (Luzon, Visayas, and Mindanao) of the Philippines. The findings showed that the overall diversity of DENV was generally low. DENV-1 was relatively more diverse than the other serotypes. Virus dispersal was evident among the three main island groups, but each island group demonstrated a distinct genotype composition. These observations suggested that the intensity of virus dispersal was not substantive enough to maintain a uniform heterogeneity among island groups so that each island group behaved as an independent epidemiological unit. The analyses suggested Luzon as one of the major sources of DENV emergence and CAR, Calabarzon, and CARAGA as important hubs of virus dispersal in the Philippines. Our findings highlight the importance of virus surveillance and molecular epidemiological analyses to gain deep insights into virus diversity, lineage dominance, and dispersal patterns that could assist in understanding the epidemiology and transmission risk of dengue in endemic regions.

Intrahost Genetic Diversity of Dengue Virus in Human Hosts and Mosquito Vectors under Natural Conditions Which Impact Replicative Fitness In Vitro

Dengue virus (DENV) is an arbovirus whose transmission cycle involves disparate hosts: humans and mosquitoes. The error-prone nature of viral RNA replication drives the high mutation rates, and the consequently high genetic diversity affects viral fitness over this transmission cycle. A few studies have been performed to investigate the intrahost genetic diversity between hosts, although their mosquito infections were performed artificially in the laboratory setting. Here, we performed whole-genome deep sequencing of DENV-1 (n = 11) and DENV-4 (n = 13) derived from clinical samples and field-caught mosquitoes from the houses of naturally infected patients, in order to analyze the intrahost genetic diversity of DENV between host types. Prominent differences in DENV intrahost diversity were observed in the viral population structure between DENV-1 and DENV-4, which appear to be associated with differing selection pressures. Interestingly, three single amino acid substitutions in the NS2A (K81R), NS3 (K107R), and NS5 (I563V) proteins in DENV-4 appear to be specifically acquired during infection in Ae. aegypti mosquitoes. Our in vitro study shows that the NS2A (K81R) mutant replicates similarly to the wild-type infectious clone-derived virus, while the NS3 (K107R), and NS5 (I563V) mutants have prolonged replication kinetics in the early phase in both Vero and C6/36 cells. These findings suggest that DENV is subjected to selection pressure in both mosquito and human hosts. The NS3 and NS5 genes may be specific targets of diversifying selection that play essential roles in early processing, RNA replication, and infectious particle production, and they are potentially adaptive at the population level during host switching.

Evolutionary dynamics of dengue virus in India

More than a hundred thousand dengue cases are diagnosed in India annually, and about half of the country's population carries dengue virus-specific antibodies. Dengue propagates and adapts to the selection pressures imposed by a multitude of factors that can lead to the emergence of new variants. Yet, there has been no systematic analysis of the evolution of the dengue virus in the country. Here, we present a comprehensive analysis of all DENV gene sequences collected between 1956 and 2018 from India. We examine the spatio-temporal dynamics of India-specific genotypes, their evolutionary relationship with global and local dengue virus strains, interserotype dynamics and their divergence from the vaccine strains. Our analysis highlights the co-circulation of all DENV serotypes in India with cyclical outbreaks every 3-4 years. Since 2000, genotype III of DENV-1, cosmopolitan genotype of DENV-2, genotype III of DENV-3 and genotype I of DENV-4 have been dominating across the country. Substitution rates are comparable across the serotypes, suggesting a lack of serotype-specific evolutionary divergence. Yet, the envelope (E) protein displays strong signatures of evolution under immune selection. Apart from drifting away from its ancestors and other contemporary serotypes in general, we find evidence for recurring interserotype drift towards each other, suggesting selection via cross-reactive antibody-dependent enhancement. We identify the emergence of the highly divergent DENV-4-Id lineage in South India, which has acquired half of all E gene mutations in the antigenic sites. Moreover, the DENV-4-Id is drifting towards DENV-1 and DENV-3 clades, suggesting the role of cross-reactive antibodies in its evolution. Due to the regional restriction of the Indian genotypes and immunity-driven virus evolution in the country, ~50% of all E gene differences with the current vaccines are focused on the antigenic sites. Our study shows how the dengue virus evolution in India is being shaped in complex ways.

Developing an appropriate evolutionary baseline model for the study of SARS-CoV-2 patient samples

Over the past 3 years, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has spread through human populations in several waves, resulting in a global health crisis. In response, genomic surveillance efforts have proliferated in the hopes of tracking and anticipating the evolution of this virus, resulting in millions of patient isolates now being available in public databases. Yet, while there is a tremendous focus on identifying newly emerging adaptive viral variants, this quantification is far from trivial. Specifically, multiple co-occurring and interacting evolutionary processes are constantly in operation and must be jointly considered and modeled in order to perform accurate inference. We here outline critical individual components of such an evolutionary baseline model-mutation rates, recombination rates, the distribution of fitness effects, infection dynamics, and compartmentalization-and describe the current state of knowledge pertaining to the related parameters of each in SARS-CoV-2. We close with a series of recommendations for future clinical sampling, model construction, and statistical analysis.

Genome evolution of dengue virus serotype 1 under selection by Wolbachia pipientis in Aedes aegypti mosquitoes

The introgression of antiviral strains of Wolbachia into Aedes aegypti mosquito populations is a public health intervention for the control of dengue. Plausibly, dengue virus (DENV) could evolve to bypass the antiviral effects of Wolbachia and undermine this approach. Here, we established a serial-passage system to investigate the evolution of DENV in Ae. aegypti mosquitoes infected with the wMel strain of Wolbachia. Using this system, we report on virus genetic outcomes after twenty passages of serotype 1 of DENV (DENV-1). An amino acid substitution, E203K, in the DENV-1 envelope protein was more frequently detected in the consensus sequence of virus populations passaged in wMel-infected Ae. aegypti than wild-type counterparts. Positive selection at residue 203 was reproducible; it occurred in passaged virus populations from independent DENV-1-infected patients and also in a second, independent experimental system. In wild-type mosquitoes and human cells, the 203K variant was rapidly replaced by the progenitor sequence. These findings provide proof of concept that wMel-associated selection of virus populations can occur in experimental conditions. Field-based studies are needed to explore whether wMel imparts selective pressure on DENV evolution in locations where wMel is established.

Contribution of phylogenetics to understanding the evolution and epidemiology of dengue virus

Dengue virus (DENV) is one of the most important arboviral pathogens in the tropics and subtropics, and nearly one‐third of the world's population is at risk of infection. The transmission of DENV involves a sylvatic cycle between nonhuman primates (NHP) and Aedes genus mosquitoes, and an endemic cycle between human hosts and predominantly Aedes aegypti. DENV belongs to the genus Flavivirus of the family Flaviviridae and consists of four antigenically distinct serotypes (DENV‐1‐4). Phylogenetic analyses of DENV have revealed its origin, epidemiology, and the drivers that determine its molecular evolution in nature. This review discusses how phylogenetic research has improved our understanding of DENV evolution and how it affects viral ecology and improved our ability to analyze and predict future DENV emergence.

Human pathogenic RNA viruses establish noncompeting lineages by occupying independent niches

Numerous pathogenic viruses are endemic in humans and cause a broad variety of diseases, but what is their potential for causing new pandemics? We show that most human pathogenic RNA viruses form multiple, cocirculating lineages with low turnover rates. These lineages appear to be largely noncompeting and occupy distinct epidemiological niches that are not regionally or seasonally defined, and their persistence appears to stem from limited outbreaks in small communities so that only a small fraction of the global susceptible population is infected at any time. However, due to globalization, interaction and competition between lineages might increase, potentially leading to increased diversification and pathogenicity. Thus, endemic viruses appear to merit global attention with respect to the prevention of future pandemics.

Many pathogenic viruses are endemic among human populations and can cause a broad variety of diseases, some potentially leading to devastating pandemics. How virus populations maintain diversity and what selective pressures drive population turnover is not thoroughly understood. We conducted a large-scale phylodynamic analysis of 27 human pathogenic RNA viruses spanning diverse life history traits, in search of unifying trends that shape virus evolution. For most virus species, we identify multiple, cocirculating lineages with low turnover rates. These lineages appear to be largely noncompeting and likely occupy semiindependent epidemiological niches that are not regionally or seasonally defined. Typically, intralineage mutational signatures are similar to interlineage signatures. The principal exception are members of the family Picornaviridae, for which mutations in capsid protein genes are primarily lineage defining. Interlineage turnover is slower than expected under a neutral model, whereas intralineage turnover is faster than the neutral expectation, further supporting the existence of independent niches. The persistence of virus lineages appears to stem from limited outbreaks within small communities, so that only a small fraction of the global susceptible population is infected at any time. As disparate communities become increasingly connected through globalization, interaction and competition between lineages might increase as well, which could result in changing selective pressures and increased diversification and/or pathogenicity. Thus, in addition to zoonotic events, ongoing surveillance of familiar, endemic viruses appears to merit global attention with respect to the prevention or mitigation of future pandemics.

Thermal stress triggers productive viral infection of a key coral reef symbiont

Climate change-driven ocean warming is increasing the frequency and severity of bleaching events, in which corals appear whitened after losing their dinoflagellate endosymbionts (family Symbiodiniaceae). Viral infections of Symbiodiniaceae may contribute to some bleaching signs, but little empirical evidence exists to support this hypothesis. We present the first temporal analysis of a lineage of Symbiodiniaceae-infecting positive-sense single-stranded RNA viruses ("dinoRNAVs") in coral colonies, which were exposed to a 5-day heat treatment (+2.1 °C). A total of 124 dinoRNAV major capsid protein gene "aminotypes" (unique amino acid sequences) were detected from five colonies of two closely related Pocillopora-Cladocopium (coral-symbiont) combinations in the experiment; most dinoRNAV aminotypes were shared between the two coral-symbiont combinations (64%) and among multiple colonies (82%). Throughout the experiment, seventeen dinoRNAV aminotypes were found only in heat-treated fragments, and 22 aminotypes were detected at higher relative abundances in heat-treated fragments. DinoRNAVs in fragments of some colonies exhibited higher alpha diversity and dispersion under heat stress. Together, these findings provide the first empirical evidence that exposure to high temperatures triggers some dinoRNAVs to switch from a persistent to a productive infection mode within heat-stressed corals. Over extended time frames, we hypothesize that cumulative dinoRNAV production in the Pocillopora-Cladocopium system could affect colony symbiotic status, for example, by decreasing Symbiodiniaceae densities within corals. This study sets the stage for reef-scale investigations of dinoRNAV dynamics during bleaching events.

Global Evolutionary History and Dynamics of Dengue Viruses Inferred from Whole Genome Sequences

Dengue is an arboviral disease caused by dengue virus (DENV), leading to approximately 25,000 deaths/year and with over 40% of the world’s population at risk. Increased international travel and trade, poorly regulated urban expansion, and warming global temperatures have expanded the geographic range and incidence of the virus in recent decades. This study used phylogenetic and selection pressure analyses to investigate trends in DENV evolution, using whole genome coding sequences from publicly available databases alongside newly sequenced isolates collected between 1963–1997 from Southeast Asia and the Pacific. Results revealed very similar phylogenetic relationships when using the envelope gene and the whole genome coding sequences. Although DENV evolution is predominantly driven by negative selection, a number of amino acid sites undergoing positive selection were found across the genome, with the majority located in the envelope and NS5 genes. Some genotypes appear to be diversifying faster than others within each serotype. The results from this research improve our understanding of DENV evolution, with implications for disease control efforts such as Wolbachia-based biocontrol and vaccine design.

Purifying Selection Determines the Short-Term Time Dependency of Evolutionary Rates in SARS-CoV-2 and pH1N1 Influenza

High-throughput sequencing enables rapid genome sequencing during infectious disease outbreaks and provides an opportunity to quantify the evolutionary dynamics of pathogens in near real-time. One difficulty of undertaking evolutionary analyses over short timescales is the dependency of the inferred evolutionary parameters on the timespan of observation. Crucially, there are an increasing number of molecular clock analyses using external evolutionary rate priors to infer evolutionary parameters. However, it is not clear which rate prior is appropriate for a given time window of observation due to the time-dependent nature of evolutionary rate estimates. Here, we characterize the molecular evolutionary dynamics of SARS-CoV-2 and 2009 pandemic H1N1 (pH1N1) influenza during the first 12 months of their respective pandemics. We use Bayesian phylogenetic methods to estimate the dates of emergence, evolutionary rates, and growth rates of SARS-CoV-2 and pH1N1 over time and investigate how varying sampling window and data set sizes affect the accuracy of parameter estimation. We further use a generalized McDonald-Kreitman test to estimate the number of segregating nonneutral sites over time. We find that the inferred evolutionary parameters for both pandemics are time dependent, and that the inferred rates of SARS-CoV-2 and pH1N1 decline by ∼50% and ∼100%, respectively, over the course of 1 year. After at least 4 months since the start of sequence sampling, inferred growth rates and emergence dates remain relatively stable and can be inferred reliably using a logistic growth coalescent model. We show that the time dependency of the mean substitution rate is due to elevated substitution rates at terminal branches which are 2-4 times higher than those of internal branches for both viruses. The elevated rate at terminal branches is strongly correlated with an increasing number of segregating nonneutral sites, demonstrating the role of purifying selection in generating the time dependency of evolutionary parameters during pandemics.

OUP accepted manuscript

The Zika virus (ZIKV) disease caused a public health emergency of international concern that started in February 2016. The overall number of ZIKV-related cases increased until November 2016, after which it declined sharply. While the evaluation of the potential risk and impact of future arbovirus epidemics remains challenging, intensified surveillance efforts along with a scale-up of ZIKV whole-genome sequencing provide an opportunity to understand the patterns of genetic diversity, evolution, and spread of ZIKV. However, a classification system that reflects the true extent of ZIKV genetic variation is lacking. Our objective was to characterize ZIKV genetic diversity and phylodynamics, identify genomic footprints of differentiation patterns, and propose a dynamic classification system that reflects its divergence levels. We analysed a curated dataset of 762 publicly available sequences spanning the full-length coding region of ZIKV from across its geographical span and collected between 1947 and 2021. The definition of genetic groups was based on comprehensive evolutionary dynamics analyses, which included recombination and phylogenetic analyses, within- and between-group pairwise genetic distances comparison, detection of selective pressure, and clustering analyses. Evidence for potential recombination events was detected in a few sequences. However, we argue that these events are likely due to sequencing errors as proposed in previous studies. There was evidence of strong purifying selection, widespread across the genome, as also detected for other arboviruses. A total of 50 sites showed evidence of positive selection, and for a few of these sites, there was amino acid (AA) differentiation between genetic clusters. Two main genetic clusters were defined, ZA and ZB, which correspond to the already characterized ‘African’ and ‘Asian’ genotypes, respectively. Within ZB, two subgroups, ZB.1 and ZB.2, represent the Asiatic and the American (and Oceania) lineages, respectively. ZB.1 is further subdivided into ZB.1.0 (a basal Malaysia sequence sampled in the 1960s and a recent Indian sequence), ZB.1.1 (South-Eastern Asia, Southern Asia, and Micronesia sequences), and ZB.1.2 (very similar sequences from the outbreak in Singapore). ZB.2 is subdivided into ZB.2.0 (basal American sequences and the sequences from French Polynesia, the putative origin of South America introduction), ZB.2.1 (Central America), and ZB.2.2 (Caribbean and North America). This classification system does not use geographical references and is flexible to accommodate potential future lineages. It will be a helpful tool for studies that involve analyses of ZIKV genomic variation and its association with pathogenicity and serve as a starting point for the public health surveillance and response to on-going and future epidemics and to outbreaks that lead to the emergence of new variants.

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.

Estimating the age of the subfamily Orthocoronavirinae using host divergence times as calibration ages at two internal nodes

Viruses of the subfamily Orthocoronavirinae can cause mild to severe disease in people, including COVID-19, MERS and SARS. Their most common natural hosts are bat and bird species, which are mostly split across four virus genera. Molecular clock analyses of orthocoronaviruses suggested the most recent common ancestor of these viruses might have emerged either around 10,000 years ago or, using models accounting for selection, many millions of years. Here, we reassess the evolutionary history of these viruses. We present time-aware phylogenetic analyses of a RNA-dependent RNA polymerase locus from 123 orthocoronaviruses isolated from birds and bats, including those in New Zealand, which were geographically isolated from other bats around 35 million years ago. We used this age, as well as the age of the avian-mammals split, to calibrate the molecular clocks, under the assumption that these ages are applicable to the analyzed viruses. We found that the time to the most recent ancestor common for all orthocoronaviruses is likely 150 or more million years, supporting clock analyses that account for selection.

Ongoing global and regional adaptive evolution of SARS-CoV-2

Understanding the trends in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) evolution is paramount to control the COVID-19 pandemic. We analyzed more than 300,000 high-quality genome sequences of SARS-CoV-2 variants available as of January 2021. The results show that the ongoing evolution of SARS-CoV-2 during the pandemic is characterized primarily by purifying selection, but a small set of sites appear to evolve under positive selection. The receptor-binding domain of the spike protein and the region of the nucleocapsid protein associated with nuclear localization signals (NLS) are enriched with positively selected amino acid replacements. These replacements form a strongly connected network of apparent epistatic interactions and are signatures of major partitions in the SARS-CoV-2 phylogeny. Virus diversity within each geographic region has been steadily growing for the entirety of the pandemic, but analysis of the phylogenetic distances between pairs of regions reveals four distinct periods based on global partitioning of the tree and the emergence of key mutations. The initial period of rapid diversification into region-specific phylogenies that ended in February 2020 was followed by a major extinction event and global homogenization concomitant with the spread of D614G in the spike protein, ending in March 2020. The NLS-associated variants across multiple partitions rose to global prominence in March to July, during a period of stasis in terms of interregional diversity. Finally, beginning in July 2020, multiple mutations, some of which have since been demonstrated to enable antibody evasion, began to emerge associated with ongoing regional diversification, which might be indicative of speciation.

Natural selection on synonymous mutations in SARS-CoV-2 and the impact on estimating divergence time

To adapt to human host environment, synonymous mutations in SARS-CoV-2 are shaped by tRNA selection, energy cost and RNA structure.

Using Wolbachia to Eliminate Dengue: Will the Virus Fight Back?

Recent field trials have demonstrated that dengue incidence can be substantially reduced by introgressing strains of the endosymbiotic bacterium Wolbachia into Aedes aegypti mosquito populations. This strategy relies on Wolbachia reducing the susceptibility of Ae. aegypti to disseminated infection by positive-sense RNA viruses like dengue. However, RNA viruses are well known to adapt to antiviral pressures. Here, we review the viral infection stages where selection for Wolbachia-resistant virus variants could occur. We also consider the genetic constraints imposed on viruses that alternate between vertebrate and invertebrate hosts, and the likely selection pressures to which dengue virus might adapt in order to be effectively transmitted by Ae. aegypti that carry Wolbachia. While there are hurdles to dengue viruses developing resistance to Wolbachia, we suggest that long-term surveillance for resistant viruses should be an integral component of Wolbachia-introgression biocontrol programs.

Complete genome characterization and evolutionary analysis of dengue viruses isolated during 2016-2017 in Pune, India

Dengue is the most common mosquito-borne viral infection in tropical and sub-tropical countries. In the recent years, frequent dengue outbreaks are being reported in many parts of India. DENV circulates as four independent serotypes posing a major public health threat around the globe. Phylogenetic and full genome sequence analyses of 19 complete DENV genome sequences presenting all the four serotypes in Pune, India (2016-2017) revealed no change in the circulating genotypes i.e., genotype V clade C (D1), genotype IVB (D2), genotype III lineage III (D3) and genotype I clade D (D4). Additionally, unique amino acid substitutions that may potentially influence viral fitness and virulence in host cells were identified. Mapping of the unique amino acid substitutions onto the T cell epitopes of the reference strains revealed that 8/10 (D1), 14/15 (D2), 3/4 (D3) and 21/74 (D4), amino acids were involved in T-cell epitope presentation for a maximum number of HLA alleles associated with disease outcome. Selection pressure analysis documented a positive selection pressure to be acting on few amino acid sites indicating continuous evolutionary changes in the viral RNA. Overall, the evolutionary and selection pressure data generated during this study may help in better understanding of DENV evolution and epidemiology.

B-Cell Memory Responses to Variant Viral Antigens

A central feature of vertebrate immune systems is the ability to form antigen-specific immune memory in response to microbial challenge and so provide protection against future infection. In conflict with this process is the ability that many viruses have to mutate their antigens to escape infection- or vaccine-induced antibody memory responses. Mutable viruses such as dengue virus, influenza virus and of course coronavirus have a major global health impact, exacerbated by this ability to evade immune responses through mutation. There have been several outstanding recent studies on B-cell memory that also shed light on the potential and limitations of antibody memory to protect against viral antigen variation, and so promise to inform new strategies for vaccine design. For the purposes of this review, the current understanding of the different memory B-cell (MBC) populations, and their potential to recognize mutant antigens, will be described prior to some examples from antibody responses against the highly mutable RNA based flaviviruses, influenza virus and SARS-CoV-2.

Ongoing Global and Regional Adaptive Evolution of SARS-CoV-2

Understanding the trends in SARS-CoV-2 evolution is paramount to control the COVID-19 pandemic. We analyzed more than 300,000 high quality genome sequences of SARS-CoV-2 variants available as of January 2021. The results show that the ongoing evolution of SARS-CoV-2 during the pandemic is characterized primarily by purifying selection, but a small set of sites appear to evolve under positive selection. The receptor-binding domain of the spike protein and the nuclear localization signal (NLS) associated region of the nucleocapsid protein are enriched with positively selected amino acid replacements. These replacements form a strongly connected network of apparent epistatic interactions and are signatures of major partitions in the SARS-CoV-2 phylogeny. Virus diversity within each geographic region has been steadily growing for the entirety of the pandemic, but analysis of the phylogenetic distances between pairs of regions reveals four distinct periods based on global partitioning of the tree and the emergence of key mutations. The initial period of rapid diversification into region-specific phylogenies that ended in February 2020 was followed by a major extinction event and global homogenization concomitant with the spread of D614G in the spike protein, ending in March 2020. The NLS associated variants across multiple partitions rose to global prominence in March-July, during a period of stasis in terms of inter-regional diversity. Finally, beginning July 2020, multiple mutations, some of which have since been demonstrated to enable antibody evasion, began to emerge associated with ongoing regional diversification, which might be indicative of speciation.

Nonsynonymous Polymorphism Counts in Bacterial Genomes: a Comparative Examination

Genomic data reveal single-nucleotide polymorphisms (SNPs) that may carry information about the evolutionary history of bacteria. However, it remains unclear what inferences about selection can be made from genomic SNP data. Bacterial species are often sampled during epidemic outbreaks or within hosts during the course of chronic infections. SNPs obtained from genomic analysis of these data are not necessarily fixed. Treating them as fixed during analysis by using measures such as the ratio of nonsynonymous to synonymous evolutionary changes (dN/dS) may lead to incorrect inferences about the strength and direction of selection. In this study, we consider data from a range of whole-genome sequencing studies of bacterial pathogens and explore patterns of nonsynonymous variation to assess whether evidence of selection can be identified by investigating SNP counts alone across multiple WGS studies. We visualize these SNP data in ways that highlight their relationship to neutral baseline expectations. These neutral expectations are based on a simple model of mutation, from which we simulate SNP accumulation to investigate how SNP counts are distributed under alternative assumptions about positive and negative selection. We compare these patterns with empirical SNP data and illustrate the general difficulty of detecting positive selection from SNP data. Finally, we consider whether SNP counts observed at the between-host population level differ from those observed at the within-host level and find some evidence that suggests that dynamics across these two scales are driven by different underlying processes.IMPORTANCE Identifying selection from SNP data obtained from whole-genome sequencing studies is challenging. Some current measures used to identify and quantify selection acting on genomes rely on fixed differences; thus, these are inappropriate for SNP data where variants are not fixed. With the increase in whole-genome sequencing studies, it is important to consider SNP data in the context of evolutionary processes. How SNPs are counted and analyzed can help in understanding mutation accumulation and trajectories of strains. We developed a tool for identifying possible evidence of selection and for comparative analysis with other SNP data. We propose a model that provides a rule-of-thumb guideline and two new visualization techniques that can be used to interpret and compare SNP data. We quantify the expected proportion of nonsynonymous SNPs in coding regions under neutrality and demonstrate its use in identifying evidence of positive and negative selection from simulations and empirical data.

Evolution and epidemiologic dynamics of dengue virus in Nicaragua during the emergence of chikungunya and Zika viruses

Arthropod-borne viruses (arboviruses) comprise a significant and ongoing threat to human health, infecting hundreds of millions annually. Three such arboviruses include circumtropical dengue, Zika, and chikungunya viruses, exhibiting continuous emergence primarily via Aedes mosquito vectors. Nicaragua has experienced endemic dengue virus (DENV) transmission involving multiple serotypes since 1985, with chikungunya virus (CHIKV) reported in 2014-2015, followed by Zika virus (ZIKV) first reported in 2016. In order to identify patterns of genetic variation and selection pressures shaping the evolution of co-circulating DENV serotypes in light of the arrival of CHIKV and ZIKV, we employed whole-genome sequencing on an Illumina MiSeq platform of random-amplified total RNA libraries to characterize 42 DENV low-passage isolates, derived from viremic patients in Nicaragua between 2013 and 2016. Our approach also revealed clinically undetected co-infections with CHIKV. Of the three DENV serotypes (1, 2, and 3) co-circulating during our study, we uncovered distinct patterns of evolution using comparative phylogenetic inference. DENV-1 genetic variation was structured into two distinct co-circulating lineages with no evidence of positive selection in the origins of either lineage, suggesting they are equally fit. In contrast, the evolutionary history of DENV-2 was marked by positive selection, and a unique, divergent lineage correlated with high epidemic potential emerged in 2015 to drive an outbreak in 2016. DENV-3 genetic variation remained unstructured into lineages throughout the period of study. Thus, this study reveals insights into evolutionary and epidemiologic trends exhibited during the circulation of multiple arboviruses in Nicaragua.

Insights on Zika virus envelope gene conservation in American outbreaks

Phylogenetic studies with Zika virus (ZIKV) have been conducted in Brazil. In this study, we sequenced 8 new sequences of the ZIKV envelope (E) gene from strains of cases from the Paraná and Mato Grosso do Sul states in 2016. A low phylogenetic signal was observed, with more than 40% of unresolved quartets, and the Maximum Likelihood Tree grouped all sequences in the Brazilian branches within the Asian genotype. In addition, a Shannon entropy analysis was conducted, showing a high stability in the E protein through the ZIKV polyprotein. Taken together, these results suggest a high degree of conservation in the ZIKV E gene from the recent American outbreaks.

Biological Characteristics and Patterns of Codon Usage Evolution for the African Genotype Zika Virus

We investigated temporal trends of codon usage changes for different host species to determine their importance in Zika virus (ZIKV) evolution. Viral spillover resulting from the potential of codon adaptation to host genome was also assessed for the African genotype ZIKV in comparison to the Asian genotype. To improve our understanding on its zoonotic maintenance, we evaluated in vitro the biological properties of the African genotype ZIKV in vertebrate and mosquito cell lines. Analyses were performed in comparison to Yellow fever virus (YFV). Despite significantly lower codon adaptation index trends than YFV, ZIKV showed evident codon adaptation to vertebrate hosts, particularly for the green African monkey Chlorocebus aethiops. PCA and CAI analyses at the individual ZIKV gene level for both human and Aedes aegypti indicated a clear distinction between the two genotypes. African ZIKV isolates showed higher virulence in mosquito cells than in vertebrate cells. Their higher replication in mosquito cells than African YFV confirmed the role of mosquitoes in the natural maintenance of the African genotype ZIKV. An analysis of individual strain growth characteristics indicated that the widely used reference strain MR766 replicates poorly in comparison to African ZIKV isolates. The recombinant African Zika virus strain ArD128000*^E/NS5 may be a good model to include in studies on the mechanism of host tropism, as it cannot replicate in the tested vertebrate cell line.

Mobilization and Activation of the Innate Immune Response to Dengue Virus

Dengue virus is an important human pathogen, infecting an estimated 400 million individuals per year and causing symptomatic disease in a subset of approximately 100 million. Much of the effort to date describing the host response to dengue has focused on the adaptive immune response, in part because of the well-established roles of antibody-dependent enhancement and T cell original sin as drivers of severe dengue upon heterotypic secondary infection. However, the innate immune system is a crucial factor in the host response to dengue, as it both governs the fate and vigor of the adaptive immune response, and mediates the acute inflammatory response in tissues. In this review, we discuss the innate inflammatory response to dengue infection, focusing on the role of evolutionarily conserved innate immune cells, their effector functions, and clinical course.

The Genetic Diversification of a Single Bluetongue Virus Strain Using an In Vitro Model of Alternating-Host Transmission

Bluetongue virus (BTV) is an arbovirus that has been associated with dramatic epizootics in both wild and domestic ruminants in recent decades. As a segmented, double-stranded RNA virus, BTV can evolve via several mechanisms due to its genomic structure. However, the effect of BTV’s alternating-host transmission cycle on the virus’s genetic diversification remains poorly understood. Whole genome sequencing approaches offer a platform for investigating the effect of host-alternation across all ten segments of BTV’s genome. To understand the role of alternating hosts in BTV’s genetic diversification, a field isolate was passaged under three different conditions: (i) serial passages in Culicoides sonorensis cells, (ii) serial passages in bovine pulmonary artery endothelial cells, or (iii) alternating passages between insect and bovine cells. Aliquots of virus were sequenced, and single nucleotide variants were identified. Measures of viral population genetics were used to quantify the genetic diversification that occurred. Two consensus variants in segments 5 and 10 occurred in virus from all three conditions. While variants arose across all passages, measures of genetic diversity remained largely similar across cell culture conditions. Despite passage in a relaxed in vitro system, we found that this BTV isolate exhibited genetic stability across passages and conditions. Our findings underscore the valuable role that whole genome sequencing may play in improving understanding of viral evolution and highlight the genetic stability of BTV.

Population genomics of louping ill virus provide new insights into the evolution of tick-borne flaviviruses

The emergence and spread of tick-borne arboviruses pose an increased challenge to human and animal health. In Europe this is demonstrated by the increasingly wide distribution of tick-borne encephalitis virus (TBEV, Flavivirus, Flaviviridae), which has recently been found in the United Kingdom (UK). However, much less is known about other tick-borne flaviviruses (TBFV), such as the closely related louping ill virus (LIV), an animal pathogen which is endemic to the UK and Ireland, but which has been detected in other parts of Europe including Scandinavia and Russia. The emergence and potential spatial overlap of these viruses necessitates improved understanding of LIV genomic diversity, geographic spread and evolutionary history. We sequenced a virus archive composed of 22 LIV isolates which had been sampled throughout the UK over a period of over 80 years. Combining this dataset with published virus sequences, we detected no sign of recombination and found low diversity and limited evidence for positive selection in the LIV genome. Phylogenetic analysis provided evidence of geographic clustering as well as long-distance movement, including movement events that appear recent. However, despite genomic data and an 80-year time span, we found that the data contained insufficient temporal signal to reliably estimate a molecular clock rate for LIV. Additional analyses revealed that this also applied to TBEV, albeit to a lesser extent, pointing to a general problem with phylogenetic dating for TBFV. The 22 LIV genomes generated during this study provide a more reliable LIV phylogeny, improving our knowledge of the evolution of tick-borne flaviviruses. Our inability to estimate a molecular clock rate for both LIV and TBEV suggests that temporal calibration of tick-borne flavivirus evolution should be interpreted with caution and highlight a unique aspect of these viruses which may be explained by their reliance on tick vectors.

Tick-borne pathogens represent a major emerging threat to public health and in recent years have been expanding into new areas. LIV is a neglected virus endemic to the UK and Ireland (though it has been detected in Scandinavia and Russia) which is closely related to the major human pathogen TBEV, but predominantly causes disease in sheep and grouse. The recent detection of TBEV in the UK, which has also emerged elsewhere in Europe, requires more detailed understanding of the spread and sequence diversity of LIV. This could be important for diagnosis and vaccination, but also to improve our understanding of the evolution and emergence of these tick-borne viruses. Here we describe the sequencing of 22 LIV isolates which have been sampled from several host species across the past century. We have utilised this dataset to investigate the evolutionary pressures that LIV is subjected to and have explored the evolution of LIV using phylogenetic analysis. Crucially we were unable to estimate a reliable molecular clock rate for LIV and found that this problem also extends to a larger phylogeny of TBEV sequences. This work highlights a previously unknown caveat of tick-borne flavivirus evolutionary analysis which may be important for understanding the evolution of these important pathogens.

Application of Next-Generation Sequencing to Reveal How Evolutionary Dynamics of Viral Population Shape Dengue Epidemiology

Dengue viral (DENV) infection results in a wide spectrum of clinical manifestations from asymptomatic, mild fever to severe hemorrhage diseases upon infection. Severe dengue is the leading cause of pediatric deaths and/or hospitalizations, which are a major public health burden in dengue-endemic or hyperendemic countries. Like other RNA viruses, DENV continues to evolve. Adaptive mutations are obscured by the major consensus sequence (so-called wild-type sequences) and can only be identified once they become the dominant viruses in the virus population, a process that can take months or years. Traditional surveillance systems still rely on Sanger consensus sequencing. However, with the recent advancement of high-throughput next-generation sequencing (NGS) technologies, the genome-wide investigation of virus population within-host and between-hosts becomes achievable. Thus, viral population sequencing by NGS can increase our understanding of the changing epidemiology and evolution of viral genomics at the molecular level. This review focuses on the studies within the recent decade utilizing NGS in different experimental and epidemiological settings to understand how the adaptive evolution of dengue variants shapes the dengue epidemic and disease severity through its transmission. We propose three types of studies that can be pursued in the future to enhance our surveillance for epidemic prediction and better medical management.

Within-Host Viral Diversity: A Window into Viral Evolution

The evolutionary dynamics of a virus can differ within hosts and across populations. Studies of within-host evolution provide an important link between experimental studies of virus evolution and large-scale phylodynamic analyses. They can determine the extent to which global processes are recapitulated on local scales and how accurately experimental infections model natural ones. They may also inform epidemiologic models of disease spread and reveal how host-level dynamics contribute to a virus's evolution at a larger scale. Over the last decade, advances in viral sequencing have enabled detailed studies of viral genetic diversity within hosts. I review how within-host diversity is sampled, measured, and expressed, and how comparative studies of viral diversity can be leveraged to elucidate a virus's evolutionary dynamics. These concepts are illustrated with detailed reviews of recent research on the within-host evolution of influenza virus, dengue virus, and cytomegalovirus.

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.

Linking influenza virus evolution within and between human hosts

Influenza viruses rapidly diversify within individual human infections. Several recent studies have deep-sequenced clinical influenza infections to identify viral variation within hosts, but it remains unclear how within-host mutations fare at the between-host scale. Here, we compare the genetic variation of H3N2 influenza within and between hosts to link viral evolutionary dynamics across scales. Synonymous sites evolve at similar rates at both scales, indicating that global evolution at these putatively neutral sites results from the accumulation of within-host variation. However, nonsynonymous mutations are depleted between hosts compared to within hosts, suggesting that selection purges many of the protein-altering changes that arise within hosts. The exception is at antigenic sites, where selection detectably favors nonsynonymous mutations at the global scale, but not within hosts. These results suggest that selection against deleterious mutations and selection for antigenic change are the main forces that act on within-host variants of influenza virus as they transmit and circulate between hosts.

"Submergence" of Western equine encephalitis virus: Evidence of positive selection argues against genetic drift and fitness reductions

Understanding the circumstances under which arboviruses emerge is critical for the development of targeted control and prevention strategies. This is highlighted by the emergence of chikungunya and Zika viruses in the New World. However, to comprehensively understand the ways in which viruses emerge and persist, factors influencing reductions in virus activity must also be understood. Western equine encephalitis virus (WEEV), which declined during the late 20th century in apparent enzootic circulation as well as equine and human disease incidence, provides a unique case study on how reductions in virus activity can be understood by studying evolutionary trends and mechanisms. Previously, we showed using phylogenetics that during this period of decline, six amino acid residues appeared to be positively selected. To assess more directly the effect of these mutations, we utilized reverse genetics and competition fitness assays in the enzootic host and vector (house sparrows and Culex tarsalis mosquitoes). We observed that the mutations contemporary with reductions in WEEV circulation and disease that were non-conserved with respect to amino acid properties had a positive effect on enzootic fitness. We also assessed the effects of these mutations on virulence in the Syrian-Golden hamster model in relation to a general trend of increased virulence in older isolates. However, no change effect on virulence was observed based on these mutations. Thus, while WEEV apparently underwent positive selection for infection of enzootic hosts, residues associated with mammalian virulence were likely eliminated from the population by genetic drift or negative selection. These findings suggest that ecologic factors rather than fitness for natural transmission likely caused decreased levels of enzootic WEEV circulation during the late 20th century.

Considering Genomic and Immunological Correlates of Protection for a Dengue Intervention

Over three billion are at risk of dengue infection with more than 100 million a year presenting with symptoms that can lead to deadly haemorrhagic disease. There are however no treatments available and the only licensed vaccine shows limited efficacy and is able to enhance the disease in some cases. These failures have mainly been due to the complex pathology and lack of understanding of the correlates of protection for dengue virus (DENV) infection. With increasing data suggesting both a protective and detrimental effect for antibodies and CD8 T-cells whilst having complex environmental dynamics. This review discusses the roles of genomic and immunological aspects of DENV infection, providing both a historical interpretation and fresh discussion on how this information can be used for the next generation of dengue interventions.

A non-parametric analytic framework for within-host viral phylogenies and a test for HIV-1 founder multiplicity

Phylogenetics is a powerful tool for understanding the diversification dynamics of viral pathogens. Here we present an extension of the spectral density profile of the modified graph Laplacian, which facilitates the characterization of within-host molecular evolution of viruses and the direct comparison of diversification dynamics between hosts. This approach is non-parametric and therefore fast and model-free. We used simulations of within-host evolutionary scenarios to evaluate the efficiency of our approach and to demonstrate the significance of interpreting a viral phylogeny by its spectral density profile in terms of diversification dynamics. The key features that are captured by the profile are positive selection on the viral gene (or genome), temporal changes in substitution rates, mutational fitness, and time between sampling. Using sequences from individuals infected with HIV-1, we showed the utility of this approach for characterizing within-host diversification dynamics, for comparing dynamics between hosts, and for charting disease progression in infected individuals sampled over multiple years. We furthermore propose a heuristic test for assessing founder heterogeneity, which allows us to classify infections with single and multiple HIV-1 founder viruses. This non-parametric approach can be a valuable complement to existing parametric approaches.

Evolutionary timescale and geographical movement of cucumber mosaic virus, with focus on Iranian strains

Cucumber mosaic virus (CMV) is a geographically widespread plant virus with a very broad host range. The virus has been detected in diverse crops all over Iran. In this study, we estimated the timescale of the evolution of CMV subgroup I and the geographical movement of the virus with a focus on Iranian strains. Analyses using the MP and CP genes and their concatenation revealed that the CMV population within subgroup I had a single ancestor dating back to about 450-550 years ago. The Iranian strains formed three clusters in a maximum-clade-credibility phylogenetic tree. It was found that the most recent common ancestor of the Iranian strains within each cluster dates back to less than 100 years ago. Our results also suggest that both short- and long-distance migration of Iranian CMV strains has occurred in the last 100 years.

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.

Chikungunya Virus Vaccine Candidates with Decreased Mutational Robustness Are Attenuated In Vivo and Have Compromised Transmissibility

Chikungunya virus (CHIKV) is a reemerged arbovirus, a member of the Togaviridae family. It circulates through mosquito vectors mainly of the Aedes family and a mammalian host. CHIKV causes chikungunya fever, a mild to severe disease characterized by arthralgia, with some fatal outcomes described. In the past years, several outbreaks mainly caused by enhanced adaptation of the virus to the vector and ineffective control of the contacts between infected mosquito populations and the human host have been reported. Vaccines represent the best solution for the control of insect-borne viruses, including CHIKV, but are often unavailable. We designed live attenuated CHIKVs by applying a rational genomic design based on multiple replacements of synonymous codons. In doing so, the virus mutational robustness (capacity to maintain phenotype despite introduction of mutations to genotype) is decreased, driving the viral population toward deleterious evolutionary trajectories. When the candidate viruses were tested in the insect and mammalian hosts, we observed overall strong attenuation in both and greatly diminished signs of disease. Moreover, we found that the vaccine candidates elicited protective immunity related to the production of neutralizing antibodies after a single dose. During an experimental transmission cycle between mosquitoes and naive mice, vaccine candidates could be transmitted by mosquito bite, leading to asymptomatic infection in mice with compromised dissemination. Using deep-sequencing technology, we observed an increase in detrimental (stop) codons, which confirmed the effectiveness of this genomic design. Because the approach involves hundreds of synonymous modifications to the genome, the reversion risk is significantly reduced, rendering the viruses promising vaccine candidates.IMPORTANCE Chikungunya fever is a debilitating disease that causes severe pain to the joints, which can compromise the patient's lifestyle for several months and even in some grave cases lead to death. The etiological agent is chikungunya virus, an alphavirus transmitted by mosquito bite. Currently, there are no approved vaccines or treatments against the disease. In our research, we developed novel live attenuated vaccine candidates against chikungunya virus by applying an innovative genomic design. When tested in the insect and mammalian host, the vaccine candidates did not cause disease, elicited strong protection against further infection, and had low risk of reversion to pathogenic phenotypes.

The evolutionary dynamics of DENV 4 genotype I over a 60-year period

Dengue virus serotype 4 (DENV 4) has had a relatively low prevalence worldwide for decades; however, likely due to data paucity, no study has investigated the epidemiology and evolutionary dynamics of DENV 4 genotype I (DENV 4-I). This study aims to understand the diversity, epidemiology and dynamics of DENV 4-I. We collected 404 full length DENV4-1 envelope (E) gene sequences from 14 countries using two sources: Yunnan Province in China (15 strains during 2013–2016) and GenBank (489 strains up to 2018-01-11). Conducting phylogenetic and phylogeographical analyses, we estimated the virus spread, population dynamics, and selection pressures using different statistical analysis methods (substitution saturation, likelihood mapping, Bayesian coalescent inference, and maximum likelihood estimation). Our results show that during the last 60 years (1956–2016), DENV 4-I was present in mainland and maritime Southeast Asia, the Indian subcontinent, the southern provinces of China, parts of Brazil and Australia. The recent spread of DENV 4-I likely originated in the Philippines and later spread to Thailand. From Thailand, it spread to adjacent countries and eventually the Indian subcontinent. Apparently diverging around years 1957, 1963, 1976 and 1990, the different Clades (Clade I-V) were defined. The mean overall evolution rate of DENV 4-I was 9.74 (95% HPD: 8.68–10.82) × 10⁻⁴ nucleotide substitutions/site/year. The most recent common ancestor for DENV 4-I traces back to 1956. While the demographic history of DENV 4-I fluctuated, peaks appeared around 1982 and 2006. While purifying selection dominated the majority of E-gene evolution of DENV 4-I, positive selection characterized Clade III (Vietnam). DENV 4-I evolved in situ in Southeast Asia and the Indian subcontinent. Thailand and Indian acted as the main and secondary virus distribution hubs globally and regionally. Our phylogenetic analysis highlights the need for strengthened regional cooperation on surveillance and sharing of sample sequences to improve global dengue control and cross-border transmission prevention efforts.

Dengue virus (DENV) can be classified into four serotypes, DENV 1, 2, 3 and 4. Although DENV 4 is the first dengue serotype to diverge in phylogenetic analyses of the genus Flavivirus, this serotype occurs at a low prevalence worldwide and spreads the least rapidly. Similar to other serotypes, DENV 4 can also cause severe dengue (SD) disease manifestations, such as dengue haemorrhagic fever (DHF) and dengue shock syndrome (DSS). To date, no study has investigated the epidemiology and dynamics of DENV 4 genotype I comprehensively. In this study, we seek to address this gap. Our study shows that the distribution of DENV 4-I is mainly restricted to Southeast Asia and the Indian subcontinent. The most recent spread of DENV 4-I likely originated from Southeast Asia–initially circulating in the Philippines, then Thailand and later on the Indian subcontinent. Viruses evolved in situ in Southeast Asia and the Indian subcontinent, respectively. Although DENV 4-I occasionally spread elsewhere, this genotype did not become widely established. The overall evolution rate of DENV 4-I was comparable with that of DENV 2–4. The nucleotide sequences indicates that the demographic history of DENV 4-I fluctuated with peaks apparent during parts of the 1980s and 2000s. Although a weak positive selection existed in Clade III -predominately in Vietnam, purifying selection dominated the E-gene evolution of DENV 4-I.

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

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

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

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

Marine RNA Virus Quasispecies Are Distributed throughout the Oceans

Very little is known about aquatic RNA virus populations and genome evolution. This is the first study that analyzes marine environmental RNA viral assemblages in an evolutionary and broad geographical context. This study contributes the largest marine RNA virus metagenomic data set to date, substantially increasing the sequencing space for RNA viruses and also providing a baseline for comparisons of marine RNA virus diversity. The new viruses discovered in this study are representative of the most abundant family of marine RNA viruses, the Marnaviridae, and expand our view of the diversity of this important group. Overall, our data and analyses provide a foundation for interpreting marine RNA virus diversity and evolution.

RNA viruses, particularly genetically diverse members of the Picornavirales, are widespread and abundant in the ocean. Gene surveys suggest that there are spatial and temporal patterns in the composition of RNA virus assemblages, but data on their diversity and genetic variability in different oceanographic settings are limited. Here, we show that specific RNA virus genomes have widespread geographic distributions and that the dominant genotypes are under purifying selection. Genomes from three previously unknown picorna-like viruses (BC-1, -2, and -3) assembled from a coastal site in British Columbia, Canada, as well as marine RNA viruses JP-A, JP-B, and Heterosigma akashiwo RNA virus exhibited different biogeographical patterns. Thus, biotic factors such as host specificity and viral life cycle, and not just abiotic processes such as dispersal, affect marine RNA virus distribution. Sequence differences relative to reference genomes imply that virus quasispecies are under purifying selection, with synonymous single-nucleotide variations dominating in genomes from geographically distinct regions resulting in conservation of amino acid sequences. Conversely, sequences from coastal South Africa that mapped to marine RNA virus JP-A exhibited more nonsynonymous mutations, probably representing amino acid changes that accumulated over a longer separation. This biogeographical analysis of marine RNA viruses demonstrates that purifying selection is occurring across oceanographic provinces. These data add to the spectrum of known marine RNA virus genomes, show the importance of dispersal and purifying selection for these viruses, and indicate that closely related RNA viruses are pathogens of eukaryotic microbes across oceans.

IMPORTANCE Very little is known about aquatic RNA virus populations and genome evolution. This is the first study that analyzes marine environmental RNA viral assemblages in an evolutionary and broad geographical context. This study contributes the largest marine RNA virus metagenomic data set to date, substantially increasing the sequencing space for RNA viruses and also providing a baseline for comparisons of marine RNA virus diversity. The new viruses discovered in this study are representative of the most abundant family of marine RNA viruses, the Marnaviridae, and expand our view of the diversity of this important group. Overall, our data and analyses provide a foundation for interpreting marine RNA virus diversity and evolution.

Tracking dengue virus type 1 genetic diversity during lineage replacement in an hyperendemic area in Colombia

Dengue virus (DENV) is a flavivirus responsible for the most common and burdensome arthropod-borne viral disease of humans[1]. DENV evolution has been extensively studied on broad geographic and time scales, using sequences from a single gene[2,3]. It is believed that DENV evolution in humans is dominated primarily by purifying selection due to the constraint of maintaining fitness in both humans and mosquitoes[4,5]. Few studies have explored DENV evolutionary dynamics using whole genome sequences, nor have they explored changes in viral diversity that occur during intra-epidemic periods. We used deep sequencing of the viral coding region to characterize DENV-1 evolution in a Colombian population sampled during two high-prevalence dengue seasons in which serotype dominance shifted. Our data demonstrate patterns of strain extinction and replacement within DENV-1 as its prevalence waned and DENV-3 became established. A comparison of whole-genome versus single-gene-based phylogenetic analyses highlights an important difference in evolutionary patterns. We report a trend of higher nonsynonymous to synonymous diversity ratios among non-structural (NS) genes, and statistically significantly higher values among these ratios in the NS1 gene after DENV-1 strain replacement. These results suggest that positive selection could be driving DENV evolution within individual communities. Signals of positive selection coming from distinct samples may be drowned out when combining multiple regions with differing patterns of endemic transmission as commonly done by large-scale geo-temporal assessments. Here, we frame our findings within a small, local transmission history which aids significance. Moreover, these data suggest that the NS1 gene, rather than the E gene, may be a target of positive selection, although not mutually exclusive, and potentially useful sentinel of adaptive changes at the population level.

Effects of life history and ecology on virus evolutionary potential

The life history traits of viruses pose many consequences for viral population structure. In turn, population structure may influence the evolutionary trajectory of a virus. Here we review factors that affect the evolutionary potential of viruses, including rates of mutation and recombination, bottlenecks, selection pressure, and ecological factors such as the requirement for hosts and vectors. Mutation, while supplying a pool of raw genetic material, also results in the generation of numerous unfit mutants. The infection of multiple host species may expand a virus' ecological niche, although it may come at a cost to genetic diversity. Vector-borne viruses often experience a diminished frequency of positive selection and exhibit little diversity, and resistance against vector-borne viruses may thus be more durable than against non-vectored viruses. Evidence indicates that adaptation to a vector is more evolutionarily difficult than adaptation to a host. Overall, a better understanding of how various factors influence viral dynamics in both plant and animal pathosystems will lead to more effective anti-viral treatments and countermeasures.

Infectivity of Dengue Virus Serotypes 1 and 2 Is Correlated with E-Protein Intrinsic Dynamics but Not to Envelope Conformations

Dengue is a mosquito-borne virus with dire health and economic impacts. Dengue is responsible for an estimated 390 million infections per year, with dengue 2 (DENV2) being the most virulent strain among the four serotypes. Interestingly, it is also in strains of this serotype that temperature-dependent large-scale morphological changes, termed "breathing," have been observed. Although the structure of these morphologies has been solved to 3.5-Å resolution, the dynamics of the viral envelope are unknown. Here, we combine fluorescence and mass spectrometry with molecular dynamics simulations to provide insights into DENV2 (NGC strain) structural dynamics in comparison with DENV1 (PVP 159). We observe hitherto unseen conformational changes and structural dynamics of the DENV2 envelope that are influenced by both temperature and divalent cations. Our results show that for DENV2 and DENV1 the intrinsic dynamics, but not the specific morphologies, are correlated with viral infectivity.

The impact of truncation selection and diffusion on cooperation in spatial games

Natural Selection is frequently modelled via proportional selection where survival is proportional to the average payoff differential. There has been little attention devoted to modelling truncation selection where replicators below a threshold are culled and survivors reproduce. Here, we systematically explore truncation selection for two strategy games in a spatial setting. We employ two variations of truncation selection: independent, where the threshold is fixed; and dependent, where the proportion culled is fixed. Further, we explore the effects of diffusion with the algorithms: contest-diffusion-offspring (CDO), and diffusion-contest-offspring (DCO). CDO and DCO frequently facilitate and diminish cooperation, respectively. For independent truncation, there are three qualitative regimes determined by the payoff threshold: cooperation decreases as the threshold rises; polymorphisms are stable; and extinction is frequent. Further, an intermediate payoff to cooperators playing defectors can maximize cooperation for the DCO algorithm with a high payoff threshold. Dependent truncation affects games differently; lower levels reduce cooperation for the Hawk Dove game and increase it for the Stag Hunt, and higher levels produce the opposite effects. Comparing these truncation methods to proportional selection, we show how they impact the prevalence of cooperation.

Persistence and Intra-Host Genetic Evolution of Zika Virus Infection in Symptomatic Adults: A Special View in the Male Reproductive System

We followed the presence of Zika virus (ZIKV) in four healthy adults (two men and two women), for periods ranging from 78 to 298 days post symptom onset. The patients were evaluated regarding the presence of the virus in different body fluids (blood, saliva, urine and semen), development of immune responses (including antibodies, cytokines and chemokines), and virus genetic variation within samples collected from semen and urine during the infection course. The analysis was focused primarily on the two male patients who shed the virus for up to 158 days after the initial symptoms. ZIKV particles were detected in the spermatozoa cytoplasm and flagella, in immature sperm cells and could also be isolated from semen in cell culture, confirming that the virus is able to preserve integrity and infectivity during replication in the male reproductive system (MRS). Despite the damage caused by ZIKV infection within the MRS, our data showed that ZIKV infection did not result in infertility at least in one of the male patients. This patient was able to conceive a child after the infection. We also detected alterations in the male genital cytokine milieu, which could play an important role in the replication and transmission of the virus which could considerably increase the risk of ZIKV sexual spread. In addition, full genome ZIKV sequences were obtained from several samples (mainly semen), which allowed us to monitor the evolution of the virus within a patient during the infection course. We observed genetic changes over time in consensus sequences and lower frequency intra-host single nucleotide variants (iSNV), that suggested independent compartmentalization of ZIKV populations in the reproductive and urinary systems. Altogether, the present observations confirm the risks associated with the long-term replication and shedding of ZIKV in the MRS and help to elucidate patterns of intra-host genetic evolution during long term replication of the virus.

Whole genome sequencing, variant analysis, phylogenetics, and deep sequencing of Zika virus strains

The recent emergence of Zika virus (ZIKV) has been concentrated in the Caribbean, Southeastern United States, and South- and Central America; resulting in travel-based cases being reported around the globe. As multi-disciplinary collaborations are combatting the ZIKV outbreak, the need to validate the sequence of existing strains has become apparent. Here, we report high-quality sequence data for multiple ZIKV strains made publicly available through the National Institutes of Health- (NIH) funded biorepository, BEI Resources (www.beiresources.org). Next-generation sequencing, 3' rapid amplification of cDNA ends (RACE), and viral genome annotation pipelines generated GenBank sequence records for 16 BEI Resources strains. Minor variants, consensus mutations, and consensus insertions/deletions were identified within the viral stocks using next-generation sequencing (NGS) and consensus changes were confirmed with Sanger sequencing. Bioinformatics analyses of the sequencing results confirm that the virus stocks available to the scientific research community through BEI Resources adequately represent the viral population diversity of ZIKV.

Inter- and intra-host sequence diversity reveal the emergence of viral variants during an overwintering epidemic caused by dengue virus serotype 2 in southern Taiwan

Purifying selection during dengue viral infection has been suggested as the driving force of viral evolution and the higher complexity of the intra-host quasi-species is thought to offer an adaptive advantage for arboviruses as they cycle between arthropod and vertebrate hosts. However, very few studies have been performed to investigate the viral genetic changes within (intra-host) and between (inter-host) humans in a spatio-temporal scale. Viruses of different serotypes from various countries imported to Taiwan cause annual outbreaks. During 2001-2003, two consecutive outbreaks were caused by dengue virus serotype 2 (DENV-2) and resulted in a larger-scale epidemic with more severe dengue cases in the following year. Phylogenetic analyses showed that the viruses from both events were similar and related to the 2001 DENV-2 isolate from the Philippines. We comprehensively analyzed viral sequences from representative dengue patients and identified three consensus genetic variants, group Ia, Ib and II, with different spatio-temporal population dynamics. The phylodynamic analysis suggested group Ib variants, characterized by lower genetic diversity, transmission rate, and intra-host variant numbers, might play the role of maintenance variants. The residential locations among the patients infected by group Ib variants were in the outer rim of case clusters throughout the 2001-2003 period whereas group Ia and II variants were located in the centers of case clusters, suggesting that group Ib viruses might serve as "sheltered overwintering" variants in an undefined ecological niche. Further deep sequencing of the viral envelope (E) gene directly from individual patient serum samples confirmed the emergence of variants belonging to three quasi-species (group Ia, Ib, and II) and the ancestral role of the viral variants in the latter phase of the 2001 outbreak contributed to the later, larger-scale epidemic beginning in 2002. These findings enhanced our understanding of increasing epidemic severity over time in the same epidemic area. It also highlights the importance of combining phylodynamic and deep sequencing analysis as surveillance tools for detecting dynamic changes in viral variants, particularly searching for and monitoring any specific viral subpopulation. Such subpopulations might have selection advantages in both fitness and transmissibility leading to increased epidemic severity.