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

Dengue Virus Surveillance in Nepal Yields the First On-Site Whole Genome Sequences of Isolates from the 2022 Outbreak

Background

The 4 serotypes of dengue virus (DENV1-4) can each cause potentially deadly dengue disease, and are spreading globally from tropical and subtropical areas to more temperate ones. Nepal provides a microcosm of this global phenomenon, having met each of these grim benchmarks. To better understand DENV transmission dynamics and spread into new areas, we chose to study dengue in Nepal and, in so doing, to build the onsite infrastructure needed to manage future, larger studies.

Methods and Results

During the 2022 dengue season, we enrolled 384 patients presenting at a hospital in Kathmandu with dengue-like symptoms; 79% of the study participants had active or recent DENV infection (NS1 antigen and IgM). To identify circulating serotypes, we screened serum from 50 of the NS1 + participants by RT-PCR and identified DENV1, 2, and 3 - with DENV1 and 3 codominant. We also performed whole-genome sequencing of DENV, for the first time in Nepal, using our new on-site capacity. Sequencing analysis demonstrated the DENV1 and 3 genomes clustered with sequences reported from India in 2019, and the DENV2 genome clustered with a sequence reported from China in 2018.

Conclusion

These findings highlight DENV's geographic expansion from neighboring countries, identify China and India as the likely origin of the 2022 DENV cases in Nepal, and demonstrate the feasibility of building onsite capacity for more rapid genomic surveillance of circulating DENV. These ongoing efforts promise to protect populations in Nepal and beyond by informing the development and deployment of DENV drugs and vaccines in real time.

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.

Intra-Host Diversity of Dengue Virus in Mosquito Vectors

Dengue virus (DENV) is the most common arbovirus, causing a significant burden on both the economy and global healthcare systems. The virus is transmitted by Aedes species of mosquitoes as a swarm of closely related virus genomes, collectively referred to as a quasispecies. The level of genomic diversity within this quasispecies varies as DENV moves through various ecological niches within its transmission cycle. Here, the factors that influence the level of DENV quasispecies diversity during the course of infection in the mosquito vectors are reviewed.

Retrospective Genomic Characterization of a 2017 Dengue Virus Outbreak, Burkina Faso

Knowledge of contemporary genetic composition of dengue virus (DENV) in Africa is lacking. By using next-generation sequencing of samples from the 2017 DENV outbreak in Burkina Faso, we isolated 29 DENV genomes (5 serotype 1, 16 serotype 2 [DENV-2], and 8 serotype 3). Phylogenetic analysis demonstrated the endemic nature of DENV-2 in Burkina Faso. We noted discordant diagnostic results, probably related to genetic divergence between these genomes and the Trioplex PCR. Forward and reverse1 primers had a single mismatch when mapped to the DENV-2 genomes, probably explaining the insensitivity of the molecular test. Although we observed considerable homogeneity between the Dengvaxia and TetraVax-DV-TV003 vaccine strains as well as B cell epitopes compared with these genomes, we noted unique divergence. Continual surveillance of dengue virus in Africa is needed to clarify the ongoing novel evolutionary dynamics of circulating virus populations and support the development of effective diagnostic, therapeutic, and preventive countermeasures.

Evolution of Viral Pathogens Follows a Linear Order

Although lessons have been learned from previous severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS) outbreaks, the rapid evolution of the viruses means that future outbreaks of a much larger scale are possible, as shown by the current coronavirus disease 2019 (COVID-19) outbreak. Therefore, it is necessary to better understand the evolution of coronaviruses as well as viruses in general. This study reports a comparative analysis of the amino acid usage within several key viral families and genera that are prone to triggering outbreaks, including coronavirus (severe acute respiratory syndrome coronavirus 2 [SARS-CoV-2], SARS-CoV, MERS-CoV, human coronavirus-HKU1 [HCoV-HKU1], HCoV-OC43, HCoV-NL63, and HCoV-229E), influenza A (H1N1 and H3N2), flavivirus (dengue virus serotypes 1 to 4 and Zika) and ebolavirus (Zaire, Sudan, and Bundibugyo ebolavirus). Our analysis reveals that the distribution of amino acid usage in the viral genome is constrained to follow a linear order, and the distribution remains closely related to the viral species within the family or genus. This constraint can be adapted to predict viral mutations and future variants of concern. By studying previous SARS and MERS outbreaks, we have adapted this naturally occurring pattern to determine that although pangolin plays a role in the outbreak of COVID-19, it may not be the sole agent as an intermediate animal. In addition to this study, our findings contribute to the understanding of viral mutations for subsequent development of vaccines and toward developing a model to determine the source of the outbreak. IMPORTANCE This study reports a comparative analysis of amino acid usage within several key viral genera that are prone to triggering outbreaks. Interestingly, there is evidence that the amino acid usage within the viral genomes is not random but in a linear order.

Flavivirus Persistence in Wildlife Populations

A substantial number of humans are at risk for infection by vector-borne flaviviruses, resulting in considerable morbidity and mortality worldwide. These viruses also infect wildlife at a considerable rate, persistently cycling between ticks/mosquitoes and small mammals and reptiles and non-human primates and humans. Substantially increasing evidence of viral persistence in wildlife continues to be reported. In addition to in humans, viral persistence has been shown to establish in mammalian, reptile, arachnid, and mosquito systems, as well as insect cell lines. Although a considerable amount of research has centered on the potential roles of defective virus particles, autophagy and/or apoptosis-induced evasion of the immune response, and the precise mechanism of these features in flavivirus persistence have yet to be elucidated. In this review, we present findings that aid in understanding how vector-borne flavivirus persistence is established in wildlife. Research studies to be discussed include determining the critical roles universal flavivirus non-structural proteins played in flaviviral persistence, the advancement of animal models of viral persistence, and studying host factors that allow vector-borne flavivirus replication without destructive effects on infected cells. These findings underscore the viral–host relationships in wildlife animals and could be used to elucidate the underlying mechanisms responsible for the establishment of viral persistence in these animals.

Attenuated dengue viruses are genetically more diverse than their respective wild-type parents

