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.     

Clade replacements in dengue virus serotypes 1 and 3 are associated with changing serotype prevalence

The evolution of dengue virus (DENV) is characterized by phylogenetic trees that have a strong temporal structure punctuated by dramatic changes in clade frequency. To determine the cause of these large-scale phylogenetic patterns, we examined the evolutionary history of DENV serotype 1 (DENV-1) and DENV-3 in Thailand, where gene sequence and epidemiological data are relatively abundant over a 30-year period. We found evidence for the turnover of viral clades in both serotypes, most notably in DENV-1, where a major clade replacement event took place in genotype I during the mid-1990s. Further, when this clade replacement event was placed in the context of changes in serotype prevalence in Thailand, a striking pattern emerged; an increase in DENV-1 clade diversity was associated with an increase in the abundance of this serotype and a concomitant decrease in DENV-4 prevalence, while clade replacement was associated with a decline in DENV-1 prevalence and a rise of DENV-4. We postulate that intraserotypic genetic diversification proceeds at times of relative serotype abundance and that replacement events can result from differential susceptibility to cross-reactive immune responses.     

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.     

Genetic specificity and potential for local adaptation between dengue viruses and mosquito vectors

Background  

Several observations support the hypothesis that vector-driven selection plays an important role in shaping dengue virus (DENV) genetic diversity. Clustering of DENV genetic diversity at a particular location may reflect underlying genetic structure of vector populations, which combined with specific vector genotype × virus genotype (G × G) interactions may promote adaptation of viral lineages to local mosquito vector genotypes. Although spatial structure of vector polymorphism at neutral genetic loci is well-documented, existence of G × G interactions between mosquito and virus genotypes has not been formally demonstrated in natural populations. Here we measure G × G interactions in a system representative of a natural situation in Thailand by challenging three isofemale families from field-derived Aedes aegypti with three contemporaneous low-passage isolates of DENV-1.     

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.     

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.     

Distribution of Fitness in Populations of Dengue Viruses

Genetically diverse RNA viruses like dengue viruses (DENVs) segregate into multiple, genetically distinct, lineages that temporally arise and disappear on a regular basis. Lineage turnover may occur through multiple processes such as, stochastic or due to variations in fitness. To determine the variation of fitness, we measured the distribution of fitness within DENV populations and correlated it with lineage extinction and replacement. The fitness of most members within a population proved lower than the aggregate fitness of populations from which they were drawn, but lineage replacement events were not associated with changes in the distribution of fitness. These data provide insights into variations in fitness of DENV populations, extending our understanding of the complexity between members of individual populations.     

Tracking Dengue Virus Intra-host Genetic Diversity during Human-to-Mosquito Transmission

Dengue virus (DENV) infection of an individual human or mosquito host produces a dynamic population of closely-related sequences. This intra-host genetic diversity is thought to offer an advantage for arboviruses to adapt as they cycle between two very different host species, but it remains poorly characterized. To track changes in viral intra-host genetic diversity during horizontal transmission, we infected Aedes aegypti mosquitoes by allowing them to feed on DENV2-infected patients. We then performed whole-genome deep-sequencing of human- and matched mosquito-derived DENV samples on the Illumina platform and used a sensitive variant-caller to detect single nucleotide variants (SNVs) within each sample. >90% of SNVs were lost upon transition from human to mosquito, as well as from mosquito abdomen to salivary glands. Levels of viral diversity were maintained, however, by the regeneration of new SNVs at each stage of transmission. We further show that SNVs maintained across transmission stages were transmitted as a unit of two at maximum, suggesting the presence of numerous variant genomes carrying only one or two SNVs each. We also present evidence for differences in selection pressures between human and mosquito hosts, particularly on the structural and NS1 genes. This analysis provides insights into how population drops during transmission shape RNA virus genetic diversity, has direct implications for virus evolution, and illustrates the value of high-coverage, whole-genome next-generation sequencing for understanding viral intra-host genetic diversity.     

Evolution of endemic and sylvatic lineages of dengue virus

