Emergence and global spread of a dengue serotype 3, subtype III virus. (Research).

Over the past two decades, dengue virus serotype 3 (DENV-3) has caused unexpected epidemics of dengue hemorrhagic fever (DHF) in Sri Lanka, East Africa, and Latin America. We used a phylogenetic approach to evaluate the roles of virus evolution and transport in the emergence of these outbreaks. Isolates from these geographically distant epidemics are closely related and belong to DENV-3, subtype Ill, which originated in the Indian subcontinent. The emergence of DHF in Sri Lanka in 1989 correlated with the appearance there of a new DENV-3, subtype III variant. This variant likely spread from the Indian subcontinent into Africa in the 1980s and from Africa into Latin America in the mid-1990s. DENV-3, subtype III isolates from mild and severe disease outbreaks formed genetically distinct groups, which suggests a role for viral genetics in DHF.

**********

Arthropod-borne viruses are responsible for the emergence of unexpected diseases in humans, as illustrated by the identification of West Nile virus encephalitis in the American hemisphere in 1999 (1). The emergence of a new disease is often attributable to the transport of a pathogen (as in the case of West Nile virus) or changes in the evolution or ecology of a native pathogen that hitherto caused mild or no disease in humans (2,3). We studied unexpected outbreaks of dengue hemorrhagic fever (DHF) in Sri Lanka, East Africa, and Latin America caused by dengue serotype 3 (DENV-3) virus.

Most persons infected with dengue viruses are asymptomatic or develop dengue fever (DF). DHF and dengue shock syndrome (DSS), which can be fatal, develop in a minority of infected persons. The pathogenesis of DHF is poorly understood, although factors such as age and previous exposure to dengue infections increase the risk for severe disease (4). Epidemiologic studies point to particular DENV strains being more virulent than others (5-8). For example, the dengue genotypes endemic to Central and South America have caused mild disease, while the Asian genotypes introduced to the region have led to DHF epidemics (9-16). Similarly, outbreaks of DHF in some Pacific islands have been traced to the introduction of Southeast Asian dengue strains (17). DENV-2 subtypes associated with mild and severe disease epidemics have distinct mutations in the E gene and 5' and 3' untranslated segments of the viral genome, although whether these mutations directly contribute to pathogenesis is unproven (18).

The distribution of DHF and DSS in Asia has been particularly puzzling. Before 1989, DHF was common in Southeast Asia but rare in the Indian subcontinent despite the circulation of all four serotypes in both regions. After 1989, this pattern of disease changed and regular epidemics of DHF were reported from several countries in the Indian subcontinent (19). Sri Lanka, in particular, experienced a dramatic and persistent increase in DHF cases (20). Epidemiologic studies of dengue in Sri Lanka have demonstrated that the intensity of virus transmission, as well as the relative abundance of each serotype, remained constant before and after the emergence of DHF (21). Thus, DHF did not emerge in Sri Lanka because of an overall increase in virus transmission or shift in serotype.

Although all four serotypes of dengue circulate in Sri Lanka, persons who have the severe form of the disease are most frequently infected with DENV-3 (20,22). Lanciotti et al. characterized the genetic relatedness of DENV-3 isolates from regions throughout the tropics and subtropics and identified four geographically distinct subtypes (23). All Sri Lankan isolates were classified as subtype III, which also includes isolates from East Africa and India, as well as recent isolates from Latin America. Because DENV-3 isolates from Sri Lanka isolated before and after 1989 (when DHF emerged) formed separate groups within subtype III, Lanciotti and colleagues postulated that a genetic shift in DENV-3 may have been responsible for the emergence of DHF (23).

In the current study, using phylogenetic methods, we analyzed DENV-3 viruses isolated from Sri Lanka for up to 10 years after the emergence of DHF to confirm the establishment of a new genotype and evaluate the roles of virus evolution and transport in establishing a new genotype. DENV-3, subtype III was introduced into Latin America in 1994 (11), and the virus has subsequently been isolated from DF and DHF outbreaks throughout Central and South America (12-16). We also examined the genetic relationships between DENV-3, subtype III isolates from Latin America, East Africa, and the Indian subcontinent. On the basis of our results, we describe the most likely scenario of events that led to the emergence of DENV-3--associated DHF in the Indian subcontinent and the Americas.

Materials and Methods

Virus Strains

The dengue virus strains sequenced for this study as well as sequences obtained from GenBank for this study are listed in Table 1. The virus isolates were obtained from the Centers for Disease Control and Prevention, Dengue Branch, Puerto Rico, and Division of Vector-Borne Infectious Diseases, Ft. Collins, Colorado; Medical Research Institute, Colombo, Sri Lanka; School of Public Health, Berkeley, California; Walter Reed Army Institute for Research, Washington, D.C.; and University of Massachusetts Medical Center, Worcester, Massachusetts.

RNA Extraction

QiaAmp Viral RNA Mini Kit (QIAGEN, Valencia, CA) was used to extract viral RNA from both the mosquito grind supernatants and infected tissue culture media following the manufacturer's protocol. Extracted RNA was stored at -70[degrees]C or immediately subjected to reverse transcription-polymerase chain reaction (RT-PCR).

RT-PCR

