Dengue-1 virus clade replacement in Thailand associated with enhanced mosquito transmission

Dengue viruses (DENV) are characterized by extensive genetic diversity and can be organized in multiple, genetically distinct lineages that arise and die out on a regular basis in regions where dengue is endemic. A fundamental question for understanding DENV evolution is the relative extent to which stochastic processes (genetic drift) and natural selection acting on fitness differences among lineages contribute to lineage diversity and turnover. Here, we used a set of recently collected and archived low-passage DENV-1 isolates from Thailand to examine the role of mosquito vector-virus interactions in DENV evolution. By comparing the ability of 23 viruses isolated on different dates between 1985 and 2009 to be transmitted by a present-day Aedes aegypti population from Thailand, we found that a major clade replacement event in the mid-1990s was associated with virus isolates exhibiting increased titers in the vector's hemocoel, which is predicted to result in a higher probability of transmission. This finding is consistent with the hypothesis that selection for enhanced transmission by mosquitoes is a possible mechanism underlying major DENV clade replacement events. There was significant variation in transmission potential among isolates within each clade, indicating that in addition to vector-driven selection, other evolutionary forces act to maintain viral genetic diversity. We conclude that occasional adaptive processes involving the mosquito vector can drive major DENV lineage replacement events.     

Nature and Extent of Genetic Diversity of Dengue Viruses Determined by 454 Pyrosequencing

Dengue virus (DENV) populations are characteristically highly diverse. Regular lineage extinction and replacement is an important dynamic DENV feature, and most DENV lineage turnover events are associated with increased incidence of disease. The role of genetic diversity in DENV lineage extinctions is not understood. We investigated the nature and extent of genetic diversity in the envelope (E) gene of DENV serotype 1 representing different lineages histories. A region of the DENV genome spanning the E gene was amplified and sequenced by Roche/454 pyrosequencing. The pyrosequencing results identified distinct sub-populations (haplotypes) for each DENV-1 E gene. A phylogenetic tree was constructed with the consensus DENV-1 E gene nucleotide sequences, and the sequences of each constructed haplotype showed that the haplotypes segregated with the Sanger consensus sequence of the population from which they were drawn. Haplotypes determined through pyrosequencing identified a recombinant DENV genome that could not be identified through Sanger sequencing. Nucleotide level sequence diversities of DENV-1 populations determined from SNP analysis were very low, estimated from 0.009–0.01. There were also no stop codon, frameshift or non-frameshift mutations observed in the E genes of any lineage. No significant correlations between the accumulation of deleterious mutations or increasing genetic diversity and lineage extinction were observed (p>0.5). Although our hypothesis that accumulation of deleterious mutations over time led to the extinction and replacement of DENV lineages was ultimately not supported by the data, our data does highlight the significant technical issues that must be resolved in the way in which population diversity is measured for DENV and other viruses. The results provide an insight into the within-population genetic structure and diversity of DENV-1 populations.     

Frequent In-Migration and Highly Focal Transmission of Dengue Viruses among Children in Kamphaeng Phet,Thailand

Revealing the patterns and determinants of the spread of dengue virus (DENV) at local scales is central to understanding the epidemiology and evolution of this major human pathogen. We performed a phylogenetic analysis of the envelope (E) genes of DENV-1, -2, -3, and -4 isolates (involving 97, 23, 5, and 74 newly collected sequences, respectively) sampled from school-based cohort and village-based cluster studies in Kamphaeng Phet, Thailand, between 2004 and 2007. With these data, we sought to describe the spatial and temporal patterns of DENV spread within a rural population where a future vaccine efficacy trial is planned. Our analysis revealed considerable genetic diversity within the study population, with multiple lineages within each serotype circulating for various lengths of time during the study period. These results suggest that DENV is frequently introduced into both semi-urban and rural areas in Kamphaeng Phet from other populations. In contrast, the persistence of viral lineages across sampling years was observed less frequently. Analysis of phylogenetic clustering indicated that DENV transmission was highly spatially and temporally focal, and that it occurred in homes rather than at school. Overall, the strength of temporal clustering suggests that seasonal bottlenecks in local DENV populations facilitate the invasion and establishment of viruses from outside of the study area, in turn reducing the extent of lineage persistence.     

