Intrahost Selection Pressures Drive Rapid Dengue Virus Microevolution in Acute Human Infections

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.

Pathophysiology of chikungunya virus infection associated with fatal outcomes

Chikungunya virus (CHIKV) is a mosquito-borne alphavirus that causes acute, subacute, and chronic human arthritogenic diseases and, in rare instances, can lead to neurological complications and death. Here, we combined epidemiological, virological, histopathological, cytokine, molecular dynamics, metabolomic, proteomic, and genomic analyses to investigate viral and host factors that contribute to chikungunya-associated (CHIK) death. Our results indicate that CHIK deaths are associated with multi-organ infection, central nervous system damage, and elevated serum levels of pro-inflammatory cytokines and chemokines compared with survivors. The histopathologic, metabolite, and proteomic signatures of CHIK deaths reveal hemodynamic disorders and dysregulated immune responses. The CHIKV East-Central-South-African lineage infecting our study population causes both fatal and survival cases. Additionally, CHIKV infection impairs the integrity of the blood-brain barrier, as evidenced by an increase in permeability and altered tight junction protein expression. Overall, our findings improve the understanding of CHIK pathophysiology and the causes of fatal infections.

Whole Genome Sequencing of DENV-2 isolated from Aedes aegypti mosquitoes in Esmeraldas, Ecuador. Genomic epidemiology of genotype III Southern Asian-American in the country

Ecuador is a tropical country reporting Dengue virus (DENV) outbreaks with areas of hyperendemic viral transmission. Entomo-virological surveillance and monitoring effort conducted in the Northwestern border province of Esmeraldas in April 2022, five pools of female Aedes aegypti mosquitoes from a rural community tested positive for DENV serotype 2 by RT-qPCR. One pool was sequenced by Illumina MiSeq, and it corresponded to genotype III Southern Asian-American. Comparison with other genomes revealed genetic similarity to a human DENV genome sequenced in 2021, also from Esmeraldas. Potential introduction events to the country could have originated from Colombia, considering the vicinity of the collection sites to the neighboring country and high human movement. The inclusion of genomic information complements entomo-virological surveillance, providing valuable insights into genetic variants. This contribution enhances our understanding of Dengue virus (DENV) epidemiology in rural areas and guides evidence-based decisions for surveillance and interventions.

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

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

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

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

Impact of Intrahost NS5 Nucleotide Variations on Dengue Virus Replication

Due to the nature of RNA viruses, their high mutation rates produce a population of closely related but genetically diverse viruses, termed quasispecies. To determine the role of quasispecies in DENV disease severity, 22 isolates (10 from mild cases, 12 from fatal cases) were obtained, amplified, and sequenced with Next Generation Sequencing using the Illumina MiSeq platform. Using variation calling, unique wildtype nucleotide positions were selected and analyzed for variant nucleotides between mild and fatal cases. The analysis of variant nucleotides between mild and fatal cases showed 6 positions with a significant difference of p < 0.05 with 1 position in the structural region, and 5 positions in the non-structural (NS) regions. All variations were found to have a higher percentage in fatal cases. To further investigate the genetic changes that affect the virus’s properties, reverse genetics (rg) viruses containing substitutions with the variations were generated and viral growth properties were examined. We found that the virus variant rgNS5-T7812G (G81G) had higher replication rates in both Baby hamster kidney cells (BHK-21) and Vero cells while rgNS5-C9420A (A617A) had a higher replication rate only in BHK-21 cells compared to wildtype virus. Both variants were considered temperature sensitive whereby the viral titers of the variants were relatively lower at 39°C, but was higher at 35 and 37°C. Additionally, the variants were thermally stable compared to wildtype at temperatures of 29, 37, and 39°C. In conclusion, viral quasispecies found in isolates from the 2015 DENV epidemic, resulted in variations with significant difference between mild and fatal cases. These variations, NS5-T7812G (G81G) and NS5-C9420A (A617A), affect viral properties which may play a role in the virulence of DENV.

Development of a Singleplex Real-Time Reverse Transcriptase PCR Assay for Pan-Dengue Virus Detection and Quantification

