Skip to main content
Advertisement
  • Loading metrics

Mapping environmental suitability of Haemagogus and Sabethes spp. mosquitoes to understand sylvatic transmission risk of yellow fever virus in Brazil

Abstract

Background

Yellow fever (YF) is an arboviral disease which is endemic to Brazil due to a sylvatic transmission cycle maintained by infected mosquito vectors, non-human primate (NHP) hosts, and humans. Despite the existence of an effective vaccine, recent sporadic YF epidemics have underscored concerns about sylvatic vector surveillance, as very little is known about their spatial distribution. Here, we model and map the environmental suitability of YF’s main vectors in Brazil, Haemagogus spp. and Sabethes spp., and use human population and NHP data to identify locations prone to transmission and spillover risk.

Methodology/Principal findings

We compiled a comprehensive set of occurrence records on Hg. janthinomys, Hg. leucocelaenus, and Sabethes spp. from 1991–2019 using primary and secondary data sources. Linking these data with selected environmental and land-cover variables, we adopted a stacked regression ensemble modelling approach (elastic-net regularized GLM, extreme gradient boosted regression trees, and random forest) to predict the environmental suitability of these species across Brazil at a 1 km x 1 km resolution. We show that while suitability for each species varies spatially, high suitability for all species was predicted in the Southeastern region where recent outbreaks have occurred. By integrating data on NHP host reservoirs and human populations, our risk maps further highlight municipalities within the region that are prone to transmission and spillover.

Conclusions/Significance

Our maps of sylvatic vector suitability can help elucidate potential locations of sylvatic reservoirs and be used as a tool to help mitigate risk of future YF outbreaks and assist in vector surveillance. Furthermore, at-risk regions identified from our work could help disease control and elucidate gaps in vaccination coverage and NHP host surveillance.

Author summary

Yellow fever virus (YFV) is an arbovirus transmitted to humans from mosquitoes and can lead to severe disease and death. Recent sporadic outbreaks coupled with low vaccination coverage have highlighted the importance of mosquito surveillance for preventing future outbreaks and potential virus spillover into dense urban areas. Yet, very little is known about the spatial distribution of mosquitoes known to transmit YFV and the factors that contribute to their environmental suitability in Brazil. We compiled an occurrence database of primary and secondary mosquito vectors belonging to Haemagogus and Sabethes species’ collected between 1991–2019 and integrated this data with environmental and land-use data to predict their spatial suitability at 1 km x 1 km resolution. Using this information, we identified suitable regions for their co-existence. We overlaid this information with human population density and locations of non-human primate host reservoirs to identify areas at risk of transmission and spillover. Our study provides high-resolution mapping tools to assist with mosquito and arbovirus surveillance which is especially useful in low-resource settings.

Introduction

In Brazil, several sylvatic mosquito species’ belonging to the genera Haemagogus and Sabethes are the primary and secondary vectors of yellow fever virus (YFV) [1,2]. YFV is an arbovirus belonging to the family Flaviviridae that causes acute infectious disease in humans. Literature states that among non-vaccinated populations from Africa and South America, about 55% (95% credible interval, CI: 37%-74%) of individuals infected with YFV are asymptomatic, while about 33% (95% CI: 13%-52%) develop mild disease. Among the approximate 12% (95% CI: 5%-26%) that develop severe disease, case fatality can reach 47% or greater (95% CI: 31%-62%) [3]. In peri-urban and urban areas, transmission of YFV can also be maintained between Aedes aegypti (Linnaeus, 1972) and humans; however, this transmission cycle has not been documented in Brazil since 1942 [4]. YFV is endemic in the Brazilian Amazon Basin and mainly remained there until the 1940s, when sporadic outbreaks were observed in the Northeast, Southeast, and Southern states [46]. Recent re-emergence of YFV in Southeastern Brazil between 2016 and 2018 caused the worst YF outbreaks in the last 70 years [7].

Despite the existence of an effective YF vaccine since 1937, eradication efforts in Brazil have been undermined by an enzootic sylvatic transmission cycle maintained in rural areas by reservoirs of Haemagogus spp. and Sabethes spp. belonging to the Culicidae family. These species feed on non-human primates (NHP) (e.g. Alouatta sp.) that dwell in high canopy forests, and occasionally on humans that live in or travel to these areas. During the 2016–2018 outbreaks, Hg. janthinomys, which can be naturally infected with YFV [8], was found to be the primary vector, while infected Hg. leucocelaenus were also abundant in outbreak areas [1]. Sa. chloropterus was found to play a local or secondary role in sylvatic transmission due to limited abundance and distribution [1]. These mosquitoes are vertically dispersed in forest areas [9] and naturally prefer to bite and breed at high canopy levels [10,11]. Due to environmental disturbances, contact between sylvatic mosquitoes and human populations have increased [12]. While the occurrence of YF cases follows a seasonal pattern, during which most cases are reported between December to May, outbreaks are still highly dependent on environmental conditions that facilitate transmission [7]. Given the frequency and severity of recent epidemics and the ubiquity of human travel to sylvatic transmission zones, there is increasing concern for a sylvatic-to-urban cycle spillover into Ae. aegypti infested areas with low vaccination coverage, especially those in near proximity to sylvatic environments. Thus, a greater understanding of the distribution of sylvatic vectors of YFV and its implications for a spillover is needed to help prevent future epidemics.

Species distribution modelling (SDM) techniques are useful for mapping the distribution of arbovirus vectors by focusing on their ecology as predicted from a combination of data on climate, environmental, and population factors [13,14]. To date, only two studies have attempted to map the distribution of sylvatic vectors of YFV. Childs et al. [15] applied SDM using a MaxEnt model to determine the combined distribution of Hg. janthinomys, Hg. leucocelaenus, and Sa. chloropterus in Brazil, by integrating data on climatic and environmental factors and occurrence data retrieved from the Global Biodiversity Information Facility (GBIF) and published literature. Similarly, Almeida et al. [16] predicted the distribution of Hg. leucocelaenus in the state of Rio Grande do Sul using occurrence data from the 2008–2009 epizootics.

However, these results alone may not accurately determine the ecological niche of each vector and their locations given use of confined occurrence data and limited presumptions about the environmental conditions that influence their distribution in Brazil. Random sampling assumptions adopted for identifying their potential absence via the MaxEnt model may reflect unbiased surveillance efforts, which is rarely the case in the field. While Childs et al. [15] adopted a pseudo-absence selection method consisting of subsampling from occurrence data on mosquito species other than Hg. janthinomys, Hg. leucocelaenus, and Sa. chloropterus, this method may incorrectly identify potential pseudo-absence locations as it classifies excluded species of Haemagogus and Sabethes sharing similar ecology with the aforementioned primary vectors as pseudo-absences. Given that little is known about these vectors and their niches, these limitations should be carefully considered. Moreover, neither study assessed the environmental suitability of Hg. janthinomys, Hg. leucocelaenus, and Sabethes spp. individually across Brazil, which is important given differences in their ecological niche and distributions across space.

Building on these previous efforts, we modelled the distribution of Hg. janthinomys, Hg. leucocelaenus, and Sabethes spp. in Brazil by compiling the most comprehensive dataset of Haemagogus and Sabethes spp. to date, which contains not only occurrence data but also records from GBIF database, published literature, and primary entomological surveys collected between 1991 and 2019. Using covariate data representative of environmental conditions (daytime land surface temperature, precipitation, relative humidity, and land-cover classes) and human activity (population density and deforestation) conducive to sylvatic mosquito activity (as confirmed by lab and field studies), we modelled the probability of environmental suitability for each species using a stacked regression ensemble approach consisting of three different models to increase predictive performance. By combining the resulting predicted species suitability from this study with modelled presence probability of non-human primate host reservoirs, within and surrounding the Southeastern Atlantic forest biome [17], we further mapped the risk of YFV transmission and spillover in highly populated areas of the region. These maps could inform policy makers the spatial distribution of current and emerging YFV risk zones, which can be targeted for transmission prevention and vaccination programmes.

Materials and methods

