Emerg Infect DisEmerging Infect. DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention15498159332027504-019710.3201/eid1009.040197ResearchResearchGenetic Divergence and Dispersal of Yellow Fever Virus, BrazilYellow Fever Virus Epidemiology, BrazilVasconcelosPedro F.C.*BryantJuliet E.Travassos da RosaAmelia P.A.TeshRobert B.RodriguesSueli G.*BarrettAlan D.T.World Health Organization Collaborating Center for Arbovirus Reference and Research, Instituto Evandro Chagas, Belém, Brazil;University of Texas Medical Branch, Galveston, Texas, USAAddress for correspondence: Pedro F.C. Vasconcelos, Arbovirus Department, Instituto Evandro Chagas, Av. Almirante Barroso, 492, 66090-000, Belém, Brazil; fax: +55-91-226-1284; email: pedrovasconcelos@iec.pa.gov.br9200410915781584

Examining viral isolates collected over 66 years shows divergence into clades and potential dispersal by human migration.

An analysis of 79 yellow fever virus (YFV) isolates collected from 1935 to 2001 in Brazil showed a single genotype (South America I) circulating in the country, with the exception of a single strain from Rondônia, which represented South America genotype II. Brazilian YFV strains have diverged into two clades; an older clade appears to have become extinct and another has become the dominant lineage in recent years. Pairwise nucleotide diversity between strains ranged from 0% to 7.4%, while amino acid divergence ranged from 0% to 4.6%. Phylogenetic analysis indicated traffic of virus variants through large geographic areas and suggested that migration of infected people may be an important mechanism of virus dispersal. Isolation of vaccine virus from a patient with a fatal case suggests that vaccine-related illness may have been misdiagnosed in the past.

Keywords: yellow fever virusBrazilevolutionepidemiologyenzootic3′NCRgenetic divergenceresearch

Yellow fever virus (YFV) is transmitted by the bite of infected mosquitoes and produces a severe hemorrhagic fever in humans. Despite a safe and effective vaccine (17D), YFV continues to be a public health problem in tropical areas of Africa and South America (13). A recent upsurge in YFV activity in Brazil and the reinfestation of urban areas with the vector mosquito Aedes aegypti have stretched disease surveillance and control resources to their limits. During 2000, YFV occurred near Brasilia, the capital city (4); in 2001, YFV spread to new areas of Minas Gerais outside the currently recognized enzootic zone (5).

In South America, YFV is maintained in enzootic cycles involving monkeys and forest-canopy mosquitoes of the genera Haemagogus and Sabethes (6,7). In Brazil, three geographic zones have been defined where YFV circulates (7) (Figure 1): 1) the region of endemicity, in which the virus is maintained in mobile monkey populations and where human cases are sporadic and rare; 2) transitional zones of emergence, in which contact is frequent between monkey and human populations (and infected mosquito vectors); and 3) regions of epidemicity where the density of susceptible human populations and competent vector species are both high, and the potential for explosive urban outbreaks is great. Brazil currently has a population of about 176 million; approximately 30 million people live in the YFV-endemic zone, 18 million live in the zone of emergence, and 128 million reside along the Atlantic coast in the YFV-free zone (8,9).

Regions where yellow fever is endemic in Brazil.

Many aspects of the molecular epidemiology and transmission cycles of YFV in the forests of South America are poorly understood, and previous studies of YFV in South America were limited to a relatively small number of isolates. We examined the genetic diversity of 79 YFV strains isolated from Brazil over 67 years and mapped the distribution of variants to investigate patterns of virus divergence and dispersal.

Material and MethodsStudy Area

Brazil is a country of enormous size and diversity, covering 8,512,000 km2. The country is divided into 26 states and the Federal District. These make up five major geographic regions, characterized by broadly different climate and vegetation zones and a highly variable distribution of the human population. The five regions (Figure 1) consist of the northern Amazon region, northeastern region (caatinga, very dry), central-western region (swamp and savannah), southeastern region (most heavily populated with an extensive system of roads and railways), and the more temperate southern region bordering Argentina and Paraguay (10).

Virus Strains

Strain-specific data such as geographic locality, passage history, source of isolate, and clinical outcome (for selected human cases) for the 79 Brazilian YFV strains used in this study are provided in Table A1. All strains were originally isolated in suckling mice and cultured once in C6/36 cells to produce seed stocks (Table A1) at the Instituto Evandro Chagas (IEC), Belém, Brazil. Methods for cell culture and virus growth have been previously described (1113), and standard laboratory precautions within biosafety level 3 facilities were undertaken to prevent cross-contamination of strains. Thirty-eight (48%) virus strains were from humans (24 from patients who died); 7 (9%) were from monkeys; and 34 (43%) were from mosquito pools, mainly Haemagogus janthinomys. Viruses represented a period of 67 years, but with unequal sample distribution: 15% of strains were obtained from 1935 to 1969, 43% were obtained from 1970 to 1989, and 42% were obtained from 1990 to 2001. The year of virus isolation is indicated in the sequence identification (e.g., Brazil35, isolated in 1935). Isolates were from 12 states: Amapá (n = 1), Bahia (n = 1), Federal District (n = 1), Goiás (n = 10), Maranhão (n = 6), Minas Gerais (n = 7), Mato Grosso (n = 4), Mato Grosso do Sul (n = 5), Pará (n = 39), Rondônia (n = 2), Roraima (n = 1), and Tocantins (n = 2).