Dengue virus (DENV) infection is a significant global health problem. There are no specific therapeutics or widely available vaccines. Early diagnosis is critical for patient management. Viral RNA detection by multiplex RT-PCR using multiple pairs of primers/probes allowing the simultaneous detection of all four DENV serotypes is commonly used. However, increasing the number of primers in the RT-PCR reaction reduces the sensitivity of detection due to the increased possibility of primer dimer formation. Here, a one tube, singleplex real-time RT-PCR specific to DENV 3′-UTR was developed for the detection and quantification of pan-DENV with no cross reactivity to other flaviviruses. The sensitivity of DENV detection was as high as 96.9% in clinical specimens collected at the first day of hospitalization. Our assay provided equivalent PCR efficiency and RNA quantification among each DENV serotype. The assay’s performance was comparable with previously established real-time RT-PCR targeting coding sequences. Using both assays on the same specimens, our results indicate the presence of defective virus particles in the circulation of patients infected with all serotypes. Dual regions targeting RT-PCR enhanced the sensitivity of viral genome detection especially during the late acute phase when viremia rapidly decline and an incomplete viral genome was clinically evident.

Assessment of the Risk of Severe Dengue Using Intrahost Viral Population in Dengue Virus Serotype 2 Patients via Machine Learning

Dengue virus, a positive-sense single-stranded RNA virus, continuously threatens human health. Although several criteria for evaluation of severe dengue have been recently established, the ability to prognose the risk of severe outcomes for dengue patients remains limited. Mutant spectra of RNA viruses, including single nucleotide variants (SNVs) and defective virus genomes (DVGs), contribute to viral virulence and growth. Here, we determine the potency of intrahost viral population in dengue patients with primary infection that progresses into severe dengue. A total of 65 dengue virus serotype 2 infected patients in primary infection including 17 severe cases were enrolled. We utilized deep sequencing to directly define the frequency of SNVs and detection times of DVGs in sera of dengue patients and analyzed their associations with severe dengue. Among the detected SNVs and DVGs, the frequencies of 9 SNVs and the detection time of 1 DVG exhibited statistically significant differences between patients with dengue fever and those with severe dengue. By utilizing the detected frequencies/times of the selected SNVs/DVG as features, the machine learning model showed high average with a value of area under the receiver operating characteristic curve (AUROC, 0.966 ± 0.064). The elevation of the frequency of SNVs at E (nucleotide position 995 and 2216), NS2A (nucleotide position 4105), NS3 (nucleotide position 4536, 4606), and NS5 protein (nucleotide position 7643 and 10067) and the detection times of the selected DVG that had a deletion junction in the E protein region (nucleotide positions of the junction: between 969 and 1022) increased the possibility of dengue patients for severe dengue. In summary, we demonstrated the detected frequencies/times of SNVs/DVG in dengue patients associated with severe disease and successfully utilized them to discriminate severe patients using machine learning algorithm. The identified SNVs and DVGs that are associated with severe dengue will expand our understanding of intrahost viral population in dengue pathogenesis.

OUP accepted manuscript

A detailed understanding of how and when severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission occurs is crucial for designing effective prevention measures. Other than contact tracing, genome sequencing provides information to help infer who infected whom. However, the effectiveness of the genomic approach in this context depends on both (high enough) mutation and (low enough) transmission rates. Today, the level of resolution that we can obtain when describing SARS-CoV-2 outbreaks using just genomic information alone remains unclear. In order to answer this question, we sequenced forty-nine SARS-CoV-2 patient samples from ten local clusters in NW Spain for which partial epidemiological information was available and inferred transmission history using genomic variants. Importantly, we obtained high-quality genomic data, sequencing each sample twice and using unique barcodes to exclude cross-sample contamination. Phylogenetic and cluster analyses showed that consensus genomes were generally sufficient to discriminate among independent transmission clusters. However, levels of intrahost variation were low, which prevented in most cases the unambiguous identification of direct transmission events. After filtering out recurrent variants across clusters, the genomic data were generally compatible with the epidemiological information but did not support specific transmission events over possible alternatives. We estimated the effective transmission bottleneck size to be one to two viral particles for sample pairs whose donor–recipient relationship was likely. Our analyses suggest that intrahost genomic variation in SARS-CoV-2 might be generally limited and that homoplasy and recurrent errors complicate identifying shared intrahost variants. Reliable reconstruction of direct SARS-CoV-2 transmission based solely on genomic data seems hindered by a slow mutation rate, potential convergent events, and technical artifacts. Detailed contact tracing seems essential in most cases to study SARS-CoV-2 transmission at high resolution.