We adopted a stacked binomial regression ensemble modelling approach to quantify the probability of environmental suitability of Hg. janthinomys, Hg. leucocelaenus, and Sabethes spp., respectively. A stacked regression ensemble is a diverse combination of base models (level-0), with potentially different assumptions about the predictions, into a meta-model (level-1), which has equal or better predictive performance than each individual model [18]. The meta-model is trained on predictions made by the base models and is adjusted based on an optimal weighted average of these models. Previously, this ensemble approach was adopted to map Plasmodium falciparum malaria parasite prevalence in sub-Saharan Africa, and was found to outperform any individual method [19]. Here we used three models in an ensemble: a lasso or elastic-net regularized generalized linear model (GLM), an extreme gradient boosted regression tree model, and a random forest model. We chose these models because they were suitable for binary prediction and have relatively good generalizability for disease mapping [20,21].

We used the following set of input variables to run the ensemble: (i) a comprehensive dataset of geocoded occurrence locations for Hg. janthinomys, Hg. leucocelaenus, and Sabethes spp. recorded from 1991–2019; (ii) selected environmental and land-cover covariate datasets that help to explain the ecological niche of each species; and (iii) a dataset outlining the geocoded locations of pseudo-absence records which refines the geographical range of species occurrence and reduces sampling bias [21]. This is described in more detail in the following sections and in S1 and S2 Figs and S1 Table, alongside full details of the data and each analysis stage. An outline of the methodology from data assembly, modelling approach, to outputs is described in Fig 1.

thumbnail
Fig 1. Schematic diagram of the modelling approach adopted to produce maps of environmental suitability for Haemagogus spp. and Sabethes spp. at 1 km x 1 km resolution.

The base layer of the map was retrieved from https://www.ibge.gov.br/geociencias/downloads-geociencias.html.

https://doi.org/10.1371/journal.pntd.0010019.g001

Presence records

To assess contemporary distributions, we extracted all records of Hg. janthinomys, Hg. leucocelaenus, and Sabethes spp. between the years 1991–2019 from our compiled dataset, which contains the geographical coordinates of species occurrence. Due to sparse records, we extracted all species of Sabethes spp. that have been confirmed of transmitting, capable of being naturally infected, or are found before and during previous and recent YFV outbreaks in Brazil. This includes Sa. chloropterus, Sa. albiprivus, Sa. soperi, Sa. glaucodaemon, Sa. cyaneus, Sa. purpureus, Sa. melanonymphe, Sa. undosus/fabrici/ignotus, Sa. intermedius, Sa. belisarioi, Sa. aurescens, Sa. identicus, Sa. xyphydes, and Sabethes sp. [1,2,8,11]. To our knowledge, our occurrence dataset (S1 Appendix) is the most comprehensive for these species in Brazil to date, which contains the geocoded information on the occurrence of pupae, larvae or eggs, and adults of each species recorded between 1922–2019 in Brazil.