Phylogenetic Analysis

Supernatants of YFV-infected Vero cells were obtained, and viral RNA was extracted by using a commercial kit (Qiagen, Valencia, CA) and processed according to the manufacturer's instructions. RNA obtained was stored at –70°C. The genomic-sense degenerate primer EMF (5´ TGGATGACSACKGARGAYAT) and genomic-complementary primer VD8 (5´ GGGTCTCCTCTAACCTCTAG) were used for reverse transcription–polymerase chain reaction amplification of a 595-bp fragment comprising 255 nucleotides of NS5 and 340 nucleotides of 3´NCR (14,15). PCR products were screened by agarose gel electrophoresis. Bands were recovered with a gel extraction kit (Qiagen) and directly sequenced with an ABI automatic sequencer at the University of Texas Medical Branch protein chemistry core facility.

Sequence editing and alignments were performed with Vector NTI (Informax, Frederick, MD), and additional manual editing of alignments was performed with the GCG Wisconsin Package Version 10.3 (Accelrys, San Diego, CA). The PAUP* program (16) was used to infer phylogenetic trees by the neighbor-joining method with Kimura 2-parameter distance corrections. Support for individual clades was determined by nonparametric bootstrapping with 1,000 replicates. The tree was rooted with a sequence of the prototype West African YFV strain Asibi (17) and the 17DD substrain vaccine virus (18).

ResultsGenetic Diversity of Brazilian YFV Strains

Sequence data were obtained for a genomic region spanning the terminal portion of NS5 and proximal region of the 3´NCR for 54 Brazilian YFV strains. In addition to the sequence data generated in this study, partial or complete 3´NCR sequences were available for an additional 25 Brazilian YFV isolates (13,14), which yielded the full dataset of 79 YFV strains from 12 states (alignment length, 576 nt). The NS5/3´NCR sequences contained 148 variable sites, of which 89 were informative. Thirty-two of the informative sites (36%) fell within the NS5 coding portion of the sequence, and 57 (64%) fell within the 3´NCR. Amino acid pairwise divergence among the partial NS5 sequences (85 amino acids in length) ranged from 0% to 4.7% (mean 2.2%).

A phylogenetic tree of the 79 Brazilian YFV sequences is shown in Figure 2. The tree is rooted with the homologous region of the Asibi prototype strain (parental virus to the 17D vaccine), isolated in Ghana in 1927 (17). Mosquito- and vertebrate-derived sequences were distributed randomly throughout the phylogenetic tree. With the exception of one 1983 isolate from Rondônia (BeH 413820) and one 1975 isolate from Aripuana in Mato Grosso (BeH 291597), all Brazilian YFV strains formed a single monophyletic clade. Two major subclades were evident: a subclade comprising isolates from Pará dating from 1954 to1968 and a subclade containing all remaining isolates from 1969 to 2001.

Brazilian NS5/3´NCR phylogeny (576 nt) based on yellow fever isolates (neighbor-joining tree, Kimura 2-parameter distance correction, midpoint rooted). Geographic origin of isolates is indicated on map. 1: North (AC, Acre; AM, Amazonas; AP, Amapá; PA, Pará; RO, Rondônia; RR, Roraima; TO, Tocantins). 2: Northeast (AL, Alagoas; BA, Bahia; CE, Ceará; MA, Maranhão; PB, Paraiba; PE, Pernambuco; PI, Piaui; RN, Rio Grande do Norte; SE, Sergipe). 3: Central West (DF, Distrito Federal; GO, Goiás; MT,Mato Grosso; MS, Mato Grosso do Sul). 4: Southeast (ES, Espírito Santo; MG, Minas Gerais; RJ, Rio de Janeiro; SP, São Paulo). 5: South (PR, Paraná; SC, Santa Catarina; RS, Rio Grande do Sul). Colors correspond to genetic clade structure. Black dots refer to isolates with unresolved phylogenetic position.

Variation within the Dominant Subclade

Genetic variation within the dominant subclade of Brazilian YFV strains showed a complex pattern of relationships that demonstrated both geographic and temporal associations. Varying levels of bootstrap support were evident for four clusters (groups 1A–1D) (Figure 2): group 1A included isolates from the 1972–1973 outbreak in Goiás and a 1984 isolate from Pará, group 1B represented sporadic human cases and samples obtained during routine surveillance in western and central regions from 1984 to 1996, group 1C consisted of isolates from sporadic cases and surveillance in eastern central states from 1971 to 1998, and group 1D represented epizootic activity from 1998 to 2001. The phylogenetic position of 12 additional isolates from Pará (10), Amapá (1), and Mato (1) was not easily resolved (Figure 2), so relationships among these strains require further sequence analysis.

The four isolates from the 1972–1973 epidemic in Goiás (group 1A, Brazil73a, 73b, 73c, and 73d) differed by 2–3 nt (0.34%–0.5%) over the length of the NS5/3´NCR fragment (595 nt). One isolate obtained 11 years later (Brazil84e) from a pool of Haemagogus janthinomys mosquitoes collected in far northern Pará (Faro), was nearly identical to the 1973 Goiás strains (2–3 nt difference).