Intra- and inter-host evolution of H9N2 influenza A virus in Japanese quail

Influenza A viruses (IAVs) are constantly evolving. Crucial steps in the infection cycle, such as sialic acid (SA) receptor binding on the host cell surface, can either promote or hamper the emergence of new variants. We previously assessed the relative fitness in Japanese quail of H9N2 variant viruses differing at a single amino acid position, residue 216 in the hemagglutinin (HA) viral surface protein. This site is known to modulate SA recognition. Our prior study generated a valuable set of longitudinal samples from quail transmission groups where the inoculum comprised different mixed populations of HA 216 variant viruses. Here, we leveraged these samples to examine the evolutionary dynamics of viral populations within and between inoculated and naïve contact quails. We found that positive selection dominated HA gene evolution, but fixation of the fittest variant depended on the competition mixture. Analysis of the whole genome revealed further evidence of positive selection acting both within and between hosts. Positive selection drove fixation of variants in non-HA segments within inoculated and contact quails. Importantly, transmission bottlenecks were modulated by the molecular signature at HA 216, revealing viral receptor usage as a determinant of transmitted diversity. Overall, we show that selection strongly shaped the evolutionary dynamics within and between quails. These findings support the notion that selective processes act effectively on IAV populations in poultry hosts, facilitating rapid viral evolution in this ecological niche.

Dynamics of SARS-CoV-2 mutations reveals regional-specificity and similar trends of N501 and high-frequency mutation N501Y in different levels of control measures

Coronavirus disease 2019 (COVID-19) is a contagious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This disease has spread globally, causing more than 161.5 million cases and 3.3 million deaths to date. Surveillance and monitoring of new mutations in the virus' genome are crucial to our understanding of the adaptation of SARS-CoV-2. Moreover, how the temporal dynamics of these mutations is influenced by control measures and non-pharmaceutical interventions (NPIs) is poorly understood. Using 1,058,020 SARS-CoV-2 from sequenced COVID-19 cases from 98 countries (totaling 714 country-month combinations), we perform a normalization by COVID-19 cases to calculate the relative frequency of SARS-CoV-2 mutations and explore their dynamics over time. We found 115 mutations estimated to be present in more than 3% of global COVID-19 cases and determined three types of mutation dynamics: high-frequency, medium-frequency, and low-frequency. Classification of mutations based on temporal dynamics enable us to examine viral adaptation and evaluate the effects of implemented control measures in virus evolution during the pandemic. We showed that medium-frequency mutations are characterized by high prevalence in specific regions and/or in constant competition with other mutations in several regions. Finally, taking N501Y mutation as representative of high-frequency mutations, we showed that level of control measure stringency negatively correlates with the effective reproduction number of SARS-CoV-2 with high-frequency or not-high-frequency and both follows similar trends in different levels of stringency.

Genetic diversity through social heterosis can increase virulence in RNA viral infections and cancer progression

In viral infections and cancer tumours, negative health outcomes often correlate with increasing genetic diversity. Possible evolutionary processes for such relationships include mutant lineages escaping host control or diversity, per se, creating too many immune system targets. Another possibility is social heterosis where mutations and replicative errors create clonal lineages varying in intrinsic capability for successful dispersal; improved environmental buffering; resource extraction or effective defence against immune systems. Rather than these capabilities existing in one genome, social heterosis proposes complementary synergies occur across lineages in close proximity. Diverse groups overcome host defences as interacting 'social genomes' with group genetic tool kits exceeding limited individual plasticity. To assess the possibility of social heterosis in viral infections and cancer progression, we conducted extensive literature searches for examples consistent with general and specific predictions from the social heterosis hypothesis. Numerous studies found supportive patterns in cancers across multiple tissues and in several families of RNA viruses. In viruses, social heterosis mechanisms probably result from long coevolutionary histories of competition between pathogen and host. Conversely, in cancers, social heterosis is a by-product of recent mutations. Investigating how social genomes arise and function in viral quasi-species swarms and cancer tumours may lead to new therapeutic approaches.

