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.

Avian malaria parasites are widespread parasites transmitted by Culicidae insects belonging to different genera. Even though several studies have been conducted recently, there is still a lack of information about potential vectors of Plasmodium parasites, especially in Neotropical regions. Former studies with free-living and captive animals in São Paulo Zoo showed the presence of several Plasmodium and Haemoproteus species. In 2015, a pilot study was conducted at the zoo to collect mosquitoes in order to find out (i) which species of Culicidae are present in the study area, (ii) what are their blood meal sources, and (iii) to which Plasmodium species might they be potential vectors. Mosquitoes were morphologically and molecularly identified. Blood meal source and haemosporidian DNA were identified using molecular protocols. A total of 25 Culicidae species were identified, and 6 of them were positive for Plasmodium/Haemoproteus DNA. Ten mosquito species had their source of blood meal identified, which were mainly birds, including some species that were positive for haemosporidian parasites in the former study mentioned. This study allowed us to expand the list of potential vectors of avian malaria parasites and to improve our knowledge of the evolutionary and ecological relationships between the highly diverse communities of birds, parasites, and vectors present at São Paulo Zoo.

Culex is the largest subgenus within the genus Culex that includes important vectors of diseases. The correct identification of mosquitoes is critical for effective control strategies. Wing geometric morphometrics (WGM) has been used to identify mosquito species alongside traditional identification methods. Here, WGM was used for eleven Culex species from São Paulo, Brazil, and one from Esquel, Argentina. Adult mosquitoes were collected using CDC (Centers for Disease Control) traps, morphologically identified and analyzed by WGM. The canonical variate analysis (CVA) was performed and a Neighbor-joining (NJ) tree was constructed to illustrate the patterns of species segregation. A cross-validated reclassification test was also carried out. From 110 comparisons in the cross-validated reclassification test, 87 yielded values higher than 70%, with 13 comparisons yielding 100% reclassification scores. Culexquinquefasciatus yielded the highest reclassification scores among the analyzed species, corroborating with the results obtained by the CVA, in which Cx. quinquefasciatus was the most distinct species. The high values obtained at the cross-validated reclassification test and in the NJ analysis as well as the segregation observed at the CVA made it possible to distinguish among Culex species with high degrees of confidence, suggesting that WGM is a reliable tool to identify Culex species of the subgenus Culex.

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.

São Paulo (SP), a densely populated state in southeast Brazil that contains one of the world’s largest urban regions, has experienced its largest yellow fever virus (YFV) outbreak in decades. Surveillance in non-human primates (NHP) is important in order to detect YFV early during an epidemic or epizootic, to quantify the magnitude of the outbreak in NHP, and to evaluate the risk of YFV spillover infection in human populations. To better understand the genetic diversity and spatial distribution of YFV during the current outbreak in southeast Brazil, we generated 46 new virus genomes from YFV positive cases identified in 18 different municipalities in SP, mostly sampled from non-human primates between April 2017 and February 2018. Our data show that most NHP cases in São Paulo state were likely caused by the introduction of a single YFV lineage from Minas Gerais to São Paulo. Phylogenetic and phylogeographic analyses of these data indicate that YFV spread southwards from Minas Gerais into São Paulo state at a typical rate of <1km per day. These results shed light in the sylvatic transmission of yellow fever in highly fragmented forested regions in São Paulo state and highlight the importance of continued operational research and surveillance of zoonotic pathogens in sentinel populations. Author’s Summary 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. YFV is not normally present in São Paulo state and therefore a large proportion of the 18 million inhabitants of the state have not been vaccinated against YFV. The presence of YFV in São Paulo state therefore represents a serious threat to public health. In Brazil, YFV typically circulates among non-human primates, with cases in humans representing isolated ‘spillover’ events from this predominantly sylvatic cycle. Understanding the epidemiological dynamics and spread of YFV in primates is therefore critical for contextualising human cases, and guiding vaccination strategies that can better protect local human populations. Here, we analyse the geographic and temporal distribution of observed cases of YFV in non-human primates in São Paulo state. We generate sequence data from 46 YFV positive cases, and perform phylogenetic and phylogeographic analyses aimed at understanding the spatial spread of YFV in São Paulo state. We show that most cases in non-human primates in the São Paulo state were likely caused by a single introduction of YFV from Minas Gerais to São Paulo. Analyses of these data indicate that YFV spread southwards from Minas Gerais into São Paulo state at a typical rate of <1 km per day, consistent with a scenario of continued spread in non-human primates and sylvatic vectors across forested patches, with occasional spillover to unvaccinated human populations. Article Summary Line Genomic surveillance of yellow fever in São Paulo during the yellow fever 2017-2018 epidemic reveals movement towards Atlantic coast.