Previous Article | Next Article ![]()
Journal of Virology, June 2008, p. 5494-5500, Vol. 82, No. 11
0022-538X/08/$08.00+0 doi:10.1128/JVI.02728-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Armed Forces Research Institute of Medical Sciences, 315/6 Rajvithi Road, Bangkok, Thailand 10400,1 Center for Infectious Disease Dynamics, Department of Biology, Mueller Laboratory, The Pennsylvania State University, University Park, Pennsylvania 16802,2 Fogarty International Center, National Institutes of Health, Bethesda, Maryland 20892,3 Center for Infectious Disease and Vaccine Research, University of Massachusetts Medical School, Worcester, Massachusetts 01655,4 State University of New York Upstate Medical University, Syracuse, New York 132105
Received 22 December 2007/ Accepted 28 February 2008
|
|
|---|
|
|
|---|
Over the last 20 years, many studies have documented the extent and structure of genetic variation in all four DENV serotypes (reviewed in reference 11). In general, these studies have revealed that each of the four serotypes contains a number of phylogenetically distinct "subtypes" (or genotypes), the genetic structures of which clearly reflect a complex pattern of viral gene flow (migration among locations) coupled with population subdivision. Most notably, some DENV genotypes appear to be restricted to specific localities, commonly South East Asia, while others have more-cosmopolitan distributions across the tropical and subtropical world. Whether these differences in geographic distribution reflect underlying differences in viral fitness, manifest as differences in transmissibility and/or virulence, is still a subject of debate.
Despite the growing database of partial and complete DENV genome sequences (17) and the increasing number of studies addressing various aspects of the molecular epidemiology of DENV infection, none to date have considered the nature of viral genetic diversity on a scale of populations separated by only a few kilometers and sampled over a time period of months. However, such studies of DENV "microevolution" are essential to obtain a full understanding of the mechanisms shaping DENV evolution. Further, given forthcoming phase III clinical trials of DENV vaccines, which may be imperfect in their protection against all four serotypes and hence shape viral evolution (7), it is critical to document the extent and structure of DENV genetic variation at vaccine trial sites and the pace of viral evolution within a limited time frame. To this end, we performed a detailed genetic analysis of the envelope (E) genes of DENV serotype 1, 2, and 3 isolates collected in Kamphaeng Phet, Thailand, during 2001 as part of a school absence-based prospective study.
|
|
|---|
38°C. Laboratory assays, including an immunoglobulin M (IgM)/IgG enzyme-linked immunosorbent assay (ELISA) and PCR, were used to diagnose dengue as previously described (5, 6).
![]() View larger version (48K): [in a new window] |
FIG. 1. Map of study area in subdistrict of Muang, Kamphaeng Phet, Thailand, showing the locations of participating schools.
|
Viral isolation and sequencing. For all suspected acute dengue cases, DENV serotype-specific PCR was performed on the first serum sample obtained for acute-illness evaluation. Reverse transcription-PCR (RT-PCR) was performed according to the protocol of Lanciotti et al. (14), with the following modifications. Avian myeloblastosis virus reverse transcriptase (Promega, Madison, WI) was used in the first-round RT-PCR. The concentrations of the primers used in the RT-PCR and nested reactions were reduced from 50 pmol to 12.5 pmol per reaction, and the number of nested PCR amplification cycles was increased to 25.
For DENV PCR-positive samples, virus isolation was performed with C6/36 cells and/or Toxorhynchites splendens mosquitoes as previously described (13, 16). DENV serotype identification was performed using a dengue serotype-specific enzyme immunoassay as previously described (18, 20). Viral RNA for sequencing was extracted from patient samples first amplified in T. splendens mosquitoes, followed by one passage in C6/36 cells by use of a QIAamp viral RNA mini kit (Qiagen, Germany) according to the manufacturer's instructions. RT was performed using random hexamer oligonucleotides with the SuperScript first-strand synthesis system (Invitrogen) according to the manufacturer's instructions. The DNA fragments of the envelope gene regions of 9 DENV-1, 40 DENV-2, and 18 DENV-3 isolates were amplified by PCR using 5 µl of cDNA in a 50-µl reaction mixture containing 0.3 mM deoxynucleoside triphosphates, 2.5 U AmpliTaq DNA polymerase (Applied Biosystems), 1x PCR buffer, 1.5 mM MgCl2, and 15 pmol of each forward and reverse primer. The PCR mixtures of DENV-1 and -3 were subjected to 1 cycle of 95°C for 5 min; 35 cycles of 94°C for 30 s, 50°C for 1 min, and 72°C for 2 min; and 1 cycle of 72°C for 15 min. The PCR mixtures of DENV-2 were subjected to the same thermal conditions as those of DENV-1 and -3, except for the annealing temperature, which was changed to 55°C. The PCR-amplified DNA fragments were purified using QIAquick PCR purification kits (Qiagen) according to the manufacturer's instructions. Purified DNA fragments were used for sequencing.
Sequencing reactions were performed by using a DYEnamic ET dye terminator sequencing kit (GE Healthcare Bio-Sciences) according to the manufacturer's instructions. The sequencing primers are available upon request. The sequencing products were cleaned by standard precipitation before sequencing with a MegaBACE 500 automated DNA sequencer (GE Healthcare Bio-Sciences). Overlapping nucleic acid sequences were combined for analysis and edited by using SEQUENCHER software (Gene Code Corporation).
Unfortunately, the number of viruses available for sequencing was a relatively small percentage of the viruses circulating in our cohort during the 2001 dengue season. In all, 244 infections, including inapparent infections, were detected in our cohort. Of the overt infections, only 86 had virus detectable by PCR, while the remaining 38 were detected based on a serological response consistent with acute dengue. Hence, we were able to sequence 78% of the PCR-positive samples but only 67 of 244 sequences (27%).
Phylogenetic analysis. The following sets of DENV E-gene sequences from Kamphaeng Phet were compiled for evolutionary analysis: for DENV-1, 9 sequences (genotype I) from three schools; for DENV-2, 40 sequences (genotype Asian I) from nine schools; and for DENV-3, 18 sequences (genotype II) from three schools. To place these isolates within the wider context of DENV evolution, particularly that from South East Asia, we also compiled, from GenBank, phylogenetically representative "background" data sets of sequences from each of the respective genotypes, many of which also came from Thailand. This resulted in final data sets of the following sizes in the phylogenetic analyses: for DENV-1, 116 sequences (1,485 nucleotides [nt]); for DENV-2, 168 sequences (1,485 nt); and for DENV-3, 124 sequences (1,479 nt).
Phylogenetic trees for all three data sets were inferred using the maximum likelihood (ML) method available in PAUP* (19), with the best-fit model of nucleotide substitution (usually GTR+I+
4) determined using the MODELTEST program (15) and employing tree bisection-reconnection branch swapping. Bootstrap resampling (1,000-replicate neighbor-joining trees under the ML substitution model) was used to determine the phylogenetic support for individual nodes. All parameter values are available from the authors on request.
Analysis of spatial and temporal structure. In the case of DENV-2, the larger number of E-gene sequences allowed us to determine whether there were distinct spatial and temporal components to the evolutionary patterns observed. This was achieved by using a parsimony approach which had previously proven highly informative in a study of DENV in the Caribbean (4). In the case of spatial structure, we tested whether there was significant clustering according to the school from which the virus was isolated. In the case of temporal structure, we tested whether viruses were clustered according to their months of isolation. In each of these analyses, DENV-2 sequences were assigned a character state dependent on the school (spatial) or month (temporal) of origin. Given the ML phylogeny for these sequences (determined as described above) and the above-defined isolate states, the minimum number of changes in character state needed to give rise to the observed distribution of states was then estimated using parsimony (with all ambiguous changes excluded). To determine the expected number of changes under the null hypothesis of complete mixing among states by space or time, the states of all isolates were randomized 500 times, and for each randomization, the number of changes in state was calculated in the manner described above. The difference between the mean numbers of observed and expected changes for each pair of states indicates the level of geographic or temporal isolation, with statistical significance calculated by comparing the total number of observed state changes to the number expected under random mixing. All these analyses were performed using PAUP*.
Analysis of selection pressures. Overall and site-specific selection pressures in the E genes of all three serotypes circulating in Kamphaeng Phet (excluding other Asian sequences) were measured as ratios of nonsynonymous (dN) to synonymous (dS) substitutions per site (dN/dS), estimated using the single-likelihood-ancestor-counting and random-effect-likelihood methods. Both methods utilized the GTR substitution model with input phylogenetic trees inferred using the neighbor-joining method available at the Datamonkey Web facility (12).
Nucleotide sequence accession numbers. All DENV gene sequences have been submitted to GenBank and assigned accession numbers EU117304 to EU117312, EU117313 to EU117352, and EU117353 to EU117370 for DENV-1, -2, and -3, respectively.
|
|
|---|
2,000 primary school age children in Kamphaeng Phet, Thailand. Overt illness was detected using school absences to identify sick children, and DENV infections were confirmed using IgM/IgG ELISA, hemagglutination inhibition, and RT-PCR (5). Inapparent infections were identified based on a fourfold rise in hemagglutination inhibition between serial serum samples. Over the 5-year study period, 312 inapparent and 304 overt infections were identified. The overt infections were detected in children attending every school in the study, and all four serotypes were identified by RT-PCR (1, 5, 6). The highest incidence of disease occurred in 2001, with 120 inapparent infections and 124 cases of overt disease (Table 1). Because of the high dengue incidence during 2001, samples were collected throughout the year, with the majority identified during the peak dengue season of June to November as part of the school absence study and the remainder identified from hospitalization of study participants (Table 2). Through ELISA typing of clinical isolates or RT-PCR of serum samples, we were able to identify the infecting DENV serotype in 86 (69%) symptomatic cases; the remaining 38 (31%) were identified as acute DENV infection by serology only. DENV-2 was the predominant serotype, identified in 52 of the 86 (60%) cases, followed by DENV-3 (25 cases; 29%) and DENV-1 (9 cases; 10%), collected in 10 of the 12 schools.
|
View this table: [in a new window] |
TABLE 1. Numbers of DENV infections in school absence-based cohort during 2001
|
|
View this table: [in a new window] |
TABLE 2. Numbers of viral isolates sequenced according to school and month of sampling
|
ML phylogenetic trees for the DENV-1, DENV-2, and DENV-3 E-gene sequences from Kamphaeng Phet, combined with the background isolates, are presented in Fig. 2, 3, and 4, respectively. In the case of DENV-1, viral isolates fall into two distinct clades, separated by other Asian viruses. That the DENV-1 isolates from school 11 fall in both clades provides compelling evidence for the independent entry of genetically distinct viruses into this spatially restricted region. A similar pattern was observed in DENV-2. Here, those viruses sampled from schools 1, 5, 8, 9, and 10 formed school-specific clusters (although only those from school 1 were clearly phylogenetically distinct from those from the other schools), while multiple genetic variants were seen in schools 2, 4, 11, and 12 (with one highly divergent variant observed in school 11). Finally, in the case of DENV-3, single clades were observed in schools 4 and 10, while two clades were again observed in school 11, although with weak bootstrap support. Notably, children from school 11 come from a densely populated urban environment at the center of the study area, close to the major transportation throughway and a public health clinic, which may have contributed to the appearance of multiple viral variants in this school.
![]() View larger version (21K): [in a new window] |
FIG. 2. ML tree of selected isolates of genotype I of DENV-1. Different schools from within Kamphaeng Phet are designated by different colors. All horizontal branch lengths are drawn to scale according to the numbers of nucleotide substitutions per site, with bootstrap support values shown next to relevant nodes. The tree is midpoint rooted for purposes of clarity only.
|
![]() View larger version (29K): [in a new window] |
FIG. 3. ML tree of selected isolates of the Asian I genotype of DENV-2. Different schools from within Kamphaeng Phet are designated by different colors. All horizontal branch lengths are drawn to scale according to the numbers of nucleotide substitutions per site, with bootstrap support values shown next to relevant nodes. The tree is midpoint rooted for purposes of clarity only.
|
![]() View larger version (30K): [in a new window] |
FIG. 4. ML tree of selected isolates of genotype II of DENV-3. Different schools from within Kamphaeng Phet are designated by different colors. All horizontal branch lengths are drawn to scale according to the numbers of nucleotide substitutions per site, with bootstrap support values shown next to relevant nodes. The tree is midpoint rooted for purposes of clarity only.
|
To determine the extent of spatial clustering more rigorously, we examined DENV-2 in more detail as this represented our best-sampled serotype. A parsimony-based analysis confirmed that there is a highly significant clustering of viral isolates by school (P < 0.005). However, despite this strong genetic segregation, there is also clear evidence (P < 0.005) of viral migration between schools 8 and 9. Notably, these two schools are geographically the closest in the study population, separated by only 5 km, and the furthest from the city, at distances of 36 and 33 km from the field site, respectively (Fig. 1). Further, the villages that feed these schools have low population densities, are highly rural, and possess poor road infrastructure, making access to this area more difficult than access to other regions of Kamphaeng Phet. As such, our data suggest that patterns of viral gene flow are determined by local geographic and economic variables. The villages that make up the Kamphaeng Phet region are diverse, with the most developed and densely populated villages lying close to the city center, where frequent travel into and out of the city is common. In contrast, villages become less developed and densely populated further from the city center, in turn restricting travel and resulting in the genetic isolation of viruses from schools 8 and 9.
In contrast, these DENV-2 sequence data exhibited no significant clustering by month of sampling, indicating that there was little replacement of viral lineages at this level of spatial and temporal resolution. Similarly, at the level of the E-gene sequences analyzed, there was little evidence for in situ evolution, such that viruses sampled from the same school over time periods ranging from days to months exhibited few mutational differences (in some cases only one mutation over a 6-month sampling period). Such relative intraschool stability sits in marked contrast to the often strong genetic differentiation among schools. This observation, coupled with the strong population subdivision exhibited by each school, indicates that the importation of DENV into the study population, rather than in situ evolution within school catchment areas, was the most important factor shaping viral genetic diversity on this spatial scale.
To determine the selection pressures acting on the DENV E-gene sequences sampled from Kamphaeng Phet, we compared the relative numbers of nonsynonymous (dN) and synonymous (dS) substitutions per site. This analysis provided no evidence for positive selection acting on any nucleotide site, although all dN/dS methods are inherently conservative. Indeed, the overall picture obtained was that of strong purifying selection, with mean dN/dS values of 0.084, 0.065, and 0.008 for DENV-1, DENV-2 and DENV-3, respectively. Such strong purifying selection has previously been observed for these serotypes in Thailand (21, 22), highlighting the predominance of genetic drift over natural selection in shaping DENV substitution dynamics in the short term. However, the reasons for the greater selective constraints in DENV-3, which is characterized by very low dN/dS values, are unclear. Further, as random genetic drift is the most common process determining substitution dynamics, allele frequencies are not expected to change greatly over the period of sampling, even given the high mutation rate of DENV (3).
Overall, these results have mixed implications for how DENV populations might respond to imperfect vaccination in the near future. DENV sampled on a short spatial and temporal scale evidently exhibits remarkably high levels of genetic diversity, thereby providing the raw material for adaptive evolution. However, this variation is more often generated by migration than by mutation accumulation, with purifying selection the dominant evolutionary force.
Funding was partially provided by the U.S. Military Infectious Disease Research Program, Ft. Detrick, MD, and the National Institutes of Health (NIH-P01-AI34533).
The opinions and assertions contained herein are the personal views of the authors and are not to be construed as official or reflecting the views of the Armed Forces Research Institute of Medical Sciences, the U.S. Department of the Army, the U.S. Department of Defense, or the National Institutes of Health.
Published ahead of print on 26 March 2008. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»