Dengue Virus Serotype 2 Intrahost Diversity in Patients with Different Clinical Outcomes

Intrahost genetic diversity is thought to facilitate arbovirus adaptation to changing environments and hosts, and it might also be linked to viral pathogenesis. Dengue virus serotype 2 (DENV-2) has circulated in Brazil since 1990 and is associated with severe disease and explosive outbreaks. Intending to shed light on the viral determinants for severe dengue pathogenesis, we sought to analyze the DENV-2 intrahost genetic diversity in 68 patient cases clinically classified as dengue fever (n = 31), dengue with warning signs (n = 19), and severe dengue (n = 18). Unlike previous DENV intrahost diversity studies whose approaches employed PCR, here we performed viral whole-genome deep sequencing from clinical samples with an amplicon-free approach, representing the real intrahost diversity scenario. Striking differences were detected in the viral population structure between the three clinical categories, which appear to be driven mainly by different infection times and selection pressures, rather than being linked with the clinical outcome itself. Diversity in the NS2B gene, however, showed to be constrained, irrespective of clinical outcome and infection time. Finally, 385 non-synonymous intrahost single-nucleotide variants located along the viral polyprotein, plus variants located in the untranslated regions, were consistently identified among the samples. Of them, 124 were exclusively or highly detected among cases with warning signs and among severe cases. However, there was no variant that by itself appeared to characterize the cases of greater severity, either due to its low intrahost frequency or the conservative effect on amino acid substitution. Although further studies are necessary to determine their real effect on viral proteins, this heightens the possibility of epistatic interactions. The present analysis represents an initial effort to correlate DENV-2 genetic diversity to its pathogenic potential and thus contribute to understanding the virus’s dynamics within its human host.

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.

Antibody-Dependent Enhancement: A Challenge for Developing a Safe Dengue Vaccine

In 2019, the United States Food and Drug Administration accorded restricted approval to Sanofi Pasteur's Dengvaxia, a live attenuated vaccine (LAV) for dengue fever, a mosquito-borne viral disease, caused by four antigenically distinct dengue virus serotypes (DENV 1-4). The reason for this limited approval is the concern that this vaccine sensitized some of the dengue-naïve recipients to severe dengue fever. Recent knowledge about the nature of the immune response elicited by DENV viruses suggests that all LAVs have inherent capacity to predominantly elicit antibodies (Abs) against the pre-membrane (prM) and fusion loop epitope (FLE) of DENV. These antibodies are generally cross-reactive among DENV serotypes carrying a higher risk of promoting Antibody-Dependent Enhancement (ADE). ADE is a phenomenon in which suboptimal neutralizing or non-neutralizing cross-reactive antibodies bind to virus and facilitate Fcγ receptor mediated enhanced entry into host cells, followed by its replication, and thus increasing the cellular viral load. On the other hand, antibody responses directed against the host-cell receptor binding domain of DENV envelope domain-III (EDIII), exhibit a higher degree of type-specificity with lower potential of ADE. The challenges associated with whole DENV-based vaccine strategies necessitate re-focusing our attention toward the designed dengue vaccine candidates, capable of inducing predominantly type-specific immune responses. If the designed vaccines elicited predominantly EDIII-directed serotype specific antibodies in the absence of prM and FLE antibodies, this could avoid the ADE phenomenon largely associated with the prM and FLE antibodies. The generation of type-specific antibodies to each of the four DENV serotypes by the designed vaccines could avoid the immune evasion mechanisms of DENVs. For the enhanced vaccine safety, all dengue vaccine candidates should be assessed for the extent of type-specific (minimal ADE) vs. cross-reactive (ADE promoting) neutralizing antibodies. The type-specific EDIII antibodies may be more directly related to protection from disease in the absence of ADE promoted by the cross-reactive antibodies.

