Skip to main content
  • Research article
  • Open access
  • Published:

Quantifying the risk of local Zika virus transmission in the contiguous US during the 2015–2016 ZIKV epidemic

Abstract

Background

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.

Methods

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.

Results

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.

Conclusions

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.

Peer Review reports

Background

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 [3]. 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 [9].

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 [23]. 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) [24]. 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.

Methods

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.

Fig. 1
figure 1

A schematic illustration of the computational framework to assess the risk of ZIKV introductions into the contiguous US. a High-resolution (0.025° × 0.025°2.5 km × 2.5 km) population density map [59] and Voronoi tessellation of the contiguous US into census areas with a major airport transportation hub at each of their centers [60]. b An example of the census area centered at Miami International Airport. c I: Travel-associated ZIKV infections entering the Miami International Airport. Location of residence of each ZIKV infection is randomly assigned with likelihood proportional to the population density within each census area. II: The probabilistic filter of the risk of exposure to mosquitoes due to socioeconomic factors such as housing conditions, sanitation, and disease awareness. III: Spatiotemporal specific ZIKV transmission dynamics are influenced by environmental factors that are temperature sensitive, including the spatial distribution of Aedes mosquitoes, seasonal mosquito abundance, and ZIKV transmissibility. d Compartmental stochastic ZIKV transmission model used to evaluate the environmental suitability of ZIKV transmission. Humans are divided into susceptible SH, exposed EH, infectious IH, and recovered RH compartments, and mosquitoes are divided into susceptible SV, exposed EV, and infectious IV compartments

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 [23]. 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 [25]. A schematic sample of the database is shown in Table 1.

Table 1 A sample of the database containing simulated travel-related ZIKV-infected individuals entering the US

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 [30], which is calibrated with historical mosquito-borne disease outbreaks in naive populations to provide a likelihood map of the individual’s exposure to mosquitoes [23]. 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:

  1. (1)

    We randomly sample one out of the simulated TCC from the statistical ensemble output of the ZIKV model [23].

  2. (2)

    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.

  3. (3)

    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.

  4. (4)

    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

$$ {r}_{tr}\left(x,\Delta t\right)=\frac{1}{N}\sum \limits_{i=1}^N{n}_i\left(x,\Delta t\right) $$
(1)

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.

Results

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.

Fig. 2
figure 2

The cumulative risk of local ZIKV transmission in the contiguous US. The cumulative risk of local ZIKV transmission at different spatial resolutions is evaluated through the full course of the simulated 2015–2016 ZIKV epidemic. a The cumulative risk map of local ZIKV transmission for each county in the contiguous US. The color scale indicates for any given county the probability of experiencing at least one ZIKV outbreak with more than 20 infections (details in Additional file 1). b High spatial resolution estimates (0.025° × 0.025°) of the cumulative risk of local ZIKV transmission through the full course of the simulated 2015–2016 ZIKV epidemic. c The complementary cumulative distribution function of the local ZIKV transmission risk for all 0.025° × 0.025°cells (on a log-log scale). The heavy tail feature of the distribution reflects strong spatial heterogeneity in terms of local ZIKV transmission risk. d The total population in the counties of the US with different risk levels of local ZIKV transmission 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 [12]. 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.

Table 2 The likelihood of a given local ZIKV transmission event in Miami-Dade, Florida, from different geographical regions (Caribbean, South America, Central America and Mexico) for the years 2015 and 2016
Fig. 3
figure 3

A breakdown of local ZIKV transmission events by the geographical origins of travel-associated ZIKV infections in Miami-Dade, Florida. ac The daily average number of ZIKV imported infections per day that trigger outbreaks with more than 20 infections, originating from the Caribbean, Central America and Mexico, and South America. d The relative contributions to the expected number of local ZIKV transmission events by different geographical regions

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 [35].

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.

Fig. 4
figure 4

Factors which co-shape the spatiotemporal risk of local ZIKV transmission in three different regions in the contiguous US. Columns from left to right represent a Miami-Dade, Florida; b Cameron, Texas; and c New York City, New York. Row 1 shows the average daily number of imported ZIKV infections. Note that for Cameron, Texas, the scale on the y-axis is different than that of Miami-Dade, Florida, and NYC, New York. Row 2 shows the average number of imported ZIKV infections that pass through the socioeconomic filter pe and reside in areas potentially exposed to mosquitoes. Row 3 shows the basic reproduction number (weekly average) calculated based on the ZIKV transmission model. Gray-shaded time windows indicate when the basic reproduction number R0 > 1 and sustainable ZIKV transmission is possible. Row 4 shows the expected daily number of ZIKV introductions with the red-shaded time window indicating the estimated time of local ZIKV transmission based on phylogenetic analysis [35]. Row 5 shows the average cumulative number of local ZIKV transmission events since January 1, 2015

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 [35]. 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 [35].

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.

Table 3 Regression analysis between log(ntr) and explanatory variables including log(Nim), log(f20° ), and log(GDP )

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.

