Fast surveillance response reveals the introduction of a new Yellow Fever Virus sub-lineage in 2021, in Minas Gerais, Brazil

In Brazil, the yellow fever virus (YFV) (family Flaviviridae, genus Flavivirus) exhibited two epidemiologically distinct transmission cycles: urban and sylvatic. In the urban cycle, which has not been recorded in Brazil since 1942, the virus is transmitted among humans by the vector Aedes aegypti.^{(1, 2)} In the sylvatic/jungle cycle, the virus is transmitted by wild mosquitoes (mainly belonging to Haemagogus and Sabethes genera) to non-human primates (NHPs) (e.g., those belonging to Alouatta and Callithrix genera) and occasionally to unvaccinated humans in close contact with forest areas.⁽³⁾ Yellow fever (YF) is endemic to the tropical rainforest of the Amazon region; from this region, waves of epidemic expansion spread toward other Brazilian regions at irregular intervals of time.⁽³⁾ During these waves, the virus usually reaches the states of Goiás (GO), in the Central-West region, and Minas Gerais (MG), in the Southeast region. Occasionally, the virus is spread to other southeast states, such as São Paulo (SP) (in 2000, 2008-2009, 2017-2018),^(4-10) Rio de Janeiro (RJ) and Espírito Santo (ES) (between 2017 and 2019).^(11-14) Even less frequently, YFV can be detected in the southernmost states of the country: Rio Grande do Sul (RS) (2001, 2008-2009, 2020-2021),^(15-17) Paraná (PR) and Santa Catarina (SC) (2018-2020).⁽¹⁸⁾ The YFV spread outside endemic areas raises at least three main concerns: 1) the possibility of re-emergence of an urban cycle; 2) the increased risk of YF in humans due to heterogeneous vaccination coverage outside the Amazon region; and 3) the extinction of threatened NHP species – especially the genus Alouatta, highly susceptible to YFV.^(5,7,19-23)

During sporadic expansion waves, MG has been an important corridor through which the virus travels before spreading to other Brazilian states.⁽¹⁴⁾ In recent times, major outbreaks occurred in MG in 2000-2003, 2010 and 2015-2018 (Fig. 1), with the latter being the largest sylvatic outbreak reported in the last 80 years.^(6,8,24-28) In the 2015-2018 outbreak, two viral sub-lineages were detected, crossing MG by different paths: the YFV_MG/SP/RS sub-lineage travelled through the west and southwest of the state, before reaching SP; the YFV_MG/ES/RJ/BA sub-lineage spread from the northwest, north, northeast and east of the state, before reaching ES and the state of Bahia (BA).^{(14, 17, 28)} Although epizootics have been reported since the end of 2020 (unfortunately, when no sampling was performed), the last reported viral detection in MG occurred in the first quarter of 2018.^{(14, 29)} Nevertheless, the virus was detected in the states of Pará (PA), Tocantins (TO), GO, SP, PR, SC and RS from the end of 2020 to June 2022 (Supplementary Fig. 1),^{(18, 30, 31)} raising awareness of the possibility of reintroduction in MG.

In this work, we describe a surveillance effort encompassing the technology of smartphone applications and the concerted action of Coordenação-Geral de Vigilância Arboviroses, Secretaria de Estado de Saúde de Minas Gerais, Laboratório Central de Saude Pública, Brazilian Ministry of Health, and members of Febre Amarela BR project (Portuguese for Yellow Fever BR) that succeeded in the first YFV confirmation in NHPs in MG in 2021, followed by genome sequencing, in only 10 days. Phylogenetic analyses revealed the introduction of a new sub-lineage in the extra-Amazonian region.

MATERIALS AND METHODS

Study area – This study was conducted in MG, which has the largest number of municipalities (853) in Brazil, the second largest population in the country (21,411,923 inhabitants, 10.1% of the Brazilian population) and the fourth largest area (586,513.99 km², larger than Spain, for example) (https://www.ibge.gov.br/cidades-e-estados/mg.html). The study took place in the north of MG, which is predominantly covered by the Cerrado (a savannah-like biome) (Fig. 1). This region has well-defined dry and rainy periods. August, the month of this investigation, is the driest month of the year, during winter. December presents the highest average rainfall, during summer.

Establishing an information network to strengthen epizootic surveillance in northern MG – In 2020, a project sponsored by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) led to the formation of a network that aggregated research institutions and public health services, aimed at strengthening the surveillance of epizootics in several Brazilian regions.^{(32, 33)}