Dengue Infections in Colombia: Epidemiological Trends of a Hyperendemic Country

Dengue is a major public health problem in hyperendemic countries like Colombia, the understanding of the epidemiological trends is important for the development of efficient public health policies. We conducted a systematic review of the epidemiologic data on dengue in Colombia from 1971 to 2020. A total of 375 relevant citations were identified, 36 of which fulfilled the inclusion criteria. The data of dengue and severe dengue cases, infection fatality rate, and serotype distribution were used to understand and identify gaps in the epidemiological knowledge in Colombia. The epidemiology of dengue in this country was characterized by five main outbreaks in 1998, 2002, 2010, 2013, and 2019 with high fatality rates in comparison with the average values reported in the Americas. The case fatality rate of severe dengue exceeded 2% and all four serotypes co-circulate throughout the country with some regional variations. Overall, the behavior of dengue in Colombia is influenced by multiple factors including seasonal temperature variation and socioeconomic conditions. Additionally, the most important barriers in the epidemiological surveillance of dengue may be due to the insufficient notification rate in some regions and the low active search for the circulation of different serotypes.

Role of cytotoxic T lymphocytes and interferon-γ in coronavirus infection: Lessons from murine coronavirus infections in mice

Murine coronavirus (CoV) is a beta-CoV that infects mice by binding to carcinoembryonic antigen-related cell adhesion molecule 1. Intraperitoneal infection with the murine CoV strain JHM (JHMV) induces acute mild hepatitis in mice. While both innate and acquired immune responses play a significant role in the protection against murine CoV infection in mice, CD8⁺ cytotoxic T lymphocytes (CTLs) and interferon-γ are essential for viral clearance in JHMV-induced hepatitis. In addition, CoVs are characterized by high diversity, caused by mutations, recombination, and gene gain/loss. 25V16G is an immune-escape JHMV variant, which lacks a dominant CTL epitope. By evading immune responses, 25V16G establishes persistent infections, leading to granulomatous serositis in interferon-γ-deficient mice. These examples of CoV-associated pathogenesis in mice might provide useful information on other CoV infections, including coronavirus disease 2019 (COVID-19).

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.

Analysis of cell-associated DENV RNA by oligo(dT) primed 5' capture scRNAseq

Dengue is one of the most widespread vector-borne viral diseases in the world. However, the size, heterogeneity, and temporal dynamics of the cell-associated viral reservoir during acute dengue virus (DENV) infection remains unclear. In this study, we analyzed cells infected in vitro with DENV and PBMC from an individual experiencing a natural DENV infection utilizing 5’ capture single cell RNA sequencing (scRNAseq). Both positive- and negative-sense DENV RNA was detected in reactions containing either an oligo(dT) primer alone, or in reactions supplemented with a DENV-specific primer. The addition of a DENV-specific primer did not increase the total amount of DENV RNA captured or the fraction of cells identified as containing DENV RNA. However, inclusion of a DENV-specific cDNA primer did increase the viral genome coverage immediately 5’ to the primer binding site. Furthermore, while the majority of intracellular DENV sequence captured in this analysis mapped to the 5’ end of the viral genome, distinct patterns of enhanced coverage within the DENV polyprotein coding region were observed. The 5’ capture scRNAseq analysis of PBMC not only recapitulated previously published reports by detecting virally infected memory and naïve B cells, but also identified cell-associated genomic variants not observed in contemporaneous serum samples. These results demonstrate that oligo(dT) primed 5’ capture scRNAseq can detect DENV RNA and quantify virus-infected cells in physiologically relevant conditions, and provides insight into viral sequence variability within infected cells.

DEN-IM: dengue virus genotyping from amplicon and shotgun metagenomic sequencing

