- Open Access
Climatic effects on mosquito abundance in Mediterranean wetlands
Parasites & Vectorsvolume 7, Article number: 333 (2014)
The impact of climate change on vector-borne diseases is highly controversial. One of the principal points of debate is whether or not climate influences mosquito abundance, a key factor in disease transmission.
To test this hypothesis, we analysed ten years of data (2003–2012) from biweekly surveys to assess inter-annual and seasonal relationships between the abundance of seven mosquito species known to be pathogen vectors (West Nile virus, Usutu virus, dirofilariasis and Plasmodium sp.) and several climatic variables in two wetlands in SW Spain.
Within-season abundance patterns were related to climatic variables (i.e. temperature, rainfall, tide heights, relative humidity and photoperiod) that varied according to the mosquito species in question. Rainfall during winter months was positively related to Culex pipiens and Ochlerotatus detritus annual abundances. Annual maximum temperatures were non-linearly related to annual Cx. pipiens abundance, while annual mean temperatures were positively related to annual Ochlerotatus caspius abundance. Finally, we modelled shifts in mosquito abundances using the A2 and B2 temperature and rainfall climate change scenarios for the period 2011–2100. While Oc. caspius, an important anthropophilic species, may increase in abundance, no changes are expected for Cx. pipiens or the salt-marsh mosquito Oc. detritus.
Our results highlight that the effects of climate are species-specific, place-specific and non-linear and that linear approaches will therefore overestimate the effect of climate change on mosquito abundances at high temperatures. Climate warming does not necessarily lead to an increase in mosquito abundance in natural Mediterranean wetlands and will affect, above all, species such as Oc. caspius whose numbers are not closely linked to rainfall and are influenced, rather, by local tidal patterns and temperatures. The final impact of changes in vector abundance on disease frequency will depend on the direct and indirect effects of climate and other parameters related to pathogen amplification and spillover on humans and other vertebrates.
Vector-borne diseases (VBD) represent a serious threat to human health and wildlife conservation. Over the last three decades and concurrent to the emergence of newly recognized diseases, the incidence and geographic distribution of VBD have increased dramatically and now account for almost one third of all cases of emerging diseases. Although VBD transmission is determined by many factors, from host immunity to pathogen circulation, in recent years great emphasis has been placed on the notion that climate change may alter patterns of VBD transmission[5–10]. Climate influences various aspects of the vector-pathogen-host system and the existence of a close relationship between climate and the force of infection (i.e. the rate at which a population of hosts is infected) has been hypothesized. Currently, our understanding of the underlying mechanisms that influence mosquito-borne disease transmission cycles is not complete and the potential impact of global warming on these diseases is today a subject of on-going debate and a source of controversy[7, 12–15].
Climate has direct effects on vectors (e.g. abundance, distribution and longevity), pathogens (e.g. incubation period, replication and lineage), hosts (e.g. abundance, distribution and behaviour) and their interactions. Mosquito abundance is an important factor that conditions vectorial capacity and Ro (the basic reproductive rate), and high abundance is often a prelude to an epidemic[7, 10]. Previous studies have concluded that high temperatures and high rainfall are positively related to mosquito abundance; even so, in some cases mosquito abundance peaks are known to be preceded by periods of drought (Additional file1).
West Nile virus (WNV), Chikungunya, Dengue and Usutu are considered to be (re)-emerging in Europe[9, 16–18]. However, only limited information regarding the influence of climatic variables on the population dynamics of vectors in Europe is available[19–21]. In particular, WNV has become widespread[22–25] and disease outbreaks have occurred when epizootic vectors transmit this virus from birds to horses or humans. Wetlands are linked to endemic WNV circulation and here we study two wetlands in SW Spain where: 1) WNV circulation since at least 2003 has been documented regularly in resident birds[28, 29], horses and humans; 2) there is a great abundance of several mosquito species potentially involved in disease transmission[32, 33] (Table 1); and 3) outbreaks of WNV in humans and horses have been reported from a neighbouring province[34–36]. In addition, the mosquito-borne flavivirus Bagaza has been detected in the surrounding area, there is a high incidence of Dirofilaria immitis and the potential malarial vector Anopheles atroparvus is abundant (Table 1).
Consequently, we investigated the relationships between several climatic variables and the female abundance of seven mosquito species on inter-annual and seasonal scales in the period 2003–2012. We also modelled the potential effects of climate change on mosquito abundance using downscaled General Circulation Model outputs under two of the emission scenarios postulated by the Spanish National Meteorological Agency (AEMET).
The studied locations lie within the Doñana National Park and Odiel Natural Park in SW Spain (Figure 1). These two sites are included on the Ramsar list of wetlands of international importance and are classified as Biosphere reserves by UNESCO. The climate is Mediterranean subtropical in type, with hot dry summers and frequent rainfall in autumn and winter. The mean, maximum and minimum temperatures reported are 18.3, 31.6 and 7.7°C, respectively, while mean annual rainfall is 516 mm (maximum monthly mean rainfall of 83 mm in December and minimum 1 mm in June). Doñana is a strongly seasonal freshwater wetland, while Odiel is very influenced by tides and has brackish wetlands. The most representative plants in the seasonal freshwater coastal areas are Scirpus maritimus L., S. littoralis Schrader, Typha latifolia L., and Phragmites australis (Cav.) Trin., while in tidal salt-marshes the vegetation is mainly composed of halophyte species such as Spartina densiflora Brongn, Salicornia ramossisima J. Woods, Sarcocornia perennis (Miller) A. J. Scott and Arthrocnemum macrostachyum (Moric.). In the Odiel salt-marshes operational mosquito control with the larvicide Bacillus thuringiensis var. israeliensis (Bti) (VectoBac 12 AS, Valent Biosciences, Libertyville, IL, USA) is carried out on a fortnightly basis every year from March to October. The avian communities of these wetlands are characterized by both migratory and sedentary species of ducks, herons, gulls, waders and other birds, with over 400 species and several millions of individuals recorded. Wild and domestic mammals (horses and cows) are also common, mainly in Doñana, and more than four million people visit the area annually, with a daily influx of more than 10,000 tourists during spring and summer.
Mosquito sampling and data
CDC (Centre for Disease Control) traps with 6-volt batteries were supplemented with a dry-ice tank of 1 kg capacity. Traps were located 2 m above ground level and were operated in 24-hour cycles. Trapping took place twice a week in February–November in 2003–2012. Traps were placed at three locations in areas of natural habitat for mosquitoes: Los Álamos (hereafter Álamos), Celestino Mutis (hereafter Mutis) and Palacio de Doñana (hereafter Palacio). Alamos and Mutis are near the tidal marshes of the Tinto and Odiel rivers, whilst Palacio is surrounded by the freshwater marshes of the Guadalquivir River (Figure 1). All three areas are environmentally protected and no landscape changes occurred in the area surrounding the traps during the study period. Mosquitoes were killed and transported to the laboratory on dry ice. Mosquitoes were counted, sexed and identified using taxonomic keys[38, 39] with a stereomicroscope and a chill table.
Temperature (mean, maximum, minimum), rainfall, relative humidity (mean, maximum and minimum), wind speed and direction, solar radiation, evapotranspiration and daylight hours during the study period were extracted from data from the closest meteorological stations to the trapping localities: Las Torres, Tomejil Moguer (http://www.mapa.es/siar/Informacion.asp) and Palacio de Doñana (http://www-rbd.ebd.csic.es/mediofisico/parametrosmeteorologicos/palaciomanual/em05.htm). Tidal data were obtained from the Anuarios de Mareas del Instituto Hidrográfico de la Marina for the ports of Mazagón (Huelva) and Bonanza (Cadiz) at the mouth of the Guadalquivir River. Data were pooled by week and the mean weekly values were computed for all variables except for rainfall, for which the total weekly precipitation was used (Figure 2).
Statistical analysis of intra-annual (seasonal) patterns
To control for temporal autocorrelation, we fitted an autoregressive time series model to the weekly female mosquito abundance data for each mosquito species, selecting the complexity with a model selection procedure based on Akaike’s Information Criterion (AIC). With the order of the autoregressive time series, we fitted an ARIMA model to each female mosquito abundance time series. We extracted the residuals of the ARIMA that were used as the response variable for the intra-annual models. The explanatory variables were mean, maximum and minimum temperature, tide height, mean and total rainfall, wind speed, wind direction, solar radiation, evapotranspiration and photoperiod (daylight hours). Given that a time lag operates for the effects of climatic variables on larval abundance (and therefore on adult mosquito abundance), it is important to use aggregated climatic variables to detect biological relationships. Therefore, a series of accumulated temperature, rainfall and tide variables were calculated for 1–4 weeks before the sampling week. The variables used in the analysis were the accumulated average temperature, accumulated rainfall and accumulated tide of 1–2, 1–3, 1–4, 2–3, 2–4 and 3–4 weeks before the sampling week. When the explanatory variables are highly correlated to each other, backward selection can select suboptimal models due to collinearity. To select a set of unrelated independent variables, we compared the Variance Inflation Factor (VIF) of covariates, then excluded the covariate with the highest VIFs, recalculated the correlation and repeated the process until all VIF values were lower than four. Based on this exploratory analysis, we selected the following explanatory variables: mean temperature during the week of capture, accumulated temperatures 1–4 weeks before capture, the accumulated tide heights two weeks before capture, accumulated tide heights 2–3 weeks before capture, mean relative humidity, wind speed, sum of rainfall during the week of capture, accumulated rainfall one week before capture, and accumulated rainfall 1–2 and 2–4 weeks before capture. Accumulated tide variables were only analysed in relation to the salt-marsh mosquitoes Ochlerotatus caspius and Oc. detritus, two species that use larval habitats that depend on tidal patterns. Given that the photoperiod is an important variable with a close collinearity to temperature, we performed a GLMM analysis for all species that only included photoperiod. For Cx. theileri, Cx. modestus and An. atroparvus we only used data from Palacio since in the two other localities these species were absent or only present in very low numbers. To model relationships between female mosquito abundance (residuals of the temporal autocorrelation model) and climatic covariates we performed a GLMM with Gaussian error. We fitted the models by pseudo-likelihood and the minimal models were chosen using AICs for each variables and a forward stepwise model selection procedure. All the independent variables were log10-transformed. Statistical analyses were performed using the R statistical package version 2.13.1 (2011 The R Foundation for Statistical computing).
Statistical analysis of inter-annual patterns
The effects of environmental variables on the estimates of annual average female mosquito abundances (2003–2012) in the commonest captured species were analyzed using Generalized Additive Models (GAMs). Given that the count data presented a significant overdispersion, the model was fitted using a negative binomial distributed error and a logarithm link. The average number of female mosquitoes was used as the response variable for each of the seven commonest species. The explanatory variables were as follows: average annual maximum temperature, annual mean temperature, average annual minimum temperature, annual rainfall and winter rainfall. Winter rainfall was calculated as the accumulated rainfall between September (week 40) and March of the following year (week 10). These variables were modelled including ‘locality’ as an independent factor and assuming that mean values for abundance may differ spatially. Model selection was performed using AICs for each variable and using a forward stepwise model selection procedure; validation was performed based on Zuur. The inter-annual analysis was only performed for Cx. pipiens, Oc. caspius and Oc. detritus because these mosquito species were the only ones regularly captured in all three localities. Statistical analyses were performed with the R statistical package version 2.13.1 (2011 The R Foundation for Statistical computing) using packages mgcv, lattice and MuMin, among others.
Climate change scenarios
The Spanish National Meteorological Agency (AEMET) has developed climate projections for the whole of the twenty-first century (http://www.aemet.es/es/elclima/cambio_climat/escenarios). We decided to use the ‘Analogue INM’ method used by the AEMET given that this statistical downscaling method, which presents the estimated data for temperature and rainfall in area E009 (Doñana) in the high-resolution reference grid for peninsular Spain and the Balearic Islands, provides the finest spatial resolution. This method provided us with the expected temperature (minimum and maximum) and rainfall values for the Doñana area using two General Circulation Models (GCMs), ECHAM (Max Planck Institute for Meteorology, Germany) and CGCM (Canadian Centre for Climate Modelling and Analysis), under the A2 and B2 emission scenarios. Scenario A2 is more extreme, and assumes self-reliance, the preservation of local identities, moderate economic development and a high growth in the global population; although energy consumption is high and changes in land use are moderately high, resources become scarce and technological change is fragmented but slower than under other scenarios. On the other hand, scenario B2 is more moderate than A2 and emphasizes environmental preservation and social equity with local solutions for economic, social and environmental sustainability; the global population is expected to increase continuously, yet more slowly than under scenario A2. B2 allows for a moderate level of economic development (like A2), but with lower energy consumption and fewer changes in land use; resources are more abundant and technological change is more diverse than under A2.
Based on the inter-annual models, we projected inter-annual variation in modelled female mosquito abundance based on the AEMET modelling and calculated the percentage shift in projected female mosquito abundance in the period 2011–2100 in relation to the reference period 1961–1990. Analyses were carried out using R software.
In total, 184 326 female mosquitoes belonging to seven species were captured (Table 2): Anopheles atroparvus Van Thiel, Culex pipiens Linnaeus, Culex theileri Theobald, Culex modestus Ficalbi, Culex perexiguus Theobald, Ochlerotatus caspius (Pallas) and Ochlerotatus detritus (Haliday). Oc. caspius was the most abundant species throughout the whole season, despite displaying a certain seasonality in its cycles (Figure 3). Cx. pipiens and Oc. detritus were abundant in all three localities in spring–summer and spring–autumn, respectively. Cx. theileri, Cx. modestus, Cx. perexiguus and An. atroparvus were especially abundant at Palacio (Table 2 and Figure 4).
Intra-annual relationships of climate and female mosquito abundance
Mean temperature was positively related to Cx. pipiens and Oc. caspius abundances (Table 3). Accumulated temperatures in the period 1–4 weeks (one month) before capture were positively related to Cx. perexiguus and Cx. modestus, and negatively to Cx. pipiens and Oc. detritus abundances. Accumulated rainfall during the period 2–4 weeks before capture was positively related to Oc. detritus abundance. Accumulated tide heights for two weeks before capture were positively related to Oc. caspius abundance and accumulated tide heights during 2–3 weeks before capture to Oc. detritus abundance. Mean relative humidity was negatively related to Cx. theileri and An. atroparvus abundances (Table 3). None of the other variables were significant (Additional file2). The photoperiod (hours of light) was positively related to Cx. perexiguus, Cx. theileri, An. atroparvus and Oc. caspius abundances but was negatively related to Oc. detritus abundance (Table 4).
Interannual relationship of climatic variables and female mosquito abundance
Cx. pipiens abundance was positively related to annual maximum temperatures and winter rainfall. Annual mean temperatures were positively related to Oc. caspius (Table 5; Figure 5), while Oc. detritus abundance was positively related to winter rainfall.
Climate change scenarios
Using the increase in temperatures and decrease in rainfall in the study area forecasted by the scenarios (Figure 6), we projected variations in the abundances of female Cx. pipiens, Oc. caspius and Oc. detritus. The other mosquito species were not represented because their inter-annual analysis did not give any significant variable. Under climate change scenarios A2 and B2 for the period 2013–2100, it is predicted that Oc. caspius abundance will considerably increase but that Oc. detritus and Cx. pipiens abundances will not change (Figure 7). Although projections do not change qualitatively between scenarios, larger quantitative changes in mosquito populations are expected under the A2 scenarios.
Mosquito abundance is related to climatic variables, although, depending on the species in question and the temporal scale used, the key variables will vary. Temperatures in general were positively related to mosquito abundance and can affect survival rates, mortality, larval productivity and other population parameters[7, 10]. Weekly temperatures were related to seasonal abundance patterns in Cx. pipiens and Oc. caspius, while accumulated temperatures (1–4 weeks before) were positively related to Cx. modestus and Cx. perexiguus abundances. These results are consistent with previous studies in USA, South Africa and Russia that link high summer temperatures to great vector abundance and WNV outbreaks[20, 47–49]. On the other hand, accumulated temperatures (1–4 weeks before) were negatively related to Cx. pipiens and Oc. detritus abundances. These negative relationships, like the negative effect of annual maximum temperatures over 40°C on Cx. pipiens abundance (Figure 5), may reflect the negative effect of high temperatures on mosquito longevity and life expectancy. These two mosquito species may be more sensitive to high temperatures, although an indirect effect may also exist given that high evapotranspiration and the absence of larval habitats coincide with the period of highest temperatures.
Accumulated rainfall (2–4 weeks before) is positively related to Oc. detritus weekly abundance. We did not detect any relationship between drought periods and mosquito abundance as has been reported in other studies[51–53]. However, both Cx. theileri and An. atroparvus within-year abundances are negatively related to relative humidity, which could be due to drought, to an effect on flight activity, or to high evapotranspiration. In general, the relationship rainfall-mosquito abundance is variable and depends on differences in the ecology and habitat selection patterns of each mosquito species.
The influence of tide heights on the abundance of salt-marsh species has previously been reported by: both Oc. caspius and Oc. detritus weekly abundances were related to tide height, albeit with a time lag of 1–2 weeks before capture in the former and of 2–3 weeks in the latter. This lag corresponds to the fortnightly seasonal rhythms of these species, which benefit from the high tides that flood the saltmarshes and create larval breeding sites. We did not detect any effect of wind on within-year abundance patterns. Wind may negatively affect captures by reducing flight activity of mosquitoes and it is important to remember that the number of captures in mosquito traps depends not only on mosquito abundance but also on mosquito activity.
Photoperiod (day-length) was positively related to the within-year abundance patterns of Cx. perexiguus, Cx. theileri, Oc. caspius and An. atroparvus, but negatively to Oc. detritus. Day-length is an important variable in spring emergence and autumn decline, especially in species that enter into diapause, and probably affects the seasonal dynamics of mosquito abundance in several species. Unlike temperature or rainfall, the seasonal pattern in day-length does not vary from year to year and provides a highly reliable anticipatory cue for future seasonal conditions. However, few studies have ever considered the photoperiod as an important variable affecting seasonal dynamics in mosquito populations. Overall, our results agree with several previous studies (Additional file1) that conclude that higher temperatures and greater accumulated temperatures, as well previous (winter) rainfall, are all important factors that help explain mosquito abundance.
Annual temperatures were related non-linearly to the inter-annual patterns of mosquito abundance in Cx. pipiens (maximum temperature) and Oc. caspius (mean temperatures). The influence of annual rainfall on mosquito abundance is a controversial question (Additional file1); the results of our study show that previous rainfall is an important variable that is positively related to Cx. pipiens and Oc. detritus annual abundances.
When analysing long-term series of mosquito abundance data it is important to consider the effects that changes in larvicide treatments may have had. In our case it is unlikely that our results were affected by changes in pesticide treatment given that no such treatment is used within 45 km of Palacio. At the other two sites, however, the tidal areas suitable for Oc. caspius and Oc. detritus larval development were treated with larvicides, but at a distance of over 1 km from the traps. The product and amount used did not change over time and treatment was applied in a fortnightly cycle after high tides. This periodicity would have hindered rather than favoured our ability to detect any relationship between tides and population size within seasons and probably had no effect on inter-annual analysis or climate change modelling. Additionally, no landscape changes that could have affected the mosquito populations were detected during the study period.
The Spanish National Meteorological Agency AEMET projections suggest that climate change will have particularly marked effects in southern Spain and will lead to higher temperatures and lower rainfall. Given the projections under A2 and B2 climate change scenarios and assuming the existence of a causal relationship between the climatic variables we studied and female mosquito abundance, we do not expect that any important shift (increase or decrease) in either Cx. pipiens (a WNV and USUV vector) or the salt-marsh mosquito Oc. detritus abundances will occur. The already abundant and anthropophilic pest mosquito Oc. caspius, a potential vector of diseases such as Tahyna, WNV and tularemia, could greatly increase in abundance, especially under the A2 climate change scenario (Figure 7), that is, more extreme estimates of change than under B2. Oc. caspius has a great resistance to high temperatures and drought, is active, bites at high temperatures (up to 36°C in adults and 33°C in larvae), and tolerates high salinities due to its great capacity for osmotic regulation. This species depends heavily on areas flooded with brackish water in estuaries during high tides in the Atlantic and, consequently, its populations are less rainfall dependent than other species. In Doñana, we have described its bloodmeal behaviour as mammalophilic, since more than 90% of its blood meals are derived from mammals, but it has little importance in the amplification of WNV in the area. In addition, in recent years the novel flavivirus Marisma Mosquito Virus has been detected in Oc. caspius pools analysed in the area, although the potential pathogenicity (if any) of this virus for vertebrates has yet to be demonstrated. We expect that the nuisance caused by Oc. caspius to local populations will increase in the future, although its implications for disease dynamics are unclear. Oc. detritus and Oc. caspius depend on tide patterns that vary in terms of sea levels and so, hypothetically, rising sea levels under global warming could lead to an increase in salt-marsh mosquito abundance. Projections of seasonal changes in rainfall patterns suggest a decrease in spring and summer precipitations and an increase in the length of drought periods. Consequently, we would also expect changes in the phenology of mosquito abundance. The reliability of these projections depends on the nature of the relationship between climate and mosquito abundance, which will have to be confirmed with manipulative experiments in addition to the observational evidence presented here since the relationship between climate and life-history processes are not likely to be linear . The next challenge is to accurately understand the contribution of multiple interacting and often non-linear underlying responses in hosts, pathogens and vectors to climate.
When we reviewed the studies analysing correlations between climatic variables and WNV human cases and infection rates, the same variables (temperature and rainfall) were related to vector mosquito abundance, infection rates and cases of infection in humans (Additional file1). Moreover, mosquito abundance has been directly related to the basic reproductive rate Ro (the total number of secondary cases arising from one infective case in a susceptible population) and the vectorial capacity C (the daily rate at which secondary cases arise from a currently infective case)[7, 10]. Additionally, temperature alters the replication of the virus in the mosquito – the extrinsic incubation period (EID) – and therefore infection rates[48, 65]. In this work, we focus on how climate change could influence transmission by altering mosquito vector abundance. An important question, given that models for vector-borne diseases use the number of vectors per host, is whether the number of vectors caught in a trap is related to the number of vectors per host. Come what may, the effects of temperature are complex due to the fact that the climate-mosquito abundance relationship is species-specific and time-dependent, and long data series are needed to continue the study of these trends in mosquito abundance. Our results highlight the fact that these effects are non-linear (see Figure 5, Cx. pipiens and maximum temperatures) and that, therefore, linear approaches will overestimate the effect of climate change on mosquito abundances at high temperatures. Without additional studies our results are not generalizable to other habitats since relationships may vary between areas and may be place-specific[49, 51]. Although it is likely that anthropogenic changes in modified wetlands will be more determinant than climate, this study was conducted in natural wetlands unaffected by any serious anthropogenic change such as population increase or agricultural and urban expansion. Thus, our results show that climate change does not necessarily lead to an increase in mosquito populations. Simplistic and generalist statements assuming that higher temperatures lead to more mosquitoes are not adequate and it is essential to carry out a careful analysis of temporal patterns in field data. In addition, climate is likely to affect other important parameters for disease transmission (such as vector survival and pathogen development rates) and, consequently, it is not possible to extrapolate from conclusions regarding climate and mosquito abundances to the risk of disease outbreaks. Other non-climatic variables (i.e. socio-economic development, circulation of pathogenic strains, pathogen replication, host ecology, anthropogenic changes, land use/land-cover and inundation patterns) may also be determinant when analysing the risk of vector-borne disease emergence. Greater understanding of the ecology of vector-borne disease is essential if we are to understand the effects of global change on vector-borne disease outbreaks, an important challenge for future decades.
Gratz NG: Emerging and resurging vector-borne diseases. Annu Rev Entomol. 1999, 44: 51-75.
Daszak P, Cunningham AA, Hyatt AD: Anthropogenic environmental change and the emergence of infectious diseases in wildlife. Acta Trop. 2001, 78 (2): 103-116.
Gubler DJ: The Global Threat of Emergent/Re-emergent Vector-Borne Diseases. Vector Biology, Ecology and Control. Edited by: Atkinson PW. 2010, Dordrecht: Springer Netherlands, 39-62.
Taylor LH, Latham SM, Woolhouse MEJ: Risk factors for human disease emergence. Philos T Roy Soc B. 2001, 356 (1411): 983-989.
Shope R: Global climate change and infectious diseases. Environ Health Persp. 1991, 96: 171-174.
Epstein PR: Global warming and vector-borne disease. Lancet. 1998, 351 (9117): 1737-
Reiter P: Climate change and mosquito-borne disease. Environ Health Persp. 2001, 109 (Suppl 1): 141-161.
Patz JA, Campbell-Lendrum D, Holloway T, Foley JA: Impact of regional climate change on human health. Nature. 2005, 438 (7066): 310-317.
Gould EA, Higgs S: Impact of climate change and other factors on emerging arbovirus diseases. Trans Roy Soc Trop Med Hyg. 2009, 103 (2): 109-121.
Lafferty KD: The ecology of climate change and infectious diseases. Ecology. 2009, 90: 888-900.
Dobson A: Climate variability, global change, immunity, and the dynamics of infectious diseases. Ecology. 2009, 90: 920-927.
Brisbois BW, Ali SH: Climate change, vector-borne disease and interdisciplinary research: social science perspectives on an environment and health controversy. EcoHealth. 2010, 7 (4): 425-438.
Rogers DJ, Randolph SE, Simon I, Hay AG, David JR: Climate Change and Vector-Borne Diseases. Global Mapping of Infectious Diseases: Methods, Examples and Emerging Applications. 2006, London: Academic Press, 345-381.
Rohr JR, Dobson AP, Johnson PTJ, Kilpatrick AM, Paull SH, Raffel TR, Ruiz-Moreno D, Thomas MB: Frontiers in climate change–disease research. Trends Ecol Evol. 2011, 26 (6): 270-277.
Patz JA, Graczyk TK, Geller N, Vittor AY: Effects of environmental change on emerging parasitic diseases. Int J Parasitol. 2000, 30 (12–13): 1395-1405.
Rezza G, Nicoletti L, Angelini R, Romi R, Finarelli AC, Panning M, Cordioli P, Fortuna C, Boros S, Magurano F: Infection with chikungunya virus in Italy: an outbreak in a temperate region. Lancet. 2007, 370 (9602): 1840-1846.
Vazquez A, Jimenez-Clavero M, Franco L, Donoso-Mantke O, Sambri V, Niedrig M, Zeller H, Tenorio A: Usutu virus: potential risk of human disease in Europe. Eurosurveillance. 2011, 16 (31): 19935-
Lindgren E, Andersson Y, Suk JE, Sudre B, Semenza JC: Monitoring EU emerging infectious disease risk due to climate change. Science. 2012, 336 (6080): 418-419.
Savage HM, Ceianu C, Nicolescu G, Karabatsos N, Lanciotti R, Vladimirescu A, Laiv L, Ungureanu A, Romanca C, Tsai TF: Entomologic and avian investigations of an epidemic of West Nile fever in Romania in 1996, with serologic and molecular characterization of a virus isolate from mosquitoes. Am J Trop Med Hyg. 1999, 61 (4): 600-611.
Platonov AE, Shipulin GA, Shipulina OY, Tyutyunnik EN, Frolochkina TI, Lanciotti RS, Yazyshina S, Platonova OV, Obukhov IL, Zhukov AN: Outbreak of West Nile virus infection, Volgograd Region, Russia, 1999. Emerg Infect Dis. 2001, 7 (1): 128-132.
Ponçon N, Toty C, L’ambert G, Le Goff G, Brengues C, Schaffner F, Fontenille D: Population dynamics of pest mosquitoes and potential malaria and West Nile virus vectors in relation to climatic factors and human activities in the Camargue, France. Med Vet Entomol. 2007, 21 (4): 350-357.
Hubálek Z, Halouzka J: West Nile fever–a reemerging mosquito-borne viral disease in Europe. Emerg Infect Dis. 1999, 5 (5): 643-650.
Calistri P, Monaco F, Savini G, Guercio A, Purpari G, Vicari D, Cascio S, Lelli R: Further spread of West Nile virus in Italy. Vet Ital. 2010, 46 (4): 467-474.
Angelini P, Tamba M, Finarelli AC, Bellini R, Albieri A, Bonilauri P, Cavrini F, Dottori M, Gaibani P, Martini E: West Nile virus circulation in Emilia-Romagna, Italy: the integrated surveillance system 2009. Eurosurveilancel. 2010, 15 (16): 19547-
Danis K, Papa A, Papanikolaou E, Dougas G, Terzaki I, Baka A, Vrioni G, Kapsimali V, Tsakris A, Kansouzidou A: Ongoing outbreak of West Nile virus infection in humans, Greece, July to August 2011. Eurosurveillance. 2011, 16 (34): 19951-
Kramer LD, Styer LM, Ebel GD: A global perspective on the epidemiology of West Nile virus. Annu Rev Entomol. 2008, 53 (1): 61-81.
Jourdain E, Gauthier-Clerc M, Bicout D, Sabatier P: Bird migration routes and risk for pathogen dispersion into Western Mediterranean Wetlands. Emerg Infect Dis. 2007, 13 (3): 365-372.
Figuerola J, Jimenez-Clavero MA, Rojo G, Gomez-Tejedor C, Soriguer R: Prevalence of West Nile virus neutralizing antibodies in colonial aquatic birds in southern Spain. Avian Pathol. 2007, 36 (3): 209-212.
Figuerola J, Jimenez-Clavero MA, Lopez G, Rubio C, Soriguer R, Gomez-Tejedor C, Tenorio A: Size matters: West Nile virus neutralizing antibodies in resident and migratory birds in Spain. Vet Microbiol. 2008, 132 (1–2): 39-46.
Jimenez-Clavero MA, Llorente F, Sotelo E, Soriguer R, Gomez-Tejedor C, Figuerola J: West Nile virus serosurveillance in horses in Donana, Spain, 2005 to 2008. Vet Rec. 2010, 167 (10): 379-380.
Bernabeu-Wittel M, Ruiz-Pérez M, del Toro MD, Aznar J, Muniain Á, de Ory F, Domingo C, Pachón J: West Nile virus past infections in the general population of Southern Spain. Enferm Infec Micr Cl. 2007, 25 (9): 561-
Muñoz J, Ruiz S, Soriguer R, Alcaide M, Viana DS, Roiz D, Vázquez A, Figuerola J: Feeding patterns of potential West Nile virus vectors in South-West Spain. PLoS ONE. 2012, 7 (6): e39549-
Vazquez A, Sanchez-Seco MP, Ruiz S, Molero F, Hernandez L, Moreno J, Magallanes A, Tejedor CG, Tenorio A: Putative new lineage of west nile virus, Spain. Emerg Infect Dis. 2010, 16 (3): 549-552.
Sotelo E, Fernández-Pinero J, Jiménez-Clavero MÁ: La fiebre/encefalitis por virus West Nile: reemergencia en Europa y situación en España. Enferm Infec Micr Cl. 2012, 30 (2): 75-83.
García-Bocanegra I, Jaén-Téllez JA, Napp S, Arenas-Montes A, Fernández-Morente M, Fernández-Molera V, Arenas A: West Nile fever outbreak in horses and humans, Spain, 2010. Emerg Infect Dis. 2011, 17 (12): 2397-2399.
Sotelo E, Gutierrez-Guzmán A, del Amo J, Llorente F, El-Harrak M, Pérez-Ramírez E, Blanco J, Höfle U, Jiménez-Clavero M: Pathogenicity of two recent Western Mediterranean West Nile virus isolates in a wild bird species indigenous to Southern Europe: the red-legged partridge. Vet Res. 2011, 42: 11-
Aguero M, Fernández-Pinero J, Buitrago D, Sánchez A, Elizalde M, San Miguel E, Villalba R, Llorente F, Jiménez-Clavero MA: Bagaza virus in partridges and pheasants, Spain, 2010. Emerg Infect Dis. 2011, 17: 1498-1501.
Becker N, Petric D, Zgomba M, Boase C, Madon M, Dahl C, Kaiser A: Mosquitoes and Their Control. 2010, Berlin: Springer
Schaffner F, Angel G, Geoffroy B, Hervy JP, Rhaiem A, Brunhes J: Les moustiques d’Europe/The mosquitoes of Europe. CD-ROM Montpellier. 2001, France: Institut de Recherche pour le Développement/EID Méditerranée
Shone SM, Ferrao PN, Lesser CR, Norris DE, Glass GE: Analysis of mosquito vector species abundances in Maryland using geographic information systems. Ann NY Acad Sci. 2006, 951: 364-368.
Zuur AF, Ieno EN, Smith GM: Analysing Ecological Data. 2007, Berlin: Springer
Zuur AF, Ieno EN, Elphick CS: A protocol for data exploration to avoid common statistical problems. Methods Ecol Evol. 2010, 1 (1): 3-14.
Bolker BM, Brooks ME, Clark CJ, Geange SW, Poulsen JR, Stevens MHH, White J: Generalized linear mixed models: a practical guide for ecology and evolution. Trends Ecol Evol. 2009, 24 (3): 127-135.
Crawley MJ: Statistical Computing: an Introduction to Data Analysis using S-Plus. 2002, Chichester: John Wiley and Sons
Ver Hoef JM, Boveng PL: Quassi-poisson vs. negative binomial regression: how should we model overdispersed count data?. Ecology. 2007, 88: 2766-2772.
Zuur AF: A Beginner’s Guide to Generalized Additive Models with R. 2014, Newburg: Highland Statistics Ltd
Jupp PG: The ecology of West Nile virus in South Africa and the occurrence of outbreaks in humans. Ann NY Acad Sci. 2001, 951 (1): 143-152.
Reisen WK, Fang Y, Martinez VM: Effects of temperature on the transmission of west nile virus by Culex tarsalis (Diptera: Culicidae). J Med Entomol. 2006, 43 (2): 309-317.
Ruiz MO, Chaves LF, Hamer GL, Sun T, Brown WM, Walker ED, Haramis L, Goldberg TL, Kitron UD: Local impact of temperature and precipitation on West Nile virus infection in Culex species mosquitoes in northeast Illinois, USA. Parasit Vectors. 2010, 3 (1): 19-
Tabachnick WJ: Challenges in predicting climate and environmental effects on vector-borne disease episystems in a changing world. J Exp Biol. 2010, 213 (6): 946-954.
Landesman WJ, Allan BF, Langerhans RB, Knight TM, Chase JM: Inter-annual associations between precipitation and human incidence of West Nile virus in the United States. Vector-Borne Zoonot. 2007, 7: 337-343.
Shaman J, Day JF, Stieglitz M: Drought-induced amplification and epidemic transmission of West Nile virus in Southern Florida. J Med Entomol. 2005, 42: 134-141.
Uejio CK, Kemp A, Comrie AC: Climatic controls on West Nile virus and Sindbis virus transmission and outbreaks in South Africa. Vector-Borne Zoonot. 2011, 12 (2): 117-125.
Leisnham PT, Sandoval-Mohapatra S: Mosquitoes associated with ditch-plugged and control tidal salt marshes on the Delmarva Peninsula. Int J Environ Res Public Health. 2011, 8 (8): 3099-3113.
Hribar LJ, DeMay DJ, Lund UJ: The association between meteorological variables and the abundance of Aedes taeniorhynchus in the Florida Keys. J Vector Ecol. 2010, 35 (2): 339-346.
Matowo NS, Moore J, Mapua S, Madumia EP, Moshi IR, Kaindoa EW, Mwangungulu SP, Kavishe DE, Sumaye RD, Lwetoijera DW, Okumu FO: Using a new odour-baited device to explore options for luring and killing outdoor-biting malaria vectors: a report on design and field evaluation of the Mosquito Landing Box. Parasites and Vectors. 2013, 6: 137-
Bradshaw WE, Holzapfel CM: Climate change. Evolutionary response to rapid climate change. Science. 2006, 312 (5779): 1477-1478.
Gong H, DeGaetano A, Harrington L: Climate-based models for West Nile Culex mosquito vectors in the Northeastern US. Int J Biometeorol. 2011, 55 (3): 435-446.
Brunet M, Casado MJ, de Castro M, Galán P, López JA, Martín JM, Pastor A, Petisco E, Ramos P, Ribalaygua J: Generación de escenarios de cambio climático regionalizados para España. 2009, Madrid: Agencia Estatal de Meteorología, Ministerio de Medio Ambiente, 158-
Bradley TJ: Physiology of osmoregulation in mosquitoes. Annu Rev Entomol. 1987, 32: 439-462.
Vázquez A, Sánchez-Seco M-P, Palacios G, Molero F, Reyes N, Ruiz S, Aranda C, Marqués E, Escosa R, Moreno J, Figuerola J, Tenorio A: Novel flaviviruses detected in different species of mosquitoes in Spain. Vector-Borne Zoonot. 2012, 12 (3): 223-229.
Ramasamy R, Surendran SN: Global climate change and its potential impact on disease transmission by salinity-tolerant mosquito vectors in coastal zones. Front Physiol. 2012, 3: 198-
de Castro M, Martín-Vide J, Alonso S: El clima de Espana: Pasado, Presente Y Escenarios de Clima para el Siglo XXI. Evaluación Preliminar de los Impactos en España por Efecto del Cambio Climático. Edited by: Moreno-Rodríguez JM. 2005, Madrid: Ministerio de Medio Ambiente, 1-64.
Altizer S, Ostfeld RS, Johnson PTJ, Kutz S, Harvell CD: Climate change and infectious diseases: from evidence to a predictive framework. Science. 2013, 341 (6145): 514-519.
Kilpatrick AM, Meola MA, Moudy RM, Kramer LD: Temperature, viral genetics, and the transmission of West Nile virus by Culex pipiens mosquitoes. PLoS Pathog. 2008, 4 (6): e1000092-
Baylis M: Research gaps in understanding how climate change will affect arboviral diseases. Anim Health Res Rev. 2013, 14 (02): 143-146.
Ponçon N, Tran A, Toty C, Luty AJF, Fontenille D: A quantitative risk assessment approach for mosquito-borne diseases: malaria re-emergence in southern France. Malaria J. 2008, 7 (1): 147-
Kilpatrick AM: Globalization, land use, and the invasion of West Nile virus. Science. 2011, 334 (6054): 323-327.
We would like to thank Juani Moreno and Antonio Magallanes from the Mosquito Control Service, Huelva, Spain, as well as Esmeralda Perez Morueta, Isidro Román and Pilar Bayon from the Estación Biológica de Doñana, Seville, Spain, for their help and technical support. This work was partially supported by projects P07-RNM-02511, RNM118 and RNM157 of the Junta de Andalucía, Plan Nacional of the Spanish Ministry of Science (CGL2009-11445), and the European Commission EDEN FP6, EuroWestnile EU grant HEALTH.2010.2.3.3-3 Project 261391 and EDENext FP7 project.
The authors declare that they have no competing interests.
SR designed and participated in field sampling; DR, SR and JF organised the analyses, results and discussion; DR analysed the data; JF and RS provided financial support. All the authors were involved in the drafting of the manuscript and gave their approval to the final version.
Electronic supplementary material
Additional file 2: Results of the models (GLMM gaussian) for all the climatic variables that are related with seasonal patterns (intra-annual) of female abundance for seven mosquito species and the estimates of the final model selected by backward selection. Significant variables are in bold. The significance (F, p) of the non-significant variables corresponds to the final model plus the non-significant variable. (DOC 118 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.