Previous Article | Next Article ![]()
Journal of Clinical Microbiology, April 2006, p. 1519-1529, Vol. 44, No. 4
0095-1137/06/$08.00+0 doi:10.1128/JCM.44.4.1519-1529.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
*
G. Palacios,2,
O. Jabado,2
N. Reyes,1
M. Niedrig,3
J. Gascón,4
M. Cabrerizo,1
W. I. Lipkin,2 and
A. Tenorio1
Laboratorio de Arbovirus y Enfermedades Víricas Importadas, Centro Nacional de Microbiología, Instituto de Salud Carlos III, Madrid, Spain,1 Jerome L. and Dawn Greene Infectious Diseases Laboratory, Columbia University, New York, New York,2 Robert Koch Institute, Berlin, Germany,3 Centro de Salud Internacional, Hospital Clinic (IDIBAPS), Barcelona, Spain4
Received 6 September 2005/ Returned for modification 7 October 2005/ Accepted 8 January 2006
|
|
|---|
|
|
|---|
The increasing incidence of dengue epidemics in the Americas, Southeast Asia, the Indian subcontinent, and the Western Pacific, along with the growing rate of DHF and DSS cases, is a global public health problem. Urbanization, overpopulation, crowding, poverty, a weakened public health infrastructure, and the globalization of trade and travel have been implicated as causes for this increase (3, 9).
Dengue viruses comprise four different serotypes, and evolutionary and epidemiological studies have facilitated clustering of dengue viruses into genotypes within serotypes (14, 30). It was posited that dengue hemorrhagic fever is a result of the immune enhancement developed after subsequent infection with two different serotypes (10). This is also a plausible mechanism to explain the observed highest risk of DHF and DSS in primary infections in infants with heterotypic immunity due to waning maternal dengue virus antibodies (17). It has also been suggested that virulent strains present an increasing replication rate in human target cells (29); therefore, the virus itself might play an important role in disease severity. Moreover, some genotypes have been associated with DHF epidemics, while others genotypes have been responsible for mild disease in certain areas while causing secondary infections (11, 37, 38). It was posited that these genotypes may present different antigenic epitopes and that these structural differences correlate with pathogenesis (18, 21). Significant differences in disease severity are observed in secondary infections, depending on the origin of the strains (4, 18, 21). Thus, the presence of a specific genotype in an area could be associated with an increased risk for the appearance of DHF and DSS cases.
Accordingly, a surveillance methodology able to identify the presence or introduction of dengue genotypes in areas of endemicity is of great interest for pathogenesis studies and for the establishment of local epidemiological control policies.
Here we describe a rapid, sensitive method for the genotyping of DENV-1, in which a short fragment in the carboxyl terminus of the E gene (E-COOH) is amplified by PCR, sequenced, and analyzed in pairwise comparisons with prototypic DENV sequences representing DENV-1 strains. Application of this method to serum samples obtained from DENV-1-infected travelers who had recently returned to Europe from India suggested the presence of two new DENV-1 lineages. Sequence analysis of the entire E gene confirmed this observation and supported the utility of the method.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. New DENV1 strains included in the study
|
The samples were also subjected to analysis by a specific DENV-1 nested RT-PCR to amplify the complete E gene. This protocol was performed by using the buffers described for the E-NS1 amplification nested RT-PCR, except for the use of specific primers (primers EGENE1-S, EGENE-R, EGENE1-SS, and EGENE/NS1-RR) that were designed based on published DENV-1 sequences by computer-assisted analysis (MACAW version 32 software, 1995; NCBI). The RT-PCRs were carried out by using an initial reverse transcription step at 48°C for 45 min, followed by a denaturation and HotStar Taq polymerase activation step (94°C, 15 min) and 40 cycles of denaturation (94°C, 30 s), primer annealing (56°C, 2 min), and primer extension (72°C, 2 min). A final incubation was carried out at 75°C for 5 min. In the second amplification reaction, the samples were subjected to a denaturation step (94°C, 5 min), followed by 40 cycles of denaturation (94°C, 30 s), primer annealing (59°C, 2 min), and primer extension (72°C, 2 min) and a further extension step at 75°C for 5 min.
Table 2 shows the sequences and the respective primer positions in the viral genome used for partial (328-nt) or whole (1,700-nt) E gene amplification. The sensitivities were 100 and 10,000 copies, respectively.
|
View this table: [in a new window] |
TABLE 2. Primers used for amplification and sequencing of the complete E gene and the E-NS1 junction of DENV-1
|
Original sequence data were first analyzed by use of CHROMAS software (version 1.3, 1996; C. McCarthy, School of Biomolecular and Biomedical Science, Faculty of Science and Technology, Griffith University, Brisbane, Queensland, Australia); the forward and reverse sequence data for each sample were aligned by using the program SEQMAN (DNASTAR Inc. Software, Madison, Wis.). The consensus sequence was compared and aligned to other samples or DNA database sequences using the program CLUSTAL X, version 1.83 (34).
Sequence analysis of amplified products. A set of 156 DENV-1 sequences (216 nt in length) comprising those from 30 new sequenced strains and 126 collected from GenBank (updated to July 2005) were used for phylogenetic analysis. To test the findings observed with the short fragment, the entire E protein-coding sequence (1,485 nt) was recovered directly from 16 clinical samples by nested RT-PCR and compared in pairwise analyses with DENV-1 complete E 116 sequences collected from GenBank.
Phylogenetic analysis. Phylogenetic analysis was performed by using the best model of nucleotide substitution (according to Modeltest [27] and Tamura and Nei [33], with correction for the proportion of invariable sites of 0.4434 and a gamma distribution of 1.4212). Programs from the MEGA package (version 3) (19) were used to produce phylogenetic trees, reconstructed by the neighbor-joining method. The statistical significance of a particular tree topology was evaluated by bootstrap resampling of the sequences 1,000 times. Pairwise comparisons of the DENV-1 database were done by global alignment by using the algorithm of Needleman and Wunsch (23), implemented by a program from EMBOSS, the European Molecular Biology Open Software Suite (28). An automated program that performs the same analysis is available at http://www.greeneidlab.columbia.edu and http://aevi.isciii.es.
Recombination event detection. Systematic screening for the presence of recombination patterns was achieved by using the nucleotide alignments and the Recombination Detection Program (RDP; Darren Martin) (22). The algorithms Bootscan (31), MaxChi (32), Chimaera (26), LARD (Likelihood Assisted Recombination Detection) (13), and Phylip Plot (6) were also used for detection of recombination.
Nucleotide sequence accession numbers. The GenBank accession numbers of the nucleotide sequences determined in this study are DQ016630 to DQ016659.
|
|
|---|
Phylogenetic analysis using the complete envelope gene. The phylogenetic tree obtained by analysis of a complete data set of all available envelope gene sequences is shown in Fig. 1. The analysis not only allowed the classification within known DENV-1 genotypes, America-Africa (AMAF), Malaysia (MAL), Thailand (THAI), Asia (ASIA), and South Pacific (SP), but also suggested the existence of lineages with distinctive geographical and temporal relationships (16, 20).
![]() View larger version (33K): [in a new window] |
FIG. 1. The complete set of nucleotide sequences of the dengue virus 1 glycoprotein E gene (1485 nt, 116 sequences) was aligned by use of the CLUSTAL W program. Phylogenetic analysis was performed by using the best model of nucleotide substitution (according to Modeltest [27] and Tamura and Nei [33], with correction for a gamma distribution of 1.4212) and the neighbor-joining method to reconstruct the phylogenetic tree (MEGA, version 3.0, software package). The statistical significance of the phylogenies constructed was estimated by bootstrap analysis with 1,000 pseudoreplicate data sets. Strains are denoted by number, country of isolation, and year of isolation; genotypes are indicated, and names are given to the right.
|
![]() View larger version (42K): [in a new window] |
FIG. 2. Phylogenetic tree of the genotype AMAF. Lineages are indicated. Significant bootstrap values ( 60) are indicated. The analysis method was the same as that reported in the legend to Fig. 1.
|
![]() View larger version (43K): [in a new window] |
FIG. 3. Phylogenetic tree of the genotype ASIA complete E gene. Lineages are indicated. Significant bootstrap values ( 60) are indicated. The analysis method was the same as that reported in the legend to Fig. 1.
|
![]() View larger version (32K): [in a new window] |
FIG. 4. Phylogenetic tree of the genotype SP complete E gene. Lineages are indicated. Significant bootstrap values ( 60) are indicated. The analysis method was the same as that reported in the legend to Fig. 1.
|
Recombination analysis. Recombination has been reported in the envelope gene in all DENV serotypes (35, 41). Because recombination events may occur at various positions in the envelope gene, phylogenetic trees based on partial envelope sequence data may be misleading. To address this concern we searched for evidence of recombination using the full envelope sequence. We detected previously reported recombination between lineages AMERICA-1 and ASIA-3 strains in S275-Singapore-1990, and between lineages AMERICA-1 and SP-3 in strains Caribbean 495-Mexico1985, Caribbean 925-Mexico1985, CV1636-Jamaica1977, and AHF82-80-Thailand1982 (35). Additionally, we found evidence of a new recombination event in the GD23-95-China-1995 sequence. This virus appears to be a recombinant of the ASIA-3 and SP-1 lineages (Fig. 5). All the recombinant strains appeared as outliers in our phylogenetic tree (Fig. 2 and 4).
![]() View larger version (23K): [in a new window] |
FIG. 5. Phylogenetic support for recombination in DENV-1. The neighbor-joining tree supporting the contrasting phylogenetic position of strain GD23-95-China 1995 in the E protein is shown. The 5' end of the sequence (until position 371 of the E gene) groups this strain with genotype ASIA (upper left panel); the 3' end sequence clustered instead with genotype SOUTH PACIFIC (upper right panel). Diversity plots for putative recombinants and their parents. (Middle left panel) Likelihood surface plot generated by LARD (Likelihood Assisted Recombination Detection). The partition that separates trees with the greatest likelihood score is the most probable recombination breakpoint. (Middle right panel) Pairwise identity plot obtained by the RDP method with parent strains A88-1988-Indonesia (SOUTH PACIFIC) and GD05_99-1995-CHI (ASIA). (Lower left panel) Chimaera results. A chi-square P value plot used to confirm the RDP result. When plotted along the length of the alignment, peaks in these 2 values indicate potential recombination breakpoints. (Lower right panel) DNA distance plot method. Pairwise distances were calculated for each window by using DNADIST (a component of the PHYLIP package), phylogenetic trees were obtained, and bootstrap support for each association is plotted against the position in the alignment of the window's center. The horizontal axis in all graphs is the nucleotide number counting from the 5' end of the E sequence. The window size and increment were 150 and 3 bp, respectively, in each case.
|
350 nt into the E gene and the 3' end breakpoint occurs in the NS1 gene, the analysis of the 3' portion of the E gene should allow detection of the inserted recombinant sequence. Consistent with this observation, strain GD23-95-China-1995, a recombinant of the ASIA-3 and SP-1 lineages, groups with the SP-1 lineage. Similarly, the cluster that includes recombinant strains from the AMERICA-1 and SP-3 lineages appears as a member of the SP-3 lineage; and the S275-Singapore strain, a recombinant between AMERICA-1 and ASIA-3 lineages, appears as a member of the AMERICA-1 lineage. Pairwise analysis by use of local alignments. An all-against-all pairwise sequence comparison was done with the DENV-1 sequences to assess the potential for the establishment of a simple program for the classification of DENV-1 genotypes and lineages similar to programs that we have created for mumps virus (25) and adenovirus (1). Sequence analysis of the 216-nt region at the carboxyl terminus amplified by PCR with consensus primers SEQAS/E/NS1-SS and SEQBR/EGENE/NS1-RR (Table 2) was sufficient to allow determination of the genotype and the lineage (Fig. 6). The validity of this method was confirmed by analysis of variance by comparing the scores of the sequence comparisons within genotypes or lineages to the scores of comparisons between genotypes. A software tool has been developed to classify DENV-1 sequences by introducing the target sequence (either the complete sequence or the carboxyl E gene sequence) and is available at http://www.greeneidlab.columbia.edu and http://aevi.isciii.es.
![]() View larger version (29K): [in a new window] |
FIG. 6. Pairwise analysis of the sequences of four dengue virus strains randomly chosen from the sequences obtained in this study. Samples are SC418VI03, SP genotype, lineages SP-1 and ASIA-4; CNM60VI05, ASIA genotype, lineage ASIA-2; WI22VI05, AMAF genotype, lineage INDIA-1; and SC416VI03, genotype AMAF, lineage INDIA-2. A pairwise Needleman-Wunsch score for the test sequence against all other members of the database was determined. The average score (with the corresponding standard deviation, maximum, and minimum values) against the members of each genotype is plotted.
|
|
|
|---|
Three (2, 15, 24) to five (7, 29) different genotypes of DENV-1 have been proposed. Our analysis of complete E gene sequences confirmed the presence of five genotypes and the clustering of sequences below the genotype level that correlate with the geographical origin and/or time of isolation and appear to represent distinct lineages with strong bootstrap support (Fig. 1).
A correlation between genotypes and pathogenesis had been suggested. It would be of interest to observe if the potential to cause DHF is a property shared at the genotype level or at the lineage level. The AMAF genotype has been described as a genotype with a high epidemiological impact due to its spread and its potential to cause DHF (29). Remarkably, two new clusters of Indian strains belonging to the genotype AMAF (lineages INDIA-1 and INDIA-2) were observed in our analysis. One lineage (INDIA-1) is related to American strains, while INDIA-2 is closer to the African and South Asiatic strains. Two different lineages are currently circulating in the subcontinent. From our data it could be inferred that this genotype is currently circulating in extensive areas both in the Americas (the AMERICA-1 lineage) and on the Indian subcontinent (the INDIA-1 and INDIA-2 lineages).
Given the geographical or temporal association among the sequences belonging to the same lineage, a classification based on lineages facilitates surveillance and tracking of dengue virus isolates. For example, in genotype ASIA, three lineages (the ASIA-1, ASIA-2, and ASIA-3 lineages) have been detected. The ASIA-2 lineage had been detected in Myanmar, Cambodia, and Thailand. We detected two members of ASIA-2 in patients returning from Kenya; this may suggest the introduction of this lineage into East Africa.
During the course of reviewing published sequence data, we discovered evidence of recombination between the ASIA-3 and the SP-1 lineages in GD23-95-China 1995. Since this virus appears to remain apart from the rest of the clusters (Fig. 5) and was isolated in 1995, it is possible that this lineage did not become established. Although recombination is accepted in Japanese encephalitis and St. Louis encephalitis viruses (8, 36), recombination in DENV has been controversial. The original reports of recombination in DENV by Holmes et al. (12) and Worobey et al. (41) have been challenged (7, 29). The only unequivocal way to resolve the scientific question with respect to GD23-95-China 1995 would be to obtain the isolate to independently confirm the results of sequence and phylogenetic analyses.
The use of a consensus PCR based on a short fragment of the E protein coupled to a similarity analysis by using a pairwise alignment method demonstrate that the DENV-1 genotype and lineage can be accurately defined. This method offers some additional advantages in microbiology laboratories, in contrast to the application of classical methods of molecular epidemiology (e.g., sequencing of the complete envelope from cell culture isolates). Among those advantages, we count the use of clinical samples as the source of viral RNA; the amplification of a short sequence to identify the serotype, the genotype, and the lineage, thus facilitating the acquisition of information in a very short time; the application of pairwise comparison analysis and the proposed online typing tool to identify DENV-1 genotypes, which minimizes the need for extensive phylogenetic analysis; and the use of an automatically updated database for comparison and a global analysis (e.g., all available DENV sequences are used each time), which minimizes the possibility of the most common mistake in phylogenetic analysis (e.g., use of a biased set of sequences for comparison). Only in the case that a new strain could not be clearly associated with a described genotype and/or lineage should the complete E gene be analyzed to identify the appearance of a new genotype. The detection and characterization of two new lineages of DENV-1 in India and the description of the spread of genotype ASIA, lineage ASIA-2, in East Africa clearly show the usefulness of this approach. This approach has been demonstrated to be equally useful for the genotyping of the four dengue virus serotypes (data not shown).
The use of this methodology for the identification of dengue virus infection could greatly contribute to the acquisition of dengue epidemiological data worldwide, detection of the appearance and spread of new genotypes and lineages, and detection of the spread and circulation of those already described and, more deeply, to a basic understanding of dengue virus pathogenesis.
M.N. and A.T. are ENIVD participants . Other ENIVD participants include P. Cassinotti and D. Schultze (Institute für Klinische Mikrobiologie und Immunologie, St. Gallen, Switzeralnd), M. Van Esbroek (Institute of Tropical Medicine, Antwerp, Belgium), A. Fomsgaard and L. Vinner (Statens Serum Institute, Copenhagen, Denmark), and M. Grandadam (Unité de Virologie Tropicale Parc du Pharo, Marseille Armées, France). J.G. is TropNetEurop participant. Other TropNetEurop participants include S. Puente (Hospital Carlos III, Madrid, Spain), O. Whichmann (Institute for Tropical Medicine, Berlin, Germany), M. Schunk (Medizinische Klinik Klinikum der Universität, Munich, Germany), and R. López-Vélez (Hospital Ramón y Cajal, Madrid, Spain). M.N. is a representative of the ENIVD Dengue Study Workgroup, and J.G. is a representative of the TropNetEurop Dengue Study Workgroup. A.T., J.G., and C.D., as well as S. Puente, are researchers of the RICET Network (Red de Investigación Cooperativa en Enfermedades Tropicales C03/04).
C.D. is under contract by an agreement between the Public Health Division of the Spanish Ministry of Health (DGSP-MSC) and the Instituto de Salud Carlos III (ISCIII) for the development of the Hemorrhagic Viral Fevers Surveillance and Control Program. This work has received financial support from the Instituto de Salud Carlos III (ISCIII) through research project grants (grants MPY 1194/02 and C03/04). G.P. and W.I.L. are supported by the Ellison Medical Foundation and the National Institutes of Health (grants AI 51292 and U54 AI57158 to W.I.L.).
Both authors contributed equally. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»