- Open Access
Ecological and environmental factors affecting transmission of sylvatic yellow fever in the 2017–2019 outbreak in the Atlantic Forest, Brazil
Parasites & Vectors volume 15, Article number: 23 (2022)
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.
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.
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.
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.
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 .
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 . 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 . 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 . In total, 27 and 262 human cases and 9 and 84 deaths were confirmed for RJ in 2017 and 2018, respectively . In 2019, there was only one confirmed NHP infection with YFV in RJ . 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 . 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.
Mosquito captures were made according to methods described elsewhere . 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 . 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 . 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)  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 . Functional groups were identified based on a functional dendrogram using the Ward hierarchical agglomerative clustering method . Principal components analysis (PCA) was performed using the package FactoMineR  to summarize the relationships among traits and identify those most commonly shared among species from the same functional group.
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 , 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.
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 . 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 .
Mosquito community diversity was characterized by two measures: (1) mosquito species richness and (2) the Shannon–Wiener index (Table 2) . The functional diversity of mosquito communities was calculated on a continuous scale using four different indices: (1) functional richness (FRic) , (2) functional evenness (FEve) , (3) functional divergence (FDiv) , and (4) functional dispersion (FDis) (Table 2) . The community-weighted means (CWM) of traits  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].
(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 . The dredge function of the MuMIn package  was used to create a set of models with all combinations of variables. The Akaike information criterion with correction for small sample size (AICc)  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 . 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.
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).
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).
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).
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).
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).
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).
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).
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 , and YFV has the potential to adapt to this species , which is why it can be considered a potential vector . 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 .
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 . 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 . 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 . 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 . 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 . 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 , 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 , nearly 6 km for Hg. leucocelaenus and 11 km for Hg. janthinomys/capricornii . 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 . 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 . 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 .
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.
Monath TP, Vasconcelos PFC. Yellow fever. J Clin Virol. 2015;64:160–73.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Franco O. História da Febre Amarela no Brasil. Ministério da Saúde. Rio de Janeiro: Departamento Nacional de Endemias Rurais; 1969.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Myneni RB, Hall FG, Sellers PJ, Marshak AL. The interpretation of spectral vegetation indexes. IEEE Trans Geosci Remote Sens. 1995;33:481–6.
Consoli RAGB, de Oliveira RL. Principais mosquitos de importância sanitária no Brasil. Rio de Janeiro: Editora Fiocruz; 1994.
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.
Forattini OP. Culicidologia médica: identificação, biologia, epidemiologia, vol. 2. São Paulo: Edusp; 2002.
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.
Murtagh F, Legendre P. Ward’s hierarchical agglomerative clustering method: which algorithms implement Ward’s criterion? J Classif. 2014;31:274–95.
Lê S, Josse J, Husson F. FactoMineR: an R package for multivariate analysis. J Stat Softw. 2008;25:1–18.
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.
Anderson MJ. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001;26:32–46.
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.
Magurran A. Measuring biologcial diversity. Oxford: Wiley-Blackwell; 2004.
Villéger S, Mason NWH, Mouillot D. New multidimensional functional diversity indices for a multifaceted framework in functional ecology. Ecology. 2008;89:2290–301.
Laliberté E, Legendre P. A distance-based framework for measuring functional diversity from multiple traits. Ecology. 2010;91:299–305.
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.
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.
Burnham K, Anderson D. Model selection and multimodel inference—a practical information-theoretic approach. 2nd ed. New York: Springer; 2002.
Grueber CE, Nakagawa S, Laws RJ, Jamieson IG. Multimodel inference in ecology and evolution: challenges and solutions. J Evol Biol. 2011;24:699–711.
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.
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.
Hurvich CM, Tsai C-L. Regression and time series model selection in small samples. Biometrika. 1989;76:297–307.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Schmidt KA, Ostfeld RS. Biodiversity and the dilution effect in disease ecology. Ecology. 2001;82:609–19.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Martins MM. Density of primates in four semi-deciduous forest fragments of São Paulo, Brazil. Biodivers Conserv. 2005;14:2321–9.
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.
Behie AM, Pavelka MS. Primates in fragments: complexity and resilience. New York: Springer; 2013.
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.
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.
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.
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.
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.
Day J. Mosquito oviposition behavior and vector control. Insects. 2016;7:65.
Arnell JH. Mosquito studies (Diptera, Culicidae) XXXII. A revision of the genus Haemagogus. Contrib Am Entomol Inst. 1973;10:1–174.
Lounibos LP, Machado-Allison CE. Oviposition and egg brooding by the mosquito Trichoprosopon digitatum in cacao husks. Ecol Entomol. 1983;8:475–8.
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.
Shannon RC. The environment and behavior of some Brazilian mosquitoes. Proc Ent Soc Wash. 1931;33:1–27.
Bates M. The natural history of mosquitoes. New York: The Macmillan Company; 1949.
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.
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.
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.
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.
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.
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.
Nielsen HT. Swarming and some other habits of Mansonia perturbans and Psorophora ferox (Diptera: Culicidae). Behaviour. 1964;24:67–88.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Chadee DD. Seasonal abundance and diel landing periodicity of Sabethes chloropterus (Diptera: Culicidae) in Trinidad, West Indies. J Med Entomol. 1990;27:1041–4.
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.
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.
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.
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.
Davis NC, Shannon RC. Further attempts to transmit yellow fever with mosquitoes of South America. Am J Epidemiol. 1931;14:715–22.
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.
Shannon RC, Whitman L, Franca M. Yellow fever virus in jungle mosquitoes. Science. 1938;88:110–1.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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).
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
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
About this article
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
- Nonhuman primate
- Functional traits