Group 1B consisted of 11 isolates. These included two isolates from patients who died (Brazil84h and Brazil91b), isolates from patients in Mato Grosso do Sul (Brazil92c) and Maranhão (Brazil93a), isolates from Haemagogus mosquitoes in western Pará (Brazil84d) and Rondônia (Brazil96a), and four isolates obtained in 1992 from mosquito pools (Brazil92a, 92b, 92d, and 92e). Although most strains in group 1B were identified in western and central regions (over a distance as large as 3,500 km2), one isolate identified in this cluster was from the northeast region in Barra do Corda, Maranhão (Brazil93a). Isolates collected from contemporaneous Haemagogus and Sabethes mosquito pools (Brazil92a, 92b, 92d, and 92e) were 100% identical.

Group 1C formed the largest cluster, with 22 strains diverging by 0% to 3.8% (0–27 nt, bootstrap 80%). These variants were distributed during an extended period from 1971 to 1994 throughout the northern Amazon region (Pará: Brazil87, 91a, 84c, 84g, 71, 78a, 84b, 94c, and 98a; Maranhão: Brazil80c, 82a, 93b, 94f, and 95), in central Goiás (Brazil88a, 80a, and 91c), and near the edge of the enzootic zone in Minas Gerais (Brazil89, 88b, 94d, 94b, and 94a). A pair of Minas Gerais isolates (Brazil94b and 94a) obtained from Haemagogus and Sabethes mosquito pools showed a 100% match across the length of the sequence fragment.

With one exception (Brazil98a), all isolates collected during 1998 to 2001 fell within a single cluster (group 1D) that had 0%–2.7% (0–16 nt) divergence. Although group 1D consistently clustered on neighbor-joining and parsimony trees, bootstrap support was low (53%). The distribution of samples collected during this period closely reflected the southward dispersal of epizootic activity. In 1998, a large number of YF cases occurred in northern Brazil, especially in Pará. In 1999 through 2000, a few cases were reported in Bahia and Tocantins, but the largest number were in central Goiás (19). In 2001, YF cases were detected in the transitional zone of Minas Gerais and further east. The isolate from Tocantins collected in 2000 (Brazil2000a) was characterized by an exceptionally large number of substitutions in the 3´NCR region. A single isolate (Brazil98a) collected from a howler monkey in Afua, Pará during the early period of the epizootic did not cluster with the other group 1D strains.

Subclade of Older Pará Strains

A group of ten isolates from Pará dating from 1954 through 1968 differed from all other Brazilian strains by 5.7% ± 1.6%. These results confirmed previous observations based on analysis of prM/E gene sequences (14), indicating 7.8% ± 2.0% divergence for this group of older Pará isolates. To date, 29 YFV strains from Pará dating from 1969 through 1999 have been examined; none of these more recent strains cluster with the earlier clade. Sequence divergence of the older Pará strains approached the threshold level used to define separate YFV genotypes (8%–9%) (20).

Alignment of previously published sequences for the structural gene region (223 codons of prM/E) (14) showed that 8 of the 10 older Pará strains had one or more mutations within the fusion peptide of the E protein (Figure 3) that is highly conserved among all flaviviruses; it is located in domain II (E98–E110) and plays a role in mediating the acid-catalyzed fusion of virions with target cell membranes (22). Three strains from 1968 (Brazil68a, 68c, and 68d) showed a C→F substitution at E105, eliminating one of the disulfide bonds that forms the structural architecture of domain II (22). Brazil68b had a D→G mutation at E98, and Brazil54 and Brazil68a shared an R→K mutation at E99. Three of the strains (Brazil55c, Brazil60, Brazil62b) had a G→S substitution at E100. The functional importance of these mutations is unknown. Of the older Pará strains, only Brazil55b and Brazil62a had the consensus sequence for the fusion peptide (Figure 3). In addition, 12 non-Pará subclade strains exhibited substitutions in the fusion peptide sequence.

Sequence alignment of the fusion peptide of the envelope (E) gene of selected yellow fever virus (YFV) strains (E98–E110). The Asibi prototype strain indicates the conserved sequence present in the majority of YFV strains and other mosquitoborne flaviviruses. A salt bridge between residues Asp E98 and Lys E110 generates the "CD loop" of residues E100–E108 (21).

Preliminary phenotypic differences among three of the older Pará strains (Brazil55b, Brazil60, and Brazil68c) in the standard mouse neuroinvasiveness model (i.e., intraperitoneal injection of 8-day-old suckling mice) (23) indicated that the strains with substitutions in the fusion peptide were less lethal (average lethal dose-50 [LD50] for Brazil60 and Brazil68c = 5.5 log10 tissue culture infective dose-50 [TCID50]) than the strain with the consensus sequence (Brazil55b) (LD50 = 0.1 log10 TCID50). Average survival times of the infected mice were 9.2, 11, and 8.8 days, respectively, for the three strains. All three strains achieved titers 6–7 log10 TCID50/mL when grown in Vero cell culture. Brazil55b had a longer passage history than other YFV strains (4 passages in suckling mouse brain) (Table A1); thus the altered phenotype may be partly attributed to selection during repeated mouse passage.

Additional Brazilian Isolates

Two Brazilian strains (Brazil83 from Rondônia and Brazil75 from Mato Grosso) failed to group within either of the two major subclades. Brazil83 appears to represent South American genotype II, as it matched 97.2% to the homologous NS5/3´NCR region of Peruvian YFV strains (14) and diverged from other Brazilian strains by 6.2% to 13.8% (mean 10.2%). Brazil75 was a vaccine virus, as it matched the 17DD vaccine virus by 99.8%. Sequence comparison of 17DD with Brazil75 revealed a single mutation (C→T) at position 16 of the 3´NCR (13).

