Skip to main content

Ecological and environmental factors affecting transmission of sylvatic yellow fever in the 2017–2019 outbreak in the Atlantic Forest, Brazil

Abstract

Background

Yellow fever virus (YFV) is an arbovirus that, despite the existence of a safe and effective vaccine, continues to cause outbreaks of varying dimensions in the Americas and Africa. Between 2017 and 2019, Brazil registered un unprecedented sylvatic YFV outbreak whose severity was the result of its spread into zones of the Atlantic Forest with no signals of viral circulation for nearly 80 years.

Methods

To investigate the influence of climatic, environmental, and ecological factors governing the dispersion and force of infection of YFV in a naïve area such as the landscape mosaic of Rio de Janeiro (RJ), we combined the analyses of a large set of data including entomological sampling performed before and during the 2017–2019 outbreak, with the geolocation of human and nonhuman primates (NHP) and mosquito infections.

Results

A greater abundance of Haemagogus mosquitoes combined with lower richness and diversity of mosquito fauna increased the probability of finding a YFV-infected mosquito. Furthermore, the analysis of functional traits showed that certain functional groups, composed mainly of Aedini mosquitoes which includes Aedes and Haemagogus mosquitoes, are also more representative in areas where infected mosquitoes were found. Human and NHP infections were more common in two types of landscapes: large and continuous forest, capable of harboring many YFV hosts, and patches of small forest fragments, where environmental imbalance can lead to a greater density of the primary vectors and high human exposure. In both, we show that most human infections (~ 62%) occurred within an 11-km radius of the finding of an infected NHP, which is in line with the flight range of the primary vectors.

Conclusions

Together, our data suggest that entomological data and landscape composition analyses may help to predict areas permissive to yellow fever outbreaks, allowing protective measures to be taken to avoid human cases.

Graphical Abstract

Background

Yellow fever is an important arboviral disease characterized by febrile and acute symptoms, with high mortality rates. Despite the existence of a safe and effective vaccine since the 1930s, yellow fever virus (YFV) continues to cause outbreaks of varying dimensions in the Americas and Africa [1].

In Brazil, the urban cycle of yellow fever, in which YFV is transmitted between humans through the bite of the domestic mosquito Aedes aegypti, has not been recorded since the early 1940s. However, the sylvatic cycle, where transmission between nonhuman primates (NHPs) occurs by the bite of arboreal mosquitoes, remains active in Brazil and other South American countries [2]. This YFV sylvatic cycle produces epizootic waves of expansion followed by retractions, threatening NHP populations at risk of extinction and affecting humans who live in the border of epizootic forests or who visit them for leisure or work [2,3,4]. The enzootic sylvatic transmission cycle seems to be perennial in the tropical rainforest of northern South America, particularly in the Amazon, from which YFV can initiate epizootic expansion waves affecting other biomes. In Brazil, the YFV expansion waves have spread towards the south and east, where the human population density is much higher than in the Amazon but vaccination coverage is frequently lower [2, 5]. The last YFV epizootic wave started in 2014/2015 in the Amazon and is still active today (2021) in Southern Brazil [6]. It is the largest recorded YFV outbreak. Between 2015 and 2019, YFV crossed the entire southeastern region of Brazil, where the largest and most populous cities infested by Ae. aegypti are located, raising concerns about the risk of re-urbanization of YFV transmission [2, 7,8,9,10].

The virus has also spread to the coastal region of Brazil covered by the Atlantic Forest that had been considered a YFV-free area for the last 80 years and whose human population therefore had even lower vaccination coverage. This was the case of the state of Rio de Janeiro (RJ), where YFV had not been detected since the late 1930s [2, 11]. The state was severely affected during the passage of the 2017–2019 epizootic wave, thus constituting an important place to study the factors involved in the YFV re-emergence. The YFV wave reached RJ in February 2017, entering from its northern portion, which borders Espírito Santo, a state that had already registered YFV infections in humans and NHPs since January 2017 [12, 13]. The virus circulated in several regions of RJ until the last detection in January 2019 [14]. In total, 27 and 262 human cases and 9 and 84 deaths were confirmed for RJ in 2017 and 2018, respectively [15]. In 2019, there was only one confirmed NHP infection with YFV in RJ [14]. Genetic and molecular studies of YFV samples demonstrated that the virus spread from the northeast to southwest of RJ through two different transmission routes, one coastal and one continental, separated by the Serra do Mar mountain chain [9, 13]. Systematic mosquito collections performed during the outbreak showed that the primary mosquito vectors in RJ and other states in southeastern Brazil were Haemagogus janthinomys/capricornii and Hg. leucocelaenus. Three other mosquito species (Sabethes chloropterus, Aedes scapularis, and Ae. taeniorhynchus) were found infected with YFV, but were considered to have a secondary role in transmission in this outbreak [16]. Concerning NHPs, 1177 deaths were recorded in RJ. Howler monkeys (Alouatta guariba clamitans) were proportionally the most highly affected species and those with the highest viral load [17, 18]. Although a considerable amount of information was available on the epidemiological aspects concerning humans and monkeys, the vectors, the virus, and the routes traveled by it in RJ, understanding the ecological and environmental factors that increase the chances of local YFV occurrence and its infection rates is key to improving yellow fever surveillance and control measures.

Indeed, the influence of climatic, environmental, and ecological factors governing the distribution of YFV primary vectors, as well as the dispersion and force of infection of YFV in a naïve area such as RJ, is still poorly known. The increase in rainfall, relative humidity of the air, temperature, and the number of local NHP species have been identified as determining factors for the occurrence of yellow fever and used as elements for the prediction of yellow fever risk in other Brazilian regions [19,20,21]. In the present work we used a multidisciplinary approach combining geographical, ecological, and entomological data to investigate the factors affecting YFV transmission and to better understand the spread of YFV in the landscape mosaic of RJ during the 2017–2019 outbreak.

Methods

Mosquito sampling

Mosquito captures were made according to methods described elsewhere [16]. Briefly, expeditions were performed between 2015 and 2019 (before and during YFV outbreak) in order to investigate fauna composition in different regions of RJ, which constitutes the country's second largest economy and the third largest population (6.75 million inhabitants), and is considered the main country gateway for foreign tourists [22]. Adult mosquitoes were captured through hand nets, mechanical aspirators, and CO2-baited BG-traps, frozen in liquid nitrogen, and subsequently tested by polymerase chain reaction (PCR) for YFV detection [16]. Each mosquito sampling point was georeferenced and classified according to the predominant environment (inside the forest—within dense forests connected to other forests; rural fragment—within forests smaller than 100 ha and surrounded by pastures; rural peri-domicile—around homesteads and country houses; urban fragment—within forests inside cities; and urban intra-domicile—within human houses inside cities) for ecological analyses. Geo-environmental characteristics including altitude, land cover/land use, forest fragment size, and Normalized Difference Vegetation Index (NDVI) [23] were also measured for each sampling point. Collected mosquitoes were classified up to the lowest taxonomic level using dichotomous keys [24,25,26], and data were tabulated according to the sampling point, date, and effort of capture (expressed in numbers of hours, traps, and people working in the capture).

Mosquito functional traits

To analyze the functional diversity of mosquito species at each sampling point, we selected nine categorical parameters (Table 1) related to behavior, physiology, habitats, and epidemiological importance, as follows: 1—oviposition preferences; 2—egg resistance to desiccation; 3—larval development speed; 4—host preference; 5—main hourly biting activity; 6—vertical distribution in the forest; 7—seasonal distribution; 8—main habitat; 9—epidemiological importance concerning YFV. We built a traits × species matrix, and the functional distances between pairs of species were computed using the Gower distance [27]. Functional groups were identified based on a functional dendrogram using the Ward hierarchical agglomerative clustering method [28]. Principal components analysis (PCA) was performed using the package FactoMineR [29] to summarize the relationships among traits and identify those most commonly shared among species from the same functional group.

