ORIGINAL RESEARCH article
Distinct YFV Lineages Co-circulated in the Central-Western and Southeastern Brazilian Regions From 2015 to 2018
- 1Laboratório de Genética Molecular de Microorganismos, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil
- 2Laboratório de Mosquitos Transmissores de Hematozoários, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil
- 3Instituto Federal do Norte de Minas Gerais, Salinas, Brazil
- 4Laboratório de Biologia Molecular de Flavivírus, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil
- 5División Biología Molecular y Genética, Departamento de Biología Molecular, Instituto de Investigaciones Biológicas Clemente Estable, Montevideo, Uruguay
- 6Laboratório Central de Saúde Pública Dr. Giovanni Cysneiros, Secretaria de Saúde de Goiás, Goiânia, Brazil
- 7Coordenação Geral de Vigilância das Doenças Transmissíveis, Departamento de Vigilância das Doenças Transmissíveis, Secretaria de Vigilância em Saúde, Ministério da Saúde, Brasília, Brazil
- 8Departamento de Moléstias Infecciosas, Hospital das Clinicas HCFMUSP, Faculdade de Medicina, Universidade de São Paulo, São Paulo, Brazil
- 9Laboratório de AIDS e Imunologia Molecular, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil
The current outbreak of yellow fever virus (YFV) that is afflicting Brazil since the end of 2016 probably originated from a re-introduction of YFV from endemic areas into the non-endemic Southeastern Brazil. However, the lack of genomic sequences from endemic regions hinders the tracking of YFV’s dissemination routes. We assessed the origin and spread of the ongoing YFV Brazilian outbreak analyzing a new set of YFV strains infecting humans, non-human primates (NHPs) and mosquitoes sampled across five Brazilian states from endemic and non-endemic regions between 2015 and 2018. We found two YFV sub-clade 1E lineages circulating in NHP from Goiás state (GO), resulting from independent viral introductions into the Araguaia tributary river basin: while one strain from 2017 clustered intermingled with Venezuelan YFV strains from 2000, the other YFV strains sampled in 2015 and 2017 clustered with sequences of the current YFV outbreak in the Brazilian Southeastern region (named YFV2015-2018 lineage), displaying the same molecular signature associated to the current YFV outbreak. After its introduction in GO at around mid-2014, the YFV2015-2018 lineage followed two paths of dissemination outside GO, originating two major YFV sub-lineages: (1) the YFVMG/ES/RJ sub-lineage spread sequentially from the eastern area of Minas Gerais state to Espírito Santo and then to Rio de Janeiro states, following the Southeast Atlantic basin; (2) the YFVMG/SP sub-lineage spread from the southwestern area of Minas Gerais to the metropolitan region of São Paulo state, following the Paraná basin. These results indicate the ongoing YFV outbreak in Southeastern Brazil originated from a dissemination event from GO almost 2 years before its recognition at the end of 2016. From GO this lineage was introduced in Minas Gerais state at least two times, originating two sub-lineages that followed different routes toward densely populated areas. The spread of YFV outside endemic regions for at least 4 years stresses the imperative importance of the continuous monitoring of YFV to aid decision-making for effective control policies aiming the increase of vaccination coverage to avoid the YFV transmission in densely populated urban centers.
In Brazil, the yellow fever virus (YFV) have been sporadically detected in human and non-human primates (NHPs) populations from the enzootic/endemic Northern (Amazon) and epidemic Central-Western regions during the second half of the 20th century (Carrington and Auguste, 2013; Monath and Vasconcelos, 2015). Since the early 2000s, the virus has progressively expanded to the Southeastern and Southern Brazilian regions and in December 2016 began the largest epizootic/epidemic of sylvatic YF registered in the country over the last 50 years (Vasconcelos, 2010; Possas et al., 2018a,b). Between December 2016 and June 2018, a total of 2,139 YF human cases were confirmed in all Southeastern Brazilian states of Minas Gerais (n = 997), São Paulo (n = 577), Rio de Janeiro (n = 307) and Espírito Santo (n = 258) with 735 deaths (case-fatality, 34%). Moreover, in 2019 YFV transmission continues in São Paulo and is emerging in the north of Paraná state from South region of Brazil (Supplementary Figure 1; Secretaria de Vigilância em Saúde, 2019).
Sequential YF outbreaks reported in the Southeastern and Southern Brazilian regions in 2000–2001, 2008–2009 and 2016–2018 were more likely caused by single independent events of re-introduction of YFV strains from endemic areas (Mir et al., 2017). A recent study speculated that the YFV strain causing the current outbreak would have been originated in the Brazilian Central-West region. This conjecture was based on the date of the most recent common ancestor of the 2016–2018 Brazilian YFV, which was estimated in a period (July 2014 to January 2016) when YFV circulation was reported in the state of Goiás (Central-Western region; Rezende et al., 2018). However, the precise routes of dissemination of YFV strains from endemic to non-endemic areas observed in Brazil in the last 15–20 years are difficult to elucidate because of the scarcity of sequences sampled from endemic regions in that period.
The spatiotemporal dynamics of dissemination of the 2016–2018 Brazilian YFV lineage within the Southeastern region also remained unclear. The first study considering full-genome YFV sequences based on samples from Espírito Santo and Rio de Janeiro states from 2017 placed the origin of the Southeastern outbreak in Espírito Santo in April 2016 (July 2015 to October 2016) and supported a rapid southward viral dissemination in direction to the great metropolitan area of Rio de Janeiro (Gomez et al., 2018). A second study also comprising full-genome YFV sequences sampled from Minas Gerais, Espírito Santo and Rio de Janeiro in 2017 supports that the outbreak arose in Minas Gerais in July 2016 (March to November 2016) and was then southerly disseminated toward Espírito Santo and Rio de Janeiro (Faria et al., 2018). A third study, based on partial genome sequences, showed that YFV strains isolated in the state of São Paulo in 2016 branched in basal position relatively to those isolated in Minas Gerais and Espírito Santo at 2017–2018 and traced the origin of current outbreak to July 2015 (July 2014 to January 2016), thus suggesting that the 2016–2018 YFV Brazilian lineage may have circulated in São Paulo state before spread to Minas Gerais and Espírito Santo (Rezende et al., 2018).
Therefore, to define in greater detail the geographic origin and subsequent dissemination routes of the 2016–2018 Brazilian YFV clade, we generated and analyzed 13 YFV new genomes obtained from humans, NHPs and mosquitoes, between 2015 and 2018, from all Southeast Brazilian states (Minas Gerais, Espírito Santo, Rio de Janeiro e São Paulo) and from Central-West region (Goiás). Interestingly, we identified two YFV lineages circulating in Central-West (Goiás) during 2015–2017, one (GO27/2015 and GO21/2017) displaying the nine amino acid signature characteristic of the 2016–2018 YFV Southeastern outbreak (Bonaldo et al., 2017) and the other (GO05/2017) isolated 2 year later that displayed an amino acid pattern typical of older YFV sub-clade 1E strains sampled in Brazil and Venezuela between 2000 and 2010 (De Souza et al., 2010). Besides, these 13 new YFV complete coding region sequences (CDS) were combined with previously described YFV CDS from Brazil (n = 67), Venezuela (n = 5) and Trinidad and Tobago (n = 1) and then subjected to phylogeographic analyses to get a better picture of the ongoing YFV epidemic in Brazil.
Materials and Methods
Viral samples from 13 infected hosts from distinct biomes and river basins in the states of Goiás (n = 3), Central-Western region, and Rio de Janeiro (n = 2), Minas Gerais (n = 3), Espírito Santo (n = 3) and São Paulo (n = 2), Southeastern region of Brazil, were analyzed (Table 1). Serum samples of human cases, liver samples of human and NHP, and homogenates of entire bodies of pooled adult female mosquitoes were collected and processed as previously described (Ferreira-De-Brito et al., 2016; Bonaldo et al., 2017; Gomez et al., 2018; Abreu et al., 2019). The YFV isolates ABR1005 and ABR1009 were obtained from the sera of a 64 and 30-year-old male patients, respectively. The ABR1005 patient was hospitalized 4 days after onset of symptoms and died due to multiple organ failure. The ABR1009 patient was hospitalized 4 days after onset of symptoms but fully recovered from the infection.
YFV Genome Detection and Nucleotide Sequencing
The sera from the individuals from São Paulo state (ABR1005 and ABR1009) and an NHP liver homogenate sample from Goiás (GO27) were employed to obtain first-passage YFV isolates by infection of monolayer cell cultures of the C6/36 clone of Aedes albopictus. Viral RNA was obtained from cell cultures or directly from samples as described elsewhere (Bonaldo et al., 2017). The set of primers utilized in PCR and sequencing procedures followed a previous report (Gomez et al., 2018). Nucleotide sequences were determined by capillary electrophoresis at the sequencing facility of Fiocruz-RJ (RPT01A – Sequenciamento de DNA – RJ). The sequences were assembled with SeqMan Pro version 8.1.5 (DNASTAR, Madison, WI, United States). The Molecular Evolutionary Genetics Analysis 7.0 program (Kumar et al., 2016) was adopted to explore the amino acid differences as well as to calculate nucleotide and amino acid distances.
Evolutionary and Phylogeographic Analyses
Complete CDS (10,239 nt in length) of the 13 newly generated YFV genomes were combined with CDS of YFV American sequences available in GenBank1 according to the following inclusion criteria: (1) link to a publication, (2) coverage of at least 99% of the viral CDS, and (3) known date and country of collection. The sequences were aligned with MAFFT (Katoh and Standley, 2013) and viral phylogenies were reconstructed by maximum likelihood (ML) analysis implemented in PhyML (Guindon et al., 2010) applying the best substitution model selected by jModelTest v1.6 (Darriba et al., 2012). The temporal signal of different combinations of sequences representing the South American genotypes I and II (SA-I+II), the South American genotype I (SA-I), and the Modern lineage of SA-I (Mir et al., 2017) were examined with TempEst v1.5.1 (Rambaut et al., 2016). The dataset with the best temporal structure (SA-I) was chosen for the subsequent time tree reconstructions by Bayesian method (Supplementary Table 1).
The time scale of the SA-I and 2016–2018 Brazilian YFV datasets were estimated using the Markov chain Monte Carlo algorithms implemented in the BEAST v1.8.4 package (Drummond and Rambaut, 2007; Drummond et al., 2012) with BEAGLE (Ayres et al., 2012) to improve running time. The evolutionary process was estimated using the best-fit nucleotide substitution model (GTR+I+Γ4), a relaxed uncorrelated lognormal molecular clock model (Drummond et al., 2006) and the non-parametric Bayesian Skyline coalescent tree prior (Drummond et al., 2005). A CTMC rate reference prior (Ferreira and Suchard, 2008) and a normal prior (mean = 4.5 × 10-4 substitution/site/year, standard deviation = 1.0 × 10-4) in the evolutionary rate were used for the analysis of the SA-I and 2016–2018 Brazilian YFV datasets, respectively.
The reconstruction of migration events throughout the phylogeny for the 2016–2018 Brazilian YFV lineage also employed the BEAST package using discrete and continuous models. The discrete phylogeographic analysis was performed using reversible (symmetric) and nonreversible (asymmetric) discrete phylogeographic models (Lemey et al., 2009), assigning discrete traits for each sequence representing the Brazilian state of isolation, except for MG, in which further geographic subdivision were employed. The spatiotemporal reconstruction in continuous space utilizing the geographic coordinates (latitude and longitude) of each YFV isolate was estimated with a homogenous Brownian diffusion (BD) model and the heterogeneous Cauchy, Gamma and Lognormal relaxed random walk (RRW) models (Lemey et al., 2010). Comparisons among the different discrete and continuous phylogeographic models were performed using the log marginal likelihood estimation (MLE) based on path sampling (PS) and stepping-stone sampling (SS) methods (Baele et al., 2012). Bayesian analyses were run for 108 generations and convergence (effective sample size > 200) was inspected using TRACER v1.7 (Rambaut et al., 2018) after discarding 10% burn-in. The maximum clade credibility (MCC) trees were summarized using TreeAnnotator v.1.8.4 (Drummond et al., 2012) and visualized with FigTree v.18.104.22.168 The viral spatio-temporal diffusion was analyzed and visualized in SPREAD (Bielejec et al., 2011) and further projected in maps generated with QGIS software3 using public access data collected from the Brazilian Institute of Geography and Statistics (Instituto Brasileiro de Geografia e Estatística [IBGE], 2019) and National Water Agency (Agência Nacional de Águas [ANA], 2019).
Identification of Two YFV Lineages Circulating in the State of Goiás During 2015–2017
Initially, we sequenced the complete genomes of YFV obtained at two Brazilian biomes to infer the origin and dissemination routes of YFV during the 2016–2018 outbreak in the country. Ten YFV strains were sampled from human, NHP and mosquitoes in all Southeastern Brazilian states: Minas Gerais, Espírito Santo, Rio de Janeiro and São Paulo, located in the Atlantic Forest biome and where massive epizootics and human cases occurred. Three other YFV genomes were recovered from human and NHP sampled in Goiás, a state of the Brazilian Central-Western region whose predominant vegetation is the Cerrado, a savanna-like biome occupying the territory between the Amazon and Atlantic rain forest of the Southeast region (Table 1). Both YFV strains from Goiás were sampled from sites located in the Tocantins-Araguaia river basin. The ML and Bayesian phylogenetic analyses placed the newly generated YFV genomes inside the sub-clade 1E (De Souza et al., 2010) of the Modern lineage of SA-I (Mir et al., 2017), with high support [aLRT/posterior probability (PP) = 1] (Figure 1 and Supplementary Figure 2). Two viral samples from Goiás (GO27/2015 and GO21/2017) infecting, respectively a capuchin-monkey from Novo Brasil on August 2015 and a person from Amarinópolis on April 2017, clustered in a highly supported (PP = 1) clade (YFV2015-2018) with all Brazilian YFV sequences from the Southeastern region associated with the current outbreak. By contrast, the other strain from Goiás (GO05/2017) that infected a howler-monkey from Nova Crixás on April 2017, was intermingled among Venezuelan YFV genomic sequences from the 2000s. The mean evolutionary rate for the YFV SA-I was estimated at 4.3 × 10-4 substitution/site/year (Supplementary Table 2), fully consistent with that previously reported (Nunes et al., 2012; Gomez et al., 2018), while the time of the most recent common ancestor (TMRCA) of the YFV2015-2018 clade was estimated on October 2013 [95% Bayesian credible interval (BCI): May 2012 to December 2014].
Figure 1. Bayesian maximum clade credibility phylogeny of the YFV South American genotype I. Only posterior probabilities (PP) of key nodes are denoted. The external node color indicates sample location (geographic region or country) accordingly to the legend. The samples obtained from NHP from Goiás state are indicated by arrows. The clade comprising the GO27/2015 and GO21/2017 samples and all YFV genomes from the ongoing outbreak (YFV2015-2018) is highlighted by a red box, and inner clades were collapsed for clarity. The node representing the most recent common ancestor of YFV2015-2018 lineage is indicated by a red circle along with its temporal origin estimate. All horizontal branch lengths are drawn to a scale of years. Sequences names are coded as accession number/date of collection/host (H – human, V – mosquito vector, M – NHP.
Analysis of amino acid signatures showed the frequent presence of the nine unique amino acid substitutions previously described (Bonaldo et al., 2017; Gomez et al., 2018) in almost all new YFV Brazilian genomes from the Southeastern region and in one sample from Goiás (GO27/2015; Supplementary Table 3 and Supplementary Figure 3). Distinctly, the other YFV strain from Goiás (GO05/2017), exhibited an amino acid pattern more similar to older YFV sub-clade 1E strains sampled in Brazil and Venezuela between 2000 and 2010 (Supplementary Table 3 and Supplementary Figure 3). These results clearly support the circulation of at least two YFV sub-clade 1E lineages in the Cerrado between 2015 and 2017 that probably resulted from independent viral introductions in the Araguaia tributary basin of the Tocantins-Araguaia primary watershed. Interestingly, the third sample from Goiás state (GO21/2017) displayed eight out of nine aminoacidic set of alterations (Supplementary Table 3). Moreover, these results revealed that the molecular signature previously associated with the 2016–2017 YFV Southeastern Brazilian strains was already present in a YFV strain isolated in the Cerrado biome in 2015.
The YFV2015-2018 Lineage Likely Arose in Goiás in 2014 and Was Disseminated to the Southeastern Region Following Two Major Routes
To determine with more precision the geographic origin and dissemination routes of the YFV2015-2018 lineage, we first applied discrete Bayesian symmetric and asymmetric phylogeographic models. The Bayes Factor test showed a strong support in favor of the discrete asymmetric phylogeographic model (Supplementary Table 4), presented in Figure 2. According to this analysis, the YFV2015-2018 lineage likely originated in the state of Goiás [posterior state probability (PSP) = 0.57] in May 2014 (95% BCI: October 2012 to July 2015). From there, it followed two paths of dissemination toward the state of Minas Gerais, originating two major YFV sub-lineages in the Southeastern region here called YFVMG/ES/RJ and YFVMG/SP.
Figure 2. Time-scaled maximum clade credibility phylogeny of the YFV2015-2018 lineage. The branches’ colors represent the most probable location of their descendent nodes as indicated in the legend. The posterior (PP) and posterior state (PSP) probabilities are denoted above and below key branches, respectively. The time for the most common ancestor of lineage YFV2015-2018 is indicated at the basal node (red circle) with the 95% BCI in parenthesis. Tips names were codified as accession number/date/region/host. The tips corresponding to the samples sequenced in this study are indicated by arrows. All horizontal branch lengths are drawn to a scale of years. BA, Bahia; ES, Espírito Santo; GO, Goiás; MG-E, Eastern Minas Gerais; MG-N, Northern Minas Gerais; MG-SE, Southeastern Minas Gerais; MG-SW, Southwest Minas Gerais; RJ, Rio de Janeiro; SP, São Paulo.
The YFVMG/ES/RJ sub-lineage probably reached initially the eastern region of Minas Gerais (MG-E, PSP = 0.51) on October 2015 (95% BCI: February 2015 to April 2016), from where it most likely spread to Espírito Santo state (PSP = 0.99) on March 2016 (95% BCI: September 2015 to August 2016) and to the northern area of Minas Gerais (MG-N, PSP = 0.47) on July 2016 (95% BCI: January 2016 to October 2016). From MG-N, the YFV spread to the south of the state of Bahia, in the Northeastern Brazilian region. From ES, the YFV was introduced in Rio de Janeiro at least four times, generating two successful intrastate transmission chains that advanced southwards: (1) one chain followed the northern side of the Serra do Mar along the Paraíba do Sul tributary basin (Supplementary Figure 4), bypassing the Rio de Janeiro metropolitan region toward the south (Valença municipality) and affecting the southeastern region of Minas Gerais (MG-SE, PSP = 0.57); (2) the other transmission chain spread through the coastal area South of Serra do Mar (Macaé tributary basin), crossing the Rio de Janeiro state metropolitan region and reaching an island located in the southern region of Rio de Janeiro (Ilha Grande).
The YFVMG/SP sub-lineage most likely spread from Goiás into the southwestern area of Minas Gerais (MG-SW, PSP = 0.99) on March 2016 (95% BCI: September 2015 to August 2016) and then disseminated from MG-SW to the metropolitan region of the state of São Paulo. Remarkably, four YFV genomes from MG-E and one from MG-SW displayed differences in the molecular amino acid signature (D15762E; R2607Q). Nevertheless, the variations were not fixed in 2018 YFV samples as identified in the position I2176V clustering in YFVMG/ES/RJ, which were collected in Paraíba tributary basin in 2017 and 2018, suggesting the maintenance of the YFV polymorphism in this region. Overall, these results expanded the geographical and temporal edge of the YFV2015-2018 lineage and further revealed a very low degree of phylogenetic intermixing of YFV2015-2018 strains from different Brazilian states during viral dissemination in the Southeastern region.
The YFV2015-2018 Lineage Was Disseminated Following Primary River Basins in the Southeastern Region
To get some insight regarding the spatiotemporal dynamics of dissemination of the YFV2015-2018 lineage, we applied different continuous phylogeographical models, assuming homogeneous (BD) and heterogeneous (RRW) dispersion rates among lineages. The RRW model with Lognormal distribution was strongly supported as the fittest diffusion model (Supplementary Table 5), indicating significant variation in the diffusion rate among the branches. The phylogeographic continuous model (Figure 3) changed slightly the epicenter of the YFV2015-2018 lineage but supports the existence of two main routes of dissemination within the Brazilian Southeastern region and few viral migrations between different states, consistent with the discrete phylogeographic reconstruction. This analysis also supports that the YFV was disseminated following major river basins.
Figure 3. Reconstructed spatiotemporal diffusion of the YFV2015-2018 lineage. Phylogeny branches were arranged in space according the internal nodes locations inferred by the continuous phylogeographic model. Branches were colored according to time as indicated by the legend. The dark gray lines represent the Brazilian states boundaries while the colored areas the different primary river basins; additional details on tributary basins are available from Instituto Brasileiro de Geografia e Estatística [IBGE] (2019). Brazilian states: DF, Distrito Federal; ES, Espírito Santo; GO, Goiás; MG, Minas Gerais; RJ, Rio de Janeiro; SP, São Paulo.
According to the continuous model, the origin of the YFV2015-2018 lineage would be the central region of Minas Gerais, within the São Francisco watershed, from where it would have independently spread to the west, east and south, reaching the Tocantins-Araguaia basin in the Goiás in 2015, the Southeast Atlantic hydrographic region in Minas Gerais at the end of 2015, and the Paraná hydrographic region of that state in the middle of 2016, respectively. From the eastern of Minas Gerais, the virus moved southward following the Southeast Atlantic watersheds distributed among the states of Minas Gerais, Espírito Santo, and Rio de Janeiro (an area covered by the Atlantic forest biome) and northward, returning to the São Francisco river basin and south of Bahia. Simultaneously, the YFV lineage showed southerly dissemination from the southwestern region of Minas Gerais toward the São Paulo state, following the Paraná basin. We estimate that YFV lineages moved, on average, 0.5 km/day (95% BCI: 0.4 to 0.7 km/day).
The current re-emergence of YFV in the Southeastern Brazilian region resulted in the largest outbreak of sylvatic YF observed in South America in the last decades. The transmission has been expanding southward in Brazil reaching sites considered YFV-free areas for 80 years, and therefore with scanty YFV vaccination coverage. As a result, 38 cases and nine deaths have been reported in January 2019, around 1,160 km from the first signal of increased incidence of YF in the Southeast (north-western region of Minas Gerais) in late 2016 (Secretaria de Vigilância em Saúde, 2019).
Previous studies pointed out that the ongoing YFV outbreak in the Brazilian Southeastern region resulted from a single introduction event of a YFV Modern lineage strain from an endemic area (Mir et al., 2017; Faria et al., 2018; Gomez et al., 2018; Rezende et al., 2018). However, the precise route of viral dissemination was not achieved due to the scarcity of Brazilian YFV sequences sampled from endemic regions over the last years. Here, the analysis of two YFV samples from Goiás (Central-Western region) obtained in 2015 and 2017 revealed that they are phylogenetically related with and carry the same amino acid signature of the YFV strains causing the current outbreak in the Southeastern region. Moreover, the discrete phylogeographic analysis showed that the YFV causing the current Brazilian outbreak probably originated in Goiás at around mid-2014, a result congruent with epidemiological reports on human and NHP infections (Supplementary Figure 1). Altogether, our data stressed the origin of the current YFV outbreak in the Central-western region and expanded the estimated TMRCA of the ongoing YFV outbreak to almost 2 years before it gained epidemiological visibility in the end of 2016 (Secretaria de Vigilância em Saúde, 2017). Interestingly, we also identified for the first time, two YFV lineages circulating in Goiás during 2015–2017, one of which followed two independent paths of dissemination toward Minas Gerais state, originating two major YFV sub-lineages in the Southeast region responsible for the severe ongoing outbreak. Of note, the molecular clock dating approach used to infer the YFV epidemic dynamics had limitations and would be undoubtedly benefited if more genome sequences (sampled from humans and reservoir populations) were available and if multiple scales of YFV evolution (between- and within-host) were considered (Frost et al., 2015).
According to the spatiotemporal epidemiologic reports of YFV infections in both humans and NHPs in the national surveillance system, after 5 years of no records of YFV cases outside the Amazon, 31 NHPs died between May and August 2014 in seven municipalities of the Tocantins state, in the Tocantins-Araguaia basin. The only corpse found still adequate to diagnosis was positive to YFV (Secretaria de Vigilância em Saúde, 2014). The Tocantins-Araguaia basin drains the territory of Tocantins and the great northern part of the neighbor state of Goiás into de Amazon river, and the gallery forests along its tributaries consist of a large network of corridors between the Amazon and Cerrado biomes. Thus, soon in early 2015, besides in the Tocantins, YFV infections were detected in Goiás, and the virus spilled over from the Tocantins-Araguaia into the São Francisco and Paraná basins, with reports scattered in northwestern Minas Gerais and southern Goiás and the Federal District (Brasília; Figure 4).
Figure 4. Municipalities reporting YFV in Brazil from late 2014 to late 2018 according to epidemiological data. Colors distinguish the semester when the first confirmed human and/or non-human YFV infection was recorded in such a municipality on a half-yearly basis during the period. The dark gray lines represent the Brazilian states boundaries while the colored areas the different primary river basins; additional details on tributary basins are available from Instituto Brasileiro de Geografia e Estatística [IBGE] (2019). Brazilian states: AC, Acre; AL, Alagoas; AP, Amapá; AM, Amazonas; BA, Bahia; CE, Ceará; DF, Distrito Federal; ES, Espírito Santo; GO, Goiás; MA, Maranhão; MT, Mato Grosso; MS, Mato Grosso do Sul; MG, Minas Gerais; PA, Pará; PB, Paraíba; PR, Paraná; PE, Pernambuco; PI, Piauí; RJ, Rio de Janeiro; RN, Rio Grande do Norte; RS, Rio Grande do Sul; RO, Rondônia; RR, Roraima; SC, Santa Catarina; SP, São Paulo; SE, Sergipe; TO, Tocantins.
Brazilian Health authorities consider this YFV reemergence as the start of the north-southeast virus spreading wave that has not stopped yet (Secretaria de Vigilância em Saúde, 2015, 2017, 2019). Accordingly, in late 2015 and early 2016, the circulation of YFV was detected in sites of Minas Gerais and Goiás states, but still essentially limited to the Cerrado biome in the three above mentioned river basins. Intriguingly, in late 2016, human cases and epizooties due to YFV were recorded simultaneously and independently in the northeast Minas Gerais as well as in southeastern Minas Gerais/northern São Paulo, where the virus has then gained the Atlantic forest ecosystem. In fact, the outbreak was officially recognized only when the epidemiological records rapidly peaked and as that the virus reached the Southeast Atlantic river basin in Espírito Santo and Rio de Janeiro, as well as continued to spread in all primary watershed of São Paulo and Minas Gerais throughout 2017. One of the causes of this delay was the introduction of Chikungunya (2014) and Zika (2015) virus in Brazil and the resulting epidemics, which probably reduced the sensitivity of surveillance and interfered with the visibility of the YFV reemergence. In 2017, scattered records of YFV circulation were also made in the states of Bahia, Goiás, Mato Grosso and even in the Amazon.
In early 2018, epidemiological data suggested that transmission was mostly concentrated in the Southeast Atlantic basin in the states of Rio de Janeiro, Minas Gerais and São Paulo and eastern Paraná basin (Figure 4). Then, the place and time of origin for the current YFV Brazilian outbreak here estimated fully agree with the above mentioned officially confirmed infections in NHP and humans reported by the Brazilian Ministry of Health (Secretaria de Vigilância em Saúde, 2015). Indeed, the Novo Brasil municipality in Goiás, from where the NHP YFV sample analyzed here was taken, is located about 100 km from the banks of the Araguaia river. The earlier detection of YFV in NHP from the Tocantins indicates that the state of Goiás probably acted as a staging post-during the dissemination of the YFV2015-2018 lineage from the North (Amazon basin) to the Southeast Brazilian regions. In this sense, the Tocantins-Araguaia watersheds may have played a major role in YFV dissemination as it extends from more than 2.500 km following two river axes (Araguaia and the Tocantins), offering a contiguous connection between the Cerrado (southward, Goiás and Tocantins states) and Amazon (northward, Pará state) biomes (Agência Nacional de Águas [ANA], 2019).
Both discrete and continuous phylogeographic models combined with the epidemiological records support of the YFV2015-2018 lineage moved toward densely populated Southeastern urban regions with low YFV vaccine coverages following major routes along different primary river basins. The phylogeographic analyses pointed out that the YFV2015-2018 lineage probably arrived in the Southeast region via the São Francisco watershed located in Minas Gerais and then moved to the Southeast Atlantic watersheds in the east and the Paraná hydrographic region in the southwest. The viral lineage that moved following the Southeast Atlantic watersheds reached the eastern and northern areas of Minas Gerais state, as well as the south of Bahia, Espírito Santo and Rio de Janeiro states. The viral lineage that followed the Paraná hydrographic region spread to the Southwest of Minas Gerais and São Paulo states. A previous study proposed that the YFV2015-2018 lineage was introduced in the southeastern region through São Paulo and then moved to other Southeastern states (Rezende et al., 2018). However, this conclusion is hampered as they used only partial genomes (with a very low number of nucleotide substitutions supporting the phylogenetic relationships) and did not conduct a formal phylogeographic analysis. Although Faria et al. (2018) already described that the YFV2015-2018 lineage displayed southward and eastward expansion from its inferred origin in Minas Gerais, our results consist of the first description of concurrent dispersion of the YFV2015-2018 lineage following two independent routes that seems to be linked to the main hydrographical basins.
Our phylogeographic analyses also support that the rapid spread of the YFV2015-2018 lineage in the Southeastern region seems to have resulted from a few successful viral disseminations events between states. Most YFV transmission in Espírito Santo was probably originated from a single successfully transmission from Minas Gerais, while most viral transmissions in Rio de Janeiro seems to have resulted from two independent introductions from Espírito Santo that subsequently spread along the coastal and northern sides of the Serra do Mar mountain system, as previously described (Gomez et al., 2018). We found that both transmission chains previously detected in Rio de Janeiro state continued to expand to beyond the metropolitan region, reaching municipalities close to the border with the state of São Paulo during 2018. The two YFV 2018 genomes from São Paulo analyzed here are the first described of the current YFV outbreak from that state and were the result of independent dissemination from the southwestern region of Minas Gerais, but more sequences from São Paulo are necessary to understand the epidemic dynamic in this state. It is unclear if most YFV infections in São Paulo resulted from a single or a few founder viral strains that spread from the southwest of Minas Gerais along the Paraná hydrographic basin, or if other viral strains may have also been disseminated from Rio de Janeiro along the Southeast Atlantic watersheds. The recent detection of YFV in NHP from the coastal area of Paraná state in 2019 (Secretaria de Vigilância em Saúde, 2019) indicates continuous dissemination of YFV into the Southern region probably following the Paraná and/or the Southeast Atlantic hydrographic basins.
We estimated that the YFV2015-2018 lineage moved on average 0.5 km/day and similar results were obtained when the YFV2015-2018’s outgroup sequences GO27/2015 and GO21/2017 were removed from the analysis. This velocity is lower than the estimates described by Faria et al. (2018) and Gomez et al. (2018) that also analyzed the dispersion of the ongoing YFV outbreak in Brazil and found dispersion rates of 4.2 and 3.4 km/day, respectively. However, those studies analyzed sequences sampled between January-April of 2017, corresponding to the wet and warmer season, when there is an increase in the density of vectors (Alencar et al., 2018) facilitating the transmission. The primary vectors in the current outbreak in Southeast Brazil are the mosquitoes Haemagogus leucocelaenus and Haemagogus janthinomys (Abreu et al., 2019), which can disperse large distances in short time (Causey et al., 1950). Our conservative velocity of YFV dispersion fully agrees with the observed distances traveled by mosquitoes, but also with howler monkeys in Southeastern Brazil (Jung et al., 2015) and also agrees with an epidemiological model based on dates and place of reported monkey deaths, which estimated YFV displacement speeds of 2.7 km/day in the warmer months and 0.5 km/day in the coldest months (Fioravanti, 2018). Thus, the YFV dispersion velocity estimated in this study would correspond to a median value between these two speeds. Besides the role of NHPs and mosquitoes, the hypothetical potential of viremic humans in the spread of sylvatic yellow fever has yet to be confirmed (Vasconcelos, 2010; Possas et al., 2018a).
Surprisingly, we found two YFV sub-clade 1E lineages circulating in the Araguaia tributary basin, indicating at least two independent introductions of YFV in that region probably from the enzootic/endemic Amazon biome in a narrow time frame. While sample GO27/2015 isolated in 2015 displays the nine unique amino acid signatures characteristic of the 2016–2018 YFV Southeastern outbreak, sample GO05/2017 isolated 2 years later at the same watershed exhibited an amino acid pattern typical of older YFV sub-clade 1E strains sampled in Brazil and Venezuela between 2000 and 2010. According to our analysis, only YFV strains related to the GO27/2015 and GO21/2017 were able to further disseminate from Goiás toward the Southeast region. We can speculate that the different pattern of molecular signatures present in the two YFV strains from Goiás modulate the spread of each viral lineage since some of them were located in key viral enzymes (Gomez et al., 2018). Alternatively, the GO05/2017 lineage may have been introduced from the Amazon into Goiás at a later time, and its dissemination toward the Southeast was hampered due to the reduction or even exhaustion of susceptible NHP hosts caused by the previous passage of the lineage that originated the YFV2015-2018 clade. Consistent with this last hypothesis, a recent ecological study concludes that dissemination of YFV in South America is not random, but it is influenced by key geo-environmental factors like diversity and number of susceptible NHP hosts (Hamrick et al., 2017). Curiously, epidemiological records showed YFV transmission in several sites in the Amazon in early 2017, including in the Tocantins-Araguaia basin in Pará state (Figure 4).
In summary, we showed that at least two different YFV lineages circulated in the Cerrado biome (Araguaia tributary basin) in a narrow time frame. One of these lineages further spread out the Cerrado biome in Goiás to the Atlantic forest biome in the Southeastern Brazilian region, originating the current Brazilian outbreak (YFV2015-2018 lineage) at around mid-2014. The ongoing YFV outbreak in Brazil disseminated in the Southeast region following two independent routes that seem to be linked to the Paraná and Southeast Atlantic hydrographic basins, comprising densely populated regions. The spread of the YFV outside the Amazon and Cerrado biomes following primary hydrographic watersheds comprising large metropolis stresses the imperative importance of the continuous monitoring of YFV coupled with in-depth phylogeographic analysis to aid decision-making Health authorities for effective prophylactic and control policies aiming the increase of vaccination coverage to avoid the YFV transmission in densely populated urban centers.
All new YFV genomes generated for this study were submitted to GenBank and their accession numbers are included in the manuscript. The information of all YFV genomes used in this study is provided in the Supplementary Table 1.
The YFV isolates ABR1005 and ABR1009 were obtained from two patients following a study protocol approved by the institutional review boards at the Hospital das Clínicas (School of Medicine, University of São Paulo) and the Infectiology Institute Emilio Ribas (CAAE: 59542216.3.1001.0068). The analysis of human samples was carried out in accordance with the recommendations of the Ethics Committees for human research at Instituto Oswaldo Cruz (CAAE 69206217.8.0000.5248), which exempted the need of a specific written informed consent from patients or their legal representatives. Capture of NHPs and mosquitoes, as well as management of NHP samples, were carried out in accordance with the Brazilian environmental authorities (SISBIO-MMA licenses 54707-6 and 52472-2, and INEA licenses 012/2016, 019/2018) and Ethics Committee of Animal use at Instituto Oswaldo Cruz (CEUA license L037/2016).
ED, GB, MB, and RL-d-O conceived the study. FdA and MN carried out the collection of biological specimens. EK provided the samples. IB and MN identified the mosquito species. AF-d-B, FdA, IB, LdS, MdC, and RdM carried out viral RNA extraction from the biological specimens and the diagnosis by RT-PCR. LdS and LR conducted the inoculation of biological specimens in cell culture. AdS, IR, and MG performed rapid viral RNA extraction and genome sequencing. AdS, MB, MG, and NF analyzed the genome sequences. ED and GB performed phylogenetic and phylogeographic analysis. ED, GB, MB, MG, NF, and RL-d-O prepared figures, tables, and supplementary material. ED, FdA, GB, IR, MB, MG, RL-d-O, EK, and AV prepared the manuscript. AR, DR, and FdA gathered, systematized, and illustrated epidemiological records. All authors critically read and approved the final version of the manuscript.
ED was financed by a Postdoctoral fellowship from the “Programa Nacional de Pós-Doutorado (PNPD)” by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior – Brazil (CAPES) – Finance Code 001. MG and IR received a Postdoctoral fellowship from the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior – Brasil (CAPES) – Finance Code 001. RL-d-O is funded by grants from Conselho Nacional Desenvolvimento Científico e Tecnológico (CNPq) (Grant Nos. 309577/2013-6 and 312446/2018-7), Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro (Grant E-26/203.064/2016), Institut Pasteur, Transversal Research Program (PTR Grant No. 528) and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (Grant No. COFECUB 799-14, AUXPE 1731/2014). MB is a recipient of CNPq fellowship for Productivity in Technological Development and Innovative Extension (Grant 309471/2016-8) and is funded by grants from Preventing and Combating the Zika Virus, MCTIC/FNDCT -CNPq/MEC-CAPES/MS-Decit. (Grants 426767/2018-7 and 88881.130684/2016-01) and INOVA-Fiocruz (Grant VPPIS-004-FIO18).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
To Vinicius Lemes da Silva, Myriam A. F Campos, Angélica Bastos, Carmen Helena, Ramos, Yulla Fernandes, Marcelo Santalucia (Secretaria de Saúde de Goiás), Romário Gabriel Aquino (Environmental Surveillance, Secretaria Municipal de Saúde de Angra dos Reis), Rodrigo F. C. Said (Secretria de Saúde de Minas Gerais), Cláudia Aarestrup, Livia Passarela, Adalberto Mitterofhe, Milton F. Castro, Adilson C. Lima (Secretarias Municipais de Saúde de Juiz de Fora e Belmiro Braga), Omar Figueiredo Neto (Environmental Surveillance, Secretaria Municipal de Saúde de Valença), Marilza L. Lange, Luciano L. Salles, Núcleo de Entomologia e Malacologia do Espírito Santo (NEMES), Tercius Barrada (Parque Estadual da Ilha Grande), Gilsa Aparecida P. Rodrigues (Secretaria de Saúde do Estado do Espírito Santo), Gilton Luiz Almada (Centro de Informação Estratégica de Vigilância em Saúde-ES), Mário Sérgio Ribeiro e Patrícia Menegueti (Secretaria de Saúde do Estado do Rio de Janeiro), Prefeitura de Macaé, Centro de Estudos Ambientais e Desenvolvimento Sustentável (CEADS/UERJ), Marcelo Celestino dos Santos, Mauro M. Muniz, Marcelo Q. Gomes, Teresa F. Silva-do-Nascimento for the help in obtaining the viral samples and support in the field work.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.01079/full#supplementary-material
Abreu, F. V. S. D., Ribeiro, I. P., Ferreira-De-Brito, A., Santos, A. A. C. D., Miranda, R. M. D., Bonelly, I. D. S., et al. (2019). Haemagogus leucocelaenus and Haemagogus janthinomys are the primary vectors in the major yellow fever outbreak in Brazil, 2016–2018. Emerg. Microb. Infect. 8, 218–231.
Agência Nacional de Águas [ANA] (2019). Divisões hidrográficas do Brasil. Available at: http://www3.ana.gov.br/portal/ANA/aguas-no-brasil/panorama-das-aguas/copy_of_divisoes-hidrograficas (accessed January 31, 2019).
Alencar, J., Mello, C. F. D., Morone, F., Albuquerque, H. G., Serra-Freire, N. M., Gleiser, R. M., et al. (2018). Distribution of haemagogus and sabethes species in relation to forest cover and climatic factors in the chapada dos guimarães National Park, State of Mato Grosso, Brazil. J. Am. Mosq. Control Assoc. 34, 85–92.
Ayres, D. L., Darling, A., Zwickl, D. J., Beerli, P., Holder, M. T., Lewis, P. O., et al. (2012). BEAGLE: an application programming interface and high-performance computing library for statistical phylogenetics. Syst. Biol. 61, 170–173. doi: 10.1093/sysbio/syr100
Baele, G., Lemey, P., Bedford, T., Rambaut, A., Suchard, M. A., and Alekseyenko, A. V. (2012). Improving the accuracy of demographic and molecular clock model comparison while accommodating phylogenetic uncertainty. Mol. Biol. Evol. 29, 2157–2167. doi: 10.1093/molbev/mss084
Bonaldo, M. C., Gomez, M. M., Dos Santos, A. A., Abreu, F. V. S., Ferreira-De-Brito, A., Miranda, R. M., et al. (2017). Genome analysis of yellow fever virus of the ongoing outbreak in Brazil reveals polymorphisms. Mem. Inst. Oswaldo Cruz 112, 447–451. doi: 10.1590/0074-02760170134
Carrington, C. V., and Auguste, A. J. (2013). Evolutionary and ecological factors underlying the tempo and distribution of yellow fever virus activity. Infect. Genet. Evol. 13, 198–210. doi: 10.1016/j.meegid.2012.08.015
De Souza, R. P., Foster, P. G., Sallum, M. A., Coimbra, T. L., Maeda, A. Y., Silveira, V. R., et al. (2010). Detection of a new yellow fever virus lineage within the South American genotype I in Brazil. J. Med. Virol. 82, 175–185. doi: 10.1002/jmv.21606
Faria, N. R., Kraemer, M. U. G., Hill, S. C., Goes De Jesus, J., Aguiar, R. S., Iani, F. C. M., et al. (2018). Genomic and epidemiological monitoring of yellow fever virus transmission potential. Science 361, 894–899. doi: 10.1126/science.aat7115
Ferreira-De-Brito, A., Ribeiro, I. P., Miranda, R. M., Fernandes, R. S., Campos, S. S., Silva, K. A., et al. (2016). First detection of natural infection of Aedes aegypti with Zika virus in Brazil and throughout South America. Mem. Inst. Oswaldo Cruz 111, 655–658. doi: 10.1590/0074-02760160332
Fioravanti, C. (2018). O Alarme Dos Macacos. Available at: http://revistapesquisa.fapesp.br/2018/01/11/o-alarme-dos-macacos/ (accessed January 31, 2019).
Gomez, M. M., Abreu, F. V. S., Santos, A., Mello, I. S., Santos, M. P., Ribeiro, I. P., et al. (2018). Genomic and structural features of the yellow fever virus from the 2016-2017 Brazilian outbreak. J. Gen. Virol. 99, 536–548. doi: 10.1099/jgv.0.001033
Guindon, S., Dufayard, J. F., Lefort, V., Anisimova, M., Hordijk, W., and Gascuel, O. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307–321. doi: 10.1093/sysbio/syq010
Hamrick, P. N., Aldighieri, S., Machado, G., Leonel, D. G., Vilca, L. M., Uriona, S., et al. (2017). Geographic patterns and environmental factors associated with human yellow fever presence in the Americas. PLoS Negl. Trop. Dis. 11:e0005897. doi: 10.1371/journal.pntd.0005897
Instituto Brasileiro de Geografia e Estatística [IBGE] (2019). Portal De mapas. Available at: https://mapas.ibge.gov.br/images/pdf/mapas/mappag99.pdf (Accessed January 31, 2019).
Jung, L., Mourthe, I., Grelle, C. E., Strier, K. B., and Boubli, J. P. (2015). Effects of local habitat variation on the behavioral ecology of two sympatric groups of brown howler monkey (Alouatta clamitans). PLoS One 10:e0129789. doi: 10.1371/journal.pone.0129789
Mir, D., Delatorre, E., Bonaldo, M., Lourenco-De-Oliveira, R., Vicente, A. C., and Bello, G. (2017). Phylodynamics of yellow fever virus in the Americas: new insights into the origin of the 2017 Brazilian outbreak. Sci. Rep. 7:7385. doi: 10.1038/s41598-017-07873-7
Nunes, M. R., Palacios, G., Cardoso, J. F., Martins, L. C., Sousa, E. C Jr., de Lima, C. P., et al. (2012). Genomic and phylogenetic characterization of Brazilian yellow fever virus strains. J. Virol. 86, 13263–13271. doi: 10.1128/JVI.00565-12
Possas, C., Lourenco-De-Oliveira, R., Tauil, P. L., Pinheiro, F. P., Pissinatti, A., Cunha, R. V. D., et al. (2018a). Yellow fever outbreak in Brazil: the puzzle of rapid viral spread and challenges for immunisation. Mem. Inst. Oswaldo Cruz 113:e180278. doi: 10.1590/0074-02760180278
Rezende, I. M., Sacchetto, L., Munhoz De Mello, E., Alves, P. A., Iani, F. C. M., Adelino, T. E. R., et al. (2018). Persistence of Yellow fever virus outside the Amazon Basin, causing epidemics in Southeast Brazil, from 2016 to 2018. PLoS Negl. Trop. Dis. 12:e0006538. doi: 10.1371/journal.pntd.0006538
Secretaria de Vigilância em Saúde (2015). Reemergência da Febre Amarela Silvestre no Brasil, 2014/2015: Situação Epidemiológica e a Importância da Vacinação Preventiva e da Vigilância Intensificada no Período Sazonal. Brazil: Ministério da Saúde.
Keywords: yellow fever virus, Brazilian outbreak, lineages, amino acid changes, phylogeography
Citation: Delatorre E, Abreu FVS, Ribeiro IP, Gómez MM, dos Santos AAC, Ferreira-de-Brito A, Neves MSAS, Bonelly I, de Miranda RM, Furtado ND, Raphael LMS, da Silva LdFF, de Castro MG, Ramos DG, Romano APM, Kallás EG, Vicente ACP, Bello G, Lourenço-de-Oliveira R and Bonaldo MC (2019) Distinct YFV Lineages Co-circulated in the Central-Western and Southeastern Brazilian Regions From 2015 to 2018. Front. Microbiol. 10:1079. doi: 10.3389/fmicb.2019.01079
Received: 20 February 2019; Accepted: 29 April 2019;
Published: 24 May 2019.
Edited by:Akio Adachi, Kansai Medical University, Japan
Reviewed by:Mike Holbrook, Battelle, United States
Tamas Bakonyi, Szent István University, Hungary
Nikos Vasilakis, The University of Texas Medical Branch at Galveston, United States
Copyright © 2019 Delatorre, Abreu, Ribeiro, Gómez, dos Santos, Ferreira-de-Brito, Neves, Bonelly, de Miranda, Furtado, Raphael, da Silva, de Castro, Ramos, Romano, Kallás, Vicente, Bello, Lourenço-de-Oliveira and Bonaldo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†Joint first authors
‡Joint last authors