- Open Access
Modeling the elevated risk of yellow fever among travelers visiting Brazil, 2018
Theoretical Biology and Medical Modellingvolume 15, Article number: 9 (2018)
Unlike the epidemic of yellow fever from 2016 to 17 in Brazil mostly restricted to the States of Minas Gerais and Espirito Santo, the epidemic from 2017 to 18 mainly involved São Paulo and Rio de Janeiro and resulted in multiple international disseminations. To understand mechanisms behind this observation, the present study analyzed the distribution of imported cases from Brazil, 2018.
A statistical model was employed to capture the risk of importing yellow fever by returning international travelers from Brazil. We estimated the relative risk of importation among travelers by the extent of wealth measured by GDP per capita and the relative risk obtained by random assignment of travelers’ destination within Brazil by the relative population size.
Upper-half wealthier countries had 2.1 to 3.4 times greater risk of importation than remainders. Even among countries with lower half of GDP per capita, the risk of importation was 2.5 to 2.8 times greater than assuming that the risk of travelers’ infection within Brazil is determined by the regional population size.
Travelers from wealthier countries were at elevated risk of yellow fever, allowing us to speculate that travelers’ local destination and behavior at high risk of infection are likely to act as a key determinant of the heterogeneous risk of importation. It is advised to inform travelers over the ongoing geographic foci of transmission, and if it appears unavoidable to visit tourist destination that has the history of producing imported cases, travelers must be strongly advised to receive vaccination in advance.
Yellow fever virus that belongs to Flavivirus causes yellow fever, transmitted via Aedes species from infected human to human . In addition to human-mosquito-human transmission cycle, nonhuman primates are also infected with the virus, and such transmission cycle has been known to be responsible for allowing continued transmission in Brazil [2,3,4]. The exposure to this virus mostly results in asymptomatic infection, but a part of cases develop fever, headache, chills, muscle pain, nausea and vomiting following an incubation period of 3–6 days . If exacerbated, the case fatality risk (CFR) given severe clinical disease is known to be 47% with the range from 40 to 80% , and there is no specific treatment. Thus, immunization is the mainstream of countermeasures , and residents of high-risk area and travelers visiting those areas are recommended to undertake vaccination.
The epidemic of yellow fever in Brazil from December 2016 to June 2017 involved 777 confirmed cases, all reported among Brazilian residents, and the CFR was estimated at 34% with 261 deaths . While the end of epidemic was declared once in September 2017, sporadic cases continued and a surge of cases started from December 2017 . On 16 January 2018, World Health Organization (WHO) recognized the ongoing epidemic and recommended vaccination among all residents of the entire State of Rio de Janeiro and São Paulo, and vaccination campaign by the government of Brazil took place with targeting States Bahia, Rio de Janeiro, São Paulo . As of 8 June 2018, 1257 confirmed cases and 394 deaths have been reported .
While no imported case was reported during the epidemic from 2016 to 17, multiple importation events have been reported during the epidemic from 2017 to 18. There have been 12 imported cases, as of 8 June 2018, notified in 8 different countries since December 2017 [9, 10]. As a possible mechanism for observing multiple importation events, the epidemic location in 2017–18 has involved Rio de Janeiro and São Paulo more than the 2016–17 epidemic . The present study aims to quantify the risk of infection among travelers visiting Brazil, 2018.
To estimate the risk of yellow fever among travelers, we analyze both the confirmed cases in Brazil from 2017 to 18  and imported cases reported abroad. As of 8 May 2018, there were 3 imported cases from Chile , 3 from Argentina , and 1 case from the Netherlands , Switzerland , France , the United Kingdom , Romania  and Germany , respectively. Nationality of all imported cases was their own country, and all these cases were regarded as importation due to travel to Brazil. Except for cases in the Netherlands and France and one case in Argentina, imported cases shared a history of visit to Ilha Grande, municipality of Angra do Reis, State of Rio de Janeiro, Brazil .
To calculate the expected risk of importation, inbound travel volume ci from each country i to Brazil was retrieved from the World Tourism Organization . In addition, we used the relative value of the Gross Domestic Product (GDP) per capita, gi of country i that was normalized by the maximum GDP in 2016 for the sake of imputation of vaccination coverage (see below) . This imputation was partly validated by statistical analysis of an association between GDP per capita and the risk of importation by country. Confirming that the variance is not significantly different between two groups by F-test, we employed Student t-test to compare the GDP per capita between countries with and without imported cases. Moreover, vaccination coverage vi of country i in 2015 was partly retrieved from a published study .
Following Dorigatti et al. , we model the expected number E(ci) of imported yellow fever cases ci in country i as
In this equation, ni is the yearly inbound number of travelers visiting Brazil from country i, vi is the vaccination coverage, qi represents the relative risk of infection among travelers from country i, which we would like to estimate through this exercise. pops and popB represent the population sizes of the three major states of 2017–18 epidemic (i.e., Minas Gerais, Rio de Janeiro and São Paulo) and the entire Brazil, respectively (pops = 81,230,574 and popB = 202,768,562 persons). The reported number of confirmed yellow fever cases in affected states is cS. We do not use the undiagnosed factor of 10, which was adopted elsewhere [5, 18], because we estimate the expected number of confirmed imported cases in the abroad. w is the mean length of stay in Brazil (w = 17 days), and fs and ps are the normalized monthly frequency of cases and the volume of travelers of month s, respectively (where fs from December 2017 to March 2018 accounted for 98.4% and ps during the same period is 47.6% of the total). Although not specified in eq. (1), month s was integrated from December 2017 up until 16 March 2018 where censoring in March was incorporated by accounting for the number of days, i.e., 16/31. μE and μI represent the mean latent and infectious period, respectively, assumed as 4.6 and 4.5 days. It should be noted that eq. (1) is intact even when we account for unascertained/asymptomatic fraction of cases. Namely, supposing that the confirmation probability among all infected individuals is α, both sides of eq. (1) is divided by α to express everything as the total number of infected individuals, and then, the constant 1/α is cancelled out from both sides.
we employ the zero-inflated Poisson distribution to describe the observed frequency of imported cases in country i, i.e.,
The mean number of imported cases from the eq. (3) is E(ci) = (1-vi)mi(qi). With respect to the vaccination coverage vi among travelers from country i, we model it as
The logit transformation was adopted, because it does not require additional parameters. The estimated vaccine coverage in a part of the routinely immunized countries was available in the published study . Countries with known vaccination coverage included Trinidad and Tobago, Panama, Argentina, Colombia, Suriname, Peru, Venezuela, Ecuador, Paraguay, Guyana, Bolivia, Angola, Nigeria, Ghana and Kenya. If the information of vaccination coverage was unavailable, we extrapolated the vaccination coverage in country i. Considering that only a small fraction of travelers is vaccinated in those countries, k was assumed to be 0.10, which is the so-called carrying capacity of the logistic distribution, and practically interpreted as theoretical possible maximum of the vaccination coverage. gi is the relative GDP per capita of country i compared to the country with the highest GDP per capita as stated above. The second case in the eq. (4) indicates that we imposed an assumption that, if no routine immunization takes place, at most 10% of those travelers visiting Brazil received vaccination, and also another assumption that the coverage follows a logit transformation of the relative GDP per capita. Because the ceiling of coverage k is a strong assumption, we varied it from 0.01 (1%) to 0.90 (90%) as part of sensitivity analysis.
We estimate qi by the list of country, i.e. by first and second half of GDP per capita, because the propensity to visit high-risk area of infection, which coincided with a resort area in the ongoing epidemic [10, 12], may vary with the extent of wealth of that particular country. qi is modelled as
where q is the relative risk of infection among travelers visiting from countries with the second half of GDP per capita as compared with an assumption that the travelers’ destination is randomly determined according to the regional population size of Brazil. a is the relative risk of importing yellow fever among wealthier countries compared to the reminders with the risk q.
Given observed counts of imported cases, c, written as a vector representing the input from all countries at risk of infection, maximum likelihood estimates of a and q were found by minimizing the negative logarithm of the following likelihood
The 95% confidence intervals (CI) were computed using the profile likelihood.
The present study analyzed data that is publicly available. As such, the datasets used in our study were de-identified and fully anonymized in advance, and the analysis of publicly available data without identity information does not require ethical approval.
Figure 1 shows the epidemic curve of 2017–18 epidemic in Brazil as a function of the week of report . The highest incidence was reported in Week 3 of 2018 followed by mostly a monotonic decline in incidence, and the incidence has greatly waned by Week 17, 2018. Figure 2 shows the comparison of GDP per capita by countries with and without imported cases. GDP per capita of countries with imported cases (n = 8) was 40,213 US dollars (95% CI: 25,614, 54,811), while that of countries without imported cases and with direct flight link from Brazil (n = 78) was 26,368 US dollars (95% CI: 21,593, 31,043). It appears that countries with imported cases have significantly greater GDP per capita than those without imported cases (t = 2.34, p = 0.04 by Student t-test).
In total, 12 imported cases were reported from 8 different countries (Fig. 3a). Other 78 countries with inbound data to Brazil were included in the following analysis. Known vaccinated fractions were as follows: Trinidad and Tobago (0.96), Panama (0.73), Argentina (0.94), Colombia (0.91), Suriname (0.89), Peru (0.90), Venezuela (0.87), Ecuador (0.78), Paraguay (0.80), Guyana (0.95), Bolivia (0.89), Angola (0.80), Nigeria (0.74), Ghana (0.89) and Kenya (0.78). Assuming k = 0.10, travelers from wealthier fraction of countries were a = 2.3 (95% CI: 0.7, 8.6) times more likely to be infected with yellow fever compared with countries below median GDP per capita. Compared with an assumption that the risk of infection was determined by the relative population size of epidemic locations (i.e., 3 states with substantial number of cases) to the entire Brazil (i.e. eq. (1) with q = 1.0), the assumption of which was employed elsewhere , even countries below median GDP per capita experienced q = 2.5 (95% CI: 0.8, 5.9) times greater risk of yellow fever.
Figure 3b shows the result from sensitivity analysis. Varying an uncertain parameter k, i.e., the possible maximum value of the vaccination coverage among travelers from countries without routine immunization against yellow fever, the estimates of qi did not vary greatly. The relative risk a of yellow fever among wealthier countries compared with countries with second half of GDP per capita ranged from 2.1 with k = 0.01 to 3.4 with k = 0.90. Also, compared with random assignment of travelers’ destination by relative population size of States, countries below median GDP per capita experienced q = 2.5 (k = 0.01) to 2.8 (k = 0.90) times greater risk of yellow fever.
Unlike the epidemic from 2016 to 17 in Brazil that was mostly restricted to the States of Minas Gerais and Espirito Santo, the epidemic from 2017 to 18 mainly involved São Paulo and Rio de Janeiro and resulted in multiple international disseminations of imported cases. To understand possible mechanisms behind this observation and also to consider possible countermeasures, the present study explored the distribution of imported cases from Brazil. Employing a statistical model, we described the risk of observing imported case, jointly estimating the relative risk of travelers by the extent of wealth (or GDP per capita) and the relative difference compared with random assignment of travelers’ destination within Brazil. As a result, it appears that wealthier travelers were at 2.1 to 3.4 times greater risk of infection than others. Moreover, even among countries with lower half of GDP per capita, the risk was 2.5 to 2.8 times greater than that with the assumption that the relative risk within Brazil is determined by regional population size.
There are two take home messages. First, we have shown that countries with wealthier GDP per capita appeared to be more often infected. The finding is in line with the fact that the imported cases arose from a holiday spot in Ilha Grande, municipality of Angra do Reis, State of Rio de Janeiro [9, 10]. It also indicates that travelers’ local destination and behavior at high risk of infection are likely to act as a key determinant of the heterogeneous risk of importing case. It is advised to well inform travelers over the ongoing geographic foci of transmission, and if it appears unavoidable to visit tourist destination that has the history of producing imported cases, travelers must be strongly advised to receive vaccination in advance.
Second, we found that even non-wealthy countries were at 2.5–2.8 times greater risk of importing yellow fever case as compared with a common modeling assumption (i.e., q = 1.0 in eq. (1)) that the destination-specific risk of infection is proportional to the relative population size of the destination to the entire country. In the case of Brazil, undoubtedly the major tourist destinations of international travelers are São Paulo and Rio de Janeiro. To precisely estimate the risk of infection among travelers, it is ideal to track down travel patterns within Brazil more in detail. A big challenge to achieve precise estimation of the risk in the future would be to quantify such risk in a finer spatial scale using limited mobility information among travelers.
Four limitations must be noted. First, the notification of yellow fever cases is undoubtedly biased by the extent of ascertainment. Thus, even though we found that travelers from countries with greater GDP per capita were at greater risk of yellow fever, the finding could partly reflect better ascertainment of cases in wealthier countries compared with the reminder. Second, we were not able to account for spatial risk of infection in a finer scale. As of 8 May 2018, the transmission has not been established within the city of Rio de Janeiro , and thus, such risk at greater precision must be communicated with risk map, as it was published elsewhere [21,22,23]. Third, the vaccination coverage among travelers from countries without routine yellow fever immunization was assumed to be proportional to GDP per capita. In the present study, our estimates were not sensitive to the ceiling of the vaccination coverage, k, but this strong assumption needs to be validated through empirical observation in the future. Fourth, the lower confidence bound of our relative risk estimates were smaller than the value of 1 (e.g., with k = 0.10, the lower 95% CI of a was 0.6), and the sample size was not substantial. This is due to limited number of countries that imported yellow fever. With greater sample size in a future follow-up study, uncertainties would be reduced, and moreover, our conclusions would be strengthened.
Apart from these future tasks for finer estimation of the risk of infection among travelers, we believe that our study successfully quantified the relative risk of infection by GDP per capita and also compared with the risk that rests on population-size specific assumption of travelers’ destination. Micro-geographic information of imported cases should be effectively shared with travelers for communication and prevention purposes.
Travelers from wealthier countries were at elevated risk of yellow fever, allowing us to speculate that travelers’ local destination and behavior at high risk of infection are likely to act as a key determinant of the heterogeneous risk of importation. As part of important messages as derived in real time studies [19, 24,25,26,27], it is advised to inform travelers over the ongoing geographic foci of transmission, and if it appears unavoidable to visit tourist destination that has the history of producing imported cases, travelers must be strongly advised to receive vaccination in advance.
Gross Domestic Product
World Health Organization
Monath TP. Yellow fever: an update. Lancet Infect Dis. 2001;1:11–20.
Dexheimer Paploski IA, Souza RL, Tauro LB, Cardoso CW, Mugabe VA, Pereira Simões Alves AB, de Jesus Gomes J, Kikuti M, Campos GS, Sardi S, Weaver SC, Reis MG, Kitron U, Ribeiro GS. Epizootic outbreak of yellow fever virus and risk for human disease in Salvador, Brazil. Ann Intern Med. 2018;168(4):301–2.
Fernandes NCCA, Cunha MS, Guerra JM, Réssio RA, Cirqueira CDS, Iglezias SD, de Carvalho J, Araujo ELL, Catão-Dias JL, Díaz-Delgado J. Outbreak of yellow fever among nonhuman Primates, Espirito Santo, Brazil, 2017. Emerg Infect Dis. 2017;23(12):2038–41.
Moreira-Soto A, Torres MC, Lima de Mendonça MC, Mares-Guia MA, Dos Santos Rodrigues CD, Fabri AA, Dos Santos CC, Machado Araújo ES, Fischer C, Ribeiro Nogueira RM, Drosten C, Sequeira PC, Drexler JF, Bispo de Filippis AM. Evidence for multiple sylvatic transmission cycles during the 2016-2017 yellow fever virus outbreak, Brazil. Clin Microbiol Infect 2018; in press. doi: https://doi.org/10.1016/j.cmi.2018.01.026.
Johansson MA, Vasconcelos PF, Staples JE. The whole iceberg: estimating the incidence of yellow fever virus infection from the number of severe cases. Trans R Soc Trop Med Hyg. 2014;108(8):482–7.
Wisseman CL, Sweet BH, Kitaoka M, Tamiya T. Immunological studies with group B arthropod-borne viruses. I. Broadened neutralizing antibody spectrum induced by strain 17D yellow fever vaccine in human subjects previously infected with Japanese encephalitis virus. Am J Trop Med Hyg. 1962;11:550–61.
Pan American Health Organization. Yellow fever. Epidemiological update on 12 January 2018. Washington D.C.: Pan American Health Organization; 2018. Available at: https://www.paho.org/hq/index.php?option=com_content&view=article&id=14043%3A12-january-2018-yellow-fever-epidemiological-update&catid=2103%3Arecent-epidemiological-alerts-updates&Itemid=42346&lang=en (Accessed: 23 Jan 2018)
World Health Organization. Disease outbreak news on 22 January 2018. Yellow fever – Brazil. Geneva: World Health Organization. Available at: http://www.who.int/csr/don/22-january-2018-yellow-fever-brazil/en/
Ministerio da Saude, Brazil. Monitoramento do Período Sazonal da Febre Amarela. Brasil – 2017/2018. 8 May 2018. INFORME N25. São Paulo: Ministerio da Saude; 2018. Available at: http://portalarquivos2.saude.gov.br/images/pdf/2018/maio/09/Informe-FA.pdf
Gossner CM, Haussig JM, de Saint Lary CB, Aaslav KK, Schlagenhauf P, Sudre B. Increased risk of yellow fever infections among unvaccinated European travellers due to ongoing outbreak in Brazil, July 2017 to March 2018. Eurosurveillance. 2018;23(11):18–00106. https://doi.org/10.2807/1560-7917.ES.2018.23.11.18-00106
Possas C, Martins RM, Oliveira RL, Homma A. Urgent call for action: avoiding spread and re-urbanisation of yellow fever in Brazil. Mem Inst Oswaldo Cruz. 2018;113:1–2.
World Health Organization. Disease outbreak news on 27 February 2018. Yellow fever – Brazil. Geneva: World Health Organization. Available at: http://www.who.int/csr/don/27-february-2018-yellow-fever-brazil/en/
World Health Organization. Disease outbreak news on 9 March 2018. Yellow fever – Brazil. Geneva: World Health Organization. Available at: http://www.who.int/csr/don/09-march-2018-yellow-fever-brazil/en/
Weigand G. YELLOW FEVER - AMERICAS (24): GERMANY ex BRAZIL (RIO DE JANEIRO). ProMED-mail, 27 March 2018. Available at: https://www.promedmail.org/post/20180327.5714395
Yearbook of Tourism Statistics dataset. Brazil. Geneva: World Tourism Organization (UNWTO); 2017. Available from http://statistics.unwto.org/publication/yearbook-tourism-statistics-2017-edition
Central Intelligence Agency. The World Factbook. Country Comparison: GDP-Per Capita (PPP). Available at: https://www.cia.gov/library/publications/the-world-factbook/rankorder/2004rank.html
Shearer FM, Moyes CL, Pigott DM, Brady OJ, Marinho F, Deshpande A, Longbottom J, Browne AJ, Kraemer MUG, O'Reilly KM, Hombach J, Yactayo S, de Araújo VEM, da Nóbrega AA, Mosser JF, Stanaway JD, Lim SS, Hay SI, Golding N, Reiner RC Jr. Global yellow fever vaccination coverage from 1970 to 2016: an adjusted retrospective analysis. Lancet Infect Dis. 2017;17(11):1209–17.
Dorigatti I, Hamlet A, Aguas R, Cattarino L, Cori A, Donnelly CA, Garske T, Imai N, Ferguson NM. International risk of yellow fever spread from the ongoing outbreak in Brazil, December 2016 to May 2017. Euro Surveill. 2017;22(28):30572.
Tsuzuki S, Lee H, Miura F, Chan YH, Jung SM, Akhmetzhanov AR, Nishiura H. Dynamics of the pneumonic plague epidemic in Madagascar, August to October 2017. Eurosruveillance. 2017;22(46) https://doi.org/10.2807/1560-7917.ES.2017.22.46.17-00710.
Couto-Lima D, Madec Y, Bersot MI, Campos SS, Motta MA, Santos FBD, Vazeille M, Vasconcelos PFDC, Lourenço-de-Oliveira R, Failloux AB. Potential risk of re-emergence of urban transmission of yellow fever virus in Brazil facilitated by competent Aedes populations. Sci Rep. 2017;7(1):4848. https://doi.org/10.1038/s41598-017-05186-3.
Shearer FM, Longbottom J, Browne AJ, Pigott DM, Brady OJ, Kraemer MUG, Marinho F, Yactayo S, de Araújo VEM, da Nóbrega AA, Fullman N, Ray SE, Mosser JF, Stanaway JD, Lim SS, Reiner RC Jr, Moyes CL, Hay SI, Golding N. Existing and potential infection risk zones of yellow fever worldwide: a modelling analysis. Lancet Global Health. 2018;6(3):e270–8.
Garske T, Van Kerkhove MD, Yactayo S, Ronveaux O, Lewis RF, Staples JE, Perea W, Ferguson NM. Yellow fever expert committee. Yellow fever in Africa: estimating the burden of disease and impact of mass vaccination from outbreak and serological data. PLoS Med. 2014;11(5):e1001638.
Hamlet A, Jean K, Perea W, Yactayo S, Biey J, Van Kerkhove M, Ferguson N, Garske T. The seasonal influence of climate and environment on yellow fever transmission across Africa. PLoS Negl Trop Dis. 2018;12(3):e0006284.
Nishiura H, Tsuzuki S, Asai Y. Forecasting the size and peak of cholera epidemic in Yemen, 2017. Future Microbiol. 2018;13:399–402.
Matsuyama R, Akhmetzhanov AR, Endo A, Lee H, Yamaguchi T, Tsuzuki S, Nishiura H. Uncertainty and sensitivity analysis of the basic reproduction number of diphtheria: a case study of a Rohingya refugee camp in Bangladesh, November-December 2017. PeerJ. 2018;6:e4583.
Nishiura H, Tsuzuki S, Yuan B, Yamaguchi T, Asai Y. Transmission dynamics of cholera in Yemen, 2017: a real time forecasting. Theor Biol Med Model. 2017;14(1):14.
Dinh L, Chowell G, Mizumoto K, Nishiura H. Estimating the subcritical transmissibility of the Zika outbreak in the State of Florida, USA, 2016. Theor Biol Med Model. 2016;13(1):20.
HN received funding support from the Japan Agency for Medical Research and Development, the Japan Science and Technology Agency (JST) CREST program (JPMJCR1413), Telecommunication Advancement Foundation, and the Japan Society for the Promotion of Science (JSPS) KAKENHI 16KT0130 and 17H04701. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
Collected datasets is publicly available and can be retrieved from Fig. 3a.
The authors are experts with interest in Infectious Disease Epidemiology and also in Theoretical Epidemiology, and the corresponding author acts as the chairperson and team leader of the Department of Hygiene, Hokkaido University Graduate School of Medicine.
Ethics approval and consent to participate
Consent for publication
The authors declare that co-author H. Nishiura is the Editor-in-Chief of Theoretical Biology and Medical Modelling. This does not alter the authors’ adherence to all the Theoretical Biology and Medical Modelling policies on sharing data and materials.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.