Table 1 List of functional traits used for each collected mosquito species

Geolocation of YFV infections in humans and NHPs

The geographical coordinates of the most probable local infections (PLI) of humans as well as places where YFV-infected NHPs were found dead were determined based on information obtained from the RJ State Department of Health or through investigations in each of the affected municipalities, supported by the local health departments. To observe the effect of the landscape on the YFV transmission, the coordinates were plotted on a map containing land use/land cover in RJ. In addition, kernel maps were generated by plotting the most probable place of YFV infections in humans, NHPs, and both to verify areas with higher infection forces. Radii of 5 and 11 km, consistent with the flight radius of Haemagogus mosquito species considered the primary vector during the outbreak [30], were plotted from each positive NHP to verify the minimum distance between an epizootic and the findings of human cases. Radii of 5 and 11 km were also plotted from the place where each YFV-positive mosquito was caught to verify the minimum distance between it and the PLI of human cases.

Data analysis

Biodiversity analysis

We built a site × species mosquito abundance matrix containing data from the 84 sampling points and 89 mosquito species sampled. A square root transformation was applied to the raw mosquito species counts to dampen the effect of dominant species, and the transformed counts were then divided by the effort of capture to control for differing sampling efforts. Analyses were performed by comparing (1) the sampling carried out before the yellow fever epidemic with that conducted during the outbreak (the latter was subdivided into sampling points with infected and uninfected mosquitoes), (2) the sampling points sampled during the yellow fever epidemic, and (3) the total period sampled, as described below.

(ii) Comparing the three different scenarios: Differences in mosquito community composition among the three epidemiological scenarios (before YFV outbreak vs. during YFV outbreak positive points vs. during YFV outbreak negative points) were analyzed using a permutational multivariate analysis of variance (PERMANOVA) based on the Bray–Curtis distance matrix [31]. SIMPER (similarity percentage) analysis was applied to assess which species were primarily responsible for the observed differences. These analyses were conducted using the vegan package in R 4.0.3 [32].

Mosquito community diversity was characterized by two measures: (1) mosquito species richness and (2) the Shannon–Wiener index (Table 2) [33]. The functional diversity of mosquito communities was calculated on a continuous scale using four different indices: (1) functional richness (FRic) [34], (2) functional evenness (FEve) [34], (3) functional divergence (FDiv) [34], and (4) functional dispersion (FDis) (Table 2) [35]. The community-weighted means (CWM) of traits [36] were calculated for each sampling point by averaging the trait expression of all mosquito species weighted by their relative abundance. Variation in CWM trait composition among habitat categories was addressed by employing canonical correspondence analysis (CCA). For the CCA, we did not include species traits related to environment (Table 1) in the CWM matrix because we used the habitat descriptions from the sampled areas as environmental factors. Functional diversity analysis and CCA analysis were performed using the FD and vegan packages in R 4.0.3, respectively [32, 37].

Table 2 Variables analyzed in the present study

(ii) Comparing only the sampling points sampled during the yellow fever epidemic (2017–2019)—infection predictors: Generalized linear mixed models (GLMMs) were used to investigate the effects of mosquito biodiversity (richness, Shannon–Wiener index, FRic, FEve, FDiv, and FDis) and Haemagogus relative abundance on (1) the Haemagogus minimum infection rate (MIR) and (2) the positivity of sampling points (Table 2). For that, only the points sampled during YFV outbreak were considered. Haemagogus MIRs were related to the predictors by fitting zero-inflated models with negative binomial errors and log link functions, and the positivity of sampling points by fitting models with binomial errors and logit link functions. A model-averaging approach was applied, which accounts for model uncertainty, increases the robustness of the parameter estimates, and assesses the relative importance of each of the predictor variables [38, 39]. Model averaging started with a global model with all the predictor variables previously described, fitted using the glmmTMB package in R [40]. The dredge function of the MuMIn package [41] was used to create a set of models with all combinations of variables. The Akaike information criterion with correction for small sample size (AICc) [42] was used to identify the best models based on the averaged model including all the equally plausible models (ΔAICc ≤ 2). Averaged parameter estimates were calculated from this set of selected models. We calculated the relative importance of each variable using the model.avg function. The relative importance was calculated using the sum of Akaike weights across all the selected models, with a weight of zero for models where a given parameter was absent [38]. In addition, we calculated the McFadden index (R2M) of the selected models as a measure of model fit.

(ii) Comparing the entire sampled period (2015–2019): As Hg. janthinomys/capricornii and Hg. leucocelaenus are the primary vectors of yellow fever, Pearson correlation tests were carried out to verify whether geo-environmental characteristics (altitude, NDVI, fragment size, land use/land cover) were related to the relative abundance of these species.

Results

Biodiversity analysis

Comparisons among the three different scenarios

Mosquito community composition varied among the three epidemiological scenarios (PERMANOVA, pseudo-F = 1.3258, P-perm = 0.04096). Aedes albopictus, Hg. leucocelaenus, and Hg. janthinomys/capricornii were the species that most contributed to overall dissimilarity between positive and negative sampling points during YFV outbreak, together accounting for 28% of the differences among scenarios.

The CCA results revealed that habitat categories (rural fragment, urban fragment, forest interior, urban intra-domicile, and rural peri-domicile) explained a significant amount of variation in mosquito community mean traits (CWM, ANOVA-like test, 999 permutations; F = 2.69, P = 0.001). The first and second CCA axis explained 30% and 17% of the variance, respectively, and together they accounted for 47% of the explained variation and 12% of total variation in the data (Fig. 1). In other words, traits related to the main yellow fever vectors (e.g. YFV natural infection, YFV experimental vector competence, primatophilic behavior, and skip oviposition) are more common in forest fragments (both rural and urban) and, secondarily, forested habitats (Fig. 1).

Fig. 1
figure 1

Canonical correspondence analysis (CCA) of mosquito trait community-weighted means (CWM) and habitat categories (forest, rural fragment, urban fragment, rural peri-domicile, and urban intra-domicile). Red rectangles highlight traits related to the main yellow fever vectors

The first two axes of the PCA explained 37.7% of the variance (Fig. 2). Altogether, we found four functional groups (Fig. 2), one of which—group 1—contained a concentration of most of the mosquitoes of tribe Aedini, including Haemagogus, primary vectors of yellow fever (Fig. 2a and b). Interestingly, the relative abundance of this functional group was increased at sampling points with YFV-positive mosquitoes (Fig. 2c).

Fig. 2
figure 2

Functional groups formed through mosquito functional trait analysis. a Mosquito species forming each functional group. b Ecological traits forming each functional group. c Relative abundance of each functional group in the three epidemiological scenarios

When comparing ecological aspects among the three epidemiological scenarios (before YFV outbreak, during YFV outbreak negative points, and during YFV outbreak positive points), it became clear that the mosquito biodiversity indicators (richness and Shannon–Wiener index) increased during the outbreak. Furthermore, sampling points with YFV-positive mosquitoes presented higher Haemagogus relative abundance and marginally higher FDis (t = −2.005, df = 25.499, P = 0.055). On the contrary, NDVI measurements were lower at sampling points with positive mosquitoes. The other biodiversity measures were not significantly different between the different scenarios (Fig. 3).

Fig. 3
figure 3

Biodiversity measures considering sampling points in three scenarios: before YFV outbreak vs. during YFV outbreak negative points vs. during YFV outbreak positive points

Comparisons between positive and negative mosquito sampling points during the yellow fever epidemic—infection predictors

