Potential for Zika virus transmission by mosquitoes in temperate climates

Mosquito-borne Zika virus (ZIKV) transmission has almost exclusively been detected in the tropics despite the distributions of its primary vectors extending farther into temperate regions. Therefore, it is unknown whether ZIKV's range has reached a temperature-dependent limit, or if it can spread into temperate climates. Using field-collected mosquitoes for biological relevance, we found that two common temperate mosquito species, Aedes albopictus and Ochlerotatus detritus, were competent for ZIKV. We orally exposed mosquitoes to ZIKV and held them at between 17 and 31°C, estimated the time required for mosquitoes to become infectious, and applied these data to a ZIKV spatial risk model. We identified a minimum temperature threshold for the transmission of ZIKV by mosquitoes between 17 and 19°C. Using these data, we generated standardized basic reproduction number R0-based risk maps and we derived estimates for the length of the transmission season for recent and future climate conditions. Our standardized R0-based risk maps show potential risk of ZIKV transmission beyond the current observed range in southern USA, southern China and southern European countries. Transmission risk is simulated to increase over southern and Eastern Europe, northern USA and temperate regions of Asia (northern China, southern Japan) in future climate scenarios.

MSCB, 0000-0002-7510-167X; SM, 0000-0002-2394-5301; GLH, 0000-0002-7567-7185 Mosquito-borne Zika virus (ZIKV) transmission has almost exclusively been detected in the tropics despite the distributions of its primary vectors extending farther into temperate regions. Therefore, it is unknown whether ZIKV's range has reached a temperature-dependent limit, or if it can spread into temperate climates. Using field-collected mosquitoes for biological relevance, we found that two common temperate mosquito species, Aedes albopictus and Ochlerotatus detritus, were competent for ZIKV. We orally exposed mosquitoes to ZIKV and held them at between 17 and 31°C, estimated the time required for mosquitoes to become infectious, and applied these data to a ZIKV spatial risk model. We identified a minimum temperature threshold for the transmission of ZIKV by mosquitoes between 17 and 19°C. Using these data, we generated standardized basic reproduction number R 0 -based risk maps and we derived estimates for the length of the transmission season for recent and future climate conditions. Our standardized R 0 -based risk maps show potential risk of ZIKV transmission beyond the current observed range in southern USA, southern China and southern European countries. Transmission risk is simulated to increase over southern and Eastern Europe, northern USA and temperate regions of Asia (northern China, southern Japan) in future climate scenarios.
The primary vectors of ZIKV are generally considered to be Aedes aegypti and Ae. albopictus [5][6][7], both of which are widespread in tropical regions where ZIKV transmission has been reported. A total of 18 tropical species have been found to be infected in the field, with eight previously confirmed in the laboratory (reviewed in [8][9][10]). While both primary vectors, especially Ae. albopictus, also occur in some temperate regions, autochthonous transmission of virus has so far been largely restricted to the warmer parts of this species range [11].
Temperature has a profound effect on the ability of insect vectors to transmit pathogens, especially viruses. Decreasing environmental temperature increases the time taken for the vector to become infectious after taking an infected blood meal (termed the extrinsic incubation period (EIP)) [12]; for example, Culex tarsalis females infected with West Nile virus (WNV) can become infectious at 5 days post-infection at 30°C, but require 36 days at 14°C [13]. At these low temperatures, there may be insufficient time for the pathogen to replicate, invade the salivary gland and become infectious within the lifespan of the vector, imposing a critical threshold temperature for virus transmission [14,15]. It is essential to understand whether this current limitation in the range of ZIKV is a result of the climate being too cold for virus transmission by vectors, or if ZIKV transmission is theoretically possible in temperate climates but has not yet occurred (e.g. because has not yet become established). With large human populations living in temperate climates it is also important to estimate the extent of ZIKV transmission in such climates, both currently and using predicted future climate scenarios.
Currently, there are no peer-reviewed data that comprehensively show at what temperature ranges ZIKV can be transmitted by temperate mosquitoes such as Ae. albopictus. Aedes albopictus is already well established in temperate Europe, the USA and Asia. Furthermore, while lowtemperature transmission has been studied for the tropical Ae. aegypti [16], no accurate minimum temperature for transmission by temperate-adapted mosquito species has been established. Such data are essential for the estimation of ZIKV transmission risk in temperate regions. Here, we test the competence of wild (or recently colonized) temperate mosquitoes Ae. albopictus and Ochlerotatus detritus and derive the EIP of ZIKV over a wide range of temperatures. We then use these experimental data from Ae. albopictus to model and estimate temperate areas at risk from the virus by combining data on reported Ae. albopictus presence [17] and current and future climate using an ensemble of calibrated general circulation models driven by the standard representative concentration pathways (RCPs) scenarios from the Intergovernmental Panel on Climate Change (IPCC).