DENV-3 RT-PCR was carried out as described by Lanciotti (23). Primers were designed to amplify and sequence a 966-bp fragment from positions 179-1,144, encompassing part of Capsid, all of PreM, and part of the E gene sequences. The reverse primer (DEN3/735) hybridized to positions 1,189-1,171 (5'-ctcctcaggcaaaaccgct-3') and the forward primer (D1 consensus) hybridized to positions 132-159 (5'-tcaatatgctgaaacgcgcgagaaaccg-3'). The reverse primer DEN3/735 was added to extracted RNA, incubated at 85[degrees]C for 90 s, and allowed to cool to room temperature. RT was carried out for 45-60 min in 20 [micro]L of reaction mix containing 25 U avian myeloblastosis virus reverse transcriptase (Roche, Nutley, NJ), deoxynucleoside tripophosphate, Mg[Cl.sub.2] and RT buffer. PCR was performed by adding a 30-[micro]L cocktail containing D1 consensus primer, PCR buffer, and EXPAND polymerase (Roche) to the 20-[micro]L RT reaction. PCR conditions were 4 min at 94[degrees]C, 30-35 cycles of 94[degrees]C for 30 s, 54[degrees]C for 30 s, and 72[degrees]C for 90 s with 5 s/cycle added to elongation step after the first 10 cycles. We separated 5 [micro]L of the reaction products on 2% agarose gels and visualized it by ethidium bromide staining. When necessary, target bands were excised and purified by using the Qiagen QIAquick Gel Extraction kit (QIAGEN) following manufacturer's instructions. All remaining PCR reaction products were purified by using the Qiagen PCR Purification kit following the manufacturer's protocol.

DNA Sequencing

Purified PCR products were sent to the automated DNA sequencing facility at the University of North Carolina, Chapel Hill, NC. The DENV-3 sequences used in this manuscript included 40 newly determined sequences, which have been submitted to GenBank (accession nos. AF547225-AF547264).

Viral Sequence Analysis