The sources of these data were from (i) The Global Biodiversity Information Facility (GBIF; https://www.gbif.org/) database, (ii) published literature (see list of papers in occurrence data set), and (iii) primary occurrence data from entomological surveys collected by the Health Secretariat of the State of Rio Grande do Sul and Projeto Cantareira as part of São Paulo State Superintendence for Endemic Diseases Control (SUCEN). To retrieve occurrences from the literature, we conducted searches on Scopus and PubMed using the search strings “haemagogus AND Brazil” and “sabethes AND Brazil” (Fig 2). We used the geographical coordinates (latitude and longitude) of each observed occurrence when available. In cases where only an area location was mentioned (e.g. name of town or city), we used Google Maps to extract the coordinates of the centroid of the cited location. The number of records retrieved by year is presented in the S1 Fig. To permit reproducibility, all occurrence data containing the aforementioned information are made freely accessible online on https://github.com/sabrinalyli/YFVMappingBrazil.git.

thumbnail
Fig 2. Number of occurrence records for Hg. janthinomys, Hg. leucocelaenus, and Sabethes spp. (1991–2019) prior and post standardization.

https://doi.org/10.1371/journal.pntd.0010019.g002

We performed spatial and temporal standardization using methods described elsewhere [14,22]. This creates a list of occurrences in which each occurrence is individually linked with a unique location and year. This ensures that occurrences from regions where multiple surveys are conducted are not overly-represented compared to occurrences from less well sampled regions. To validate the geo-positioning of occurrence locations, we overlaid the occurrence points with a map of Brazil to ensure that points were geocoded to land only. Point locations that were incorrectly geocoded to a location over water were adjusted using an algorithm described elsewhere [23] or manually checked for validity.

Environment and land-use variables

Data on the spatial and temporal distribution of environmental conditions are extracted from high-resolution satellite raster imagery. We standardized all covariates to a 1 km x 1 km spatial resolution by resampling each raster using the bilinear method [24] and have outlined the rationale for inclusion in our models below. A summary of metadata for our covariate datasets along with details of covariate sensitivity analysis can be found in S1 Table and S3 Fig. All data processing was done in R (Vers. 4.0.2) using the rgdal, raster, and sf packages.

Temperature, precipitation, and relative humidity

Climatic factors such as temperature, relative humidity, and precipitation influence the abundance and activity of Haemagogus and Sabethes spp. [2527]. These mosquitoes exhibit daytime biting behaviour [28,29]; in particular, Hg. leucocelaenus only travel closer to the ground to feed during the day [30]. To account for temperature in our model, we adopted mean annual daytime land surface temperature layers extracted from MODIS v6 time series processed using the Temporal Fourier Analysis (TFA) described elsewhere [31]. For precipitation, we calculated yearly mean total precipitation, yearly minimum precipitation, and yearly maximum precipitation, averaged between years 2001–2019 using data retrieved from ERA5 dataset provided by the European Centre for Medium-Range Weather Forecasts (ECMWF). Minimum annual relative humidity from 2015 was determined as a percentage of water vapour needed to saturate the air given surrounding temperature [32], and was extracted for each 1 km x 1 km pixel.

Enhanced Vegetation Index (EVI) and vegetation type

Haemagogous and Sabethes species have an affinity for forest canopies and inhabit tree holes, bamboo nodes, bromelias, and coconut husks in vegetated areas [33,34]. Many complete their larval development in temporarily flooded forest areas [35]. Haemagogus and Sabethes spp. are primarily found in rain forests and deciduous forests [25,36]. To characterize vegetation cover, we used EVI retrieved from MODIS, a measure of vegetation canopy greenness ranging from 0–1, where values near 1 indicate lush and dense vegetation. We also included consensus land-cover layers processed by Tuanmu et al. [37] detailing the prevalence of evergreen and deciduous trees.

Forest cover and deforestation

Sylvatic mosquito abundance could be influenced by increasing deforestation activities as a result of regular landscape transformation to accommodate agricultural activities and cattle grazing in recent years [4,38]. To characterize forest cover and deforestation, we adopted gross forest cover loss for the years 2000–2019, using data processed by Hansen et al. [39]. We also included a consensus land cover on cultivated and managed vegetation to account for land that has been anthropogenically modified for agricultural purposes.

Elevation

While Haemagogus and Sabethes spp. are known to prefer canopy habitats within the forest strata, they can also be found near the ground depending on the time of the year and climate conditions [11,40,41]. Thus, to account for vertical oviposition, we extracted minimum elevation in meters above sea level from a raster processed by Amatulli et al. [42].

Population density

Higher population density in cities has the potential to increase physical contact between people and YFV mosquito vectors, which could potentially instigate and maintain another urban YFV transmission pathway in Brazil [43]. Furthermore, including population density in our model allows us to account for the possible occurrence of Hg. janthinomys and Hg. leucocelaenus in forest-city ecotones. To account for population density, we included population density estimates adjusted to match country estimates reported by the United Nations from WorldPop (https://www.worldpop.org/) calculated using methods described elsewhere [44].

Non-human primate reservoirs

We retrieved the presence probability of NHP reservoirs predicted at 1 km x 1 km resolution for six species within the Southeastern Atlantic forest in Brazil, which was determined using occurrence records from 1992–2016 retrieved from Culot et al. [17]. Methods for producing these rasters are described elsewhere [45]. The six species of monkeys were of the following genera: Alouatta sp., Brachyteles sp., Callicebus sp., Callithrix sp., Leontopithecus sp., and Sapajus sp. These species have been found to host and transmit YFV in Brazil [46,47]. In particular, Howler monkeys (Alouatta sp.) are extremely susceptible to YFV and can develop fatal disease [48].

Modelling approach

Our stacked binomial regression ensemble consists of an elastic-net regularized GLM, extreme gradient boosted regression trees, and random forest. An elastic-net regularized GLM can be used as a binary classifier that adopts regularization to shrink the coefficients of the linear model to reduce overlearning and to allow for better generalization [49]. Both extreme gradient boosted regression trees and random forest use an ensemble of regression trees to increase prediction accuracy [50]. The method of gradient boosting builds trees in a stagewise process, where new trees help to correct errors created by the previous tree, while random forest trains each tree independently using a randomly sampled training dataset that makes the model more robust and less prone to overfitting [50]. Boosted regression trees, first described for disease mapping in Elith et al. [51], have been adopted to map several arboviral diseases and vectors [5254]. The method has also been previously adopted to estimate environmental suitability of Aedes vectors known to facilitate arbovirus transmission [14,55]. Within the ensemble, we adopted extreme gradient boosting, which is designed to be computationally efficient with better model performance [56].

From the ensemble, we retrieved a predicted probability of suitability at 1 km x 1 km spatial resolution. To obtain uncertainty estimates around the prediction, we fitted an ensemble of 100 sub-models by applying a separate bootstrap sampling of the dataset for each model. We then determined the distribution of pixel-based uncertainty by calculating the 95% confidence interval derived from the 100 sub-models. All data analysis were conducted using raster, gbm, dismo, and seegSDM packages in R (Vers. 4.0.2). All models were fitted using the Superlearner, glmnet, xgboost, and ranger packages in R.

Pseudo-absence records for removing sample selection bias

A sample of background data defined by absence records is required by BRT models to understand conditions available to the modelled species in an area. Phillips et al. [21] proposed a pseudo-absence selection approach using background records that share the same environmental characteristics as the presence records. This minimizes sample selection bias while highlighting the differentiation between the distribution of presence records and that of the background. We adopted this method by compiling a pseudo-absence dataset containing records of similar acrodendrophilic mosquitoes that belong to the same Culicidae family and are captured using the same methods and equipment in the field. A total number of 1,056 records including Culex spp., Limatus spp., Wyeomyia spp., Psorophora spp., Orthopodomyia spp., Mansonia spp., and Coquillettidia spp. observed during 1992–2020 were extracted from the GBIF database (S1 Appendix). The locations of these records are shown in S2 Fig.

Mapping suitable areas, risk of transmission, and risk of spillover

To better visualize areas that are suitable for vector species, we mapped predicted suitability above a certain threshold. This was done by highlighting pixels where predicted suitability was greater than 0.5 for each species and for all species. We also conducted a sensitivity analysis by changing the threshold to 0.25 and 0.75. Using these data, we mapped the risk of transmission and risk of spillover from sylvatic environments by determining areas of overlap between two to three variables: high environmental suitability of vector species, high predicted presence probability of NHPs, and presence of human population (yes/no) based on population density. High environmental suitability was defined by suitable areas filtered by the threshold. Likewise, high predicted presence probability of NHPs was defined by a probability greater than 0.5.

We defined the risk of transmission by the presence of two variables, which interact to maintain transmission. Firstly, a transmission cycle can be maintained by the presence of sylvatic mosquito vectors and NHP hosts. Secondly, although YF transmission can be maintained by Aedes spp., it was not observed during the 2016–2018 outbreaks. Nonetheless, increasing overlap between human populations and NHP host reservoirs could instigate an urban cycle. Lastly, the presence of human populations in areas where sylvatic vectors are present can also foster transmission. Similarly, we defined the risk of spillover by the presence of all three variables in an area, where local transmission is maintained by the presence of mosquito vectors and NHP hosts that feed on nearby human populations. Finally, areas where only one variable is present are defined as low risk. Despite the lack of overlap, transmission could still occur if humans, NHP hosts, and mosquito vectors are assumed to be mobile within given suitable environmental conditions and reasonable travel behaviour.

Model validation

We calculated the overall predictive accuracy of the model ensemble by quantifying the area under the curve (AUC) statistic and its 95% confidence interval. This was done by running our model ensemble on a randomly sampled subsection of the data referred to as the training set of the data, which contains 80% of the original data. Cross-validation on each model was performed by evaluating the mean area under the curve (AUC) statistic from 10 separate cross-validation folds. This model was then applied for prediction on the holdout set (20% of the original data). Using 2000 bootstraps, we computed the 95% confidence interval. All analyses were performed using the ROCR [57] and pROC [58] packages in R.

All covariates used in the model were assessed for their variable importance using the permutation method [59,60], which is measured as a change in the model error. Values of each covariate in the data were randomly permuted and assessed for its effect on the model performance using the difference between a baseline performance measure (e.g. RMSE) before and after permutation. This was performed on the stacked binomial regression ensemble using the vip package in R [61].

Results

To characterize the presence of vectors we retrieved a total of 1,091 occurrences (Hg. janthinomys = 176; Hg. leucocelaenus = 539; Sabethes spp. = 376) over the time period 1991–2019 after standardization. Out of all the occurrences recorded for Hg. janthinomys, most were recorded in the states of Rio de Janeiro (RJ: 27.3%), São Paulo (SP: 17.6%), and Bahia (BA:15.9%) (Fig 3A and 3B) likely due to higher surveillance in recent years as a result of recent outbreaks. For Hg. leucocelaenus, 55% of all occurrences were recorded in Rio Grande do Sul (RS), while 20.5% and 17.6% were recorded in São Paulo and Rio de Janeiro, respectively (Fig 3C). The distribution of Sabethes spp. occurrences adopted in our study were found in São Paulo (27.8%), Rio de Janeiro (19.4%), Rio Grande do Sul (14.4%) and Bahia (8.2%), (Fig 3D).

thumbnail
Fig 3. Environmental suitability and occurrence records of Haemagogus and Sabethes spp.

A, Map of Brazil by 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) and regions (colour-coded). B-D, Maps indicating the point locations (red) of standardized occurrence records of each species collected between 1991–2019 and predicted environmental suitability for Hg. janthinomys, Hg. leucocelaenus, and Sabethes spp. The base layer of the map was retrieved from https://www.ibge.gov.br/geociencias/downloads-geociencias.html.

https://doi.org/10.1371/journal.pntd.0010019.g003

Fig 3B and 3D also shows the predicted suitability of Hg. janthinomys, Hg. leucocelaeunus, and Sabethes spp. To evaluate prediction accuracy of the outputs, we determined model performance using the AUC (area under the curve) on our holdout set and found high accuracy (>0.97). While these species are distributed similarly, there are some differences. The Southeastern region part of the Atlantic forest biome encompassing the states of Minas Gerais (MG), Rio de Janeiro (RJ), Paraná (PA), and parts of São Paulo (SP) were predicted to be the most suitable locations for Hg. janthinomys and Hg. leucocelaenus, where the species have been recorded in the last three decades. While few species were recorded in these regions, the Central-West Brazilian states incorporating the Cerrado biome of Goiás (GO) and Mato Grosso do Sul (MS), along with some areas within Tocantins (TO), Federal District (DF), and Mato Grosso (MT) were also found to be suitable. Highly suitable areas were also predicted in the Northeast state of Bahia and Piauí (PI), and nearby adjacent areas in the states of Ceará (CE), Paraíba (PB), and Pernambuco (PE), where the species have yet to be reported. The Southern states of Rio Grande do Sul (RS), Paraná (PR), and Santa Catarina (SC) were predicted to be highly suitable for both Hg. leucocelaenus and Sabethes spp. For all species, suitability was predicted to be low across the Northern region along the Amazon river basin, where contemporary records of species are sparse. Maps summarizing uncertainty associated with predicted estimates are shown in S4 Fig.

These patterns are further reflected in the maps showing suitable areas (suitability probability ≥ 0.5) for each species (Fig 4A and 4C) and for all species combined (Fig 4D). Our sensitivity analysis shows that 0.5 is an appropriate- threshold to avoid underestimation and overestimation of suitable zones S5 Fig. States that are suitable for all species are mainly concentrated in highly populated areas of the Southeastern region and along part of Northeast and Central-West.

thumbnail
Fig 4. Suitable regions for Haemagogus and Sabethes spp.

A-C, Maps of areas identified as suitable with environmental suitability ≥ 0.5. D, Map showing areas suitable for only one species, two species, or all species analyzed in this study. The base layer of the map was retrieved from https://www.ibge.gov.br/geociencias/downloads-geociencias.html.

https://doi.org/10.1371/journal.pntd.0010019.g004

Likewise, our risk map based on the presence of mosquitoes, NHP host reservoirs, and human populations (Fig 5) shows that areas at risk of YFV spillover are concentrated along highly populated municipalities of Paraná, the cities of São José do Rio Preto and São Paulo, as well as municipalities in Rio de Janeiro and Espírito Santo states. Areas at risk of spillover also include the Federal District and Goiás. These areas are surrounded by comparatively less populated municipalities which are at the risk of transmission. The regions between Goiás and Minas Gerais have the lowest risk given the lack of NHP host reservoirs. We superimposed our spillover risk maps with data on confirmed (PCR, immunohistochemistry and/or IgM ELISAs) YFV cases among humans and non-human primates between 2015–2020 [62] retrieved from the Brazilian Ministry of Health National Notifiable Disease System (SINAN; Sistema de Informação de Agravos de Notificação). This shows that previous human and NHP cases clustered mainly in spillover zones. However, the original site of human infection may be different than the site of reporting attributed to passive surveillance coupled with human movement.

thumbnail
Fig 5. Risk of transmission and spillover along the Southeastern Atlantic Forest Biome.

Map showing the distribution in risk status in the Southeastern municipalities of São Paulo, Minas Gerais, Rio de Janeiro, Espírito Santo, and Southern state of Paraná and adjacent municipalities in the Central-West, including the Federal District, Goiás, and Mato Grosso. Risk is defined by the presence of sylvatic mosquito vectors, non-human primate reservoirs, and human population density. Low risk refers to the presence of only one of these factors. Risk of transmission is defined by the overlapping presence of two of these factors, while risk of spillover is defined by the overlapping presence of all three factors. Blue points show the locations of confirmed cases among non-human primates (NHP) recorded between years 2015–2020 while green triangles indicate the locations of confirmed human cases recorded between 2016–2019. The base layer of the map was retrieved from https://www.ibge.gov.br/geociencias/downloads-geociencias.html.

https://doi.org/10.1371/journal.pntd.0010019.g005

Predictors with the strongest variable importance on the model ensemble for Hg. janthinomys included population density (permuted model error: 0.106±0.001), which shares a negative relationship with species distribution, the prevalence of deciduous tree cover (0.1±0.004), and elevation (0.075±0.001). For modelling Hg. leucocelaenus, the most influential predictors were elevation (0.097±0.002) and minimum daytime land surface temperature (0.069±0.003), followed by population density (0.067±0.002). Similar results were found for modelling Sabethes spp., where the prevalence of deciduous tree cover (0.091±0.002), population density (0.084±0.002), and elevation (0.078±0.002) had the strongest relative influence. Relative to all covariates, mean forest cover loss had consistently one of the least variable importance for all models. Variable importance along with partial dependence plots for Hg. janthinomys, Hg. leucocelaenus, and Sabethes spp can be viewed interactively from https://github.com/sabrinalyli/YFVMappingBrazil. Overall predictive performance was high for all model ensembles based on computed AUC (S6 Fig).

Discussion

To understand the spatial distribution and patterns of YFV vectors and their risk to humans, we used a stacked regression ensemble modelling approach to map the environmental suitability of main YFV vectors from the Haemagogus and Sabethes genera at a high spatial resolution of 1 km x 1 km. This was done using occurrence records reported between 1991–2019 combined with associated environmental and land-cover covariates. Using predicted estimates, we identified suitable regions for each vector and all vectors. We also mapped the distribution of risk of transmission and spillover in regions that overlapped human settlements and the presence probability of NHP host reservoirs.

Previous studies by Childs et al. [15] and Almeida et al. [16] used MaxEnt models with different assumptions about the selection of pseudo-absences. Using a stacked binomial regression model ensemble and pseudo-absence selection based on background records that share the same environmental characteristics as the presence records, we found higher suitability in the Northeast and lower suitability in the North, which is the opposite of the findings presented by Childs et al. [15]. Similar to Almeida et al. [16], our model ensemble also predicted high suitability for Hg. leucocelaenus in Rio Grande do Sul. Compared to Almeida et al. [16], the integration of additional occurrence points in our model ensemble allowed the prediction of higher suitability for the south of Rio Grande do Sul. Whereas temperature did not play an important role in model prediction in Childs et al. [15], we found it an important covariate in ranking after elevation, population density, and deciduous forest cover.

Despite the existence of effective vaccines, Brazil experienced the largest sporadic YF outbreaks of recent years during 2016–2018, which resulted in more than 2,153 diagnosed human cases and 744 deaths in highly populated urban areas of Southeast Brazil that were previously considered low risk [63]. These outbreaks, which were identified to be of sylvatic origin [1], spilled into the transition zone between the Cerrado biome and the Atlantic Forest in Minas Gerais and moved towards the Southeastern states of Minas Gerais, Espírito Santo, Rio de Janeiro, and São Paulo, (61), which corresponds to predicted areas of high environmental suitability for primary and secondary sylvatic vectors in our study. We also identified highly suitable areas in the South and Central-West regions, including Santa Catarina and Rio Grande do Sul, where outbreaks were predicted to occur every six to ten years [4] and have been recently reported for 2019 and 2021 [64,65]. These findings are consistent with recent spatial epidemiological patterns of sylvatic YF outbreaks in the Central-West and Southeastern regions [7,66]. Low predicted suitability in the Amazon river basin is likely due to decreased entomological reporting in recent years, mainly due to lack of surveillance of sylvatic vectors in the Amazon region. These maps can be used as a public health tool and help guide the surveillance of YFV vectors, especially in areas lacking entomological data but were identified as environmentally suitable.

Our understanding of each species’ ecology is confirmed by the relative influence of environmental covariates on the models. The distribution of each species may be dependent on the height of treetop canopies, but further investigation is required. Here we used elevation to characterize vertical oviposition activity. A previous study found that the risk of YFV cases at elevations between 300–800 MASL was six times higher compared to other elevations [67], indicating a potential linkage to sylvatic mosquito activity. Population density strongly influenced species’ distribution via an inverse relationship, corroborating species’ preference for non-urbanized settings. The prevalence of deciduous tree cover had a strong influence on shaping the distribution of both Hg. janthinomys and Hg. leucocelaenus. Minimum day-time temperature had a strong influence on both Hg. leucocelaenus and Sabethes spp. species. In particular, higher suitability was detected when temperature was above 27°C, which is consistent with recent findings on temperature-population dynamics of Hg. leucocelaenus [25,68].

Forest loss was found to have an overall weak influence on the models for all three species. Evidence on the relationship between deforestation and mosquito abundance is scarce and remains unclear. In Brazil, Hg. janthinomys and Hg. leucocelaenus are highly dispersed between forest areas, and can be found in small forest patches (e.g. 7 hectares) next to urban areas and isolated high-income neighbourhoods located in forest valleys [1], In particular, Hg. leucocelaenus are often found in modified forest patches near human activities, including isolated forest fragments in cultivated areas and in pasture grasslands generally associated with small rivers and creeks as well as urban green areas [30,69,70]. A systematic review that investigated the effects of deforestation on mosquitoes found that half of mosquito species favoured deforested areas [71]. However the loss of mosquito habitats from the progression of deforestation may also reduce transmission risk [72]. Moreover, some considerations have also been given to the changing agricultural ecosystem in Brazil, which is a main driver of deforestation. Seasonally varying agricultural activities such as planting and harvesting to accommodate the growth of cash crops have been found to capture seasonal YFV reporting in humans and NHPs [38], but little is known about the impacts of regular land conversion through agriculture on mosquito habitats and transmission activity. Nonetheless, more field observations and investigations are needed to understand the effects of agriculture and deforestation on mosquito abundance.

By combining predicted estimates of vector suitability, predicted presence of NHP host reservoirs, and population density, we further characterized and mapped the distribution in risk of YFV transmission and spillover within and surrounding the Southeastern Atlantic Forest biome. Our detection of spillover and transmission risk in and around Minas Gerais, Rio de Janeiro, Espírito Santo, and São Paulo, where previous spillover events have been recorded (62), suggests that disease control and prevention should continue to focus in those areas. Furthermore, our map highlighting spillover risk in Paraná also coincides with the occurrence of epizootics in the state between 2019 and early 2021 [73]. Given that YFV cases among NHPs are recorded through passive surveillance and are thus often based in areas with human settlements, we cannot rule out the fact that varying levels of risk may subsist in non-spillover zones that are not identified on our map. Despite this, these results do not reflect the risk of YFV outbreaks as population level vaccination coverage is not considered. Even if vaccination coverage is considered, the heterogeneity of current vaccination data by municipality may not accurately portray the coverage gradient across space at the local level. Furthermore, we also did not consider the dynamics of human mobility, which also contribute to outbreak risk. Nonetheless, our maps can help guide public health officials in detecting at-risk municipalities, where further investigation can be conducted by comparing vaccination coverage and effectiveness.

In September 2018, only a handful of at-risk municipalities within the states in our study region reached 95% vaccination coverage; this included 13% and 19% of municipalities in Paraná (13%) São Paulo, respectively [74]. However as of May 2020, coverage in the majority of the municipalities in São Paulo, Paraná, Santa Catarina, Espírito Santo, and Minas Gerais was still less than 60% for individuals aged 15 years and older, and less than 80% for individuals aged 14 years and younger [75]. Despite improved coverage, the lack of seroprevalence data for YF presents uncertainty in the disease burden experienced by at-risk populations. A recent seroepidemiological study focusing on the Serro region in Minas Gerais found that about 26% of urban and rural populations did not present neutralizing antibodies against YFV, including some that were vaccinated [76]. In particular, children were less likely to present YFV-neutralizing antibodies and were less likely to have received a booster dose. Furthermore, there has been a reduction in childhood immunization in Brazil from 2017 to 2019 [77]. This highlights potential underestimation of disease burden and the likelihood of new outbreaks originating among vulnerable populations across all ages living in at-risk areas. Moreover, it is uncertain whether current vaccination efforts have been hindered by the co-existence of the COVID-19 pandemic, which may have redirected outbreak surveillance efforts and the allocation of health resources [78].

Our study has a few limitations. It does not account for seasonal patterns of species distribution, which could influence the distribution in environmental suitability of each species throughout the year. Given that the abundance of Haemagogus and Sabethes spp. is highest between December and May during the rainy season, the risk of transmission and spillover may also vary throughout the year [79]. Despite assembling the most comprehensive database to date, our data comes from a few specific geographic areas and thus may not cover the full range of ecological conditions that influence each species. This could be addressed by limiting predictions to locations outside of model training regions [80] or using a temperature suitability index [81]. Due to limited lab and field evidence on how temperature affects the extrinsic incubation period of these species, their gonotrophic cycle, and oviposition behaviour, future entomological surveys could focus on these observations in various locations to update our knowledge which is mainly based on lab evidence presented from a few decades ago [8284]. Similarly, while our pseudo-absence records on Culicidae species are assembled from GBIF only, future works could explore the addition of records from literature. While doing so may not overcome the issue of differential surveillance efforts conducted by various research teams, it could potentially minimise selection bias and improve modelled outputs.

Our finding of restricted spatial distribution in environmental suitability for Sabethes spp. concur with previous field observations. Sabethes spp. was observed in low abundance, distribution, and infection rates and played a local or secondary role during the 2016–2018 outbreaks [1]. Previous studies have noted the competence of Sa. albiprivus, Sa. petrocchiae, and Sa. chloropterus in transmitting YFV [85], however some species are limited to certain ecological conditions. Among the Sabethes genera, Sa. chloropterus was found mainly in forest areas and while considered the primary YFV sylvatic vector [2], was found with low infection rates during previous outbreaks [1,86].

While beyond the scope of this paper, the role of Aedes spp. mosquitoes and its interactions with the sylvatic cycle in facilitating YFV transmission should be further investigated given conflicting findings. In Brazil, evidence of YFV transmission from Ae. aegypti and Ae. albopictus, which are sensitive to YFV infection, has not been observed since 1942 [1] but other species such as Ae. scapularis and Ae. taeniorhynchus were found to be naturally infected during the 2016–2018 outbreaks. While these occasional findings do not define their role as YFV vectors, more investigation is needed on their potential vectorial capacity. In addition, considering human mobility in future studies could also help to identify changes in disease spread facilitated by viremic humans, including those that are asymptomatic, which can be applied to infer the geographical shift in risk over time. Furthermore, current uncertainties related to our knowledge of NHP host species, their viremia, and the spatial distributions of their reservoirs [87] should be further studied. While most NHP exist in sylvatic environments, some are capable of surviving in degraded environments [88]. NHP behaviour of various species across Brazil should be closely monitored to identify reservoir establishment and for epizootics to delimit transmission. In this aspect, the promising application of phylogenetic data and genomic analysis could help to identify the roles of NHP species in facilitating transmission, the variation of YFV infection rates across NHP species, and trends in epizootic cases [89].

In conclusion, our findings expand on our limited understanding on the geographic distribution of potential sylvatic vector reservoirs in Brazil. This work contributes to the growing body of research on vector and vector-borne disease mapping at high spatial resolution in low-and middle-income countries that lack comprehensive surveillance coverage. Our identification of suitable regions for main YFV vectors nationally and areas at risk of YFV spillover in the Southeastern region of Brazil could help guide knowledge users and policy makers working in vector surveillance and vaccine allocation.

Supporting information

S1 Fig.

Temporal distribution of occurrence records for Hg. janthinomys, Hg. leucocelaenus, and Sabethes spp.

https://doi.org/10.1371/journal.pntd.0010019.s001

(TIFF)

S2 Fig.

Distribution of pseudo-absence points. The base layer of the map was retrieved from https://www.ibge.gov.br/geociencias/downloads-geociencias.html.

https://doi.org/10.1371/journal.pntd.0010019.s002

(TIFF)

S3 Fig.

Maps illustrating the distribution of environmental covariate data. The base layer of the map was retrieved from https://www.ibge.gov.br/geociencias/downloads-geociencias.html.

https://doi.org/10.1371/journal.pntd.0010019.s003

(TIFF)

S4 Fig.

(A) Maps of prediction uncertainty calculated based on the 95% confidence interval and (B) frequency of presence occurrence points falling in each predicted range of suitability. The base layer of the map was retrieved from https://www.ibge.gov.br/geociencias/downloads-geociencias.html.

https://doi.org/10.1371/journal.pntd.0010019.s004

(TIFF)

S5 Fig.

Maps showing suitable regions for Haemagogus and Sabethes spp. when environmental suitability threshold is set to (A) ≥ 0.25 and (B) ≥ 0.75. The base layer of the map was retrieved from https://www.ibge.gov.br/geociencias/downloads-geociencias.html.

https://doi.org/10.1371/journal.pntd.0010019.s005

(TIFF)

S6 Fig.

AUC and 95% confidence interval generated for the holdout set for Hg. janthinomys, Hg. leucocelaenus, and Sabethes spp.

https://doi.org/10.1371/journal.pntd.0010019.s006

(TIFF)

S1 Table.

Metadata of environmental and land-cover covariates used in stacked binomial regression model.

https://doi.org/10.1371/journal.pntd.0010019.s007

(DOCX)

S1 Appendix.

List of occurrence data downloads (with DOI) from Global Diversity Information Facility (GBIF).

https://doi.org/10.1371/journal.pntd.0010019.s008

(DOCX)

References

  1. 1. de Abreu FVS, Ribeiro IP, Ferreira-de-Brito A, dos Santos AAC, de Miranda RM, Bonelly I de S, et al. Haemagogus leucocelaenus and Haemagogus janthinomys are the primary vectors in the major yellow fever outbreak in Brazil, 2016–2018. Emerg Microbes Infect. 2019;8: 218–231. pmid:30866775
  2. 2. Silva NIO, Sacchetto L, de Rezende IM, Trindade G de S, LaBeaud AD, de Thoisy B, et al. Recent sylvatic yellow fever virus transmission in Brazil: the news from an old disease. Virol J. 2020;17: 9. pmid:31973727
  3. 3. Johansson MA, Vasconcelos PFC, Staples JE. The whole iceberg: estimating the incidence of yellow fever virus infection from the number of severe cases. Trans R Soc Trop Med Hyg. 2014;108: 482–487. pmid:24980556
  4. 4. Monath TP, Vasconcelos PFC. Yellow fever. J Clin Virol. 2015;64: 160–173. pmid:25453327
  5. 5. Benchimol JL, editor. Febre amarela: a doença e a vacina, uma história inacabada [online]. Rio de Janeiro: Editora FIOCRUZ, 2001. 470 p. online. Rio de Janeiro: Fiocruz; 2001.
  6. 6. Taylor RM. Epidemiology. In: Strode George K., Bugher John C., editors. Yellow Fever. New York: McGraw-Hill; 1951. pp. 427–459.
  7. 7. de Oliveira Figueiredo P, Stoffella-Dutra AG, Barbosa Costa G, Silva de Oliveira J, Dourado Amaral C, Duarte Santos J, et al. Re-Emergence of Yellow Fever in Brazil during 2016–2019: Challenges, Lessons Learned, and Perspectives. Viruses. 2020;12. pmid:33143114
  8. 8. Pinheiro GG, Rocha MN, de Oliveira MA, Moreira LA, Andrade Filho JD. Detection of Yellow Fever Virus in Sylvatic Mosquitoes during Disease Outbreaks of 2017–2018 in Minas Gerais State, Brazil. Insects. 2019;10: 136. pmid:31083286
  9. 9. Lourenço-de-Oliveira R, Failloux A-B. High risk for chikungunya virus to initiate an enzootic sylvatic cycle in the tropical Americas. PLoS Negl Trop Dis. 2017;11: e0005698. pmid:28662031
  10. 10. Alencar J, Morone F, De Mello CF, Dégallier N, Lucio PS, Da Serra-Freire NM, et al. Flight Height Preference for Oviposition of Mosquito (Diptera: Culicidae) Vectors of Sylvatic Yellow Fever Virus Near the Hydroelectric Reservoir of Simplício, Minas Gerais, Brazil. J Med Entomol. 2013;50: 791–795. pmid:23926776
  11. 11. Pinto CS, Confalonieri UEC, Mascarenhas BM. Ecology of Haemagogus sp. and Sabethes sp. (Diptera: Culicidae) in relation to the microclimates of the Caxiuanã National Forest, Pará, Brazil. Mem Inst Oswaldo Cruz. 2009;104: 592–598. pmid:19722082
  12. 12. Warren AJ. Landmarks in the conquest of yellow fever. In: Strode GK, Bugher JC, editors. Yellow fever. New York: McGraw-Hill; 1951. pp. 5–37.
  13. 13. Elith J, Leathwick JR. Species Distribution Models: Ecological Explanation and Prediction Across Space and Time. Annu Rev Ecol Evol Syst. 2009;40: 677–697.
  14. 14. Kraemer MUG, Sinka ME, Duda KA, Mylne AQN, Shearer FM, Barker CM, et al. The global distribution of the arbovirus vectors Aedes aegypti and Ae. albopictus. eLife. 2015;4: e08347. pmid:26126267
  15. 15. Childs ML, Nova N, Colvin J, Mordecai EA. Mosquito and primate ecology predict human risk of yellow fever virus spillover in Brazil. Philos Trans R Soc B Biol Sci. 2019;374: 20180335. pmid:31401964
  16. 16. Almeida MAB, Santos E dos, Cardoso J da C, Silva LG da, Rabelo RM, Bicca-Marques JC. Predicting Yellow Fever Through Species Distribution Modeling of Virus, Vector, and Monkeys. EcoHealth. 2019;16: 95–108. pmid:30560394
  17. 17. Culot L, Pereira LA, Agostini I, Almeida MAB de, Alves RSC, Aximoff I, et al. ATLANTIC-PRIMATES: a dataset of communities and occurrences of primates in the Atlantic Forests of South America. Ecology. 2019;100: e02525. pmid:30317556
  18. 18. Laan MJ van der, Polley EC, Hubbard AE. Super Learner. Stat Appl Genet Mol Biol. 2007;6. pmid:17910531
  19. 19. Bhatt S, Cameron E, Flaxman SR, Weiss DJ, Smith DL, Gething PW. Improved prediction accuracy for disease risk mapping using Gaussian process stacked generalization. J R Soc Interface. 2017;14: 20170520. pmid:28931634
  20. 20. Hastie T, Tibshirani R, Friedman JH. The Elements of Statistical Learning: Data Mining, Inference, and Prediction. Springer; 2009.
  21. 21. Phillips SJ, Dudík M, Elith J, Graham CH, Lehmann A, Leathwick J, et al. Sample selection bias and presence-only distribution models: implications for background and pseudo-absence data. Ecol Appl. 2009;19: 181–197. pmid:19323182
  22. 22. Messina JP, Brady OJ, Pigott DM, Brownstein JS, Hoen AG, Hay SI. A global compendium of human dengue virus occurrence. Sci Data. 2014;1: 140004. pmid:25977762
  23. 23. Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013;496: 504–507. pmid:23563266
  24. 24. Hijmans RJ, van Etten J. raster: Geographic data analysis and modeling. R Package Version. 2016;2.
  25. 25. Alencar J, Mello CF de, Morone F, Albuquerque HG, Serra-Freire NM, Gleiser RM, et al. 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. 2018;34: 85–92. pmid:31442164
  26. 26. de Araújo WS, Vieira TM, de Souza GA, Bezerra IC, Corgosinho PHC, Borges MAZ. Nocturnal Mosquitoes of Pará State in the Brazilian Amazon: Species Composition, Habitat Segregation, and Seasonal Variation. J Med Entomol. 2020 [cited 2 Oct 2020]. pmid:32484514
  27. 27. La Corte R, Maia PCR, Dolabella SS, Cruz DER, Marteis LS. Mosquitoes of the Caatinga. III. Larval Habitats, Frequency, and Dynamics of Immature and Adult Stages in a Dry Brazilian Forest. J Med Entomol. 2019;56: 120–128. pmid:30247710
  28. 28. Marteis LS, Natal D, Sallum MAM, Medeiros-Sousa AR, Oliveira TMP, La Corte R. Mosquitoes of the Caatinga: 1. Adults stage survey and the emerge of seven news species endemic of a dry tropical forest in Brazil. Acta Trop. 2017;166: 193–201. pmid:27876648
  29. 29. Santos EB, Favretto MA, Müller GA. When and what time? On the seasonal and daily patterns of mosquitoes (Diptera: Culicidae) in an Atlantic Forest remnant from Southern Brazil. Austral Entomol. 2020;59: 337–344.
  30. 30. Cardoso J da C, de Almeida MAB, dos Santos E, da Fonseca DF, Sallum MAM, Noll CA, et al. Yellow Fever Virus in Haemagogus leucocelaenus and Aedes serratus Mosquitoes, Southern Brazil, 2008. Emerg Infect Dis. 2010;16: 1918–1924. pmid:21122222
  31. 31. Scharlemann JPW, Benz D, Hay SI, Purse BV, Tatem AJ, Wint GRW, et al. Global Data for Ecology and Epidemiology: A Novel Algorithm for Temporal Fourier Processing MODIS Data. PLoS ONE. 2008;3. pmid:18183289
  32. 32. Linacre ET. A simple formula for estimating evaporation rates in various climates, using temperature data alone. Agric Meteorol. 1977;18: 409–424.
  33. 33. Forattini OP. Entomologia médica: 4o volume: psychodidae. phlebotominae. leishmanioses. bartonelose. Entomol Médica 4o Vol Psychodidae Phlebotominae Leishmanioses Bartonelose. 1973; 658–658.
  34. 34. Silva-Inacio CL, Paiva AAP de, Araújo JMG de, Ximenes M de FF de M, Silva-Inacio CL, Paiva AAP de, et al. Ecological relationships of Haemagogus spegazzinii (Diptera: Culicidae) in a semiarid area of Brazil. Rev Soc Bras Med Trop. 2020;53. pmid:33263687
  35. 35. Tátila-Ferreira A, Maia D de A, de Abreu FVS, Rodrigues WC, Alencar J. Oviposition behavior of Haemagogus leucocelaenus (Diptera: culicidae), a vector of wild yellow fever in Brazil. Rev Inst Med Trop São Paulo. 2017;59. pmid:28793027
  36. 36. Mucci LF, Júnior RPC, de Paula MB, Scandar SAS, Pacchioni ML, Fernandes A, et al. Feeding habits of mosquitoes (Diptera: Culicidae) in an area of sylvatic transmission of yellow fever in the state of São Paulo, Brazil. J Venom Anim Toxins Trop Dis. 2015;21: 6. pmid:25810711
  37. 37. Tuanmu M-N, Jetz W. A global 1-km consensus land-cover product for biodiversity and ecosystem modelling. Glob Ecol Biogeogr. 2014;23: 1031–1045. https://doi.org/10.1111/geb.12182
  38. 38. Hamlet A, Ramos DG, Gaythorpe KAM, Romano APM, Garske T, Ferguson NM. Seasonality of agricultural exposure as an important predictor of seasonal yellow fever spillover in Brazil. Nat Commun. 2021;12: 3647. pmid:34131128
  39. 39. Hansen MC, Potapov PV, Moore R, Hancher M, Turubanova SA, Tyukavina A, et al. High-Resolution Global Maps of 21st-Century Forest Cover Change. Science. 2013;342: 850–853. pmid:24233722
  40. 40. Alencar J, Mello CF de, Gil-Santana HR, Guimarães AÉ, Almeida SAS de, Gleiser RM. Vertical oviposition activity of mosquitoes in the Atlantic Forest of Brazil with emphasis on the sylvan vector, Haemagogus leucocelaenus (Diptera: Culicidae). J Vector Ecol. 2016;41: 18–26. pmid:27232120
  41. 41. Hendy A, Hernandez-Acosta E, Valério D, Mendonça C, Costa ER, Júnior JTA, et al. The vertical stratification of potential bridge vectors of mosquito-borne viruses in a central Amazonian forest bordering Manaus, Brazil. Sci Rep. 2020;10: 18254. pmid:33106507
  42. 42. Amatulli G, Domisch S, Tuanmu M-N, Parmentier B, Ranipeta A, Malczyk J, et al. A suite of global, cross-scale topographic variables for environmental and biodiversity modeling. Sci Data. 2018;5: 180040. pmid:29557978
  43. 43. Possas C, Lourenço-de-Oliveira R, Tauil PL, Pinheiro F de P, Pissinatti A, da Cunha RV, et al. Yellow fever outbreak in Brazil: the puzzle of rapid viral spread and challenges for immunisation. Mem Inst Oswaldo Cruz. 2018;113. pmid:30427974
  44. 44. Lloyd CT, Chamberlain H, Kerr D, Yetman G, Pistolesi L, Stevens FR, et al. Global spatio-temporally harmonised datasets for producing high-resolution gridded population distribution datasets. Big Earth Data. 2019;3: 108–139. pmid:31565697
  45. 45. Hamlet ATP. Yellow fever in South America: The role of environment and host on transmission dynamics. 2020 [cited 2 Dec 2021].
  46. 46. Dietz JM, Hankerson SJ, Alexandre BR, Henry MD, Martins AF, Ferraz LP, et al. Yellow fever in Brazil threatens successful recovery of endangered golden lion tamarins. Sci Rep. 2019;9: 12926. pmid:31506447
  47. 47. Sacchetto L, Silva NIO, Rezende IM de, Arruda MS, Costa TA, de Mello ÉM, et al. Neighbor danger: Yellow fever virus epizootics in urban and urban-rural transition areas of Minas Gerais state, during 2017–2018 yellow fever outbreaks in Brazil. PLoS Negl Trop Dis. 2020;14: e0008658. pmid:33017419
  48. 48. Almeida MAB, Cardoso J da C, dos Santos E, da Fonseca DF, Cruz LL, Faraco FJC, et al. Surveillance for Yellow Fever Virus in Non-Human Primates in Southern Brazil, 2001–2011: A Tool for Prioritizing Human Populations for Vaccination. Gubler DJ, editor. PLoS Negl Trop Dis. 2014;8: e2741. pmid:24625681
  49. 49. Zou H, Hastie T. Regularization and variable selection via the elastic net. J R Stat Soc Ser B Stat Methodol. 2005;67: 301–320.
  50. 50. Breiman L, Friedman JH, Olshen RA, Stone CJ. Classification And Regression Trees. Boca Raton: Routledge; 1984. https://doi.org/10.1201/9781315139470
  51. 51. Elith J, Leathwick JR, Hastie T. A working guide to boosted regression trees. J Anim Ecol. 2008;77: 802–813. pmid:18397250
  52. 52. Longbottom J, Browne AJ, Pigott DM, Sinka ME, Golding N, Hay SI, et al. Mapping the spatial distribution of the Japanese encephalitis vector, Culex tritaeniorhynchus Giles, 1901 (Diptera: Culicidae) within areas of Japanese encephalitis risk. Parasit Vectors. 2017;10: 148. pmid:28302156
  53. 53. Messina JP, Kraemer MU, Brady OJ, Pigott DM, Shearer FM, Weiss DJ, et al. Mapping global environmental suitability for Zika virus. eLife. 2016;5. pmid:27090089
  54. 54. Messina JP, Brady OJ, Golding N, Kraemer MUG, Wint GRW, Ray SE, et al. The current and future global distribution and population at risk of dengue. Nat Microbiol. 2019;4: 1508–1515. pmid:31182801
  55. 55. Kraemer MUG, Reiner RC, Brady OJ, Messina JP, Gilbert M, Pigott DM, et al. Past and future spread of the arbovirus vectors Aedes aegypti and Aedes albopictus. Nat Microbiol. 2019;4: 854–863. pmid:30833735
  56. 56. Chen T, Guestrin C. XGBoost: A Scalable Tree Boosting System. Proc 22nd ACM SIGKDD Int Conf Knowl Discov Data Min. 2016; 785–794.
  57. 57. Sing T, Sander O, Beerenwinkel N, Lengauer T. ROCR: visualizing classifier performance in R. Bioinforma Oxf Engl. 2005;21: 3940–3941. pmid:16096348
  58. 58. Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez J-C, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011;12: 77. pmid:21414208
  59. 59. Breiman L. Random Forests. Mach Learn. 2001;45: 5–32.
  60. 60. Fisher A, Rudin C, Dominici F. All Models are Wrong, but Many are Useful: Learning a Variable’s Importance by Studying an Entire Class of Prediction Models Simultaneously. ArXiv180101489 Stat. 2019 [cited 16 Jul 2021]. Available: http://arxiv.org/abs/1801.01489 pmid:34335110
  61. 61. Greenwell BM, Boehmke BC. Variable Importance Plots—An Introduction to the vip Package. R J. 2020;12: 343–366.
  62. 62. Hill SC, Souza R de, Thézé J, Claro I, Aguiar RS, Abade L, et al. Genomic Surveillance of Yellow Fever Virus Epizootic in São Paulo, Brazil, 2016–2018. PLOS Pathog. 2020;16: e1008699. pmid:32764827
  63. 63. Giovanetti M, Mendonça MCL de, Fonseca V, Mares-Guia MA, Fabri A, Xavier J, et al. Yellow Fever Virus Reemergence and Spread in Southeast Brazil, 2016–2019. J Virol. 2019;94. pmid:31597773
  64. 64. Andrade MS, Campos FS, Campos A a. S, Abreu FVS, Melo FL, Cardoso JC, et al. Real-time genomic surveillance during the 2021 re-emergence of the yellow fever virus in Rio Grande do Sul State, Brazil. 2021 Jul p. 2021.07.16.452727.
  65. 65. Ministério da Saúde. Saúde Brasil 2019: 14—Reemergência e manutenção extra-Amazônica da febre amarela no Brasil, 2014 a 2019: principais desafios para a vigilância, a prevenção e o controle. Brasília–DF; 2019 pp. 305–327.
  66. 66. de Rezende IM, Sacchetto L, Munhoz de Mello É, Alves PA, Iani FC de M, Adelino TÉR, et al. Persistence of Yellow fever virus outside the Amazon Basin, causing epidemics in Southeast Brazil, from 2016 to 2018. PLoS Negl Trop Dis. 2018;12. pmid:29864115
  67. 67. Hamrick PN, Aldighieri S, Machado G, Leonel DG, Vilca LM, Uriona S, et al. Geographic patterns and environmental factors associated with human yellow fever presence in the Americas. Barker CM, editor. PLoS Negl Trop Dis. 2017;11: e0005897. pmid:28886023
  68. 68. Couto-Lima D, Andreazzi CS, Leite PJ, Bersot MIL, Alencar J, Lourenço-de-Oliveira R. Seasonal population dynamics of the primary yellow fever vector Haemagogus leucocelaenus (Dyar & Shannon) (Diptera: Culicidae) is mainly influenced by temperature in the Atlantic Forest, southeast Brazil. Mem Inst Oswaldo Cruz. 2020;115. pmid:32696917
  69. 69. Camargo-Neves VLF de, Poletto DW, Rodas LAC, Pachioli ML, Cardoso RP, Scandar SAS, et al. Entomological investigation of a sylvatic yellow fever area in São Paulo State, Brazil. Cad Saúde Pública. 2005;21: 1278–1286. pmid:16021266
  70. 70. Montagner FRG, Silva OS, Jahnke SM. Mosquito species occurrence in association with landscape composition in green urban areas. Braz J Biol. 2018;78: 233–239. pmid:28793030
  71. 71. Burkett-Cadena ND, Vittor AY. Deforestation and vector-borne disease: Forest conversion favors important mosquito vectors of human pathogens. Basic Appl Ecol. 2018;26: 101–110. pmid:34290566
  72. 72. Barros FSM, Honório NA. Deforestation and Malaria on the Amazon Frontier: Larval Clustering of Anopheles darlingi (Diptera: Culicidae) Determines Focal Distribution of Malaria. Am J Trop Med Hyg. 2015;93: 939–953. pmid:26416110
  73. 73. Coordenação-Geral do Programa Nacional de Imunizações do Departamento de Imunização e Doenças Transmissíveis (CGPNI/DEIDT/SVS), Coordenação-Geral de Laboratórios de Saúde Pública do Departamento de Articulação Estratégica de Vigilância em Saúde (CGLAB/DAEVS/SVS). Boletim Epidemiológico: Informe quinzenal sarampo–Brasil, semanas epidemiológicas 43 de 2020 a 1 de 2021. Brasília—DF: Ministério da Saúde: Secretaria de Vigilância em Saúde; 2021 Feb. Available: https://www.gov.br/saude/pt-br/media/pdf/2021/fevereiro/11/boletim_epidemiologico_svs_4.pdf
  74. 74. WHO. Yellow fever–Brazil. In: WHO [Internet]. World Health Organization; 11 Feb 2019 [cited 15 Mar 2021]. Available: http://www.who.int/csr/don/11-february-2019-yellow-fever-brazil/en/
  75. 75. Coordenação-Geral de Vigilância de Agravos e Doenças Não Transmissíveis (CGDANT/DASNT/SVS). Boletim Epidemiológico: Monitoramento dos casos de arboviroses urbanas transmitidas pelo Aedes Aegypti (dengue, chikungunya e zika), Semanas Epidemiológicas 1 a 19, 2020. Brasília/DF: Ministério da Saúde: Secretaria de Vigilância em Saúde; 2020 May. Available: https://antigo.saude.gov.br/images/pdf/2020/May/20/Boletim-epidemiologico-SVS-20-aa.pdf
  76. 76. Stoffella-Dutra AG, Silva de Oliveira J, Barbosa Costa G, Geessien Kroon E, Santos Abrahão J, Desiree LaBeaud A, et al. Absence of YF-neutralizing antibodies in vulnerable populations of Brazil: A warning for epidemiological surveillance and the potential risks for future outbreaks. Vaccine. 2020;38: 6592–6599. pmid:32788140
  77. 77. Césare N, Mota TF, Lopes FFL, Lima ACM, Luzardo R, Quintanilha LF, et al. Longitudinal profiling of the vaccination coverage in Brazil reveals a recent change in the patterns hallmarked by differential reduction across regions. Int J Infect Dis. 2020;98: 275–280. pmid:32619762
  78. 78. Ribeiro VST, Telles JP, Tuon FF. Arboviral diseases and COVID-19 in Brazil: Concerns regarding climatic, sanitation, and endemic scenario. J Med Virol. 2020;92: 2390–2391. pmid:32462677
  79. 79. Cavalcante KRLJ Tauil PL. Epidemiological characteristics of yellow fever in Brazil, 2000–2012. Epidemiol E Serv Saude Rev Sist Unico Saude Bras. 2016;25: 11–20. pmid:27861674
  80. 80. Shearer FM, Longbottom J, Browne AJ, Pigott DM, Brady OJ, Kraemer MUG, et al. Existing and potential infection risk zones of yellow fever worldwide: a modelling analysis. Lancet Glob Health. 2018;6: e270–e278. pmid:29398634
  81. 81. Brady OJ, Golding N, Pigott DM, Kraemer MUG, Messina JP, Reiner RC Jr, et al. Global temperature constraints on Aedes aegypti and Ae. albopictus persistence and competence for dengue virus transmission. Parasit Vectors. 2014;7: 338. pmid:25052008
  82. 82. Dégallier N, SÁFilho GC, Monteiro HAO, Castro FC, Vaz Da Silva O, Brandão RCF, et al. Release–Recapture Experiments with Canopy Mosquitoes in the Genera Haemagogus and Sabeihes (Diptera: Culicidae) in Brazilian Amazonia. J Med Entomol. 1998;35: 931–936. pmid:9835682
  83. 83. Galindo P. Bionomics of Sabethes chloropterus Humboldt, a vector of sylvan yellow fever in Middle America. Am J Trop Med Hyg. 1958;7: 429–440. pmid:13559598
  84. 84. Mondet B. [Conditions for the laboratory survival of Haemagogus janthinomys dyar, 1921 (Diptera:Culicidae)]. Rev Soc Bras Med Trop. 1997;30: 11–14. pmid:9026827
  85. 85. Couto-Lima D, Madec Y, Bersot MI, Campos SS, Motta M de A, Santos FB dos, et al. Potential risk of re-emergence of urban transmission of Yellow Fever virus in Brazil facilitated by competent Aedes populations. Sci Rep. 2017;7. pmid:28687779
  86. 86. Mascheretti M, Tengan CH, Sato HK, Suzuki A, Souza RP de, Maeda M, et al. Febre amarela silvestre: reemergência de transmissão no estado de São Paulo, Brasil, 2009. Rev Saúde Pública. 2013;47: 881–889. pmid:24626492
  87. 87. Mares-Guia MAM de M, Horta MA, Romano A, Rodrigues CDS, Mendonça MCL, dos Santos CC, et al. Yellow fever epizootics in non-human primates, Southeast and Northeast Brazil (2017 and 2018). Parasit Vectors. 2020;13: 90. pmid:32075684
  88. 88. Teixeira B, Hirsch A, Goulart VDLR, Passos L, Teixeira CP, James P, et al. Good neighbours: distribution of black-tufted marmoset (Callithrix penicillata) in an urban environment. Wildl Res. 2015;42: 579–589.
  89. 89. Jesus JG de, Gräf T, Giovanetti M, Mares-Guia MA, Xavier J, Maia ML, et al. Yellow fever transmission in non-human primates, Bahia, Northeastern Brazil. PLoS Negl Trop Dis. 2020;14: e0008405. pmid:32780745