(b) Virus/blood meal preparation
At 7 days post-emergence, female adults were removed, transferred and counted into 1 l cylindrical polypropylene DISPO-SAFE containers (Microbiological Supply Company, UK), with a fine mesh covering the container opening and stored for 24 h without access to sugar or water. Blood meals containing ZIKV were provided for 3 h using a Hemotek feeding system and 'Hemotek feeding membrane' (Hemotek Ltd, UK) odorized by rubbing against human skin (greater than 0.4 log decrease in viral titre was observed over this period). The PE243 isolate of ZIKV from Brazil was used for all infection experiments, virus was cultured in Vero E6 cells, using Dulbecco's Modified Eagle Medium (Thermo Fisher Scientific) and 10% fetal bovine serum (Thermo Fisher Scientific), and prepared to a final titre of 1 × 10 6 PFU ml −1 in 50% human blood. Blood was obtained from NHS transfusion service in Speke, Merseyside, UK. Blood was heparinized to prevent coagulation. Unfed adults were removed from the cage, and the fed mosquitoes were incubated at the appropriate temperature and 70% RH for the respective time period with a 12 : 12 light : dark cycle.

(c) Incubation and saliva extraction
Mortality and competence of wild-obtained Oc. detritus and colony Ae. albopictus were tested at six different temperatures: 17°C, 19°C, 21°C, 24°C, 27°C and 31°C. Adult females were sacrificed at eight time points: 0, 5, 7, 10, 14, 17, 21 and 28 days postinfection. As numbers of mosquitoes available were much lower, Cx. pipiens pipiens and Cs. annulata were only incubated at 21°C and tested at 17 days post-infection. Saliva was extracted using the methods described in [18], mosquitoes were allowed to salivate for 15 min into 10 µl light mineral oil (NF/FCC) (Thermo Fisher Scientific), and expelled directly into TRIzol (Thermo Fisher Scientific); saliva secretion was observed under a dissection microscope from all individuals. Once saliva had been extracted, the entire body was also submerged in TRIzol (Thermo Fisher Scientific). Samples were stored at −80°C until RNA extraction.

(d) Detection of virus
RNA was extracted according to TRIzol (Thermo Fisher Scientific) instructions, and cDNA was generated using SuperScript VILO (Thermo Fisher Scientific). Taqman (Thermo Fisher Scientific) qRT-PCR was used to detect the presence of ZIKV RNA in the samples. The primers and probes used were: sense AAR TACACATACCARAACAAAGTGGT, antisense TCCRCTCCCY CTYTGGTCTTG, probe 6Fam-CTYAGACCAGCTGAAR-Tamra; amplification regime, 5 min 95°C, 40 cycles of 15 s at 95°C and 60 s at 60°C [21]. The limit of detection of the analysis was measured by a known titre of live virus submitted to RNA royalsocietypublishing.org/journal/rspb Proc. R. Soc. B 287: 20200119 extraction, cDNA generation and qRT-PCR as per methods described here. Live virus titre was calculated from PFU ml −1 of stock virus, briefly, Vero E6 cells were infected with serial dilutions of ZIKV and incubated under an overlay of DMEM with 2% FCS and 0.6% Avicel (FMC BioPolymer) at 37°C for 5-7 days. Cells were fixed with 4% formaldehyde then stained with Giemsa for visualization. The limit of detection is described as being 10 0.6 PFU, the lowest value where all three triplicates were positive.
To demonstrate that ZIKV is infectious after incubation at 19°C, an additional 30 mosquitoes were fed with a ZIKV containing blood meal using the same methods described earlier. After the blood meal, the mosquitoes were incubated at 19°C for 21 days, before their saliva was extracted. Cytopathic effect (CPE) assays were used to test saliva samples collected for infectious virus. Salivary collections were pooled (into 10 groups of three) and added to Vero E6 cells maintained in DMEM with 2% FBS and 0.05 mg ml −1 gentamycin at 37°C and 5% CO 2 . Cells were monitored daily for CPE for 7 days.