Dengue poses a significant burden of individual health, health systems and the economy in dengue endemic regions. As such, dengue vaccine development has been an active area of research. Previous studies selected attenuated vaccine candidates based on plaque size. However, these candidates led to mixed safety outcome in clinical trials, suggesting it is insufficiently informative as an indicator of dengue virus (DENV) attenuation. In this study, we examined the genome diversity of wild-type DENVs and their attenuated derivatives developed by Mahidol University and tested in phase 1 clinical trials. We found that the attenuated DENVs, in particular the strain under clinical development by Takeda Vaccines, DENV2 PDK53, showed significantly higher genome diversity than its wild-type parent, DENV2 16681. The determinant of genomic diversity was intrinsic to the PDK53 genome as infectious clone of PDK53 showed greater genomic diversity after a single in vitro passage compared to 16681 infectious clone. Similar trends were observed with attenuated DENV1 and DENV4, both of which were shown to be attenuated clinically, but not DENV3 that was not adequately attenuated clinically. Taken together, evidence presented here suggests that genome diversity could be developed into a marker of DENV attenuation.

Alternation between taxonomically divergent hosts is not the major determinant of flavivirus evolution

Flaviviruses display diverse epidemiological and ecological features. Tick-borne and mosquito-borne flaviviruses (TBFV and MBFV, respectively) are important human pathogens that alternate replication in invertebrate vectors and vertebrate hosts. The Flavivirus genus also includes insect-specific viruses (ISFVs) and viruses with unknown invertebrate hosts. It is generally accepted that viruses that alternate between taxonomically different hosts evolve slowly and that the evolution of MBFVs and TBFVs is dominated by strong constraints, with limited episodes of positive selection. We exploited the availability of flavivirus genomes to test these hypotheses and to compare their rates and patterns of evolution. We estimated the substitution rates of CFAV and CxFV (two ISFVs) and, by taking into account the time-frame of measurement, compared them with those of other flaviviruses. Results indicated that CFAV and CxFV display relatively different substitution rates. However, these data, together with estimates for single-host members of the Flaviviridae family, indicated that MBFVs do not display relatively slower evolution. Conversely, TBFVs displayed some of lowest substitution rates among flaviviruses. Analysis of selective patterns over longer evolutionary time-frames confirmed that MBFVs evolve under strong purifying selection. Interestingly, TBFVs and ISFVs did not show extremely different levels of constraint, although TBFVs alternate among hosts, whereas ISFVs do not. Additional results showed that episodic positive selection drove the evolution of MBFVs, despite their high constraint. Positive selection was also detected on two branches of the TBFVs phylogeny that define the seabird clade. Thus, positive selection was much more common during the evolution of arthropod-borne flaviviruses than previously thought. Overall, our data indicate that flavivirus evolutionary patterns are complex and most likely determined by multiple factors, not limited to the alternation between taxonomically divergent hosts. The frequency of both positive and purifying selection, especially in MBFVs, suggests that a minority of sites in the viral polyprotein experience weak constraint and can evolve to generate new viral phenotypes and possibly promote adaptation to new hosts.

Epidemiological and genetic aspects of the largest dengue outbreak recorded in 2015 in Southeastern Brazil

Dengue virus, the etiological agent of dengue fever (DF) occurs in four genetically distinct serotypes (DENV1-4), being transmitted by female Aedes mosquitoes. DF incidence is increasing in Brazil, following vector dispersal, proliferation and DENV serotypes introduction, co-circulation and substitution. Medium- and small-sized cities in Sao Paulo State, such as Marilia (Midwest region), have been affected by huge epidemics. To understand the evolution of DENV epidemics in medium-sized cities, in this study a historical data on DENV incidence (2000-2015) in Marilia, was evaluated. Previous studies disclosed regional and specific DF outcomes associated with 2007 outbreak in that city, when co-circulating DENV1 and DENV3 presented different hematological profiles. In this study, characteristics of 2007 DENV epidemics were compared to the epidemiological, hematological and demographic outlines of the major outbreak of DENV1 in Marilia in 2015. DENV1 genetic diversity was assessed through capsid and pre-membrane junction encoding gene (CprM) sequencing. The results revealed circulation of DENV1 serotype from 2007 to 2015, with epidemics occurring every three-years until 2013 and then, increasing yearly. There were significant differences in hematological profiles of DENV1 patients between 2015 and 2007. CprM showed DENV1 genetic variability in 2015, contrasting with the unique sequence pattern in 2007. These results reinforce the regional and temporal characteristics of DENV epidemics that need local public health research to improve care for people and to limit the spread of new serotypes/genotypes to uninfected areas.

Genetic Variation in the Domain II, 3' Untranslated Region of Human and Mosquito Derived Dengue Virus Strains in Sri Lanka

Genetic variations in dengue virus (DENV) play a distinct role in epidemic emergence. The DENV 3′ UTR has become a recent interest in research. The objective of the study was to examine the genetic variation in the domain II, 3′ UTR region of human and mosquito-derived DENV. DENV-infected human sera were orally infected to laboratory reared Aedes aegypti mosquitoes. The domain II, 3′ UTR of each human- and mosquito-derived sample was amplified. The nucleotide sequence variation, phylogenetic and secondary structure analysis was carried out incorporating respective regions of so far recorded Sri Lankan and the reference genotype strains of the DENV3 and DENV1 serotypes. The human- and mosquito-derived domain II, 3′ UTR were identical in nucleotide sequences within the serotypes isolated, indicating the conserved nature of the region during host switch. The sequence analysis revealed distinct variations in study isolates compared to so far recorded Sri Lankan isolates. However, despite single nucleotide variations, the maintenance of structural integrity was evident in related strains within the serotypes in the secondary structure analysis. The phylogenetic analysis revealed distinct clade segregation of the study sequences from so far reported Sri Lankan isolates and illustrated the phylogenetic relations of the study sequences to the available global isolates of respective serotypes.

Effect of low-passage number on dengue consensus genomes and intra-host variant frequencies