Phylogeography of Recently Emerged DENV-2 in Southern Viet Nam

Revealing the dispersal of dengue viruses (DENV) in time and space is central to understanding their epidemiology. However, the processes that shape DENV transmission patterns at the scale of local populations are not well understood, particularly the impact of such factors as human population movement and urbanization. Herein, we investigated trends in the spatial dynamics of DENV-2 transmission in the highly endemic setting of southern Viet Nam. Through a phylogeographic analysis of 168 full-length DENV-2 genome sequences obtained from hospitalized dengue cases from 10 provinces in southern Viet Nam, we reveal substantial genetic diversity in both urban and rural areas, with multiple lineages identified in individual provinces within a single season, and indicative of frequent viral migration among communities. Focusing on the recently introduced Asian I genotype, we observed particularly high rates of viral exchange between adjacent geographic areas, and between Ho Chi Minh City, the primary urban center of this region, and populations across southern Viet Nam. Within Ho Chi Minh City, patterns of DENV movement appear consistent with a gravity model of virus dispersal, with viruses traveling across a gradient of population density. Overall, our analysis suggests that Ho Chi Minh City may act as a source population for the dispersal of DENV across southern Viet Nam, and provides further evidence that urban areas of Southeast Asia play a primary role in DENV transmission. However, these data also indicate that more rural areas are also capable of maintaining virus populations and hence fueling DENV evolution over multiple seasons.     

Successful propagation of flavivirus infectious cDNAs by a novel method to reduce the cryptic bacterial promoter activity of virus genomes

Reverse genetics is a powerful tool to study single-stranded RNA viruses. Despite tremendous efforts having been made to improve the methodology for constructing flavivirus cDNAs, the cause of toxicity of flavivirus cDNAs in bacteria remains unknown. Here we performed mutational analysis studies to identify Escherichia coli promoter (ECP) sequences within nucleotides (nt) 1 to 3000 of the dengue virus type 2 (DENV2) and Japanese encephalitis virus (JEV) genomes. Eight and four active ECPs were demonstrated within nt 1 to 3000 of the DENV2 and JEV genomes, respectively, using fusion constructs containing DENV2 or JEV segments and empty vector reporter gene Renilla luciferase. Full-length DENV2 and JEV cDNAs were obtained by inserting mutations reducing their ECP activity in bacteria without altering amino acid sequences. A severe cytopathic effect occurred when BHK21 cells were transfected with in vitro-transcribed RNAs from either a DENV2 cDNA clone with multiple silent mutations within the prM-E-NS1 region of dengue genome or a JEV cDNA clone with an A-to-C mutation at nt 90 of the JEV genome. The virions derived from the DENV2 or JEV cDNA clone exhibited infectivities similar to those of their parental viruses in C6/36 and BHK21 cells. A cis-acting element essential for virus replication was revealed by introducing silent mutations into the central portion (nt 160 to 243) of the core gene of DENV2 infectious cDNA or a subgenomic DENV2 replicon clone. This novel strategy of constructing DENV2 and JEV infectious clones could be applied to other flaviviruses or pathogenic RNA viruses to facilitate research in virology, viral pathogenesis, and vaccine development.     

Emergence of a New Lineage of Dengue Virus Type 2 Identified in Travelers Entering Western Australia from Indonesia, 2010-2012