(e) Model construction and assumptions
First, we used a logistic regression model to estimate the EIP10, e.g. the EIP measured as the time until 10% of infected mosquitoes become infectious for a given temperature, based on our experimental data available at 17°C, 19°C, 24°C, 27°C and 31°C. See section 'Models of extrinsic incubation period (EIP10)' in the electronic supplementary material for further details.
Second, we include this newly obtained EIP10(T°) scheme for ZIKV and Ae. albopictus into a basic reproduction number R 0 model for ZIKV (one vector-one host model). Given recently published evidence, the model has been slightly modified.
R 0 (T ) is given by The biting (a), mortality (µ) and extrinsic incubation (ν = 1/ EIP10(T°)) rates depend dynamically on temperature. All other epidemiological parameters are fixed to a constant value (electronic supplementary material, table S1). A rainfall criterion is applied to mask desertic regions following [22]. The final R 0 (T ) estimate is standardized to range between 0 and 1, to be consistent with previously published models and make our results directly comparable with such studies [23]. We also calculate the potential length of the ZIKV transmission season (LTS in months) based on standardized R 0 (T ). For further details about the standardized R 0 (T ) model and LTS, see sections 'Derivation of R 0 (T ) for ZIKV-Ae. albopictus model' and 'Derivation of the length of the ZIKV transmission season (LTS)-Ae. albopictus model' in the electronic supplementary material.
Finally, we integrate the R 0 (T ) model with observed gridded climate data for the recent period and with an ensemble of climate change projections for the future to investigate regions at risk of ZIKV transmission by Ae. albopictus. See section 'Observed climate datasets and climate change scenarios input data' in the electronic supplementary material for further details.