Discussion

Brazil accounts for approximately 25% of all YF cases reported from South America (4). YFV activity has been reported from each of the five regions in the country (2,24). From 1950 to 2003, no single region consistently produced the highest number of YF cases; however, YFV activity in the dry northeast was rare (Figure 4). Overall, the largest number of cases was reported from Goiás (central-western) and Pará (northern) regions.

Yellow fever incidence in Brazil by region, 1950–2003.

Previous studies of the genetic relationships among global variants of YFV have indicated that divergence across the length of the ≈11-kb genome is relatively uniform, and the 3´NCR contains useful markers for subtype-specific distinctions (13,25). Our data showed a single genotype of YFV circulating in Brazil (South American genotype I), with the exception of a single strain from Rondônia (South America genotype II). The Brazilian strains made up two subclades: a group of older strains from Pará from 1954 to 1968 and a larger group dating from 1969 to the present. The clear genetic distinction between strains isolated before 1968 (10 in total) and all subsequent strains suggests that the older Pará strains have been replaced by a dominant new lineage. Several members of the older Pará subclade had one or more mutations within the fusion peptide region of the E protein, including substitutions that disrupted conserved disulfide bridges (Figure 3). Although the functional importance of these substitutions is unknown, changes in the fusion peptide would affect protein folding and conformation and binding interactions with target membranes (22). Continued sampling and surveillance of YFV strains from Pará are necessary to confirm whether the variants represented in this subclade have truly become extinct or whether they are being maintained in cycles that have yet to be identified.

At least three articles in the past 50 years have reported YFV epizootics that began in northern and western Amazonian regions and then spread south into Paraná, Santa Catarina, and Rio Grande do Sul. The first was a large epidemic that swept south from Goiás in the 1930s and 1940s (26). The second was in 1963, when an epizootic started in Mato Grosso and extended eventually as far south as Missiones and Corrientes Provinces in northern Argentina by 1966 (27). In 1972 to 1973, an outbreak occurred in Goiás, and although investigations at the time suggested the epizootic was highly localized (28), cases reported from Paraguay in the following years were attributed to spread of viruses from Goiás. Virus maintenance and dispersal have presumably involved sequential infections of migrating groups of monkeys (2932). Four isolates collected during the 1972–1973 Goiás outbreak were included in this study; these isolates (group 1A, Brazil73a–d) had nearly identical (99.6%) NS5/3´NCR sequences, and a fifth isolate with a highly similar sequence was obtained 10 years later from Faro, in northern Pará. Other than the single Faro isolate (Brazil84e), no additional descendents have been identified of those outbreak strains.

YFV transmission was particularly active during the rainy season in Maranhão in 1993 to 1994 (33). Serologic studies of rural and urban populations at the time indicated an overall attack rate of 75 per 1,000; incidence of clinical disease was 3.5 per 1,000 persons in urban areas and 4.2 per 1,000 in rural areas. Five of the six Maranhão isolates from this time were in group 1C (Figure 2), whereas one isolate from a patient with a fatal case (Brazil93a) appeared to be related to group 1B (Figure 2). In addition to the 1993–1994 outbreak in Maranhão, the subclade of group 1C strains was also associated with sporadic cases throughout 1971 to 1998 in Pará, Minas Gerais, and Goiás.

The most recent increase in epizootic YFV activity in Brazil occurred from 1998 to 2001, with cases distributed over a large region covering eight states (4). Several human cases were reported to have originated from the Chapada dos Veadeiros National Park, a tourist canyon located near Brasilia, Goiás. The proximity to major cities raised alarm, and reports at the time designated Goiás as the epicenter of the outbreak (4). The phylogenetic evidence presented here indicates that YFV activity in 1998 in Pará, in 2000 in Goiás, and in 2001 in Minas Gerais was all part of one continuous epizootic, which dispersed genetic variants from the 1D group of viruses (Figure 2).

Investigations in 2000 and 2001 indicated that YFV activity had expanded beyond the typical borders of the enzootic zone, into areas where the virus had not been reported for >100 years. The appearance of nearly identical variants across very large distances over short periods (e.g., group 1D spanning >3,000 km within 1 month, group 1B strains with isolates >2,000 km apart within 1 year) suggests that humans, rather than other primates or mosquitoes, may be partly responsible for the spread of YFV variants. Because Haemagogus mosquitoes are forest species and are relatively fragile, patterns of traffic and commerce would not likely lead to translocation of infected mosquitoes. From 1998 to 2001, the virus may have been transported from Pará (34) to Goiás> (4) and then to Minas Gerais by asymptomatic carriers or viremic persons in the prodromal phase. The general trend in the southward movement of the virus (group 1D, Figure 2) could be interpreted to reflect the pattern of labor migration from less populated areas of the northern Amazon to cities of the southeast. However, humans are not believed to play a major role in either virus maintenance or dispersal within South America, because human contact with the forest species Haemagogus janthinomys is infrequent, and the length of the viremic period in infected humans is brief (3–4 days).

