São Paulo, a densely inhabited state in southeast Brazil that contains the fourth most populated city in the world, recently experienced its largest yellow fever virus (YFV) outbreak in decades. YFV does not normally circulate extensively in São Paulo, so most people were unvaccinated when the outbreak began. Surveillance in non-human primates (NHPs) is important for determining the magnitude and geographic extent of an epizootic, thereby helping to evaluate the risk of YFV spillover to humans. Data from infected NHPs can give more accurate insights into YFV spread than when using data from human cases alone. To contextualise human cases, identify epizootic foci and uncover the rate and direction of YFV spread in São Paulo, we generated and analysed virus genomic data and epizootic case data from NHPs in São Paulo. We report the occurrence of three spatiotemporally distinct phases of the outbreak in São Paulo prior to February 2018. We generated 51 new virus genomes from YFV positive cases identified in 23 different municipalities in São Paulo, mostly sampled from NHPs between October 2016 and January 2018. Although we observe substantial heterogeneity in lineage dispersal velocities between phylogenetic branches, continuous phylogeographic analyses of generated YFV genomes suggest that YFV lineages spread in São Paulo at a mean rate of approximately 1km per day during all phases of the outbreak. Viral lineages from the first epizootic phase in northern São Paulo subsequently dispersed towards the south of the state to cause the second and third epizootic phases there. This alters our understanding of how YFV was introduced into the densely populated south of São Paulo state. Our results shed light on the sylvatic transmission of YFV in highly fragmented forested regions in São Paulo state and highlight the importance of continued surveillance of zoonotic pathogens in sentinel species. Since July 2016, the southeast region of Brazil has experienced the largest yellow fever virus (YFV) outbreak in decades. São Paulo is the most densely populated state in southeast Brazil. The outbreak has caused serious public health concern in the state, as YFV does not normally circulate widely there and most of the 21 million inhabitants were correspondingly unvaccinated against YFV when the outbreak began. In Brazil, YFV typically circulates among non-human primates (NHPs), and human cases represent isolated spill over events from this predominantly sylvatic cycle. Understanding the epidemiological dynamics and spread of YFV in NHPs is therefore critical for contextualising human cases, and guiding vaccination strategies that can better protect local human populations. Here, we aim to contextualise human cases, identify epizootic foci and uncover the rate and direction of YFV spread in São Paulo. We analyse the geographic and temporal distribution of observed cases of YFV in NHPs in São Paulo state, and identify three distinct phases of the epizootic. We generate virus genome sequences from 51 YFV-positive cases and perform phylogenetic and phylogeographic analyses aimed at understanding the spatial spread of YFV in São Paulo state. Analyses of these data indicate that YFV spread from the north of São Paulo state into more densely populated southern regions. Although we observe substantial heterogeneity in the rate at which different sampled YFV lineages spread, the typical rate of spread was low with a mean rate of ~1 km per day. This is consistent with a scenario in which the majority of transmission events occurred between NHPs primates and sylvatic vectors across forested patches.

After detecting the death of Howlers monkeys (genus Alouatta) and isolation of yellow fever virus (YFV) in Buri county, São Paulo, Brazil, an entomological research study in the field was started. A YFV strain was isolated from newborn Swiss mice and cultured cells of Aedes albopictus - C6/36, from a pool of six Haemagogus (Conopostegus) leucocelaenus (Hg. leucocelaenus) mosquitoes (Dyar & Shannon) collected at the study site. Virus RNA fragment was amplified by RT-PCR and sequenced. The MCC Tree generated showed that the isolated strain is related to the South American I genotype, in a monophyletic clade containing isolates from recent 2008-2010 epidemics and epizootics in Brazil. Statistical analysis commonly used were calculated to characterize the sample in relation to diversity and dominance and indicated a pattern of dominance of one or a few species. Hg. leucocelaenus was found infected in Rio Grande do Sul State as well. In São Paulo State, this is the first detection of YFV in Hg. leucocelaenus.

In 2000, an outbreak of sylvatic yellow fever possibly occurred in gallery forests of the Grande river in the Paraná basin in the northwestern region of São Paulo state. The aim of this study was to obtain information on the bionomics of Haemagogus and other mosquitoes inside tree holes in that area. Eighteen open tree holes were sampled for immature specimens. Adults were collected twice a month in the forest in Santa Albertina county from July 2000 to June 2001. The seasonal frequency of fourth instars was obtained by the Williams geometric mean (Mw), while the adult frequency was estimated either by hourly arithmetic or the Williams' means. Cole's index was applied to evaluate larval inter-specific associations. Among the ten mosquito species identified, the most abundant was Aedes terrens Walker followed by Sabethes tridentatus Cerqueira and Haemagogus janthinomys Dyar. Larval and adult abundance of these species was higher in summer than in winter. Although larval abundance of Hg. janthinomys peaked in the rainy season, correlation with rainfall was not significant. Six groups of larval associations were distinguished, one of which the most positively stable. The Hg. janthinomys and Ae. terrens association was significant, and Limatus durhamii Theobald was the species with most negative associations.

Yellow fever virus (YFV) is an endemic arbovirus in parts of Africa and the Americas. In Brazil, following the eradication of the urban transmission cycle, YFV is maintained in a sylvatic cycle involving several species of neotropical primates and mosquitoes of the genera Haemagogus and Sabethes, which serve as primary and secondary vectors, respectively. During the 2016–2019 outbreak in São Paulo State, a total of 3731 mosquito pools were collected from sites with ongoing epizootic events in 192 municipalities. The RT-qPCR analysis detected YFV in 46 pools (1.4%) across nine mosquito species, including both primary and secondary vectors, as well as species from the genera Aedes and Psorophora. Differences in viral loads were observed among species. While Aedes aegypti was not found to be positive, the detection of natural YFV infection in other Aedes species raises concerns about potential virus reurbanization. Further studies are needed to clarify the role of additional mosquito species in YFV transmission in Brazil.