Dengue virus (DENV) transmission is ubiquitous throughout the tropics. More than 70% of the current global dengue disease burden is borne by people who live in the Asia-Pacific region. We sequenced the E gene of DENV isolated from travellers entering Western Australia between 2010–2012, most of whom visited Indonesia, and identified a diverse array of DENV1-4, including multiple co-circulating viral lineages. Most viruses were closely related to lineages known to have circulated in Indonesia for some time, indicating that this geographic region serves as a major hub for dengue genetic diversity. Most notably, we identified a new lineage of DENV-2 (Cosmopolitan genotype) that emerged in Bali in 2011–2012. The spread of this lineage should clearly be monitored. Surveillance of symptomatic returned travellers provides important and timely information on circulating DENV serotypes and genotypes, and can reveal the herald wave of dengue and other emerging infectious diseases.     

Genome-wide patterns of intrahuman dengue virus diversity reveal associations with viral phylogenetic clade and interhost diversity

Analogous to observations in RNA viruses such as human immunodeficiency virus, genetic variation associated with intrahost dengue virus (DENV) populations has been postulated to influence viral fitness and disease pathogenesis. Previous attempts to investigate intrahost genetic variation in DENV characterized only a few viral genes or a limited number of full-length genomes. We developed a whole-genome amplification approach coupled with deep sequencing to capture intrahost diversity across the entire coding region of DENV-2. Using this approach, we sequenced DENV-2 genomes from the serum of 22 Nicaraguan individuals with secondary DENV infection and captured ～75% of the DENV genome in each sample (range, 40 to 98%). We identified and quantified variants using a highly sensitive and specific method and determined that the extent of diversity was considerably lower than previous estimates. Significant differences in intrahost diversity were detected between genes and also between antigenically distinct domains of the Envelope gene. Interestingly, a strong association was discerned between the extent of intrahost diversity in a few genes and viral clade identity. Additionally, the abundance of viral variants within a host, as well as the impact of viral mutations on amino acid encoding and predicted protein function, determined whether intrahost variants were observed at the interhost level in circulating Nicaraguan DENV-2 populations, strongly suggestive of purifying selection across transmission events. Our data illustrate the value of high-coverage genome-wide analysis of intrahost diversity for high-resolution mapping of the relationship between intrahost diversity and clinical, epidemiological, and virological parameters of viral infection.     

Microevolution of Dengue viruses circulating among primary school children in Kamphaeng Phet, Thailand

To determine the extent and structure of genetic variation in dengue viruses (DENV) on a restricted spatial and temporal scale, we sequenced the E (envelope) genes of DENV-1, -2, and -3 isolates collected in 2001 from children enrolled in a prospective school-based study in Kamphaeng Phet, Thailand, and diagnosed with dengue disease. Our analysis revealed substantial viral genetic variation in both time and space, with multiple viral lineages circulating within individual schools, suggesting the frequent gene flow of DENV into this microenvironment. More-detailed analyses of DENV-2 samples revealed strong clustering of viral isolates within individual schools and evidence of more-frequent viral gene flow among schools closely related in space. Conversely, we observed little evolutionary change in those viral isolates sampled over multiple time points within individual schools, indicating a low rate of mutation fixation. These results suggest that frequent viral migration into Kamphaeng Phet, coupled with population (school) subdivision, shapes the genetic diversity of DENV on a local scale, more so than in situ evolution within school catchment areas.     

Dengue Virus in Sub-tropical Northern and Central Viet Nam: Population Immunity and Climate Shape Patterns of Viral Invasion and Maintenance

Dengue virus transmission occurs in both epidemic and endemic cycles across tropical and sub-tropical regions of the world. Incidence is particularly high in much of Southeast Asia, where hyperendemic transmission plagues both urban and rural populations. However, endemicity has not been established in some areas with climates that may not support year-round viral transmission. An understanding of how dengue viruses (DENV) enter these environments and whether the viruses persist in inapparent local transmission cycles is central to understanding how dengue emerges in areas at the margins of endemic transmission. Dengue is highly endemic in tropical southern Vietnam, while increasingly large seasonal epidemics have occurred in northern Viet Nam over the last decade. We have investigated the spread of DENV-1 throughout Vietnam to determine the routes by which the virus enters northern and central regions of the country. Phylogeographic analysis of 1,765 envelope (E) gene sequences from Southeast Asia revealed frequent movement of DENV between neighboring human populations and strong local clustering of viral lineages. Long-distance migration of DENV between human population centers also occurred regularly and on short time-scales, indicating human-mediated viral invasion into northern Vietnam. Human populations in southern Vietnam were found to be the primary source of DENV circulating throughout the country, while central and northern Vietnam acted as sink populations, likely due to reduced connectedness to other populations in the case of the central regions and to the influence of temperature variability on DENV replication and vector survival and competence in the north. Finally, phylogeographic analyses suggested that viral movement follows a gravity model and indicates that population immunity and physical and economic connections between populations may play important roles in shaping patterns of DENV transmission.     