The identification of one strain from western Rondônia (Brazil83) with high identity to South American genotype II strains (from Peru and Bolivia) suggests that the two South American YFV genotypes cocirculate in regions of western Brazil. The genetic divergence and distribution of YFV are similar that of Oropouche virus (an Orthobunyavirus transmitted by midges). Like YFV, Oropouche virus has split into two distinct genotypes in South America with separate eastern and western regions of circulation, and the only region where the two genotypes have been found to cocirculate is in western Brazil, in Rondônia (35). The ecologic correlates for these patterns remain uncertain; however, they suggest that surveys of the border regions between Brazil, Peru, and Bolivia may identify additional YFV genotypes.

One unanticipated result of this study was identifying a vaccine virus (Brazil75) from what had been presumed to be a case of natural exposure to wild-type YFV. This isolate was obtained from the blood of a patient who died in Aripuanã, Mato Grosso; this patient had been vaccinated 5 days before becoming ill and died 9 days postvaccination. Serious adverse effects resembling wild-type YF (viscerotropic disease) have only recently been reported (19,36,37). The complete genomic sequence of Brazil75 is currently being sought to address the question of vaccine reversion. Although no controlled studies have examined the safety or efficacy of YF vaccination in immunosuppressed patients, these findings underscore the importance of evaluating patient history before delivering vaccine and also suggest that other cases of vaccine-related illness may have been misdiagnosed in the past.

In summary, we describe considerable genetic variability among YFV variants circulating in Brazil and identify clusters of strains associated with epizootics in different geographic regions. Brazilian YFV strains have diverged into two subclades, one of which has become the dominant lineage in recent years. We suggest a potential role for human migration in mediating virus dispersal. Expansion of YFV outside the enzootic zone presents an ongoing risk for reintroduction of the virus to urban areas and highlights the need for continued surveillance and control.

Suggested citation for this article: Vasconcelos PFC, Bryant JE, Travassos da Rosa APA, Tesh RB, Rodrigues SG, Barrett ADT. Genetic divergence and dispersal of yellow fever virus, Brazil. Emerg Infect Dis [serial on the Internet]. 2004 Sep [date cited]. http://dx.doi.org/10.3201/eid1009.040197

Acknowledgments

We thank Zouraide Guerra Costa for sharing important unpublished data.

This study was financially supported by CNPq (Brazilian Agency for Scientific and Technologic Development, process 302770/02-0) and the National Institutes of Health (grant AI 50175). P.F.C. Vasconcelos was supported by The Lancet Publishing Group through the 2002 Lancet International Fellowship Award at Galveston. The contributions of J.E. Bryant were supported by the Centers for Disease Control and Prevention Fellowship Training Program in Vector-Borne Infectious Diseases, grant no. T01/CCT622892.