Results
Zika RNA was detected in the saliva of both Ae. albopictus (Verano colony) and Oc. detritus at all temperatures from 19°C to 31°C, but not at 17°C (figure 1). The presence of ZIKV in the saliva of Ae. albopictus and Oc. detritus was recorded from or before day 17, with the earliest time point being day 7 at 31°C for Ae. albopictus (figure 1a) and day 10 at 27°C and 31°C for Oc. detritus (figure 1b). The titre of ZIKV in the saliva of Ae. albopictus is 3.8× higher than in the saliva of Oc. detritus ( p < 0.00001) across all saliva-positive individuals (electronic supplementary material, tables S2). The bodies of all females were positive at day 0 (immediately after feeding) demonstrating that virus was ingested. On day 5, all bodies were negative for virus demonstrating the eclipse phase of virus infection with undetectably low titres. From day 7 onwards virus was detected in the bodies of both species (for numbers of mosquito bodies and saliva tested at each temperature and time point see electronic supplementary material, table S3). Mortality of Ae. albopictus and Oc. detritus was similar at all temperatures up to 24°C; above this, the mortality rate of the more temperate Oc. detritus was significantly higher (electronic supplementary material, figure S1). The negative competence results for smaller numbers of tested Cx. pipiens pipiens and Cs. annulata can be found in the electronic supplementary material, table S4. Due to small numbers they are not discussed further in the main text.
In order to demonstrate that ZIKV in the saliva of the mosquitoes incubated at the lowest positive temperature (19°C) is viable, we tested the ability of virus extracted from saliva to infect VERO cells via a cytopathic effect (CPE) assay. Saliva from 30 G 11 Ae. albopictus (Verano colony) females was extracted 21 days after ZIKV challenge and were pooled into 10 lots of three samples. Four of these 10 pools were positive for viable virus, indicating at least 13% of saliva samples were positive (consistent with our results for Ae. albopictus at day 21 and 19°C; see figure 1a). This confirmed that our minimum temperature for demonstrable transmission does produce infectious virus and hence does present a transmission risk. Pictures of these cells can be found in the supplementary materials (see electronic supplementary material, table S5 and figure S2).
Using our experimental data, table 1 gives maximumlikelihood estimates and associated 95% confidence intervals for EIP 10 , e.g. the time in days to reach 10% saliva-positive for Ae. albopictus and Oc. detritus.
Model 1, which represents the estimated proportion of mosquitoes susceptible to saliva infection, includes a linear effect of temperature and separate estimates of infection plateaux (δs) for the two species. Model 1 was identified as the preferred model (see electronic supplementary material, text, tables S1 and S6, and figure for further details about model design, selection and validation).
The employed logistic regression model (Model 1) identified a negative linear relationship between the duration of the EIP and temperature for both Ae. albopictus and Oc. detritus within the temperature ranges used in this study (electronic supplementary material, figure S3b). Additionally, the plateau value representing the proportion of mosquitoes susceptible to saliva infection (δs) for Ae. albopictus (0.27) was significantly greater than that of Oc. detritus (0.20) (electronic supplementary material, table S7 and figure S4), indicating that Ae. albopictus is more competent for ZIKV at all temperatures. The resulting standardized R 0 (T) function peaks around 29.2°C (electronic supplementary material, figure S3d). Our standardized R 0 (T) peak estimate for Ae. albopictus potential to transmit ZIKV (electronic supplementary material, figure S3d, solid red line) occurs at slightly higher temperature than highlighted previously for Ae. aegypti [24] and is very similar to a published estimate [25] (electronic supplementary material, figure S3d, solid black line). Note that if we fitted an exponential scheme instead of a linear scheme for our EIP(T) estimates, this does not significantly impact the shape of the resulting R 0 (T) function as mortality (U-shape curve, electronic supplementary material, figure S3c) tends to determine the upper and lower suitability boundaries for standardized R 0 (T ) (see electronic supplementary material, figure S5). Based on our laboratory experiments, the minimum temperature threshold for ZIKV transmission lies between 17°C and 19°C (figure 1a); this temperature range corresponds to standardized R 0 (T ) ranging between 0.201 and 0.295; these values are highlighted in beige on figure 2 and electronic supplementary material, figure S3d. Above 19°C, we know that ZIKV transmission by Ae. albopictus occurs in the laboratory (figure 1a); these values are highlighted in orange, red and dark red on figure 2 and electronic supplementary material, figure S3d. In the following, we focus on regions where annual standardized R 0 (T ) exceeds these critical thresholds, using the same colour code and we also discuss potential changes in the simulated length of the ZIKV transmission season (LTS thereafter) at global scale.
Annual standardized R 0 (T ) maps reveal that ZIKV transmission risk is largest over the tropics but a moderate risk also extends into temperate regions of North America, Europe and Asia      royalsocietypublishing.org/journal/rspb Proc. R. Soc. B 287: 20200119 scenario for the 2080s, a large increase in standardized R 0 is simulated over the south of France, the eastern coast of Spain, Italy, the Adriatic coasts, and parts of Hungary, Serbia, Croatia, Moldova, Romania and Bulgaria (figure 3h). The ZIKV transmission season could extend to six to seven months in future over a large area covering the south of France, Italy, the Balkans and regions surrounding the Black Sea (electronic supplementary material, figure S8). Interestingly, no ZIKV transmission could occur over the UK based on recent climate conditions (electronic supplementary material, figure S6e), but short seasonal transmission (one to three months) might occur over southern UK during the 2080s based on the highest emission scenario (electronic supplementary material, figure S8 h). ZIKV transmission could also extend farther north in the USA, such as into Nebraska, West Virginia and Pennsylvania (figure 3; electronic supplementary material, figure S9). In temperate Asia, the southeastern half of China and Yunnan province, Taiwan, southern Japan, and to a lesser extent the coasts of South Korea could be at risk of ZIKV transmission by temperate mosquitoes by the 2080s (figure 3; electronic supplementary material, figure S10). LTS could range between five and six months over temperate regions in Asia by the 2080s (electronic supplementary material, figure S10 h). While Ae. albopictus is not believed to be the primary vector of ZIKV in Africa or South America, there is a clear increase in simulated standardized R 0 over high altitude regions ( plateaux of Ethiopia, western Kenya, plateaux of Madagascar, see figure 3; electronic supplementary material, figure S11). A large part of South America is still simulated to be climatically suitable for the transmission of ZIKV in future (figure 3; electronic supplementary material, figure S12). However, there is a decrease in simulated R 0 (T) (figure 3) and a shortening of the ZIKV transmission season over Brazil in future for the 2080s (electronic supplementary material, figure S12 h), due to high temperatures.