High mating frequency and variation with lineage ratio in dependent-lineage harvester ants

Explaining the evolution of multiple mating is a challenge because of the associated costs. For social insects, mating frequency may have fitness consequences due to effects on social interactions or genetic diversity within colonies. Here, we investigated the evolution of mating frequency in a social insect species with a unique genetic system that requires multiple mating. In certain populations of Pogonomyrmex harvester ants, there are two interbreeding yet genetically distinct mitochondrial lineages. Queens must mate with males of the opposite lineage to produce workers and with males of the same lineage to produce reproductive females. We expected queens of the dependent-lineage system to exhibit high mating frequencies relative to other social insects. Furthermore, we expected queens from populations of highly asymmetric lineage ratios to exhibit even higher mating frequencies, to adequately sample the population and successfully mate with males of the less common lineage. To test these predictions, we estimated the mating frequency of 16 P. barbatus queens, and compared these mating frequencies between two populations, one with relatively equal lineage ratio (60:40) and a second with a highly asymmetrical lineage ratio (96:4). Overall, observed mating frequency exceeded 10, which is high in comparison to other social insects, and our estimates of effective mating frequency were among the highest of Pogonomyrmex species. Mating frequency at the site with the asymmetrical lineage ratio was also significantly higher than the site with the more even ratio. Our results suggest that obligate multiple mating as well as lineage ratio contribute to the evolution of high mating frequency in dependent-lineage populations.     

Fitness Impact and Stability of a Transgene Conferring Resistance to Dengue-2 Virus following Introgression into a Genetically Diverse Aedes aegypti Strain

In 2006, we reported a mariner (Mos1)-transformed Aedes aegypti line, Carb77, which was highly resistant to dengue-2 virus (DENV2). Carb77 mosquitoes expressed a DENV2-specific inverted-repeat (IR) RNA in midgut epithelial cells after ingesting an infectious bloodmeal. The IR-RNA formed double-stranded DENV2-derived RNA, initiating an intracellular antiviral RNA interference (RNAi) response. However, Carb77 mosquitoes stopped expressing the IR-RNA after 17 generations in culture and lost their DENV2-refractory phenotype. In the current study, we generated new transgenic lines having the identical transgene as Carb77. One of these lines, Carb109M, has been genetically stable and refractory to DENV2 for >33 generations. Southern blot analysis identified two transgene integration sites in Carb109M. Northern blot analysis detected abundant, transient expression of the IR-RNA 24 h after a bloodmeal. Carb109M mosquitoes were refractory to different DENV2 genotypes but not to other DENV serotypes. To further test fitness and stability, we introgressed the Carb109M transgene into a genetically diverse laboratory strain (GDLS) by backcrossing for five generations and selecting individuals expressing the transgene''s EGFP marker in each generation. Comparison of transgene stability in replicate backcross 5 (BC₅) lines versus BC₁ control lines demonstrated that backcrossing dramatically increased transgene stability. We subjected six BC₅ lines to five generations of selection based on EGFP marker expression to increase the frequency of the transgene prior to final family selection. Comparison of the observed transgene frequencies in the six replicate lines relative to expectations from Fisher''s selection model demonstrated lingering fitness costs associated with either the transgene or linked deleterious genes. Although minimal fitness loss (relative to GDLS) was manifest in the final family selection stage, we were able to select homozygotes for the transgene in one family, Carb109M/GDLS.BC₅.HZ. This family has been genetically stable and DENV2 refractory for multiple generations. Carb109M/GDLS.BC₅.HZ represents an important line for testing proof-of-principle vector population replacement.     