Virus isolates used in this study
IsolateStrain IDPassage levelaSourcebStatecCommunityBiotoped/Climatee
JSSBrazil35 (38)Mosq1, SM8HumanunkMTNo dataWC/EHW
Be H111Brazil54 (39)c6/36, SM10Human *PAOribocaEA/EHW
BeAr 162Brazil55BSM4, C6/36#1Hg. janthinomysPAPirelli MaritubaEA/EHW
BeAr 189Brazil55C (39)sm1 c6/36#1Sabethes sp.PAPirelli MaritubaEA/EHW
BeAn 23536Brazil602SM1 c6/36#1Cebus sp.PABelém Brasília Km94EA/EHW
BeAr 46299Brazil62A (39)c6/36Haemagogus sp.PABelém Brasília Km94EA/EHW
BeAr 44824Brazil62B (39)SM1 c6/36#1Haemagogus sp.PABelém Brasília km87EA/EHW
BeAn142028Brazil68A2c6/36#1Saguinus sp.PAAbaetetubaEA/EHW
BeH141816Brazil68BSM3, c6/36#1HumanunkPAAbaetetubaEA/EHW
BeAr142658Brazil68C (39)SM2, c6/36#1Haemagogus sp.PABarcarenaEA/EHW
BeAn142027Brazil68D (39)OriginalSaguinus midasPAAbaetetubaEA/EHW
BeH171995Brazil69ASM1 c6/36#1Human +PABelémEA/EHW
BeH203410Brazil71 (39)SM2, c6/36#1HumanunkPAPeixe BoiEA/EHW
BeAr233164Brazil73A (39)Mosq 4Hg. janthinomysGOGoiásWC/TMW
BeAr232869Brazil73B (39)Mosq 1, SM2Haemagogus sp.GOGoiásWC/TMW
BeAr233436Brazil73C (39)OriginalHaemagogus sp.GOBela VistaWC/TMW
BeH233393Brazil73DSM1 c6/36#1Human+GOMara RosaWC/TMW
BeH291597Brazil75SM2, Mosq 1Human +MTAripuanãWC/EHW
BeAr301129Brazil76ASM1 c6/36#1Haemagogus sp.PAAlenquerEA/EHW
BeH301035Brazil76BSM1 c6/36#1Human *PAAlenquerEA/EHW
BeH324213Brazil77SM1 c6/36#1Human *PAAltamiraEA/EHW
BeH350698Brazil78A (38,39)SM2, Mosq 1Human unkPATomé AcuEA/EHW
BeAr350397Brazil78BSM2, Mosq 1Haemagogus sp.PAFaz. Cangalha FormosaEA/EHW
BeH340824Brazil78CSM1 c6/36#1Human +PABelterra-SantarémEA/EHW
BeH371220Brazil79ASM1 c6/36#1Human +PAAltamiraEA/EHW
BeAr378600Brazil80ASM2, C6/36#1Haemagogus sp.GOUruacuWC/TMW
BeH385780Brazil80B (39)SM1 c6/36#1Human +PAAltamiraEA/EHW
BeH379501Brazil80CSM2Human +MAMontes AltosNEW/EHS
BeH394880Brazil81Original, C6/36#1Human *PAConceição AraguaiaEA/EHW
BeH403366Brazil82ASM1 c6/36#1Human +MAAmarante do MaranhãoNEW/EHS
BeH405954Brazil82BSM1 c6/36#1Human *PAAltamiraEA/EHW
BeH 413820Brazil83SM1 c6/36#1Human +ROPorto VelhoWA/EHW
BeH425381Brazil84A (39)c6/36Human +APTribo OyampiNA/EHW
BeAn 424208Brazil84BOriginal, C6/36#1Chiropotes satanasPATucuruíEA/EHW
BeAr424719Brazil84CSM1 c6/36#1Haemagogus sp.PAS. DomingosEA/EHW
BeAr424083Brazil84DSM1 c6/36#1Hg. albomaculatusPAMonte AlegreEA/EHW
BeAr424492Brazil84ESM1 c6/36#1Hg. janthinomysPAFaroEA/EHW
BeH422255Brazil84FSM1 c6/36#1Human +PAFaroEA/EHW
BeH423602Brazil84GSM1 c6/36#1Human +PAS. DomingosEA/EHW
BeH422312Brazil84HSM1Human +PAMonte AlegreEA/EHW
BeH436823Brazil85ASM1 c6/36#1Human +MTSinopWC/EHW
BeAr437159Brazil85BSM1 c6/36#1Hg. janthinomysMTSinopWC/EHW
BeH463676Brazil87SM1 c6/36#1Human +PABrevesEA/EHW
BeH474245Brazil88ASM1 c6/36#1Human +GOCaiapôniaWC/TMW
BeH474297Brazil88BSM1 c6/36#1Human +MGParacatuSE/TSS
BeH485717Brazil89SM1 c6/36#1Human +MGFrancisco DumontSE/TSS
BeAr 511437Brazil91A2SM1, C30c6/36#1Hg. janthinomysPABarcarenaEA/EHW
BeH511843Brazil91BSM1, C6/36#2Human +RRTribo YanomamyWA/EHS
BeAn 510268Brazil91CSM1 c6/36#1Alouatta sp.GOGoiásWC/TMW
BeAr 512943Brazil92A (39)SM1 C6/36-1 Vero2Hg. janthinomysMSSidrolândiaWC/TMW
BeAr513008Brazil92B (39)SM1 c6/36#1Sabethes sp.MSSidrolândiaWC/TMW
BeH512772Brazil92C (39)SM2, C6/36#1Human unkMSCampo GrandeWC/TMW
BeAr513060Brazil92DSM1 c6/36#1Haemagogus sp.MSCampo GrandeWC/TMW
BeAr 513292Brazil92E (39)SM1 C6/36-1 Vero1Sabethes cloropterusMSJaraguaríWC/TMW
BeH520988Brazil93ASM1 c6/36#1Human +MABarra do CordaNEW/EHS
BeH521244Brazil93BSM1 c6/36#1Human *MAMiradorNEW/EHS
BeAr 527785Brazil94A (39)SM1 c6/36#1Sabethes cloropterusMGArinosSE/TSS
BeAr527198Brazil94B (39)SM1 c6/36#1Haemagogus sp.MGArinosSE/TSS
BeAr527547Brazil94CSM2, vero1Haemagogus sp.PADon EliseuEA/EHW
BeH526722Brazil94DSM2, c6/36#1Human +MGArinosSE/TSS
BeAr528057Brazil94EOriginal, C6/36#1Hg. janthinomysMGArinosSE/TSS
BeAr527410Brazil94FOriginal, C6/36#1Hg. janthinomysMAPastos BonsNEW/EHS
BeH535010Brazil95SM2, c6/36#1Human +MAPastos BonsNEW/EHS
BeAr 544276Brazil96A (39)SM1 c6/36#1Hg. janthinomysROCabixiWA/EHW
BeAn 604552Brazil98AOriginal, C6/36#1Alouatta belzebulPAAfuaNA/EHW
BeAr603401Brazil98BOriginal, C6/36#1Hg. janthinomysPAAltamiraEA/EHW
BeAr605158Brazil98CSM1 c6/36#1Hg. janthinomysPAAfuáNA/EHW
BeH603325Brazil98DOriginal, C6/36#1Human unkPAAfuáEA/EHW
BeH605427Brazil98ESM1 c6/36#1Human +PACarajasEA/EHW
BeAr614320Brazil99BSM2Haemagogus sp.PABrevesNA/EHW
BeAr617127Brazil99DSM1 c6/36#1Hg. janthinomysTOMonte do CarmoEA/EHS
BeH613582Brazil99ESM1 c6/36#1Human *PABrevesNA/EHW
BeAr628124Brazil00A (39)SM1 c6/36#1Hg. janthinomysTOParanãEA/EHS
BeH 622491Brazil00BSM1Human +DFBrasiliaWC/THA
BeAn625923Brazil00COriginal, C6/36#1Alouatta sp.GOGoiásWC/TMW
BeH 622205Brazil0DSM1 c6/36#1Human *GOGoiásWC/TMW
BeAr630768Brazil01ASM1Hg. janthinomysGOAlto ParaisoWC/TMW
BeAr631464Brazil01BSM1 c6/36#1Sa. chloropterusBAJaborandiNED/THA
BeAr645693Brazil01DSM1 c6/36#1HaemagogusMGUberabaSE/TSS