Intra-host single nucleotide variants (iSNVs) have been increasingly used in genomic epidemiology to increase phylogenetic resolution and reconstruct fine-scale outbreak dynamics. These analyses are preferably done on sequence data from direct clinical samples, but in many cases due to low viral loads, there might not be enough genetic material for deep sequencing and iSNV determination. Isolation of the virus from clinical samples with low-passage number increases viral load, but few studies have investigated how dengue virus (DENV) culture isolation from a clinical sample impacts the consensus sequence and the intra-host virus population frequencies. In this study, we investigate consensus and iSNV frequency differences between DENV sequenced directly from clinical samples and their corresponding low-passage isolates. Twenty five DENV1 and DENV2 positive sera and their corresponding viral isolates (T. splendens inoculation and C6/36 passage) were obtained from a prospective cohort study in the Philippines. These were sequenced on MiSeq with minimum nucleotide depth of coverage of 500×, and iSNVs were detected using LoFreq. For both DENV1 and DENV2, we found a maximum of one consensus nucleotide difference between clinical sample and isolate. Interestingly, we found that iSNVs with frequencies ≥5 % were often preserved between the samples, and that the number of iSNV positions, and sample diversity, at this frequency cutoff did not differ significantly between the sample pairs (clinical sample and isolate) in either DENV1 or DENV2 data. Our results show that low-passage DENV isolate consensus genomes are largely representative of their direct sample parental viruses, and that low-passage isolates often mirror high frequency within-host variants from direct samples.

Full-Length Genome and Partial Viral Genes Phylogenetic and Geographical Analysis of Dengue Serotype 3 Isolates

Dengue fever is among the most common vector-borne diseases. Dengue virus (DENV), responsible for dengue fever as well as dengue hemorrhagic fever, belongs to the genus flavivirus and family Flaviviridae. Flaviviruses infect various vertebrate species and arthropods and are also responsible for diseases in birds, wild animals, and primates. DENV consists of a single-stranded, positive-sense RNA genome ~11 kb in size. Complete genome and partial gene sequences of geographically distinct DENV-3 strains were retrieved from the GenBank database. The evolutionary divergence of the 33 whole-genome and individual gene sequences of the nucleotides and amino acids of DENV-3 strains were generated with the maximum likelihood (ML) and Bayesian phylogenetic study (BEAST) methods using the MEGA 7 software. The genome size varied from 10,484 to 10,724 nucleotides among the strains with distinct geographical backgrounds belonging to Central America, South-Central Asia, and Eastern Asia. A phylogenetic analysis of the nucleotide and amino acid sequences of these DENV-3 isolates revealed extensive differences in the topologies due to PrM/M, NS1, NS2B, and NS3 genes. These results suggest substantial variation in the evolutionary pathways of the studied genes and genomes.

Principles of dengue virus evolvability derived from genotype-fitness maps in human and mosquito cells

Dengue virus (DENV) cycles between mosquito and mammalian hosts. To examine how DENV populations adapt to these different host environments, we used serial passage in human and mosquito cell lines and estimated fitness effects for all single-nucleotide variants in these populations using ultra-deep sequencing. This allowed us to determine the contributions of beneficial and deleterious mutations to the collective fitness of the population. Our analysis revealed that the continuous influx of a large burden of deleterious mutations counterbalances the effect of rare, host-specific beneficial mutations to shape the path of adaptation. Beneficial mutations preferentially map to intrinsically disordered domains in the viral proteome and cluster to defined regions in the genome. These phenotypically redundant adaptive alleles may facilitate host-specific DENV adaptation. Importantly, the evolutionary constraints described in our simple system mirror trends observed across DENV and Zika strains, indicating it recapitulates key biophysical and biological constraints shaping long-term viral evolution.

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.

Whole genome sequencing and genetic variations in several dengue virus type 1 strains from unusual dengue epidemic of 2017 in Vietnam

Background

Dengue hemorrhagic fever is an acute viral infection transmitted by mosquitoes. In the 2017, a dengue epidemic occurred in Hanoi in a short time interval and many cases were serious with associated mortality. This was the largest and unusual dengue fever outbreak in the North of Vietnam over the past 20 years. The objective of the present study was to understand the genetic characteristics of the DENV-1 strain in the 2017 epidemic and its relationship with previous viruses in Vietnam and the rest of the world.

Methods

Complete genomes of 72 DENV-1 from patients in the 2017 epidemic were sequenced using NGS. The full genome sequences were then analyzed to find out the genetic variants in the groups of 72 strains, followed by their comparison with other strains that caused disease in Vietnam previously and several other regions of the world, revealing a genetic relationship between them.

Results

The complete genome sequence of 72 DENV-1 strains comprised 10,697 nucleotides with an open reading frame coding for 3392 amino acids. The genomic analysis revealed different amino acid substitutions in all genes, especially varying at position S75 (Capsid), M125 (PrM), D54 (E), T147, V180 (NS1), G45, Y126, I154 (NS2A), A94 (NS2B), M298 (NS3), K47, V68 (NS4A), I29 (NS4B), and R166, E536, G614, T821 (NS5). The genetic analysis suggested that the viruses were most closely related to the causative virus of the dengue outbreak in Vietnam and Cambodia from 2006 to 2008. These results indicated that DENV-1 from the dengue epidemic 2017 in Northern Vietnam originated from the virus that caused the dengue outbreak during the 2007 to 2008 period in Vietnam.

Conclusion

The present study is the first of its kind to describe complete genome sequence as well as genetic variants and phylogenetic analysis of DENV-1 associated with the unusual dengue epidemic of 2017 in northern Vietnam. These results provide detailed evidence to elucidate the origin, circulation, and genetic evolution of DENV in Vietnam.

Evolution of Subgenomic RNA Shapes Dengue Virus Adaptation and Epidemiological Fitness

Changes in dengue virus (DENV) genome affect viral fitness both clinically and epidemiologically. Even in the 3′ untranslated region (3′ UTR), mutations could affect subgenomic flaviviral RNA (sfRNA) production and its affinity for host proteins, which are necessary for successful viral replication. Indeed, we recently showed that mutations in DENV2 3′ UTR of epidemic strains increased sfRNA ability to bind host proteins and reduce interferon expression. However, whether 3′ UTR differences shape the overall DENV evolution remains incompletely understood. Herein, we combined RNA phylogeny with phylogenetics to gain insights on sfRNA evolution. We found that sfRNA structures are under purifying selection and highly conserved despite sequence divergence. Only the second flaviviral nuclease-resistant RNA (fNR2) structure of DENV2 sfRNA has undergone strong positive selection. Epidemiological reports suggest that substitutions in fNR2 may drive DENV2 epidemiological fitness, possibly through sfRNA-protein interactions. Collectively, our findings indicate that 3′ UTRs are important determinants of DENV fitness in human-mosquito cycles.

•

Dengue viruses (DENVs) preserve RNA elements in their 3′ untranslated region (UTR).

•

Quantification of natural selection revealed positive selection on DENV2 sfRNA

•

Flaviviral nuclease-resistant RNAs (fNR) in the 3′ UTRs contribute to DENV speciation