Establishment and lineage replacement of H6 influenza viruses in domestic ducks in southern China

Domestic ducks in southern China act as an important reservoir for influenza viruses and have also facilitated the establishment of multiple H6 influenza virus lineages. To understand the continuing evolution of these established lineages, 297 H6 viruses isolated from domestic ducks during 2006 and 2007 were genetically and antigenically analyzed. Phylogenetic analyses showed that group II duck H6 viruses had replaced the previously predominant group I lineage and extended their geographic distribution from coastal to inland regions. Group II H6 virus showed that the genesis and development of multiple types of deletions in the neuraminidase (NA) stalk region could occur in the influenza viruses from domestic ducks. A gradual replacement of the N2 NA subtype with N6 was observed. Significant antigenic changes occurred within group II H6 viruses so that they became antigenically distinguishable from group I and gene pool viruses. Gene exchange between group II H6 viruses and the established H5N1, H9N2, or H6N1 virus lineages in poultry in the region was very limited. These findings suggest that domestic ducks can facilitate significant genetic and antigenic changes in viruses established in this host and highlight gaps in our knowledge of influenza virus ecology and even the evolutionary behavior of this virus family in its aquatic avian reservoirs.     

Fitness declines in Tobacco etch virus upon serial bottleneck transfers

It has been well established that populations of RNA viruses transmitted throughout serial bottlenecks suffer from significant fitness declines as a consequence of the accumulation of deleterious mutations by the onset of Muller's ratchet. Bottlenecks are unavoidably linked to different steps of the infectious cycle of most plant RNA viruses, such as vector-mediated transmissions and systemic colonization of new leaves. Here we report evidence for fitness declines by the accumulation of deleterious mutations in the potyvirus Tobacco etch virus (TEV). TEV was inoculated into the nonsystemic host Chenopodium quinoa, and local lesions were isolated and used to initiate 20 independent mutation accumulation lineages. Weekly, a random lesion from each lineage was isolated and used to inoculate the next set of plants. At each transfer, the Malthusian growth rate was estimated. After 11 consecutive transfers, all lineages suffered significant fitness losses, and one even became extinct. The average rate of fitness decline was 5% per day. The average pattern of fitness decline was consistent with antagonistic epistasis between deleterious mutations, as postulated for antiredundant genomes. Temporal fitness fluctuations were not explained by random noise but reflected more complex underlying processes related to emergence and self-organization phenomena.     

The evolutionary dynamics of influenza A virus adaptation to mammalian hosts

Few questions on infectious disease are more important than understanding how and why avian influenza A viruses successfully emerge in mammalian populations, yet little is known about the rate and nature of the virus’ genetic adaptation in new hosts. Here, we measure, for the first time, the genomic rate of adaptive evolution of swine influenza viruses (SwIV) that originated in birds. By using a curated dataset of more than 24 000 human and swine influenza gene sequences, including 41 newly characterized genomes, we reconstructed the adaptive dynamics of three major SwIV lineages (Eurasian, EA; classical swine, CS; triple reassortant, TR). We found that, following the transfer of the EA lineage from birds to swine in the late 1970s, EA virus genes have undergone substantially faster adaptive evolution than those of the CS lineage, which had circulated among swine for decades. Further, the adaptation rates of the EA lineage antigenic haemagglutinin and neuraminidase genes were unexpectedly high and similar to those observed in human influenza A. We show that the successful establishment of avian influenza viruses in swine is associated with raised adaptive evolution across the entire genome for many years after zoonosis, reflecting the contribution of multiple mutations to the coordinated optimization of viral fitness in a new environment. This dynamics is replicated independently in the polymerase genes of the TR lineage, which established in swine following separate transmission from non-swine hosts.     

