ABSTRACT
Although the United States (U.S.) declared measles eliminated in 2000, so long as measles circulates globally, continued importations will trigger sporadic U.S. outbreaks. The United States is highly connected to the rest of the world, and importations occur largely as a result of air travel by U.S. residents and international visitors. We developed a model to assess the risk of measles virus importation from Europe, Asia, and Africa through air travel. We projected 308 (95% prediction interval, 151-518) total measles cases imported from 66 countries during 2006-2015; 290 importations were reported. The model projected a spike in importation risk from Israel during October 2018, coincident with large, importation-triggered outbreaks in New York and New Jersey. Our model shows U.S. importation risk primarily originates from European and Asian countries. Models such as this inform preemptive actions to mitigate infectious disease threats, but only if we prioritize global surveillance and data sharing.
INTRODUCTION
Approximately 78 million tourists1 and 44 million residents2 enter the United States (U.S.) from abroad via commercial flights every year – equivalent to over 300,000 travelers daily. With each traveler comes the potential for an infectious pathogen to be brought into the United States. The rapid, long-distance transmission potential posed by global air travel has become more evident during the last two decades, particularly in the context of emerging pathogens. In 2002, SARS coronavirus emerged in China. It spread rapidly to 36 countries, transmitted largely by air travelers, causing 8,098 cases and 774 deaths worldwide3. Six years later, the 2009 H1N1 influenza pandemic spread quickly across the world from Mexico; it caused 284,000 deaths worldwide and stoked widespread fear of a repeat 1918 influenza pandemic4. Recent deadly outbreaks of Middle East respiratory syndrome (MERS) and Ebola have prompted travel advisories, airport screening, and travel restrictions for exposed individuals5–10.
While these outbreaks highlight concerns about the spread of emerging pathogens, global travel is similarly critical to the spread and reintroduction of infectious diseases that countries have struggled to eliminate, particularly measles. In 2000, the United States declared measles elimination11. However, since 2000, the United States has reported more than 100 measles outbreaks, over 2,000 reported cases, and 3 measles-related deaths, all arising from infections in travelers12–14. With increasing globalization, the threat posed by geographically distant outbreaks will continue to grow. Achieving elimination of a highly contagious disease like measles requires immense efforts and resources, and because of the constant threat of reintroduction, maintaining elimination can be challenging and expensive.
Having a model to help predict when and where measles importation may occur, therefore, can assist in planning for resource allocation. Relying on global measles incidence data and 2006-2015 airline itinerary information for international arrivals into the United States from other countries, we quantified the potential for measles importation using a framework previously developed to model arbovirus spread15,16. For simplicity, we broadly define this potential as the projected “risk of importation” and quantify it with crude estimates of the absolute number of importations per time period and the ratio of these importations between locations and time periods. Because a multitude of factors can modify how this potential for importation translates into actual imported cases, we focus on this projected potential and not the absolute estimates. We validated our approach with data on reported measles importations. This model provides the risk of measles importation, both where importations are most likely to occur (i.e., into which U.S. cities) and where they are most likely to originate (i.e., from which countries). This work provides a platform that public health officials can use for real-time risk estimates and demonstrates the value of global disease surveillance efforts.
RESULTS
We estimated the monthly risk of measles importation into the United States from 2006 through 2015 from 66 countries in Europe, Africa, and Asia for which data were available on measles incidence and air travel to the United States. For some countries measles surveillance data were not available for all periods. Hence, we restricted our analysis for each country to months when data were available. The final data set comprised 7,260 country-months (92% of those possible over the study period). During this time, the U.S. Centers for Disease Control and Prevention (CDC) reported 412 confirmed imported measles cases. Of these, 376 (91%) originated in one of the 66 countries included in our analysis, and 290 (70% of all importations) were reported during months when measles surveillance data were available from origin countries. We estimated importation risk for all 7,260 country-months with surveillance data and compared model estimates to the 290 imported cases reported during those months.
Between 2006-2015, the 66 included countries reported 1,320,171 confirmed cases of measles, and travelers made an estimated 340 million trips to the United States from these countries. Global incidence of measles was characterized by annual periodicity, with several distinct outbreaks in Europe, Asia, and Africa (Figure 1). In Africa, large outbreaks occurred in Malawi (134,000 cases; 2010), Nigeria (58,000 cases; 2013), Burkina Faso (54,000 cases; 2009), and Zambia (38,000 cases; 2010-11)17–20. In Asia, a large measles outbreak occurred in the Philippines (59,000 cases; 2014), and both India and China reported substantial average annual case counts (34,000 and 47,000, respectively)21. In Europe, notable outbreaks occurred in Bulgaria (23,500 cases; 2009-10), France (16,000 cases; 2010-11), and Ukraine (2005-06: 42,000 cases; 2011-12: 14,000 cases)21–23. Overall, most travel to the United States originated from Europe (68%) and Asia (30%), with only 2% from Africa. In the most recent years, 2010-2015, travel volume increased consistently, by 4%, 6%, and 10% per year from Europe, Asia, and Africa, respectively (Figure 2).
Reported measles cases and importations
Reported monthly measles case counts in source countries correlated moderately with reported monthly measles importations into the United States (r=0.52; Figure 1). Correlation was highest for Asia (r=0.67), where 802,535 reported measles cases resulted in 135 reported importations (Figure 1B, Table1). The number of importations from Europe was similar (n=136), though reported cases (162,779) was lower, as was correlation (r=0.39). Despite 362,182 total cases reported from twelve African countries, only 19 importations into the United States were reported, and correlation between reported cases and importations from Africa was low (r=0.23).
Estimated risk of importation into the United States, by source region
To examine contributions by region and country, we compared mean annual estimates of importation and relative risk. Mean annual estimates were used rather than cumulative estimates as data were not available during part of the study period for several high-risk countries (e.g., data from China were not available for 2006-2010, nor were data available from India for 2006-2013, during which 10 and 50 importations occurred respectively). On average, estimated measles importations predominantly originated from Asia and Europe (60% and 35%, respectively), with only 5% reported from Africa annually. Taken alone, neither case nor travel data fully captures of the risk of measles importation: among the three regions, Asia and Africa reported the majority of measles cases (61% and 27%), with only 12% occurring in Europe. In contrast, 68% of travel volume to the United States originated from Europe, compared to 30% from Asia and only 2% from Africa. Using this model and accounts for reported cases, travel volume, and source country population, we estimated an average of 22.9 (95% CI, 1.8-76) measles cases imported from Europe annually and 19.7 (95% CI, 6.5-64) from Asia (Table 1), with Africa contributing just over one measles importation each year (mean=1.3; 95% CI, 0.1-5.1). These are comparable to average annual reported importations of 17.2 (95% CI, 0.5-68) for Europe, 26.1 (95% CI, 8.4-68) for Asia, and 2.1 (95% CI, 0-8.8) for Africa.
Estimated risk of importation into the United States, by country
The five countries with the highest average annual estimated risk of measles case exportation to the United States were India, the United Kingdom (U.K.), France, the Philippines, and China (Table 2, Figure 3). Risk was greatest for countries with high travel volume to the United States and either persistent measles cases (e.g., the U.K., India, China) or large outbreaks (e.g., France, the Philippines). Yearly country-specific mean importation risk correlated moderately with mean travel volume (r=0.42) and reported cases (r=0.43), and poorly with population size (r=0.20).
Outbreaks and associated importation risk
The periods of greatest estimated risk during 2006-2015 occurred during the first halves of 2008, 2011, and 2014, corresponding to major outbreaks in Japan and Switzerland (2008), France (2011), and the Philippines (2014) (Figure 3–4). Major outbreaks do not always constitute high risk: while the 2011 outbreak in France (16,350 cases reported) produced 17 reported importations, the larger 2010 outbreak in Bulgaria (23,500 cases reported) produced none (Table S6). While Bulgaria had more cases, annual traveler volume from Bulgaria to the United States was only 3% of that from France. Accounting for these differences in cases, travel, and population, we estimated the importation risk from the French outbreak was 5.2 (95% CI, 1.8-15.2) times as high as from the Bulgarian outbreak, with 23 (95% CI, 2033) versus 6 (95% CI, 4-11) expected importations into the United States, respectively.
Seasonality and timing of importation risk
Global measles cases (Figure 1) and travel (Figure 2) fluctuate seasonally, driving seasonality in importation risk into the United States, with elevated risk typically occurring in April–July (Figures 3, 5). Measles incidence typically peaks between March and May for Europe and Asia but is more irregular in Africa (Figure S2). In contrast, travel to the United States from all regions peaks in July and August (Figure S3). The temporal alignment of incidence and travel is critical to importation risk. Were measles cases in Europe to peak three months later (in July), when travel between the United States and Europe is at its height, we estimate the number of importations from Europe would increase by approximately 15% (95% CI, 8-32%). For Asia and Africa, this shift in measles timing would increase importations by 4% (95% CI, 3-5%) and 11% (95% CI, 2-12%), respectively.
Estimated risk of importation, by state
Our model predicted that New York/New Jersey/Connecticut (NY/NJ/CT, combined because of shared airports) and California had the highest risks of importation among U.S. states, experiencing 10.4 (95% CI, 8.5-12.3) and 8.4 (95% CI, 6.7-10.2) times the mean risk across all states during the study period (Figure 5A-B, Table 3). We estimated there were 7.8 (95% CI, 0-40) and 6.2 (95% CI, 0-35) annual importations on average, into these locales, respectively. Risk correlates highly with mean annual state travel volume and population size (Pearson’s r=0.88, r=0.89). Several states remained at significantly elevated risk after accounting for population size (NY/NJ/CT, California, Hawaii) and travel volume (California, Hawaii; Table 3). Hawaii and West Coast states experienced greater importation risk from Asian countries, while states on the East Coast were more at risk from Europe (Figure 5C). In Florida, 82% of estimated annual importations were associated with travelers originating from European countries, contrasting with only 40% for Washington. For Hawaii, in one year alone (2008), Japan contributed 53% (95% CI, 28-100%) of the risk during the ten-year study period.
Model Accuracy and Validation
The model’s 95% prediction interval estimates successfully captured the reported number of measles cases imported into the United States for 99% percent of the 7,260 country-months under study (Figure 6). We also assessed its ability to predict the occurrence of at least one importation, which we defined as the probability of at least one imported case being greater than 50%. Sensitivity of this metric was low for both source countries (18%) and U.S. states (18%) and specificity was high (99% for both), indicating that many low-probability importations occur. The positive predictive value was 40% for sources and 49% for states indicating that the 50% probability threshold is a well-calibrated and reliable indicator of risk. Negative predictive value was very high for both (greater than 97%). Model estimates of imported cases correlated highly with numbers of reported imported cases, for both annual and mean country-specific imported cases (ρ =0.71 and ρ =0.83), while correlating poorly or moderately with travel volume, measles incidence, and population (ρ=0.31-0.53, Figure S5-S6). We also compared the full model (with both travel and incidence data) to travel-only or incidence-only models, finding that the full model was significantly more likely to predict the reported importations (likelihood ratio [LR]=44 compared to incidence-only and LR=22 compared to travel-only; Table S5). It should be noted that reported importation cases were not used in any way to fit or calibrate the model, and this performance was obtained completely from reasoning from first principles.
Forecasting measles importations
We conducted a secondary analysis predicting measles importations using forecasted travel data for 20162019. Predicted importations were not validated against reported importations, as full details on imported cases were not yet available for this time period. The greatest importation risk occurred in NY/NJ/CT, California, and Florida, associated with travelers originating from both Asia and Europe. The periods of greatest importation risk occurred during March-May 2018, with the Philippines, Italy, France, India, UK, Greece, and Ukraine contributing the most risk, and October 2018-February 2019, from Israel, Ukraine, and Philippines (Figure 7). During this latter period, we estimated that 14.3 (95% CI, 7-22) importations into the United States occurred from Israel, with 5.8 (95% CI, 2-11) of these into NY/NJ/CT. This is a significant divergence from previous risk from Israel, from which there were only 5 reported and 7 (95% CI, 2-12) estimated importations for 2006-2015.
DISCUSSION
Two decades after declaring measles elimination, the United States continues to experience measles outbreaks, reporting a record 1,172 cases during January-August 201911,24–26. This measles incidence is not only a public health challenge for the U.S, but a global problem, with numerous countries at or close to elimination recently experiencing significant outbreaks. Increasing measles incidence and outbreaks can be attributed to insufficient vaccination coverage, at least partly caused by vaccination hesitancy and apathy27. With increasing global incidence, the potential for reintroduction of measles virus into countries where the disease had been eliminated, particularly through air travel, poses a significant threat. During 2000-2016, a quarter of measles cases were directly imported into the United States, and the rest resulted from those importations28.
As in the United States, recent outbreaks in the Republic of Korea29, Israel30, Japan31, New Zealand32, and other countries demonstrate importation of measles virus as a major challenge for control. While herd immunity limits endemic transmission and large-scale outbreaks in the United States, without complete immunity (United States was estimated in 2010 to have 92% overall immunity), potential for transmission from importations remains. This becomes particularly problematic when infected travelers live or interact within clusters of susceptibility, producing even greater potential to spark local transmission33–35. To mitigate the risk posed by frequent reintroductions, it becomes critical that we understand when and from where infected travelers will arrive, and where they will go, enabling public health agencies to potentially limit their numbers through public awareness and vaccination and limit transmission through prioritized resources and rapid detection and response.
Previous work has used case and travel data to extrapolate the risk of measles importation. Work by Sarkar et al. recently highlighted the potential to assign risk based on source case and travel data, and how importation can compound outbreak risk produced by non-medical exemptions36. Massad et al. estimated the number of imported measles cases into Saudi Arabia during Hajj, emphasizing the risk of transmission and global dissemination following the annual pilgrimage37. While these previous studies have demonstrated the value of accounting for both travel and case data, our study is the first to estimate monthly measles importation risk and numbers from international sources into the United States and validate them with detailed data on reported importations collected by the CDC. As we demonstrate, because measles incidence and travel fluctuate over time, correctly accounting for their synchrony through highly time-resolved estimation is critical to appropriately attribute risk. This time-resolution, coupled with measles-specific clinical characteristics (i.e., incubation time, time to characteristic rash), enabled us to accurately project the monthly number of measles importations from specific source countries to specific states in the United States
As expected, we found that importation risk typically concentrates in states/districts with large populations or those with airports that serve as common ports of entry (e.g., New York, California, District of Columbia, and Massachusetts), though some exceptions were found (e.g., Hawaii). We found that two factors (cases of disease in the country of origin and air travel volume) combined to drive importation risk. Periods of elevated importation risk were generally the consequence of large-scale outbreaks (e.g., in France or the Philippines). However, the minimal importation risk produced by large outbreaks in Bulgaria and Malawi demonstrate the necessity of sufficient travel. Overall, both Asia and Europe contributed similar overall importation risk to the United States, with Asia’s greater measles cases countering Europe’s travel volume. Despite massive measles outbreaks and high incidence, Africa contributes only minimal estimated and reported risk, highlighting the critical importance of accounting not only for cases, but also connectivity.
Real-time application of this model demonstrates enormous potential for identifying and responding to emerging risk (Figure 7). In October 2018, eight imported cases from Israel were identified in New York and New Jersey, resulting in outbreaks totaling at least 275 cases largely concentrated among the orthodox Jewish community38. Historically, Israel has presented a minimal importation risk to the United States, with only 5 total importations reported from 2006-2015 (we estimated 7). Using projected travel data and real-time surveillance data, our model estimated a dramatic increase in risk, peaking in October 2018, coincident with an outbreak of 4,250 reported in Israel. Of interest, the Israeli outbreak resulted from international travelers importing measles to that country30.
The model framework presented here is a parsimonious representation of a much more complex system, and discrepancies between reported and estimated importations partially reflect the model’s simplicity. We do not account for individual-level characteristics, such as vaccination status, or the correlation between travel and infection probabilities—both likely invalid in many populations. Age, for example, is typically associated with infection risk (through susceptibility and contact) and individual likelihood of traveling (e.g. travel purpose, destination)39,40. Traveler’s countries of residence, which was not available at the granularity needed and thus not accounted for, may impact the probability of infection: U.S. residents constituted 40% of international travel to the United States, but 62% of imported measles cases (during 2001-2016)28. However, the travel data do account for the actual origin country and U.S. destination state/state group, giving confidence to our attribution of risk to locations. Furthermore, while we assume people do not travel once rash has begun, it is possible that some infected individuals do travel with rash or develop rash during travel.
While random individual variability should have minimal impact on accuracy, population-wide characteristics can have a large impact. We saw this during the 2010 Bulgarian outbreak where 90% of cases were among the Roma population and, despite 6 estimated importations, none were reported22,41. Similarly, the 2013 Dutch outbreak, which was concentrated (>84%) among spatially-isolated orthodox Protestants, produced only 2 reported importations compared with 7 estimated42. Both populations have well-documented histories of insufficient vaccination, are typically socially and spatially isolated, and may be less likely to travel internationally or interact with travelers41–43.
This model is also sensitive to data limitations, particularly with reported measles cases (often underestimated and underreported44–47) and reported importation cases (required for validation). For an imported case to be reported, the person must either (1) seek medical care, be diagnosed, then be reported by the provider, or (2) be linked to secondary cases who sought care. Even with highly characteristic measles symptomology and essentially no asymptomatic infections, historically, measles reporting has been incomplete, even in the United States46–49. Moreover, foreign visitors might be less likely to seek care, particularly if they depart from the United States early during the infectious period. Importation estimates during a large measles outbreak in Japan during 2007-2008 highlight the challenges of using reported data. While an estimated 4,000-18,000 cases occurred in 200750,51, World Health Organization (WHO) data on reported measles cases include only five cases. This significant underreporting of cases in the country of origin (i.e., Japan) led our model to estimate no importations when ten were reported. In contrast, following national policy changes to require reporting of all measles cases, Japan reported more than 17,000 cases in 2008, from which our model estimated 16 importations. Interestingly, no importations into the United States from Japan were reported during 2008, potentially resulting from changes in vaccination policies, travel advisories, behavior changes, or underreporting of importations50. Similarly for Pakistan, reported cases were consistently lower than those found in published sources, contributing to fewer estimated importations than reported52.
Despite clear limitations, our simple model does surprisingly well at estimating the number of cases imported into the United States based solely on reasoning from first principles. Similar models simulating yellow fever epidemics15, chikungunya virus spread in the Americas16, and dengue importation in China53 have also demonstrated high accuracy, indicating that this model framework emulates the key processes of disease importation. This accuracy also suggests that some of the more difficult to measure components (e.g., reporting rates, traveler characteristics, population mixing) likely are balanced across most of these countries, most of the time. Nonetheless, a model like this is unlikely to predict rare, but impactful importation events such as that which caused the 2014 outbreak in an Ohio Amish community or the 2017 outbreak in the Minnesota Somali community34,54.
For the United States and other governments concerned with importation of diseases like measles, this work highlights two important priorities. First, reducing the risk contributed by resident travelers is essential, accomplished through increased risk awareness55 and vaccination coverage56 among travelers (particularly high-risk travelers and those from communities with heterogeneous vaccination coverage57). Second, this work stresses the tremendous value of information to both domestic biosecurity and global disease control. Our world is highly connected, and infectious diseases are global problems that require global cooperation and efforts, including global disease surveillance systems. The better global systems are at collecting and sharing data, the more timely and accurate those data will be, and the more likely it is that they will be useful in preparedness efforts to prevent outbreaks and protect people. Here we focus on measles importations into the United States, but our methods can readily be applied to other countries and other diseases. Virulent pathogens, like arboviruses, zoonotic coronaviruses, Nipah virus, and highly antibiotic-resistant bacteria, continually emerge around the world and, with the right conditions, can rapidly spread across local and global populations. When this occurs, we need accurate case information and effective tools to guide efforts to prevent and control these disease threats.
CONCLUSIONS
Until we achieve global eradication, measles virus importations will continue. These importations present a continual risk to spark measles outbreaks in the United States and elsewhere, even after national elimination. With effective tools, like this model, and timely surveillance and travel data, we can identify high-risk destination locations and high-risk source countries in real-time, potentially enabling us to mitigate future biosecurity threats from infectious disease importations through preemptive actions. Increasing and maintaining high routine vaccination coverage, in both the United States and abroad, is the only foolproof means of minimizing the number and the subsequent impact of measles importations.
METHODS
Model
Using monthly measles case counts and air travel volume from 66 countries in Europe, Asia, and Africa, we estimated the expected numbers of imported measles cases and the relative risk of measles importation into forty-two U.S. destination states or regions. Our approach was based on previous models of importation for arboviruses.15,16 We assumed that importation primarily occurs through air travel, and that it is directly proportional to the number of travelers and the number of cases of infection in the source population. We estimated the monthly number of infected travelers arriving at destination d as a binomial process dependent on the number of infected individuals (Is,m) in source country s during month m, the duration of infection during which an individual could travel (D), and the probability that an individual travels from s to d and is detected in d during month m (ps,d,m). This is estimated as follows:
Using the monthly number of travelers from s to d (Ts,d,m), the total population size of s (Ns,m), and the number of days per month (δm), we approximated ps,d,m:
We performed 10,000 stochastic simulations, sampling D from a distribution for each simulation, and calculating ns,d,m for each s, d, and m combination.
Data
Source country measles case and population data
Monthly reported numbers of measles cases are publicly available from the World Health Organization and the European Centers for Disease Control. These data were checked against country and regional reports where possible to verify accuracy. Annual case counts from monthly data were checked against reported annual measles cases available from the WHO21. Where annual totals differed by more than 25%, monthly data were considered invalid and excluded. Monthly case data were missing for multiple years for several countries. Where data were deemed invalid or missing, these country/year combinations were excluded from the analysis (e.g., India 2006-2013). Population data are publicly available from the United Nations World Population Prospect data58.
Travel data
Summary travel data are proprietary data from Diio and OAG in two batches59, 60. First, data were collected for 28 European countries for 2006-2015 from Diio in October 2016. Data for the additional 38 countries – 21 Asian, 12 African, and 5 European countries (Belarus, Bosnia-Herzegovina, Croatia, Switzerland, Ukraine) for 2010-2015 were collected from OAG in July 2017. These data are modeled estimates of trip numbers based on ticket sales and reporting from carriers and include locations of initial departure and final destination. As international airports like O’Hare International Airport in Chicago or John F. Kennedy International Airport in New York also serve nearby states, we grouped several of these states together to limit bias, though this likely does not fully resolve this bias. To provide model results for all 66 countries for a ten-year span from January 2006 – December 2015, we imputed missing travel data from OAG using the available data and regionally-aggregated monthly travel data from the U.S. National Travel and Tourism Office61 (Table S8). The imputation model was validated using the 28-country European travel data set, using leave-one-out methods (see Supplement for more details).
Measles importations to the United States
To validate the model, we use data on the reported number of cases of measles virus infection imported into the United States, which were provided by the U.S. Centers for Disease Control and Prevention. These data include the month and state of detection, and country of origin or travel. Reported importations that originated from countries other than the 66 included or during time periods where travel or measles case data were not available were excluded to ensure comparability with model results (n=112). Among the excluded reported importations, 50 originated in India, occurring during 2006-2013, during which monthly measles case data are not available for India. Because this is such a large portion of the reported importations (12%), we opted to examine the annual means of importation and risk, rather than the ten-year totals for the full study period, whereby annual means were estimated across years for which data were complete. Because India reported relatively consistent annual measles case counts (median=32,546; IQR=27,440-43,480), the estimated mean annual cases and importations for India should not change substantially were data for 2006-2013 included. We also should note that the numbers of importations that we report do not fully concur with those previously published28, as individuals reporting multiple potential countries of infection (n=22) were counted as an importation from each reported country.
Disease Parameters
The duration during which an infected individual could travel post-infection (D). For measles, we assumed individuals do not travel once the rash develops. Based on previous estimates, we quantify D as the sum of a normally-distributed incubation period (mean=12.5 days, 95% CI=11.8- 13.3) and a uniformly-distributed time to rash (2-3 days)62.
Primary analysis data were restricted to January 2006 through December 2015, the period for which data were complete for all sources.
Forecasting Analysis
To demonstrate our ability to forecast the risk of measles importation, especially in real-time, we repeated the primary analysis for January 2016 through February 2019, using a combination of reported measles case data and forecasted travel data for all 66 countries. These forecasted travel data were generated using the same model constructed to impute missing travel volume data.
Model validation
Model performance was evaluated using Pearson’s correlation between the estimated importations and risk, and reported importation cases, by year and by month. We further validated the model through comparison of fits to reported importation data with those of probability models that excluded data (i.e., excluded travel data, case data, and both; see Supplement).
Relative risk of importation
Using the simulation results for the absolute number of importations, we also estimated the relative risks of importation. There are a multitude of factors that can directly influence the absolute magnitude of importations, including individual-level characteristics of travelers and the accuracy of reported measles cases. While we cannot directly account for all these factors, assuming they are relatively consistent across destinations and sources, we can use the relative risk (RR) of importation as a metric to compare risk between destinations and sources. We estimate the RR of importation as the ratio of the estimated number of imported cases to the mean number of importations for that simulation across all time periods and sources or destinations as the reference. Thus, for each simulation, a single reference value is established across all countries and time periods. The mean RR and confidence intervals across the 10,000 simulations are estimated for each destination/time or source/time combination.
All analyses were performed using R 3.5.2 and RStudio 1.1.463.
Data Availability
These analyses require three primary sources of data: reported measles cases from source locations, travel volume data, and population data. Reported monthly measles cases are available publicly from the World Health Organization (https://www.who.int/immunization/monitoring_surveillance/burden/vpd/surveillance_type/active/measles_monthlydata/en/). Population data are publicly available from the United Nations World Population Prospects data (https://population.un.org/wpp/). Monthly air travel are primarily from OAG (OAG Aviation Worldwide, https://www.oag.com/). These data were provided under contract and cannot be provided by the authors. Additionally, data on reported measles importation were used for validation. These data were provided by the U.S. Centers for Disease Control and Prevention and are not currently available to the public.
Acknowledgments
This work was supported by the Bill & Melinda Gated Foundation [S.A.T. and J.L., grant OPP1094793]. The air passenger data used in this study were purchased from Diio, LLC and OAG Aviation Worldwide Ltd. These data were used under license for the current study and so are not publicly available. The authors are available to share the air passenger data upon reasonable request from OAG and with the permission of OAG Aviation Worldwide Ltd. The authors no longer have access to Diio, LLC data. For researchers, the opportunity to access the data would be through the procurement of Diio Mi FMg. Please contact Diio’s Senior Vice President, Jordan Kayloe, at jordan.kayloe{at}diio.net.
The opinions expressed by authors contributing to this journal do not necessarily reflect the opinions of the Centers for Disease Control and Prevention or the institutions with which the authors are affiliated.Author information
Affiliations:
Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, 21205, USA
Shaun A. Truelove, Justin Lessler
Biomedical Advanced Research and Development Authority, U.S. Department of Health & Human Services, Washington, DC, 20201, USA
Luis Mier-y-Teran-Romero
Division of Vector-Borne Diseases, Centers for Disease Control and Prevention, Atlanta, GA, 30333, USA
Allison Taylor Walker, Andre Berro
Division of Global Migration and Quarantine, Centers for Disease Control and Prevention, Atlanta, GA, 30333, USA
Paul Gastanaduy
Division of Viral Diseases, Centers for Disease Control and Prevention, San Juan, Puerto Rico, 00920, USA
Michael A. Johansson
Center for Communicable Disease Dynamics, Harvard TH Chan School of Public Health, Boston, MA, 02115, USA
Michael A. Johansson
Contributions:
S.A.T., L.M.T.R, J.L., and M.A.J. conceived and designed the project. S.A.T. drafted the manuscript.
S.A.T. and L.M.T.R conducted the analyses. S.A.T., L.M.T.R, A.T.W., P.G., A.B., and M.A.J. contributed to data acquisition. All authors discussed the results, contributed critical revisions, and approved the final manuscript.