•

A highly evolving fNR structure appears to increase DENV2 epidemiological fitness

Human antibody response to Zika targets type-specific quaternary structure epitopes

The recent Zika virus (ZIKV) epidemic in the Americas has revealed rare but serious manifestations of infection. ZIKV has emerged in regions endemic for dengue virus (DENV), a closely related mosquito-borne flavivirus. Cross-reactive antibodies confound studies of ZIKV epidemiology and pathogenesis. The immune responses to ZIKV may be different in people, depending on their DENV immune status. Here, we focus on the human B cell and antibody response to ZIKV as a primary flavivirus infection to define the properties of neutralizing and protective antibodies generated in the absence of preexisting immunity to DENV. The plasma antibody and memory B cell response is highly ZIKV type-specific, and ZIKV-neutralizing antibodies mainly target quaternary structure epitopes on the viral envelope. To map viral epitopes targeted by protective antibodies, we isolated 2 type-specific monoclonal antibodies (mAbs) from a ZIKV case. Both mAbs were strongly neutralizing in vitro and protective in vivo. The mAbs recognize distinct epitopes centered on domains I and II of the envelope protein. We also demonstrate that the epitopes of these mAbs define antigenic regions commonly targeted by plasma antibodies in individuals from endemic and nonendemic regions who have recovered from ZIKV infections.

Complete genome of DENV2 isolated from mosquitoes in Mexico

Dengue virus is the most prevalent arbovirus in Mexico, and although the diversity of this virus has been studied, the vast majority of sequences have been derived from viruses isolated from the human host. In this work, we aimed to sequence and to analyze DENVs derived from wild mosquitoes captured in Acapulco Guerrero, Mexico. We succeeded in determining three full genome sequences of such viruses and were able to compare them with other reported sequences from human and mosquito-derived DENVs. We found 15 nonsynonymous and 88 synonymous substitutions that were present more frequently in mosquito viruses than what would be expected by chance, although the limited number of genomes reported so far puts a constraint on the conclusions that can be derived from these analyses. Also, given the high depth of coverage attained in one of the genomes a variant analysis was carried out, finding 68 polymorphic sites in this genome. Interestingly, six of them corresponded to SNV that were detected as potentially differential between mosquitoes and humans, indicating that a that at least some positions may be maintained as polymorphic, which may facilitate host transmission.

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.

Dynamic interactions between RNA viruses and human hosts unravelled by a decade of next generation sequencing

BACKGROUND

Next generation sequencing (NGS) methods have significantly contributed to a paradigm shift in genomic research for nearly a decade now. These methods have been useful in studying the dynamic interactions between RNA viruses and human hosts.

SCOPE OF THE REVIEW

In this review, we summarise and discuss key applications of NGS in studying the host - pathogen interactions in RNA viral infections of humans with examples.

MAJOR CONCLUSIONS

Use of NGS to study globally relevant RNA viral infections have revolutionized our understanding of the within host and between host evolution of these viruses. These methods have also been useful in clinical decision-making and in guiding biomedical research on vaccine design.

GENERAL SIGNIFICANCE

NGS has been instrumental in viral genomic studies in resolving within-host viral genomic variants and the distribution of nucleotide polymorphisms along the full-length of viral genomes in a high throughput, cost effective manner. In the future, novel advances such as long read, single molecule sequencing of viral genomes and simultaneous sequencing of host and pathogens may become the standard of practice in research and clinical settings. This will also bring on new challenges in big data analysis.

Structural and Nonstructural Genes Contribute to the Genetic Diversity of RNA Viruses

With the advent of advanced sequencing technology, studies of RNA viruses have shown that genetic diversity can contribute to both attenuation and virulence and the paradigm is that this is controlled by the error-prone RNA-dependent RNA polymerase (RdRp). Since wild-type yellow fever virus (YFV) strain Asibi has genetic diversity typical of a wild-type RNA virus, while 17D virus vaccine has limited diversity, it provides a unique opportunity to investigate RNA population theory in the context of a well-characterized live attenuated vaccine. Utilizing infectious clone-derived viruses, we show that genetic diversity of RNA viruses is complex and that multiple genes, including structural genes and NS2B and NS4B genes also contribute to genetic diversity. We suggest that the replication complex as a whole, rather than only RdRp, drives genetic diversity, at least for YFV.

One paradigm to explain the complexity of viral RNA populations is that the low fidelity of the RNA-dependent RNA polymerase (RdRp) drives high mutation rates and consequently genetic diversity. Like most RNA viruses, wild-type yellow fever virus (YFV) replication is error-prone due to the lack of proofreading by the virus-encoded RdRp. However, there is evidence that replication of the live attenuated YF vaccine virus 17D, derived from wild-type strain Asibi, is less error-prone than wild-type RNA viruses. Recent studies comparing the genetic diversity of wild-type Asibi and 17D vaccine virus found that wild-type Asibi has the typical heterogeneous population of an RNA virus, while there is limited intra- and interpopulation variability of 17D vaccine virus. Utilizing chimeric and mutant infectious clone-derived viruses, we show that high and low genetic diversity profiles of wild-type Asibi virus and vaccine virus 17D, respectively, are multigenic. Introduction of either structural (pre-membrane and envelope) genes or NS2B or NS4B substitutions into the Asibi and 17D backbone resulted in altered variant population, nucleotide diversity, and mutation frequency compared to the parental viruses. Additionally, changes in genetic diversity of the chimeric and mutant viruses correlated with the phenotype of multiplication kinetics in human alveolar A549 cells. Overall, the paradigm that only the error-prone RdRp controls genetic diversity needs to be expanded to address the role of other genes in genetic diversity, and we hypothesize that it is the replication complex as a whole and not the RdRp alone that controls genetic diversity.

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.

Using bioinformatics tools for the discovery of Dengue RNA-dependent RNA polymerase inhibitors

Background

Dengue fever has rapidly manifested into a serious global health concern. The emergence of various viral serotypes has prompted the urgent need for innovative drug design techniques. Of the viral non-structural enzymes, the NS5 RNA-dependent RNA polymerase has been established as a promising target due to its lack of an enzymatic counterpart in mammalian cells and its conserved structure amongst all serotypes. The onus is now on scientists to probe further into understanding this enzyme and its mechanism of action. The field of bioinformatics has evolved greatly over recent decades, with updated drug design tools now being publically available.

Methods

