Volume 18, Number 11—November 2012
Phylogeography of Dengue Virus Serotype 4, Brazil, 2010–2011
Dengue virus serotype 4 (DENV-4) reemerged in Roraima State, Brazil, 28 years after it was last detected in the country in 1982. To study the origin and evolution of this reemergence, full-length sequences were obtained for 16 DENV-4 isolates from northern (Roraima, Amazonas, Pará States) and northeastern (Bahia State) Brazil during the 2010 and 2011 dengue virus seasons and for an isolate from the 1982 epidemic in Roraima. Spatiotemporal dynamics of DENV-4 introductions in Brazil were applied to envelope genes and full genomes by using Bayesian phylogeographic analyses. An introduction of genotype I into Brazil from Southeast Asia was confirmed, and full genome phylogeographic analyses revealed multiple introductions of DENV-4 genotype II in Brazil, providing evidence for >3 introductions of this genotype within the last decade: 2 from Venezuela to Roraima and 1 from Colombia to Amazonas. The phylogeographic analysis of full genome data has demonstrated the origins of DENV-4 throughout Brazil.
Dengue virus (DENV), a widespread arthropod-borne virus that commonly affects humans, belongs to the family Flaviviridae, genus Flavivirus, and is classified into 4 distinct serotypes (DENV 1–4). DENV is most prevalent in tropical and subtropical areas, where eco-epidemiologic conditions appear to sustain the virus. In particular, these regions harbor 2 competent DENV vectors (Aedes aegypti and A. albopictus mosquitoes) and have environments favorable for DENV (1). In past decades, the number of countries reporting DENV cases and those with endemic DENV has increased dramatically. These increases reflect the expanding habitat of the Aedes spp. mosquito vectors, the poorly planned urbanization of many cities in developing countries, an increased number of susceptible human hosts, and the rapid spread of DENV serotypes through global human travel networks (2–4). According to the World Health Organization, ≈3 billion persons living in >100 countries are at risk of being infected at least once annually by 1 of the 4 DENV serotypes (4).
In Brazil, DENV-1–3 were responsible for ≈5 million cases of DENV infection during 1990–2009, resulting in >15,000 reported cases of dengue hemorrhagic fever and ≈1,000 DENV-related deaths (5–7). DENV-4 reemerged in Brazil in 2010, 28 years after it was last detected in the country; the site of the reemergence was Roraima State, northern Brazil (8), the same state in which DENV-4 had last been detected in 1982 (9). Brazilian Ministry of Health data for 2010 and 2011 show there were 1,666,208 cases of DENV infection, including 26,659 severe cases and 1,097 associated deaths (10).
We describe the genetic characterization and spatiotemporal patterns of spread for DENV-4 strains isolated from 4 Brazilian states: Pará, Amazonas, and Roraima in northern Brazil and Bahia in northeastern Brazil. To characterize the origins of DENV-4 reemergence, we performed discrete Bayesian phylogeographic analysis on 98 full-length DENV-4 genomes and compared the results with those of a similar analysis on 314 envelope gene sequences.
We included 16 DENV-4 isolates in this study (Table). The viruses corresponded to low-passage virus strains (passage no. 1) from C6/36 cells obtained from the Department of Arbovirology and Hemorrhagic Fevers, Instituto Evandro Chagas, Brazilian Ministry of Health (Ananindeua, Brazil).
Sequencing, Assembly, and Accession Numbers
Nearly complete genome sequences were obtained by using high-throughput sequencing on a GS FLX+ System (454 Life Sciences, Branford, CT, USA) as described (11). The sequences of the 5′ and 3′ termini were obtained by using 5′ and 3′ rapid amplification of cDNA ends systems (Invitrogen, Carlsbad, CA, USA), according to the manufacturer’s instructions, with the following specific sets of primers: DENV-4 5′SP1 (5′-AKCCCTGTCTTGGGTCCAGC-3′), DENV-4 5′SP2 (5′- TTGAACGCCTCTTGAAGGTC-3′), DENV-4 5′SP3 (CTCTTGAAGGTCCAGGTCTA), DENV-4 3′SP1 (5′-ATATCTGAATGGCAGCCATC-3′), and DENV-4 3′SP2 (5′-TCACTGGCTGTTTCTTCTGCT-3′). All 5′ and 3′ rapid amplification of cDNA ends amplicons were cloned into a plasmid bacterial system by using the TOPO TA Cloning Kit (Invitrogen) and directly sequenced (in both directions) by using the ABI Prism BigDye Terminator v1.1 Cycle Sequencing Kit on an ABI Prism 3130 DNA analyzer (Applied Biosystems, Foster City, CA, USA) with the M13F/M13R set of primers (Invitrogen).
Data generated by the GS FLX+ System and ABI 3130 platforms were assembled by using the mapping reference method as implemented by the GS Reference Mapper Program (available in Newbler v.2.6 software, http://454.com/products/analysis-software/index.asp). The mapper program was used to rearrange the reads against a given reference sequence by using the following default parameters: input = 20 bp, all contig threshold = 100, large contig threshold = 200, minimum overlap length = 40, minimum overlap identity = 70%, k-mer = 12 (seed step), and k-mer = 16 (seed length). A total of 16 new full-length DENV-4 sequences were obtained and deposited in GenBank (accession nos. JQ513330–JQ513345).
Nucleotide Data Compilation and Alignment
We complemented our data with available full genome and envelope protein sequences for a total of 98 full genome (10,624 nt) and 314 full-length DENV-4 envelope protein (1,485 nt) sequences. These included 18 sequences from Brazil (3 from Belém, Pará State; 1 from Santarém, Pará State; 5 from Manaus, Amazonas State; 8 from Boa Vista, Roraima State; and 1 from Salvador, Bahia State); 26 full genome and 75 envelope protein sequences from the Caribbean region (Puerto Rico, Dominica, Dominican Republic, Jamaica, Puerto Rico, Trinidad and Tobago, Barbados, Montserrat, Martinique, and Bahamas); 8 full genome and 10 envelope protein sequences from Colombia; 34 full genome and 48 envelope protein sequences from Venezuela; 9 full genome and 88 envelope protein sequences from Mainland Southeast Asia (Thailand-Bangkok, Vietnam, Cambodia, Myanmar, and Peninsular Malaysia); and 3 full genome and 22 envelope protein sequences from Maritime Southeast Asia (Philippines, Indonesia, Singapore, and East Malaysia).
We aligned the sequences by using MAFFT (12). After we manually edited the resulting alignment by using Se-Al (13), the total lengths of the full genome and envelope protein alignments consisted of 10,624 bp and 1,485 bp, respectively. All sequences were screened for recombination by using the Phi-test (14), which is available in the SplitsTree4 program (15).
Evolutionary Reconstruction of DENV-4 Dispersal
The parameters of a full probabilistic model of evolutionary history, including timed sequence evolution and spatial dispersal, were estimated by using a discrete Bayesian asymmetric diffusion approach (16) implemented in the BEAST software package (13,17). We used the Bayesian skyride model (18) as a flexible tree prior and a general time-reversible model with a discretized gamma distribution (general time-reversible + 4Γ) to account for among-site rate variation. To calibrate the time scale of the trees, we obtained isolation dates (in years) from the GenBank annotations, and to accommodate rate variation among lineages, we used a lognormal relaxed molecular clock approach (19).
This approach enabled us to estimate ancestral spatial locations throughout the phylogenetic history while accounting for uncertainty in the phylogenetic and diffusion process (17). For the full genome and envelope protein datasets, 3 Markov-chain Monte Carlo analyses were run for 50 million states and sampled once every 10,000 states. We used the BEAGLE library (13) together with BEAST (13) to augment the computational speed. After we removed 10% of the burn-in, we combined the runs by using LogCombiner (www.molecularevolution.org/software/phylogenetics/beast). Maximum clade credibility (MCC) trees were summarized by using TreeAnnotator and visualized by using FigTree (12). We use the SPREAD application (20) to visualize and convert the estimated divergence times and spatial estimates annotated in the MCC trees to a keyhole markup language file (the files are available from authors upon request). All evolutionary parameters are reported as posterior means along with their 95% Bayesian credibility intervals.
To obtain the expectations for the location state transitions, we estimated Markov jump counts (18,21) along the branches of the posterior tree distribution (22). We applied kernel density estimation in R to summarize Markov jump densities through time from particular locations.
Asian DENV-4 Genotype I in Northern Brazil
The analysis (Phi-test implemented in SplitsTree) of the 98 DENV-4 full genomes did not provide support for recombination. Therefore, all 98 complete nucleotide (10,624 nt) and 314 (1,485 nt) envelope protein sequences were included in the Bayesian phylogeographic analyses by using a discrete asymmetric diffusion model and a flexible demographic prior (Technical Appendix Table).
The complete genome and envelope protein MCC genealogies confirmed that 3 major phylogenetic groups were compatible with the genotypes currently established for DENV-4: genotypes I–III (Figure 1). In particular, genotype I (n = 5 full genome, n = 85 envelope protein sequences) included strains that were isolated in Southeast Asia and a unique strain (BHI 3581 Be H 781363) that was isolated from an autochthonous febrile patient in 2011 in Salvador, Bahia State. The dense sampling in the envelope protein dataset made it possible for us to establish a relatively recent common ancestor with an Asian genotype I dating back to 2004 (95% Bayesian credible interval [BCI] 2002–2007). In particular, results from the genomic and envelope protein phylogeographic analyses confirmed with high posterior probability (>0.99) that the Brazilian strain from Bahia was introduced from Mainland Southeast Asia.
Multiple Introductions of DENV-4 Genotype II in Brazil
Genotype II (n = 86 full genome, n = 223 envelope protein sequences) included the other 17 strains sampled in Brazil (Municipalities of Manaus, Belém, Santarém, and Boa Vista) and in Venezuela, Colombia, Central America, the Caribbean, and Southeast Asia. Genotype III (n = 2 full genome, 4 envelope protein sequences) included strains that were exclusively isolated from Southeast Asia (Thailand). On the basis of the envelope gene analysis, we reconstructed 3 introductions of genotype II DENV-4 in Brazil (Figure 1, left panel). First, the envelope protein MCC tree indicated a Brazilian cluster of strains identified in Boa Vista and Manaus as part of the major South American genotype II cluster. Second, the full genome MCC tree broke the cluster in Boa Vista and Manaus into 3 different lineages in the same genotype II cluster. Third, the large South American genotype II cluster, as part of the full genome MCC genealogy (Figure 1, right panel), indicates additional introductions of genotype II into Brazil. Thus, we identified a total of 6 introductions of DENV-4 in Brazil.
The most recent common ancestor of genotypes I–III is believed to have a Southeast Asian origin (7). In line with recent findings (23), we estimate that genotype II, the predominant genotype, was likely first introduced into the Americas (Caribbean region) around 1978 (95% BCI 1976–1979) (Figure 1; Technical Appendix Figure). Following a period of circulation in the Caribbean region, this genotype independently entered South American countries >4 times, as inferred respectively by the full genome and envelope protein analyses. Most of these introductions into South America seem to have led to relatively shallow phylogenetic clusters that probably reflect short-term survival of the virus in South America.
Nevertheless, our data provide support for a major genotype II cluster that reflects the successful establishment of this genotype in South American countries neighboring Brazil in or around 1994 (95% BCI 1992–1995). From Venezuela and Colombia, the genotype spread to Brazil >3 times in the last decade, as inferred by the full genome phylogeographic analyses (Figure 1, right panel). In contrast, the DENV-4 outbreak in Pará State (Santarém and Belém) was identified around 2008 (95% BCI 2006–2009) and is directly linked to strains from the Caribbean region. Overall, we found >3 distinct DENV-4 introductions into Roraima State: 1 from the Caribbean region ≈3 decades ago (indicated as 1 in Figure 1, left panel), and 2 from the strain that has become endemic in Venezuela within the last decade (indicated as 3 and 5 in Figure 1, right panel). One of these strains resulted in a well-supported cluster that dated back to about 2006 (95% BCI 2005–2007) and comprised 6 strains.
Global Patterns of DENV-4 Dissemination
To identify the major export and import locations for DENV-4 at a global scale, we summarized the number of introductions from and to each of the locations that were estimated by using Markov jump count approaches as part of the full genome analysis (Figure 2, panel A). Our results pinpoint Maritime Southeast Asia, the Caribbean region, and Venezuela as the 3 most prominent sources of DENV-4 export. In addition, to investigate the temporal dynamics of viral dispersal from the 3 locations in the evolutionary history of DENV-4, we estimated the Markov jump density from these locations through time (Figure 2, panel B). We distinguished different peaks, which appear to reflect 3 distinct waves of DENV-4 migration; the most recent wave corresponded to exportations from South American countries bordering Brazil, particularly Venezuela, where DENV-4 has been established longer.
After an initial establishment in the Caribbean region, genotype II most likely first dispersed to Brazil (Roraima State), Venezuela, Colombia, and, more recently, to Brazil (Belém, Pará State) again (Figure 3). Locations such as Santarém, Pará State, and Manaus, Amazonas State, do not have a direct relation to the Caribbean locations; instead, the strain found in Santarém is linked to a cluster of sequences from other locations in Pará State. This finding further indicates that there is ongoing circulation of DENV-4 within Brazil. Nevertheless, most recent introductions of genotype II in Brazil in the last decade seem to have been fueled by strains originating from bordering countries.
DENV-4 was not detected in Brazil for 28 years after the first clinical and laboratory reports of dengue fever cases in Roraima State during 1981–1982 (9). Since then, no additional cases were reported in the country until the reemergence of DENV-4 in Boa Vista, Roraima State, in 2010 (8). Following those first recent detections of DENV-4 in Brazil, the virus was identified in 2011 in the northern Brazilian states of Amazonas, Amapá, and Pará. Furthermore, DENV-4 was serologically detected in persons in several other states, demonstrating the potential for this virus to spread quickly to several regions (8).
Since the first detection of DENV-4 in Brazil in 1982 (8,9), partial genomic studies have confirmed that the original virus was directly associated with the Caribbean strains, which are known as genotype II. Genotype II has been responsible for several outbreaks in many countries in the Caribbean and South America (8,24–28). Phylogenetic analyses of different strains demonstrated the presence of 2 distinct genotypes (genotypes I and II) of DENV-4 in Brazil. The introduction of a new serotype and a distinct DENV-4 genotype (Asian genotype I) into Brazil after nearly 3 decades of no circulation highlights the potential for future outbreaks of genotype I in this country. The presence of a serotype against which the population is not immunized indicates that the country is at risk for a sharp increase in the number of dengue virus infections, including severe cases (8).
The phylogeographic analyses of full genomes and envelope genes confirmed the co-circulation of 2 distinct DENV-4 genotypes (I and II) in Brazil. Genotype II is most common in South America and the Caribbean region, and genotype I was represented by 1 strain isolated from Bahia State in northeastern Brazil. The analyses also confirmed the introduction of DENV-4 into Pará State from the Caribbean region and suggest that this introduction is very recent (Technical Appendix Figure).
Previous phylogenetic analyses have been performed on partial and a few complete DENV-4 genomes (8,29–31). The results we obtained by using the new set of complete genomes suggest that DENV-4 genotype II emerged and reemerged in Brazil from >3 distinct origins (Southeast Asia, the Caribbean region, and Venezuela), which demonstrates a dispersal pattern in Brazil that is far more complex than expected from standard epidemiologic data.
We demonstrate that full genome analysis complements the analysis of more widely available envelope protein sequences. Although a dense sampling may have its own particular benefits for divergence time and spatial diffusion estimation, full genomes can provide more phylogenetic resolution to unravel individual introductions from populations of closely related strains.
Our results indicate that DENV-4 was exported multiple times from the Caribbean region to northern South American locations before the virus became established in South America. Overall, the spatiotemporal patterns of DENV-4 revealed in South America over a long time scale are reminiscent of a source-sink model of virus dispersal (32), in which the location of primary source virus populations shifted from Maritime Southeast Asia to the Caribbean region and, most recently, to northern South America, albeit with co-circulation in each location. This dispersal pattern differs from the source-sink dynamics of seasonal influenza, in which each epidemic wave seems to be caused by emergence from the same source (33). Active molecular epidemiologic surveillance will be essential for better characterizing local source populations of DENV-4.
In conclusion, we found that 2 distinct genotypes (I and II) of DENV-4 are circulating in Brazil, and we provide insight into the origin and dispersal of the DENV throughout northern Brazil and areas of several South American countries. Further studies are needed to analyze complete genomes from other countries to which DENV-4 is endemic; such studies will more fully elucidate the geographic dispersal dynamics of DENV-4 in regions of the Americas where it is endemic.
Dr Nunes is a molecular biologist at the Instituto Evandro Chagas, Ananindeua, Brazil. His primary research interests are molecular epidemiology and phylogenetic and evolutionary aspects of arboviruses, hantaviruses, and other microbial agents.
This study was supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) grants (302987/2008-8 and 301641/2010-2) and a CNPq/Coordenação de Aperfeiçoamento de Pessoal de Nível Superior/Fundação de Amparo a Pesquisa do Estado do Pará grant (573739/2008-0). N.R.F. is supported by Fundação para a Ciência e Tecnologia (grant SFRH/BD/64530/2009). M.A.S. is partially supported by the National Science Foundation (grant DMS0856099. A.B.A. is supported by Fundação para a Ciência e Tecnologia (grant SFRH/BPD/65605/2009). M.A.S. and P.L. had financial support from the European Research Council under the European Community's Seventh Framework Program (FP7/2007-2013; under grant 278433-PREDEMICS and ERC grant 260864). The collaboration between M.R.T.N., N.R.F., and A.B.A. was supported by the Network of the Portuguese-Speaking Countries Community, Instituto de Higiene e Medicina Tropical, Universidade Nova de Lisboa, Lisboa, Portugal.
- Gubler DJ. Dengue and dengue hemorrhagic fever: its history and resurgence as a global public health problem. In: Gubler DJ, Kuno G, eds. Dengue and dengue hemorrhagic Fever. New York: CAB International; 1997. p. 1–22.
- Gibbons RV, Vaughn DW. Dengue: an escalating problem. BMJ. 2002;324:1563–6.
- Lindbäck H, Lindbäck J, Tegnell A, Janzon R, Vene S, Ekdahl K. Dengue fever in travelers to the tropics, 1998 and 1999. Emerg Infect Dis. 2003;9:438–42.
- World Health Organization. Dengue fever and dengue hemorrhagic fever. Geneva: the Organization; 2009.
- Ministério da Saúde/Secretara de Vigilância em Saúde MS/SVS. Dengue cases Technical report; 2009 [in Portuguese] [cited 2011 Dec 17]. http://portal.saude.gov.br/portal/arquivos/pdf/boletim_dengue_informe11maio.pdf
- Teixeira MG, Barreto ML, Costa MCN, Ferreira LDA, Vasconcelos PFC, Cairncross S. Dynamics of dengue virus circulation: a silent epidemic in a complex urban área. Trop Med Int Health. 2002;7:757–62.
- Dussart P, Baril L, Petit L, Beniguel L, Quang LC, Ly S, Clinical and virological study of dengue cases and the members of their households: the multinational DENFRAME project. PLoS Negl Trop Dis. 2012;6:e1482.
- Temporão JG, Penna GO, Carmo EH, Coelho GE, do Socorro Silva Azevedo R, Teixeira Nunes MR, Dengue virus serotype 4, Roraima State, Brazil. Emerg Infect Dis. 2011;17:938–40.
- Osanai CH, Travassos da Rosa AP, Tang AT, do Amaral RS, Passos AD, Tauil PL. Dengue outbreak in Boa Vista, Roraima [in Portuguese]. Rev Inst Med Trop Sao Paulo. 1983;25:53–4.
- Ministério da Saúde/Secretaria de Vigilância em Saúde MS/SVS. Sempre é hora de combater a dengue [cited 2011 Dec 17]. http://www.combatadengue.com.br/
- Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005;437:376–80.
- Katoh K, Kuma K, Toh T, Miyata T. MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005;33:511–8.
- Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.
- Bruen TC, Philippe H, Bryant D. A simple and robust statistical test for detecting the presence of recombination. Genetics. 2006;172:2665–81.
- Huson DH, Bryant D. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006;23:254–67.
- Edwards CJ, Suchard MA, Lemey P, Welch JJ, Barnes I, Fulton TL, Ancient hybridization and an Irish origin for the modern polar bear matriline. Curr Biol. 2011;21:1251–8.
- Lemey P, Rambaut A, Drummond AJ, Suchard MA. Bayesian phylogeography finds its roots. PLOS Comput Biol. 2009;5:e1000520.
- Minin VN, Suchard MA. Counting labeled transitions in continuous-time Markov models of evolution. J Math Biol. 2008;56:391–412.
- Drummond AJ, Ho SY, Phillips MJ, Rambaut A. Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006;4:e88.
- Bielejec F, Rambaut A, Suchard MA, Lemey P. SPREAD: spatial phylogenetic reconstruction of evolutionary dynamics. Bioinformatics. 2011;27:2910–2.
- Talbi C, Lemey P, Suchard MA, Abdelatif E, Elharrak M, Nourlil J, Phylodynamics and human-mediated dispersal of a zoonotic virus. PLoS Pathog. 2010;6:e1001166.
- O’Brien JD, Minin VN, Suchard MA. Learning to count: robust estimates for labeled distances between molecular sequences. Mol Biol Evol. 2009;26:801–14.
- Allicock OM, Lemey P, Tatem AJ, Pybus OG, Bennett SN, Mueller BA, Phylogeography and population dynamics of dengue viruses in the Americas. Mol Biol Evol. 2012;29:1533–43.
- Centers for Disease Control and Prevention. Dengue type 4 infections in US travelers to the Caribbean. MMWR Morb Mortal Wkly Rep. 1981;30:249–50.
- Forshey BM, Morrison AC, Cruz C, Rocha C, Vilcarromero S, Guevara C, Dengue virus serotype 4, northeastern Peru, 2008. Emerg Infect Dis. 2009;15:1815–8.
- Fernández J, Vera L, Tognarelli J, Fasce R, Araya P, Villagra E, Detection of dengue virus type 4 in Easter Island, Chile. Arch Virol. 2011;156:1865–8.
- Nogueira RM, Eppinghaus AL. Dengue virus type 4 arrives in the state of Rio de Janeiro: a challenge for epidemiological surveillance and control. Mem Inst Oswaldo Cruz. 2011;106:255–6.
- Shu PY, Su CL, Liao TL, Yang CF, Chang SF, Lin CC, Molecular characterization of dengue viruses imported into Taiwan during 2003–2007: geographic distribution and genotype shift. Am J Trop Med Hyg. 2009;80:1039–46.
- Lanciotti RS, Gubler DJ, Trent DW. Molecular evolution and phylogeny of dengue-4 viruses. J Gen Virol. 1997;78:2279–84.
- Dussart P, Lavergne A, Lagathu G, Lacoste V, Martial J, Morvan J, Reemergence of dengue virus type 4, French Antilles and French Guiana, 2004–2005. Emerg Infect Dis. 2006;12:1748–51.
- de Souza RP, Rocco IM, Maeda AY, Spenassatto C, Bisordi I, Suzuki A, Dengue virus type 4 phylogenetics in Brazil 2011: looking beyond the veil. PLoS Negl Trop Dis. 2011;5:e1439.
- Holmes EC. The evolutionary history and phylogeography of human viruses. Annu Rev Microbiol. 2008;62:307–28.
- Rambaut A, Pybus OG, Nelson MI, Viboud C, Taubenberger JK, Holmes EC. The genomic and epidemiological dynamics of human influenza A virus. Nature. 2008;453:615–9.