Discussion
In this study, we assessed the competence of common species from both wild populations of Oc. detritus (UK), and newly colonized Ae. albopictus (Verano colony, Italy) and found that they are competent to transmit viable ZIKV at temperatures as low as 19°C. We also presented the most comprehensive analysis of the effect of temperature on ZIKV incubation to date, establishing that the minimum temperature for ZIKV transmission is between 17°C and 19°C in two different competent temperate vectors, and provide estimates of the time required for 10% to become infectious (as measured by detection of virus in saliva in infected females). Using these estimates, we modelled risk across temperate regions from Ae. albopictus, both in current climatic conditions and in future projected conditions. We used the recently isolated PE243 strain of ZIKV from Brazil, and upper estimates of field-representative titres of ZIKV for mosquito infections [32]. This integration of laboratory studies on field-derived or recently field-derived mosquitoes with modelling makes this study highly applicable to the ongoing risk of ZIKV spread.
Our data show that Ae. albopictus is a more competent laboratory vector than Oc. detritus, in terms of the proportion of infectious individuals and the titre of the virus in the saliva, being significantly higher in Ae. albopictus. However, there is a considerable proportion of infectious Oc. detritus, suggesting that there may be other temperate species which are also highly competent for ZIKV and present a transmission risk (for a detailed list of competent vectors for ZIKV, excluding Oc. detritus, see [10]).
Our finding of a minimum temperature for transmission of ZIKV in both Ae. albopictus and Oc. detritus of between 17°C and 19°C is supported by the findings of Heitmann et al. [33] showing ZIKV transmission in Ae. albopictus at 27°C but none at 18°C after 21 days. Interestingly, this minimum temperature is substantially higher than the minimum recorded temperature for transmission of the closely related WNV at 14°C [13], indicating that ZIKV may not present as much risk to temperate regions as WNV. ZIKV appears to develop at lower temperatures than DENV, for which most studies have shown no or excessively slow viral development at temperatures lower than 20°C [34], indicating that ZIKV presents a greater risk to temperate regions. Previous studies of ZIKV competence by Ae. albopictus have almost exclusively used only 'standard' tropical insectary temperatures of 27 ± 1°C or with a single additional lower temperature [33]. To our knowledge, no other study has examined ZIKV development in a temperate-adapted vector at multiple temperatures or investigated the effect of temperatures as low as 17°C, giving our work unique applicability to temperate regions. Indeed, it is widely believed that climate change may lead to the spread of many vectors and vector-borne diseases to more extreme latitudes and to higher altitudes [23], with strong evidence that in some cases this has already happened, driven by the changing climate [35,36]. Predicting extensions of vector-borne disease distributions into cooler, more temperate regions requires studies at these low temperatures, rather than extrapolation from higher temperature studies, so that the lower temperature threshold can be determined with an acceptable accuracy. Disease model parameterization often relies on extrapolation and our investigation of mosquito infections at low temperatures allows us to better estimate transmission potential, a critical factor to understanding the impact of climate change on vector-borne diseases [15].
Other studies have investigated the importance of wildcaught temperate vectors. Gendernalik et al. [37] demonstrated 80% competence in a common worldwide vector Ae. vexans (caught in the USA) incubated at 27°C. Studies on wild-caught Ae. japonicus highlight the importance of assessing the EIP at temperate-realistic temperatures, as Jansen et al., 2018 [38] showed competence at 27°C but very little at 21°C. Combined, these studies show significant potential risk from temperate vectors, but highlight the need to assess transmission at lower temperatures to determine realistic risk in these regions.
Our R 0 risk maps indicate a significant current risk over large areas of southeast USA. These areas spread considerably further than the limit of where autochthonous transmission of ZIKV has occurred to date, which has been restricted almost exclusively to Florida and, to a lesser extent, Texas [27,28]. Aedes aegypti is thought to be the main ZIKV vector in these regions, however, our results suggest that Ae. albopictus could theoretically sustain ZIKV transmission in these states. Furthermore, our model's prediction of large risk of Ae. albopictus transmission of ZIKV over southern provinces of China, southern Japan, and Taiwan is supported by autochthonous royalsocietypublishing.org/journal/rspb Proc. R. Soc. B 287: 20200119 transmission of dengue virus. These data are reported on a yearly basis for the past 25 years over Guangdong, Guangxi and Hunan provinces, with the largest outbreak reported in 2014 [39], and are consistent with former published estimates [25]. This strongly supports concerns about the risk posed by travellers coming from ZIKV endemic countries to these southern provinces [40]. In contrast to the maps of USA and temperate Asia, our risk maps shows that despite Ae. albopictus presence in large areas of Europe [41], the recent climatic conditions do not result in large R 0 values. This suggests that the range of this important vector extends beyond the minimum temperature limit of ZIKV and consequently the presence of vector alone is not a good indicator of present day risk of sustained Zika transmission. However, seasonal transmission of ZIKV by Ae. albopictus is simulated over southern Europe. This finding is consistent with recently reported autochthonous cases of Zika in Hyères, near to Marseille in southeastern France [42]. We show that ZIKV transmission risk could increase significantly in the future over southern and eastern Europe, extend farther north into the USA and into temperate regions of Asia; and decrease over South America and West Africa, if we consider the highest emission scenario for the 2080s. ZIKV transmission by Ae. albopictus could increase over high-altitude regions; this finding is similar to expected impacts of climate change on malaria [44]. Our results are also consistent with other studies that show likely trends in future [23,45]; however our modelling framework also allows us to explore the impact of rainfall changes on simulated R 0 . Aedes albopictus can flourish in urban settings, using man-made containers to lay its eggs. Consequently, rainfall might not be a limiting factor in densely populated areas or on irrigated lands.
We acknowledge certain limitations in our experimental design, as follows.
(1) Our study focuses on the risk in temperate regions and hence does not account for transmission by the primary ZIKV vector Ae. aegypti, which is important in transmission in some of the warmer areas of regions we investigated. We don't consider this a limitation per se, but rather a focus on different regions. (2) Our analysis is based solely on a temperate strain of Ae.
albopictus; our finding that Oc. detritus is also a highly competent temperate vector, as well as other published data showing other competent vectors (reviewed in [10]), suggests that there may be a risk of transmission in temperate regions outside of the range of Ae. albopictus.
For example, Oc. detritus larvae mature in brackish water, often in estuaries or expansive wetland areas [46,47] that are very different to the small volume containers, tyres or natural pools that are preferred by Ae. albopictus and Ae. aegypti. Locally, they can both achieve very high densities and are voracious human feeders; such species could present a major additional vector control challenge in the event of a ZIKV outbreak.
(3) Our experimental work focuses on fixed temperature incubation, even though diurnal temperature variability has been shown to be important for many of the parameters in the R 0 estimation [48,49]. The breadth of different temperature cycling regimens that would need to be tested make such work unfeasible in a study of this scale. Further studies isolating a single average temperature and testing multiple temperature regimens would be needed to assess this effect fully. (4) Our work demonstrates that the minimum temperature for transmission of viable ZIKV is between 17°C and 19°C. Future experiments, including EIP and virus infectivity at these temperatures, are planned to narrow this range at low temperatures and therefore produce more accurate estimations of the absolute limit of risk.
Our estimates of the potential length of the ZIKV transmission season by Ae. albopictus match recently published findings. In our model, seasonal ZIKV transmission could occur during four to five months over southern Europe; this finding is consistent with another R 0 model estimate (about four months transmission shown in [45]). The pattern and magnitude of the simulated length of the ZIKV transmission season over northern and southern America is also consistent with [16,23], who also showed potential year-round ZIKV transmission over Central America, the Caribbean and the northern half of South America. Simulated LTS (about seven to eight months) over the southern provinces of China (Guangdong and Guangxi) also corresponds to observed seasonal transmission of dengue fever over these regions [50,51]. Our model only considers the impact of climate on the mosquito vector potential to transmit ZIKV. Such a model cannot account for other important factors including the number of infected travellers arriving, and the immunological history of the local population [52,53]. Herd immunity is an important factor in risk and is expected to be low in naive European, North American and temperate Asian populations. Conversely, antibody-dependent enhancement (ADE) of infection, most commonly associated with DENV, may enhance ZIKV viremia and therefore potentially increase the likelihood of mosquito acquisition of the virus [54]. Such ADE may have been a significant driving force in non-flavivirus-naive populations, but is likely to be reduced in temperate populations. The interplay between climatic and other external factors will likely have a large effect on the risk to a population. Using field applicable experimental data, we present a comprehensive analysis of the relationship between ZIKV risk and temperature. Our results show that the risk of autochthonous ZIKV transmission extends beyond the regions in which it has been reported thus far, especially in Asia and the USA. They also show in Europe that vectors can be present in regions that are too cold for sustained ZIKV transmission. The comprehensive analysis of EIP and minimum temperature requirements presented here offers a route to enable informed risk management and outbreak preparedness in more specific temperate situations.