For these purposes, preliminary field expeditions to 12 municipalities in MG were conducted. In each municipality, meetings and lectures were organised with public health officials and environmental surveillance officials to educate them on the importance of epizootic surveillance as a prompt sign of YFV in circulation. The agents were trained to notify epizootics using the Sistema de Informação em Saúde Silvestre Georreferenciado (SISS-Geo) web app⁽³⁴⁾ (or georeferenced wildlife health information system) developed by Fundação Oswaldo Cruz-Fiocruz, attached to the Brazilian Ministry of Health. This platform works as a public repository for collaborative monitoring and disease surveillance of wild animals in Brazil. For monitoring, any person can send pictures of sylvatic animals and mortality (epizootic) in wildlife, as well as metadata. The app collects real-time geographic coordinates and sends information to state and federal surveillance centres. We also organised fieldwork to train teams on vector and NHP sample collection. Finally, several WhatsApp application groups were created with all trainees, constituting an information network whose main aim was to distribute educational material and exchange news about epizootics.^{(32, 33)} An organisation chart showing the steps involved in organising the network is shown in Fig. 2. The Secretaria Estadual de Saúde supported all the field expeditions.

Investigation of the epizootics in northern MG – On August 24, 2021, the information network alerted us about the occurrence of epizootics in two municipalities, Icaraí de Minas and Ubaí, located in the northern region of MG, 595 km from Belo Horizonte, the capital of the state. Both municipalities remained unaffected during the 2015-2018 outbreak. On August 25, 2021, a team comprising members of the Febre Amarela BR project, members of the Secretaria Estadual de Saúde and officials from the affected municipalities were rapidly assembled and dispatched. On August 25 and 26, 2021, the team surveyed the riparian forests of those municipalities, searching for NHP carcasses to collect biological samples. All detected epizootics were registered in the SISS-Geo app and the Sistema de Informação de Agravos de Notificação (SINAN). Samples were collected following safety protocols,⁽³⁵⁾ preserved in liquid nitrogen (-196ºC) and sent to the reference laboratory of the Sistema Único de Saúde (SUS) located at the Fundação Ezequiel Dias (FUNED) in MG, as well as to the laboratory of Febre Amarela BR project for viral diagnosis.

YFV detection – Samples (different fragments of viscera) were lysed with TRIzol™, stored in liquid nitrogen and sent to the sequencing laboratory on dry ice (-80ºC). Total RNA extraction with TRIzol was performed according to the manufacturer’s instructions. YFV RNA was detected using a previously published by reverse transcription-quantitative polymerase chain reaction (RT-qPCR) protocol.⁽³⁶⁾

Genome sequencing – Samples were subjected to cDNA synthesis using the LunaScript™ RT Super-Mix Kit (NEB) following the manufacturer’s instructions. Then, multiplex tiling PCR was performed using previously published YFV primers⁽³⁷⁾ with 40 amplification cycles (denaturation: 95 °C for 15 s and annealing/extension: 65 °C for 5 min) of PCR using Q5 high-fidelity DNA polymerase (NEB). A volume of 7.5 µL of the PCR product was used for DNA library preparation using the Rapid Barcoding Kit 96 SQK-RBK110.96 (Oxford Nanopore Technologies, Oxford, UK) following the manufacturer’s instructions. The sequencing library was loaded onto an R9.4 flow cell (Oxford Nanopore Technologies) and sequenced for 6-18 h using the MiNKOW software. The resulting Fast5 files were base called and demultiplexed using Guppy (Version 4.4.2, Oxford Nanopore Technologies). Consensus sequences were generated by de novo assembly using Genome Detective.⁽³⁸⁾

Phylogenetic analyses – To perform phylogenetic analyses, we selected all near-complete YFV sequences of South American genotype I, available at National Center for Biotechnology Information (NCBI) (n = 315), excluding sequences < 8 kb and those of vaccine and patent-related viruses. Metadata such as sample collection dates and geographic positioning were retrieved from GenBank files or genome-associated publications (manual curation). The genome named MG66-L (reported here; details below) and 315 genomes from NCBI were aligned using MAFFT v.7.480.⁽³⁹⁾ The subsequent alignment was then used to infer a maximum likelihood tree topology in IQTREE version 2,⁽⁴⁰⁾ using GTR with 1,000 bootstrap replicates. The transfer bootstrap support for splits in the topology was inferred using Booster.⁽⁴¹⁾ The new genome sequences were sent to the NCBI GenBank database under accession numbers OL519587 to OL519589.

RESULTS