For Haemagogus MIRs, model selection indicated that FDis, FRic, richness, and Shannon–Wiener index were the most important predictors, as they were included in all the top-ranked candidate models (proportion of explained deviance R2M around 0.38 and relative importance of overall predictor = 1.00, Additional file 1: Tables S1 and S3). Model-averaged coefficients (based on a 95% confidence interval that excluded 0) for MIRs indicated that an increase in FDis and decrease in FRic, species richness, and Shannon-Weiner diversity led to an increase in the Haemagogus MIR (Table 3). FDis and relative abundance of Haemagogus were also important predictors of the positivity of sampling points, as they were included in most of the top-ranked candidate models (proportion of explained deviance R2M around 0.20 and relative importance of overall predictor > 0.6, Additional file 1: Tables S2 and S4), but their model-averaged coefficients overlapped 0 (Table 4), evidencing a weak effect. Although included among the top-ranked candidate models, FDiv and relative abundance of Haemagogus, and species richness and FEve had a weak effect on Haemagogus MIRs and the positivity of sampling points, respectively (model-averaged coefficients for all covariates overlapped 0, Tables 3 and 4).

Table 3 Model-averaged standardized coefficients (based on models summarized in Additional file 1: Table S1), unconditional standard errors, 95% confidence intervals, and relative importance of biodiversity predictors of Haemagogus MIR in the sampling points during the 2017–2019 YFV outbreak in Rio de Janeiro, Brazil
Table 4 Model-averaged standardized coefficients (based on models summarized in Additional file 1: Table S2), unconditional standard errors, 95% confidence intervals, and relative importance of biodiversity predictors of the positivity of the sampling points during YFV outbreak in Rio de Janeiro, Brazil

Comparisons over the entire sampled period

Among the tested geo-environmental factors (altitude, NDVI, fragment size, land use/land cover), NDVI was the only one that correlated with the abundance of Hg. janthinomys/capricornii (Pearson correlation r = 0.256, P = 0.02) and Hg. leucocelaenus (Pearson correlation r = 0.269, P = 0.014). Altitude, fragment size, and land use/land cover had no direct influence on the abundance of Hg. janthinomys/capricornii (Pearson correlations r = 0.017, P = 0.878; r = −0.007, P = 0.949; r = −0.175, P = 0.113, respectively) or Hg. leucocelaenus (Pearson correlations r = −0.177, P = 0.109; r = −0.001, P = 0.996; r = −0.215, P = 0.051, respectively).

Geospatial analysis

The geographical coordinates of the most probable places of YFV infection for 65 NHPs and 269 humans as well as the 81 mosquito sampling points in RJ were plotted on maps illustrating the predominant land use and vegetation cover throughout the state (Fig. 4). YFV-infected mosquitoes were found in six out of the 26 sampling points sampled during the yellow fever outbreak (Fig. 4). As expected, human and NHP YFV infections were abundant in areas with forest cover and forest fragments, but practically absent in areas predominantly covered by pastures and highly dense urbanized areas (Fig. 3). The kernel maps (Fig. 5) showed that, despite YFV infections in both NHPs (Fig. 5a) and humans (Fig. 5b) being distributed throughout almost half of the state, it is possible to identify at least three clusters with a greater number of records (Fig. 5c). Interestingly, two of these clusters (rectangles d and e) coincide with the most heavily forested areas, located in the Bahia da Ilha Grande and Serrana regions, and the third (rectangle f) is composed of small, highly fragmented forest, located in the Médio Paraíba and Centro-Sul regions (Fig. 4 and 5c).

Fig. 4
figure 4

Rio de Janeiro state map showing land use/land cover along with mosquito sampling points and distribution of most probable places of infection for NHPs and humans infected with YFV

Fig. 5
figure 5

Kernel maps showing YFV infection density in Rio de Janeiro state, Brazil in 2015–2019. Areas covered by forests are also shown. a YFV infection density in NHP. b YFV infection density in humans. c YFV infection density in both NHP and humans. Rectangles d, e and f highlight three YFV infection clusters. Numbers indicate the regions of Rio de Janeiro: 1—Bahia da Ilha Grande, 2—Médio Paraíba, 3—Metropolitana I, 4—Centro-Sul, 5—Serrana, 6—Metroplotana II, 7—Baixada Litorânea, 8—Norte, 9—Noroeste

We found that most of the YFV confirmed human cases were concentrated within radii of 5 and 11 km (compatible with the flight range of the primary vectors) around the PLI of YFV infection in NHPs (Fig. 6a) or in mosquitoes (Fig. 6b). In fact, 62% (n = 166) of the human cases were within a radius of 11 km of a YFV-positive NHP (Fig. 7a) and 90% were within 25 km (Fig. 7b).

Fig. 6
figure 6

Maps showing 5- and 11-km radius from each point where an infected NHP (a) or mosquito (b) was found in Rio de Janeiro state, Brazil in 2015–2019. Confirmed human cases are shown as black dots

Fig. 7
figure 7

a Number of YFV human cases within a distance of 11 km from a YFV-infected NHP. b Cumulative percentage of YFV human cases as a function of the distance from a YFV-infected NHP in Rio de Janeiro state, Brazil in 2015–2019

Discussion

Despite several studies investigating different aspects of the recent YFV outbreak in southeast Brazil, this is the first time that ecological and geographical analysis tools were used to shed light on the determinants of the spread and transmission strength of the virus in RJ, a state without YFV circulation for around 80 years. From the standpoint of mosquito biodiversity, we demonstrate that the increase in Haemagogus abundance and functional dispersion, together with the reduction in mosquito species richness and FEve associated with the increased abundance of functional group 1, seems to increase the risk of YFV infection in mosquitoes. From a geographical perspective, we show that YFV transmission was more intense in two different types of environments: in continuous forest areas and in mosaics composed of forest fragments surrounded by pastures. Furthermore, we show that the majority of PLI in humans are within a radius of 11 km from confirmed YFV epizootics in NHPs. Together, these data help us understand the complexity of the factors determining YFV dispersion and provide support for further prediction risk analyses and preventive measures which could protect human populations.

When we compared the three epidemiological scenarios (before YFV outbreak, during YFV outbreak positive points, and during YFV outbreak negative points), we found an increase in mosquito species richness, Shannon-Winner index, and Haemagogus abundance during the yellow fever outbreak. In other words, the increase in species richness offset the increase in Haemagogus abundance and increased the Shannon–Wiener index. Furthermore, the increased presence and abundance of the two species considered primary vectors of YFV in the southeastern outbreak (Hg. janthinomys/capricornii and Hg. leucocelaenus) [16, 43], combined with those of the potential vector Ae. albopictus, constituted the main differences between areas with and without detection of yellow fever in mosquitoes. Natural YFV infection has never been confirmed through viral isolation or complete genome sequencing in Ae. albopictus in Brazil, and Brazilian populations of this species have low vector competence for YFV [18, 44]. However, it can bite YFV-infected NHP on the forest ground layer [14], and YFV has the potential to adapt to this species [45], which is why it can be considered a potential vector [46]. Haemagogus janthinomys/capricornii and Hg. leucocelaenus, as well as Ae. albopictus, are tree-hole breeders, and lay eggs capable of resisting months of desiccation, hatching after one (for Ae. albopictus) or several (for Haemagogus) immersions in water [47,48,49]. Therefore, they are species that depend on and respond to similar environmental conditions. Concerning climate, their population peak is recorded in the rainy summer characterized by increased rainfall and temperature [24, 50], which also determines the seasonal period of yellow fever transmission in Brazil [51].

