http://orcid.org/0000-0001-9716-6899
Figures
Abstract
Background
From the end of 2016 until the beginning of 2019, Brazil faced a massive sylvatic yellow fever (YF) outbreak. The 2016–2019 YF epidemics affected densely populated areas, especially the Southeast region, causing thousands of deaths of humans and non-human primates (NHP).
Methodology/Principal findings
We conducted a molecular investigation of yellow fever virus (YFV) RNA in 781 NHP carcasses collected in the urban, urban-rural interface, and rural areas of Minas Gerais state, from January 2017 to December 2018. Samples were analyzed according to the period of sampling, NHP genera, sampling areas, and sampling areas/NHP genera to compare the proportions of YFV-positive carcasses and the estimated YFV genomic loads. YFV infection was confirmed in 38.1% of NHP carcasses (including specimens of the genera Alouatta, Callicebus, Callithrix, and Sapajus), from the urban, urban-rural interface, and rural areas. YFV RNA detection was positively associated with epidemic periods (especially from December to March) and the rural environment. Higher median viral genomic loads (one million times) were estimated in carcasses collected in rural areas compared to urban ones.
Conclusions/Significance
The results showed the wide occurrence of YF in Minas Gerais in epidemic and non-epidemic periods. According to the sylvatic pattern of YF, a gradient of viral dissemination from rural towards urban areas was observed. A high YF positivity was observed for NHP carcasses collected in urban areas with a widespread occurrence in 67 municipalities of Minas Gerais, including large urban centers. Although there was no documented case of urban/Aedes YFV transmission to humans in Brazil during the 2016–2019 outbreaks, YFV-infected NHP in urban areas with high infestation by Aedes aegypti poses risks for YFV urban/Aedes transmission and urbanization.
Author summary
Brazil faced the most massive sylvatic yellow fever (YF) outbreak in 2016–2019. The outbreak affected highly densely populated areas, and Minas Gerais was the most affected state with thousands of deaths of human and non-human primates (NHP). We investigated the yellow fever virus (YFV) RNA in NHP carcasses collected throughout Minas Gerais in 2017 and 2018. We demonstrated the wide occurrence of YFV-infected NHP, including the viral persistence during the non-epidemic dry season of 2017. YFV RNA was detected in NHP carcasses in the urban, urban-rural interface and rural areas. We have also detected new YF cases in 49 municipalities where YF cases have not been previously detected during the outbreaks. Estimates of YFV genomic load in naturally infected NHP carcasses showed high and similar loads in specimens (Alouatta, Callithrix, and Callicebus) collected in rural areas and lower genomic loads in the urban-rural interface and urban Callithrix specimens. The presence of YFV inside urban areas poses an imminent risk, although no human case was epidemiologically linked to urban/Aedes transmission during the last outbreaks in Brazil.
Citation: Sacchetto L, Silva NIO, Rezende IMd, Arruda MS, Costa TA, de Mello ÉM, et al. (2020) Neighbor danger: Yellow fever virus epizootics in urban and urban-rural transition areas of Minas Gerais state, during 2017-2018 yellow fever outbreaks in Brazil. PLoS Negl Trop Dis 14(10): e0008658. https://doi.org/10.1371/journal.pntd.0008658
Editor: Marilia Sá Carvalho, Fundacao Oswaldo Cruz, BRAZIL
Received: January 23, 2020; Accepted: July 28, 2020; Published: October 5, 2020
Copyright: © 2020 Sacchetto et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the manuscript and its Supporting Information files.
Funding: This work was supported by different agencies. BPD was supported by Fundação de Amparo à Pesquisa do Estado de Minas Gerais (FAPEMIG), under Grant APQ-01574-17, and Secretaria de Estado de Saúde de Minas Gerais (SES-MG)/Secretaria de Estado de Planejamento de Minas Gerais (SEPLAG-MG)/Instituto René Rachou/Fundação Oswaldo Cruz (IRR/FIOCRUZ) under Grant Yellow fever. GST was supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) under grant 440593/2016-6 and Decit/SCTIE/MoH under grant 440593/2016-6. APP was supported by FAPEMIG and Coodernação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES). BdT was supported by the RESERVOIRS program supported by European (ERDF/FEDER) funds and assistance from Collectivité Territoriale de la Guyane and Direction Régionale pour la Recherche et la Technologie, and the MicroBIOME project granted by the Laboratoire d’Excellence CEBA “Investissement d’Avenir” under Grant ANR-10-LABX-2501. PAA was supported by FAPEMIG under grant APQ-01989-18. MLN was supported by Fundação de Amparo à Pesquisa do estado de São Paulo (FAPESP) under grant 13/21719-3. LS and IMR were supported by Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)/Brazil under graduate scholarship grant (0001) and NIOS was supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) under graduate scholarship grant, VEM was supported by Pró-Reitoria de Graduação-Universidade Federal de Minas Gerais (UFMG) under graduate scholarship grant. MAS, TAC, and GFGO were supported by CNPq FAPEMIG, and CAPES undergraduate scholarship grants. BPD, APP, GST, EGK, and MLN are CNPq Research Fellows. CNPq: http://www.cnpq.br/ CAPES: https://www.capes.gov.br/ FAPEMIG: https://fapemig.br/pt/ FAPESP: http://www.fapesp.br/ Pró-Reitoria de Graduação-Universidade Federal de Minas Gerais: https://www2.ufmg.br/prograd/ The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Yellow fever virus (YFV) (family Flaviviridae, genus Flavivirus) is endemic in tropical and subtropical regions of Africa, Central and South America. Yellow fever (YF) presents a broad spectrum of severity, with clinical manifestations ranging from self-limited febrile to fatal disease in humans [1,2]. There is no specific treatment for YF, and vaccination is recommended for people who live in or travel to at-risk areas [2].
Urban YF, transmitted by Aedes aegypti vector, was a severe threat to Brazil’s human health until 1942 when the last case was registered in Acre state [1]. Nowadays, the virus is maintained in the sylvatic cycle involving Neotropical non-human primates (NHP) and sylvatic mosquitoes (Haemagogus spp. and Sabethes spp.) [2]. Although these mosquitoes mainly feed on NHP, they can incidentally transmit the virus to humans [2]. The reemergence of sylvatic YF outside the Brazilian Amazon basin has been reported, especially during the rainy season, called the epidemic period [3].
Since the eradication of urban YF, Brazil faced the largest sylvatic YF outbreak from the end of 2016 until the beginning of 2019 [4]. YFV was introduced into Southeast [4–8], from the Midwest region [8] spreading southward, causing thousands of NHP and human cases and deaths [4,9–11]. From July 2016 until June 2019, 2,240 human cases and 760 deaths were confirmed, mostly in the Southeast region (99.2%), and over 2,590 YF epizootics were confirmed in Brazil [4,10–13]. Minas Gerais, Southeast Brazil, was one of the most affected states during the YF outbreaks. Between July 2016 and June 2018, the Secretary of Health of Minas Gerais confirmed 1,002 human cases and 340 deaths, and more than 490 epizootics in NHP [14]. Minas Gerais is in the recommended area for YF vaccination since 2008, but until the end of 2016, the state presented low average vaccination coverage (57.26%), reaching 92.71% by the end of 2018 [15].
The rapid spread and the extent of the latest recent YF outbreaks in the Southeast region, including highly densely populated urban centers in the states of Minas Gerais, São Paulo, Espírito Santo, and Rio de Janeiro [10,11], have raised concerns regarding YFV reintroduction in epidemiologically receptive areas, with millions of naïve people for YFV. NHP is the primary host of the YF sylvatic cycle, and even though NHP are mostly observed in rural/sylvatic zones, these animals are commonly found in urban-rural transition and urban areas in Brazil. Given the high numbers of epizootics during YF outbreaks, this study aimed at the molecular investigation of YFV in NHP from Minas Gerais. Our results showed the widespread occurrence of YFV in Minas Gerais and, hastily, YFV in NHP carcasses collected inside urban matrixes and urban-rural transition areas, including large urban centers.
Material and methods
Study area and samples
During YF outbreaks in Brazil (January 2017 to December 2018), we analyzed 771 liver samples from carcasses of free-living NHP from different areas of Minas Gerais state. Data regarding ten NHP carcasses previously screened for the YFV RNA [7] were included in data analysis, totalizing 781 individuals analyzed here. The NHP carcasses consisted of a convenient sampling, collected in the context of the Brazilian Yellow Fever Surveillance Program [16], of epizootics that occurred in 2017 and 2018. Since biological samples were obtained from carcasses, blood sampling was not possible, and liver samples were selected for YFV investigation, due to the hepatotropic of YFV in humans and NHP [1,4, 16, 17]. Minas Gerais is in Southeast Brazil (Fig 1), with 586,521.121 km2 and 21,040,662 inhabitants. The state has 853 municipalities grouped into 12 mesoregions. Three main biomes cover the state: Cerrado, Atlantic Forest, and Caatinga [18] (Fig 1).
(A) Geopolitical map of Brazil with biomes distribution. (B) Map of Minas Gerais state indicating municipalities from where NHP carcasses were collected and investigated for YFV RNA. Municipalities are colored as follows: in grey are 120 municipalities without detection of YFV in NHP (in this study); in light red are 110 municipalities with detection of YFV in NHP (this study and official bulletins); and in blue are 49 municipalities with detection of YFV in NHP (only in this study) where YFV circulation has not been described, during 2016–2018 outbreaks. Based on official bulletins, YFV was detected (humans or NHP) in another 31 municipalities of Minas Gerais by the State Surveillance Secretary during 2016–2018 outbreaks [13,21,22]. Abbreviations referring to Federal District and states as follows: AC: Acre; AL: Alagoas; AP: Amapá; AM: Amazonas; BA: Bahia; CE: Ceará; DF: Distrito Federal; ES: Espírito Santo; GO: Goiás; MA: Maranhão; MT: Mato Grosso; MS: Mato Grosso do Sul; MG: Minas Gerais; PA: Pará; PB: Paraíba; PR: Paraná; PE: Pernambuco; PI: Piauí; RR: Roraima; RO: Rondônia; RJ: Rio de Janeiro; RN: Rio Grande do Norte; RS: Rio Grande do Sul; SC: Santa Catarina; SP: São Paulo; SE: Sergipe; and TO: Tocantins. Numbers indicate the mesoregions of Minas Gerais as follows: 1. Triângulo Mineiro e Alto Paranaíba; 2. Northwest; 3. North; 4. Jequitinhonha; 5. Vale do Mucuri; 6. Vale do Rio Doce; 7. Zona da Mata; 8. South/Southwest; 9. West; 10. Central; 11. Metropolitan region; 12. Campo das Vertentes. Maps were created using the QGIS software version 3.8.2.
The NHP carcasses were collected in the field by health surveillance agents (Brazilian Yellow Fever Surveillance Program [16], frozen (-20 °C), and transported to the Laboratory of Zoonosis of Belo Horizonte/Minas Gerais (LZOON-BH). At LZOON-BH, NHP carcasses were defrosted and visually inspected regarding the stage of decomposition, taking into account different tissue aspects: (i) brightness (bright or opaque); (ii) color (regular aspect, greenish, greyish or whitish); (iii) consistency (if the tissue was firm or softened); (iv) smell (whether the carcass was smelling rotten or not); and (v) presence or absence of parasites (larvae). Based on that, liver samples were categorized according to their preservation status: good (n = 465; 59.5%), intermediate (n = 211; 27.0%), or bad (n = 59; 7.5%). After defrosting, liver fragments (0.5 cm) were collected from the carcasses and placed into microtubes containing RNA later solution (Ambion, USA) by the LZOON team. The microtubes were kept at 4 °C for 16 hours and then stored at -20 °C until used for RNA extraction at Laboratório de Vírus/Universidade Federal de Minas Gerais (UFMG).
Using the information on the address or geographic coordinates where carcasses were collected, each carcass was classified based on land use: urban (built-up) areas, urban-rural transition areas, and rural/sylvatic areas. Urban areas (built-up) corresponded to cities, towns, and villages or isolated urban areas, characterized by transformations resulting from urban development (e.g., paved streets, electric lighting, sewerage, buildings, and intense human occupation) [19]. A buffer zone (2.0 kilometers) was established on the outskirts of urban areas, delineating an urban-rural transition area (interface). The rural/sylvatic areas covered the entire area located outside these limits. Each localization of carcasses was checked using satellite imagery (Google Maps, available at https://www.google.com/maps) to determine the boundaries of the urban areas [19] urban-rural interface areas and rural/sylvatic areas as stated above. NHP carcasses analyzed here were collected in Cerrado and Atlantic Forest areas, but not in Caatinga regions.
NHP carcasses were identified at species or genus level, by the Veterinary team at LZOON or by the Zoology/Mastozoology teams at UFMG, using morphological criteria [20–22]. Out of 781 NHP analyzed in this study, 360 carcasses (collected in 2017) were forward to UFMG laboratories, and some of them were taxidermized and deposited in the Center of Taxonomic Collections of Universidade Federal de Minas Gerais (CCT-UFMG), Brazil [20–22]. The Committee on Ethics in the Use of Animals of UFMG approved the study (CEUA 98/2017).
YFV molecular investigation
Total RNA was extracted from approximately 30 mg of the liver preserved in RNA later solution (Ambion, USA), using RNeasy Minikit (Qiagen, USA). RNA extraction was performed in batches of 11 samples plus one negative control (nuclease-free water). RNA samples were submitted to a one-step real-time polymerase chain reaction (RT-qPCR), with primers targeting the gene coding for β-actin [23], using GoTaq 1-Step RT-qPCR System (Promega, USA). All samples were positive, showing the suitability of samples for RT-qPCR. RNA samples and the negative extraction controls were screened in duplicate for YFV RNA using RT-qPCR [24]. Non-template (nuclease-free water) and YFV 17DD RNA (provided by Bio-Manguinhos/FIOCRUZ) were used as negative and positive controls, respectively. Samples were considered positive when they presented amplification in duplicate, considering the threshold for cycle quantification value (Cq) of 37. Since between the range of 37–40 Cq indicate minimal quantities of DNA, with low confidence for the results, Cqs>37 or undetermined were considered negative. The results of RT-qPCR runs were manually inspected for the correction of baseline and threshold parameters whenever necessary due to heterogeneity in the amount of input RNA among different samples.
The Cqs are inversely proportional to the input target nucleic acid in the sample, which could indirectly reflect the YFV genomic load in the samples. In that way, we performed a quantitative RT-qPCR (Bio Gene Research Febre Amarela PCR kit—Bioclin-Quibasa, Brazil) using RNA obtained from 28 YFV-positive samples to have an estimative of YFV genomic load in the liver. We determined the number of genomic copies per gram (gc/g) of the liver, and we estimated the genomic loads in all YFV-positive carcasses using linear regression and the observed Cqs.
From 11 YFV-positive samples, a partial sequence of the envelope gene (primers 5-left and 6-right) was obtained [5], using the GoTaq Hot Start Colorless Master Mix (Promega, USA). The amplicons were purified and sequenced with the ABI3130 platform (Applied Biosystems, USA). The consensus sequences were generated (Geneious v.9.1.6—Biomatters, NZL; https://www.geneious.com/) and used for phylogenetic reconstruction, using the maximum likelihood method implemented in MEGA7 (https://www.megasoftware.net/). The final tree was edited and visualized in FigTree v.1.4.4 (http://tree.bio.ed.ac.uk/software/figtree/).
Statistical analyses
For statistical analyses, samples were subdivided into different groups according to the period of sampling: (i) epidemic (December-May) and non-epidemic (June-November) periods; (ii) bimester within epidemic (December/January, February/March, and April/May) and within non-epidemic periods (June/July, August/September, and October/November); (iii) NHP genera: Callithrix, Callicebus, and Alouatta (Sapajus was not included due to the low number of specimens); (iv) environment: rural, urban-rural interface, and rural areas; and (v) environment together with NHP genera: urban Callithrix, urban-rural Callithrix, rural Callithrix, rural Callicebus, and rural Alouatta (given the low number of specimens, Alouatta and Callicebus from urban and urban-rural interface were not included in the analyses). The comparison of proportions of YFV-positive carcasses within groups was performed using chi-square and chi-square for trend. The Kruskal-Wallis test and multiple pairwise-comparisons between groups by Wilcoxon rank-sum test were used to compare the (median) Cq values observed for the positive samples within each group. Median Cq values were compared only using data from carcasses in good and intermediate conservation status (see Results). Statistical analyses were run in R v.3.6.0 (https://www.r-project.org/) and Epitools (http://epitools.ausvet.com.au/). Bonferroni correction for multiple comparisons was applied, with differences considered statistically significant at p< 0.0125.
Results
We analyzed data from 781 NHP carcasses collected from January 2017 to December 2018 (2017: 50.8%, and 2018: 49.1%), during the epidemic (75.8%) and non-epidemic periods (24.2%), in urban (57.7%), urban-rural interface (6.9%), and rural areas (35.3%) (S1 Fig). The identified NHP carcasses belonged to different families and genera: Atelidae (Alouatta guariba, Alouatta caraya, and Alouatta sp.) (6.6%); Pitheciidae (Callicebus personatus, Callicebus nigrifrons, and Callicebus sp.) (4.1%); Cebidae (Sapajus sp. and Sapajus nigritus) (0.2%); and Callitrichidae (Callithrix geoffroyi, Callithrix aurita, Callithrix penicillata, Callithrix sp., and hybrid marmosets) (84.3%) (Table 1).
YFV infection was confirmed in 298 out of 781 carcasses (38.1%, 95% CI = 34.7–41.6%). Regardless the preservation status of the carcasses, all samples were analyzed, and YFV RNA was detected in the liver of carcasses in bad (16/59, 27.1%, 95% CI = 17.4–39.6%), intermediate (85/211, 40.2%, 95% CI = 33.9–47.0%) or good (187/465, 40.2%, 95% CI = 35.6–44.7%) preservation status. No association between the preservation status of carcasses and YFV RNA detection was observed (x2 = 3.91, p = 0.14).
YFV RNA was detected in specimens of Alouatta, Callicebus, Callithrix (Table 1), and Sapajus. The investigation of the association between the positivity and NHP genera considered only data regarding Callithrix, Callicebus, and Alouatta (n = 743). An association between the YFV positivity and NHP genera (x2 = 60.51, p <0.0001) was observed, with a stronger positive association between the detection of YFV RNA and specimens of Alouatta (78.8%) and Callicebus (68.7%) compared to Callithrix (31.8%) (Fig 2A, Table 2). No difference in positivity for YFV was observed between Alouatta and Callicebus (x2 = 0.60, p = 0.44), and a higher odds ratio (OR) for YFV detection was observed for both genera (OR = 7.96 and OR = 4.70, respectively), compared to Callithrix (Table 2).
YFV ND: yellow fever virus RNA was not detected. YFV pos: yellow fever virus RNA was detected. Pearson residuals (standardized) were extracted from the chi-square function and plotted. In each cell, the size of the circle is proportional to the amount of cell contribution, and colors indicate positive residuals (blue) or negative residuals (red). The vertical bars indicate Pearson residuals values. Analyses were run in R software v.3.6.0.
The detection of YFV RNA was negatively associated with the non-epidemic periods while positively associated with epidemic periods (x2 = 85.03, p< 0.0001) (Fig 2B), especially the first two bimesters of the epidemic period (Dec/Jan and Feb/Mar) (S2 Fig). An extremely high OR (OR = 14.63, 95% CI = 5.80–36.91, p< 0.0001) was observed for YFV detection in December /January followed by February/March (OR = 14.63, 95% CI = 2.00–12.35, p< 0.0001). The detection of YFV RNA was negatively associated with the last bimester of the epidemic period (Apr/May) (S2 Fig). No difference was observed between the proportions of YFV-infected carcasses collected in 2017 (n = 155/397; 39.0%, 95% CI = 34.3–43.9%) and 2018 (n = 126/384; 32.8%; 95% CI = 28.3–37.6%) (x2 = 0.19, p = 0.65). Throughout 2017, YFV was detected with a higher proportion of YFV-infected carcasses during the epidemic period (x2 = 17.37, OR = 3.28, 95% CI = 1.87–5.77, p< 0.001). In 2018, YFV was detected only in the epidemic period (S1 Fig).
NHP carcasses were collected in 279 municipalities located in 12 mesoregions of Minas Gerais state, covered by Cerrado and Atlantic Forest. YFV was detected in 159 municipalities (11 mesoregions) (Fig 1B), within urban (n = 67), urban-rural interface (n = 20), and rural (n = 106) environments. In 2017 and 2018, the highest numbers of YFV-positive carcasses were from Metropolitan (42.6%), Zona da Mata (51.0%), and the South/Southwest (43.9%) regions (S1 Table).
Regarding the environment, differences were observed among the proportions of YFV-positive carcasses collected in urban (28.7%) urban-rural interface (35.8%), and rural areas (53.6%) (x2 = 44.65, p< 0.0001) (Table 2). Analyses indicated a linear trend (x2 = 44.27, p< 0.0001) and higher odds ratio for detection of YFV from urban towards urban-rural interface (OR = 1.40) and rural areas (OR = 2.85) (Table 2). The detection of YFV had a strong positive association with the rural environment and a strong negative association with the urban environment (Fig 2C).
Next, we detected an association of YFV-positive carcasses with the environment and NHP genera altogether (x2 = 60.51, p< 0.0001) (Fig 2D). The proportions of YFV-positive carcasses of rural Alouatta (76.7%), rural Callicebus (69.2%), and rural Callithrix (42.2%) were higher than urban-rural Callithrix (33.3%) (x2≥ 7.70, p≤ 0.005) and urban Callithrix (27.5%) (x2≥ 18.34, p≤ 0.0001). A difference between the proportion of urban and rural YFV-infected Callithrix was observed (x2 = 11.76, p< 0.001), but no difference was observed between rural and urban-rural Callithrix (x2 = 2.78, p = 0.95) nor between urban and urban-rural Callithrix (x2 = 1.10, p = 0.29). Lower linear trend (x2 = 60.15, p< 0.0001) and odds ratio for YFV detection were observed, from rural Alouatta (OR = 8.66), rural Callicebus (OR = 5.90) and rural Callithrix (OR = 1.92) to urban-rural Callithrix (OR = 1.27) and finally urban Callithrix (Table 2).
Since Cqs reflect the YFV genomic load in liver samples, we analyzed the median Cq values for YFV-positive samples. There was no difference in median Cqs observed in carcasses in good versus intermediate conservation status (x2 = 9.08, p< 0.010), carcasses in bad conservation status were discarded from this analysis, and only results obtained from carcasses in good and intermediate conservation status were considered. A negative linear regression was observed, and a significant equation was obtained [10E (-0.2986 x Cq + 13.5094)] (R2 = 0.99 and p< 2.2E-16) to estimate the YFV genomic load using the Cqs in liver samples under the experimental conditions used here.
The median Cqs differed regarding the environment where carcasses were collected (x2 = 57.19, p< 0.0001) (Table 3). Lower median Cqs (higher estimated YFV genomic loads) were observed for carcasses collected in rural areas (median Cq = 13.0) compared to carcasses from urban-rural interface (median Cq = 31.5) (Table 3), and urban areas (median Cq = 32.6) (p = 0.0026 and p< 0.0001, respectively). Regarding the genera, analyses of median Cqs among Alouatta (median Cq = 13.0), Callicebus (median Cq = 8.5), and Callithrix (median Cq = 32.0) indicated a significant difference (x2 = 31.86, p< 0.0001) (Table 3). Differences in median Cqs were observed for total Callithrix compared to Alouatta (p< 0.001) and to Callicebus (p = 0.0045). Finally, for one specimen of Sapajus nigritus, Cq of 30.2 (estimate of 3.10E+04 gc/g liver) was observed.
When we analyzed median Cqs considering both NHP genera and environments, differences were observed (x2 = 51.58, p< 0.0001) (Table 3). Rural Alouatta (median Cq = 14.1), rural Callicebus (median Cq = 13.9), and rural Callithrix (median Cq = 15.0) presented lower median Cqs (higher estimated YFV genomic loads) compared to urban Callithrix (median Cq = 33.0) (p≤ 0.0014). No difference was observed in median Cqs when Alouatta, Callicebus, and Callithrix collected in rural areas were compared to each other (p> 0.20). No difference in median Cqs was observed between urban-rural (median Cq = 33.8) and urban Callithrix (p = 0.62).
Partial nucleotide sequences of the envelope gene of 11 samples were obtained (Genbank accession numbers: MN517211-MN517221) (S2 Table), confirming wild-type YFV in NHP carcasses from urban and rural areas of different municipalities of Minas Gerais state. Phylogenetic analysis, involving 112 nucleotide (nt) sequences (621 nt) (S3 Table), indicated that YFV isolates obtained here formed a monophyletic cluster together with other YFV isolates from 2016 to 2018, within South America genotype I (S3 Fig).
Discussion
From late 2016 to 2019, Southeast Brazil faced large and extensive sylvatic YF outbreaks, with thousands of human cases and epizootics, especially in MG state [4,10–13]. Based on demand for YFV investigation, we received liver samples belonging to NHP carcasses from January/2017 to December/2018. NHP are the primary hosts of the sylvatic YF cycle, and to better understand viral circulation in Minas Gerais state during recent outbreaks, we investigated YFV RNA in NHP carcasses. YFV was detected with high overall positivity of 38% in specimens of Alouatta, Callicebus, Sapajus, and Callithrix. We observed the widespread circulation of YFV in Minas Gerais, detecting the virus in 49 municipalities where YF was not previously reported during the recent YF outbreaks [14,25,26]. In Minas Gerais state, the mesoregions with the highest numbers of YFV-positive carcasses were the same ones with the highest records of human YF cases from 2016 to 2018 [4, 14, 23, 24]. Our phylogenetic analyses confirmed the YFV in different urban areas of Minas Gerais, close to human populations. Besides, the YFV detected in urban and rural environments belonged to the same lineage, indicating the dissemination of YFV from rural/sylvatic towards the urban areas during the recent YF outbreaks.
A gradient of YFV dissemination was observed from rural towards urban-rural interface and urban areas, in conformity with the sylvatic pattern of YF [5], when one would expect a higher exposure to the virus in rural environments. According to the seasonal pattern, a higher positivity for YFV infection was observed during the epidemic periods, when a more extensive transmission of the virus and a higher number of cases would likely occur [9–11]. Our data indicate the period from December to March with the highest chances of YFV detection in NHP carcasses, what should be considered by the YF surveillance system, and strengthening the surveillance, especially at the beginning of the epidemic period. Nevertheless, the YFV circulation during the non-epidemic period of 2017 was confirmed here and by official bulletins [10,14] confirming the persistence of YFV in Minas Gerais, from 2017 to 2018 [7], followed by a decline by mid-2018 until 2019 [15].
The NHP genera and species studied here are distributed along with the Atlantic Forest and Cerrado biomes [4,27], the biomes from which NHP carcasses were sampled. Alouatta spp. has been described as extremely sensitive to YFV infection, developing severe and fatal disease [1,4]. However, specimens of Alouatta presenting antibodies against YFV shows that howler monkeys can survive from YF [4,28,29]. Indeed, high rates of howler monkeys presenting protective or neutralizing antibodies against YFV were described in studies conducted in different regions of Panama (50.7% out of 203 individuals) [29] and Brazil (25.5% out of 192 specimens) [28]. There is little information regarding the susceptibility of Callicebus spp. to YFV infection. However, humoral immunity against YFV was demonstrated in 28.6% of 14 specimens investigated by Kumm and Laemmert (1950) [28], in Brazil. In the present study, Callicebus and Alouatta presented similar proportions of YFV-infected carcasses with similar low median Cqs in RT-qPCR. These results might indicate the high sensitivity to YFV of Callicebus, as C. nigrifrons and C. personatus. The YF surveillance and animal conservation programs should consider this information since some of those species are threatened (C. nigrifrons) and vulnerable (C. personatus) with decreasing population trends [30]. Similarly, species within Sapajus genus present population decreasing trends with some already considered critically endangered [30]. Here, only two specimens were sampled, and one carcass was YFV-positive. Although Sapajus specimens may develop severe disease, they have been considered less susceptible to YFV [4]. Our results and previous studies [5,7,31] showed YF outbreaks causing thousands of deaths of NHP, which lives in one of the most threatened ecosystems of the world, the Atlantic Forest [4,18].
Despite many marmoset carcasses received (mostly from urban areas), the overall positivity for YFV in Callithrix was lower than Alouatta and Callicebus. Due to their ability to adapt and survive in degraded and secondary habitats [32], free-living marmosets frequently inhabit urban green areas inside Brazilian cities, close to human populations. The proximity of Callithrix to the anthropic environment and its wide occurrence in urban areas could increase the chances of dead marmosets to be found in a good preservation status. However, we observed a strong negative correlation for detection of YFV in urban marmosets and higher median Cqs compared to Callicebus and Alouatta, suggesting lower genomic viral loads in those urban marmosets. Previous studies showed similar data when Callithrix presented higher Cqs (lower viral loads) than Alouatta during the recent YFV epizootics in Brazil [33]. On the other hand, we did not observe a difference among rural NHP (Callithrix, Alouatta, and Callicebus), which presented median estimates of genomic viral loads of 10E+6 times higher than Callithrix from urban and urban-rural transition areas. When only marmosets were evaluated, we detected a higher positivity in rural areas and a decreasing gradient towards urban areas, consistent with the sylvatic pattern of YF. In that way, the lower overall positivity for YFV could result from the lower exposure to YFV-infected sylvatic vectors in urban areas. Haemagogus spp. is highly competent and considered the primary vector of YFV in Brazil [1,4,34], including the recent outbreaks in the country [34,35]. These mosquitoes are primatophilic species, being found in the tree canopy [36]. However, some species as Haemagogus leucocelaenus and Haemagogus janthynomis can disseminate through great distances [34] and may adapt to modified environments, as urban-rural interface [34,35] and urban green areas [37].
In the present study, we observed YFV-positive carcasses collected in urban areas boundaries and different neighborhoods, including downtown areas or near to the central area of different cities. For example, 25-YFV positive marmosets were collected in urban areas of Belo Horizonte (2.5 million inhabitants), the state capital, and a total of 76 YFV-positive marmosets were sampled in the Metropolitan Region of Minas Gerais state (5.0 million inhabitants). Due to some carcasses location (for example, near the downtown area of Belo Horizonte), the land use, and marmosets territorial behavior [32], some of these NHP would have been likely infected in urban matrixes.
Although NHP have a crucial role in the epidemiology of YF and despite YFV-infected NHP in urban areas, there was no documented case of urban/Aedes transmission of YFV in Brazil since 1942 [1,4,9–11]. The lower viral load (represented by the estimated viral genomic load) in urban NHP could lower the chances for further viral transmission to vectors. Besides, the closure of some parks in urban centers, the vaccination of humans, and the distinct behavior patterns and habitats observed for hosts (NHP and humans) and vectors in urban areas might have played a role in preventing the urban/Aedes transmission of YFV to humans. However, the high infestation of Brazilian urban areas by Aedes aegypti poses a risk for YFV urban/Aedes transmission [4,34,37,38]. Concerning Aedes, previous studies in Brazil have confirmed the competence of Aedes aegypti and Aedes albopictus, although the latter has shown low transmission competency to YFV [33].
Here we demonstrated the wide occurrence of YFV-infected NHP in Minas Gerais state, including several urban areas, during the epidemic and non-epidemic periods. Although the number of YFV-infected NHP and humans dramatically decreased from the mid-2018, YFV persisted in the Southeast region [5,7,8], at least during three consecutive transmission seasons [39]. These data reinforce that YFV has suitable ecological and climate conditions for its maintenance in Southeast Brazil. Further studies to investigate the sensitivity of urban NHP and the competence of YFV vectors found in urban and urban-rural transition areas should be conducted to understand better the dynamics of YF and the risks for the occurrence of YF in urban centers.
Supporting information
S1 Fig. Characteristics of non-human primate carcasses used for yellow fever virus investigation.
(A) Balloon plot of observed frequencies (freq) of NHP carcasses related to the year (2017 and 2018) and the sampling period. Epi: epidemic period (December to May), non-epi: non-epidemic period (June to November). (B) Balloon plot of observed frequencies (freq) of NHP carcasses related to the genera Alouatta, Callicebus, Callithrix, or non-identified specimens (Non-ID), and the environment rural, urban-rural interface (UR-I), or urban areas. Balloon plots were created using Rand; each cell contains a dot whose size and color reflect the relative magnitude of the corresponding component, according to the frequencies presented. (C) Distribution of NHP carcasses according to the epidemiological weeks, in 2017 (left) and 2018 (right). YFV-positive samples are shown in red, and samples non-detectable for YFV RNA by RT-qPCR [24] are shown in blue. Epidemic periods, in 2017 and 2018, correspond to epidemiological weeks 1 to 22 and 49 to 52 of each year.
https://doi.org/10.1371/journal.pntd.0008658.s001
(TIF)
S2 Fig. Association between detection of yellow fever virus (YFV) in the liver of non-human primate carcasses with each bimester within the yellow fever epidemic and non-epidemic periods.
YFV ND: yellow fever virus RNA non-detectable. YFV pos: yellow fever virus RNA detectable. Pearson residuals (standardized) were extracted from the chi-square function and plotted. In each cell, the circle’s size is proportional to the amount of cell contribution, and colors indicate positive residuals (blue) or negative residuals (red). The vertical bars indicate Pearson residuals values. Analyses were run in R software v.3.6.0.
https://doi.org/10.1371/journal.pntd.0008658.s002
(TIF)
S3 Fig. Phylogenetic tree of yellow fever virus based on partial envelope gene sequence.
The phylogenetic tree was reconstructed using the Maximum Likelihood method based on the General Time Reversible model, using Gamma distribution (5 categories). A total of 1000 bootstrap replicates were run, and circles in the respective nodes represent the results. The analysis involved 112 nucleotide (nt) sequences (621 nt). Branch lengths are drawn to a scale of nucleotide substitutions per site according to the scale. Sequences are identified by Genbank accession. The analysis was conducted in MEGA7, and the final tree was edited and visualized in FigTree v.1.4.4.
https://doi.org/10.1371/journal.pntd.0008658.s003
(TIF)
S1 Table. Non-human primate (NHP) carcasses collected and tested for yellow fever virus (YFV) RNA, according to the mesoregions of Minas Gerais state, Brazil.
YFV: yellow fever virus. RT-qPCR: one-step real-time polymerase chain reaction. NHP: non-human primate. Liver samples of NHP carcasses collected in Minas Gerais state, Brazil (January 2017–December 2018), were tested for YFV RNA using the RT-qPCR [24]. (+) detection of YFV RNA by RT-qPCR. (-) non-detection of YFV RNA by RT-qPCR.
https://doi.org/10.1371/journal.pntd.0008658.s004
(DOC)
S2 Table. Information regarding non-human primate (NHP) carcasses from which the yellow fever virus sequences were obtained.
aAll municipalities are in Minas Gerais state, Southeast of Brazil. NHP: non-human primate. ID: identification. Non-ID: non-identified specimens. Feb.: February. Apr.: April. Jan.: January.
https://doi.org/10.1371/journal.pntd.0008658.s005
(DOC)
S3 Table. Information of yellow fever virus (YFV) sequences included in the dataset.
ID: identification. YFV: yellow fever virus. NHP: non-human primate. NA: not available. This dataset included 112 YFV nucleotide (nt) sequences, spanning 621 nt (from the nucleotide 1,766 to 2,386 of MK333808.1 sequence) from South American and African genotypes.
https://doi.org/10.1371/journal.pntd.0008658.s006
(DOC)
References
- 1. da C Vasconcelos PF. Yellow Fever. Rev Soc Bras Med Trop. 2003;36: 275–93. pmid:12806465
- 2. Douam F, Ploss A. Yellow Fever Virus: Knowledge Gaps Impeding the Fight Against an Old Foe. Trends Microbiol. 2018;26: 913–28. pmid:29933925
- 3. da C Vasconcelos PF. Yellow fever in Brazil: thoughts and hypotheses on the emergence in previously free areas. Rev Saude Publica. 2010;44: 1144–9. pmid:21109907
- 4. Silva NIO, Sacchetto L, de Rezende IM, de S Trindade G, LaBeaud AD, de Thoisy Benoit and Drumond BP. Recent sylvatic yellow fever virus transmission in Brazil: the news from an old disease. Virol J. 2020; https://doi.org/10.1186/s12985-019-1277-7 pmid:31973727
- 5. Faria NR, Kraemer MUG, Hill SC, De Jesus JG, Aguiar RS, Iani FCM, et al. Genomic and epidemiological monitoring of yellow fever virus transmission potential. Science (80-). 2018;361: 894–99. pmid:30139911
- 6. Gómez MM, de Abreu FVS, Dos Santos AAC, de Mello IS, Santos MP, Ribeiro IP, et al. Genomic and structural features of the yellow fever virus from the 2016–2017 Brazilian outbreak. J Gen Virol. 2018;99: 536–548. pmid:29469689
- 7. de Rezende IM, Sacchetto L, Munhoz de Mello É, Alves PA, de M Iani FC, Adelino TÉR, et al. Persistence of Yellow fever virus outside the Amazon Basin, causing epidemics in Southeast Brazil, from 2016 to 2018. PLoS Negl Trop Dis. 2018;12: e0006538. pmid:29864115
- 8. Delatorre E, Santos De Abreu FV, Ribeiro IP, Gómez MM, Cunha Dos Santos AA, Ferreira-De-Brito A, et al. Distinct YFV Lineages Co-circulated in the Central-Western and Southeastern Brazilian Regions from 2015 to 2018. Front Microbiol. 2019;10. pmid:31178835
- 9.
SVS-MS. Centro de operações de emergência em saúde pública sobre febre amarela Informe No 43 2017. In: Secretaria de Vigilância em Saúde—Ministério da SaúdeMinistério da Saúde [Internet]. 2017 [cited 7 Jul 2019] pp. 1–7. http://portalarquivos.saude.gov.br/images/pdf/2017/junho/02/COES-FEBRE-AMARELA—INFORME-43—Atualiza—-o-em-31maio2017.pdf
- 10.
SVS-MS. Monitoramento do Período Sazonal da Febre Amarela Brasil 2017/2018 Informe N° 27. In: Secretaria de Vigilância em Saúde—Ministério da Saúde [Internet]. 2018 [cited 6 Jul 2019] pp. 1–12. https://portalarquivos2.saude.gov.br/images/pdf/2018/outubro/08/Informe-FA.pdf
- 11.
SVS-MS. Monitoramento de Febre Amarela Brasil 2019 Informe No 18. In: Secretaria de Vigilância em Saúde—Ministério da Saúde [Internet]. 2019 [cited 6 Jul 2019] pp. 1–8. https://portalarquivos2.saude.gov.br/images/pdf/2019/junho/13/Informe-de-Monitoramento-de-Febre-Amarela-Brasil—n-18.pdf
- 12.
MS-BR. Monitoramento do período sazonal da Febre Amarela Brasil 2018/2019 Informe No 3. In: Ministério da Saúde do Brasil [Internet]. 2019 [cited 5 Jul 2019] pp. 1–12. https://portalarquivos2.saude.gov.br/images/pdf/2019/janeiro/28/informe-FA-n.3-21jan19.pdf
- 13.
MS-BR. Monitoramento do período sazonal da Febre Amarela Brasil 2017/2018 Informe No 1 2017/2018. In: Ministério da Saúde do Brasil [Internet]. 2017 [cited 5 Jul 2019] pp. 1–12. https://portalarquivos2.saude.gov.br/images/PDF/2017/novembro/14/Informe-FA-14-11-17.pdf
- 14.
SES-MG. Febre Amarela Silvestre em Minas Gerais 18/10/2018. In: Secretaria do Estado de Saúde de Minas Gerais. [Internet]. 2018 [cited 6 Jul 2019] pp. 1–12. http://www.saude.mg.gov.br/images/noticias_e_eventos/000_2018/BoletinsEpidemiologicos/Boletim_-_Febre_Amarela_2018_-_Consolidado_V3.pdf
- 15.
SES-MG. Febre Amarela Silvestre em Minas Gerais—19/02/2019. In: Secretaria de Estado de Saúde de Minas Gerais [Internet]. 2019 [cited 5 Jul 2019] pp. 1–12. http://www.saude.mg.gov.br/images/noticias_e_eventos/000_2019/jane_fev_mar/Febre_Amarela/Boletim_atualização_FA_12-02-2019.pdf
- 16.
SVS-MS. Febre Amarela Febre do Nilo Ocidental Febre Maculosa Brasileira e Outras Riquetsioses. Guia de Vigilância em Saúde. 4a edição. Brasília: Secretaria de Vigilância em Saúde, Ministério da Saúde; 2019. p. 727.
- 17. Monath TP, Vasconcelos PFC. Yellow fever. J Clin Virol. Elsevier B.V.; 2015;64: 160–173. pmid:25453327
- 18. Stan K, Sanchez-Azofeifa A, Espírito-Santo M, Portillo-Quintero C. Simulating deforestation in Minas Gerais, Brazil, under changing government policies and socioeconomic conditions. PLoS One. 2015;10: e0137911. pmid:26371876
- 19.
IBGE. Sinopse Preliminar do Censo Demográfico. In: Instituto Brasileiro de Geografia e Estatística [Internet]. 2000 [cited 2 Jun 2020] p. 415. https://biblioteca.ibge.gov.br/visualizacao/periodicos/308/cd_2000_v7.pdf
- 20. Gregorin R. Taxonomia e variação geográfica das espécies do gênero Alouatta Lacépède (Primates, Atelidae) no Brasil. Rev Bras Zool. 2006;23: 64–144.
- 21.
Groves C. Smithsonian series in comparative evolutionary biology. Primate taxonomy. Washington, DC: Smithsonian Institution Press; 2001. pp. 1–350.
- 22. van Roosmalen MGM, van Roosmalen T, Mittermeier R. A taxonomic review of the titi monkeys, genus Callicebus Thomas, 1903, with the description of two new species, Callicebus bernhardi and Callicebus stephennashi, from Brazilian Amazonia. Neotrop Primates. 2002;10: 1–52.
- 23. de Rezende IM, Alves PA, Arruda MS, Gonçalves AP, Oliveira GFG, Pereira LS, et al. Yellow fever virus genotyping tool and investigation of suspected adverse events following yellow fever vaccination. Vaccines. 2019;7. pmid:31817103
- 24. Domingo C, Patel P, Yillah J, Weidmann M, Méndez JA, Nakouné ER, et al. Advanced yellow fever virus genome detection in point-of-care facilities and reference laboratories. J Clin Microbiol. 2012;50: 4054–60. pmid:23052311
- 25.
SES-MG. Febre Amarela Silvestre em Minas Gerais—20/06/2018. In: Secretaria de Estado de Saúde de Minas Gerais [Internet]. 2018 [cited 5 Jul 2019] pp. 1–12. http://www.saude.mg.gov.br/images/documentos/Boletim_Febre%20Amarela_21.06.2018_atualizada.pdf
- 26.
SES-MG. Situação epidemiológica da febre amarela silvestre em Minas Gerais, 2017–21/08/2017. In: Secretaria de Estado de Saúde de Minas Gerais [Internet]. 2017 [cited 5 Jul 2019] pp. 1–29. http://www.saude.mg.gov.br/images/editais/AtualizaoFA-22%20de%20agosto%202017.pdf
- 27.
Paglia AP, Fonseca GAB, Rylands AB, Herrmann G, Aguiar LMS, Chiarello AG, et al. Lista Anotada dos Mamíferos do Brasil/Annotated Checklist of Brazilian Mammals. (2a Edição/2nd Edition). Occasional Papers in Conservation Biology. Arlington, VA: Conservation International; 2012.
- 28. Kumm HW, Laemmert HW. The geographical distribution of immunity to yellow fever among the primates of Brazil. Am J Trop Med Hyg. 1950;30: 733–48. pmid:14771397
- 29. Clark HC. Endemic yellow fever in Panama and neighboring areas. Am J Trop Med Hyg. 1952;1: 78–86. pmid:14903438
- 30.
IUCN. The IUCN Red List of Threatened Species. In: International Union for Conservation of Nature and Natural Resources [Internet]. 2019 [cited 5 Jul 2019]. https://www.iucnredlist.org/
- 31. Cunha MS, da Costa AC, de Azevedo Fernandes NCC, Guerra JM, dos Santos FCP, Nogueira JS, et al. Epizootics due to Yellow Fever Virus in São Paulo State, Brazil: viral dissemination to new areas (2016–2017). Sci Rep. 2019;9. pmid:30940867
- 32. Teixeira B, Hirsch A, Goulart VDLR, Passos L, Teixeira CP, James P, et al. Good neighbours: Distribution of black-tufted marmoset (Callithrix penicillata) in an urban environment. Wildl Res. 2015;42: 579–89.
- 33. Mares-Guia MAMDM, Horta MA, Romano A, Rodrigues CDS, Mendonça MCL, Dos Santos CC, et al. Yellow fever epizootics in non-human primates, Southeast and Northeast Brazil (2017 and 2018). Parasites and Vectors. 2020;13: 90. pmid:32075684
- 34. de Abreu FVS, Ribeiro IP, Ferreira-de-Brito A, dos Santos AAC, de Miranda RM, de S Bonelly I, et al. Haemagogus leucocelaenus and Haemagogus janthinomys are the primary vectors in the major yellow fever outbreak in Brazil, 2016–2018. Emerg Microbes Infect. 2019;8: 218–31. pmid:30866775
- 35. Pinheiro GG, Rocha MN, de Oliveira MA, Moreira LA, Filho JDA. Detection of yellow fever virus in sylvatic mosquitoes during disease outbreaks of 2017–2018 in minas gerais state, brazil. Insects. 2019;10: e136. pmid:31083286
- 36. Causey OR, Kumm HW, Laemmert HW. Dispersion of forest mosquitoes in Brazil: Further studies. Am J Trop Med Hyg. 1950;s1-30: 301–12. pmid:15419408
- 37. Mucci LF, Medeiros-Sousa AR, Ceretti-Júnior W, Fernandes A, Camargo AA, Evangelista E, et al. Haemagogus leucocelaenus and Other Mosquitoes Potentially Associated With Sylvatic Yellow Fever In Cantareira State Park In the São Paulo Metropolitan Area, Brazil. J Am Mosq Control Assoc. 2016;32: 329–32. pmid:28206862
- 38. Mir D, Delatorre E, Bonaldo M, Lourenço-De-Oliveira R, Vicente AC, Bello G. Phylodynamics of Yellow Fever Virus in the Americas: New insights into the origin of the 2017 Brazilian outbreak. Sci Rep. 2017;7: 7385. pmid:28785067
- 39. De Abreu FVS, Delatorre E, Dos Santos AAC, Ferreira-De-Brito A, De Castro MG, Ribeiro IP, et al. Combination of surveillance tools reveals that yellow fever virus can remain in the same atlantic forest area at least for three transmission seasons. Mem Inst Oswaldo Cruz. 2019;114: e190076. pmid:31038550