Dengue virus (DENV) represents a public health threat and economic burden in affected countries. The availability of genomic data is key to understanding viral evolution and dynamics, supporting improved control strategies. Currently, the use of high-throughput sequencing (HTS) technologies, which can be applied both directly to patient samples (shotgun metagenomics) and to PCR-amplified viral sequences (amplicon sequencing), is potentially the most informative approach to monitor viral dissemination and genetic diversity by providing, in a single methodological step, identification and characterization of the whole viral genome at the nucleotide level. Despite many advantages, these technologies require bioinformatics expertise and appropriate infrastructure for the analysis and interpretation of the resulting data. In addition, the many software solutions available can hamper the reproducibility and comparison of results. Here we present DEN-IM, a one-stop, user-friendly, containerized and reproducible workflow for the analysis of DENV short-read sequencing data from both amplicon and shotgun metagenomics approaches. It is able to infer the DENV coding sequence (CDS), identify the serotype and genotype, and generate a phylogenetic tree. It can easily be run on any UNIX-like system, from local machines to high-performance computing clusters, performing a comprehensive analysis without the requirement for extensive bioinformatics expertise. Using DEN-IM, we successfully analysed two types of DENV datasets. The first comprised 25 shotgun metagenomic sequencing samples from patients with variable serotypes and genotypes, including an in vitro spiked sample containing the four known serotypes. The second consisted of 106 paired-end and 76 single-end amplicon sequences of DENV 3 genotype III and DENV 1 genotype I, respectively, where DEN-IM allowed detection of the intra-genotype diversity. The DEN-IM workflow, parameters and execution configuration files, and documentation are freely available at https://github.com/B-UMMI/DEN-IM).

Differential replicative fitness of the four dengue virus serotypes circulating in Colombia in human liver Huh7 cells

Dengue has been a significant public health problem in Colombia since the simultaneous circulation of the four dengue virus serotypes. The replicative fitness of dengue is a biological feature important for virus evolution and contributes to elucidating the behavior of virus populations and viral pathogenesis. However, it has not yet been studied in Colombian isolates. This study aimed to compare the replicative fitness of the four dengue virus serotypes and understand the association between the serotypes, their in vitro infection ability, and their replication in target cells. We used three isolates of each DENV serotype to infect Huh-7 cells at an MOI of 0.5. The percentage of infected cells was evaluated by flow cytometry, cell viability was evaluated by MTT assay, and the pathogenicity index was calculated as a ratio of both parameters. The replicative fitness was measured by the number of viral genome copies produced using quantitative PCR and the production of infectious viral progeny was measured by plaque assay. We showed that Huh-7 cells were susceptible to infection with all the different strain isolates. Nevertheless, the biological characteristics, such as infectious ability and cell viability, were strain-dependent. We also found different degrees of pathogenicity between strains of the four serotypes, representative of the heterogeneity displayed in the circulating population. When we analyzed the replicative fitness using the mean values obtained from RT-qPCR and plaque assay for the different strains, we found serotype-dependent behavior. The highest mean values of replicative fitness were obtained for DENV-1 (log 4.9 PFU/ml) and DENV-4 (log 5.28 PFU/ml), followed by DENV-2 (log 3.9 PFU/ml) and DENV-3 (log 4.31 PFU/ml). The internal heterogeneity of the replicative fitness within each serotype could explain the simultaneous circulation of the four DENV serotypes in Colombia.

Impact of RNA Virus Evolution on Quasispecies Formation and Virulence

RNA viruses are known to replicate by low fidelity polymerases and have high mutation rates whereby the resulting virus population tends to exist as a distribution of mutants. In this review, we aim to explore how genetic events such as spontaneous mutations could alter the genomic organization of RNA viruses in such a way that they impact virus replications and plaque morphology. The phenomenon of quasispecies within a viral population is also discussed to reflect virulence and its implications for RNA viruses. An understanding of how such events occur will provide further evidence about whether there are molecular determinants for plaque morphology of RNA viruses or whether different plaque phenotypes arise due to the presence of quasispecies within a population. Ultimately this review gives an insight into whether the intrinsically high error rates due to the low fidelity of RNA polymerases is responsible for the variation in plaque morphology and diversity in virulence. This can be a useful tool in characterizing mechanisms that facilitate virus adaptation and evolution.

Evolution of Attenuation and Risk of Reversal in Peste des Petits Ruminants Vaccine Strain Nigeria 75/1