Recent disease outbreaks have raised awareness of tropical pathogens, especially mosquito-borne viruses. Dengue virus (DENV) is a widely studied mammalian pathogen transmitted by various species of mosquito in the genus Aedes, especially Aedes aegypti and Aedes albopictus. The prevailing view of the research community is that endemic viral lineages that cause epidemics of DENV in humans have emerged over time from sylvatic viral lineages, which persist in wild, non-human primates. These notions have been examined by researchers through phylogenetic analyses of the envelope gene (E) from the four serotypes of DENV (serotypes DENV-1 to DENV-4). In these previous reports, researchers used visual inspection of a phylogeny in order to assert that sylvatic lineages lead to endemic clades. In making this assertion, these researchers also reasserted the model of periodic sylvatic to endemic disease outbreaks. Since that study, there has been a significant increase in data both in terms of metadata (e.g., place and host of isolation) and genetic sequences of DENV. Here, we re-examine the model of sylvatic to endemic shifts in viral lineages through a phylogenetic tree search and character evolution study of metadata on the tree. We built a dataset of nucleotide sequences for 188 isolates of DENV that have metadata on sylvatic or endemic sampling along with three orthologous sequences from West Nile virus as the outgroup for the phylogenetic analysis. In contrast to previous research, we find that there are several shifts from endemic to sylvatic lineages as well as sylvatic to endemic lineages, indicating a much more dynamic model of evolution. We propose that a model that allows oscillation between sylvatic and endemic hosts better captures the dynamics of DENV transmission.     

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.     

Identification of the progenitors of Indonesian and Vietnamese avian influenza A (H5N1) viruses from southern China

The transmission of highly pathogenic avian influenza H5N1 virus to Southeast Asian countries triggered the first major outbreak and transmission wave in late 2003, accelerating the pandemic threat to the world. Due to the lack of influenza surveillance prior to these outbreaks, the genetic diversity and the transmission pathways of H5N1 viruses from this period remain undefined. To determine the possible source of the wave 1 H5N1 viruses, we recently conducted further sequencing and analysis of samples collected in live-poultry markets from Guangdong, Hunan, and Yunnan in southern China from 2001 to 2004. Phylogenetic analysis of the hemagglutinin and neuraminidase genes of 73 H5N1 isolates from this period revealed a greater genetic diversity in southern China than previously reported. Moreover, results show that eight viruses isolated from Yunnan in 2002 and 2003 were most closely related to the clade 1 virus sublineage from Vietnam, Thailand, and Malaysia, while two viruses from Hunan in 2002 and 2003 were most closely related to viruses from Indonesia (clade 2.1). Further phylogenetic analyses of the six internal genes showed that all 10 of those viruses maintained similar phylogenetic relationships as the surface genes. The 10 progenitor viruses were genotype Z and shared high similarity (>/=99%) with their corresponding descendant viruses in most gene segments. These results suggest a direct transmission link for H5N1 viruses between Yunnan and Vietnam and also between Hunan and Indonesia during 2002 and 2003. Poultry trade may be responsible for virus introduction to Vietnam, while the transmission route from Hunan to Indonesia remains unclear.     

Using RNA interference to develop dengue virus resistance in genetically modified Aedes aegypti

Diseases caused by arthropod-borne viruses are significant public health problems, and novel methods are needed to control pathogen transmission. We hypothesize that genetic manipulation of Aedes aegypti mosquitoes can profoundly and permanently reduce vector competence and subsequent transmission of dengue viruses (DENV) to human hosts. We have identified RNA interference (RNAi) as a potential anti-viral, intracellular pathway in the vector that can be triggered by expression of virus-specific, double stranded RNAs (dsRNAs) to reduce vector competence to DENV. We identified DENV-derived RNA segments using recombinant Sindbis viruses to trigger RNAi, that when expressed in mosquitoes ablate homologous DENV replication and transmission. We also demonstrated that heritable expression of DENV-derived dsRNA in cultured mosquito cells can silence virus replication. We now have developed a number of transgenic mosquito lines that transcribe the effector dsRNA from constitutive promoters such as immediate early 1 (baculovirus) and polyubiquitin (Drosophila melanogaster). We have detected DENV-specific small interfering RNAs, the hallmark of RNAi, in at least one of these lines. Surprisingly, none of these lines expressed dsRNA in relevant tissues (e.g., midguts) that will ultimately affect transmission. A major challenge now is to express the effector dsRNA from tissue-specific promoters to allow RNAi to silence virus replication at critical sites in the vector such as midguts and salivary glands. If successful, this strategy has the advantage of harnessing a naturally occurring vector response to block DENV infection in a mosquito vector and profoundly affect virus transmission.     

