Beginning in late 2016 Brazil faced the worst outbreak of Yellow Fever in recent decades, mainly located in southeastern rural regions of the country. In the present study we characterize the Yellow Fever Virus (YFV) associated with this outbreak in São Paulo State, Brazil. Blood or tissues collected from 430 dead monkeys and 1030 pools containing a total of 5,518 mosquitoes were tested for YFV by quantitative RT-PCR, immunohistochemistry (IHC) and indirect immunofluorescence. A total of 67 monkeys were YFV-positive and 3 pools yielded YFV following culture in a C6/36 cell line. Analysis of five nearly full length genomes of YFV from collected samples was consistent with evidence that the virus associated with the São Paulo outbreak originated in Minas Gerais. The phylogenetic analysis also showed that strains involved in the 2016–2017 outbreak in distinct Brazilian states (i.e., Minas Gerais, Rio de Janeiro, Espirito Santo) intermingled in maximum-likelihood and Bayesian trees. Conversely, the strains detected in São Paulo formed a monophyletic cluster, suggesting that they were local-adapted. The finding of YFV by RT-PCR in five Callithrix monkeys who were all YFV-negative by histopathology or immunohistochemistry suggests that this YFV lineage circulating in Sao Paulo is associated with different outcomes in Callithrix when compared to other monkeys.

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.

Yellow Fever (YF) is a severe disease caused by Yellow Fever Virus (YFV), endemic in some parts of Africa and America. In Brazil, YFV is maintained by a sylvatic transmission cycle involving non-human primates (NHP) and forest canopy-dwelling mosquitoes, mainly Haemagogus-spp and Sabethes-spp. Beginning in 2016, Brazil faced one of the largest Yellow Fever (YF) outbreaks in recent decades, mainly in the southeastern region. In São Paulo city, YFV was detected in October 2017 in Aloutta monkeys in an Atlantic Forest area. From 542 NHP, a total of 162 NHP were YFV positive by RT-qPCR and/or immunohistochemistry, being 22 Callithrix-spp. most from urban areas. Entomological collections executed did not detect the presence of strictly sylvatic mosquitoes. Three mosquito pools were positive for YFV, 2 Haemagogus leucocelaenus, and 1 Aedes scapularis. In summary, YFV in the São Paulo urban area was detected mainly in resident marmosets, and synanthropic mosquitoes were likely involved in viral transmission.

Classical insect-flaviviruses (cISFVs) and dual host-related insect-specific flavivirus (dISFV) are within the major group of insect-specific flavivirus. Remarkably dISFV are evolutionarily related to some of the pathogenic flavivirus, such as Zika and dengue viruses. The Evolutionary relatedness of dISFV to flavivirus allowed us to investigate the evolutionary principle of host adaptation. Additionally, dISFV can be used for the development of flavivirus vaccines and to explore underlying principles of mammalian pathogenicity. Here we describe the genetic characterization of a novel putative dISFV, termed Guapiaçu virus (GUAPV). Distinct strains of GUAPV were isolated from pools of Aedes terrens and Aedes scapularis mosquitoes. Additionally, we also detected viral GUAPV RNA in a plasma sample of an individual febrile from the Amazon region (North of Brazil). Although GUAPV did not replicate in tested mammalian cells, 3′UTR secondary structures duplication and codon usage index were similar to pathogenic flavivirus.