Peste des Petits Ruminants (PPR) is a highly infectious disease caused by a virus of the Morbillivirus genus. The current PPR eradication effort relies mainly on the implementation of massive vaccination campaigns. One of the most widely used PPR vaccines is the Nigeria 75/1 strain obtained after attenuation by 75 serial passages of the wild type isolate in cell cultures. Here we use high throughput deep sequencing of the historical passages that led to the Nigeria 75/1 attenuated strain to understand the evolution of PPRV attenuation and to assess the risk of reversal in different cell types. Comparison of the consensus sequences of the wild type and vaccine strains showed that only 18 fixed mutations separate the two strains. At the earliest attenuation passage at our disposal (passage 47), 12 out of the 18 mutations were already present at a frequency of 100%. Low-frequency variants were identified along the genome in all passages. Sequencing of passages after the vaccine strain showed evidence of genetic drift during cell passages, especially in cells expressing the SLAM receptor targeted by PPRV. However, 15 out of the 18 mutations related to attenuation remained fixed in the population. In vitro experiments suggest that one mutation in the leader region of the PPRV genome affects virus replication. Our results suggest that only a few mutations can have a serious impact on the pathogenicity of PPRV. Risk of reversion to virulence of the attenuated PPRV strain Nigeria 75/1 during serial passages in cell cultures seems low but limiting the number of passages during vaccine production is recommended.

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

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

Intrahost Selection Pressure Drives Equine Arteritis Virus Evolution during Persistent Infection in the Stallion Reproductive Tract

Equine arteritis virus (EAV) is the causative agent of equine viral arteritis (EVA), a reproductive and respiratory disease of horses. Following natural infection, 10 to 70% of infected stallions can become carriers of EAV and continue to shed virus in the semen. In this study, sequential viruses isolated from nasal secretions, buffy coat cells, and semen of seven experimentally infected and two naturally infected EAV carrier stallions were deep sequenced to elucidate the intrahost microevolutionary process after a single transmission event. Analysis of variants from nasal secretions and buffy coat cells lacked extensive positive selection; however, characteristics of the mutant spectra were different in the two sample types. In contrast, the initial semen virus populations during acute infection have undergone a selective bottleneck, as reflected by the reduction in population size and diversifying selection at multiple sites in the viral genome. Furthermore, during persistent infection, extensive genome-wide purifying selection shaped variant diversity in the stallion reproductive tract. Overall, the nonstochastic nature of EAV evolution during persistent infection was driven by active intrahost selection pressure. Among the open reading frames within the viral genome, ORF3, ORF5, and the nsp2-coding region of ORF1a accumulated the majority of nucleotide substitutions during persistence, with ORF3 and ORF5 having the highest intrahost evolutionary rates. The findings presented here provide a novel insight into the evolutionary mechanisms of EAV and identified critical regions of the viral genome likely associated with the establishment and maintenance of persistent infection in the stallion reproductive tract.IMPORTANCE EAV can persist in the reproductive tract of infected stallions, and consequently, long-term carrier stallions constitute its sole natural reservoir. Previous studies demonstrated that the ampullae of the vas deferens are the primary site of viral persistence in the stallion reproductive tract and the persistence is associated with a significant inflammatory response that is unable to clear the infection. This is the first study that describes EAV full-length genomic evolution during acute and long-term persistent infection in the stallion reproductive tract using next-generation sequencing and contemporary sequence analysis techniques. The data provide novel insight into the intrahost evolution of EAV during acute and persistent infection and demonstrate that persistent infection is characterized by extensive genome-wide purifying selection and a nonstochastic evolutionary pattern mediated by intrahost selective pressure, with important nucleotide substitutions occurring in ORF1a (region encoding nsp2), ORF3, and ORF5.

Application of deep sequencing methods for inferring viral population diversity