Regarding the functional diversity of mosquitoes, it is interesting to note that some characteristics are more abundant in forest environments (e.g. breeding sites in plant axils and in shallow pools on the ground, as well as solitary larvae breeders) and others in urban environments (e.g. artificial breeding sites). A greater number of traits are also found in forested areas, indicating greater diversity of ecological niches. Importantly, traits related to the biology of YFV vectors (e.g. vector competence, natural infection, skip oviposition) are striking among forest environments and in rural and urban forest fragments, which reflects the capacity of these environments to support populations of these vectors. Currently, the expansion of cities into forest fragments for estate speculation and/or the use of these fragments for leisure activities has been observed, which increases the chances of human infections in the sylvatic cycle and raises concerns regarding the risk of re-urbanization of yellow fever, especially in large cities surrounded by or interspersed with green spots [2]. São Paulo, Goiânia, and Nova Iguaçu are examples of large Brazilian cities with recent detection of yellow fever in their nearby forest areas, which increases the need for constant surveillance [6, 52, 53]. From the functional diversity analyses, it was possible to describe, for the first time, the existence of four main functional groups within the mosquito communities. The increased abundance of group 1, formed mainly by Aedini mosquitoes including Haemagogus species, primary vectors of yellow fever, in areas where positive mosquitoes were found may be an important indicator of areas receptive to viral circulation. Therefore, it will be important to carry out studies on the abundance of these functional groups to detect new areas receptive to YFV.

Curiously, although an increase in mosquito species richness and diversity was detected at the points sampled during the yellow fever outbreak, the Haemagogus MIRs showed the opposite behavior. That is, when comparing positive areas, we found that Haemagogus infection rates increased when mosquito diversity decreased and when relative abundance of Haemagogus increased, which could be explained by the concept of the “dilution effect.” The dilution effect describes the idea that the presence of different host species would increase the evenness of species abundance and dilute the chances of pathogen transmission, due to the presence of several incompetent reservoirs [54]. In this sense, in a more diverse region (greater number of vertebrate host and mosquito species), the probability of transmission of YFV is reduced due to greater availability of hosts—many of which are incapable of amplifying the virus—and competition for resources among vectors, which can be eclectic concerning blood-feeding sources [55]. Interestingly, we observed that the reduction in mosquito diversity, expressed by the Shannon–Wiener, richness, and FRic indices, positively influenced the Haemagogus MIR, corroborating the dilution effect theory. Similar results have been found for other mosquito-borne diseases. Decreasing bird and mosquito diversity, for example, has enhanced the incidence of human West Nile fever cases in the United States, and reduced functional bird diversity has increased the prevalence of avian malaria and the diversity of Plasmodium strains in the Brazilian Atlantic Forest [56,57,58]. In conclusion, the most important predictors of infection were the high relative abundance of Haemagogus, and consequently lower mosquito species richness and diversity. Areas and situations with these characteristics should be monitored by propensity for viral circulation. Nevertheless, we could not establish a Haemagogus population threshold that indicated major or imminent risk of YFV circulation in RJ, where the virus caused an outbreak and has disappeared since 2019. Cross-sectional studies in enzootic/endemic areas, such as the Amazon region, may help to unravel this issue.

The type of land cover/land use is highlighted as a determining factor for the occurrence of sylvatic yellow fever. Forested areas are more susceptible to viral circulation due to the higher probability of occurrence of NHP species and primary vector mosquito species that use water in tree holes for breeding [19, 21, 59, 60]. In fact, we found two important clusters of confirmed YFV infections in humans and NHPs in the areas most densely covered by forests, which in RJ correspond to the large forest continuum distributed along the mountain chain of Serra do Mar (Fig. 4 and 5c, rectangles d and e). However, a third cluster was verified in the Médio Paraíba and Centro-Sul regions of RJ (Fig. 5c, rectangle f). It is a mosaic composed of small forest fragments mostly surrounded by pastures (Fig. 5c, rectangle f). This type of landscape, characterized by intermediate levels of forest cover with numerous fragments producing higher availability of forest edges, was recently identified as more prone to the occurrence of human cases in a municipality-level analysis [61]. In this type of environment, the spread of the virus between the fragments possibly occurs through wind currents that facilitate and enhance the flight of vector mosquitoes [20, 30, 62]. Furthermore, humans become infected when approaching the fragments to rest, plant, harvest, or extract other forest resources. Recently, a correlation between the seasonality of agriculture and the appearance of YFV human cases in Brazil was demonstrated [63]. During planting and harvesting times, people working in agriculture and extractivism, who represent about 45% of the reported YFV cases, are more exposed to contact with the wild environment [63], which together with the dilution effect would help to explain the infection force in these regions with small forest fragments. Therefore, it appears that there are two distinct socio-ecological contexts, both favorable to the circulation of the YFV: The first is continuous forest environments that, due to the diversity and abundance of ecological niches, have a great support capacity for harboring many specimens of different species of NHP and of mosquitoes. In this context, humans would be exposed when entering the woods for tourism, leisure, or in search of bucolic moments nearby [2, 64], but would be partially protected by the dilution effect. The second is fragmented forest mosaics surrounded by pastures or plantations. Although small, these forest fragments can support groups of Alouatta and Haemagogus mosquitoes, the main vertebrate and invertebrate hosts, respectively. [16, 65,66,67,68]. In this context, humans, especially rural workers, would become infected when approaching the fragments during their work routine.

To the best of our knowledge, this is the first time that the minimum distance has been measured between the identification of a YFV-infected NHP and the occurrence of human infections. The 11-km radius, where 62% of human cases are concentrated, matches the flight capacity of the main vectors identified in RJ [16], nearly 6 km for Hg. leucocelaenus and 11 km for Hg. janthinomys/capricornii [30]. It is important to consider that finding a dead monkey and confirming it for YFV is not trivial, as it depends on the capacity of surveillance combined with the opportunity to collect viable samples for diagnosis [69]. The molecular analysis of a large set of YFV samples obtained from the 2017 and 2019 outbreak [9, 13, 70] has proved that the virus reached RJ from southern Espírito Santo. However, the large territorial gap in reports of YFV-infected NHPs between the southernmost border of Espírito Santo and the nearest NHP and human depicted in Fig. 3 highlights the failure of YFV epizootic surveillance. Therefore, it is possible that most of the YFV epizootics have not been detected and reported and/or laboratory-confirmed, which probably contributed to the increase in the distances measured. Even so, according to our data, strategic vaccination of the population within a radius of up to 25 km from a confirmed epizooty disease would prevent more than 90% of human cases, even in areas without a record of epizootic disease, considered epidemiological silent areas. Timely vaccination within a 2-km radius of a dead monkey helped to prevent several human cases in the 2008–2009 outbreak in southern Brazil [71]. The strategy of rapid and mobile vaccination campaigns to reach the more vulnerable rural populations within a radius of up to 25 km from the encounter of an infected mosquito may also be applied, especially when the surveillance of epizooties is failing. Although the opportunity for finding infected mosquito collections is even rarer, entomological and primatological surveys are complementary. In areas such as the city of Teresópolis, for example, where 22 human cases were confirmed and georeferenced, only positive mosquitoes were found, due to a lack of timely collection of NHP [15].

Conclusions

Together, our results add new pieces to the puzzle towards understanding the epidemiology of sylvatic yellow fever in the Brazilian Atlantic Forest from data obtained prior to and during the outbreak in RJ. Importantly, we showed that there are at least two sets of favorable conditions for the circulation of yellow fever: areas with large forest continuums, capable of harboring many hosts of the virus, and areas with small forest fragments, where environmental imbalance can lead to a greater density of primary vectors, linking the loss of biodiversity as a risk factor. For both sets, changes in the relative abundance of functional groups and species composition of vector communities seem to determine not only the possibility of YFV circulation, but also the rate of vector infection and the occurrence of human cases. In this way, future analyses of vector biodiversity and landscape may help to predict areas permissive to yellow fever outbreaks.