aPassage history of seed strain in collection. SM, suckling mouse.
bOutcome of human cases, unk = unknown, + = fatal, * = recovered.
cAP, Amapa; BA, Bahia; DF, Federal District; GO, Goiás; MA, Maranhão; MG, Minas Gerais; MS, Mato Grosso do Sul; MT, Mato Grosso; PA, Pará; RO, Rondônia; RR, Roraima; TO, Tocantins.
dEA, East Amazon; NA, Northern Amazon; NED, Northeast Dry; NEW, Northeast Wet; S, South; SE, Southeast; WA, West Amazon; WC, West-Central.
eEHW, equatorial, hot and wet; EHS, equatorial, hot and semi-wet; THA, tropical, hot and semi-arid; THW, tropical, semi-hot and wet; TMW, tropical, mesothermal, wet; TSW, tropical, semi-hot, wet; TSS, tropical, semi-hot and semi-wet.

Dr. Vasconcelos is a physician and chief of the Arbovirus Department at Evandro Chagas Institute, National Reference Center for Arbovirus, Ministry of Health, Belém, Brazil. His research interests focus on arthropodborne and rodentborne viruses, especially epidemiology, molecular biology, and diagnosis.

ReferencesMonath TP Epidemiology of yellow fever: current status and speculations on future trends. In: Saluzzo JF, Dodet B, editors. Amsterdam: Elsevier; 1997 p. 143–65.Robertson SE, Hull BP, Tomori O, Bele O, LeDuc JW, Esteves K Yellow fever: a decade of reemergence. JAMA. 1996;276:115762 10.1001/jama.1996.035401400450258827969Vainio J, Cutts F Yellow fever. Global programme for vaccines and immunization. WHO/EPI/GEN/98.11. Geneva: World Health Organization; 1998Vasconcelos PFC, Costa ZG, Travassos da Rosa ES, Luna E, Rodrigues SG, Barros VLRS, Epidemic of jungle yellow fever in Brazil, 2000: implications of climatic alterations in disease spread. J Med Virol. 2001;65:598604 10.1002/jmv.207811596099Filippis AM, Schatzmayr HG, Nicolai C, Baran M, Miagostovich MP, Sequeira PC, Jungle yellow fever, Rio de Janeiro. Emerg Infect Dis. 2001;7:484511384541Dégallier N, Travassos da Rosa AP, Hervé JP, Travassos da Rosa JFS, Vasconcelos PFC, Silva CJM, A comparative study of yellow fever in Africa and South America. J Braz Assoc Advanc Sci. 1992;44:14351Mondet B, Travassos da Rosa APA, Vasconcelos PFC The risk of urban yellow fever outbreaks in Brazil by dengue vectors Aedes aegypti and Aedes albopictus. Bull Soc Pathol Exot. 1996;89:107138924767Fundação Nacional de Saúde (FUNASA) Distribuição de casos confirmados de febre amarela notificados por unidade federativa no Brasil, 1981–2002. Brasília: Ministério da Saúde; 2003Massad E, Burattini MN, Coutinho FA, Lopez LF Dengue and the risk of urban yellow fever reintroduction in Sao Paulo State, Brazil. Rev Saude Publica. 2003;37:47784 10.1590/S0034-8910200300040001312937709Figueiredo LT The Brazilian flaviviruses. Microbes Infect. 2000;2:16439 10.1016/S1286-4579(00)01320-411113383Beaty BJ, Calisher CH, Shope RE Arboviruses. In: Lennette EH, Schmidt NJ, editors. Diagnostic procedures for viral, rickettsial, and chlamydial infections. 6th ed. Washington: American Public Health Association; 1989 p. 797–855.Kuno G, Chang GJ, Tsuchiya KR, Karabatsos N, Cropp CB Phylogeny of the genus Flavivirus. J Virol. 1998;72:73839420202Wang E, Weaver SC, Shope RE, Tesh RB, Watts DM, Barrett AD Genetic variation in yellow fever virus: duplication in the 3´ noncoding region of strains from Africa. Virology. 1996;225:27481 10.1006/viro.1996.06018918913Bryant JE, Barrett ADT Comparative phylogenies of yellow fever isolates from Peru and Brazil. FEMS Immunol Med Microbiol. 2003;39:10318 10.1016/S0928-8244(03)00238-414625093Deubel V, Digoutte JP, Monath TP, Girard M Genetic heterogeneity of yellow fever virus strains from Africa and the Americas. J Gen Virol. 1986;67:20913 10.1099/0022-1317-67-1-2093944583Swofford DL PAUP*. Phylogenetic analysis using parsimony (* and other methods). Sunderland (MA): Sinauer Associates; 1998Hahn CS, Dalrymple JM, Strauss JH, Rice CM Comparison of the virulent Asibi strain of yellow fever virus with the 17D vaccine strain derived from it. Proc Natl Acad Sci U S A. 1987;84:201923 10.1073/pnas.84.7.20193470774Duarte dos Santos CN, Post PR, Carvalho R, Ferreira II, Rice CM, Galler R Complete nucleotide sequence of yellow fever virus vaccine strains 17DD and 17D-213. Virus Res. 1995;35:3541 10.1016/0168-1702(94)00076-O7754673Vasconcelos PFC, Luna EJ, Galler R, Silva LJ, Coimbra TL, Barros VLRS, Serious adverse events associated with yellow fever 17DD vaccine in Brazil: a report of two cases [errata in Lancet. 2001;358:336; Lancet. 2001;358:1018] Lancet. 2001;358:917 10.1016/S0140-6736(01)05326-011463409Mutebi JP, Wang H, Li L, Bryant JE, Barrett AD Phylogenetic and evolutionary relationships among yellow fever virus isolates in Africa. J Virol. 2001;75:69997008 10.1128/JVI.75.15.6999-7008.200111435580Rey FA, Heinz FX, Mandl C, Kunz C, Harrison SC The envelope glycoprotein from tick-borne encephalitis virus at 2 angstrom resolution. Nature. 1995;375:2918 10.1038/375291a07753193Allison SL, Schalich J, Stiasny K, Mandl CW, Heinz FX Mutational evidence for an internal fusion peptide in flavivirus envelope protein E. J Virol. 2001;75:426875 10.1128/JVI.75.9.4268-4275.200111287576Fitzgeorge R, Bradish CJ The in vitro differentiation of strains of yellow fever virus in mice. J Gen Virol. 1980;46:113 10.1099/0022-1317-46-1-16766176Monath TP Milestones in the conquest of yellow fever. In: Koprowski H, Oldstone MBA, editors. Microbes and man—then and now. Bloomington (IN): Medi-Ed Press; 1996 p. 95–111.Heraud JM, Hommel D, Hulin A, Deubel V, Poveda JD, Sarthon JL, First case of yellow fever in French Guiana since 1902. Emerg Infect Dis. 1999;5:42932 10.3201/eid0503.99031410341180Soper FL Ventures in world health. Scientific publication no. 355. Washington: Pan American Health Organization; 1977Monath TP Yellow fever: Victor, Victoria? Conqueror, conquest? Epidemics and research in the last 40 years and prospects for the future. Am J Trop Med Hyg. 1991;45:1431867343Pinheiro FP, Travassos da Rosa APA, Moraes MAP, Almeida Neto JC, Camargo S, Filgueiras JP An epidemic of yellow fever in central Brazil, 1972–1973. I. Epidemiological studies. Am J Trop Med Hyg. 1978;27:12532626268Chippaux A, Deubel V, Moreau JP, Reynes JM Current situation of yellow fever in Latin America. Bull Soc Pathol Exot. 1993;86:46047819801Kiple KF, Cooper DB Yellow fever. In: Kiple KF, Graham RR, editors. Cambridge: Cambridge University Press; 1993 p. 1100–7.Monath TP Yellow fever. In: Monath TP, editor. The arboviruses—epidemiology and ecology. Vol. 5. Boca Raton (FL): CRC Press; 1989 p. 139–231.Vasconcelos PFC, Travassos da Rosa APA, Pinheiro FP, Rodrigues SG, Travassos da Rosa ES, Cruz ACR, Aedes aegypti, dengue, and re-urbanization of yellow fever in Brazil and other South American countries: past and present situation and future perspectives. Dengue Bull. 1999;23:5566Vasconcelos PFC, Rodrigues SG, Dégallier N, Moraes MAP, Travassos da Rosa JFS, Travassos da Rosa ES, Epidemic of sylvatic yellow fever in southeast region of Maranhão State, Brazil, 1993–1994. Epidemiological and entomological findings. Am J Trop Med Hyg. 1997;57:13279288803Vasconcelos PFC, Travassos da Rosa APA, Rodrigues SG, Travassos da Rosa ES, Montiero HÃO, Cruz ACR, Yellow fever in Pará State, Amazon region of Brazil, 1998–1999. Entomological and epidemiological findings. Emerg Infect Dis. 2001;7:5659 10.3201/eid0703.01033811485676Saeed MF, Wang H, Nunes M, Vasconcelos PFC, Weaver SC, Shope RE, Nucleotide sequences and phylogeny of the nucleocapsid gene of Oropouche virus. J Gen Virol. 2000;81:743810675412Chan RC, Penney DJ, Little D, Carter IW, Roberts JA, Rawlinson WD Hepatitis and death following vaccination with 17D-204 yellow fever vaccine. Lancet. 2001;358:1212 10.1016/S0140-6736(01)05341-711463415Martin M, Tsai TF, Cropp B, Chang GJJ, Holmes DA, Tseng J, Fever and multisystem organ failure associated with 17D-204 yellow fever vaccination: a report of four cases. Lancet. 2001;358:98104 10.1016/S0140-6736(01)05327-211463410Wang E, Weaver SC, Shope RE, Tesh RB, Watts DM, Barrett AD Genetic variation in yellow fever virus: duplication in the 3´noncoding region of strains from Africa. Virology. 1996;225:27481 10.1006/viro.1996.06018918913Bryant JE, Barrett ADT Comparative phylogenies of yellow fever isolates from Peru and Brazil. FEMS Immunol Med Microbiol.