The first deep sequencing method was announced in 2005. Due to an increasing number of sequencing data and a reduction in the costs of each sequencing dataset, this innovative technique was soon applied to genetic investigations of viral genome diversity in various viruses, particularly RNA viruses. These deep sequencing findings documented viral epidemiology and evolution and provided high-resolution data on the genetic changes in viral populations. Here, we review deep sequencing platforms that have been applied in viral quasispecies studies. Further, we discuss recent deep sequencing studies on viral inter- and intrahost evolution, drug resistance, and humoral immune selection, especially in emerging and re-emerging viruses. Deep sequencing methods are becoming the standard for providing comprehensive results of viral population diversity, and their applications are discussed.

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

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

Evolutionary Virology at 40

RNA viruses are diverse, abundant, and rapidly evolving. Genetic data have been generated from virus populations since the late 1970s and used to understand their evolution, emergence, and spread, culminating in the generation and analysis of many thousands of viral genome sequences. Despite this wealth of data, evolutionary genetics has played a surprisingly small role in our understanding of virus evolution. Instead, studies of RNA virus evolution have been dominated by two very different perspectives, the experimental and the comparative, that have largely been conducted independently and sometimes antagonistically. Here, we review the insights that these two approaches have provided over the last 40 years. We show that experimental approaches using in vitro and in vivo laboratory models are largely focused on short-term intrahost evolutionary mechanisms, and may not always be relevant to natural systems. In contrast, the comparative approach relies on the phylogenetic analysis of natural virus populations, usually considering data collected over multiple cycles of virus-host transmission, but is divorced from the causative evolutionary processes. To truly understand RNA virus evolution it is necessary to meld experimental and comparative approaches within a single evolutionary genetic framework, and to link viral evolution at the intrahost scale with that which occurs over both epidemiological and geological timescales. We suggest that the impetus for this new synthesis may come from methodological advances in next-generation sequencing and metagenomics.

Mechanisms and Concepts in RNA Virus Population Dynamics and Evolution

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

Complete Genome Characterization of the 2017 Dengue Outbreak in Xishuangbanna, a Border City of China, Burma and Laos

A dengue outbreak abruptly occurred at the border of China, Myanmar, and Laos in June 2017. By November 3rd 2017, 1184 infected individuals were confirmed as NS1-positivein Xishuangbanna, a city located at the border. To verify the causative agent, complete genome information was obtained through PCR and sequencing based on the viral RNAs extracted from patient samples. Phylogenetic trees were constructed by the maximum likelihood method (MEGA 6.0). Nucleotide and amino acid substitutions were analyzed by BioEdit, followed by RNA secondary structure prediction of untranslated regions (UTRs) and protein secondary structure prediction in coding sequences (CDSs). Strains YN2, YN17741, and YN176272 were isolated from local residents. Stains MY21 and MY22 were isolated from Burmese travelers. The complete genome sequences of the five isolates were 10,735 nucleotides in length. Phylogenetic analysis classified all five isolates as genotype I of DENV-1, while isolates of local residents and Burmese travelers belonged to different branches. The three locally isolates were most similar to the Dongguan strain in 2011, and the other two isolates from Burmese travelers were most similar to the Laos strain in 2008. Twenty-four amino acid substitutions were important in eight evolutionary tree branches. Comparison with DENV-1SS revealed 658 base substitutions in the local isolates, except for two mutations exclusive to YN17741, resulting in 87 synonymous mutations. Compared with the local isolates, 52 amino acid mutations occurred in the CDS of two isolates from Burmese travelers. Comparing MY21 with MY22, 17 amino acid mutations were observed, all these mutations occurred in the CDS of non-structured proteins (two in NS1, 10 in NS2, two in NS3, three in NS5). Secondary structure prediction revealed 46 changes in the potential nucleotide and protein binding sites of the CDSs in local isolates. RNA secondary structure prediction also showed base changes in the 3′UTR of local isolates, leading to two significant changes in the RNA secondary structure. To our knowledge, this study is the first complete genome analysis of isolates from the 2017 dengue outbreak that occurred at the border areas of China, Burma, and Laos.

Drivers of Dengue Intrahost Evolution

RNA viruses circulate as rapidly evolving swarms of related variants, and their evolutionary dynamics within hosts may be key to understanding virus emergence. In this issue of Cell Host & Microbe, Parameswaran et al. (2017) combine next-generation sequencing and functional studies to characterize viral populations in acute dengue infections.