Research article | Open | Open Peer Review | Published:
Quantifying the risk of local Zika virus transmission in the contiguous US during the 2015–2016 ZIKV epidemic
BMC Medicinevolume 16, Article number: 195 (2018)
Local mosquito-borne Zika virus (ZIKV) transmission has been reported in two counties in the contiguous United States (US), prompting the issuance of travel, prevention, and testing guidance across the contiguous US. Large uncertainty, however, surrounds the quantification of the actual risk of ZIKV introduction and autochthonous transmission across different areas of the US.
We present a framework for the projection of ZIKV autochthonous transmission in the contiguous US during the 2015–2016 epidemic using a data-driven stochastic and spatial epidemic model accounting for seasonal, environmental, and detailed population data. The model generates an ensemble of travel-related case counts and simulates their potential to have triggered local transmission at the individual level in the 2015–2016 ZIKV epidemic.
We estimate the risk of ZIKV introduction and local transmission at the county level and at the 0.025° × 0.025° cell level across the contiguous US. We provide a risk measure based on the probability of observing local transmission in a specific location during a ZIKV epidemic modeled after the epidemic observed during the years 2015–2016. The high spatial and temporal resolution of the model allows us to generate statistical estimates of the number of ZIKV introductions leading to local transmission in each location. We find that the risk was spatially heterogeneously distributed and concentrated in a few specific areas that account for less than 1% of the contiguous US population. Locations in Texas and Florida that have actually experienced local ZIKV transmission were among the places at highest risk according to our results. We also provide an analysis of the key determinants for local transmission and identify the key introduction routes and their contributions to ZIKV transmission in the contiguous US.
This framework provides quantitative risk estimates, fully captures the stochasticity of ZIKV introduction events, and is not biased by the under-ascertainment of cases due to asymptomatic cases. It provides general information on key risk determinants and data with potential uses in defining public health recommendations and guidance about ZIKV risk in the US.
From 2015 to 2016, the Zika virus (ZIKV) epidemic spread across most countries in the Americas, including the United States (US) [1,2,3]. As of July 3, 2018, three US territories, including Puerto Rico, have reported 37,255 ZIKV cases mostly due to widespread local transmission [3, 4]. Laboratory evidence of possible ZIKV infections has been found in 4900 pregnant women from US territories, 167 of whom have had pregnancy outcomes with ZIKV-related birth defects [3, 5, 6]. The US states and District of Columbia have reported 5710 travel-associated ZIKV cases, including 2474 pregnant women with evidence of ZIKV infection and 116 ZIKV-related birth defects . Two geographical locations have experienced local transmission of ZIKV in the contiguous US: Miami-Dade County, in Florida, and Cameron County, in Texas [7, 8]. While the outbreaks in Florida and Texas were limited, the indirect impact on the local economy has been remarkable .
Concerns have been raised that several other locations in the contiguous US were at risk of ZIKV transmission, thus triggering a number of studies aimed at identifying populations at highest risk of local transmission [10,11,12,13,14,15,16,17,18,19,20]. In particular, detailed studies based on environmental suitability, epidemiological factors, and travel-related case importations have been used to estimate the risk for specific counties in the US [21, 22]. In this study, we quantify the risk of local ZIKV transmission by using a data-driven stochastic and spatial epidemic model accounting for seasonal, environmental, and detailed population data. The model also accounts for the association between socioeconomic status and the risk of exposure to mosquitoes, and it has been previously used to estimate the introduction of Zika in the Americas and the spatial and temporal dynamics of the epidemic . By using an extensive likelihood analysis with data from places with a reliable epidemiological surveillance system, the model generates a stochastic ensemble of simulations estimating the place and time of introduction of ZIKV in Brazil and the unfolding of the epidemic in the Americas. For each simulation, the individual-level scale of the model allows for the construction of daily travel-related case counts (TCCs) tracking the number of infected individuals in the contiguous US at the county level and at the finer spatial resolution of 0.025° × 0.025° corresponding approximately to 2.5 km × 2.5 km cells, comparable in size to the ZIKV active transmission areas identified in Florida by the Centers for Disease Control and Prevention (CDC) . Using the time series of county-specific TCCs and the mechanistic transmission model, it is possible to estimate the probability that a specific location would experience local ZIKV transmission during the 2015–2016 time window. The methodology proposed here provides a statistical estimate of ZIKV transmission risk that is not biased by the under-ascertainment of infections and the single historical occurrence of the case importation timeline that fail to account for the full stochasticity of transmission. The TCC database also allows us to identify key sources and routes of ZIKV introductions. Results from our study can provide guidance to public health agencies in their efforts to identify populations and seasons at high risk of ZIKV transmission, so that resources towards outbreak prevention and response can be allocated more efficiently.
We consider three major factors associated with local ZIKV transmission in the contiguous US: the intensity of travel-related infection importations, the environmental suitability for ZIKV transmission, and the socioeconomic risk of exposure to mosquitoes. In this study, we develop a data-driven computational framework (Fig. 1) to quantitatively account for these three factors and to evaluate their impact on ZIKV transmission. Based on this framework, we assess the risk of local ZIKV transmission across the contiguous US through the full course of the 2015–2016 ZIKV epidemic.
The starting point of our methodology is the construction of a synthetic database of TCC entering the US through airport transportation hubs. The database is generated from simulations based on a large-scale spatial model simulating the 2015–2016 ZIKV epidemics, where both symptomatic and asymptomatic ZIKV infections are considered . The synthetic database of TCC contains for each infected individual the time of arrival, stage of ZIKV infection, airports of origin and arrival, and location of residence in the contiguous USFootnote 1 . A schematic sample of the database is shown in Table 1.
Each infected individual’s likelihood of exposure to mosquito bites and his/her capability of triggering local ZIKV transmission is affected by the ecological presence of mosquitoes in his/her location of residence. Indeed, our model integrates mosquito abundance data (Ae. aegypti and Ae. albopictus) [26, 27] that takes into account for temperature suitability, precipitation, vegetation, and urbanization and considers seasonal variations in the mosquito density determined by daily temperature. The individual’s socioeconomic status, which is strongly associated with factors such as sanitation conditions, accessibility to air conditioning, and level of disease awareness, also affects the likelihood of exposure to mosquitoes [14, 28, 29]. Our computational framework considers a data layer based on global socioeconomic indicators , which is calibrated with historical mosquito-borne disease outbreaks in naive populations to provide a likelihood map of the individual’s exposure to mosquitoes . This map serves as a spatial filter (Fig. 1c-II) that probabilistically selects individuals exposed to mosquito bites down to the resolution of a 0.25° × 0.25° cell containing his/her location of residence. Each of the exposed individuals can potentially trigger detectable local ZIKV transmissions (Fig. 1c-III, d), according to the stochastic mechanistic ZIKV transmission model that takes into account mosquito abundance, the current temperature in the area, and the transmission dynamics of ZIKV (see Additional file 1: Supplementary Information). We define a detectable local transmission as the generation of 20 or more autochthonous transmission infections triggered by a single ZIKV infection introduction. Smaller outbreaks would likely go unnoticed assuming a 5% to 10% detection rate of infections due to the large proportion of asymptomatic cases [31,32,33]. Due to fine spatial and temporal resolution, the transmission model is able to account for the significant variability in the ZIKV basic reproduction number (R0) across locations, as well as the variability within the same location at different times. These differences in R0 are driven by temperature and the mosquito abundance, among other variables. The details of the mechanistic model and the calculation of the socioeconomic risk of exposure to mosquitoes are reported in Additional file 1. More technically, we can define the following procedure:
We randomly sample one out of the simulated TCC from the statistical ensemble output of the ZIKV model .
For each infected individual in the TCC, we stochastically determine whether he/she is potentially exposed to mosquito bites based on the probability of exposure pe at the location of residence x. pe is calibrated based on socioeconomic indicators and x identifies a specific county or spatial cell. In each location x, these individuals could potentially trigger local transmission.
Based on the individual’s stage of infection (exposed or infectious), time of introduction, and location of residence (at 0.025° × 0.025° resolution), we simulate local ZIKV transmission with the same stochastic transmission model used in the global model (described in Additional file 1: Supplementary Information) with the specific parameters calibrated to each 0.25° × 0.25° cell in the US.
For each simulated TCC, the above procedure identifies all the infections triggering detectable local transmission. For every time interval ∆t and geographical area x of interest, we can associate variable n(x, ∆t) = 1 if there is at least one imported infection from the TCC that triggers detectable local transmission, and n(x, ∆t) = 0 otherwise.
In order to provide a probabilistic risk measurement, we execute N = 106 resamplings from the ensemble of simulated TCC generated by the model and repeat the above procedure. The resampling procedure accounts for the many possible TCCs compatible with the observed ZIKV epidemic and stochastic effects in the local transmission. This is because not all case importations will result in local outbreaks, even in areas where transmission is favored. The risk of local ZIKV transmission for area x during time window ∆t can be thus defined as
where i indexes the 106 outcomes from the resampled TCCs. This definition of the risk can be aggregated at various spatial (0.025° × 0.025°) and temporal resolutions (≥ 1 day), and it can be used to generate risk maps of ZIKV introduction across the contiguous US. Unless otherwise specified, we consider in this study the local transmission risk rtr(x) that is defined on the ∆t referring to the time window spanning from January 1, 2015 to December 31, 2016. This definition of risk can be interpreted as the probability of observing a detectable local transmission in a specific area per ZIKV epidemic.
By using the methodology outlined in the previous section, we provide quantitative estimates of rtr(x) both at the county level and at 0.025° × 0.025° cell resolution. Figure 2a shows the risk of ZIKV introduction at county level in the contiguous US through the full course of the simulated 2015–2016 ZIKV epidemics. We consider four main brackets for the risk and the associated population sizes. At the county level, the highest risk bracket rtr(x) > 0.5 includes only 0.71% of the total population in the contiguous US. In these areas, one would expect to observe detectable local transmission events with a probability above 50% during the simulated 2015–2016 ZIKV epidemic. Even when we extend the high-risk bracket to include counties with rtr(x) > 1/8, this includes only 2.56% of the total population in the contiguous US. Thus, the risk of local transmission is extremely concentrated to specific geographical locations. Figure 2d shows the population living in counties with different risk brackets of ZIKV introduction and their percentage with respect to the total population in the contiguous US.
The counties of Miami-Dade, Florida, and Cameron, Texas, where local transmission was observed in the year 2016, were both estimated to be high-risk locations (risk bracket, greater than 1/4). Densely populated areas along the Gulf Coast also show up as high-risk locations, in agreement with estimates from other models . The risk of ZIKV introduction and local transmission rtr(x) is highly spatially heterogeneous (Fig. 2a, b). This heterogeneity persists even within the state of Florida, where most areas are estimated to be environmentally suitable for ZIKV transmission all year long [12, 34]. This is mostly because of socioeconomic and local climate heterogeneities. At a spatial granularity of 0.025° × 0.025°, it is possible to perform a statistical analysis of the risk distribution. In Fig. 2c, we report the distribution of cell-specific risks rtr(x). The distribution has a very right-skewed heavy tail extending over more than four orders of magnitude, a clear signature of the large heterogeneity of the risk in the contiguous US.
It is worth stressing that the source of ZIKV introductions in each location is time-dependent, since the TCC is determined by both the magnitude of the epidemic in the regions of the Americas affected by ZIKV and travel patterns from these areas. Our model explicitly simulates individual ZIKV-infected travelers, with detailed information about the traveler’s origin and destination at the daily scale. This allows us to decompose the relative contribution of potential ZIKV introductions from different epidemic regions and to identify routes of high risk with high spatiotemporal resolution. In Table 2, we report the likelihood of local ZIKV transmission in Miami-Dade, Florida, for the year 2015 and 2016 triggered by infection importations from the Caribbean, Central America and Mexico, and South America. The likelihood accounts for intensity of ZIKV transmission in epidemic regions, the travel volume between the source regions and Miami-Dade, and the time-dependent environmental suitability of local transmission in Miami-Dade. In Fig. 3, we report the daily risk of ZIKV infections in Miami-Dade from different geographical regions as well as the time-dependent relative contributions of different regions to the risk throughout the years 2015 and 2016.
As shown in both Table 2 and Fig. 3, in 2015, countries in the Caribbean and South America were major contributors to ZIKV introduction risk in Miami-Dade. On the other hand, countries in Central America and Mexico became major contributors in 2016. This reflects the fact that the ZIKV epidemic started earlier in South American countries, including Brazil and Colombia, and later on spread to countries in Central America and Mexico. Caribbean countries, however, remained a major source of infection importation in both 2015 and 2016. This is possibly due to the high travel volumes between Florida and the Caribbean, as well as high incidence rate and weak seasonality of ZIKV transmission in that region. This is in line with epidemiological data from Florida’s Department of Health, as well as phylogenetic analysis based on sequenced ZIKV genomes from both infected humans and mosquitoes in Florida .
In Fig. 4, we zoom in on three representative areas to disentangle the key determinants shaping the spatiotemporal risk of local ZIKV transmission. Panels a, b, and c in Fig. 4 represent geographical areas covering Miami-Dade, Florida; Cameron, Texas; and New York City, New York. Both Miami-Dade and New York City experienced a high volume of ZIKV infection importations due to high population density and close proximity to major international transportation hubs. Cameron, Texas, on the other hand, had far fewer ZIKV infection importations. However, due to socioeconomic factor (among other factors), the population in Cameron, Texas, is more likely to be exposed to mosquitoes than the populations of Miami-Dade and New York City. Consequently, the volume of Cameron’s imported infections that were exposed to mosquito bites is comparable to those of Miami-Dade and New York City.
The environmental suitability of ZIKV transmission in the three areas is remarkably different. The basic reproduction number R0 is above the epidemic threshold (R0 > 1) in Miami-Dade throughout the year, indicating ZIKV transmission is environmentally suitable all year long. Cameron, Texas, has moderate environmental suitability, where R0 drops below the threshold in winter seasons. New York City is far less environmentally suitable for ZIKV transmission, with a narrow time window of approximately 2 months during summer when R0 is larger than 1.
Given the individual-level resolution of the model, we can focus on the daily average number of travel-associated ZIKV infections leading to local transmission. This is a different indicator than risk. The latter is defined as the probability of observing at least one event of detectable local transmission in the area, thus overlooking the number of different introduction events that trigger local transmission. The profile of daily ZIKV introductions that would lead to local transmission (Fig. 4, row 4) is jointly shaped by ZIKV infection importations, socioeconomic risk of exposure to mosquitoes, and the environmental suitability of ZIKV transmission. The cumulative number of ZIKV introductions leading to local transmission was high in both Miami-Dade, Florida, and Cameron, Texas, where local transmission occurred in the year 2016. The time of ZIKV introduction in Miami-Dade, Florida, is estimated to have occurred between January and May 2016 based on phylogenetic analysis of sequenced ZIKV genomes from infected patients and Ae. aegypti mosquitoes . Our model suggests (Fig. 4, row 4) high risk of ZIKV introduction during the same time window, despite relatively low environmental suitability. The high risk of introduction in Miami-Dade between January and May 2016 was mainly driven by a high influx of imported ZIKV infections. Based on our simulations, Miami-Dade county has on average 1.29 cumulative introductions leading to local transmission events (95%CI (0–9)) throughout 2015 and 2016 (Fig. 4, row 5, insert). However, the distribution of the number of introductions is positively skewed (skewness γ1 = 4.40), with a maximum of 55 introductions. This indicates the possibility of multiple introductions during the ZIKV outbreak in Miami-Dade, Florida, in line with estimates from phylogenetic analysis .
To investigate to what extent the spatial variation of local ZIKV transmission is driven by key socioeconomic and environmental determinants, we first consider a regression model exploring the relation between the average number of local ZIKV transmissions (log(ntr) is the dependent variable) and three key determinants: the number of ZIKV importations, average temperature, and the GDP per capita. Specifically, the explanatory variables include:
log(Nim), the logarithm of the cumulative average number of TCC for each 0.25° × 0.25° cell from January 1, 2015, to December 31, 2016.
log(f20°), the logarithm of the fraction of days over the year with an average temperature larger than 20 °C for each 0.25° × 0.25° cell.
log(GDP), the gross domestic product per capita in terms of purchasing power parity for each 0.25° × 0.25° cell.
In Table 3, we show that if all three explanatory variables are included in the regression (model 1), the model can explain 73.9% of the variance in the number of average introductions leading to local transmission in each cell x. While only considering log(Nim) and log(f20° ) (model 2), we can explain 56.2% of the variance, and using log(Nim) (model 3) alone can explain 47.5% of the variance. It is worth remarking that such a simple statistical analysis cannot fully explain the variance of log(ntr) due to the nonlinear dependency between ZIKV transmission, vector population dynamics, and temperature. It is also due to the highly nonlinear nature of the disease transmission dynamics captured by the epidemic threshold (where the basic reproduction number (R0) needs to be larger than one to be able to spread in a population). In addition, more than 90% of the geographical areas in the contiguous US are not included in the regression because the simulations project no local transmission events in those areas. However, 77% (in terms of areas) of these “risk-free” areas are not environmentally suitable for ZIKV transmission according to our model.
To better illustrate the role of the three main drivers of Zika transmission, we conduct a sensitivity analysis considering three counterfactual scenarios. In each counterfactual scenario, we modify one of the three drivers across the contiguous US to uniformly mimic the conditions in Miami-Dade, Florida, while keeping the other two drivers intact. Specifically:
In counterfactual scenario 1, the environmental suitability (the temperature and thus all temperature-modulated disease parameters) and socioeconomic risk of exposure remain the same, while for all airports in the US, the ZIKV infection importations are set to be the same as those of the airport in Miami-Dade, Florida.
In counterfactual scenario 2, the ZIKV infection importations and the socioeconomic risk of exposure to mosquitoes remain the same. However, in this scenario, the temperature and consequently all temperature-modulated parameters of ZIKV transmission model across the contiguous US are set to be the same as those in Miami-Dade, Florida.
In counterfactual scenario 3, the ZIKV infection importations and the environmental suitability are kept intact, while the socioeconomic risks of exposure to mosquitoes across the contiguous US are set to be the same as that in Miami-Dade, Florida.
For each of the three counterfactual scenarios, we repeat the analysis performed with the real data and generate the cumulative county-level risk map of local ZIKV transmission during the years 2015–2016 (see Additional file 1: Supplementary Information, Section 4). All three risk maps of the counterfactual scenarios are distinctly different from the risk map of Fig. 2a. Particularly, in counterfactual scenario 1, under unrealistic high intensity of ZIKV infection importations, all areas with overlapping favorable environmental and socioeconomic determinants are at high risk of local ZIKV transmission. In counterfactual scenario 2, with unrealistic favorable environmental suitabilities of ZIKV across the US, the areas at high risk are no longer restricted to the proximity of the US southern border. Many counties with low average yearly temperature and absence of Aedes mosquitos in the real world present significant risks of local ZIKV transmission. In counterfactual scenario 3, with high socioeconomic status equivalent to Miami-Dade, Florida, both southern Texas and populated areas along the Gulf Coast were relieved from high probability of encountering Zika, leaving southern Florida as the only focus of high risk. Thus, all three drivers are necessarily required to evaluate the risk of local ZIKV transmission in the contiguous US.
A prominent feature of our findings is the spatiotemporal heterogeneity of ZIKV transmission risk across the contiguous US. Spatially, our model estimates that approximately 68.9% of the people in the contiguous US live in areas that are environmentally suitable for ZIKV transmission, in line with other models’ estimates . However, taking all ZIKV introduction and transmission determinants into consideration, areas with non-negligible risk (greater than 1/8) are concentrated in densely populated areas along the Gulf Coast, capturing 2.56% of the US population. From a temporal perspective, certain areas experience strong seasonality of ZIKV environmental suitability, with a narrow time window when ZIKV transmission is possible. Given limited resources, identifying seasons and regions of high risk may help guide resource allocation for high-risk population screening, intervention, and vector control. Our model is also able to identify the high-risk routes of ZIKV importations through air travel. Imported infections originating from Caribbean countries served as a major contributor to trigger local ZIKV transmission in Florida. Although it has the highest number of estimated ZIKV infections among all countries, Brazil is not a major contributor overall (5.75% of potential introductions leading to local transmission across the contiguous US). This is due to Rio de Janeiro and Sao Paulo, two of the largest transportation hubs in Brazil which make up 65% of the international travel to US from Brazil, being located in the Southern region where ZIKV transmission activity is relatively low. In addition, Rio de Janeiro and Sao Paulo have the opposite seasonality compared to the contiguous US. When it is environmentally suitable for ZIKV transmission in Rio de Janerio and Sao Paulo, it is not suitable for ZIKV transmission in most of the US. Thus, imported ZIKV infections from Brazil were less likely to fuel potential transmissions in the US.
Our model also suggests that in Miami-Dade, Florida, the overall risk of ZIKV introduction in 2015 is comparable to that in 2016, while local transmission is only observed in 2016. This could be explained by the stochasticity of transmission events. Another possibility is that because of the high asymptomatic rate of ZIKV infections, limited local transmission events occurred in 2015 without being picked up by the surveillance system. Awareness of ZIKV was low in 2015 as the World Health Organization declared ZIKV as a Public Health Emergency of International Concerns only in early 2016. Around the same time, the CDC announced a Health Alert Network advisory for Zika virus , marking the start of active monitoring of ZIKV activities in the US.
The proposed model has several limitations. The high volume of cruise ship stops along coastal areas of Florida to the Caribbean may elevate the risk of ZIKV transmissions beyond what is estimated in our model. Sexual transmission and transmission through other routes, not considered by our model, may facilitate the risk of local transmission even further. From January 1, 2015, to August 9, 2017, there were 49 reported ZIKV cases in the contiguous US acquired through other routes, including sexual transmission [3, 37,38,39]. This indicates that a larger population may be affected by ZIKV [40,41,42]. In addition, ZIKV RNA was detected in semen as long as 92 days after symptom onset and is able to be sexually transmitted 31–42 days after symptom onset . ZIKV’s ability to persist in infected males and the potential to infect through sexual transmission long after symptom onset are troublesome. However, the specific risk through sexual transmission or other transmission routes are not well understood, and the overall impact of ZIKV infections acquired through other routes remains unclear. As such, we do not include them in our study . Risk of exposure to mosquitoes associated with socioeconomic factors is widely recognized but poorly quantified. In our model, we utilize seroprevalence studies from nine chikungunya outbreaks on confined, naive populations to estimate this association, in line with other approaches used to estimate the ZIKV attack rate . Further studies however are needed to advance our understanding of the association between risk of exposure to mosquitoes and socioeconomic status.
Our model assumes the mosquito abundance is explicitly modulated by temperature, since many studies suggest that temperature is the main driver of the seasonal variation of mosquito abundance [45,46,47,48]. The effect of rainfall as an environmental driver is indirectly included into our model through incorporating the mosquito presence data created by Kraemer et al. . The study suggests that for both Ae. aegypti and Ae. albopictus, maximum and minimum precipitation make significant contributions to explain the spatial distribution of Aedes mosquitoes, consequently affecting the environmental suitability of local ZIKV transmission. However, a full mechanistic modeling of the influence of rainfall (i.e., daily timescale) on the mosquito lifecycle, while interesting, is still out of reach on a global scale. Along with rainfall, human water supplies may also affect the availability of stagnant water, especially in urban settings [45, 49]. Without controlling for the effect of human water supplies, the effect of precipitation could be positive [50,51,52], negative , or no effect at all [54, 55]. In Additional file 1: Supplementary Information, we provide a figure illustrating the seasonal abundance provided by our model.
In our model, we consider both Ae. aegpyti and Ae. albopictus as competent vectors to transmit the ZIKV. However, the competence of Ae. albopictus to transmit ZIKV is debated, and the notable differences in the spatial distributions of Ae. aegpyti and Ae. albopictus make it crucial for evaluating the global risk of ZIKV [27, 56]. However, these differences are less relevant when limiting the risk assessment within the spatial range of the contiguous US. This is because the geographical distribution of the environmental suitability of Ae. aegpyti and Ae. albopictus is largely overlapping within the contiguous US, based on the studies by Johnson et al. . The environmental suitability distribution of Ae. albopictus extends a bit further north when compared to that of Ae. aegpyti. In the areas where only Ae. albopictus are present, the overall environmental suitability of ZIKV transmission is very low due to the presence of strong seasonality, and our model estimates that those areas would have minimal risk of experiencing local ZIKV transmission in the years 2015–2016 (Fig. 2).
In 2017–2018, ZIKV transmission activities in most countries throughout the Americas has plummeted , in agreement with model estimates [23, 58]. The risk of ZIKV introduction in the contiguous US would be expected to be negligible as imported infections triggering the local transmission would be drastically reduced. However, one should exercise caution as vector-transmitted diseases are known to show strong spatial heterogeneity and seasonality and are affected by socioeconomic factors. The stochastic nature of ZIKV transmission could leave a considerable amount of naive populations living in regions at risk of ZIKV transmission. Furthermore, expansion of the Aedes mosquito distribution, human migration, and shifts in socioeconomic status could lead to more populations being at risk for local ZIKV transmission. It is more likely that ZIKV transmission activities in the future may resemble the current situation of chikungunya, where transmission activities could flare up sporadically. The possible sporadic outbreaks of ZIKV would continue to pose a risk to the contiguous US, where most of the population is naive to the virus and a large fraction live in areas environmentally suitable for ZIKV transmission.
In this study, we show that the overall risk of ZIKV introduction and local transmission during the 2015–2016 outbreak is jointly determined by the intensity of ZIKV importations, environmental suitability for ZIKV transmissions, and the socioeconomic risk of exposure to mosquitoes. Our estimates suggest that the risk of ZIKV introductions has a very strong spatial and temporal heterogeneity. The areas in the contiguous US at non-negligible risk (that is, greater than 1/8) only account for 2.6% of the total population in the contiguous US. The model is able to identify the hotspots for ZIKV introductions, and it reveals the relative contributions of ZIKV introductions from different geographical regions over time. The results of our study have the potential to guide the development of ZIKV prevention and response strategies in the contiguous US.
Although there has been reporting in the media about the traffic to and from Latin and Caribbean countries, airline traffic in 2016 has been stable with a mere 4.4% increase.
Centers for Disease Control and Prevention
Gross domestic product
Purchasing power parity
Travel-related case counts
World Health Organization. Zika virus and complications: 2016 Public Health Emergency of International Concern: World Health Organization; 2017. http://www.who.int/emergencies/zika-virus/en/ Accessed 17 Sept 2017
Pan American Health Organization. Zika virus infection: Pan American Health Organization; 2017. https://goo.gl/oABMtf Accessed 17 Sept 2017
Centers for Disease Control and Prevention. Zika virus: Centers for Disease Control and Prevention; 2017. https://www.cdc.gov/zika/index.html Accessed 17 Sept 2017
Walker WL. Zika virus disease cases—50 states and the District of Columbia, January 1–July 31, 2016. Morb Mortal Wkly Rep. 2016;65(36):983.
Reynolds MR, Jones AM, Petersen EE, Lee EH, Rice ME, Bingham A, Ellington SR, Evert N, Reagan-Steiner S, Oduyebo T, et al. Vital signs: update on Zika virus-associated birth defects and evaluation of all US infants with congenital Zika virus exposure-US Zika Pregnancy Registry, 2016. Morb Mortal Wkly Rep. 2017;66(13):366–73.
Honein MA, Dawson AL, Petersen EE, Jones AM, Lee EH, Yazdy MM, Ahmad N, Macdonald J, Evert N, Bingham A, et al. Birth defects among fetuses and infants of US women with evidence of possible Zika virus infection during pregnancy. JAMA. 2017;317(1):59–68.
Centers for Disease Control and Prevention. Advice for people living in or traveling to Brownsville, Texas: Centers for Disease Control and Prevention; 2017. https://www.cdc.gov/zika/intheus/texas-update.html Accessed 17 Sept 2017
Centers for Disease Control and Prevention. Advice for people living in or traveling to South Florida: Centers for Disease Control and Prevention; 2017. https://www.cdc.gov/zika/intheus/florida-update.html Accessed 17 Sept 2017
Lee BY, Alfaro-Murillo JA, Parpia AS, Asti L, Wedlock PT, Hotez PJ, Galvani AP. The potential economic burden of Zika in the continental United States. PLoS Negl Trop Dis. 2017;11(4):0005531.
Keegan LT, Lessler J, Johansson MA. Quantifying Zika: advancing the epidemiology of Zika with quantitative models. J Infect Dis. 2017;216(suppl 10):884–90.
Monaghan AJ, Morin CW, Steinhoff DF, Wilhelmi O, Hayden M, Quattrochi DA, Reiskind M, Lloyd AL, Smith K, Schmidt CA, et al. On the seasonal occurrence and abundance of the Zika virus vector mosquito Aedes aegypti in the contiguous United States. PLoS Curr. 2016;8. https://doi.org/10.1371/currents.outbreaks.50dfc7f46798675fc63e7d7da563da76.
Messina JP, Kraemer MU, Brady OJ, Pigott DM, Shearer FM, Weiss DJ, Golding N, Ruktanonchai CW, Gething PW, Cohn E, et al. Mapping global environmental suitability for Zika virus. Elife. 2016;5:15272.
Bogoch II, Brady OJ, Kraemer MU, German M, Creatore MI, Brent S, Watts AG, Hay SI, Kulkarni MA, Brownstein JS, et al. Potential for Zika virus introduction and transmission in resource-limited countries in Africa and the Asia-Pacific region: a modelling study. Lancet Infect Dis. 2016;16(11):1237–45.
Perkins TA, Siraj AS, Ruktanonchai CW, Kraemer MU, Tatem AJ. Model-based projections of Zika virus infections in childbearing women in the Americas. Nat Microbiol. 2016;1:16126.
Alfaro-Murillo JA, Parpia AS, Fitzpatrick MC, Tamagnan JA, Medlock J, Ndeffo-Mbah ML, Fish D, Ávila-Agüero ML, Marín R, Ko AI, et al. A cost-effectiveness tool for informing policies on Zika virus control. PLoS Negl Trop Dis. 2016;10(5):0004743.
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.
Rocklöv J, Quam MB, Sudre B, German M, Kraemer MU, Brady O, Bogoch II, Liu-Helmersson J, Wilder-Smith A, Semenza JC, et al. Assessing seasonal risks for the introduction and mosquito-borne spread of Zika virus in Europe. EBioMedicine. 2016;9:250–6.
Lourenço J, de Lima MM, Faria NR, Walker A, Kraemer MU, Villabona-Arenas CJ, Lambert B, de Cerqueira EM, Pybus OG, Alcantara LC, et al. Epidemiological and ecological determinants of Zika virus transmission in an urban setting. eLife. 2017;6:e29820.
Ajelli M. Modeling mosquito-borne diseases in complex urban environments. Acta Trop. 2017;176:332–4.
Ajelli M, Moise IK, Hutchings TCS, Brown SC, Kumar N, Johnson NF, Beier JC. Host outdoor exposure variability affects the transmission and spread of Zika virus: insights for epidemic control. PLoS Negl Trop Dis. 2017;11(9):0005851.
Castro LA, Fox SJ, Chen X, Liu K, Bellan SE, Dimitrov NB, Galvani AP, Meyers LA. Assessing real-time Zika risk in the United States. BMC Infect Dis. 2017;17(1):284.
Fox SJ, Bellan SE, Perkins TA, Johansson MA, Meyers LA. Downgrading disease transmission risk estimates using terminal importations. bioRxiv. 2018:265942. https://doi.org/10.1101/265942.
Zhang Q, Sun K, Chinazzi M, y Piontti AP, Dean NE, Rojas DP, Merler S, Mistry D, Poletti P, Rossi L, et al. Spread of Zika virus in the Americas. Proc Natl Acad Sci. 2017;114(22):4334–43.
Centers for Disease Control and Prevention. Guidance for areas with local Zika virus transmission in the continental United States and Hawaii: Centers for Disease Control and Prevention; 2017. https://www.cdc.gov/zika/geo/domestic-guidance.html Accessed 17 Sept 2017
The International Civil Aviation Organization (ICAO): The world of air transport in 2016 (2016). https://www.icao.int/annual-report-2016/Pages/the-world-of-air-transport-in-2016.aspx Accessed 17 Sept 2017.
Centers for Disease Control and Prevention. Estimated range of Aedes aegypti and Aedes albopictus in the United States, 2017: Centers for Disease Control and Prevention; 2017. https://www.cdc.gov/zika/vector/range.html Accessed 17 Sept 2017
Kraemer MU, Sinka ME, Duda KA, Mylne AQ, Shearer FM, Barker CM, Moore CG, Carvalho RG, Coelho GE, Van Bortel W, et al. The global distribution of the arbovirus vectors Aedes aegypti and Ae. albopictus. Elife. 2015;4:08347.
Sissoko D, Moendandze A, Malvy D, Giry C, Ezzedine K, Solet JL, Pierre V. Seroprevalence and risk factors of chikungunya virus infection in Mayotte Indian Ocean, 2005-2006: a population-based survey. PLoS One. 2008;3(8):3066. https://doi.org/10.1371/journal.pone.0003066.
Reiter P, Lathrop S, Bunning M, Biggerstaff B, Singer D, Tiwari T, Baber L, Amador M, Thirion J, Hayes J, et al. Texas lifestyle limits transmission of dengue virus. Emerg Infect Dis. 2003;9(1):86.
Nordhaus, W.D., Chen, X.: Global gridded geographically based economic data (G-econ), version 4 (2016).
Duffy MR, Chen T-H, Hancock WT, Powers AM, Kool JL, Lanciotti RS, Pretrick M, Marfel M, Holzbauer S, Dubray C, et al. Zika virus outbreak on Yap Island, federated states of Micronesia. N Engl J Med. 2009;360(24):2536–43.
Russell SP, Ryff KR, Gould CV, Martin SW, Johansson MA. Detecting local Zika virus transmission in the continental United States: a comparison of surveillance strategies. PLOS Currents Outbreaks. 2017:145102. Edition 1. https://doi.org/10.1371/currents.outbreaks.cd76717676629d47704170ecbdb5f820.
Moghadas SM, Shoukat A, Espindola AL, Pereira RS, Abdirizak F, Laskowski M, Viboud C, Chowell G. Asymptomatic transmission and the dynamics of Zika infection. Sci Rep. 2017;7(1):5829.
Marini G, Guzzetta G, Rosà R, Merler S. First outbreak of Zika virus in the continental United States: a modelling analysis. Euro Surveill. 2017;22(37).
Grubaugh ND, Ladner JT, Kraemer MU, Dudas G, Tan AL, Gangavarapu K, Wiley MR, White S, Thézé J, Magnani DM, et al. Genomic epidemiology reveals multiple introductions of Zika virus into the United States. Nature. 2017;546(7658):401–5.
Bogoch II, Brady OJ, Kraemer M, German M, Creatore MI, Kulkarni MA, Brownstein JS, Mekaru SR, Hay SI, Groot E, et al. Anticipating the international spread of Zika virus from Brazil. Lancet. 2016;387(10016):335–6.
Foy BD, Kobylinski KC, Foy JLC, Blitvich BJ, da Rosa AT, Haddow AD, Lanciotti RS, Tesh RB. Probable non–vector-borne transmission of Zika virus, Colorado, USA. Emerg Infect Dis. 2011;17(5):880.
Russell K, Hills SL, Oster AM, Porse CC, Danyluk G, Cone M, Brooks R, Scotland S, Schiffman E, Fredette C, et al. Male-to-female sexual transmission of Zika virus—United States, January–April 2016. Clin Infect Dis. 2016;64(2):211–3.
McCarthy M. Zika virus was transmitted by sexual contact in Texas, health officials report. BMJ. 2016;352:i720.
Gao D, Lou Y, He D, Porco TC, Kuang Y, Chowell G, Ruan S. Prevention and control of Zika as a mosquito-borne and sexually transmitted disease: a mathematical modeling analysis. Sci Rep. 2016;6:28070.
Allard A, Althouse BM, Hébert-Dufresne L, Scarpino SV. The risk of sustained sexual transmission of Zika is underestimated. PLoS Pathog. 13(9):1006633–2017.
Yakob L, Kucharski A, Hue S, Edmunds WJ. Low risk of a sexually-transmitted Zika virus outbreak. Lancet Infect Dis. 2016;16(10):1100–2.
Gaskell KM, Houlihan C, Nastouli E, Checkley AM. Persistent Zika virus detection in semen in a traveler returning to the United Kingdom from Brazil, 2016. Emerg Infect Dis. 2017;23(1):137.
Kim CR, Counotte M, Bernstein K, Deal C, Mayaud P, Low N, Broutet N, et al. Investigating the sexual transmission of Zika virus. Lancet Glob Health. 2018;6(1):24–5.
Tran A, L’Ambert G, Lacour G, Benoît R, Demarchi M, Cros M, Cailly P, Aubry-Kientz M, Balenghien T, Ezanno P. A rainfall- and temperature-driven abundance model for Aedes albopictus populations. Int J Environ Res Public Health. 2013;10(5):1698–719. https://doi.org/10.3390/ijerph10051698.
Alto BW, Juliano SA. Precipitation and temperature effects on populations of Aedes albopictus (Diptera: Culicidae): implications for range expansion. J Med Entomol. 2001;38(5):646–56.
Gomes AF, Nobre AA, Cruz OG. Temporal analysis of the relationship between dengue and meteorological variables in the city of Rio de Janeiro, Brazil, 2001-2009. Cad Saude Publica. 2012;28(11):2189–97.
Xu L, Stige LC, Chan K-S, Zhou J, Yang J, Sang S, Wang M, Yang Z, Yan Z, Jiang T, et al. Climate variation drives dengue dynamics. Proc Natl Acad Sci. 2017;114(1):113–8.
Barrera R, Amador M, Mackay AJ. Population dynamics of Aedes aegypti and dengue as influenced by weather and human behavior. PLoS Negl Trop Dis. 2011;5(12):1378. https://doi.org/10.1371/journal.pntd.0001378.
Lourenço-de-Oliveira R, Castro MG, Braks MAH, Lounibos LP. The invasion of urban forest by dengue vectors in Rio de Janeiro. J Vector Ecol. 2004;29:94–100.
Reiskind M, Lounibos L. Spatial and temporal patterns of abundance of Aedes aegypti L.(Stegomyia aegypti) and Aedes albopictus (Skuse)[Stegomyia albopictus (Skuse)] in southern Florida. Med Vet Entomol. 2013;27(4):421–9.
Li M-T, Sun G-Q, Yakob L, Zhu H-P, Jin Z, Zhang W-Y. The driving force for 2014 dengue outbreak in Guangdong, China. PLoS One. 2016;11(11):0166211.
Roiz D, Rosà R, Arnoldi D, Rizzoli A. Effects of temperature and rainfall on the activity and dynamics of host-seeking Aedes albopictus females in northern Italy. Vector Borne Zoonotic Dis. 2010. https://doi.org/10.1089/vbz.2009.0098.
Luciano T, Severini IF, Di Luca IM, Bella IA, ryP Roberto R. Seasonal patterns of oviposition and egg hatching rate of Aedes albopictus in Rome. J Am Mosq Control Assoc. 2003;19(1):100.
Azil AH, Long SA, Ritchie SA, Williams CR. The development of predictive tools for pre-emptive dengue vector control: a study of Aedes aegypti abundance and meteorological variables in North Queensland, Australia. Tropical Med Int Health. 2010;15(10):1190–7.
Gardner LM, Chen N, Sarkar S. Global risk of Zika virus depends critically on vector status of Aedes albopictus. Lancet Infect Dis. 2016;16(5):522–3.
Johnson TL, Haque U, Monaghan AJ, Eisen L, Hahn MB, Hayden MH, Savage HM, McAllister J, Mutebi J-P, Eisen RJ. Modeling the environmental suitability for Aedes (Stegomyia) aegypti and Aedes (Stegomyia) albopictus (Diptera: Culicidae) in the contiguous United States. J Med Entomol. 2017;54(6):1605–14.
Ferguson NM, Cucunubá ZM, Dorigatti I, Nedjati-Gilani GL, Donnelly CA, Basáñez M-G, Nouvellet P, Lessler J. Countering the Zika epidemic in Latin America. Science. 2016;353(6297):353–4.
CIESIN-Columbia University: Global population count grid time series estimates (2017). https://doi.org/10.7927/H4CC0XNV Accessed 31 Dec 2017.
Balcan D, Gonçalves B, Hu H, Ramasco JJ, Colizza V, Vespignani A. Modeling the spatial spread of infectious diseases: the GLobal Epidemic and Mobility computational model. J Comput Sci. 2010;1(3):132–45.
This work was supported by the Models of Infectious Disease Agent Study, National Institute of General Medical Sciences Grant U54GM111274 and the NIH supplement Grant R01 AI102939-05.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
AV has received research support unrelated to this paper (through his employer Northeastern University) from Metabiota Inc. The other authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information. (PDF 5211 kb)