Genetic characterization of rabies field isolates from Thailand

We sequenced 512 nucleotides in two variable regions of the N gene of 23 rabies isolates from the northeastern part of Thailand by direct sequencing of PCR-amplified products. The sequencing data revealed two new lineages in these rabies isolates. Based on the results of this study together with the findings of our earlier study, the rabies isolates in Thailand were divided into two genogroups, designated as T1 and T2, which were predominantly localized in the northern and northeastern areas, respectively. Each of these two genogroups consisted of four lineages. There was a correlation between these eight lineages and the geographical origins of the isolates. Two lineages belonging to the T2 genogroup from the northeastern area of Thailand were newly identified in this study. The isolates in both genogroups were also prevalent in the central area of Thailand. Each lineage in the T1 and T2 genogroups was found independently in dogs in the upper and lower southern areas of Thailand, respectively. These genetic data and the historical background suggest that rabies viruses belonging to the T2 genogroups were prevalent many years ago in the central and northeastern areas of Thailand and were later transferred to the lower southern parts of Thailand.     

Evolutionary potential of the extrinsic incubation period of dengue virus in Aedes aegypti

Dengue fever is the most common arboviral disease worldwide. It is caused by dengue viruses (DENV) and the mosquito Aedes aegypti is its primary vector. One of the most powerful determinants of a mosquito's ability to transmit DENV is the length of the extrinsic incubation period (EIP), the time it takes for a virus to be transmitted by a mosquito after consuming an infected blood meal. Here, we repeatedly measured DENV load in the saliva of individual mosquitoes over their lifetime and used this in combination with a breeding design to determine the extent to which EIP might respond to the evolutionary forces of drift and selection. We demonstrated that genetic variation among mosquitoes contributes significantly to transmission potential and length of EIP. We reveal that shorter EIP is genetically correlated with reduced mosquito lifespan, highlighting negative life‐history consequences for virus‐infected mosquitoes. This work highlights the capacity for local genetic variation in mosquito populations to evolve and to dramatically affect the nature of human outbreaks. It also provides the impetus for isolating mosquito genes that determine EIP. More broadly, our dual experimental approach offers new opportunities for studying the evolutionary potential of transmission traits in other vector/pathogen systems.     

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.     

Changes in population dynamics during long-term evolution of sabin type 1 poliovirus in an immunodeficient patient

The evolution of the Sabin strain of type 1 poliovirus in a hypogammaglobulinemia patient for a period of 649 days is described. Twelve poliovirus isolates from sequential stool samples encompassing days 21 to 649 after vaccination with Sabin 1 were characterized in terms of their antigenic properties, virulence in transgenic mice, sensitivity for growth at high temperatures, and differences in nucleotide sequence from the Sabin 1 strain. Poliovirus isolates from the immunodeficient patient evolved gradually toward non-temperature-sensitive and neurovirulent phenotypes, accumulating mutations at key nucleotide positions that correlated with the observed reversion to biological properties typical of wild polioviruses. Analysis of plaque-purified viruses from stool samples revealed complex genetic and evolutionary relationships between the poliovirus strains. The generation of various coevolving genetic lineages incorporating different mutations was observed at early stages of virus excretion. The main driving force for genetic diversity appeared to be the selection of mutations at attenuation sites, particularly in the 5' noncoding region and the VP1 BC loop. Recombination between virus strains from the two main lineages was observed between days 63 and 88. Genetic heterogeneity among plaque-purified viruses at each time point seemed to decrease with time, and only viruses belonging to a unique genotypic lineage were seen from day 105 after vaccination. The relevance of vaccine-derived poliovirus strains for disease surveillance and future polio immunization policies is discussed in the context of the Global Polio Eradication Initiative.     