In this study, bioinformatics tools were used to provide a comprehensive sequence and structural analysis of the two most prominent serotypes of Dengue RNA-dependent RNA polymerase. A list of popular flavivirus inhibitors were also chosen to dock to the active site of the enzyme. The best docked compound was then used as a template to generate a pharmacophore model that may assist in the design of target-specific Dengue virus inhibitors.

Results

Comparative sequence alignment exhibited similarity between all three domains of serotype 2 and 3.Sequence analysis revealed highly conserved regions at residues Meth530, Thr543 Asp597, Glu616, Arg659 and Pro671. Mapping of the active site demonstrated two highly conserved residues: Ser710 and Arg729. Of the active site interacting residues, Ser796 was common amongst all ten docked compounds, indicating its importance in the drug design process. Of the ten docked flavivirus inhibitors, NITD-203 showed the best binding affinity to the active site. Further pharmacophore modeling of NITD-203 depicted significant pharmacophoric elements that are necessary for stable binding to the active site.

Discussion

This study utilized publically available bioinformatics tools to provide a comprehensive framework on Dengue RNA-dependent RNA polymerase. Based on docking studies, a pharmacophore model was also designed to unveil the crucial pharmacophoric elements that are required when constructing an efficacious DENV inhibitor. We believe that this study will be a cornerstone in paving the road toward the design of target-specific inhibitors against DENV RdRp.

Genomic variations in the 3'-termini of Rice stripe virus in the rotation between vector insect and host plant

A large number of plant RNA viruses circulate between plants and insects. For RNA viruses, host alternations may impose a differential selective pressure on viral populations and induce variations in viral genomes. Here, we report the variations in the 3'-terminal regions of the multiple-segment RNA virus Rice stripe virus (RSV) that were discovered through de novo assembly of the genome using RNA sequencing data from infected host plants and vector insects. The newly assembled RSV genome contained 16- and 15-nt extensions at the 3'-termini of two genome segments compared with the published reference RSV genome. Our study demonstrated that these extensional sequences were consistently observed in two RSV isolates belonging to distinct genetic subtypes in RSV-infected rice, wheat and tobacco. Moreover, the de novo assembled genome of Southern rice black-streaked dwarf virus also contained 3'-terminal extensions in five RNA segments compared with the reference genome. Time course experiments confirmed that the 3'-terminal extensions of RSV were enriched in the vector insects, were gradually eliminated in the host plant and potentially affected viral replication. These findings indicate that variations in the 3'-termini of viral genomes may be different adaptive strategies for plant RNA viruses in insects and plants.

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

Dengue, caused by four dengue virus serotypes (DENV-1 to DENV-4), is a highly prevalent mosquito-borne viral disease in humans. Yet, selection pressures driving DENV microevolution within human hosts (intrahost) remain unknown. We employed a whole-genome segmented amplification approach coupled with deep sequencing to profile DENV-3 intrahost diversity in peripheral blood mononuclear cell (PBMC) and plasma samples from 77 dengue patients. DENV-3 intrahost diversity appears to be driven by immune pressures as well as replicative success in PBMCs and potentially other replication sites. Hotspots for intrahost variation were detected in 59%-78% of patients in the viral Envelope and pre-Membrane/Membrane proteins, which together form the virion surface. Dominant variants at the hotspots arose via convergent microevolution, appear to be immune-escape variants, and were evolutionarily constrained at the macro level due to viral replication defects. Dengue is thus an example of an acute infection in which selection pressures within infected individuals drive rapid intrahost virus microevolution.

Understanding dengue virus evolution to support epidemic surveillance and counter-measure development

Dengue virus (DENV) causes a profound burden of morbidity and mortality, and its global burden is rising due to the co-circulation of four divergent DENV serotypes in the ecological context of globalization, travel, climate change, urbanization, and expansion of the geographic range of the Ae.aegypti and Ae.albopictus vectors. Understanding DENV evolution offers valuable opportunities to enhance surveillance and response to DENV epidemics via advances in RNA virus sequencing, bioinformatics, phylogenetic and other computational biology methods. Here we provide a scoping overview of the evolution and molecular epidemiology of DENV and the range of ways that evolutionary analyses can be applied as a public health tool against this arboviral pathogen.

Identification of infection by Chikungunya, Zika, and Dengue in an area of the Peruvian coast. Molecular diagnosis and clinical characteristics

Objective

To assess the presence of Dengue, Chikungunya, and Zika in serum samples of patients with acute febrile illness in Piura, Peru and describe the most common clinical features.

Results

Dengue was the most common arbovirus detected in 170/496 (34.3%), followed by Zika in 39/496 (7.9%) and Chikungunya in 23/496 (4.6%). Among the 170 samples positive for Dengue, serotype 2 was the most predominant type present in 97/170 (57.1%) of samples, followed by the serotype 3 in 9/170 (5.3%). Headaches, muscle pain, and joint pain were the most common symptoms associated with fever in patients with Dengue and Zika. No symptoms predominance was observed in patients with Chikungunya.Dengue is considered the most frequent arbovirus in Peru and the number of cases has increased dramatically in the last 5 years. However, it is not the only arbovirus that circulates along the northern coast of Peru. It has also been determined the presence of Zika and Chikungunya in our population, which may suggest the circulation of other arboviruses that have not been detected.

Mosquitoes Transmit Unique West Nile Virus Populations during Each Feeding Episode

Arthropod-borne viruses (arboviruses), such as Zika virus, chikungunya virus, and West Nile virus (WNV), pose continuous threats to emerge and cause large epidemics. Often, these events are associated with novel virus variants optimized for local transmission that first arise as minorities within a host. Thus, the conditions that regulate the frequency of intrahost variants are important determinants of emergence. Here, we describe the dynamics of WNV genetic diversity during its transmission cycle. By temporally sampling saliva from individual mosquitoes, we demonstrate that virus populations expectorated by mosquitoes are highly diverse and unique to each feeding episode. After transmission to birds, however, most genetic diversity is removed by strong purifying selection. Further, transmission of potentially mosquito-adaptive WNV variants is strongly influenced by genetic drift in mosquitoes. These results highlight the complex evolutionary forces a novel virus variant must overcome to alter infection phenotypes at the population level.

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