On August 25, 2021 (day 1), following the alert of our information network, members of the Febre Amarela BR project found and notified an epizootic affecting six black-and-gold howler monkeys (Alouatta caraya) in Icaraí de Minas (MG), of which only one carcass was suitable for sampling (liver, spleen and kidney samples were collected and named MG66-L, MG66-S and MG66-K, respectively). On August 26 (day 2), another epizootic was detected, affecting six black-tufted marmosets (Callithrix penicillata) in Ubaí, and only one carcass met the conditions for sampling collection (from the liver, spleen, kidney and brain, named MG67-L, MG67-S, MG67-k and MG67-B, respectively). Samples were sent to the laboratory (September 2, day 9) and tested (September 3, day 10). All sampled tissues were positive for YFV by RT-qPCR (Table).

The interval between the alert of the information network and YFV diagnosis was 10 days (Fig. 3). After confirming the virus detection, health authorities were immediately notified. Control measures were intensified, and further investigations were conducted in the affected and surrounding areas. Following diagnosis, a sequencing library was prepared, and sequencing was performed on the same day. There was only an 8-hour interval between RNA extraction and generation of the near-whole genome sequence.

Three near-complete YFV genomes were generated from three different fragments of viscera (liver, spleen and kidney) from the same carcass of Alouatta caraya (named MG66). A comparison of the three YFV genomes revealed no intra-host sequence differences. Therefore, the MG66-L genome was used for the phylogenetic analysis.

Such analyses revealed that the YFV genome generated here (MG66 – OL519587|Brazil|MG|Monkey|2021-08-25) was related to the South American genotype I, as expected, and clustered with an isolate (MF370546|Brazil|PA|Monkey|2017-04-05) sampled from an Alouatta caraya identified in PA (Amazonian region) in 2017 (Fig. 4A and Supplementary Fig. 2), indicating the introduction of a new YFV sub-lineage in MG, now called YFV_PA/MG. Moreover, they were not related to the viruses detected in the same year (2021) in RS, the southernmost state of Brazil. Further analysis of the polyprotein encoded by genomes MF370546 (PA-2017) and OL519587 (MG-2021) revealed only five amino acid differences (T128I, X108I, N808D, V3060I and A3149V), whereas comparison of sequences OL519587 and MZ712133 (RS-2021) revealed nine amino acid differences. These findings confirm the expansion and circulation of at least two YFV sub-lineages (YFV_PA/MG and YFV_MG/SP/RS) in the extra-Amazon region in 2021.⁽¹⁷⁾ The nucleotide difference matrix between the selected genomes is shown in Fig. 4B.

DISCUSSION

MG has historically been affected by several YFV introductions. However, owing to limited logistics, it has been difficult to conduct surveillance and collect biological samples promptly, both of which require a robust field and laboratory surveillance infrastructure. In this study, we describe the results of setting up an inter-institutional network capable of collaborative actions to investigate epizootics, which resulted in the first confirmation of YFV circulation in MG, Brazil, in 2021.

MG has played a role as an YFV dissemination route during expansion waves from the Amazon region. Despite the limited resources available in the past, studies on viral dispersion routes in the region have been conducted since 1938 (Supplementary Fig. 3).^{(24, 42)} In the last 2 decades, outbreaks that resulted in human cases were recorded in 2001 (32 human cases), 2002-2003 (63 human cases) and two isolated cases between 2008-2009. The largest outbreak occurred between 2016-2018, with 1,006 human cases and 448 confirmed epizootic cases.^{(26, 27, 43, 44)} Importantly, there were no previous records of viral circulation in the two affected municipalities, suggesting the presence of a naive NHP population.^{(45, 46)}

Since 2019, to the best of our knowledge, despite the notification of several epizootics, no new reports on YFV circulation have been published in the north of MG.⁽²⁹⁾ The wide territorial extension, with huge rural areas far from urban and research centres, combined with socioeconomic inequality and the lack of trained professionals in several municipalities, made it difficult to implement appropriate surveillance of epizootics, thus significantly reducing the chances of timely sample collection for virus detection. These challenges reaffirm the need to introduce new tools and strategies to strengthen surveillance. Among them, the following points are worth mentioning. a) The use of smartphone applications to notify epizootics and collect the geographic coordinates of the occurrence in real time. In this sense, the SISS-Geo app, which has been progressively implemented in Brazil, speeds up the arrival of information to involved institutions and provides accurate geographic location and streamlines outbreak response.⁽³⁴⁾ b) The use of instant messaging apps increases the number of “watchers,” providing information through real-time text transmission. Creating WhatsApp groups with a local or regional reach, comprising different members of the population (health agents, cyclists, hikers, students, rural workers, etc.), increases the chances of detecting epizootics and speeding up the dissemination of important information such as vaccination campaigns, NHPs and ecosystem conservation efforts.^{(32, 33)} Ultimately, this network of professional and citizen-driven wildlife disease surveillance serves both public health and science popularisation, the latter being a crucial yet neglected part in epidemic management. c) The creation of multi-institutional networks to facilitate and speed up sample collection and diagnosis. In the present study, research institutions linked to the Febre Amarela BR project provided an integrated, multidisciplinary team, aggregating staff from universities, students and governmental organs at municipal, state and federal health secretariats whose coordinates and combined actions allowed the achievement of the results reported here.