Assessment of Genetic Diversity in Thai Isolates of Pyricularia grisea by Random Amplification of Polymorphic DNA

One hundred and seventy‐four isolates of Pyricularia grisea were collected from various hosts such as barley, rice, weed and wild rice in Thailand. Seven arbitrary decamer primers from the set of University of British Columbia were employed and nine lineages were classified. Lineages B, C and H were predominant, contributing up to 70% of total pathogens in this study. Analysis showed that the distribution of each lineage differs from the predominant lineages across Thailand in such that other lineages were restricted in particular area. For instance, lineage A was limited only in southern Thailand, whereas wide distribution of lineages B and C reflected an influence of both biological and physical effects on pathogen variation. Principal component analysis resulted in a total of four groups of blast pathogen with small distinctions between barley‐, rice‐, weed‐ and wild rice‐infected blast. Bridging relationships occurred among border isolates of weed and rice blast suggesting a chance of migrations between hosts. Higher diversity was observed in northern, north‐eastern and central Thailand while eastern and southern parts were rather low. Genetic diversity indices elucidated an abundance of pathogen lineages existing in northern Thailand suggesting that it should be the centre of diversity.     

Molecular epidemiology of rabies in Thailand.

For the purpose of making clear the dynamics of rabies viruses that are prevalent among dogs in Asia, especially Thailand, nucleoprotein (N) genes of isolates derived from Thailand were partially sequenced, and a phylogenetic analysis was performed on the basis of the sequencing data. Firstly, all 27 isolates from Thailand belonged to one group that was distantly related to an isolate from China and was separated into at least six lineages. On the other hand, the isolate from Japan was related to viruses from the Arctic. Secondly, in order to analyze the diversity of the N gene more conveniently, restriction fragment length polymorphism (RFLP) analysis was performed on the N gene of 27 isolates from Thailand. The RFLP analysis could distinguish the lineages of each isolate, and the lineages of additional 34 isolates were deduced by this method. On examination of the geographical distribution of the six lineages, based on the results of phylogenetic and RFLP analyses, it was clear that infection cycles of the rabies virus in Thailand have tended to be maintained endemically.     

Molecular Evolution of Human H1N1 and H3N2 Influenza A Virus in Thailand, 2006–2009

Background

Annual seasonal influenza outbreaks are associated with high morbidity and mortality.

Objective

To index and document evolutionary changes among influenza A H1N1 and H3N2 viruses isolated from Thailand during 2006–2009, using complete genome sequences.

Methods

Nasopharyngeal aspirates were collected from patients diagnosed with respiratory illness in Thailand during 2006–2009. All samples were screened for Influenza A virus. A total of 13 H1N1 and 21 H3N2 were confirmed and whole genome sequenced for the evolutionary analysis using standard phylogenetic approaches.

Results

Phylogenetic analysis of HA revealed a clear diversification of seasonal from vaccine strain lineages. H3N2 seasonal clusters were closely related to the WHO recommended vaccine strains in each season. Most H1N1 isolates could be differentiated into 3 lineages. The A/Brisbane/59/2007 lineage, a vaccine strain for H1N1 since 2008, is closely related with the H1N1 subtypes circulating in 2009. HA sequences were conserved at the receptor-binding site. Amino acid variations in the antigenic site resulted in a possible N-linked glycosylation motif. Recent H3N2 isolates had higher genetic variations compared to H1N1 isolates. Most substitutions in the NP protein were clustered in the T-cell recognition domains.

Conclusion

In this study we performed evolutionary genetic analysis of influenza A viruses in Thailand between 2006–2009. Although the current vaccine strain is efficient for controlling the circulating outbreak subtypes, surveillance is necessary to provide unambiguous information on emergent viruses. In summary, the findings of this study contribute the understanding of evolution in influenza A viruses in humans and is useful for routine surveillance and vaccine strain selection.