Arthropod-borne viruses are among the most genetically constrained RNA viruses, yet they have a remarkable propensity to adapt and emerge. We studied wild birds and mosquitoes naturally infected with West Nile virus (WNV) in a 'hot spot' of virus transmission in Chicago, IL, USA. We generated full coding WNV genome sequences from spatiotemporally matched bird and mosquito samples using high-throughput sequencing, allowing a molecular evolutionary assessment with deep coverage. Mean FST among samples was 0.66 (±0.02 SE) and was bimodal, with mean nucleotide diversity being higher between samples (interhost πN = 0.001; πS = 0.024) than within them (intrahost πN < 0.0001; πS < 0.001). Eight genomic sites with FST > 1.01 (in the PrM, NS2a, NS3, NS4b, and 5'-noncoding genomic regions) showed bird versus mosquito variant frequency differences of >30 per cent and/or polymorphisms fixed in ≥5 host or vector individuals, suggesting host tropism for these variants. However, phylogenetic analyses demonstrated a lack of grouping by bird or mosquito, most inter-sample differences were synonymous (mean interhost πN/πS = 0.04), and there was no significant difference between hosts and vectors in either their nucleotide diversities or levels of purifying selection (mean intrahost πN/πS = 0.28 in birds and πN/πS = 0.21 in mosquitoes). This finding contrasts with the 'trade-off' and 'selective sieve' hypotheses that have been proposed and tested in the laboratory, which predict strong host versus vector effects on WNV genetic variation, with heightened selective constraint in birds alternating with heightened viral diversity in mosquitoes. Overall, our data show WNV to be highly selectively constrained within and between both hosts and vectors but still able to vary at a limited number of sites across the genome. Such site-specific plasticity in the face of overall selective constraint may offer a mechanism whereby highly constrained viruses such as WNV and its relatives can still adapt and emerge.

Phylogenetic analyses of DENV-3 isolated from field-caught mosquitoes in Thailand

Dengue virus serotype 3 (DENV-3) can cause all forms of dengue diseases and is a predominant serotype in many countries. This serotype is classified into five genotypes: I-V. Genotypes I-III have widely spread throughout the world, whereas genotypes IV and V are rare. Despite the impact on the spread of dengue diseases, only a few studies have reported the characteristics of DENV present in mosquito vectors. Hence, this study aimed to identify DENV-3 genotypes and reveal genetic variation of this virus presented in field-caught mosquitoes collected from endemic areas in Thailand during 2011-2015. First, we examined the effectiveness of the E gene sequence on DENV-3 genotyping, with results supporting the use of this gene for genotype identification. Then, we sequenced this gene in ten DENV-3 strains isolated from mosquitoes. The results showed that eight and two samples were genotypes III and V, respectively, and that they are closely related to DENV-3 isolated from Southeast and East Asian samples. The translated E gene sequences showed 25 unique amino acid (AA) residues located at 23 positions. Eight out of 25 residues have different chemical properties compared to the conserved AAs that are distributed across the three domains functioning in virus-host interaction. Hence, our study reports the first DENV-3 genotype V in Thailand, with these viruses potentially influencing both the disease severity and epidemic potential of DENV-3.

Full-genome dengue virus sequencing in mosquito saliva shows lack of convergent positive selection during transmission by Aedes aegypti

Like other pathogens with high mutation and replication rates, within-host dengue virus (DENV) populations evolve during infection of their main mosquito vector, Aedes aegypti. Within-host DENV evolution during transmission provides opportunities for adaptation and emergence of novel virus variants. Recent studies of DENV genetic diversity failed to detect convergent evolution of adaptive mutations in mosquito tissues such as midgut and salivary glands, suggesting that convergent positive selection is not a major driver of within-host DENV evolution in the vector. However, it is unknown whether this conclusion extends to the transmitted viral subpopulation because it is technically difficult to sequence DENV genomes in mosquito saliva. Here, we achieved DENV full-genome sequencing by pooling saliva samples collected non-sacrificially from 49 to 163 individual Ae. aegypti mosquitoes previously infected with one of two DENV-1 genotypes. We compared the transmitted viral subpopulations found in the pooled saliva samples collected in time series with the input viral population present in the infectious blood meal. In all pooled saliva samples examined, the full-genome consensus sequence of the input viral population was unchanged. Although the pooling strategy prevents analysis of individual saliva samples, our results demonstrate the lack of strong convergent positive selection during a single round of DENV transmission by Ae. aegypti. This finding reinforces the idea that genetic drift and purifying selection are the dominant evolutionary forces shaping within-host DENV genetic diversity during transmission by mosquitoes.

Complete coding sequence of dengue virus serotype 4 isolated from field-caught mosquitoes in Thailand

This report is the first to characterise the complete coding sequence of a dengue virus serotype 4 (DENV-4) genotype I that was isolated from field-caught mosquitoes from an endemic area in Thailand in June 2013. The sequence was assembled from high-throughput sequencing reads generated by Illumina HiSeq. Three out of four observed intra-sample variants caused an amino acid variation in C, NS2B, and NS5 genes. The C4279T variant located in the NS2B gene can indirectly affect the proteolytic activity of the NS3 protein. The sequence provided in this study might be useful for the epidemiological study of DENV-4.

Dengue virus genomic variation associated with mosquito adaptation defines the pattern of viral non-coding RNAs and fitness in human cells

The Flavivirus genus includes a large number of medically relevant pathogens that cycle between humans and arthropods. This host alternation imposes a selective pressure on the viral population. Here, we found that dengue virus, the most important viral human pathogen transmitted by insects, evolved a mechanism to differentially regulate the production of viral non-coding RNAs in mosquitos and humans, with a significant impact on viral fitness in each host. Flavivirus infections accumulate non-coding RNAs derived from the viral 3’UTRs (known as sfRNAs), relevant in viral pathogenesis and immune evasion. We found that dengue virus host adaptation leads to the accumulation of different species of sfRNAs in vertebrate and invertebrate cells. This process does not depend on differences in the host machinery; but it was found to be dependent on the selection of specific mutations in the viral 3’UTR. Dissecting the viral population and studying phenotypes of cloned variants, the molecular determinants for the switch in the sfRNA pattern during host change were mapped to a single RNA structure. Point mutations selected in mosquito cells were sufficient to change the pattern of sfRNAs, induce higher type I interferon responses and reduce viral fitness in human cells, explaining the rapid clearance of certain viral variants after host change. In addition, using epidemic and pre-epidemic Zika viruses, similar patterns of sfRNAs were observed in mosquito and human infected cells, but they were different from those observed during dengue virus infections, indicating that distinct selective pressures act on the 3’UTR of these closely related viruses. In summary, we present a novel mechanism by which dengue virus evolved an RNA structure that is under strong selective pressure in the two hosts, as regulator of non-coding RNA accumulation and viral fitness. This work provides new ideas about the impact of host adaptation on the variability and evolution of flavivirus 3’UTRs with possible implications in virulence and viral transmission.

Flaviviruses constitute the most important and diverse group of arthropod-transmitted viruses, including relevant human pathogens such as dengue, Zika, yellow fever, and West Nile viruses. The natural alternation of these viruses between vertebrate and invertebrate hosts imposes a selective pressure on the viral population, with potential epidemiological implications. However, the selective forces and mechanisms that act on the viral RNA during host adaptation are largely unknown. Using dengue virus, we found that viral replication in mosquito or human cells leads to the accumulation of different patterns of viral non-coding RNAs that differently regulate viral fitness in each host. Carrying out this process requires changes of the structure of the viral 3’UTR that is under strong selective pressure in the two hosts. Our findings provide the first link between flavivirus host adaptation, fitness, and the production of viral non-coding RNAs, and support a model in which opposite selective pressures in the two hosts drive flavivirus 3’UTR evolution.

Extended Evaluation of Virological, Immunological and Pharmacokinetic Endpoints of CELADEN: A Randomized, Placebo-Controlled Trial of Celgosivir in Dengue Fever Patients

CELADEN was a randomized placebo-controlled trial of 50 patients with confirmed dengue fever to evaluate the efficacy and safety of celgosivir (A study registered at ClinicalTrials.gov, number NCT01619969). Celgosivir was given as a 400 mg loading dose and 200 mg bid (twice a day) over 5 days. Replication competent virus was measured by plaque assay and compared to reverse transcription quantitative PCR (qPCR) of viral RNA. Pharmacokinetics (PK) correlations with viremia, immunological profiling, next generation sequence (NGS) analysis and hematological data were evaluated as exploratory endpoints here to identify possible signals of pharmacological activity. Viremia by plaque assay strongly correlated with qPCR during the first four days. Immunological profiling demonstrated a qualitative shift in T helper cell profile during the course of infection. NGS analysis did not reveal any prominent signature that could be associated with drug treatment; however the phylogenetic spread of patients’ isolates underlines the importance of strain variability that may potentially confound interpretation of dengue drug trials conducted during different outbreaks and in different countries. Celgosivir rapidly converted to castanospermine (Cast) with mean peak and trough concentrations of 5727 ng/mL (30.2 μM) and 430 ng/mL (2.3 μM), respectively and cleared with a half-life of 2.5 (± 0.6) hr. Mean viral log reduction between day 2 and 4 (VLR2-4) was significantly greater in secondary dengue than primary dengue (p = 0.002). VLR2-4 did not correlate with drug AUC but showed a trend of greater response with increasing Cmin. PK modeling identified dosing regimens predicted to achieve 2.4 to 4.5 times higher Cmin. than in the CELADEN trial for only 13% to 33% increase in overall dose. A small, non-statistical trend towards better outcome on platelet nadir and difference between maximum and minimum hematocrit was observed in celgosivir-treated patients with secondary dengue infection. Optimization of the dosing regimen and patient stratification may enhance the ability of a clinical trial to demonstrate celgosivir activity in treating dengue fever based on hematological endpoints. A new clinical trial with a revised dosing regimen is slated to start in 2016 (NCT02569827). Furthermore celgosivir’s potential value for treatment of other flaviruses such as Zika virus should be investigated urgently.

Trial Registration: ClinicalTrials.gov NCT01619969

Dengue virus is currently threatening 40% of the world’s population. An approximately 60% efficacious vaccine has been registered for use in Mexico, Brazil, the Philippines, Paraguay and El Salvador, but there are no approved antiviral treatments available. We have shown that celgosivir, an endoplasmic reticulum alpha glucosidase inhibitor, has submicromolar activity against all 4 serotypes of dengue virus (DENV) and also efficacious in mouse model of infection. The strong preclinical pharmacology motivated the conduct of an investigator-initiated, Phase 1b randomized, double-blind, placebo-controlled trial of celgosivir in 50 adult dengue patients. Although the trial did not meet the primary endpoints of lowering viremia or fever, the safety profile of the drug prompted extended hematological, pharmacokinetic, immunological and viral sequence profiling. Here we report several non-significant trends of pharmacological effect of celgosivir on platelet count, hematocrit, and NS1 clearance in secondary dengue patients. In addition, pharmacokinetic modeling identified an alternate dosing regimen that is predicted to achieve a 4.5-fold increase in minimum drug concentrations during treatment (C_min) with only a modest increase in overall dose. A new Phase 2a clinical trial with an optimized dosing regimen of celgosivir (ClinicalTrials.gov number NCT02569827) is scheduled to start in the latter part of 2016.

Vectors as Epidemiological Sentinels: Patterns of Within-Tick Borrelia burgdorferi Diversity

Hosts including humans, other vertebrates, and arthropods, are frequently infected with heterogeneous populations of pathogens. Within-host pathogen diversity has major implications for human health, epidemiology, and pathogen evolution. However, pathogen diversity within-hosts is difficult to characterize and little is known about the levels and sources of within-host diversity maintained in natural populations of disease vectors. Here, we examine genomic variation of the Lyme disease bacteria, Borrelia burgdorferi (Bb), in 98 individual field-collected tick vectors as a model for study of within-host processes. Deep population sequencing reveals extensive and previously undocumented levels of Bb variation: the majority (~70%) of ticks harbor mixed strain infections, which we define as levels Bb diversity pre-existing in a diverse inoculum. Within-tick diversity is thus a sample of the variation present within vertebrate hosts. Within individual ticks, we detect signatures of positive selection. Genes most commonly under positive selection across ticks include those involved in dissemination in vertebrate hosts and evasion of the vertebrate immune complement. By focusing on tick-borne Bb, we show that vectors can serve as epidemiological and evolutionary sentinels: within-vector pathogen diversity can be a useful and unbiased way to survey circulating pathogen diversity and identify evolutionary processes occurring in natural transmission cycles.

Lyme disease, caused by a bacteria carried by deer ticks, is the most common vector-borne disease in North America and over 30,000 cases are reported each year in the United States. Ticks may be infected with multiple strains of the Lyme disease bacteria, which differ in transmissibility and the harm they pose to humans. In this study, we collected 98 infected deer ticks from across the United States and southern Canada. We used genetic techniques to investigate the diversity of the Lyme disease bacteria infecting each individual tick. We find that 70% of ticks are infected with multiple strains of the Lyme disease bacteria, indicating that humans may be exposed to and infected with multiple bacterial strains from a single tick bite. We also find evidence that the Lyme disease bacteria is evolving in response to the immune defenses of its natural hosts (including rodents and birds). Our study shows that individual ticks and other disease vectors can be studied as epidemiological sentinels, which reveal the extensive diversity of pathogens circulating in natural disease cycles and how they are evolving.

Genetic Drift, Purifying Selection and Vector Genotype Shape Dengue Virus Intra-host Genetic Diversity in Mosquitoes

Due to their error-prone replication, RNA viruses typically exist as a diverse population of closely related genomes, which is considered critical for their fitness and adaptive potential. Intra-host demographic fluctuations that stochastically reduce the effective size of viral populations are a challenge to maintaining genetic diversity during systemic host infection. Arthropod-borne viruses (arboviruses) traverse several anatomical barriers during infection of their arthropod vectors that are believed to impose population bottlenecks. These anatomical barriers have been associated with both maintenance of arboviral genetic diversity and alteration of the variant repertoire. Whether these patterns result from stochastic sampling (genetic drift) rather than natural selection, and/or from the influence of vector genetic heterogeneity has not been elucidated. Here, we used deep sequencing of full-length viral genomes to monitor the intra-host evolution of a wild-type dengue virus isolate during infection of several mosquito genetic backgrounds. We estimated a bottleneck size ranging from 5 to 42 founding viral genomes at initial midgut infection, irrespective of mosquito genotype, resulting in stochastic reshuffling of the variant repertoire. The observed level of genetic diversity increased following initial midgut infection but significantly differed between mosquito genetic backgrounds despite a similar initial bottleneck size. Natural selection was predominantly negative (purifying) during viral population expansion. Taken together, our results indicate that dengue virus intra-host genetic diversity in the mosquito vector is shaped by genetic drift and purifying selection, and point to a novel role for vector genetic factors in the genetic breadth of virus populations during infection. Identifying the evolutionary forces acting on arboviral populations within their arthropod vector provides novel insights into arbovirus evolution.

During infection of their arthropod vectors, arthropod-borne viruses (arboviruses) such as dengue viruses traverse several anatomical barriers that are believed to cause dramatic reductions in population size. Such population bottlenecks challenge the maintenance of viral genetic diversity, which is considered critical for fitness and adaptability of arboviruses. Anatomical barriers in the vector were previously associated with both maintenance of arboviral genetic diversity and alteration of the variant repertoire. However, the relative role of random processes and natural selection, and the influence of vector genetic heterogeneity have not been elucidated. In this study, we used high-throughput sequencing to monitor dengue virus genetic diversity during infection of several genetic backgrounds of their mosquito vector. Our results show that initial infection of the vector is randomly founded by only a few tens of individual virus genomes. The overall level of viral genetic diversity generated during infection was predominantly under purifying selection but differed significantly between mosquito genetic backgrounds. Thus, in addition to random evolutionary forces and the purging of deleterious mutations that shape dengue virus genetic diversity during vector infection, our results also point to a novel role for vector genetic factors in the genetic breadth of virus populations.

Increasing Clinical Severity during a Dengue Virus Type 3 Cuban Epidemic: Deep Sequencing of Evolving Viral Populations

During the dengue virus type 3 (DENV-3) epidemic that occurred in Havana in 2001 to 2002, severe disease was associated with the infection sequence DENV-1 followed by DENV-3 (DENV-1/DENV-3), while the sequence DENV-2/DENV-3 was associated with mild/asymptomatic infections. To determine the role of the virus in the increasing severity demonstrated during the epidemic, serum samples collected at different time points were studied. A total of 22 full-length sequences were obtained using a deep-sequencing approach. Bayesian phylogenetic analysis of consensus sequences revealed that two DENV-3 lineages were circulating in Havana at that time, both grouped within genotype III. The predominant lineage is closely related to Peruvian and Ecuadorian strains, while the minor lineage is related to Venezuelan strains. According to consensus sequences, relatively few nonsynonymous mutations were observed; only one was fixed during the epidemic at position 4380 in the NS2B gene. Intrahost genetic analysis indicated that a significant minor population was selected and became predominant toward the end of the epidemic. In conclusion, greater variability was detected during the epidemic's progression in terms of significant minority variants, particularly in the nonstructural genes. An increasing trend of genetic diversity toward the end of the epidemic was observed only for synonymous variant allele rates, with higher variability in secondary cases. Remarkably, significant intrahost genetic variation was demonstrated within the same patient during the course of secondary infection with DENV-1/DENV-3, including changes in the structural proteins premembrane (PrM) and envelope (E). Therefore, the dynamic of evolving viral populations in the context of heterotypic antibodies could be related to the increasing clinical severity observed during the epidemic.

IMPORTANCE

Based on the evidence that DENV fitness is context dependent, our research has focused on the study of viral factors associated with intraepidemic increasing severity in a unique epidemiological setting. Here, we investigated the intrahost genetic diversity in acute human samples collected at different time points during the DENV-3 epidemic that occurred in Cuba in 2001 to 2002 using a deep-sequencing approach. We concluded that greater variability in significant minor populations occurred as the epidemic progressed, particularly in the nonstructural genes, with higher variability observed in secondary infection cases. Remarkably, for the first time significant intrahost genetic variation was demonstrated within the same patient during the course of secondary infection with DENV-1/DENV-3, including changes in structural proteins. These findings indicate that high-resolution approaches are needed to unravel molecular mechanisms involved in dengue pathogenesis.

Genomic approaches for understanding dengue: insights from the virus, vector, and host

The incidence and geographic range of dengue have increased dramatically in recent decades. Climate change, rapid urbanization and increased global travel have facilitated the spread of both efficient mosquito vectors and the four dengue virus serotypes between population centers. At the same time, significant advances in genomics approaches have provided insights into host–pathogen interactions, immunogenetics, and viral evolution in both humans and mosquitoes. Here, we review these advances and the innovative treatment and control strategies that they are inspiring.

True Blood: dengue virus evolution

This month's Genome Watch highlights insights into different selection pressures in dengue virus evolution.