Overlapping individual nucleic acid sequences were assembled with the aid of VECTOR NTI ContigExpress (InforMax, Inc., Bethesda, MD). Sequences were aligned and analyzed by using the following software: Clustal X (available from: URL: http://inn-prot.weizmann.ac.il/software/ClustalX.html), PAUP* (available from: URL: http://www.sinauer.com), PHYLIP (available from: URL: http://evolution.genetics.washington.edu/phylip.html), and MEGA II (available from: URL: http://www.megasoftware.net). Genetic distances were calculated by using Tamura-Nei distance algorithm with 1,000 bootstrap replicates; the trees were generated by using the Minimum Evolution method. The phylogenetic tree in Figure 1 is based on a 708-base segment, positions 437-1,145, spanning pre-M/M and a portion of the E gene. The phylogenetic tree presented in Figure 2 is based on 966-base region spanning positions 179-1,145 on the viral genome, capturing a portion of the C gene, all of pre-M/M gene, and a portion of the E gene.

[FIGURE 1-2 OMITTED]

Results and Discussion

Many investigators have used viral nucleotide sequence data and phylogenetic methods to understand genetic relationships between viruses, as well as the epidemiology of viral disease. Phylogenetic studies have shown that dengue viruses can move long distances between continents (24) as well as short distances between neighboring countries (25). Our goal was to use a phylogenetic approach to understand recent DHF outbreaks caused by DENV-3 infections in the Indian subcontinent and Latin America.

Previous phylogenetic analysis of DENV-3 has principally relied on complete or partial sequences of the pre-M/M and E genes (13,16,17,23). Our analysis used a 708-base segment, positions 437-1,145, spanning pre-M/M and a portion of the E gene, coding for 236 amino acids. This region was selected because it both conserved the original phylogenetic relationship identified by Lanciotti et al. and, in preliminary analysis with previously established sequences, captured 44% of the variable sites within DENV-3, subtype III Sri Lankan sequences. No insertion/deletion mutations and no hypervariable regions were detected in this span.

A total of 40 DENV-3 sequences, including 21 sequences available from GenBank and 19 newly determined Indian subcontinent and African sequences (Table 1) were compared. Dates of isolation ranged from 1963 to 1998. With the exception of 63PuertoR, all sequences were from low-passage (<4) virus cultures. Several approaches to phylogenetic analysis, including maximum likelihood, parsimony, and distance methods, were compared. All approaches yielded identical or nearly identical topologies. Results presented here used the Tamura-Nei algorithm to calculate genetic distances and the minimum evolution method to create the trees (Figure 1). The tree identifies four distinct lineages that correspond to the region of isolation, reproducing the same evolutionary relationship first described by Lanciotti, et al. (23). Subtype I includes isolates from Southeast Asia and the South Pacific islands; subtype II consists of isolates from Thailand; subtype III is comprised of isolates from the Indian subcontinent, East Africa, and a single isolate from Samoa; and subtype IV includes Puerto Rico and Tahiti. Similarity within subtypes was high, with subtype III showing the greatest mean similarity (98.4%), followed by subtypes I, II, and IV (Table 4).

All 24 Sri Lankan, Indian, and East African strains fell into subtype III (Figure 1). The circulating virus genotypes within this region have remained closely related over the relatively long period of 18 years (1981-1998), indicating that countries bordering the western Indian Ocean form a geographically distinct region with regard to DENV-3 viruses. DENV-2 viruses in the regions also form a subtype with a similar geographic distribution (26,27). Frequent trade between East Africa, Western Indian Ocean islands, and the Indian subcontinent may have been responsible for the movement of dengue viruses throughout the region (26,28). Rico-Hesse, for example, demonstrated the introduction of DENV-2 to Africa from islands in the Indian Ocean (26). The earliest subtype III virus on record is an isolate from India in 1966; this virus occupies a node that is ancestral to all the subsequent Asian and African isolates (R.S. Lanciotti, pers. comm.), suggesting that the DENV-3, subtype III viruses have their origin in the Indian subcontinent and have subsequently spread out of the region.

In Sri Lanka, regular epidemics of DHF have been observed only after 1988. DENV-3 is responsible for many of the infections that progress to DHF (20,22). DENV-3 isolates obtained before and after the emergence of DHF are very closely related and belong to subtype III, indicating that the emergence of DHF on the island is not due to the introduction of a new subtype from outside the region. However, within subtype III, most Sri Lankan isolates (except for 93SriLan1) from before and after the emergence of DHF segregated into two distinct clades, designated groups A and B (Figure 1). Group A, with nine isolates from 1981 to 1989, consists of viruses collected up to the year epidemic DHF emerged in Sri Lanka but contains no isolates from later than 1989. Group B includes eight isolates from 1989 to 1998 but none from before 1989. Temporally, the two groups are continuous, by virtue of sharing isolates in 1989. Group A includes isolate 89SriLan2, while group B contains 89SriLan1, 89SriLan3, and 89SriLanA. However, the groups do not form a continuous lineage; they share a common ancestor only at the node for subtype III (Figure 1). Group B shares ancestral nodes with isolates from India and East Africa. Because the Indian and East African isolates overlap temporally with group A (all isolates are from the 1980s), group A and group B lineages likely diverged sometime before 1981 and followed distinct evolutionary pathways.

We propose two likely scenarios that led to the emergence of group B viruses in Sri Lanka. One possibility is that the group B viruses were introduced from India or East Africa into Sri Lanka because Indian and East African isolates from the mid-1980s are closely related to Sri Lankan group B viruses (Figure 1). Of the two regions, India is the more likely source because of geographic proximity to Sri Lanka, although the East African viruses could be the direct ancestors of the group B viruses. Another possibility is that both groups co-circulated in Sri Lanka in the early 1980s, with group B being a minor population. Some selective force operating in the late 1980s may have shifted the balance in favor of group B viruses. In either case, group B viruses emerged in Sri Lanka because a subtype III variant already established in the greater region became more common in Sri Lanka and not because a novel virus evolved and emerged de novo on the island.

DENV-3, subtype III was detected in the Americas during DF and DHF outbreaks in Nicaragua and Panama in 1994 (11). Subsequently, the virus has spread to many countries in Latin America, and DENV-3-associated DHF was confirmed in several countries (13,14,16,29-31) (Figure 3). To establish the relationship of recent Latin American DENV-3 isolates to each other and to the previously identified Indian subcontinent and East African subtype III isolates, we sequenced and analyzed a 966-base region spanning positions 179-1,145 on the viral genome, capturing a portion of the C gene, all of pre-M/M gene, and a portion of the E gene. This region adds 288 positions to the 5' end of the sequences initially presented in this study.

[FIGURE 3 OMITTED]

Forty-three isolates were sequenced (21 from Mexico and Central and South America, 16 from Sri Lanka, 1 from India, and 5 from East Africa) (Table 1). The D3H-87 belonging to DENV-3, subtype 1 was used as an outgroup. Years of isolation ranged from 1983 to 2001, an 18-year span. Except for the D3H-87 outgroup, all other isolates were low-passage clinical isolates. Most of the nucleotide mutations were silent: only 12 amino acid positions showed any variability and only 2 positions showed variability in more than one isolate. Consequently, the evolutionary relationships observed in this analysis likely reflect the results of genetic drift and are unlikely to have been influenced by host-specific selection events on this portion of the genome (26).

All sequences included in this analysis fell within subtype III (data not shown). Several approaches to phylogenetic analysis were compared, and all approaches yielded identical or nearly identical topologies. We used the Tamura-Nei algorithm to calculate genetic distances and the minimum evolution method to create the trees (Figure 2). Bootstrap values are shown at critical nodes. Despite the high overall similarity of the isolates in this analysis, geographically and temporally distinct groups formed separate lineages. Generally, two separate lineages formed within subtype III. The first consists of group A viruses isolated from 1981 to 1989 in Sri Lanka. These viruses have been associated only with DF. The second is composed of Sri Lankan group B, Indian, East African, and all of the isolates from Mexico and Central and South America.

Within group A, members are closely related, with a nucleotide mean similarity of 99.4% (Table 5). Within the expanded group B and related viruses, three distinct clades exist: a group of closely related Sri Lankan isolates from 1989 to 1998, 5 East African isolates from 1985 to 1993, and 21 isolates from 1994 to 2001 from Latin America. Isolates 84India and 93SriLan1 are less closely related to the other geographically distinct isolates in the larger second lineage.

The isolates from Latin America all emerge from a common node on the tree, suggesting a single introduction of a virus and the subsequent diversification of the virus population from the founding strain. The DENV-3, subtype III isolates from Nicaragua, Panama, and Costa Rica are closest to the Latin American group's originating node, with the more recent isolates found farther from that node, reflecting the viral population's ongoing evolution after the point source introduction.

The internal branch from the Latin American group shares a common node with the isolates from East Africa. The common hypothetical ancestor for Latin America and East Africa then shares a common node with the Sri Lankan group B virus isolates. Both on the phylogenetic tree and in pair-wise comparisons (Table 5), the Latin American group was more closely related to the isolates from East Africa than to the group B Sri Lankan isolates. Furthermore, the East African isolates pre-date the earliest Latin American isolates by 9 years, while the less closely related Sri Lankan group B isolates are nearly contemporaneous with the Latin American isolates. Therefore, the point source DENV-3 introduction into Latin America is most likely to have its origins in East Africa and not the Indian subcontinent (Figure 4).

[FIGURE 4 OMITTED]

Little is known about dengue activity in Africa, particularly DENV-3 (32). DENV-3 was first detected on the African continent in 1984 to 1985 during an outbreak in Mozambique (32). Later studies of U.S. troops in Africa and the Persian Gulf suggested that DENV-3 is endemic in those regions but largely undetected (33). Our results show that all East African DENV-3 isolates belong to subtype III. The fact that DENV-3 was only first isolated from East Africa in 1985, whereas the viruses were present in the Indian subcontinent at least as far back as 1966 (R.S. Lanciotti, pers. comm.), suggests that DENV-3, subtype III was introduced from the Indian subcontinent into East Africa in or before 1984 (Figure 4). This introduction led to the establishment of a stable East African group of DENV-3, subtype III because all the isolates from Mozambique, Kenya, and Somalia isolated from 1985 to 1993 form a distinct clade within subtype III (Figure 3).

The DENV-3, subtype III viruses introduced into Latin America are most closely related to subtype III viruses in East Africa (Figure 3). Although we can only speculate about the exact mode of transport of DENV-3 into Latin America, we propose that Panama, with its canal that attracts goods as well as civilians and military personnel from other parts of the world, may have been the point of introduction of subtype III into the Americas. Similarly, the introduction of DENV-2 in 1981 into Cuba may be attributable to Cuban military personnel traveling between Southeast Asia and Cuba (24,34).

Epidemiologic and clinical studies on dengue in Indonesia in the 1970s pointed to strain differences between DENV-3 viruses contributing to transmission and disease severity (35,36). Despite their overall similarity at the nucleotide level, the DENV-3, subtype III isolates examined in this study have been associated with severe or mild disease outbreaks (Figure 3). Sri Lankan group A viruses were isolated during a time of little to no DHF, while group B viruses were isolated after the emergence of DHF in Sri Lanka. The emergence of DHF in Sri Lanka was not accompanied by a change in dengue transmission or the abundance of any particular serotype (21). Implicating DENV-3 directly as the cause of DHF in Sri Lanka has been difficult because few virus isolates are available from DHF patients in Sri Lanka. However, during dengue surveillance studies in 1997, only DENV-3 was isolated from hospitalized dengue cases, whereas DENV-1, DENV-2, and DENV-3 were isolated from patients visiting outpatient clinics (22). These observations suggest that DENV-3 is responsible for severe dengue disease in Sri Lanka. Further studies are required to better establish the relative contribution of DENV-3 to severe disease in Sri Lanka.

The current studies support a viral genetic basis for severe and mild disease outbreaks. We found that the population of DENV-3 viruses associated with DHF in Sri Lanka did not appear to be direct descendants of the group A viruses that were circulating before DHF emerged in that country. The Sri Lankan 1989-1997 isolates are more closely related to the isolates from East Africa and the isolates from the Americas than they are to the isolates from 1981 to 1989 in Sri Lanka (Figure 3). All three groups of subtype III viruses (Sri Lankan group B, East African group, and Latin American group) associated with DHF are more closely related to each other than they are to the pre-DHF group A viruses from Sri Lanka (Figure 3). Thus, all the viruses within subtype III are closely related (mean 98.4% identity at the nucleotide level), yet they form distinct phylogenetic groups associated with mild or severe disease.

The Sri Lankan group B viruses may be associated with severe disease unlike group A viruses because the group B viruses are inherently more virulent. Alternatively, the ability of preexisting dengue antibody to neutralize group A viruses and enhance group B viruses may account for the observed associations with severe and mild disease. In a recent study, antibodies against American DENV-1 viruses neutralized the Native American DENV-2 genotype better than the Southeast Asian DENV-2 genotype that is currently circulating in the Americas and causing DHF (37). This study lends support to the idea that Asian DENV-2 may produce a more severe disease not because of inherent virulence properties but because persons with previous primary DENV-1 infections may enhance infection with this genotype and neutralize infections with the Native American DENV-2 genotype. Similarly, DENV-2 and -3 are the common serotypes in Sri Lanka, and persons with previous primary DENV-2 infections could neutralize the DENV-3 group A viruses better than the group B viruses. This difference may explain the unexpected emergence of DHF associated with group B. Further comparative studies with group A and B viruses are needed to understand their association with mild and severe disease, respectively.

Table 1. Dengue virus type 3 sequences used (a)

Strain       Y        Location      Name                   Subtype

D1266       1983      Sri Lanka     83SriLan1              III
D1306       1983      Sri Lanka     83SriLan2              III
D1307       1983      Sri Lanka     83SriLan3              III
D1336       1983      Sri Lanka     83SriLan4              III
D1440       1984      Sri Lanka     84SriLan1              III
073         1985      Sri Lanka     85SriLan               III
D2783       1989      Sri Lanka     89SriLan1              III
D2863       1989      Sri Lanka     89SriLan2              III
D2803       1989      Sri Lanka     89SriLan3              III
D3197       1990      Sri Lanka     90SriLan1              III
D5231       1993      Sri Lanka     93SriLan1              III
D9397       1994      Sri Lanka     94SriLan1              III
L57         1997      Sri Lanka     97SriLan1              III
K1          1998      Sri Lanka     98SriLan               III
1557        1985      Mozambique    85Mozamb1              III
1558        1985      Mozambique    85Mozamb2              III
1559        1985      Mozambique    85Mozamb3              III
251991      1991      Kenya         91Kenya                III
SOM079      1993      Somalia       93Somalia              III
32267       1994      Nicaragua     94Nicara1              III
6845        1998      Nicaragua     98Nicara1              III
7431        1998      Nicaragua     98Nicara2              III
7071        1998      Nicaragua     98Nicara3              III
BC 96/94    1994      Panama        94Panama1              III
032231      1994      Panama        94Panama2              III
BC 13/96    1994      Panama        94Panama3              III
BC 20/97    1996      Mexico        96Mexico1              III
BC 172/97   1996      Mexico        96Mexico2              III
BC 184/97   1996      Mexico        96Mexico3              III
BC 173/97   1996      Mexico        96Mexico4              III
17605       1995      Costa Rica    95CostaR1              III
17608       1995      Costa Rica    95CostaR2              III
322473      1995      Costa Rica    95CostaR3              III
322488      1995      Costa Rica    95CostaR4              III
20/8        1997      Guatemala     97Guatem1              III
366-781     1998      Puerto Rico   98PuertorR1            III
400-996     2000      Puerto Rico   00PuertoR1             III
MK          1998      El Salvador   98ElSalv1              III
612210      2001      Venezuela     01Venezue1             III
VEN03       2001      Venezuela     01VEN03                III
Ref. 18     1981-91   Sri Lanka     81,85,89,91 SriLanA    III
Ref. 18     1984      India         84IndiaA               III
Ref. 18     1986      Samoa         86Samoa                III
Ref. 18     1962-86   Thailand      62,73,86,86 Thailand   II
Ref. 18     1983      Philippines   83Philipp              I
Ref. 18     1989      Tahiti        89Tahiti               I
Ref. 18     1992      Fiji          92Fiji                 I
Ref. 18     1973-85   Indonesia     73,78,85 Indones       I
Ref. 18     1974-81   Malaysia      74,81 Malaysi          I
Ref. 18     1956      Philippines   D3H-87                 I
Ref. 18     1963-77   Puerto Rico   63,77 PuertoR          IV
Ref. 18     1965      Tahiti        65 Tahiti              IV

Strain      Sequence source   GenBank accession no.

D1266       This study        AF547225
D1306       This study        AF547226
D1307       This study        AF547227
D1336       This study        AF547228
D1440       This study        AF547229
073         This study        AF547241
D2783       This study        AF547230
D2863       This study        AF547231
D2803       This study        AF547232
D3197       This study        AF547233
D5231       This study        AF547234
D9397       This study        AF547235
L57         This study        AF547242
K1          This study        AF547243
1557        This study        AF547236
1558        This study        AF547237
1559        This study        AF547238
251991      This study        AF547239
SOM079      This study        AF547240
32267       This study        AF547244
6845        This study        AF547245
7431        This study        AF547246
7071        This study        AF547262
BC 96/94    This study        AF547247
032231      This study        AF547248
BC 13/96    This study        AF547249
BC 20/97    This study        AF547250
BC 172/97   This study        AF547251
BC 184/97   This study        AF547252
BC 173/97   This study        AF547253
17605       This study        AF547254
17608       This study        AF547255
322473      This study        AF547256
322488      This study        AF547257
20/8        This study        AF547263
366-781     This study        AF547258
400-996     This study        AF547264
MK          This study        AF547259
612210      This study        AF547260
VEN03       This study        AF547261
Ref. 18     GenBank           L11431,L11436-L11438
Ref. 18     GenBank           L11424
Ref. 18     GenBank           L11435
Ref. 18     GenBank           L11440-L11442,L11620
Ref. 18     GenBank           L11432
Ref. 18     GenBank           L111619
Ref  18     GenBank           L11422
Ref. 18     GenBank           L11425,L11426,L11428
Ref. 18     GenBank           L11429,L11427
Ref. 18     GenBank           L11423
Ref. 18     GenBank           L11433,L11434
Ref. 18     GenBank           L11439

(a) Includes original identifier for strain, year of isolation, taxa
name used in this paper, dengue virus 3 subtype, source of viral
sequence, and GenBank accession numbers.

Table 2. Amino acid substitutions conserved within each dengue virus 3
subtype for the isolates used to create the phylogenetic tree in
Figure 1 (a)
                                        Position

Subtype    Name        31   55   57   128   135   148   188   234

Outgroup   56Philipp    I    H    T     L     I     L     D     I
I          73Indones   --    L   --     F    --    --    --     V
I          74Malaysi   --    L   --     F    --    --    --     V
I          78Indones   --    L   --     F    --    --    --     V
I          81Malaysi   --    L   --     F    --    --    --     V
I          83Philipp   --    L   --     F    --    --    --     V
I          85Indones   --    L   --     F    --    --    --     V
I          89Tahiti    --    L   --     F    --    --    --     V
I          92Fiji      --    L   --     F    --    --    --     V
II         62Thailan   --   --    A    --    --     W    --    --
II         73Thailan   --   --    A    --    --    --    --    --
II         86Thailan   --    L    A    --    --    --    --    --
II         87Thailan   --    L    A    --    --    --    --    --
III        85Mozamb1   --   --   --    --    --    --    --    --
III        85Mozamb2   --   --   --    --    --    --    --    --
III        85Mozamb3   --   --   --    --    --    --    --    --
III        84IndiaA    --   --   --    --    --    --    --    --
III        91 Kenya    --   --   --    --    --    --    --    --
III        93Somolio   --   --   --    --    --    --    --    --
III        81SriLanA   --   --   --    --    --    --    --    --
III        83SriLan1   --   --   --    --    --    --    --    --
III        83SriLan2   --   --   --    --    --    --    --    --
III        83SriLan3   --   --   --    --    --    --    --    --
III        83SriLan4   --   --   --    --    --    --    --    --
III        84SriLan1   --   --   --    --    --    --    --    --
III        85SriLanA   --   --   --    --    --    --    --    --
III        85SriLan    --   --   --    --    --    --    --    --
III        89SriLan2   --   --   --    --    --    --    --    --
III        89SriLanA   --   --   --    --    --    --    --    --
III        89SriLan1   --   --   --    --    --    --    --    --
III        89SriLan3   --   --   --    --    --    --    --    --
III        90SriLan1   --   --   --    --    --    --    --    --
III        91SriLanA   --   --   --    --    --    --    --    --
III        93SriLan1   --   --   --    --    --    --    --    --
III        94SriLan1   --   --   --    --    --    --    --    --
III        97SriLan1   --   --   --    --    --    --    --    --
III        98SriLan1   --   --   --    --    --    --    --    --
III        86Samoa     --   --   --    --    --    --    --    --
IV         63PuertoR    T   --   --     F     L     M     E    --
IV         65Tahiti     T   --   --     F     L     M     E    --
IV         77PuertoR    T   --   --     F     L     M     E    --

(a) Reference strain 56Philipp is the highly passaged laboratory strain
H87. Positions are numbered sequentially from the first position in the
pre-M/M protein

Table 3. Nucleotide substitutions conserved within dengue virus 3,
subtype III groups

                                           Position

Group          Strain       338   429   503   566   653   686   695

Reference      H87            A     A     C     G     G     T     T
Group A        83SriLan1     --     G     T     A    --    --     C
Group A        83SriLan2     --     G     T     A    --    --     C
Group A        83SriLan3     --     G     T     A    --    --     C
Group A        83SriLan4     --     G     T     A    --    --     C
Group A        84SriLan1     --     G    --    --    --    --    --
Group A        85SriLan1     --     G     T     A    --    --     C
Group A        89SriLan2     --     G     T     A    --    --     C
unclassified   84India        G     G    --    --     A    --    --
unclassified   93SriLan1      G     G    --    --     A    --    --
East Africa    85Mozamb1      G    --    --    --     A    --    --
East Africa    85Mozamb2      G    --    --    --     A    --    --
East Africa    85Mozamb3      G    --    --    --     A    --    --
East Africa    91Kenya        G    --    --    --     A    --    --
East Africa    93Somalia      G    --    --    --     A    --    --
Group B        89SriLan1      G     G    --    --     A    --    --
Group B        89SriLan3     --     G    --    --     A    --    --
Group B        90SriLan1      G     G    --    --     A    --    --
Group B        94SriLan1     --     G    --    --     A    --    --
Group B        97SriLan1      G     G    --    --     A    --    --
Group B        98SriLan1      G     G    --    --     A    --    --
L. America     94Nicara1      G     G    --    --     A     C    --
L. America     94Panama1      G     G    --    --     A     C    --
L. America     94Panama2      G     G    --    --     A     C    --
L. America     94Panama3      G     G    --    --     A     C    --
L. America     CostaRica1     G     G    --    --     A     C    --
L. America     CostaRica2     G     G    --    --     A     C    --
L. America     CostaRica3     G     G    --    --     A     C    --
L. America     CostaRica4     G     G    --    --     A     C    --
L. America     96Mexico1      G     G    --    --     A     C    --
L. America     96Mexico2      G     G    --    --     A     C    --
L. America     96Mexico3      G     G    --    --     A     C    --
L. America     96Mexico4      G     G    --    --     A     C    --
L. America     97Guatem1      G     G    --    --     A     C    --
L. America     98Nicara1      G     G    --    --     A     C    --
L. America     98Nicara2      G     G    --    --     A     C    --
L. America     98Nicara3      G     G    --    --     A     C    --
L. America     98ElSalv1      G     G    --    --     A     C    --
L. America     01Venezue1     G     G    --    --     A     C    --
L. America     01VEN03        G     G    --    --     A     C    --
L. America     98PuertoR1     G     G    --    --     A     C    --
L. America     00PuertoR1     G     G    --    --     A     C    --

Group          Strain       707   728   734   749   791   821   866

Reference      H87            G     C     C     C     A     T     C
Group A        83SriLan1     --    --     T     T    --    --     T
Group A        83SriLan2     --    --     T     T    --    --     T
Group A        83SriLan3     --    --    --     T    --    --     T
Group A        83SriLan4     --    --     T     T    --    --     T
Group A        84SriLan1     --    --    --     T    --    --    --
Group A        85SriLan1     --    --     T     T    --    --     T
Group A        89SriLan2     --    --     T     T    --    --     T
unclassified   84India       --    --    --     T    --    --    --
unclassified   93SriLan1     --    --    --    --    --    --    --
East Africa    85Mozamb1     --    --    --    --    --    --    --
East Africa    85Mozamb2     --    --    --    --    --    --    --
East Africa    85Mozamb3     --    --    --    --    --    --    --
East Africa    91Kenya       --    --    --    --    --    --    --
East Africa    93Somalia     --    --    --    --    --    --    --
Group B        89SriLan1      A     T    --    --    --     C    --
Group B        89SriLan3      A     T    --    --    --     C    --
Group B        90SriLan1      A     T    --    --    --     C    --
Group B        94SriLan1      A     T    --    --    --     C    --
Group B        97SriLan1      A     T     T    --    --     C    --
Group B        98SriLan1      A     T    --    --    --     C    --
L. America     94Nicara1     --    --    --    --     G     C    --
L. America     94Panama1     --    --    --    --     G     C    --
L. America     94Panama2     --    --    --    --     G     C    --
L. America     94Panama3     --    --    --    --     G     C    --
L. America     CostaRica1    --    --    --    --     G     C    --
L. America     CostaRica2    --    --    --    --     G     C    --
L. America     CostaRica3    --    --    --    --     G     C    --
L. America     CostaRica4    --    --    --    --     G     C    --
L. America     96Mexico1     --    --    --    --     G     C    --
L. America     96Mexico2     --    --    --    --     G     C    --
L. America     96Mexico3     --    --    --    --     G     C    --
L. America     96Mexico4     --    --    --    --     G     C    --
L. America     97Guatem1     --    --    --    --     G     C    --
L. America     98Nicara1     --    --    --    --     G     C    --
L. America     98Nicara2     --    --    --    --     G     C    --
L. America     98Nicara3     --    --    --    --     G     C    --
L. America     98ElSalv1     --    --    --    --     G     C    --
L. America     01Venezue1    --    --    --    --     G     C    --
L. America     01VEN03       --    --    --    --     G     C    --
L. America     98PuertoR1    --    --    --    --     G     C    --
L. America     00PuertoR1    --    --    --    --     G     C    --

Group          Strain       896   902   978   1010   1019   1056

Reference      H87            T     T     C     C      T      C
Group A        83SriLan1      C     C     T     T     --     --
Group A        83SriLan2      C     C     T     T     --     --
Group A        83SriLan3     --     C    --     T     --     --
Group A        83SriLan4      C     C     T     T     --     --
Group A        84SriLan1      C    --    --    --     --     --
Group A        85SriLan1      C     C     T     T     --     --
Group A        89SriLan2      C     C     T     T     --     --
unclassified   84India       --    --    --    --      G     --
unclassified   93SriLan1     --    --    --    --      G     --
East Africa    85Mozamb1     --    --    --    --      G     --
East Africa    85Mozamb2     --    --    --    --      G     --
East Africa    85Mozamb3     --    --    --    --      G     --
East Africa    91Kenya       --    --    --    --      G     --
East Africa    93Somalia     --    --    --    --      G     --
Group B        89SriLan1     --    --    --    --      G     --
Group B        89SriLan3     --    --    --    --      G     --
Group B        90SriLan1     --     C    --    --      G     --
Group B        94SriLan1     --    --    --    --      G     --
Group B        97SriLan1      C    --    --    --      G     --
Group B        98SriLan1      C    --    --    --      G     --
L. America     94Nicara1     --    --    --    --      G      T
L. America     94Panama1     --    --    --    --      G      T
L. America     94Panama2     --    --    --    --      G      T
L. America     94Panama3     --    --    --    --      G      T
L. America     CostaRica1    --    --    --    --      G      T
L. America     CostaRica2    --    --    --    --      G      T
L. America     CostaRica3    --    --    --    --      G      T
L. America     CostaRica4    --    --    --    --      G      T
L. America     96Mexico1     --    --    --    --      G      T
L. America     96Mexico2     --    --    --    --      G      T
L. America     96Mexico3     --    --    --    --      G      T
L. America     96Mexico4     --    --    --    --      G      T
L. America     97Guatem1     --    --    --    --      G      T
L. America     98Nicara1     --    --    --    --      G      T
L. America     98Nicara2     --    --    --    --      G      T
L. America     98Nicara3     --    --    --    --      G      T
L. America     98ElSalv1     --    --    --    --      G      T
L. America     01Venezue1    --    --    --    --      G      T
L. America     01VEN03       --    --    --    --      G      T
L. America     98PuertoR1    --    --    --    --      G      T
L. America     00PuertoR1    --    --    --    --      G      T

(a) Reference strain is the highly passaged laboratory strain H87.
Positions are numbered sequentially from the first nucleotide position
at the 5' end of the genome.

Table 4. Summary of within- and between-subtype nucleotide
mean similarity for dengue virus 3 isolates shown in Figure 1 (a)

                             Between subtype
                              similarity (%)
          Within subtype
Subtype   similarity (%)    I      II     III

I              98.1
II             97.7        94.9
III            98.4        95.6   96.3
IV             97.6        92.3   92.5   92.7

(a) Mean similarities were calculated with the Tamura-Nei distance
algorithm.

Table 5. Summary of within- and between-group nucleotide mean
similarity for the dengue virus 3, subtype III virus isolates shown in
Figure 2 (a)
                                     Between group similarity
                Within-group
Subgroup         similarity    Subgroup A   East Africa   Subgroup B

Subgroup A         99.4%
East Africa        99.5%          98.2%
Subgroup B         98.8%          97.9%        98.7%
Latin America      99.5%          98.0%        99.0%         98.5%

(a) Mean similarities were calculated with the Tamura-Nei distance
algorithm.

Acknowledgments

We thank Irene Bosch, Vance Vorndam, Niranjan Kanesathasan, and Eric Wagar for virus isolates and Tissa Vitarana, Gaya Colombage, Nalini Withana, and staff members in the Virology Department at the Medical Research Institute in Colombo for their help.

This work was supported by a Junior Faculty Development Award and other funds from the University of North Carolina at Chapel Hill (A.M.S.)

References

(1.) Briese T, Jia XY, Huang C, Grady LJ, Lipkin WI. Identification of a Kunjin/West Nile-like flavivirus in brains of patients with New York encephalitis. Lancet 1999;354:1261-2.

(2.) Lanciotti RS, Roehrig JT, Deubel V, Smith J, Parker M, Steele K, et al. Origin of the West Nile virus responsible for an outbreak of encephalitis in the northeastern United States. Science 1999:286:2333-7.

(3.) Morse SS. Emerging viruses. New York: Oxford University Press: 1993.

(4.) Halstead SB. Epidemiology of dengue and dengue hemorrhagic fever. In: Kuno DJGaG, editor. Dengue and dengue hemorrhagic fever. New York: CAB International; 1997. p. 23-44.

(5.) Gubler DJ. Dengue and dengue hemorrhagic lever. Clin Microbiol Rev 1998;11:480-96.

(6.) Rosen L. The emperor's new clothes revisited, or reflections on the pathogenesis of dengue hemorrhagic fever. Am J Trop Med Hyg 1977;26:337-43.

(7.) Gubler DJ, Reed D, Rosen L, Hitchcock JR Jr. Epidemiologic, clinical, and virologic observations on dengue in the Kingdom of Tonga. Am J Trop Mod Hyg 19787;27:581-9.

(8.) Gubler DJ, Suharyono W, Lubis I, Eram S, Gunarso S. Epidemic dengue 3 in central Java, associated with low viremia in man. Am J Trop Med Hyg 1981;30:1094-9.

(9.) Guzman MG, Kouri GP, Bravo J, Soler M, Velasquez S, Morier L. Dengue haemorrhagic fever in Cuba. A retrospective seroepidemiologic study. Am J Trop Med Hyg 1990;42:179-84.

(10.) Rico-Hesse R, Harrison LM, Salas RA, Tovar D, Nisalak A, Ramos C, et al. Origins of dengue type 2 viruses associated with increased pathogenicity in the Americas. Virology 1997;230:244-51.

(11.) Dengue type 3 infection. Nicaragua and Panama, October-November 1994. Wkly Epidemiol Rec 1995;70:41-3.

(12.) Gubler DJ. The global pandemic of dengue/dengue haemorrhagic fever: current status and prospects for the future. Ann Acad Med Singapore 1998;27:227-34.

(13.) Guzman M, Huelva G, Saenz E, Quiroz E, De los Reyes J, Balmaseda A. Reintroduccion del dengue 3 en las Americas: 1994-1996. Archivos Venezolanos de Medicina Tropial 1998;2:8-19.

(14.) Balmaseda A, Sandoval E, Perez L, Gutierrez CM, Harris E. Application of molecular typing techniques m the 1998 dengue epidemic in Nicaragua. Am J Trop Med Hyg 1999;61:893-7.

(15.) Nogueira RM, Miagostovich MP, de Filippis AM, Pereira MA, Schatzmayr HG. Dengue virus type 3 in Rio de Janeiro, Brazil. Mem Inst Oswaldo Cruz 2001;96:925-6.

(16.) Usuku S, Castillo L, Sugimoto C, Noguchi Y, Yogo Y, Kobayashi N. Phylogenetic analysis of dengue-3 viruses prevalent in Guatemala during 1996-1998. Arch Virol 2001;146:1381-90.

(17.) Chungue E, Deubel V, Cassar O, Laille M, Martin PM. Molecular epidemiology of dengue 3 viruses and genetic relatedness among dengue 3 strains isolated from patients with mild or severe form of dengue fever in French Polynesia. J Gen Virol 1993;74:2765-70.

(18.) Leitmeyer KC, Vaughn DW, Watts DM, Salas RA, de Chacon IV, Ramos C, et al. Dengue virus structural differences that correlate with pathogencsis. J Virol 1999;73:4738-47.

(19.) Srivastava VK, Suri S, Bhasin A, Srivastava L, Bharadwaj M. An epidemic of dengue haemorrhagic fever and dengue shock syndrome in Delhi: a clinical study. Ann Trop Paediatr 1990;10:329-34.

(20.) Vitarana UT, Jayasekera N, Withane N, Gubler DJ. Finding the cause of dengue hemorrhagic fever outbreaks in Sri Lanka. Arbovirus Research in Australia 1993:6:125-9.

(21.) Messer WB, Sivananthan K, Elvitigala J, Preethimala LD, Ramesh R, Withana N, et al. Epidemiology of dengue in Sri Lanka before and after the emergence of epidemic dengue hemorrhagic fever. Am J Trop Med Hyg 2002;66:765-73.

(22.) de Silva AM, Sivananthan K, Withana N, Vorndam V, Gubler DJ. Dengue 3 virus is responsible for recent epidemics of dengue hemorrhagic fever in Sri Lanka. In: Annual Meeting of the American Society for Tropical Medicine and Hygiene. San Juan, Puerto Rico; 1998.

(23.) Lanciotti RS, Lewis JG, Gubler DJ, Trent DW. Molecular evolution and epidemiology of dengue-3 viruses. J Gen Virol 1994:75:65-75.

(24.) Rico-Hesse R, Harrison LM, Salas RA, Tovar D, Nisalak A, Ramos C, et al. Origins of dengue type 2 viruses associated with increased pathogenicity in the Americas. Virology 1997;230:244-51.

(25.) Kobayashi N, Thayan R, Sugimoto C, Oda K, Saat Z, Vijayamalar B, et al. Type-3 dengue viruses responsible for the dengue epidemic in Malaysia during 1993-1994. Am J Trop Med Hyg 1999;60:904-9.

(26.) Rico-Hesse R. Molecular evolution and distribution of dengue viruses type 1 and 2 in nature. Virology 1990;174:479-93.

(27.) Lewis JA, Chang GJ, Lanciotti RS, Kinney RM, Mayer LW, Trent DW. Phylogenetic relationships of dengue-2 viruses. Virology 1993;197:216-24.

(28.) Wang E, Ni H, Xu R, Barrett AD, Watowich SJ, Gubler DJ, et al. Evolutionary relationships of endemic/epidemic and sylvatic dengue viruses. J Virol 2000;74:3227-34.

(29.) Gubler D. Dengue and dengue hemorrhagic fever: its history and resurgence as a global public health problem. In: Kuno DGaG, editor. Dengue and dengue hemorrhagic fever. New York: CAB International; 1997.

(30.) Isturiz RE, Gubler DJ, Brea del Castillo J. Dengue and dengue hemorrhagic fever in Latin America and the Caribbean. Inflect Dis Clin North Am 2000; 14:121-40.

(31.) Briseno-Garcia B, Gomez-Dantes H, Argott-Ramirez E, Montesano R, Vazquez-Martinez AL, Ibanez-Bernal S, et al. Potential risk for dengue hemorrhagic fever: the isolation of serotype dengue-3 in Mexico. Emerg Infect Dis 1996;2:133-5.

(32.) Gubler DJ, Sather GE, Kuno G, Cabral JR. Dengue 3 virus transmission in Africa. Am J Trop Med Hyg 1986;35:1280-4.

(33.) Kanesa-thasan N, Chang GJ, Smoak BL, Magill A, Burrous MJ, Hoke CH Jr. Molecular and epidemiologic analysis of dengue virus isolates from Somalia. Emerg Infect Dis 1998;4:299-303.

(34.) Gubler DJ, Trent DW. Emergence of epidemic dengue/dengue hemorrhagic fever as a public health problem in the Americas. Infect Agents Dis 1993;2:383-93.

(35.) Gubler DJ, Suharyono W, Lubis I, Eram S, Sulianti Saroso J. Epidemic dengue hemorrhagic lever in rural Indonesia. I. Virological and epidemiological studies. Am J Trop Med Hyg 1979;28:701-10.

(36.) Lee E, Gubler DJ, Weir RC, Dalgamo L. Genetic and biological differentiation of dengue 3 isolates obtained from clinical cases in Java, Indonesia, 1976-1978. Arch Virol 1993;133:113-25.

(37.) Kochel TJ, Watts DM, Halstead SB, Hayes CG, Espinoza A, Felices V. Effect of dengue-1 antibodies on American dengue-2 viral infection and dengue haemorrhagic tbver. Lancet 2002;360:310-2.

Address for correspondence: Aravinda de Silva, Department of Microbiology and Immunology, CB#7290, University of North Carolina, Chapel Hill, NC 27599, USA; fax: [phone]; email: [e-mail]

William B. Messer, * Duane J. Gubler, ([dagger]) Eva Harris, ([double dagger]) Kamalanayani Sivananthan, ([section]) and Aravinda M. de Silva *

* University of North Carolina, Chapel Hill, North Carolina, USA; ([dagger]) Centers for Disease Control and Prevention, Fort Collins, Colorado, USA; ([double dagger]) University of California, Berkeley, California, USA; and ([section]) Medical Research Institute, Colombo, Sri Lanka

Dr. Messer earned his Ph.D. in ecology from the University of North Carolina at Chapel Hill. His broad research interests include the emergence, movement, and evolution of human pathogens. For his doctoral research, he evaluated host, viral, and environmental factors that may have contributed to the sudden emergence of dengue hemorrhagic lever in Sri Lanka. He attends medical school at the University of North Carolina.