Genomic surveillance efforts have shown that the virus circulating in MG was related to the virus circulating in PA, in the Amazon region, in 2017 (sub-lineage YFV_PA/MG). This reveals a new wave of viral expansion in MG, 6 years after the outbreak associated with the YFV_MG/SP/RS and YFV_MG/ES/RJ/BA sub-lineages, between 2015 and 2018. Interestingly, in 2021, at least two different viral sub-lineages were found circulating at the same time in the extra-Amazon region, represented by YFV_MG/SP/RS, detected in RS,⁽¹⁷⁾ and YFV_PA/MG. Notably, the latter was detected at the height of the dry season in one of the driest regions of MG, which makes it more difficult for the main species of mosquito vectors (i.e., Haemagogus janthinomys and Hg. leucocelaenus) to survive.⁽¹³⁾

In this study, the importance of coordinated efforts between local populations and different institutions is emphasised to increase surveillance capacity. Through this system, it was possible to confirm YFV occurrence, which enabled the establishment of timely response measures. Genomic surveillance has allowed the identification of a new YFV expansion wave in MG, which needs to be monitored to ensure that effective preventive measures are implemented in appropriate time and place. These include an increase in vaccine coverage in humans and improved communication regarding the risk of YFV to health professionals and the general population to avoid human cases. This study also highlights the need for awareness and public information on the role of NHPs as sentinels for the occurrence of YFV. Finally, this study highlights the importance of collaborative, integrated surveillance, allowing prompt action in the event of YF suspicion.

ACKNOWLEDGEMENTS

To the contributions of the Divisão de Vigilância Sanitária e Ambiental from MG, to the effort of the Grupo de Vigilância da Febre Amarela of MG that was at the forefront of the preparation for virus outbreak response and field investigations, to Beatrízio Rodrigues Almeida, Hermes Almeida and Joselio Ribeiro Paraíso, for their valuable help during the field work, to Aline Tátila Ferreira, for her help in creating images, to Fernando Lucas Melo, for their valuable help with genome sequencing, to countless colleagues from municipalities’ health departments, who conducted the investigation of epizootics collecting samples in the field and to the Coordenação-Geral de Vigilância Arboviroses, Brazilian Ministry of Health. The Febre Amarela BR project (https://www.febreamarelabr.com.br/) is supported by grants from Conselho Nacional de Desenvolvimento Científico e Tecnológico/Departamento de Ciência e Tecnologia/Secretaria de Ciência, Tecnologia e Insumos Estratégicos/Brazilian Ministry of Health (CNPq/Decit/SCTIE/MS) (grant number 443215/2019-7). MSA is granted a post-doctoral scholarship (DTI-A) from CNPq. ACF, FSC, BMR and PMR are CNPq research fellows. FVSA is on receipt of a grant (APQ-01403-21) from Fundação de Amparo à Pesquisa do Estado de Minas Gerais (FAPEMIG) and (21/2020) from Instituto Federal do Norte de Minas Gerais (IFNMG).

AUTHORS’ CONTRIBUTION

Conceived and designed the experiments: FVSA, MSA, FSC, AASC, JCC, ES, DST, GRA, ACF, BMR, PMR and MABA. Performed the experiments: MSA, VF, MAP, LOM, SMAT, MEGS, JCC, ES, LCB, CMDS, NFDM, CHO, AJJS, SBV, MAMG, GGM, LOM and MABA. Analysed the data: MSA, FSC, VF, AASC, FVSA, APS, JCC, ES, DST, APMR, BMR, AODT, PMR and MABA. Contributed reagents/materials/analysis tools: ACF, BMR, PMR, AODT, DCCC and MABA. Contributed to the writing of the manuscript: MSA, FSC, VF, AASC, MABA, APS, JCC, ES, AODT, DST, SBV, APMR, BMR, PMR and FVSA. All authors have read and agreed to the published version of the manuscript. We declare that no conflict of interest exists.