Determinants of the population growth of the West Nile virus mosquito vector Culex pipiens in a repeatedly affected area in Italy
- Paolo Mulatti^{1},
- Heather M Ferguson^{2}Email author,
- Lebana Bonfanti^{1},
- Fabrizio Montarsi^{1},
- Gioia Capelli^{1} and
- Stefano Marangon^{1}
https://doi.org/10.1186/1756-3305-7-26
© Mulatti et al.; licensee BioMed Central Ltd. 2014
Received: 21 August 2013
Accepted: 8 January 2014
Published: 15 January 2014
Abstract
Background
The recent spread of West Nile Virus in temperate countries has raised concern. Predicting the likelihood of transmission is crucial to ascertain the threat to Public and Veterinary Health. However, accurate models of West Nile Virus (WNV) expansion in Europe may be hampered by limited understanding of the population dynamics of their primary mosquito vectors and their response to environmental changes.
Methods
We used data collected in north-eastern Italy (2009–2011) to analyze the determinants of the population growth rate of the primary WNV vector Culex pipiens. A series of alternative growth models were fitted to longitudinal data on mosquito abundance to evaluate the strength of evidence for regulation by intrinsic density-dependent and/or extrinsic environmental factors. Model-averaging algorithms were then used to estimate the relative importance of intrinsic and extrinsic variables in describing the variations of per-capita growth rates.
Results
Results indicate a much greater contribution of density-dependence in regulating vector population growth rates than of any environmental factor on its own. Analysis of an average model of Cx. pipiens growth revealed that the most significant predictors of their population dynamics was the length of daylight, estimated population size and temperature conditions in the 15 day period prior to sampling. Other extrinsic variables (including measures of precipitation, number of rainy days, and humidity) had only a minor influence on Cx. pipiens growth rates.
Conclusions
These results indicate the need to incorporate density dependence in combination with key environmental factors for robust prediction of Cx. pipiens population expansion and WNV transmission risk. We hypothesize that detailed analysis of the determinants of mosquito vector growth rate as conducted here can help identify when and where an increase in vector population size and associated WNV transmission risk should be expected.
Keywords
Culex pipiens Culicidae Mosquito population dynamics Density dependence West nile virus Disease vectorsBackground
In the last decade, changes in global climate have raised concern as being a potential trigger for the expansion of infectious diseases, and in particular those borne by arthropod vectors (Vector Borne Diseases – VBDs), into previously unaffected areas [1]. The emergence and recurrence of VBDs such as Blue Tongue, West Nile Disease, and Chikungunya into areas of Europe formerly thought not to be environmentally suitable for their transmission have been considered as a confirmation of this threat within temperate zones [2–5]. Nevertheless, climate change alone may not represent a sufficient cause to explain the behaviour or the expansion of VBDs in newly affected areas [6]. Other factors, including socio-economic and land-use factors, may nullify or counter any impacts of climate change in enhancing conditions for VBD establishment and transmission (e.g. as has been hypothesized for malaria in Camargue and in south-eastern England [6]). However, other areas have reported an expansion of VBDs in areas where environmental conditions were not thought to favour transmission, such as West Nile in USA [6]. The variable response of VBD transmission to environmental change [6, 7] indicates that prediction of the risk of emergence and spread into new areas is complex, and likely requires detailed knowledge of the specific pathogen-vector system and their environmental dependencies. We hypothesize that a thorough understanding of the core ecological determinants of the population dynamics of arthropod vectors is fundamentally required to improve predictions of the potential extent and rate of VBD expansion and risk.
West Nile virus (WNV) has become one of the most widely distributed arboviruses in the world [8]. The virus is normally maintained and spread through a bird-mosquito cycle involving ornithophilic Culex spp. mosquitoes [8, 9]. However, the virus can also be spread to a wide range of incidental hosts including humans and horses, via mammophilic and/or anthropophilic mosquito species (including Aedes spp. and Ochlerotatus spp.), and it has the potential to cause severe illness characterised by neurological disorders [10]. The presence of WNV has also been sporadically recorded in hard (Ixodus spp. and Amblyomma spp.) and soft ticks (Argas spp.), with some cases of transmission to other hosts [9]. Various non-vectorial modes of transmission have also been observed in a variety of mammals, including humans, suggesting the possibility that a WNV cycle could be sustained in the absence of mosquito vectors [11]. Nevertheless, non-vectorial routes, and also non-mosquito vectors, appear to play a very minor role in the maintenance of the WNV cycle in nature [12].
West Nile Virus has been circulating in Europe for more than 30 years, primarily in central and Mediterranean countries, with cases being recurrently reported both in humans and in horses [13]. The first appearance of WNV in Italy occurred in 1998 in Tuscany, with 14 infected horses and no human cases [14]. No further clinical manifestations of WNV were reported until 2008, when 251 horse premises tested positive for WNV in the Po Valley of north-eastern Italy [15]. Since then WNV has been recurrently detected in horses, wild birds, and humans [16, 17], suggesting the likely endemisation of the disease [18]. In 2010, WNV was also identified in Culex pipiens mosquitoes in the neighbouring regions of Veneto and Emilia Romagna. Since then the virus has been isolated from Cx. pipiens only, suggesting these mosquitoes are the primary WNV vectors in Italy [19].
Despite the widespread distribution of Cx. pipiens, understanding of the determinants of their population dynamics both within Italian WNV transmission areas, and other parts of Europe is rather limited. Previously it has been assumed that environmental factors such as temperature and rainfall are the overwhelming drivers of mosquito population growth [20]. However, recent investigations of other tropical and temperate mosquito species indicate that endogenous density dependence may have a strong influence on their population growth, as arising primarily from intra-specific competition for resources within their aquatic larval habitats [21–23]. The primary implication of this phenomenon is that intra-specific competition may limit the size of vector populations at levels below which would otherwise be expected from abiotic conditions alone. As such, existence of strong density-dependent regulation within vector populations may limit the potential for WNV transmission to below that estimated from environmental conditions on their own.
Here the existence of such a phenomenon within the putative main WNV vector in Italy was investigated for the first time. Using longitudinal data obtained over three years from a WNV surveillance programme implemented in north-eastern Italy, the seasonal abundance of Cx. pipiens in north-eastern Italy was analyzed to quantify the relative contribution of density dependence and exogenous environmental factors to their population growth. This information will provide a more realistic understanding of the constraints on WNV vector populations within Europe, and facilitate prediction of their spatial and/or temporal expansion and associated epidemiological impact. Ultimately, such detailed knowledge of vector population determinants could provide a more reliable framework for development of risk-based WNV surveillance and early detection programs.
Methods
Mosquito collection
Only Cx. pipiens has been considered in the analyses, since it has been acknowledged as being the most important WNV vector in Italy [19]. Moreover, as the two subspecies belonging to Cx. pipiens (Cx. pipiens pipiens and Cx. pipiens molestus) are morphologically indistinguishable, the analyses were conducted on Cx pipiens complex as a whole.
Environmental monitoring
Environmental data were recorded at weather stations located in the closest proximity to mosquito traps (mean distance: 5,095 m, 95% CI 4,507 – 5,683 m). Collected data included daily measurements of the maximum, minimum and mean temperatures, minimum, maximum, and average relative humidity rates, and cumulative daily precipitation.
where T_{ mean } is the daily average temperature, and L.T. and U.T. are the lower and upper temperature thresholds respectively. In the present study, values of 13°C and 33°C were used for these lower and upper thresholds respectively [30, 31]. The daily GDDs were then summed to obtain the cumulative GDD for the 15 days preceding the mosquito capture.
Culex pipiens mosquitoes can undergo a diapause, defined as a period when development and activity is suppressed [32]. Entrance into the diapause stage is strongly dependent on photoperiod, with the proportion of mosquitoes in diapause rapidly decreasing as the length of daylight increases [32]. Consequently the length of daylight was recorded for each capture day to test its association with mosquito population dynamics (data collected from http://aa.usno.navy.mil/).
Environmental variables used in the generalized linear mixed-effects models
Variable | Description |
---|---|
PREC | Daily cumulative precipitation (mm) |
GDD | Growing degree days |
HMN | Daily minimum relative humidity (%) |
PREC.15d | Cumulative precipitation in 15 days prior to capture(mm) |
DPREC.15d | Number of rainy days in 15 days prior to capture |
HAV.15d | Average relative humidity in 15 days prior to capture (%) |
GDD.15d | Growing degree days in 15 days prior to capture |
LDM | Length of daylight (minutes) |
Statistical analyses
where N_{ t } is the population density at time t, N_{ t+1 } is the mosquito abundance at the next capture.
where, r_{ m } is the maximum intrinsic growth rate, assumed in our models to be equal to the highest observed growth rate per capture site over the whole observation period, N is the abundance of mosquitoes per capture, K is the carrying capacity (which was assumed to vary for the three study years, and was set as the maximum number of mosquitoes recorded per capture site in each year of observation), and Env is the environmental variable assumed to be a proxy of critical habitat. We assumed that environmental factors likely related to the availability of suitable larval habitat would have the strongest effects on determining carrying capacity. Therefore, we considered the average humidity and cumulative precipitation in the 15 days prior to adult capture as the environmental factors most likely to modify the carrying capacity of mosquito populations [33, 34].
All of the statistical models were fitted using Maximum Likelihood mixed-effects linear regression (Generalized Linear Mixed-effects Model, GLMM) [35]. Grouping variables related to the year of observation (2009, 2010 or 2011) and site of capture were considered as crossed random effects. For all of the 15 alternative models the Akaike’s information criterion corrected for small sample size was calculated, and the strength of support was evaluated through calculation of their Akaike weights (w AIC_{c} ), which may be interpreted as the probability that a model is the most likely [36]. For each of the tested GLMMs, a likelihood-ratio based pseudo-R^{2} (R^{ 2 }_{ LR }) was calculated as a measure of goodness of fit [37, 38]. The model with the strongest degree of statistical support was selected as the most appropriate base model for further investigations of the role of additional environmental factors on mosquito populations.
After the best base model for population growth was identified, additional model selection procedures were conducted to assess the further contribution of other environmental factors. Model selection was based on Information Theoretic methodologies based on the corrected Akaike Information Criterion (IT-AIC_{c}) [39–41]. A full initial model was created by including all of the environmental variables considered (Table 1), and a set of models containing the terms included in base population model and subsets of terms of the global model was generated. All of the models were compared through w AIC_{c} to assess whether there was a single model with high statistical support, or many models with similar values of w AIC_{c}. In the case of uncertainty in model selection, a model averaging algorithm was conducted which allowed combination of the parameter estimates from a selected set of models, considering the contribution of each model as proportional to its likelihood weight [39, 41]. The set of models whose total cumulative w AIC_{c} was at least equal to 0.95 were selected for model averaging. In doing this, it can be inferred that the selected set of models is expected to contain the AIC-best model with a probability of 0.95 if further data were added [41]. The importance of each single variable in the model with averaged coefficients, were calculated as the cumulative w AIC_{c} for the models in which the variables were included as predictors. This value, may also be interpreted as the probability that a particular variable is included in the best AIC_{c} model [41].
The generalized mixed-effect models were fitted using the nlme package [42], the IT-AIC_{c} approach was performed through the package MuMIn, in the R statistical software version 2.15.1 [43].
Results
Annual trends in mosquito distribution and abundance
Number of mosquitoes collected during the sampling periods 2009-2011
2009 | 2010 | 2011* | |
---|---|---|---|
Sampling period (first capture – last capture) | 27 May – 11 Nov | 3 May – 26 Oct | 3 May – 25 Oct |
Total collected mosquitoes | 35, 129 | 137, 897 | 85, 136 |
Total number of captures | 222 | 536 | 733 |
Avg. no. mosquitoes per capture | 158. 24 (S.E.: 17. 77) | 257. 27 (S.E.: 24. 64) | 116. 14 (S.E.: 7.85) |
Total collected Cx. pipiens | 27, 721 (78. 91%) | 119, 847 (86. 91%) | 68, 932 (80. 97%) |
Avg. no. Cx. pipiens per capture | 124.87 (S.E.: 15. 01) | 223. 60 (S.E.: 22. 40) | 94. 04 (S.E.: 6. 74) |
Total no. of identified species | 12 | 16 | 14 |
Predictors of mosquito population growth rate
Statistical models built to define the baseline population dynamics of Cx. pipiens
Model | ΔAIC_{c} | w AIC_{c} | R ^{ 2 } _{ LR } |
---|---|---|---|
Gompertz | 0.00 | > 0.999 | 0.59 |
Gompertz mod. by humidity | 33.26 | < 0.001 | 0.58 |
Ricker | 192.99 | < 0.001 | 0.52 |
Ricker mod. by humidity | 195.03 | < 0.001 | 0.52 |
Growing degree-day (15 dd) | 209.23 | < 0.001 | 0.52 |
Growing degree-day | 334.21 | < 0.001 | 0.47 |
Gompertz mod. by by precipitation | 354.26 | < 0.001 | 0.46 |
Rainy days (15 dd) | 355.79 | < 0.001 | 0.46 |
Precipitation (15 dd) | 357.21 | < 0.001 | 0.46 |
Reference model | 359.02 | < 0.001 | 0.46 |
Exponential | 359.02 | < 0.001 | 0.46 |
Ricker mod. by precipitation | 360.05 | < 0.001 | 0.46 |
Minimum daily humidity | 360.62 | < 0.001 | 0.46 |
Precipitation | 360.77 | < 0.001 | 0.46 |
Average daily humidity (15 dd) | 361.02 | < 0.001 | 0.46 |
Both Ricker- and Gompertz-logistic models performed better than any single density-independent environmental factor. Within alternative density-dependent models, those in which the mosquito carrying capacity was modified by proxy indices of habitat availability (i.e. average humidity and cumulative precipitation in the 15 days prior to mosquito capture), had lower explanatory power than base Ricker- and Gompertz-logistic alternatives (Table 3). Overall, the Gompertz-logistic model provided a substantially better representation of the seasonal growth dynamics of Cx. pipiens populations than any of the other base models considered, with an overwhelming w AIC_{c} support greater than 0.999 (Table 3). Only 6 of the 15 alternative growth models resulted explaining more variation than the reference model: the base Gompertz- and Ricker-logistic models (R^{ 2 }_{ LR } = 0.59, and R^{ 2 }_{ LR } = 0.55, respectively), their version with carrying capacity modified by humidity (R^{ 2 }_{ LR } =0.58 and R^{ 2 }_{ LR } = 0.55, for the modified Gompertz and Ricker models respectively), and the density-independent models with GDD and the cumulative GDD in 15 days as predictor (R^{ 2 }_{ LR } = 0.47, and R^{ 2 }_{ LR } = 0.52, respectively).
Coefficients of the fixed effects for average model for Cx. pipiens population dynamics
Coefficient | Estimate | 95% C.I. | S.E. | Importance |
---|---|---|---|---|
(Intercept) | -0.234 | -0.270; -0.197 | 0.019 | -- |
LDM | 0.789 | 0.740; 0.839 | 0.025 | 1.00 |
Ln(N)/Ln(K) | -0.421 | -0.471; -0.371 | 0.026 | 1.00 |
GDD | 0.029 | -0.007; 0.064 | 0.018 | 0.57 |
GDD.15d | -0.106 | -0.151; -0.061 | 0.023 | 1.00 |
PREC | 0.004 | -0.023; 0.031 | 0.014 | 0.25 |
PREC.15d | 0.002 | -0.032; 0.036 | 0.017 | 0.26 |
DPREC.15d | -0.021 | -0.055; 0.013 | 0.018 | 0.42 |
HMN | 0.000 | -0.028; 0.028 | 0.014 | 0.24 |
HAV.15d | -0.014 | -0.050; 0.021 | 0.018 | 0.32 |
Comparisons between the relative effects of the different factors were facilitated by scaling and centering the environmental variables before model fitting. Therefore, it was possible to make inferences simply considering the magnitude of coefficients, with higher absolute values indicating stronger effects on the population growth rates [44], despite the environmental factors being recorded with different units of measure. The highest absolute value for the coefficients was associated with the length of daylight. This association with daylight is hypothesized to be driven by the known influence that the shortening of photoperiod has on determining the beginning of the diapause phase in Cx. pipiens[32]. The positive sign of the coefficient indicates that mosquito reproduction capacity was higher when daylight lasted longer (i.e. in summer), while it dropped with a decrease of the photoperiod. The form of density-dependence acting on Cx. pipiens was predicted to be strongly negative, indicating that when the population density approaches the carrying capacity, their growth rate decreases. The impact of the cumulative GDD was also found to be negative (Table 4), which could be indicative of negative effects occurring over prolonged periods where temperatures occur near the maximum tolerance for Cx. pipiens, and/or may have been related to limited ability of larvae to acclimate above specific temperature thresholds (e.g. 25°C for Cx. pipiens) [45]. Alternatively, this apparently negative impact of GDD may be a product of its correlation with another environmental variable not measured here that has a directly negative association on GDD (e.g. predator abundance, etc). Further research into this potential relationship with GDD is required to clarify its impact on mosquito population growth.
Overall, the combination of length of daylight, density-dependence, and cumulative GDD in the 15 days prior to capture generated the best description of Cx. pipiens growth dynamics. All three of these variables were included in the top 48 models (out of a set of 128) that accounted for a cumulative w AIC_{c} ≤0.95. In contrast, the other environmental factors considered (daily GDD, daily precipitation, cumulative precipitation and number of rainy days within 15 days prior to mosquito collections, and minimum daily humidity and average humidity in 15 days before mosquito captures) were of relatively low significance (Table 4), as interpreted from their lower probability of being included in the best AIC_{c} model [41].
Discussion
West Nile Virus mosquito vector species ecology and dynamics are commonly considered to depend mostly on external environmental factors, such as rainfall or temperature [20, 46, 47]. While all mosquito populations are likely regulated by such environmental factors to some extent, recent evidence from tropical species indicates that density-dependence can also play a large role in determining population growth rates [21, 22]. The existence of density-dependence has not been widely investigated in temperate mosquito species, nor in vectors of WNV, therefore, we developed and tested a series of statistical models, to investigate the relative importance of density-dependent versus independent factors in predicting the population growth rates of the WNV vector Cx. pipiens, in a recurrently affected area of Italy.
Our analyses showed that density-dependent models provide a much better description of Cx. pipiens growth rates during the study period than any extrinsic environmental variable on its own. This suggests that mosquito density has a more significant effect on regulating the dynamics of Cx. pipiens than previously thought. Moreover, little or no interactions between environmental factors we assumed to be strictly correlated to the habitat of mosquitoes (average humidity rate and cumulative precipitation, both related to the environmental situation within 15 days prior to mosquito captures) and the mosquito carrying capacity resulted in our study. This indicates that these climatic factors may be not critical determinants of the carrying capacity during the surveillance season (May-November) in the study area.
The inclusion of further environmental factors to the base density-dependent model (Gompertz-logistic) through the IT-AIC_{c} approach allowed definition of the relative importance of each of both density-dependent factors and all of the environmental variables. The length of the daylight, the density-dependent factor included in the Gompertz-logistic model, and the cumulative GDD resulted in being the most important factor variables in describing variation in Cx. pipiens per capita growth rates. In contrast, the other environmental predictors considered here (daily GDD, daily precipitation, cumulative precipitation and number of rainy days within 15 days prior to mosquito collections, and minimum daily humidity and average humidity in 15 days before mosquito captures) were judged to have relatively less impact of mosquito population dynamics, on their limited importance as defined through the IT-AIC_{c} approach. We caution that these findings should not be interpreted as dismissing the role of abiotic factors on Cx. pipiens dynamics, but as evidence that the upper limit of mosquito population growth rates may be set just as much or more by intra-specific competition than environmental factors.
The IT-AIC_{c} approach used here indicated that a model set with a cumulative w AIC_{c} of 0.95, interpreted as being likely to contain the best AIC-model of mosquito population growth, was achieved from the sum weightings of the highest ranked 48 out of the 128 model set. This can be interpreted as reflecting a high level of uncertainty in defining the single best model (with several performing equally well). The relatively similar performance of several alternative models tested here could arise because of correlations between environmental variables, which would consequently limit the amount of additional information derived from varying the specific combinations of variables tested. However, in the model with averaged coefficients, only the three factors with the higher probability to be included in the best AIC-model resulted in having a significant coefficient. The respectively positive and negative sign of the coefficients for the length of daylight and for the density-dependent factors are easily interpretable: Cx. pipiens population growth is predicted to go up as day length lengthens, and fall as mosquito density rises. However, interpretation of the negative coefficient for the cumulative growth degree-days in 15 days prior to the capture is less clear. This predicted decrease in population growth rate as the number of growing days during the larval development period increases appears counterintuitive, but may be a product of a non-linear relationship between temperature and growth rate; and/or negative correlations between GDD and other key mosquito resources (e.g. evaporation of larval habitat/reduced rainfall associated with higher GDD). This finding appears to contrast with previous studies on WNV mosquitoes in Italy, which described the population density as positively correlated to the temperature [20, 47, 48]. As GDD is a composite measure that incorporates not just temperature but the length of time at conditions deemed optimal for growth, it is conceivable that these measures of environmental conditions have differing implications for population growth. Further study is required to investigate this apparent discrepancy in the nature of temperature effects.
While the modelling investigation here has been useful to highlight the contribution of density-dependent as well as environmental factors to Cx. pipiens population dynamics; we caution that a substantial amount of variation in mosquito population growth rate remained unexplained by the statistical models tested. Notably, the overall goodness of fit of the 48 models used in model averaging was marginally higher than that of the base population dynamics model (which had a rank of 65 out of the 128 tested). Thus further more detailed analyses including data from a longer time series and/or incorporation of additional environmental variables from what was measured here, and their interactions will be required to obtain a model that can predict Cx. pipiens population growth with a high degree of quantitative accuracy; as opposed to revealing the relative strength and direction with which environmental factors and density-dependence influence these dynamics, as done here.
In seeking to estimate the impact of environmental factors in driving this process, a frequently-used approach (and one used here) is to test for the impact of specific environmental factors while controlling for the random effect of site. A possible limitation of this approach is that some environmental factors may be relatively static across time. This could make it difficult to distinguish the relative impact of ‘site’ from ‘site-specific’ environmental factors, and potentially lead to an underestimation of the relative importance of environmental factors considered independently of site. While we acknowledge this as a potential limitation of this approach, we do not feel it could have created a bias in our results as all the environmental variables we considered (rainfall, temperature and humidity-related) varied very extensively within a site over all years, and had little collinearity with ‘site’.
Whilst the statistical models tested here were subject to some limitations, including the assumption of a fixed length of the mosquito life cycle and larval development throughout the study period, our results provide useful insights into the population dynamics of Cx. pipiens in Western Europe. Most notably, the inclusion of density-dependence yielded better prediction of their population growth rates. Consequently, failure to account for density dependence when analysing mosquito population growth rate data may lead to biased inference, especially when the population density becomes remarkably high. Mosquito populations may respond in different ways to fluctuations of environmental conditions according to their abundances. Specifically, at the start of the mosquito growing season when numbers are low, small increases in favourable environmental conditions may lead to rapid increase in population growth. In contrast, when mosquito population density is very high and intra-specific competition high, further improvements in environmental conditions may yield no change as population growth is suppressed by competition [21–23, 49]. Density-dependent models allow for these intrinsic regulatory effects to be incorporated, and thus can provide information on how populations may respond to environmental change that would be missed in studies that consider only the linear effects of environmental factors.
Although the inclusion of density-dependence as shown here can significantly improve our understanding of mosquito population dynamics, we caution that further investigation is required in a wider range of areas that differ in macro-ecological and socio-economic factors to confirm the generality of this conclusion. Furthermore, it could be that wider consideration of other environmental and demographic data that were not included here could yield significantly better predictions of mosquito and population dynamics, and thus reduce the relative importance of density –dependence from what has been described here. However, density dependence is known to play an important role in regulating animal populations in many other ecological systems [50], and we hypothesize it is likely to remain important in this and other mosquito species even in other environmental conditions; albeit its relative importance over climatic factors may be reduced or reversed in areas with more extreme environmental variations in time.
Mosquito species collected within the 2-year surveillance (2009–2011) conducted here included not only Cx. pipiens, but many others that may act as vectors for several VBDs. In addition to WNV, a number of newly introduced arboviruses have been detected in the past five years in Italy. Besides WNV which now appears to be endemic in north-eastern Italy, other viruses including Usutu and Chikungunya viruses have also been isolated from mosquitoes collected in this study, or in adjoining regions [51–53]. These events, together with the introduction of exotic mosquito species that could be related to both climate changes and increased global trade [54], make knowledge of mosquito population dynamics of paramount importance to better understand mechanisms of VBD dissemination and persistence.
Conclusions
Our study provides a first attempt to gain a deeper insight into WNV vector dynamics in Italy. The inclusion of growth rate fluctuations, and of density dependent effects, could enhance models to predict the trend of mosquito populations in areas considered at risk of WNV introduction and spread. Furthermore, the trend of growth rate could give a better idea on when and where to expect an increase in vector population density, and therefore an amplified risk of WNV transmission. Based on this information, appropriate intervention measures such as the implementation of risk-based surveillance programmes could be envisaged, in order to better assess the risk of human infections.
Declarations
Acknowledgments
The research was funded by the Regional Prevention Plan entitled “Surveillance of vector-borne zoonoses and evaluation of the risk of disease spread in Veneto Region” (Regional Committee resolution n. 1094 of 26 July 2011).
The authors wish to thank Paul Johnson and Dan Haydon, University of Glasgow, for their valuable advice in compiling the scripts and interpreting the results of the models used. We wish also to thank Matteo Mazzucato and the Parasitology Laboratory of the IZSVe, for support in managing and elaborating the entomological data analysed in the present study, and the Regional Agencies for Environmental Protection (ARPA) of Veneto and FVG for providing the environmental data.
Authors’ Affiliations
References
- Rogers DJ, Randolph SE: Climate change and vector-borne diseases. Adv Parasitol. 2006, 62: 345-381.View ArticlePubMedGoogle Scholar
- Medlock JM, Hansford KM, Schaffner F, Versteirt V, Hendrickx G, Zeller H, Van Bortel W: A review of the invasive mosquitoes in Europe: ecology, public health risks, and control options. Vector Borne Zoonotic Dis. 2012, 12: 435-447. 10.1089/vbz.2011.0814.PubMed CentralView ArticlePubMedGoogle Scholar
- Bezirtzoglou C, Dekas K, Charvalos E: Climate changes, environment and infection: facts, scenarios and growing awareness from the public health community within Europe. Anaerobe. 2011, 17: 337-340. 10.1016/j.anaerobe.2011.05.016.View ArticlePubMedGoogle Scholar
- Purse BV, Mellor PS, Rogers DJ, Samuel AR, Mertens PPC, Baylis M: Climate change and the recent emergence of bluetongue in Europe. Nat Rev Microbiol. 2005, 3: 171-181. 10.1038/nrmicro1090.View ArticlePubMedGoogle Scholar
- Brugger K, Rubel F: Characterizing the species composition of European Culicoides vectors by means of the Köppen-Geiger climate classification. Parasit Vectors. 2013, 6: 333-10.1186/1756-3305-6-333.PubMed CentralView ArticlePubMedGoogle Scholar
- Randolph SE, Rogers DJ: The arrival, establishment and spread of exotic diseases: patterns and predictions. Nat Rev Microbiol. 2010, 8: 361-371. 10.1038/nrmicro2336.View ArticlePubMedGoogle Scholar
- Randolph SE: To what extent has climate change contributed to the recent epidemiology of tick-borne diseases?. Vet Parasitol. 2010, 167: 92-94. 10.1016/j.vetpar.2009.09.011.View ArticlePubMedGoogle Scholar
- Kramer LD, Styer LM, Ebel GD: A global perspective on the epidemiology of West Nile virus. Annu Rev Entomol. 2008, 53: 61-81. 10.1146/annurev.ento.53.103106.093258.View ArticlePubMedGoogle Scholar
- Hubálek Z, Halouzka J: West Nile fever: a reemerging mosquito-borne viral disease in Europe. Emerg Infect Dis. 1999, 5: 643-650. 10.3201/eid0505.990505.PubMed CentralView ArticlePubMedGoogle Scholar
- Mackenzie JS, Gubler DJ, Petersen LR: Emerging flaviviruses: the spread and resurgence of Japanese encephalitis, West Nile and dengue viruses. Nat Med. 2004, 10 (12 Suppl): S98-S109.View ArticlePubMedGoogle Scholar
- Blitvich BJ: Transmission dynamics and changing epidemiology of West Nile virus. Anim Heal Res Rev. 2008, 9: 71-86. 10.1017/S1466252307001430.View ArticleGoogle Scholar
- Petersen LR, Brault AC, Nasci RS: West Nile virus: review of the literature. JAMA. 2013, 310: 308-315. 10.1001/jama.2013.8042.PubMed CentralView ArticlePubMedGoogle Scholar
- Calistri P, Giovannini A, Hubalek Z, Ionescu A, Monaco F, Savini G, Lelli R: Epidemiology of west nile in europe and in the mediterranean basin. Open Virol J. 2010, 4: 29-37.PubMed CentralPubMedGoogle Scholar
- Autorino GL, Battisti A, Deubel V, Ferrari G, Forletta R, Giovannini A, Lelli R, Murri S, Scicluna MT: West Nile virus epidemic in horses, Tuscany Region, Italy. Emerg Infect Dis. 2002, 8: 1372-1378. 10.3201/eid0812.020234.PubMed CentralView ArticlePubMedGoogle Scholar
- Calistri P, Giovannini A, Savini G, Monaco F, Bonfanti L, Ceolin C, Terregino C, Tamba M, Cordioli P, Lelli R: West Nile virus transmission in 2008 in north-eastern Italy. Zoonoses Public Health. 2010, 57: 211-219. 10.1111/j.1863-2378.2009.01303.x.View ArticlePubMedGoogle Scholar
- Busani L, Capelli G, Cecchinato M, Lorenzetto M, Savini G, Terregino C, Vio P, Bonfanti L, Pozza MD, Marangon S: West Nile virus circulation in Veneto region in 2008–2009. Epidemiol Infect. 2011, 139: 818-825. 10.1017/S0950268810001871.View ArticlePubMedGoogle Scholar
- Savini G, Capelli G, Monaco F, Polci A, Russo F, Di Gennaro A, Marini V, Teodori L, Montarsi F, Pinoni C, Pisciella M, Terregino C, Marangon S, Capua I, Lelli R: Evidence of West Nile virus lineage 2 circulation in Northern Italy. Vet Microbiol. 2012, 158: 267-273. 10.1016/j.vetmic.2012.02.018.View ArticlePubMedGoogle Scholar
- Monaco F, Savini G, Calistri P, Polci A, Pinoni C, Bruno R, Lelli R: West Nile disease epidemic in Italy: first evidence of overwintering in Western Europe?. Res Vet Sci. 2009, 2011 (91): 321-326.Google Scholar
- Romi R, Pontuale G, CIufolini MG, Fiorentini G, Marchi A, Nicoletti L, Cocchi M, Tamburro A: Potential vectors of West Nile virus following an equine disease outbreak in Italy. Med Vet Entomol. 2004, 18: 14-19. 10.1111/j.1365-2915.2004.0478.x.View ArticlePubMedGoogle Scholar
- Bisanzio D, Giacobini M, Bertolotti L, Mosca A, Balbo L, Kitron U, Vazquez-Prokopec GM: Spatio-temporal patterns of distribution of West Nile virus vectors in eastern Piedmont Region, Italy. Parasit Vectors. 2011, 4: 230-10.1186/1756-3305-4-230.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang G-J, Bradshaw CJA, Whelan PI, Brook BW: Importance of endogenous feedback controlling the long-term abundance of tropical mosquito species. Popul Ecol. 2008, 50: 293-305. 10.1007/s10144-008-0082-8.View ArticleGoogle Scholar
- Yang G-J, Brook BW, Whelan PI, Cleland S, Bradshaw CJA: Endogenous and exogenous factors controlling temporal abundance patterns of tropical mosquitoes. Ecol Appl. 2008, 18: 2028-2040. 10.1890/07-1209.1.View ArticlePubMedGoogle Scholar
- Russell TL, Lwetoijera DW, Knols BGJ, Takken W, Killeen GF, Ferguson HM: Linking individual phenotype to density-dependent population growth: the influence of body size on the population dynamics of malaria vectors. Proc Biol Sci. 2011, 278: 3142-3151. 10.1098/rspb.2011.0153.PubMed CentralView ArticlePubMedGoogle Scholar
- Roiz D, Roussel M, Muñoz J, Ruiz S, Soriguer R, Figuerola J: Efficacy of mosquito traps for collecting potential West Nile mosquito vectors in a natural Mediterranean wetland. Am J Trop Med Hyg. 2012, 86: 642-648. 10.4269/ajtmh.2012.11-0326.PubMed CentralView ArticlePubMedGoogle Scholar
- Horsfall WR: Mosquitoes: Their Bionomics and Relation to Disease. 1972, New York: Haftner Pub. Co, 723-Google Scholar
- Henn JB, Metzger ME, Kwan JA, Harbison JE, Fritz CL, Riggs-Nagy J, Shindelbower M, Kramer VL: Development time of Culex mosquitoes in stormwater management structures in California. J Am Mosq Control Assoc. 2008, 24: 90-97. 10.2987/5644.1.View ArticlePubMedGoogle Scholar
- Benincasa F, Maracchi G, Rossi P: Agrometeorologia. 1991, Bologna: Patron EditoreGoogle Scholar
- Baker D, Sharratt B, Chiang H, Zandlo J, Ruschy D: Base temperature selection for the prediction of European corn borer instars by the growing degree-day method. Agric For Meteorol. 1984, 32: 55-60. 10.1016/0168-1923(84)90028-5.View ArticleGoogle Scholar
- Wang J, Ogden NH, Zhu H: The impact of weather conditions on Culex pipiens and Culex restuans (Diptera: Culicidae) abundance: a case study in Peel region. J Med Entomol. 2011, 48: 468-475. 10.1603/ME10117.View ArticlePubMedGoogle Scholar
- Loetti V, Schweigmann N, Burroni N: Development rates, larval survivorship and wing length of Culex pipiens (Diptera: Culicidae) at constant temperatures. J Nat Hist. 2011, 45: 2203-2213. 10.1080/00222933.2011.590946.View ArticleGoogle Scholar
- Clements AN: Biology of Mosquitoes. Volume 1: Development, Nutrition and Reproduction. 1992, London: Chapman & Hall, 509-Google Scholar
- Vinogradova EB: Diapause in aquatic insects, with emphasis on mosquitoes. Diapause Aquat Invertebr Theory Hum Use. Edited by: Alekseev VR, de Stasio BT, Gilbert JJ. 2007, Dordrecht, The Netherlands: Springer, 83-113.View ArticleGoogle Scholar
- Vinogradova EB, Shaikevich EV: Morphometric, physiological and molecular characteristics of underground populations of the urban mosquito Culex pipiens Linnaeus f. molestus Forskål (Diptera: Culicidae) from several areas of Russia. Eur Mosq Bull. 2007, 22: 17-24.Google Scholar
- Savage HM, Anderson M, Gordon E, McMillen L, Colton L, Charnetzky D, Delorey M, Aspen S, Burkhalter K, Biggerstaff BJ, Godsey M: Oviposition activity patterns and West Nile virus infection rates for members of the Culex pipiens complex at different habitat types within the hybrid zone, Shelby County, TN, 2002 (Diptera: Culicidae). J Med Entomol. 2006, 43: 1227-1238. 10.1603/0022-2585(2006)43[1227:OAPAWN]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Pinheiro JC, Bates DM: Mixed-Effects Models in S and S-PLUS. 2000, New York: Springer, 537-View ArticleGoogle Scholar
- Burnham KP, Anderson DR: Multimodel inference: understanding AIC and BIC in model selection. Sociol Methods Res. 2004, 33: 261-304. 10.1177/0049124104268644.View ArticleGoogle Scholar
- Magee L: R2 Measures based on Wald and likelihood ratio joint significance tests. Am Stat. 1990, 44: 250-253.Google Scholar
- Vonesh EF, Chinchilli VM, Pu K: Goodness-of-fit in generalized nonlinear mixed-effects models. Biometrics. 1996, 52: 572-587. 10.2307/2532896.View ArticlePubMedGoogle Scholar
- Burnham KP, Anderson DR: Model Selection and Multi-Model Inference: a Practical Information-Theoretic Approach. 2002, New York: Springer, 488-2Google Scholar
- Johnson JB, Omland KS: Model selection in ecology and evolution. Trends Ecol Evol. 2004, 19: 101-108. 10.1016/j.tree.2003.10.013.View ArticlePubMedGoogle Scholar
- Whittingham MJ, Stephens PA, Bradbury RB, Freckleton RP: Why do we still use stepwise modelling in ecology and behaviour?. J Anim Ecol. 2006, 75: 1182-1189. 10.1111/j.1365-2656.2006.01141.x.View ArticlePubMedGoogle Scholar
- Pinheiro J, Bates D, DebRoy S, Sarkar D, R Development Core Team: nlme: Linear and Nonlinear Mixed Effects Models. 2012Google Scholar
- R Development Core Team: R: A Language and Environment for Statistical Computing. 2012, Vienna: R Foundation for Statistical ComputingGoogle Scholar
- Seber GAF, Lee AJ: Linear Regression Analysis. 2003, New York: John Wiley and Sons, Ltd, 557-View ArticleGoogle Scholar
- Buffington JD: Temperature acclimation of respiration in Culex pipiens pipiens (Diptera: Culicidae) and the influence of seasonal selection. Comp Biochem Physiol. 1969, 30: 865-878. 10.1016/0010-406X(69)90041-3.View ArticleGoogle Scholar
- Kingsolver JG: Weather and the population dynamics of insects: integrating physiological and population ecology. Physiol Zool. 1989, 62: 314-334.Google Scholar
- Lebl K, Brugger K, Rubel F: Predicting Culex pipiens/restuans population dynamics by interval lagged weather data. Parasit Vectors. 2013, 6: 129-10.1186/1756-3305-6-129.PubMed CentralView ArticlePubMedGoogle Scholar
- Roiz D, Vazquez A, Rosà R, Muñoz J, Arnoldi D, Rosso F, Figuerola J, Tenorio A, Rizzoli A: Blood meal analysis, flavivirus screening, and influence of meteorological variables on the dynamics of potential mosquito vectors of West Nile virus in northern Italy. J Vector Ecol. 2012, 37: 20-28. 10.1111/j.1948-7134.2012.00196.x.View ArticlePubMedGoogle Scholar
- Agnew P, Haussy C, Michalakis Y: Effects of density and larval competition on selected life history traits of Culex pipiens quinquefasciatus (Diptera: Culicidae). J Med Entomol. 2000, 37: 732-735. 10.1603/0022-2585-37.5.732.View ArticlePubMedGoogle Scholar
- Brook BW, Bradshaw CJ: Strength of evidence for density dependence in abundance time series of 1198 species. Ecology. 2006, 87: 1445-1451. 10.1890/0012-9658(2006)87[1445:SOEFDD]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Tamba M, Bonilauri P, Bellini R, Calzolari M, Albieri A, Sambri V, Dottori M, Angelini P: Detection of Usutu virus within a West Nile virus surveillance program in Northern Italy. Vector Borne Zoonotic Dis. 2011, 11: 551-557. 10.1089/vbz.2010.0055.View ArticlePubMedGoogle Scholar
- Bonilauri P, Bellini R, Calzolari M, Angelini R, Venturi L, Fallacara F, Cordioli P, Angelini P, Venturelli C, Merialdi G, Dottori M: Chikungunya virus in Aedes albopictus, Italy. Emerg Infect Dis. 2008, 14: 852-854. 10.3201/eid1405.071144.PubMed CentralView ArticlePubMedGoogle Scholar
- Roiz D, Vázquez A, Rosso F, Arnoldi D, Girardi M, Cuevas L, Perez-Pastrana E, Sánchez-Seco MP, Tenorio A, Rizzoli A: Detection of a new insect flavivirus and isolation of Aedes flavivirus in Northern Italy. Parasit Vectors. 2012, 5: 223-10.1186/1756-3305-5-223.PubMed CentralView ArticlePubMedGoogle Scholar
- Mills JN, Gage KL, Khan AS: Potential influence of climate change on vector-borne and zoonotic diseases: a review and proposed research plan. Environ Health Perspect. 2010, 118: 1507-1514. 10.1289/ehp.0901389.PubMed CentralView ArticlePubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.