Discussion

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 [36]. 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 [3], 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 [43]. 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 [44]. 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 [14]. 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. [27]. 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 [53], 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. [57]. 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 [2], 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.

Conclusion

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.

Notes

  1. 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.

Abbreviations

CDC:

Centers for Disease Control and Prevention

GDP:

Gross domestic product

PPP:

Purchasing power parity

TCCs:

Travel-related case counts

US:

United States

ZIKV:

Zika virus

References

  1. 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

  2. Pan American Health Organization. Zika virus infection: Pan American Health Organization; 2017. https://goo.gl/oABMtf Accessed 17 Sept 2017

  3. 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

  4. 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.

    Article  Google Scholar 

  5. 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.

    Article  Google Scholar 

  6. 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.

    Article  Google Scholar 

  7. 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

  8. 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

  9. 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.

    Google Scholar 

  10. 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.

    Article  Google Scholar 

  11. 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://0-doi-org.brum.beds.ac.uk/10.1371/currents.outbreaks.50dfc7f46798675fc63e7d7da563da76.

  12. 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.

    Article  Google Scholar 

  13. 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.

    Article  Google Scholar 

  14. 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.

    Article  CAS  Google Scholar 

  15. 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.

    Article  Google Scholar 

  16. 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.

    Article  Google Scholar 

  17. 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.

    Article  Google Scholar 

  18. 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.

  19. Ajelli M. Modeling mosquito-borne diseases in complex urban environments. Acta Trop. 2017;176:332–4.

    Article  Google Scholar 

  20. 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.

    Article  Google Scholar 

  21. 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.

    Article  Google Scholar 

  22. Fox SJ, Bellan SE, Perkins TA, Johansson MA, Meyers LA. Downgrading disease transmission risk estimates using terminal importations. bioRxiv. 2018:265942. https://0-doi-org.brum.beds.ac.uk/10.1101/265942.

  23. 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.

    Article  Google Scholar 

  24. 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

  25. 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.

  26. 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

  27. 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.

    Article  Google Scholar 

  28. 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://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0003066.

    Article  CAS  Google Scholar 

  29. 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.

    Article  Google Scholar 

  30. Nordhaus, W.D., Chen, X.: Global gridded geographically based economic data (G-econ), version 4 (2016).

    Google Scholar 

  31. 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.

    Article  CAS  Google Scholar 

  32. 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://0-doi-org.brum.beds.ac.uk/10.1371/currents.outbreaks.cd76717676629d47704170ecbdb5f820.

  33. 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.

    Article  Google Scholar 

  34. 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).

  35. 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.

    Article  CAS  Google Scholar 

  36. 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.

    Article  Google Scholar 

  37. 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.

    Article  Google Scholar 

  38. 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.

    Article  Google Scholar 

  39. McCarthy M. Zika virus was transmitted by sexual contact in Texas, health officials report. BMJ. 2016;352:i720.

  40. 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.

    Article  CAS  Google Scholar 

  41. 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.

  42. 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.

    Article  Google Scholar 

  43. 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.

    Article  Google Scholar 

  44. 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.

    Article  Google Scholar 

  45. 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://0-doi-org.brum.beds.ac.uk/10.3390/ijerph10051698.

    Article  PubMed  PubMed Central  Google Scholar 

  46. 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.

    Article  CAS  Google Scholar 

  47. 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.

    Article  Google Scholar 

  48. 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.

    Article  CAS  Google Scholar 

  49. 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://0-doi-org.brum.beds.ac.uk/10.1371/journal.pntd.0001378.

    Article  Google Scholar 

  50. 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.

  51. 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.

    Article  CAS  Google Scholar 

  52. 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.

    Google Scholar 

  53. 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://0-doi-org.brum.beds.ac.uk/10.1089/vbz.2009.0098.

    Article  Google Scholar 

  54. 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.

    Google Scholar 

  55. 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.

    Article  Google Scholar 

  56. 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.

    Article  Google Scholar 

  57. 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.

    Article  Google Scholar 

  58. 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.

    Article  CAS  Google Scholar 

  59. CIESIN-Columbia University: Global population count grid time series estimates (2017). https://0-doi-org.brum.beds.ac.uk/10.7927/H4CC0XNV Accessed 31 Dec 2017.

  60. 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.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable

Funding

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.

Author information

Authors and Affiliations

Authors

Contributions

KS, MEH, IML, and AV designed the research. KS, QZ, MC, APyP, SM, DM, PP, LR, and AV performed the research. KS, QZ, MC, APyP, DM, NED, DPR, and AV analyzed the data. KS, QZ, MC, APyP, NED, DPR, SM, DM, PP, LR, MEH, IML, and AV wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Alessandro Vespignani.

Ethics declarations

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Competing interests

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.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Supplementary information. (PDF 5211 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sun, K., Zhang, Q., Pastore-Piontti, A. et al. Quantifying the risk of local Zika virus transmission in the contiguous US during the 2015–2016 ZIKV epidemic. BMC Med 16, 195 (2018). https://0-doi-org.brum.beds.ac.uk/10.1186/s12916-018-1185-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12916-018-1185-5

Keywords