Availability of data and materials

The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.

References

  1. 1.

    Monath TP, Vasconcelos PFC. Yellow fever. J Clin Virol. 2015;64:160–73.

    PubMed  Google Scholar 

  2. 2.

    Possas C, Lourenço-de-Oliveira R, Tauil PL, Pinheiro FDP, Pissinatti A, Venâncio R, et al. Yellow fever outbreak in Brazil : the puzzle of rapid viral spread and challenges for immunisation. Mem Inst Oswaldo Cruz. 2018;113:1–12.

    PubMed  Google Scholar 

  3. 3.

    Romano APM, Costa ZGA, Ramos DG, Andrade MA, de Jayme VS, de Almeida MAB, et al. Yellow fever outbreaks in unvaccinated populations, Brazil, 2008–2009. PLoS Negl Trop Dis. 2014;8:e2740.

    PubMed  PubMed Central  Google Scholar 

  4. 4.

    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(1):12926.

    PubMed  PubMed Central  Google Scholar 

  5. 5.

    Vasconcelos PFC, Costa ZG, Travassos da Rosa ES, Luna E, Rodrigues SG, Barros VLRS, et al. Epidemic of jungle yellow fever in Brazil, 2000: implications of climatic alterations in disease spread. J Med Virol. 2001;65:598–604.

    CAS  PubMed  Google Scholar 

  6. 6.

    Brasil. Boletim Epidemiológico—Situação epidemiológica da febre amarela—Monitoramento 2020/2021. 2021. https://www.gov.br/saude/pt-br/media/pdf/2021/fevereiro/11/boletim_epidemiologico_svs_4.pdf. Acessed 10 Oct 2021.

  7. 7.

    Silva NIO, Sacchetto L, De Rezende IM, Trindade GDS, 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.

    PubMed  PubMed Central  Google Scholar 

  8. 8.

    Faria NR, Kraemer MUG, Hill SC, de Jesus JG, Aguiar RS, Iani FCM, et al. Genomic and epidemiological monitoring of yellow fever virus transmission potential. Science. 2018;361:894–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Gómez MM, de Abreu FVS, dos Santos AAC, de Mello IS, Santos MP, Ribeiro IP, et al. Genomic and structural features of the yellow fever virus from the 2016–2017 Brazilian outbreak. J Gen Virol. 2018;99:536–48.

    PubMed  Google Scholar 

  10. 10.

    Cunha MS, da Costa AC, de Azevedo Fernandes NCC, Guerra JM, dos Santos FCP, Nogueira JS, et al. Epizootics due to yellow fever virus in São Paulo state, Brazil: viral dissemination to new areas (2016–2017). Sci Rep. 2019;9:5474.

    PubMed  PubMed Central  Google Scholar 

  11. 11.

    Franco O. História da Febre Amarela no Brasil. Ministério da Saúde. Rio de Janeiro: Departamento Nacional de Endemias Rurais; 1969.

    Google Scholar 

  12. 12.

    Bonaldo MC, Gómez MM, dos Santos AA, de Abreu FVS, Ferreira-de-Brito A, de Miranda RM, et al. Genome analysis of yellow fever virus of the ongoing outbreak in Brazil reveals polymorphisms. Mem Inst Oswaldo Cruz. 2017;112:447–51.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Delatorre E, de Abreu FVS, Ribeiro IP, Gómez MM, dos Santos AAC, Ferreira-de-Brito A, et al. Distinct YFV lineages co-circulated in the Central-Western and Southeastern Brazilian regions from 2015 to 2018. Front Microbiol. 2019;10:1–12.

    Google Scholar 

  14. 14.

    de Abreu FVS, Delatorre E, dos Santos AAC, Ferreira-de-Brito A, de Castro MG, Ribeiro IP, et al. Combination of surveillance tools reveals that yellow fever virus can remain in the same Atlantic Forest area at least for three transmission seasons. Mem Inst Oswaldo Cruz. 2019;114:e190076.

    PubMed  PubMed Central  Google Scholar 

  15. 15.

    Secretaria de Saúde de Estado do Rio de Janeiro. Dados—Febre Amarela RJ. 2021. http://www.febreamarela.rj.gov.br/dados. Acessed 3 May 2021.

  16. 16.

    de Abreu FVS, Ribeiro IP, Ferreira-de-Brito A, dos Santos AAC, de Miranda RM, de Bonelly IS, 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–31.

    PubMed  PubMed Central  Google Scholar 

  17. 17.

    de Mares-Guia MAMM, 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.

    PubMed  PubMed Central  Google Scholar 

  18. 18.

    de Abreu FVS, Ferreira-de-Brito A, Azevedo ADS, Linhares JHR, de Oliveira SV, Hime Miranda E, et al. Survey on non-human primates and mosquitoes does not provide evidences of spillover/spillback between the urban and sylvatic cycles of yellow fever and zika viruses following severe outbreaks in Southeast Brazil. Viruses. 2020;12:364.

    PubMed Central  Google Scholar 

  19. 19.

    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. PLoS Negl Trop Dis. 2017;11:e0005897.

    PubMed  PubMed Central  Google Scholar 

  20. 20.

    de Almeida MAB, dos Santos E, da Cardoso JC, da Silva LG, Rabelo RM, Bicca-Marques JC. Predicting yellow fever through species distribution modeling of virus, vector, and monkeys. EcoHealth. 2018;16:95–108.

    PubMed  Google Scholar 

  21. 21.

    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.

    Google Scholar 

  22. 22.

    IBGE. IBGE divulga estimativa da população dos municípios para 2020. 2020 https://agenciadenoticias.ibge.gov.br/agencia-sala-de-imprensa/2013-agencia-de-noticias/releases/28668-ibge-divulga-estimativa-da-populacao-dos-municipios-para-2020. Acessed 4 May 2021.

  23. 23.

    Myneni RB, Hall FG, Sellers PJ, Marshak AL. The interpretation of spectral vegetation indexes. IEEE Trans Geosci Remote Sens. 1995;33:481–6.

    Google Scholar 

  24. 24.

    Consoli RAGB, de Oliveira RL. Principais mosquitos de importância sanitária no Brasil. Rio de Janeiro: Editora Fiocruz; 1994.

    Google Scholar 

  25. 25.

    Lane J, Cerqueira NL. Os Sabetíneos da América, (Diptera, Culicidae). Arquivos de Zoologia do Estado de São Paulo. 1942;3:473–849.

    Google Scholar 

  26. 26.

    Forattini OP. Culicidologia médica: identificação, biologia, epidemiologia, vol. 2. São Paulo: Edusp; 2002.

    Google Scholar 

  27. 27.

    Pavoine S, Vallet J, Dufour A-B, Gachet S, Daniel H. On the challenge of treating various types of variables: application for improving the measurement of functional diversity. Oikos. 2009;118:391–402.

    Google Scholar 

  28. 28.

    Murtagh F, Legendre P. Ward’s hierarchical agglomerative clustering method: which algorithms implement Ward’s criterion? J Classif. 2014;31:274–95.

    Google Scholar 

  29. 29.

    Lê S, Josse J, Husson F. FactoMineR: an R package for multivariate analysis. J Stat Softw. 2008;25:1–18.

    Google Scholar 

  30. 30.

    Causey OR, Laemmert HW, Kumm HW. Dispersion of forest mosquitoes in Brazil: further studies 1. Am J Trop Med Hyg. 1950;s1–30:301–12.

    Google Scholar 

  31. 31.

    Anderson MJ. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001;26:32–46.

    Google Scholar 

  32. 32.

    Oksanen AJ, Blanchet FG, Friendly M, Kindt R, Legendre P, Mcglinn D, et al. Package vegan: community ecology package. R package version 4.0.3. 2019. https://cran.r-project.org/web/packages/vegan/vegan.pdf. Acessed 10 Apr 2021.

  33. 33.

    Magurran A. Measuring biologcial diversity. Oxford: Wiley-Blackwell; 2004.

    Google Scholar 

  34. 34.

    Villéger S, Mason NWH, Mouillot D. New multidimensional functional diversity indices for a multifaceted framework in functional ecology. Ecology. 2008;89:2290–301.

    PubMed  Google Scholar 

  35. 35.

    Laliberté E, Legendre P. A distance-based framework for measuring functional diversity from multiple traits. Ecology. 2010;91:299–305.

    PubMed  Google Scholar 

  36. 36.

    Kleyer M, Dray S, Bello F, Lepš J, Pakeman RJ, Strauss B, et al. Assessing species and community functional responses to environmental gradients: which multivariate methods? J Veg Sci. 2012;23:805–21.

    Google Scholar 

  37. 37.

    Laliberté, E., Legendre, P. and BS. FD: measuring functional diversity from multiple traits, and other tools for functional ecology. R package version 1.0-12. 2014. https://cran.r-project.org/web/packages/FD/FD.pdf. Acessed 10 Apr 2021.

  38. 38.

    Burnham K, Anderson D. Model selection and multimodel inference—a practical information-theoretic approach. 2nd ed. New York: Springer; 2002.

    Google Scholar 

  39. 39.

    Grueber CE, Nakagawa S, Laws RJ, Jamieson IG. Multimodel inference in ecology and evolution: challenges and solutions. J Evol Biol. 2011;24:699–711.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Brooks ME, Kristensen K, van Benthem KJ, Magnusson A, Berg CW, Nielsen A, et al. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R J. 2017;9:378–400.

    Google Scholar 

  41. 41.

    Barton MK. MuMIn: Multi-Model Inference. R package version 1.43.17. 2020. https://cran.r-project.org/web/packages/MuMIn/MuMIn.pdf. Acessed 10 Apr 2021.

  42. 42.

    Hurvich CM, Tsai C-L. Regression and time series model selection in small samples. Biometrika. 1989;76:297–307.

    Google Scholar 

  43. 43.

    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.

    PubMed Central  Google Scholar 

  44. 44.

    Couto-Lima D, Madec Y, Bersot MI, Campos SS, de Motta MA, dos Santos FB, 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:4848.

    PubMed  PubMed Central  Google Scholar 

  45. 45.

    Amraoui F, Pain A, Piorkowski G, Vazeille M, Couto-Lima D, de Lamballerie X, et al. Experimental adaptation of the yellow fever virus to the mosquito Aedes albopictus and potential risk of urban epidemics in Brazil, South America. Sci Rep. 2018;8:14337.

    PubMed  PubMed Central  Google Scholar 

  46. 46.

    Pereira dos Santos T, Roiz D, Santos de Abreu FV, Luz SLB, Santalucia M, Jiolle D, et al. Potential of Aedes albopictus as a bridge vector for enzootic pathogens at the urban-forest interface in Brazil. Emerg Microbes Infect. 2018;7:191.

    PubMed  PubMed Central  Google Scholar 

  47. 47.

    Tátila-Ferreira A, de Maia DA, Alencar J. Development of preimaginal stages of Haemagogus leucocelaenus (Diptera: Culicidae) in laboratory conditions. Entomol News. 2017;127:142–50.

    Google Scholar 

  48. 48.

    Alencar J, Gleiser RM, Morone F, de Mello CF, dos Silva JS, Serra-Freire NM, et al. A comparative study of the effect of multiple immersions on Aedini (Diptera: Culicidae) mosquito eggs with emphasis on sylvan vectors of yellow fever virus. Mem Inst Oswaldo Cruz. 2014;109:114–7.

    PubMed  Google Scholar 

  49. 49.

    Alencar J, de Almeida HM, Marcondes CB, Guimarães AÉ. Effect of multiple immersions on eggs and development of immature forms of Haemagogus janthinomys from South-Eastern Brazil (Diptera: Culicidae). Entomol News. 2008;119:239–45.

    Google Scholar 

  50. 50.

    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:e200218.

    PubMed  PubMed Central  Google Scholar 

  51. 51.

    Romano APM, Ramos DG, Araújo FAA, de Siqueira GAM, Ribeiro MPD, Leal SG, et al. Febre amarela no Brasil: recomendações para a vigilância, prevenção e controle. Epidemiol e Serviços Saúde. 2011;20:101–6.

    Google Scholar 

  52. 52.

    Mucci LF, Medeiros-Sousa AR, Ceretti-Júnior W, Fernandes A, Camargo AA, Evangelista E, et al. Haemagogus leucocelaenus and other mosquitoes potentially associated with sylvatic yellow fever in Cantareira State Park in the São Paulo metropolitan area, Brazil. J Am Mosq Control Assoc. 2016;32:329–32.

    PubMed  Google Scholar 

  53. 53.

    Cunha MS, Tubaki RM, de Menezes RMT, Pereira M, Caleiro GS, Coelho E, et al. Possible non-sylvatic transmission of yellow fever between non-human primates in São Paulo city, Brazil, 2017–2018. Sci Rep. 2020;10:15751.

    CAS  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Schmidt KA, Ostfeld RS. Biodiversity and the dilution effect in disease ecology. Ecology. 2001;82:609–19.

    Google Scholar 

  55. 55.

    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 Incl Trop Dis. 2015;21:6.

    PubMed  PubMed Central  Google Scholar 

  56. 56.

    Swaddle JP, Calos SE. Increased avian diversity is associated with lower incidence of human West Nile infection: observation of the dilution effect. PLoS ONE. 2008;3:e2488.

    PubMed  PubMed Central  Google Scholar 

  57. 57.

    Fecchio A, Lima MR, Bell JA, Schunck F, Corrêa AH, Beco R, et al. Loss of forest cover and host functional diversity increases prevalence of avian malaria parasites in the Atlantic Forest. Int J Parasitol. 2021;51:719–28.

    PubMed  Google Scholar 

  58. 58.

    Johnson BJ, Munafo K, Shappell L, Tsipoura N, Robson M, Ehrenfeld J, et al. The roles of mosquito and bird communities on the prevalence of West Nile virus in urban wetland and residential habitats. Urban Ecosyst. 2012;15:513–31.

    PubMed  PubMed Central  Google Scholar 

  59. 59.

    De Paiva CA, Oliveira APDS, Muniz SS, Calijuri ML, Dos Santos VJ, Alves SDC. Determination of the spatial susceptibility to yellow fever using a multicriteria analysis. Mem Inst Oswaldo Cruz. 2019;114:e180509.

    PubMed  PubMed Central  Google Scholar 

  60. 60.

    Wilk-da-Silva R, Mucci LF, Ceretti-Junior W, de Duarte AMRC, Marrelli MT, Medeiros-Sousa AR. Influence of landscape composition and configuration on the richness and abundance of potential sylvatic yellow fever vectors in a remnant of Atlantic Forest in the city of São Paulo, Brazil. Acta Trop. 2020;204:105385.

    PubMed  Google Scholar 

  61. 61.

    Ilacqua RC, Medeiros-Sousa AR, Ramos DG, Obara MT, Ceretti-Junior W, Mucci LF, et al. Reemergence of yellow fever in Brazil: the role of distinct landscape fragmentation thresholds. J Environ Public Health. 2021;2021:8230789.

    PubMed  PubMed Central  Google Scholar 

  62. 62.

    Causey OR, Hughes TP, Laemmert HW. The invasion of small forests by yellow fever virus as indicated by immunity in Cebus monkeys 1. Am J Trop Med Hyg. 1949;S1-29:555–65.

    Google Scholar 

  63. 63.

    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.

    CAS  PubMed  PubMed Central  Google Scholar 

  64. 64.

    Hamer DH, Angelo K, Caumes E, van Genderen PJJ, Florescu SA, Popescu CP, et al. Fatal yellow fever in travelers to Brazil, 2018. Morb Mortal Wkly Rep. 2018;67:340–1.

    Google Scholar 

  65. 65.

    Chiarello AG, de Melo FR. Primate population densities and sizes in Atlantic forest remnants of northern Espirito Santo, Brazil. Int J Primatol. 2001;22:379–96.

    Google Scholar 

  66. 66.

    Martins MM. Density of primates in four semi-deciduous forest fragments of São Paulo, Brazil. Biodivers Conserv. 2005;14:2321–9.

    Google Scholar 

  67. 67.

    Abreu FVS, dos Santos E, Gomes MQ, Vargas WP, Oliveira Passos PH, Nunes e Silva C, et al. Capture of Alouatta guariba clamitans for the surveillance of sylvatic yellow fever and zoonotic malaria: which is the best strategy in the tropical Atlantic Forest? Am J Primatol. 2019;81:e23000.

    PubMed  Google Scholar 

  68. 68.

    Behie AM, Pavelka MS. Primates in fragments: complexity and resilience. New York: Springer; 2013.

    Google Scholar 

  69. 69.

    Brasil. Guia de vigilância de epizootias em primatas não humanos e entomologia aplicada à vigilância da febre amarela. 2nd ed. Brasília: Ministério da Saúde; 2017.

  70. 70.

    Giovanetti M, de Mendonça MCL, 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:e01623-e1719.

    PubMed  PubMed Central  Google Scholar 

  71. 71.

    Almeida MAB, da Cardoso JC, 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. PLoS Negl Trop Dis. 2014;8:e2741.

    PubMed  PubMed Central  Google Scholar 

  72. 72.

    Galindo P. Bionomics of Sabethes chloropterus Humboldt, a vector of sylvan yellow fever in Middle America. Am J Trop Med Hyg. 1958;7:429–40.

    CAS  PubMed  Google Scholar 

  73. 73.

    Farnesi LC, Menna-Barreto RFS, Martins AJ, Valle D, Rezende GL. Physical features and chitin content of eggs from the mosquito vectors Aedes aegypti, Anopheles aquasalis and Culex quinquefasciatus: connection with distinct levels of resistance to desiccation. J Insect Physiol. 2015;83:43–52.

    CAS  PubMed  Google Scholar 

  74. 74.

    Day J. Mosquito oviposition behavior and vector control. Insects. 2016;7:65.

    PubMed Central  Google Scholar 

  75. 75.

    Arnell JH. Mosquito studies (Diptera, Culicidae) XXXII. A revision of the genus Haemagogus. Contrib Am Entomol Inst. 1973;10:1–174.

    Google Scholar 

  76. 76.

    Lounibos LP, Machado-Allison CE. Oviposition and egg brooding by the mosquito Trichoprosopon digitatum in cacao husks. Ecol Entomol. 1983;8:475–8.

    Google Scholar 

  77. 77.

    Tátila-Ferreira A, de Maia DA, 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 Sao Paulo. 2017;59:e60.

    PubMed  PubMed Central  Google Scholar 

  78. 78.

    Shannon RC. The environment and behavior of some Brazilian mosquitoes. Proc Ent Soc Wash. 1931;33:1–27.

    Google Scholar 

  79. 79.

    Bates M. The natural history of mosquitoes. New York: The Macmillan Company; 1949.

    Google Scholar 

  80. 80.

    Bradshaw WE, Holzapfel CM, Kleckner CA, Hard JJ. Heritability of development time and protandry in the pitcher-plant mosquito Wyeomyia smithii. Ecology. 1997;78:969–76.

    Google Scholar 

  81. 81.

    Lima MIB, Motta MA, Lourenço De Oliveira R. Aspectos bionômicos de Sabethes (Peytonulus) identicus Dyar & Knab, 1907 (Diptera: Culicidae) em laboratório. 2014. https://www.arca.fiocruz.br/bitstream/icict/34787/2/maria_bersot_ioc_espec_2014.pdf. Acessed 11 Jun 2021.

  82. 82.

    Lourenço-de-Oliveira R. Alguns aspectos da ecologia dos mosquitos (Diptera: Culicidae) de uma área de planície (Granjas Calábria), em Jacarepaguá, Rio de Janeiro. I. frequência comparativa das espécies em diferentes ambientes e métodos de coleta. Mem Inst Oswaldo Cruz. 1984;79:479–90.

    Google Scholar 

  83. 83.

    Lourenço-de-Oliveira R, Heyden R. Alguns aspectos da ecologia dos mosquitos (Diptera: Culicidae) de uma área de planície (Granjas Calábria) em Jacarepaguá, Rio de Janeiro: IV. Preferências alimentares quanto ao hospedeiro e freqüência domiciliar. Mem Inst Oswaldo Cruz. 1986;81:15–27.

    Google Scholar 

  84. 84.

    Guimarães AÉ, Arlé M, Machado RNM. Mosquitos no Parque Nacional da Serra dos Orgãos, estado do Rio de Janeiro, Brasil: IV. Preferência alimentar. Mem Inst Oswaldo Cruz. 1987;82:277–85.

    Google Scholar 

  85. 85.

    Alencar J, Lorosa ES, Dégallier N, Serra-Freire NM, Pacheco JB, Guimarães AÉ. Feeding patterns of Haemagogus janthinomys (Diptera: Culicidae) in different regions of Brazil. J Med Entomol. 2005;42:981–5.

    PubMed  Google Scholar 

  86. 86.

    Nielsen HT. Swarming and some other habits of Mansonia perturbans and Psorophora ferox (Diptera: Culicidae). Behaviour. 1964;24:67–88.

    CAS  PubMed  Google Scholar 

  87. 87.

    Lourenço-de-Oliveira R, da Silva TF, Heyden R. Alguns aspectos da ecologia dos mosquitos (Diptera: Culicidae) de uma área de planície (Granjas Calábria), em Jacarepaguá, Rio de Janeiro: II. Frequência mensal e no ciclo lunar. Mem Inst Oswaldo Cruz. 1985;80:123–33.

    Google Scholar 

  88. 88.

    Deane LM, Ferreira Neto JA, Lima MM. The vertical dispersion of Anopheles (Kerteszia) cruzi in a forest in Southern Brazil suggests that human cases of malaria of simian origin might be expected. Mem Inst Oswaldo Cruz. 1984;79:461–3.

    CAS  PubMed  Google Scholar 

  89. 89.

    Guimarães AE, Arlé M, Machado RNM. Mosquitos no Parque Nacional da Serra dos Órgãos, estado do Rio de Janeiro, Brasil. II. Distribuição vertical. Mem Inst Oswaldo Cruz. 1985;80:171–85.

    Google Scholar 

  90. 90.

    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 Sabethes (Diptera: Culicidae) in Brazilian Amazonia. J Med Entomol. 1998;35:931–6.

    PubMed  Google Scholar 

  91. 91.

    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–5.

    PubMed  Google Scholar 

  92. 92.

    Lira-Vieira AR, Gurgel-Goncalves R, Moreira IM, Yoshizawa MAC, Coutinho ML, Prado PS, et al. Ecological aspects of mosquitoes (Diptera: Culicidae) in the gallery forest of Brasilia National Park, Brazil, with an emphasis on potential vectors of yellow fever. Rev Soc Bras Med Trop. 2013;46:566–74.

    PubMed  Google Scholar 

  93. 93.

    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.

    CAS  PubMed  PubMed Central  Google Scholar 

  94. 94.

    Pereira-Silva JW, Ríos-Velásquez CM, de Lima GR, MarialvadosSantos EF, Belchior HCM, Luz SLB, et al. Distribution and diversity of mosquitoes and Oropouche-like virus infection rates in an Amazonian rural settlement. PLoS ONE. 2021;16:e0246932.

    CAS  PubMed  PubMed Central  Google Scholar 

  95. 95.

    Galindo P, Carpenter SJ, Trapido H. A contribution to the ecology and biology of tree hole breeding mosquitoes of Panama1. Ann Entomol Soc Am. 1955;48:158–64.

    Google Scholar 

  96. 96.

    Galindo P, Trapido H, Carpenter SJ, Blanton FS. The abundance cycles of arboreal mosquitoes during six years at a sylvan yellow fever locality in Panama. Ann Entomol Soc Am. 1956;49:543–7.

    Google Scholar 

  97. 97.

    Chadee DD. Seasonal abundance and diel landing periodicity of Sabethes chloropterus (Diptera: Culicidae) in Trinidad, West Indies. J Med Entomol. 1990;27:1041–4.

    CAS  PubMed  Google Scholar 

  98. 98.

    Lourenço-de-Oliveira R, Castro MG, Braks MAH, Lounibos LP. The invasion of urban forest by dengue vectors in Rio de Janeiro. J Vector Ecol. 2004;29:94–100.

    PubMed  Google Scholar 

  99. 99.

    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–44.

    Google Scholar 

  100. 100.

    Hendy A, Hernandez-Acosta E, Chaves BA, Fé NF, Valério D, Mendonça C, et al. Into the woods: Changes in mosquito community composition and presence of key vectors at increasing distances from the urban edge in urban forest parks in Manaus, Brazil. Acta Trop. 2020;206:105441.

    CAS  PubMed  PubMed Central  Google Scholar 

  101. 101.

    Davis NC, Shannon RC. Studies on yellow fever in South America : V. Transmission experiments with certain species of Culex and Aedes. J Exp Med. 1929;50:803–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  102. 102.

    Davis NC, Shannon RC. Further attempts to transmit yellow fever with mosquitoes of South America. Am J Epidemiol. 1931;14:715–22.

    Google Scholar 

  103. 103.

    Whitman L, Antunes PCA. Studies on the capacity of mosquitoes of the genus Haemagogus to transmit yellow fever 1. Am J Trop Med Hyg. 1937;17(suppl. 1):825–31.

    Google Scholar 

  104. 104.

    Shannon RC, Whitman L, Franca M. Yellow fever virus in jungle mosquitoes. Science. 1938;88:110–1.

    CAS  PubMed  Google Scholar 

  105. 105.

    Dégallier N, Travassos da Rosa APA, Vasconcelos PFC, Travassos da Rosa ES, Rodrigues SG, Sa FGC, et al. New entomological and virological data on the vectors of sylvatic yellow fever in Brazil. J Braz Assoc Advanc Sci. 1992;44:136–42.

    Google Scholar 

  106. 106.

    Vasconcelos PFC, Rodrigues SG, Degallier N, Moraes MAP, Travassos Da Rosa JFS, Travassos Da Rosa ES, et al. An epidemic of sylvatic yellow fever in the southeast region of Maranhao State, Brazil, 1993–1994: epidemiologic and entomologic findings. Am J Trop Med Hyg. 1997;57:132–7.

    CAS  PubMed  Google Scholar 

  107. 107.

    Vasconcelos PFC, Rosa APAT, Rodrigues SG, Rosa EST, Monteiro HAO, Cruz ACR, et al. Yellow fever in Pará State, Amazon region of Brazil, 1998–1999: entomologic and epidemiologic findings. Emerg Infect Dis. 2001;7:565–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  108. 108.

    Lourenço-de-Oliveira R, Vazeille M, de Filippis AMB, Failloux A-B. Large genetic differentiation and low variation in vector competence for dengue and Yellow Fever viruses of Aedes albopictus from Brazil, the United States, and the Cayman Islands. Am J Trop Med Hyg. 2003;69:105–14.

    PubMed  Google Scholar 

  109. 109.

    da Cardoso JC, 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–24.

    PubMed Central  Google Scholar 

  110. 110.

    Moreno ES, Rocco IM, Bergo ES, Brasil RA, Siciliano MM, Suzuki A, et al. Reemergência de febre amarela: Detecção de transmissão no estado de São Paulo, Brasil, 2008. Rev Soc Bras Med Trop. 2011;44:290–6.

    PubMed  Google Scholar 

  111. 111.

    Goenaga S, Fabbri C, Dueñas JCR, Gardenal CN, Rossi GC, Calderon G, et al. Isolation of yellow fever virus from mosquitoes in Misiones province, Argentina. Vector-Borne Zoonotic Dis. 2012;12:986–93.

    PubMed  Google Scholar 

  112. 112.

    Cunha MS, Faria NR, Caleiro GS, Candido DS, Hill SC, Claro IM, et al. Genomic evidence of yellow fever virus in Aedes scapularis, southeastern Brazil, 2016. Acta Trop. 2020;205:105390.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

To Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ), and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) and Instituto Serrapilheira for funding this work. To Marcelo Gomes Quintela, Marcelo Celestino dos Santos and Renato Carvalho for their help in the fieldwork. To Agostinho C. N. Pereira for the help on mosquito identification. To Instituto Estadual do Ambiente (INEA-RJ) for their valuable support during the fieldwork. To Grupo Técnico de Vigilância de Arboviroses, Ministério da Saúde for the access to epidemiological data and support for the fieldwork.

Funding

Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro (Grant no. E-26/010.001537/2014). Conselho Nacional de Desenvolvimento Científico e Tecnológico (Grant no. 430808/2018-8). Serrapilheira Institute (Grant no. 1912-32354).

Author information

Affiliations

Authors

Contributions

RLO and FVSA conceived the study. FVSA, MSASN, and MAM were involved in mosquito sampling. MSASN and MAM identified mosquitoes. PSM, MSR, CMGD, FVSA, and RLO obtained geographical coordinates. ARR, MAFMM, CSA, and FVSA organized the data bank. CSA and MAFMM performed ecological analyses. ARR, MAFMM, and CB performed geographical analyses and prepared the maps. CSA, FVSA, and RLO wrote the manuscript. All authors reviewed the paper.

Corresponding authors

Correspondence to Filipe Vieira Santos de Abreu or Ricardo Lourenço-de-Oliveira.

Ethics declarations

Ethics approval and consent to participate

Methods were approved by the Brazilian Ministry of Environment (SISBIO Nº 52472-2 and 54707-5) and Rio de Janeiro’s Environment Agency (INEA 012/2016 and 019/2018) as well as by the Institutional Research Ethics Committee (CAAE: 12437919.1.0000.5248; CEP-Instituto Oswaldo Cruz).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1

. Top-ranked candidate models explaining variation in Haemagogus MIR in the sampling points during YFV outbreak in Rio de Janeiro, RJ, Brazil. Table S2. Top-ranked candidate models explaining variation in the positivity of sampling points during YFV outbreak in Rio de Janeiro, RJ, Brazil.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Abreu, F.V.S., de Andreazzi, C.S., Neves, M.S.A.S. et al. Ecological and environmental factors affecting transmission of sylvatic yellow fever in the 2017–2019 outbreak in the Atlantic Forest, Brazil. Parasites Vectors 15, 23 (2022). https://doi.org/10.1186/s13071-021-05143-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-021-05143-0

Keywords

  • Haemagogus
  • Mosquito
  • Nonhuman primate
  • Functional traits