Abstract
Despite the considerable morbidity and mortality of Yellow fever virus (YFV) infections in Brazil, as well as its widespread presence in non-human primate host, our understanding of disease outbreaks is hampered by limited viral genomic data. Determining the timing and spatial corridors of YFV spread, as well as the geographic hotspots that link the endemic north of the country with epidemic extra-Amazonian regions, are central to predicting and preventing future outbreak events and epidemics. Here, we tracked the recent spread of the virus by integrating genome sequences of new YFV infections sampled from infected non-human primates and humans with both epidemiological and vector data. Through a combination of phylogenetic and epidemiological models we reconstructed the recent transmission history of YFV within different epidemic seasons in Brazil. A suitability index based on the highly domesticated Aedes aegypti was able to capture the seasonality of reported human infections. Spatial modelling revealed spatial hotspots with both past reporting and low vaccination coverage, which coincided with many of the largest urban centres in the Southeast. Phylodynamic analysis unravelled the circulation of three distinct YFV lineages, and provided proof of the directionality of a known spatial corridor of viral spread that connects the endemic North with the extra-Amazonian basin. This study illustrates that genomics linked with field sampling of animals and humans within a One Health framework can provide new insights into the landscape of YFV transmission, augmenting traditional approaches to infectious disease surveillance and control.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
This work was supported in part through National Institutes of Health USA grant U01 AI151698 for the United World Arbovirus Research Network (UWARN) and the Brazilian Ministry of Health (SCON2021-00180). MG is funded by PON Ricerca e Innovazione 2014-2020. JL is funded by BioISI (Biosystems and Integrative Sciences Institute), Faculdade de Ciencias da Universidade de Lisboa (UIDP/4046/2020). FP is funded by the One Health Poultry Hub, a Global Challenges Research Fund (GCRF) and United Kingdom Research and Innovation (UKRI) initiative. ECH is funded by an Australian Research Council Australian Laureate Fellowship (FL170100022). CNDS was founded by CNPq (307176-2018-5). Authors additionally thanks Fiocruz/CVSLR for the logistic support.
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
The project was supported by the Pan American World Health Organization (PAHO) and the Brazilian Ministry of Health (MoH) as part of the arboviral genomic surveillance efforts within the terms of Resolution 510/2016 of CONEP (Comissao Nacional de Etica em Pesquisa, Ministerio da Saude; National Ethical Committee for Research, Ministry of Health).
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.
Yes
Data Availability
New sequences generated as part of this study have been deposited in GenBank under accession numbers OP508570-OP508716 (Table 1). Aggregate estimates of YFV temporal suitability (index P) and its input (humidity, temperature, rainfall) for the South and Southeast between 2014 and 2021 are provided as a comma separated values (CSV) file (Supplementary Data File 1). Estimates of the probability of reporting human cases per municipality for the South and Southeast (output of spatial regression model) are provided as a comma separated values (CSV) file (Supplementary Data File 2) geocoded with GEOCMU7 (a seven digit unique identifier).
