Figures
Abstract
With approximately half of the world’s population at risk of contracting dengue, this mosquito-borne disease is of global concern. International travellers significantly contribute to dengue’s rapid and large-scale spread by importing the disease from endemic into non-endemic countries. To prevent future outbreaks and dengue from establishing in non-endemic countries, knowledge about the arrival time and location of infected travellers is crucial. We propose a network model that predicts the monthly number of dengue-infected air passengers arriving at any given airport. We consider international air travel volumes to construct weighted networks, representing passenger flows between airports. We further calculate the probability of passengers, who travel through the international air transport network, being infected with dengue. The probability of being infected depends on the destination, duration and timing of travel. Our findings shed light onto dengue importation routes and reveal country-specific reporting rates that have been until now largely unknown. This paper provides important new knowledge about the spreading dynamics of dengue that is highly beneficial for public health authorities to strategically allocate the often limited resources to more efficiently prevent the spread of dengue.
Citation: Liebig J, Jansen C, Paini D, Gardner L, Jurdak R (2019) A global model for predicting the arrival of imported dengue infections. PLoS ONE 14(12): e0225193. https://doi.org/10.1371/journal.pone.0225193
Editor: Pierre Roques, CEA, FRANCE
Received: May 14, 2019; Accepted: October 30, 2019; Published: December 4, 2019
Copyright: © 2019 Liebig 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: The air travel data used in this study are owned by a third party, and were licensed for use under contract by International Air Travel Association (IATA)- Passenger Intelligence Services (PaxIS): http://www.iata.org/services/statistics/intelligence/paxis/Pages/index.aspx. The same data can be purchased for use by any other researcher by contacting: Phil GENNAOUI Regional Manager - Aviation Solutions (Asia Pacific) Tel: +65 6499 2314 | Mob: +65 9827 0414 [email protected] | www.iata.org Monthly dengue incidence rates for all countries for the years 2011 and 2015 are available as supplementary information files. Median population age data is publicly available from the Central Intelligence Agency at https://www.cia.gov/library/publications/resources/the-world-factbook/fields/343.html. International tourism arrival data used in this study are owned by a third party and cannot be shared publicly. The data is available for purchase from the World Tourism Organisation at https://www.e-unwto.org/action/doSearch?ConceptID=2445&target=topic. Dengue notification data from Australia is publicly available from the Australian Department of Health at http://www9.health.gov.au/cda/source/rpt_1_sel.cfm. Dengue notification data from Europe is publicly available from the European Centre for Disease Prevention and Control at http://atlas.ecdc.europa.eu/public/index.aspx. Case-based dengue notifications for Queensland cannot be shared as it contains confidential information. The authors gained access to this data in accordance with Section 284 of the Public Health Act 2005.
Funding: The authors received internal CSIRO funding. 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
The well connected structure of the global air transportation network and the steadily increasing volume of international travel has a vast impact on the rapid, large-scale spread of arboviral and other diseases [1–7]. A recent example of disease introduction to a novel region is the spread of the Zika virus from Brazil to Europe, the United States and other countries, which prompted the World Health Organisation (WHO) to announce a public health emergency of international concern in early 2016. Investigations confirmed that international viraemic travellers were a major contributing factor to the rapid spread [8].
With an estimated 50-100 million symptomatic infections each year [9, 10], dengue is ranked the most important mosquito-borne disease [11, 12]. The rapid geographic spread is, to a great extent, driven by the increase in international air travel [13, 14]. In addition, dengue is severely under-reported, making it extremely challenging to monitor and prevent the spread of the disease. Presumably, 92% of symptomatic infections are not reported to health authorities [10]. Low reporting rates can have many reasons, including low awareness levels and misdiagnosis [9, 15].
Due to the rapid global spread of dengue as well as severe under-reporting, many countries are facing the threat of ongoing local transmission in the near future [11]. In non-endemic countries, local outbreaks are usually triggered by an imported case [16], a person who acquired the disease overseas and transmitted the virus to local mosquitoes. To prevent ongoing dengue transmission in non-endemic countries, it is critical to forecast the importation of disease cases into these areas and move from responsive containment of dengue outbreaks to proactive outbreak mitigation measures.
The majority of existing models forecast relative rather than absolute risk of dengue importation and are unable to predict the total number of imported disease cases [13, 17, 18]. The few models that can predict absolute numbers are region-specific rather than global [19–21]. The most recently proposed model estimates the total number of imported dengue cases for 27 European countries [21], however, the model has several limitations: (i) Monthly incidence rates were based on dengue cases reported to the World Health Organisation (WHO) despite dengue being under-reported and the general consensus that the actual number of cases is much higher than the figures published by the WHO [9, 10]; (ii) Only 16 countries were considered as possible sources of importation. The authors reason that these 16 countries contribute 95% of all global dengue cases, referring to numbers published by the WHO. Since African countries do not report to the WHO, and dengue remains an under-reported disease in many other countries [22–25], it is likely that the percentage contribution to the number of global dengue cases by the 16 selected countries is strongly biased; (iii) Seasonal distributions of dengue cases were inferred based on information from only two source countries (Latin American countries were assumed to have similar seasonalities to Brazil, while Thailand served as a proxy for countries in South-East Asia). The assertion that all countries within a given global region experience similar seasonal fluctuations in dengue infections is likely inaccurate. For example, dengue notifications peak between April and December in Thailand, while Indonesia reports the highest number of dengue cases from November to April [26].
The contribution of this paper is twofold. First, we develop a network model that overcomes the limitations of previous models by employing global air passenger volumes, country-specific dengue incidence rates and country-specific temporal infection patterns. We construct weighted directed networks, using data collected by the International Air Transportation Association (IATA) to capture the movement of air passengers. We calculate monthly, country-specific dengue incidence rates by combining data from the Global Health Data Exchange [27], the most comprehensive health database, and known seasonal patterns in reported dengue infections [26]. Further, we distinguish between two categories of travellers: returning residents and visitors. The number of days people from these two categories spend in an endemic country, and therefore the risk of being infectious on arrival, vary greatly. The model predicts the number of imported dengue cases per month for any given airport and can be applied with relative ease to other vector-borne diseases of global concern, such as malaria, Zika or chikungunya.
Second, we apply the model to infer time-varying, region-specific reporting rates, defined as the ratio of reported to actual infections. Dengue reporting rates vary greatly across space and time, often by several orders of magnitude, and hence are difficult to determine [10]. The usual approach towards estimating country-specific reporting rates is to carry out cohort or capture-recapture studies that can be costly, are time consuming and may be biased [28]. Consequently, dengue reporting-rates remain unknown for most countries [10].
In this paper we focus on those countries that are most at risk of dengue introduction, i.e. non-endemic countries with vector presence. These countries will have the greatest benefit from our model as knowledge about the likely arrival times and places of infected people is crucial to prevent local outbreaks.
Materials and methods
CSIRO’s human research ethics committee CSSHREC has approved this study (approval number: Ethics Clearance 142/16). All data were analysed anonymously and individuals cannot be identified.
IATA data
The International Air Transportation Association (IATA) has approximately 280 airline members who together contribute to approximately 83% of all air traffic. Data is collected in form of travel routes, detailing the origin, destination and stopover airports. It contains over 10,000 airports in 227 different countries and dependencies. For each route the total number of passengers per month is given. We do not have any information on stopover times and whether passengers are leaving the airport during their stopover and therefore assume that all passengers continue their journey to the final destination instantly. Table A in S1 File lists the IATA 3-Letter Codes used to abbreviate airports in the main manuscript. As the recorded itineraries do not include any travel on chartered flights, we compare the IATA passenger volumes to official airport passenger statistics [29–47] to quantify the potential discrepancies between actual travel patterns and that reported by IATA. Table B in S1 File lists the countries where the difference in passenger numbers is greater than 15% (at country level) and countries where airport statistics were not available and the tourist data suggests inaccuracies in the IATA data (i.e. the number of tourists arriving in a particular country is larger than the total number of passengers arriving). We also excluded Singapore as a source of importation for Australia for the following reason: The Department of Home Affairs publishes Arrival Card data [48] that can be used to validate the IATA data. A comparison of the monthly travel volume from Singapore to Australia revealed that the IATA data overestimates travel volumes by approximately 112% on average in 2011 and 2015. This may be due to individuals who travel from other countries to Singapore and then directly continue to Australia and do not book their entire trip in one itinerary (this would be recorded as two separate trips in the IATA data that cannot be linked to each other). Due to this large discrepancy in the travel data we believe that our model will significantly overestimate the number of dengue infections imported from Singapore, and therefore exclude it as a source country for Australia.
The air transportation network
We begin by constructing twelve weighted, directed networks, using IATA data, to represent the monthly movement of air passengers during a given year. The networks are denoted , with m = 1,…,12 indicating the month of the year. The node set V comprises more than 10,000 airports recorded by IATA. To distinguish the travellers by their country of embarkation, we represent the edges of the network as ordered triples, (i, j, ωi,j(c, k)) ∈ E, where i, j ∈ V and ωi,j(c, k) is a function that outputs the number of passengers who initially embarked in country c with final destination airport k and travel from airport i to airport j as part of their journey.
Incidence rates and seasonal distributions
Calculating the number of infected passengers requires daily infection probabilities. We derive these from country-level yearly estimates of symptomatic dengue incidence rates that are published together with their 95% confidence intervals by the Global Health Data Exchange [27]. The estimates are obtained using the model published in [10] and account for under-reporting.
We first deduce monthly incidence rates using information on dengue seasonality published by the International Association for Medical Assistance to Travellers [26]. To do so we associate a weight with each month that indicates the intensity of transmission. To assign the weights we use a modified cosine function with altered period that matches the length of the peak-transmission season. The function is shifted and its amplitude adjusted so that its maximum occurs midway through the peak-season with value equal to the length of the peak-season divided by 2π. The months outside the peak-season receive a weight of one if dengue transmission occurs year around and a weight of zero if dengue transmission ceases outside the peak-season. The weights are then normalised and multiplied by the yearly incidence rate for the corresponding country. Normalising the weights ensures that the sum of the monthly incidence rates is equal to the yearly incidence rate. To calculate the lower and upper bounds of the monthly incidence rates, we multiply the normalised weights by the lower and upper bounds of the 95% confidence interval given for the yearly incidence rates.
The average probability, βc,m, of a person becoming infected on any given day during month m in country c is then given by
(1)
where γc,m is the monthly dengue incidence rate in country c during month m and dm is number of days in month m. Note that Eq (1) converts the daily incidence rate into the probability of a single person becoming infected with dengue on any given day during month m.
Inferring the number of infected passengers
Next, we present a mathematical model that approximates the number of dengue-infected people for each edge in the network . The time between being bitten by an infectious mosquito and the onset of symptoms is called the intrinsic incubation period (IIP). This period closely aligns with the latent period, after which dengue can be transmitted to mosquitoes [49]. The IIP lasts between 3 and 14 days (on average 5.5 days) and was shown to follow a gamma distribution of shape 53.8 and scale equal to 0.1 [50]. After completion of the IIP a person is infectious for approximately 2 to 10 days (on average 5 days) [50, 51]. The length of the infectious period was shown to follow a gamma distribution of shape 25 and scale equal to 0.2 [50]. We denote the sum of the IIP and the infectious period by n, which is rounded to the nearest integer after the summation. For travellers to import the infection from country c into a new location r they must have been infected with dengue within the last n − 1 days of their stay in country c. We now consider the following two cases: tc ≥ n − 1 and tc < n − 1, where tc is number of days spent in country c before arriving in region r. Since we do not know the exact date of arrival for travellers, we assume that arrival and departure dates fall within the same month and hence βc,m is the same for every day during the travel period.
If tc ≥ n − 1, that is the individual spent more time in country c than the sum of the lengths of the IIP and the infectious period, the probability of not being infected on return is equal to . The first term covers the possibility that the individual did not get infected whilst staying in country c and the second term covers the possibility that the individual got infected and recovered before arriving at a given airport (see Fig A in S1 File). Hence, the probability of a person, who arrives at a given airport from country c during month m, being infected with dengue is given by
(2)
If tc < n − 1, that is the individual spent less time in country c than the sum of the lengths of the IIP and the infectious period, the probability of not being infected on return is equal to , which covers the possibility that the individual did not get infected whilst staying in country c. Since tc < n − 1, the probability of recovery before arriving at a given airport is zero. Hence, the probability of a person, who arrives from country c at a given airport during month m, being infected with dengue is given by
(3)
We distinguish between two different types of travellers arriving at a given airport of region r: returning residents and visitors. We define a returning resident as a traveller who resides in region r and a visitor as a traveller who resides in country c and visits region r. Returning residents are expected to have stayed a couple of weeks in the endemic country, while visitors may have spent their whole life in the country.
Since we lack information on how long each individual spent in country c before arriving at an airport of region r, we substitute parameter tc by if the person is a returning resident,
being the average number of days a returning resident spends in country c before returning home. If the person is a visitor, parameter tc is substituted by
, the average number of days a visitor spends in country c before arriving at an airport of region r. We distinguish between returning residents and visitors since
.
We assume that the length of stay for returning residents follows a normal distribution with mean equal to 15 days and standard deviation of 2, i.e. . A previous study has shown that employees around the world are on average entitled to approximately 15 days of annual leave [52]. On the other hand, visitors likely spent all their lives in the endemic country. We assume that
, where μvis is equal to c’s median population age. Median population ages by country are published in the World Factbook by the Central Intelligence Agency [53].
For simplicity we do not take immunity to the different dengue strains into consideration.
Proportion of returning residents and visitors
Lastly, we need to infer the proportions of returning residents and visitors. As this information is not contained in the IATA itineraries, we use international tourism arrival data from the World Tourism Organisation [54]. The data contains the yearly number of international tourist arrivals by air for each destination country. From the IATA data we calculate the total number of arrivals per year for each country and hence can infer the ratio of visitors to returning residents. As we lack sufficient data, we assume that the ratio of visitors to residents is the same for each month.
Calculating the absolute number of infected passengers
Given the above, we can now determine the number of infected passengers Ik,m arriving at airport k during month m as follows:
(4)
where q is the proportion of residents inferred from the international tourism arrival data,
(5)
and
(6)
Evaluation of the model’s uncertainty
We performed a thousand runs of the model for each edge in the network, drawing the parameters from their respective distributions, to calculate the mean and standard deviation of dengue-infected passengers. In addition, we have conducted a global sensitivity analysis to identify the model parameters with the greatest influence. We used Sobol’s method [55] with 100,000 samples to carry out the sensitivity analysis. The parameter ranges are shown in Table 1. The analysis was done with SALib [56], an open-source Python library.
https://doi.org/10.1371/journal.pone.0225193.t001
Results
We run our model for two different years to explore the robustness of the proposed methodology. Specifically, the analysis is conducted for 2011 and 2015. The results for the year 2015 are presented in the main manuscript, while the results for 2011 are presented in the supplementary material. Fig 1 shows the number of predicted imported dengue infections per airport for August 2015, where the area of a node increases with the number of dengue cases imported through the corresponding airport. The map clearly shows that many non-endemic regions where the dengue-transmitting vectors Aedes aegypti or Aedes albopictus are present (coloured in light grey) have airports that are predicted to receive a high number of dengue infections. For a list of dengue endemic and non-endemic countries see Table C in S1 File. As resources for the control and prevention of dengue are often limited [57], these countries face a high risk of future endemicity.
The map shows the output of our model for August 2015.The area of a node increases with the number of dengue cases imported through the corresponding airport. Airports that are predicted to not receive any infections are not shown on the map. Endemic countries are coloured dark grey. Countries that are non-endemic and where dengue vectors Aedes aegypti and/or Aedes albopictus are present are coloured in light grey. The blue circles correspond to the top ten airports identified in Fig 2. The map was created with the Python GeoPandas package and publicly available shapefiles from Natural Earth (http://www.naturalearthdata.com/).
https://doi.org/10.1371/journal.pone.0225193.g001
In Fig 2 and Fig B in S1 File we plot the number of predicted dengue importations over time for the ten airports that receive the highest number of cases, lie in non-endemic regions with vector presence and where local cases have been reported in the past (more detailed plots with confidence intervals are shown in Fig C in S1 File). While the majority of airports listed in Fig 2 and Fig B in S1 File are predicted to receive between 50 and 150 cases each month, Miami International Airport (MIA) is estimated to receive between 146 and 309 cases each month during both years. With Orlando International Airport (MCO) and Fort Lauderdale–Hollywood International Airport (FLL) also represented amongst the airports with the highest number of imported cases, Florida faces a high risk of local dengue outbreaks. Los Angeles International Airport (LAX) is predicted to receive the second highest number of imported cases. In 2011 its monthly predictions vary between 97 and 205 cases and in 2015 between 113 and 253 cases. The remaining airports listed in Fig 2 and Fig B in S1 File are located in France, Germany, the Netherlands, Texas, and Queensland, Australia. A full ranking of all airports located in non-endemic countries with vector presence can be found in Table D in S1 File.
The number of predicted imported dengue infections for the top ten airports in non-endemic countries/states with vector presence for each month in 2015. A break in a line indicates that the corresponding airport was not amongst the top ten during the respective month. Airports are abbreviated using the corresponding IATA code. A full list of abbreviations can be found in the supplementary material (see Table A in S1 File).
https://doi.org/10.1371/journal.pone.0225193.g002
In addition to calculating the number of imported dengue infections per airport, the model further provides the number of infected passengers travelling between any two airports, thus revealing common importation routes. Table 2 and Table E in S1 File list the routes that carry the highest number of infected passengers whose final destinations lie in non-endemic countries with vector presence. Table F in S1 File lists the routes that carry the highest number of infected passengers whose final destinations lie in non-endemic countries irrespective of whether vectors are present. For example, the route between Denpasar and Perth is ranked third in 2011 in Table F in S1 File, but it is not considered in the ranking shown in Table E in S1 File, as there are no vectors in Perth. Fig D in S1 File shows a map of all importation routes into non-endemic countries with vector presence.
https://doi.org/10.1371/journal.pone.0225193.t002
In both years the highest predicted number of infected passengers are recorded during the northern hemisphere’s summer. The route between São Paulo International Airport (GRU) and Miami International Airport (MIA) is the exception, where the highest number of infected passengers is predicted during April. The routes with the highest estimated number of dengue-infected passengers terminate at airports in countries that are non-endemic and where dengue-transmitting vectors are present.
Returning residents and visitors
Next, we aggregate airports by country/state to predict the number of imported dengue infections on a coarser level. For non-endemic countries that cover an area larger than 5,000,000 km2 and where dengue vectors are present we aggregate airports by state. These countries are Russia, the United States of America and Australia. The comparison between passenger volumes recorded by IATA and official airport statistics indicated that the IATA data for Russia may be inaccurate, i.e. the difference in passenger numbers is larger than 15% (see Material and Methods). Hence, we did not perform a state-level analysis for this country. In Australia vectors are present only in Queensland [58]. While vectors have been observed in more than 40 different US states, autochthonous cases have been reported only in California, Florida, Hawaii and Texas [59].
Our model separately calculates the number of dengue-infected people amongst returning residents and visitors and hence we can identify which of these groups is more likely to import the disease into a given country or state. Fig 3 and Fig E in S1 File show the results for six non-endemic countries/states with vector presence that are predicted to receive the highest number of dengue importations each month. Results for the remaining countries and states are shown in Figs F–K in S1 File. We observe that the contributions of returning residents and visitors to the total number of imported dengue infections is predicted to vary greatly between the different countries and states. In Florida and Queensland returning residents are predicted to be the main source of dengue importation. In France and Italy approximately one third of all dengue infections are predicted to be imported by visitors while in Spain visitors import around 75% of all imported cases. For Switzerland we do not have any information about the ratio of returning residents to visitors. For the United States there is evidence in the form of surveillance reports that returning residents are indeed the main contributors to dengue importations [60]. For Queensland we predict that 95% and 94% of infections were imported by returning residents in 2011 and 2015, respectively. Our predictions are supported by Queensland’s dengue notification data (provided by Queensland Health), showing that 97% and 92% of all dengue importations in 2011 and 2015, respectively, were imported by returning residents.
Here we show the results for non-endemic countries/states with vector presence with the highest number of predicted imported dengue cases in 2015. The bars are stacked to distinguish between returning residents (green) and visitors (blue). The blue solid line corresponds to the total number of imported cases. The error bars correspond to the model’s coefficient of variation (see Material and methods). The six countries were selected because they are predicted to receive the highest number of dengue importations, are non-endemic and dengue vectors are established.
https://doi.org/10.1371/journal.pone.0225193.g003
Countries of acquisition
In addition to being able to distinguish between returning residents and visitors, the model also divides the imported cases according to their places of acquisition. Fig 4 and Fig L in S1 File show the model’s estimated percentage contribution of dengue importations by source country.
The predicted percentage contribution by source country and month in 2015. The size and colour of the circles indicate the percentage contribution of the corresponding country to the total number of imported cases. The y-labels indicate the yearly percentage contribution of the corresponding source country.
https://doi.org/10.1371/journal.pone.0225193.g004
Florida is predicted to import most infections from the Caribbean and Latin America, with infections acquired in Puerto Rico (PRI) predicted to peak during June and July and infections acquired in Brazil predicted to peak between January and April. We hypothesise that Florida receives such a high number of imported dengue cases due to its close proximity to the Caribbean, which has been endemic since the 1970s [61]. France is predicted to receive many infections from the Caribbean, in particular from Martinique and Guadeloupe which are French overseas regions and hence a high volume of air traffic from these regions to metropolitan France is expected. These predictions align with the fact that outbreaks of dengue in France coincide with outbreaks in the French West Indies, where most reported cases are acquired [62, 63]. In Italy the model predicts that the most common countries of acquisition are India and Brazil. India and Brazil are also the most common countries of acquisition for Switzerland in 2011. In 2015 Switzerland is predicted to receive most of their dengue importations from India and Thailand. Spain is predicted to import the majority of infections from Latin America and the Caribbean. For Queensland the model predicts that imported cases are acquired mostly in South-East Asia with Indonesia being the largest source. This is in agreement with previous studies [64] and the dengue case data that was provided by Queensland Health. In addition, we performed a rank-based validation of these results.
We obtained dengue case data from Queensland Health, which records the places of acquisition for each dengue case reported in Queensland. We rank the countries of acquisition by the total number of predicted and reported dengue-infected people who arrive in Queensland. We then plot the reported ranking against the predicted ranking. In addition, we plot the absolute number of reported importations against the absolute number of predicted importations and calculate Spearman’s rank correlation coefficient. Fig 5 and Fig M in S1 File show the results.
(A) Countries are ranked by the total number of predicted and reported imported dengue cases. The reported ranking is then plotted against the predicted ranking. Countries that were ranked by the model, but did not appear in the dataset receive a rank of i + 1, were i is the number of unique importation sources according to the dengue case data. Similarly, countries that appeared in the data and were not ranked by the model receive a rank of i + 1. For circles that lie on the x = y line (grey solid line) the predicted and reported rankings are equal. Circles that lie between the two dashed lines correspond to countries with a difference in ranking that is less than or equal to five. The circle areas are scaled proportionally to the number of reported cases that were imported from the corresponding country. Spearman’s rank correlation coefficient between the absolute numbers of reported and predicted importations is equal to 0.6. (B) The absolute number of reported dengue importations are plotted against the absolute number of predicted importations.
https://doi.org/10.1371/journal.pone.0225193.g005
The rank-based validation of our model demonstrates that overall, the model captures the different importation sources well. It does particularly well for the countries from which Queensland receives the most infections. Spearman’s rank correlation coefficient is equal to 0.6 for the year 2015 and equal to 0.58 for the year 2011. Below we explain some of the differences between the data and the model output.
For the rank-based validation the two largest outliers in both years are Fiji and Taiwan. The predicted ranking for Fiji in 2011 is 2, while the reported ranking is 10. In 2015 we estimate Fiji to be ranked fifth, however no cases were reported in 2015 and hence Fiji is ranked last amongst the reported cases. According to the Fijian government tourists are less likely to contract the disease than local residents as they tend to stay in areas that are not infested by Aedes aegypti mosquitoes [65] or where there is likely considerable control effort undertaken by tourism accommodation operators. Since the incidence rates incorporated into our model do not distinguish between different regions of a source country, the model is unable to account for such nuances. In 2011 and 2015 we estimate Taiwan to be ranked seventh and eighth, respectively, however no cases were reported in both years. This result is surprising as dengue occurs year-round in Taiwan [26] and approximately 44,000 and 16,000 Queensland residents travelled to Taiwan in 2011 and 2015, respectively.
Some of the differences between the observed percentages and the predicted percentages can be explained by under-reporting. It is possible that dengue awareness among travellers to one country is greater than the awareness amongst travellers to another country. Travellers with higher awareness levels are more likely to report to a doctor if feeling unwell after their return.
Country-specific reporting rates
The reporting rate of a disease is defined as the ratio of reported infections to actual infections. Dengue reporting rates vary greatly across space and time and are difficult to determine [10]. The usual approach to estimating country-specific reporting rates is to carry out cohort or capture-recapture studies that can be costly, are time consuming and may be biased [28].
We utilised our model to infer country- and state-specific reporting rates of imported cases by performing a least squares linear regression without intercept.
Table 3 and Table G in S1 File show the estimated yearly and seasonal reporting rates of imported cases for Queensland, Florida, France, Italy and Spain. To distinguish locally acquired and imported cases in Queensland, we use case-based data from Queensland Health where the country of acquisition is recorded. Travel-related dengue cases reported in Europe are published by the European Centre for Disease Prevention and Control (http://ghdx.healthdata.org/gbd-results-tool). Data for Florida is available from the Florida Department of Health (http://www.floridahealth.gov/diseases-and-conditions/mosquito-borne-diseases/surveillance.html).
https://doi.org/10.1371/journal.pone.0225193.t003
The results show that estimated reporting rates of imported cases are highest in Queensland, in particular during autumn. This is expected as dengue awareness campaigns are intensified between November and April [66]. In contrast, Florida has the lowest dengue reporting rate (1.3% in 2011 and 1.4% in 2015). This finding is supported by a previous study which found that awareness levels in Florida are extremely low [67]. The estimated reporting rates for the European countries are also low; however, the model predicts a substantial increase from 2011 to 2015. The question why reporting rates in Queensland are higher is challenging to answer, as we do not have any information about the true number of imported cases. However, Queensland has one of the best dengue prevention programs in the world. According to Queensland Health, other states and countries frequently ask for training and advice regarding surveillance and awareness campaigns.
Model uncertainty
We found that the average coefficient of variation of our importation model is 19.5% across both years. That is, the model’s standard deviation is on average equal to 19.5% of its mean. Fig N in S1 File shows the distribution of the coefficient of variation for several destinations.
The results from the global sensitivity analysis show that tc is the most important of the three model parameters with a total-order index of 0.94 (see Fig O in S1 File). The different values of the first-order and total-order indices indicate interaction between the model parameters. The second-order indices show that there is significant interaction between parameters tc and βc,m with a second-order index of 0.19, as well as between parameters tc and n with a second-order index of 0.1.
Since the range of parameter tc is large ([1, 29200] days), we performed the sensitivity analysis again for a shorter range of values ([1, 30] days) that is more realistic for returning residents who spend their holidays in an endemic country. In this case, parameter βc,m, with a total-order index of 0.6, is more important than tc, which has a total-order index of 0.35 (see Fig O in S1 File). The second-order indices show that there is still significant interaction between parameters tc and βc,m with a second-order index of 0.06, and between parameters tc and n with a second-order index of 0.07.
Discussion
To mitigate the risk of outbreaks from importation of dengue into non-endemic regions it is critical to predict the arrival time and location of infected individuals. We modelled the number of dengue infections arriving each month at any given airport, which enabled us to estimate the number of infections that are imported into different countries and states each month. In addition, the model determines the countries of acquisition and hence is able to uncover the routes along which dengue is most likely imported. Our results can also be used to estimate country- and state-specific reporting rates of imported cases.
Such knowledge can inform surveillance, education and risk mitigation campaigns to better target travellers along high risk importation routes at the most appropriate times. It will also help authorities to more efficiently surveil those airports with the highest risk of receiving dengue-infected passengers.
The model proposed here overcomes many of the shortcomings of previous models, however, it is not without limitations. Validation through comparison of reported cases to predicted cases is infeasible due to the high degree of under-reporting. However, we demonstrate that the coefficient of variation of the model with 19.5% on average is low (see Material and Methods). A rank-based validation for Queensland confirmed that the different importation sources are accurately predicted.
Incidence rates may vary considerably from region to region within the same country [65] and higher resolution data could improve the model’s predictions, as it would better reflect the export of dengue cases from the individual regions. Region-specific incidence rates can, for instance, be combined with spatial patterns of the visiting frequency of travellers to determine the likelihood of travellers to export dengue out of endemic countries. Additional data on individuals’ travel behaviour may also be beneficial, as it can be analysed to improve the estimation of the average time that a person has spent in a specific country before arriving at a given airport. Our assumption that returning residents and visitors are exposed to the same daily incidence rates is a simplification. Further details on the types of accommodation, for example, resorts vs local housing, could also be used to inform the daily incidence rates, due to variations in vector control. The global sensitivity analysis has revealed that tc, the number of days a traveller has spent in country c, is the most important model parameter. Hence, additional data on individuals’ travel behaviour may substantially improve the model. Knowledge about the exact age of visitors who reside in non-endemic countries would also improve the model. Currently, we assume that the age of a visitor is equal to the median age of the population of the country in which the visitor resides. In reality, the age of air passengers may differ from the median age, especially for developing countries.
In temperate regions local conditions may not allow for dengue to be transmitted during the winter months. Thus, even a large number of imported cases during those months would not trigger local outbreaks. Variable seasonality patterns due to El Niño Southern Oscillation can affect the spread of dengue in tropical and subtropical regions. An interesting direction for future research is to combine the here proposed model with knowledge of local conditions and weather phenomena like El Niño Southern Oscillation to evaluate the risk of local outbreaks. In this work we studied dengue importation via air travel. In future, we will also consider other modes of transportation to develop a more comprehensive model.
Supporting information
S1 File. Supplementary material.
This file contains all supplementary figures and tables.
https://doi.org/10.1371/journal.pone.0225193.s001
(PDF)
S2 File. Dengue incidence rates for 2011.
https://doi.org/10.1371/journal.pone.0225193.s002
(CSV)
S3 File. Dengue incidence rates for 2015.
https://doi.org/10.1371/journal.pone.0225193.s003
(CSV)
Acknowledgments
We would like to thank Frank de Hoog and Simon Dunstall for their constructive feedback which helped us to improve the model. We would also like to thank Queensland Health for providing dengue outbreak data. This work is part of the DiNeMo project.
References
- 1. Bloom DE, Black S, Rappuoli R. Emerging infectious diseases: A proactive approach. Proc Natl Acad Sci. 2017;114(16):4055–4059. pmid:28396438
- 2. Brockmann D, Helbing D. The hidden geometry of complex, network-driven contagion phenomena. Science. 2013;342(6164):1337–1342. pmid:24337289
- 3. Dorigatti I, Hamlet A, Aguas R, Cattarino L, Cori A, Donnelly CA, et al. International risk of yellow fever spread from the ongoing outbreak in Brazil, December 2016 to May 2017. Eurosurveillance. 2017;22(28):30572. pmid:28749337
- 4. Guimera R, Mossa S, Turtschi A, Amaral LAN. The worldwide air transportation network: Anomalous centrality, community structure, and cities’ global roles. Proc Natl Acad Sci. 2005;102(22):7794–7799.
- 5. Huang Z, Tatem AJ. Global malaria connectivity through air travel. Malar J. 2013;12(1):269. pmid:23914776
- 6. Hufnagel L, Brockmann D, Geisel T. Forecast and control of epidemics in a globalized world. Proc Natl Acad Sci. 2004;101(42):15124–15129. pmid:15477600
- 7. Tatem AJ, Huang Z, Das A, Qi Q, Roth J, Qiu Y. Air travel and vector-borne disease movement. Parasitology. 2012;139(14):1816–1830. pmid:22444826
- 8. Bogoch II, Brady OJ, Kraemer MUG, German M, Creatore MI, Kulkarni MA, et al. Anticipating the international spread of Zika virus from Brazil. Lancet. 2016;387(10016):335–336. pmid:26777915
- 9. Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013;496(7446):504–507. pmid:23563266
- 10. Stanaway JD, Shepard DS, Undurraga EA, Halasa YA, Coffeng LE, Brady OJ, et al. The global burden of dengue: An analysis from the Global Burden of Disease Study 2013. Lancet Infect Dis. 2016;16(6):712–723. pmid:26874619
- 11. Murray NEA, Quam M, Wilder-Smith A. Epidemiology of dengue: Past, present and future prospects. Clin Epidemiol. 2013;5:299–309. pmid:23990732
- 12.
World Health Organization. WHO heralds “new phase” in the fight against neglected tropical diseases. 2013. Available from: http://www.who.int/mediacentre/news/releases/2013/ntds_report_20130116/en/.
- 13. Gardner L, Sarkar S. A global airport-based risk model for the spread of dengue infection via the air transport network. PLOS One. 2013;8(8):e72129. pmid:24009672
- 14. Mackenzie JS, Gubler DJ, Petersen LR. Emerging flaviviruses: The spread and resurgence of Japanese encephalitis, West Nile and dengue viruses. Nat Med. 2004;10(12):S98–S109. pmid:15577938
- 15. Simmons CP, Farrar JJ, van Vinh Chau N, Wills B. Dengue. N Engl J Med. 2012;366(15):1423–1432. pmid:22494122
- 16. Chang F-S, Tseng Y-T, Hsu P-S, Chen C-D, Lian I-B, Chao D-Y. Re-assess Vector Indices Threshold as an Early Warning Tool for Predicting Dengue Epidemic in a Dengue Non-endemic Country. PLOS Negl Trop Dis. 2015;9(9):e0004043. pmid:26366874
- 17. Semenza JC, Sudre B, Miniota J, Rossi M, Hu W, Kossowsky D, et al. International dispersal of dengue through air travel: Importation risk for Europe. PLOS Negl Trop Dis. 2014;8(12):e3278. pmid:25474491
- 18. Gardner LM, Sarkar S. Risk of dengue spread from the Philippines through international air travel. Transp Res Rec J Transp Res Board. 2015;2501:25–30.
- 19. Wilder-Smith A, Quam M, Sessions O, Rocklöv J, Liu-Helmersson J, Franco L, et al. The 2012 dengue outbreak in Madeira: Exploring the origins. Eurosurveillance. 2014;19(8). pmid:24602277
- 20. Quam MB, Khan K, Sears J, Hu W, Rocklöv J, Wilder-Smith A. Estimating air travel–associated importations of dengue virus into Italy. J Travel Med. 2015;22(3):186–193.
- 21. Massad E, Amaku M, Coutinho FAB, Struchiner CJ, Burattini MN, Khan K, et al. Estimating the probability of dengue virus introduction and secondary autochthonous cases in Europe. Sci Rep. 2018;8(1):4629. pmid:29545610
- 22. Standish K, Kuan G, Avilés W, Balmaseda A, Harris E. High dengue case capture rate in four years of a cohort study in Nicaragua compared to national surveillance data. PLOS Negl Trop Dis. 2010;4(3):e633. pmid:20300515
- 23. Kakkar M. Dengue fever is massively under-reported in India, hampering our response. BMJ. 2012;345(17):e8574. pmid:23255584
- 24. Vong S, Goyet S, Ly S, Ngan C, Huy R, Duong V, et al. Under-recognition and reporting of dengue in Cambodia: A capture–recapture analysis of the National Dengue Surveillance System. Epidemiol Infect. 2012;140(03):491–499.
- 25. Wahyono TYM, Nealon J, Beucher S, Prayitno A, Moureau A, Nawawi S, et al. Indonesian dengue burden estimates: Review of evidence by an expert panel. Epidemiol Infect. 2017;145(11):2324–2329. pmid:28545598
- 26.
International Association for Medical Assistence to Travellers (IAMAT). 2018 [cited 6/11/18]. Available from: https://www.iamat.org/.
- 27.
Global Burden of Disease Collaborative Network. Global Burden of Disease Study 2016 (GBD 2016) Results. 2017 [cited 9/11/18]. Available from: http://ghdx.healthdata.org/gbd-results-tool.
- 28. Toan NT, Rossi S, Prisco G, Nante N, Viviani S. Dengue epidemiology in selected endemic countries: Factors influencing expansion factors as estimates of underreporting. Trop Med Int Heal. 2015;20(7):840–863.
- 29.
Airline Network News and Analysis. Aviation database. 2018 [cited 1/12/18]. Available from https://www.anna.aero/databases/.
- 30.
Civil Aviation Administration of China. 2015 Civil Aviation Airport Throughput Ranking. 2016 [cited 1/12/18]. Available from http://www.caac.gov.cn/XXGK/XXGK/TJSJ/201603/t20160331_30105.html.
- 31.
Ministry of Land Infrastructure Transport and Tourism. Ministry of Land, Infrastructure, Transport and Tourism,. 2016. Available from https://web.archive.org/web/20161021205147/http://www.mlit.go.jp/common/001141840.pdf.
- 32.
Ministry of Transportation Republic of Indonesia. Air transport traffic. 2018 [cited 1/12/18]. Available from http://hubud.dephub.go.id/?en/llu/index/filter:bulan,0.
- 33.
Airports Council International. North America airport traffic annual reports. 2015 [cited 1/12/18]. Available from https://airportscouncil.org/intelligence/north-american-airport-traffic-reports/north-america-airport-traffic-annual-reports/.
- 34.
Directorate General of Civil Aviation—Kuwait. Statistics. 2015 [cited 1/12/18]. Available from https://www.dgca.gov.kw/en/civil-aviation/media-and-info/statistics.
- 35.
EANA—Navegación Aérea Argentina. Statistical reports. 2016 [cited 1/12/18]. Available from https://www.eana.com.ar/estadisticas#node-150.
- 36.
Angkasa Pura II. Innovating beyond excellence—2015 annual report. 2015 [cited 1/12/18]. Available from https://cms.angkasapura2.co.id/NUWEB_PUBLIC_FILES/angkasapura2/Annual_25_07_2016__09_04_23.pdf.
- 37.
Civil Aviation Authority of the Philippines. Passenger movement. 2015 [cited 1/12/18]. Available from https://data.gov.ph/dataset/civil-aviation-authority-philippines-passenger-movement.
- 38.
Instituto de Estadísticas de Puerto Rico. Carga y passajeros aéreos y marítimos. 2015 [cited 1/12/18]. Available from http://www.estadisticas.gobierno.pr/iepr/Estadisticas/InventariodeEstad%C3%ADsticas/tabid/186/ctl/view_detail/mid/775/report_id/9485fbd8-efa5-4583-a8a4-a39c74eb846f/Default.aspx.
- 39.
Malaysia Airports. Connecting you seamlessly to the world—Annual report 2015. 2015 [cited 1/12/18]. Available from http://mahb.listedcompany.com/misc/ar/ar2015.pdf.
- 40.
Zambia: Transport Data Portal. Air transport statistics Zambia. 2016 [cited 1/12/18]. Available from http://zambiamtc.opendataforafrica.org/ZMATS2015/air-transport-statistics-zambia?airports=1000000-kenneth-kaunda-international-airport.
- 41.
California Department of Transportation. 2015:2014 air passenger and air cargo traffic activity report. 2016 [cited 1/12/18]. Available from http://www.dot.ca.gov/hq/planning/aeronaut/documents/statistics/15trafficgrowthrpt.pdf.
- 42.
Airports Authoriy of India. Passengers. 2016 [cited 1/12/18]. Available from https://www.aai.aero/sites/default/files/traffic-news/Mar2k15annex3.pdf.
- 43.
Associazione Italiana Gestori Aeroporti. Statistiche. 2016 [cited 1/12/18]. Available from http://www.assaeroporti.com/statistiche/.
- 44.
Office National des Aéroports. Passengers. 2017 [cited 1/12/18]. Available from http://www.onda.ma/en/I-am-a-Professional/Companies/statistics/(offset)/30.
- 45.
Airports Company South Africa. Statistcs. 2018 [cite 1/12/18]. Available from http://www.airports.co.za/news/statistics.
- 46.
Egyptian Holding Company for Airports and Air Navigation. Statistics. 2015 [cite 1/12/18]. Available from http://www.ehcaan.com/statistics.aspx.
- 47.
Airports Council International. Passenger traffic 2015 final (annual). 2018 [cited 1/12/18]. Available from https://aci.aero/data-centre/annual-traffic-data/passengers/2015-final-summary/.
- 48.
Department of Home Affairs. Overseas Arrivals and Departures. 2017 [cited 1/12/18]. Available from: https://data.gov.au/dataset/overseas-arrivals-and-departures.
- 49. Chan M, Johansson MA. The incubation periods of dengue viruses. PLOS One. 2012;7(11):e50972. pmid:23226436
- 50. Chowell G, Diaz-Dueñas P, Miller JC, Alcazar-Velazco A, Hyman JM, Fenimore PW, et al. Estimation of the reproduction number of dengue fever from spatial epidemic data. Math Biosci. 2007;208(2):571–589. pmid:17303188
- 51. Gubler DJ, Clark GG. Dengue/dengue hemorrhagic fever: The emergence of a global health problem. Emerg Infect Dis. 1995;1(2):55–57. pmid:8903160
- 52.
Messenger JC, Lee S, McCann D. Working time around the world: Trends in working hours, laws, and policies in a global comparative perspective. 1st ed. London: Routledge; 2007.
- 53.
Central Intelligence Agency. The World Factbook. 2017 [cited 12/11/18] Available from: https://www.cia.gov/library/publications/resources/the-world-factbook/fields/343.html.
- 54.
World Tourism Organisation. All Countries: Inbound Tourism: Arrivals by mode of transport 1995—2017. 2018 [cited 12/11/18] Available from: https://www.e-unwto.org/action/doSearch?ConceptID=2445&target=topic.
- 55. Sobol IM. Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates. Math Comput Simulat. 2001;55:271–280.
- 56. Herman J, Usher W. SALib: An open-source Python library for Sensitivity Analysis. The Journal of Open Source Software. 2017;2(9):97.
- 57. Morrison AC, Zielinski-Gutierrez E, Scott TW, Rosenberg R. Defining challenges and proposing solutions for control of the virus vector Aedes aegypti. PLOS Med. 2008;5(3):e68. pmid:18351798
- 58. Beebe NW, Cooper RD, Mottram P, Sweeney AW. Australia’s dengue risk driven by human adaptation to climate change. PLOS Negl Trop Dis. 2009;3(5):e429.
- 59. Hahn MB, Eisen RJ, Eisen L, Boegler KA, Moore CG, McAllister J, et al. Reported Distribution of Aedes (Stegomyia) aegypti and Aedes (Stegomyia) albopictus in the United States, 1995-2016 (Diptera: Culicidae). J Med Entomol. 2016;53(5):1169–1175.
- 60. Van Dodewaard CAM, Richards SL. Trends in dengue cases imported into the United States from Pan America 2001–2012. Environ Health Insights. 2015;9:33–40.
- 61. Amarakoon Dharmaratne, Chen A, Rawlins S, Chadee DD, Taylor M, et al. Dengue epidemics in the Caribbean-temperature indices to gauge the potential for onset of dengue. Mitig Adapt Strateg Glob Chang. 2008;13(4):341–357.
- 62. Vasquez V, Haddad E, Perignon A, Jaureguiberry S, Brichler S, Leparc-Goffart I, et al. Dengue, chikungunya, and Zika virus infections imported to Paris between 2009 and 2016: Characteristics and correlation with outbreaks in the French overseas territories of Guadeloupe and Martinique. Int J Infect Dis. 2018;72:34–39. pmid:29782922
- 63. La Ruche G, Dejour-Salamanca D, Bernillon P, Leparc-Goffart I, Ledrans M, Armengaud A, et al. Capture–recapture method for estimating annual incidence of imported dengue, France, 2007–2010. Emerg Infect Dis. 2013;19(11):1740–1748.
- 64. Warrilow D, Northill JA, Pyke AT. Sources of Dengue Viruses Imported into Queensland, Australia, 2002-2010. Emerg Infect Dis. 2012;18(11):1850–1857. pmid:23092682
- 65.
The Fijian Government. Minimal risk of dengue fever in major tourism areas of Fiji. 2014 [cited 1/12/18]. Available from: http://www.fiji.gov.fj/Media-Center/Press-Releases/MINIMAL-RISK-OF-DENGUE-FEVER-IN-MAJOR-TOURISM-AREA.aspx.
- 66.
Queensland Health. Queensland dengue management plan 2015-2020. 2015. Available from: https://www.health.qld.gov.au/__data/assets/pdf_file/0022/444433/dengue-mgt-plan.pdf.
- 67. Hayden MH, Cavanaugh JL, Tittel C, Butterworth M, Haenchen S, Dickinson K, et al. Post outbreak review: Dengue preparedness and response in Key West, Florida. Am J Trop Med Hyg. 2015;93(2):397–400. pmid:26078319