Bet hedging via multiple mating: A meta‐analysis

Polyandry has been hypothesized to allow females to “bet hedge” against mating only with unsuitable mates, reducing variance in offspring fitness between members of a polyandrous lineage relative to a single‐mating one. Theoretically, this reduction in fitness variance could select for polyandrous genotypes even when polyandry carries a direct cost, especially in small populations. However, this hypothesis is controversial and difficult to test empirically. Here, I apply a novel simulation model to 49 published empirical datasets, and quantify the potential selective advantage of multiple mating via reduced offspring fitness variance. For a wide range of assumptions, including those that most favor the evolution of bet hedging, I show that any fitness gains are meager. The variance in offspring quality caused by mate identity does not appear to be high enough for bet hedging to drive the evolution of polyandry.     

The kinase inhibitor SFV785 dislocates dengue virus envelope protein from the replication complex and blocks virus assembly

Dengue virus (DENV) is the etiologic agent for dengue fever, for which there is no approved vaccine or specific anti-viral drug. As a remedy for this, we explored the use of compounds that interfere with the action of required host factors and describe here the characterization of a kinase inhibitor (SFV785), which has selective effects on NTRK1 and MAPKAPK5 kinase activity, and anti-viral activity on Hepatitis C, DENV and yellow fever viruses. SFV785 inhibited DENV propagation without inhibiting DENV RNA synthesis or translation. The compound did not cause any changes in the cellular distribution of non-structural 3, a protein critical for DENV RNA synthesis, but altered the distribution of the structural envelope protein from a reticulate network to enlarged discrete vesicles, which altered the co-localization with the DENV replication complex. Ultrastructural electron microscopy analyses of DENV-infected SFV785-treated cells showed the presence of viral particles that were distinctly different from viable enveloped virions within enlarged ER cisternae. These viral particles were devoid of the dense nucleocapsid. The secretion of the viral particles was not inhibited by SFV785, however a reduction in the amount of secreted infectious virions, DENV RNA and capsid were observed. Collectively, these observations suggest that SFV785 inhibited the recruitment and assembly of the nucleocapsid in specific ER compartments during the DENV assembly process and hence the production of infectious DENV. SFV785 and derivative compounds could be useful biochemical probes to explore the DENV lifecycle and could also represent a new class of anti-virals.     

The human immunodeficiency virus type 1 envelope confers higher rates of replicative fitness to perinatally transmitted viruses than to nontransmitted viruses

Selection of a minor viral genotype during perinatal transmission of human Immunodeficiency virus type 1 (HIV-1) has been observed, but there is a lack of information on the correlation of the restrictive transmission with biological properties of the virus, such as replicative fitness. Recombinant viruses expressing the enhanced green fluorescent protein or the Discosoma sp. red fluorescent (DsRed2) protein carrying the V1 to V5 regions of env from seven mother-infant pairs (MIPs) infected by subtype C HIV-1 were constructed, and competition assays were carried out to compare the fitness between the transmitted and nontransmitted viruses. Flow cytometry was used to quantify the frequency of infected cells, and the replicative fitness was determined based on a calculation that takes into account replication of competing viruses in a single infection versus dual infections. Transmitted viruses from five MIPs with the mothers chronically infected showed a restrictive env genotype, and all the recombinant viruses carrying the infants' Env had higher replicative fitness than those carrying the Env from the mothers. This growth fitness is lineage specific and can be observed only within the same MIP. In contrast, in two MIPs where the mothers had undergone recent acute infection, the viral Env sequences were similar between the mothers and infants and showed no further restriction in quasispecies during perinatal transmission. The recombinant viruses carrying the Env from the infants' viruses also showed replication fitness similar to those carrying the mothers' Env proteins. Our results suggest that newly transmitted viruses from chronically infected mothers have been selected to have higher replicative fitness to favor transmission, and this advantage is conferred by the V1 to V5 region of Env of the transmitted viruses. This finding has important implications for vaccine design or development of strategies to prevent HIV-1 transmission.