Figures
Abstract
Background
Mosquito feeding behaviour determines the degree of vector–host contact and may have a serious impact on the risk of West Nile virus (WNV) epidemics. Feeding behaviour also interacts with other biotic and abiotic factors that affect virus amplification and transmission.
Methodology/Principal Findings
We identified the origin of blood meals in five mosquito species from three different wetlands in SW Spain. All mosquito species analysed fed with different frequencies on birds, mammals and reptiles. Both ‘mosquito species’ and ‘locality’ explained a similar amount of variance in the occurrence of avian blood meals. However, ‘season of year’ was the main factor explaining the presence of human blood meals. The differences in diet resulted in a marked spatial heterogeneity in the estimated WNV transmission risk. Culex perexiguus, Cx. modestus and Cx. pipiens were the main mosquito species involved in WNV enzootic circulation since they feed mainly on birds, were abundant in a number of localities and had high vector competence. Cx. perexiguus may also be important for WNV transmission to horses, as are Cx. pipiens and Cx. theileri in transmission to humans. Estimates of the WNV transmission risk based on mosquito diet, abundance and vector competence matched the results of previous WNV monitoring programs in the area. Our sensitivity analyses suggested that mosquito diet, followed by mosquito abundance and vector competence, are all relevant factors in understanding virus amplification and transmission risk in the studied wild ecosystems. At some of the studied localities, the risk of enzootic circulation of WNV was relatively high, even if the risk of transmission to humans and horses was less.
Citation: Muñoz J, Ruiz S, Soriguer R, Alcaide M, Viana DS, Roiz D, et al. (2012) Feeding Patterns of Potential West Nile Virus Vectors in South-West Spain. PLoS ONE 7(6): e39549. https://doi.org/10.1371/journal.pone.0039549
Editor: Bradley S. Schneider, Global Viral Forecasting Initiative, United States of America
Received: January 30, 2012; Accepted: May 21, 2012; Published: June 22, 2012
Copyright: © 2012 Muñoz et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was partially supported by projects RNM118, RNM157 and P07-RNM-02511 of the Junta de Andalucía, the European Union FP7 grants HEALTH.2010.2.3.3-3 Project 261391 EuroWestNile, and HEALTH.2010.Single-stage Project 261504 EDENEXT, and the Spanish Ministry of Science project (CGL2009-11445). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. No additional external funding received for this study.
Competing interests: The authors have declared that no competing interests exist.
Introduction
West Nile virus (WNV) is a zoonotic mosquito-transmitted arbovirus whose enzootic cycle is maintained by birds and mosquitoes [1]. Virus amplification and maintained circulation require successive feeding events in viraemic birds but transmission to humans and/or horses (i.e. dead-end hosts) can cause disease outbreaks [2]. It is well known that virus amplification is affected by the characteristics of the mosquito species (i.e. its vector competence [3]), vertebrate host characteristics (i.e. host competence [4]) and abiotic factors (i.e. temperature and rainfall [5]). One potentially important variable in the understanding of virus amplification and transmission is the degree of vertebrate-mosquito contact or, in other words, the degree of interaction between mosquitoes and their hosts, since not all vertebrates are competent hosts for this virus or are susceptible to this disease [4].
The use of serological and molecular methods has greatly improved our capacity for describing feeding patterns of haematophagous arthropods in the wild [6], [7]. These approaches have allowed us to detect blood from a large variety of vertebrates in mosquitoes and to identify important differences in diet between seasons and/or localities [8], [9]. Previous studies have focused on mosquito feeding behaviour in rural, urban and peri-urban zones (see [10]-[13]). Nevertheless, natural ecosystems that contain diverse and abundant mosquito communities, as well as a wide range of vertebrate hosts, have to date been poorly studied (e.g. [14]). Aside from mosquito feeding preferences, vector competence and mosquito abundance are also significant parameters in the amplification of WNV [15] and could be related to virus epidemics.
Currently, most of our knowledge of WNV ecology and dynamics comes from North America as a result of the introduction of this arbovirus in 1999. Since its rapid expansion throughout this continent, WNV has caused important human and wildlife health problems [16]. Nevertheless, in Europe the patterns of WNV circulation are still unclear. WNV outbreaks occur nearly every year but in different and often widely separated regions. A number of studies have proposed that this pattern is the result of infected migratory birds arriving from Africa and seeding the virus in different areas [17]. However, other results suggest that the virus remains endemic all year in some areas of Europe, but only becomes epidemic in horses or humans under certain conditions [18]. Seroconversions in wild animals and virus isolation in mosquitoes suggest that WNV has been circulating in Spain regularly at least since 2003 [19]-[21], during which time cases of illness have been reported sporadically in birds, horses and humans [22] (ProMED-mail, Archive Numbers: 20101119.4203 and 20100925.3478). Despite the potential importance for understanding patterns of WNV (and flavivirus in general) amplification and transmission, very little is known about mosquito feeding patterns in Spain and in Europe in general [23]-[27].
In the present study, first of all we use molecular methods to describe feeding patterns in five mosquito species captured in three different wetlands in SW Spain (Doñana, Odiel and La Janda) from 2007 to 2009. We then quantify the relative influence of inter-specific, spatial and temporal variation on mosquito feeding patterns. Finally, we estimate the transmission risk of WNV in birds, horses and humans for different mosquito species and localities, and evaluate the sensitivity of these estimates to the natural range of variation in three important parameters for the transmission of WNV.
Methods
Study Area and Mosquito Samples
Our study area consists of a ‘multi-ecosystem’ landscape (e.g. forests, dunes, beaches, scrubland, marshland and rural settlements). The Doñana Natural Space occupies over 1,060 km2 [28] and receives around 350,000 visitors a year (data from last three years), with peaks during spring and summer (mainly in March–May and July–August). This area is considered one of the most important and outstanding protected spaces in Europe and is a major wintering and stopover area for migratory birds. Previous studies performed in this area have documented the circulation of WNV between resident birds [19], [29] and horses [30]. The presence of WNV, Usutu virus (USUV) and other flaviviruses in various mosquito species were detected during the sampling period for the present study [20], [21], although no case of WNV-related illness was reported. The estuary of the Odiel river (Huelva province) covers 10,000 ha of tidal marshland in close proximity to a coastal area frequented by large numbers of tourists in summer. La Janda (Cádiz province) consists of 3,000 ha of rice paddies on the site of a drained lake.
CDC light-traps baited with CO2 were placed at 12 sampling sites (Fig. 1) and operated for 24 hours every 1–2 weeks between March and November for three years (2007–2009). Trapping of mosquitoes was done with all the necessary permits from the regional Department of the Environment (Consejería de Medio Ambiente, Junta de Andalucía). Samples were preserved in dry ice and then transported to the laboratory for identification of species, gender and feeding status. Specimens were separated using Petri plates with white filter paper on a chill table and identified following the appropriate identification keys [31], [32] using a stereomicroscope. Culex mosquitoes belonging to the Univittatus complex were named as Culex perexiguus based on male genitalia following Harbach [33]. Individual engorged females, identified visually by abdominal blood contents, were stored at −20°C for subsequent blood meal analyses. The number of mosquito females captured per sampling site ranged from 2,490 to 61,984, while the number of females captured per species ranged from 12,846 (Cx. pipiens) to 99,191 (Cx. theileri), see Table 1.
Molecular Identification of Blood Meals
Individual mosquito abdomens were separated from the rest of the insects’ bodies and put into individual PCR tubes using two disposable pipette tips per specimen. DNA was isolated from the abdominal contents using the HotSHOT protocol as described in [6]. DNA extracts from the mosquito blood meals were used for a first standard Polymerase Chain Reaction (PCR) assay, whose PCR products were subsequently used for a nested PCR to amplify a fragment of vertebrate Cytochrome c Oxidase Subunit I (COI) mitochondrial gene. For both PCR assays we used previously described primers and thermal-cycling conditions [6] with a PTC-100 Programmable Thermal Controller (MJ Research). Amplicons of the appropriate size, as determined by agarose gel electrophoresis were purified using the commercial ExoSAP-IT reagent (GE Healthcare Life Sciences). Sequencing reactions were performed according to BigDye 1.1 technology (Applied Biosystems) using BCV-RV2 primer. Labelled DNA fragments were resolved through an ABI 3130xl automated sequencer (Applied Biosystems). Sequences were analysed and annotated using the software Sequencher™ v.4.5 (Gene Codes Corp., © 1991–2005) and identified by comparison with the GenBank DNA sequence database (National Center for Biotechnology Information, 2008) and the BOLD Systems platform (http://www.boldsystems.org/views/login.php) to assign unknown COI sequences to particular species. Positive identification and host species assignment were made when exact or nearly exact matches (>98%) were obtained and differentiated from the nearest most similar species by around 2%.
Statistical Analyses
To explain the presence of bird, human, horse and mammal blood in mosquitoes, we fitted Generalized Linear Mixed Models (GLMM) with binomial error distribution using ‘mosquito species’, ‘wetland’, ‘locality’ nested on ‘wetland’, ‘season’ and ‘year’ as random factors. Variance components were used to estimate the relative (%) contribution of each of these random factors to diet variation and the statistical significance determined by Likelihood Ratio tests. Four different models were built to analyse the factors related to the presence of bird, human, horse and mammal blood (including horse and human blood meals). For example, for the analyses of avian blood meals the dependent variable was scored as 1 when avian blood was detected in a female mosquito and 0 if only blood from other vertebrate groups was detected. Models were fitted with the Procedure GLIMMIX implemented in SAS 9.2. (SAS-Institute, Cary, NC, USA) using Laplace approximation, since it is considered to be more adequate than Pseudo-likelihood when fitting mixed models with binary response variables [34].
Rarefaction-based estimates on host richness were computed using the software Analytic Rarefaction v.2.0 (http://www.huntmountainsoftware.com/html/rarefaction.html) for each mosquito species.
Estimation of WNV Transmission Risk
The relative risk that a species of mosquito might infect humans, horses or birds with WNV was estimated using Kilpatrick’s formula [15], which we modified by not including the Minimum Infection Rate (MIR) values for each mosquito due to the lack of information on WNV infection rates for most of the species and localities analysed in our study. In our formula (see below) we used the proportion of avian blood meals as an indicator of the exposure of each mosquito species to WNV, since only in some bird species viraemia is high enough to infect the mosquitoes that feed on them. Unfortunately, host competence values for avian species in Europe do not exist for as many species as in North America and so for the moment it is not possible to conduct more detailed analyses that divide avian species into competent and non-competent hosts. Although high viraemias have been reported in several mammal species [35], none of these species occur in our study areas. Our risk transmission estimation is a relative measure of the number of WNV-infectious bites that a host species or group (e.g. humans, horses or birds) receives from a mosquito species in relation to mosquito abundance, its vector competence and its feeding pattern. We applied the following formulas for each mosquito species and locality:where A is the abundance of mosquitoes estimated as the mean number of females captured per trap/night, Fa is the fraction of blood meals taken from avian hosts, FH the fraction of blood meals taken from humans, and Fh the fraction of blood meals taken from horses. Cv is a measurement of the vector competence estimated as the proportion of bites from infected mosquitoes that transmit the virus. Given that no data on vector competence from Spain are available, we used previously reported data from France and South Africa [3], [36]. Our results are given as an index representing the potential total risk to humans, horses and birds for each mosquito species (see [15]).
To determine the importance of natural variation in each parameter for potential risk of WNV transmission, we performed a sensitivity analysis using the R package [37], [38]. We estimated the Sobol Index values according to the method described in [39], [40] for each parameter. These values represent the amount of variation in transmission risk that can be attributed to each parameter, both due to its pure effect (i.e. first-order index) and to interactions with other parameters (i.e. total index). In order to increase the coverage of the sampled distributions (i.e. uniform) and all possible combinations of the values from all parameters, a Latin hypercube sampling (LHS) technique was applied [41]. Briefly, two 10,000 (number of samples) by six (parameters) pseudo-random matrices were created to allow parameters to vary uniformly within their real ranges at the sampling sites. The range of variation for each parameter used were: blood meal component (Fa ranging between 0 and 1; FH ranging between 0 and 0.2; and Fh ranging between 0 and 1), mosquito abundance (A, ranging between 0.21 and 326), and virus competence (Cv, ranging between 0.1 and 0.9). Each combination corresponding to the blood meal fractions was forced to add up to one, but no restrictions were placed on the values of the other two parameters.
Results
Mosquito and Host Identification
We identified ten mosquito species from a total of 221,819 samples, of which only five (Culex modestus, Cx. perexiguus, Cx. pipiens, Cx. theileri and Ochlerotatus caspius) provided enough blood meals to perform statistical analyses. All of these five species are considered potential WNV vectors [3], [31]. A total of 911 blood meal DNA extracts were amplified via PCR, and then sequenced and identified to species level (see Table 1).
We identified 61 vertebrate host species, which included 42 birds, 17 mammals and two reptiles (see Table 1 and Fig. 2); only Culex modestus and Cx. pipiens were found to feed on reptiles. Culex perexiguus fed on 29 host species (22 birds and seven mammals). The number of host species identified for the other four mosquito species ranged between 20 and 22. Rarefaction curves (Fig. 3) grew rapidly at first, although in three Culex species (Cx. modestus, Cx. perexiguus and Cx. pipiens) the number of host species detected did not reach a plateau.
Only the main bird and mammal species are shown. Numbers above bars indicate the number of bloodmeals per mosquito species obtained for each sampled wetland. Percentage of blood meals for each vertebrate species/group are give as numbers inside the bars when higher than 5%. Only species/wetlands with more than 10 bloodmeals were included in the figure.
Factors Driving Mosquito Feeding Patterns
Patterns of variation in the prevalence of avian, mammal and horse blood meals were similar (see Table 2). Only ‘mosquito species’ (around 49–65% of explained variance) and ‘locality’ (around 30–50% of explained variance) were significantly associated to blood meal origin. ‘Season’ and ‘study year’ also explained a significant amount of variance in the presence of horse blood meals, but none in the case of avian and mammal blood meals. However, a very different pattern was found in the case of human blood meals, which were concentrated in summer, the only season that was significantly related to human blood presence (Table 2). ‘Mosquito species’, ‘locality’ and ‘year’ were unrelated to human blood presence. ‘Wetland’ did not explain a significant amount of variance in mosquito feeding on any of the vertebrate groups or species considered.
Potential Transmission Risk to Avian Species, Horses and Humans
Culex perexiguus had the highest potential for enzootic virus transmission, with a mean bird-to-bird estimated transmission risk of 31.52±62.4 (S.D.), followed by Cx. modestus (19.9±39.7) and Cx. pipiens (4.6±4.9). Culex perexiguus also had the highest transmission risk to horses (1.99±3.4), followed by Cx. theileri (0.1±0.2), although this risk was very low in comparison to birds. The potential transmission risk to humans was low for Cx. pipiens (0.025±0.05), Cx. theileri (0.02±0.05) and Oc. caspius (0.006±0.01) (see Fig. 4).
For localities with the same estimates dots have been slightly displaced horizontally to avoid complete symbol overlapping. Solid circles indicate those mosquito species and sampling sites in which WNV has been reported in mosquitoes (Vázquez et al. [20], [21]). Estimates are reported for 4 localities for Cx. modestus, 4 for Cx. perexiguus, 5 for Cx. pipiens, 10 for Cx. theileri, and 10 for Oc. caspius. One was added to all estimates to allow plotting on a log10 scale.
The sensitivity analysis suggested that ‘blood meal origin’ makes a greater contribution to transmission risk than either ‘mosquito abundance’ or ‘vector competence’ (Fig. 5A). The estimated values of enzootic transmission risk are large in some of our study areas when compared with theoretical values (Fig. 5B), but tend towards a lower range of possible values for potential transmission risk in humans and horses.
A. First order (circles) and total (triangles) sensitivity Sobol Index values for birds, humans and horses are given. Bars denote 95% confidence intervals. B. Relationship between the blood meal component and the infection risk estimated in the sensitivity analysis; i.e. Fig. A (lower panel; grey dots). The blood meal component is calculated as the product of mosquito blood meals on birds and mosquito blood meals on the group in question (birds, humans or horses). The black line is a lowess (locally weighted scatterplot smoothing) smoother that helps visualize the positive tendency between the parameter and the risk of infection; black diamonds correspond to risk values at the sampling locations.
Discussion
Our results showed that the studied mosquito species fed on a wide range of vertebrates (e.g. birds, mammals and reptiles) and thus all of them can be considered as generalists (see Fig. 2), as has been reported in other regions of the world [31], [42]. However, important inter-specfic differences in the occurrence of avian and mammal blood meals were detected: Culex modestus, Cx. perexiguus and Cx. pipiens (i.e. the main WNV vectors, see below) feed mainly on birds (78–85% of blood meals), while Cx. theileri and Oc. caspius feed mainly on mammals (81–87%) (Fig. 2). Interestingly, the rarefaction curves of the former three Culex species had the steepest gradients since they fed on a large number of hosts, which translates into a potentially more complex WNV enzootic cycle involving a larger number of avian host species. The techniques used to trap the mosquitoes may also potentially affect the characteristics of the sampled fraction of the population. CDC traps baited with CO2 capture preferentially host-seeking females or females that have not been able to successfully complete a blood meal (i.e. because they have fed on hosts with strong anti-mosquito behaviour) [43]. However, it is important to note that we found no differences in the vertebrate origin of blood meals between female mosquitoes captured with CDC light traps lured with CO2 and those captured in resting boxes (BG traps without any attractant) [44]. Additionally, this potential bias may affect the absolute values of the parameters (i.e. the proportion of mammal blood meals) but is not likely to affect the conclusions of variance decomposition analyses or the relative estimates of transmission risk.
Variance Decomposition Analyses
Our analyses indicate that the major factors involved in the feeding behaviour of these arthropods are the species and the sampling site where adult females are captured. Recently, Chaves et al. [10] concluded that mosquito feeding patterns were unrelated to vertebrate class (i.e. there were no inter-specific differences in preferences for mammals or birds). However, our results suggest that these authors’ conclusions are the product of an artefact caused by having ignored interaction strengths (i.e. the feeding frequency in each clade). Our results, nevertheless, highlight the need to include interaction strength (i.e. frequency) and the serious impact of differences in vertebrate-mosquito contact in the amplification and transmission risk potential of pathogens. Our rarefaction analyses indicate that larger sampling sizes than those currently used in most mosquito feeding/behaviour studies are also essential if accurate descriptions of the vertebrate species present in mosquito diets are to be obtained, above all if we aim to infer anything from the absence of a particular species from mosquito diet. Additionally, our variance decomposition analyses indicate that, in terms of the presence of vertebrate classes in mosquito diets (birds or mammals), large spatial variation exists, although ‘mosquito species’ does still explain a large amount of variance.
In addition to ‘mosquito species’, our results showed that ‘locality’ (i.e. spatial heterogeneity) is a factor that affects the origin of blood meals, although this effect occurred at the level of trap location rather than at the ‘wetland-complex’ level. It is important to note that Culex females do not fly long distances [31] but feed near the area of capture; consequently, vertebrate availability at such local scales will determine the origin of their blood meals. Hamer et al. [45] observed that the composition of vertebrate communities determines vector feeding patterns. Nevertheless, the nature of the relationship between vertebrate abundance and mosquito feeding is still poorly understood [10], [42], [45], [46]; for example, mosquitoes do not necessarily feed on the most abundant prey species [45], [47]. Vertebrate species may differ widely in their host competence (i.e. the sum of the probability that an infected host will transmit virus to a biting mosquito on each of the seven days following infection [4], [48]) and important differences may exist even between different species of birds [4]. Consequently, different vertebrate communities may vary dramatically in their capacity for pathogen amplification, circulation and transmission to particular host species. Unfortunately, not enough information on the host competence of European avian species exists [49] to be able to refine our analyses by dividing avian species into host-competent and non-competent groups. Thus, our results should be taken as a generalized approach in the absence of further detailed information.
In our study, the frequency of feeding on humans was only affected by the season. The low number of human blood meals reported in our study is related to the fact that we worked in natural areas with very low anthropic presence. Although the small sample size may have reduced our capacity to detect other effects (i.e. mosquito species) as statistically significant, the amount of variance explained (Table 2) does support the importance of seasonal differences in the occurrence of human blood meals. Seasonal differences in the frequency of human blood meals have been associated with the occurrence of WNV epidemics in Colorado and California, U.S.A. [50], [51], where in an urban environment the increase in feeding on humans was linked to a decrease in the abundance of the preferred avian hosts [8]. In our case, in addition to potential changes in the abundance of other preferred vertebrates, the increased presence of human blood meals during the summer season could be a reflection of the increased number of visits by humans that make them more accessible to mosquitoes. Human behaviour might have an impact on the transmission of some vector-borne virus in urban zones [52] and the same could be applicable to WNV in non-urban areas.
Mosquito Species Relevant for WNV Transmission in SW Spain
The greatest potential risk of WNV transmission in the studied wetlands in SW Spain occurs between birds (i.e. enzootic cycle). It is fifty times less between birds and horses and two hundred times less between birds and humans (Figs. 4 and 5), although risk values were found to be very heterogeneous between localities. In fact, the WNV antibody prevalences found in Doñana for birds (10.4% [29]), horses (5.3% [30]) and humans (0.5% [53]) follow the same patterns of variation as our estimated transmission risks. However, these differences are also affected by the different age structure of the populations of these organisms and by the different duration of detectable antibodies in the bloodstream after infection by WNV.
In our study area, our results for Cx. perexiguus show a high potential risk of WNV transmission from bird to bird and from bird to horse, this mosquito species being the primary candidate for enzootic (bird-to-bird) transmission and for epizootic transmission to horses. Culex pipiens is also an enzootic vector, but has only a low risk for epizootic transmission to humans. There was also a potentially high risk of transmission for Cx. modestus from bird to bird; despite being a secondary enzootic vector, this species was not implicated in epizootic cycles, unlike in the Camargue in France [24], [54]. Finally, Cx. theileri and Oc. caspius can be considered secondary epizootic vectors. Interestingly, the Cx. perexiguus and Cx. pipiens positive pools for WNV detected in Spain [20], [21] were collected in the two combinations of ‘locality-mosquitoes’ ranked among the highest transmission risk values (black dots in Fig. 4). This suggests that patterns of variation in estimated WNV transmission risk may be related to real WNV transmission patterns; nevertheless, a larger sample size (i.e. number of localities with and without WNV mosquito positive pools) is needed if we are to test this relationship statistically. The rapid accumulation of detailed studies regarding the origin of blood meals should allow us to investigate in the future how these parameters differ whenever and wherever epidemics among horses and humans occur.
Sensitivity Analyses for WNV Transmission Risk
Our sensitivity analysis ranked the blood feeding pattern (diet) as the main factor related to virus transmission, followed by mosquito abundance and vector competence. This conclusion reflects the great variability in diet composition reported in natural mosquito populations and is not an artefact caused by including two blood meal related variables as the components of diet in the formula. Note that diet is analysed in the sensitivity analyses as a simple parameter and that analyses of sensitivity conducted using the same levels of variation for the three parameters analysed (diet, mosquito abundance and vector competence) would assign the same importance to each factor. We believe that it is important to emphasize the fact that the importance of vector competence may have been underestimated in these analyses because vector competence can vary geographically [55] (but also with temperature and virus characteristics), and we have used the same values (estimated from mosquitoes Camargue in France for 3 species and from South Africa for other two species) for all the localities. However, to the best of our knowledge, no information of the potential variability in vector competence at the spatial scale analysed in our study (150 km) exists. Thus, our analyses suggest that changes in the relative frequency of blood meals in certain vertebrate groups might lead to changes in the transmission risk for potential host vertebrate species. This result reinforces the idea that broad-based surveys on host communities are needed if we are to accurately estimate transmission risks.
One lesson learnt from WNV outbreaks is that a particular outbreak is never the same as the previous one, and that the mosquito species involved in the new outbreak may differ [54]. Not surprisingly, the estimated risk of transmission for a single species varied strongly between localities (Fig. 4), thereby illustrating the potential for spatial heterogeneity in the relevance of different mosquito species for WNV circulation. Recently, Simpson et al. [51] concluded that seasonal changes in blood meal origin in Cx. pipiens was the more influential parameter explaining the intensity and timing of WNV transmission in different localties of North America. In our area, Cx. perexiguus (the European species in the Univittatus complex [33]), morphologically similar to Cx. univitattus (a major WNV vector in Africa and in the Middle East [16], [56]), is the primary enzootic vector and a potential epizootic vector for horses. These two species have been recorded as the principal vectors in several WNV outbreaks in regions such as Egypt [57], South Africa [36], [58], Kenya [59], Israel [60]; and Portugal [61]. Therefore, we should not underestimate the role of Cx. perexiguus as an enzootic and epizootic vector of WNV in SW Spain and Europe.
In conclusion, our study represents the first quantification of the relative importance of mosquito species, location and season in explaining differences in blood meal origin in mosquitoes from a natural wetland ecosystem and the consequences that this finding may have for WNV dynamics in Europe. The presence of human blood meals was only related to season, while the presence of avian or mammal blood meals depends on mosquito species and locality. Estimates of transmission risks had important spatial and inters-pecific differences and closely match what little data is available on WNV detection in the area. Diet composition had a stronger impact than either mosquito abundance or vector competence on estimated transmission risks due to high spatial and inter-specific differences in diet composition. Our findings provide invaluable information regarding vector surveillance for WNV activity and highlight the potential risk of WNV transmission to horses and between birds by a poorly known Culex species, Cx. perexiguus.
Acknowledgments
Esmeralda Perez, Oscar Gonzalez, Isabel Martín, Juani Moreno, Antonio Magallanes and Cristina Pérez participated in the collection, identification and blood meal analyses.
Author Contributions
Conceived and designed the experiments: RS SR JF. Performed the experiments: JM SR MA. Analyzed the data: DR DV JF. Contributed reagents/materials/analysis tools: AV JM SR JF DR MA DV. Wrote the paper: JM SR RS MA DV DR AV JF.
References
- 1.
Ostlund EN (2008) West Nile Fever, pp 377–385. In: OIE (Ed.), Manual of diagnostic tests and vaccines for terrestrial animals. Chapter 2.1.20. Available: http://www.oie.int/en/international-standard-setting/terrestrial-manual/access-online.
- 2. Kramer LD, Styer LM, Ebel GD (2008) A global perspective on the epidemiology of West Nile virus. Annu Rev Entom 53: 61–81.
- 3. Balenghien T, Vazeille M, Grandadam M, Schaffner F, Zeller H, et al. (2008) Vector competence of some French Culex and Aedes mosquitoes for West Nile virus. Vector Borne Zoonotic Dis 8: 589–595.
- 4. Komar N, Langevin S, Hinten S, Nemeth N, Edwards E, et al. (2003) Experimental infection of North American birds with the New York 1999 strain of West Nile virus. Emerg Infect Dis 9: 311–322.
- 5. Ruiz MO, Chaves LF, Hamer GL, Sun T, Brown WM, et al. (2010) Local impact of temperature and precipitation on West Nile virus infection in Culex species mosquitoes in northeast Illinois, USA. Parasit Vectors 3: 19.
- 6. Alcaide M, Rico C, Ruiz S, Soriguer R, Muñoz J, et al. (2009) Disentangling vector-borne transmission networks: A universal DNA Barcoding method to identify vertebrate hosts from Arthropod bloodmeals. Plos One 4: e7092.
- 7. Kent RJ (2009) Molecular methods for arthropod bloodmeal identification and applications to ecological and vector-borne disease studies. Mol Ecol Resour 9: 4–18.
- 8. Kilpatrick AM, Kramen LD, Jones MJ, Marra PP, Daszak P (2006) West Nile virus epidemics in North America are driven by shifts in mosquito feeding behavior. Plos Biology 4: e82.
- 9. Gómez-Díaz E, Figuerola J (2010) New perspectives in tracing vector-borne interaction networks. Trends Parasitol 26: 470–476.
- 10. Chaves LF, Harrington LC, Keogh CL, Nguyen AM, Kitron UD (2010) Blood feeding patterns of mosquitoes: random or structured? Front Zool 7: 1–11.
- 11. Apperson CS, Hassan HK, Harrison BA, Savage HM, Aspen SE, et al. (2004) Host feeding patterns of established and potential mosquito vectors of West Nile virus in the Eastern United States. Vector Borne Zoonotic Dis 4: 71–82.
- 12. Molaei G, Farajollahi A, Armstrong PM, Oliver J, Howard JJ, et al. (2009) Identification of bloodmeals in Anopheles quadrimaculatus and Anopheles punctipennis from eastern equine encephalitis virus foci in northeastern U.S.A. Med Vet Entomol 23: 350–356.
- 13. Murdock CC, Olival KJ, Perkins SL (2010) Molecular identification of host feeding patterns of snow-melt mosquitoes (Diptera: Culicidae): Potential implications for the transmission ecology of Jamestown Canyon virus. J Med Entomol 47: 226–229.
- 14. Cupp EW, Zhang D, Yue X, Cupp MS, Guyer C, et al. (2004) Identification of reptilian and amphibian blood meals from mosquitoes in an Eastern equine encephalomyelitis virus focus in Central Alabama. Am J Trop Med Hyg 71: 272–276.
- 15. Kilpatrick AM, Kramer LD, Campbell SR, Alleyne EO, Dobson AP, et al. (2005) West Nile virus risk assessment and the bridge vector paradigm. Emerg Infect Dis 11: 425–429.
- 16. Kilpatrick AM (2011) Globalization, land use, and the invasion of West Nile virus. Science 334: 323–327.
- 17. Jourdain E, Toussaint Y, Leblond A, Bicout DJ, Sabatier P, et al. (2007) Bird species potentially involved in introduction, amplification, and spread of West Nile virus in a Mediterranean wetland, The Camargue (Southern France). Vector Borne Zoonotic Dis 7: 15–33.
- 18. Hubálek Z (2000) European experience with the West Nile virus ecology and epidemiology: could it be relevant for the New World? Viral Immunol 13: 415–426.
- 19. Figuerola J, Soriguer R, Rojo G, Gómez-Tejedor C, Jimenez-Clavero MA (2007) Seroconversion in wild birds and local circulation of West Nile virus, Spain. Emerg Infect Dis 13: 1915–1917.
- 20. Vázquez A, Sánchez-Seco MP, Ruiz S, Molero F, Hernández L, et al. (2010) Putative new lineage of West Nile virus, Spain. Emerg Infect Dis 16: 549–552.
- 21. Vázquez A, Ruiz S, Herrero L, Moreno J, Molero F, et al. (2011) West Nile and Usutu viruses in mosquitoes in Spain, 2008–2009. Am J Trop Med Hyg 85: 178–181.
- 22. Hofle U, Blanco JM, Crespo E, Naranjo V, Jimenez-Clavero MA, et al. (2008) West Nile virus in the endangered Spanish imperial eagle. Vet Microbiol 129: 171–178.
- 23. Fyodorova MV, Savage HM, Lopatina JV, Bulgakova TA, Ivanitsky AV, et al. (2006) Evaluation of potential West Nile virus vectors in Volgograd region, Russia, 2003 (Diptera: Culicidae): species composition, bloodmeal host utilization, and virus infection rates of mosquitoes. J Med Entomol 43: 552–563.
- 24. Balenghien T, Fouque F, Sabatier P, Bicout DJ (2006) Horse-, bird-, and human-seeking behavior and seasonal abundance of mosquitoes in a West Nile virus focus of Southern France. J Med Entomol 43: 936–946.
- 25. Muñoz J, Eritja R, Alcaide M, Montalvo T, Soriguer R, et al. (2011) Host-feeding patterns of native Culex pipiens and invasive Aedes albopictus mosquitoes (Diptera: Culicidae) in urban zones from Barcelona, Spain. J Med Entomol 48: 956–960.
- 26. Roiz D, Vázquez A, Rosà R, Muñoz J, Arnoldi D, et al. (2012) 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 37: 1–9.
- 27. Valerio L, Marini F, Bongiorno G, Facchinelli L, Pombi M, et al. (2010) Host-feeding patterns of Aedes albopictus (Diptera: Culicidae) in urban and rural contexts within Rome province, Italy. Vector Borne Zoonotic Dis 10: 291–294.
- 28.
Díaz-Delgado R (2010) An integrated monitoring programme for Doñana Natural Space. In: Hurford C, Schneider M, Cowx I, editors. pp. 325–337. Dordrecht: Springer.
- 29. Figuerola J, Jimenez-Clavero MA, Lopez G, Rubio C, Soriguer R, et al. (2008) Size matters: West Nile virus neutralizing antibodies in resident and migratory birds in Spain. Vet Microbiol 132: 39–46.
- 30. Jiménez-Clavero MA, Llorente F, Sotelo E, Soriguer R, Gómez-Tejedor C, et al. (2010) West Nile virus serosurveillance in horses in Doñana, Spain, 2005 to 2008. Vet Rec 167: 379–380.
- 31.
Becker N, Petric D, Zgomba M, Boase C, Madon M, et al. (2010) Mosquitoes and their control. Berlin: Springer-Verlag. 577 p.
- 32.
Schaffner F, Angel G, Geoffrey B, Hervy J, Rhaiem A, et al. (2001) The mosquitoes of Europe, an identification and training programme. Montpellier, France: IRD.
- 33. Harbach (1999) The identity of Culex perexiguus Theobald versus ex. Univittatus Theobald in southern Europe. European Mosquito Bulletin 4: 7.
- 34. Bolker BM, Brooks ME, Clark CJ, Geange SW, Poulsen JR, et al. (2008) Generalized linear mixed models: a practical guide for ecology and evolution. Trends Ecol Evol 24: 127–135.
- 35. Platt KB, Tucker BJ, Halbur PG, Tiawsirisup S, Blitvich BJ, et al. (2007) West Nile virus viremia in Eastern chipmunks (Tamias striatus) sufficient for infecting different mosquitoes. Emerg Infect Dis 13: 831–837.
- 36. Jupp PG, McIntosh BM (1970) Quantitative experiments on the vector capability of Culex (Culex) univittatus Theobald with West Nile and Sindbis viruses. J Med Entomol 7: 371–373.
- 37.
R Development Core Team (2009) A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Available: http://www.R-project.org.
- 38.
Pujol G (2008) Sensitivity: Sensitivity analysis. R package version 1.4–0. Available: http://www.R-project.org.
- 39. Saltelli A (2002) Sensitivity analysis for importance assessment. Risk Anal 22: 579–590.
- 40.
Saltelli A, Ratto M, Andres T, Campolongo F, Cariboni J, et al. (2008) Global sensitivity analysis. West Sussex, UK: John Wiley and Sons. 292 p.
- 41.
Carnell R (2009) lhs: Latin Hypercube Samples. R package version 0.5. Available: cran.r-project.org/web/packages/lhs/index.html.
- 42.
Silver JB (2008) Mosquito Ecology: field sampling methods. Third Edition. London: Springer. 1477 p.
- 43. Kent RJ, Gonzalez-Reiche AS, Morales-Betoulle ME, Komar N (2010) Comparasion of engorged Culex quinquefasciatus collection and blood-feeding pattern among four mosquito collection methods in Puerto Barrios, Guatemala, 2007. J Am Mosq Control Assoc 26: 332–336.
- 44. Roiz D, Roussel M, Muñoz J, Ruiz S, Soriguer R, et al. (2012) Efficacy of mosquito traps for collecting potential West Nile mosquito vectors in a natural Mediterranean wetland. Am J Trop Med Hyg 86: 642–648.
- 45. Hamer GL, Kitron UD, Goldberg TL, Brawn JD, Loss SR, et al. (2009) Host selection by Culex pipiens mosquitoes and West Nile virus amplification. Am J Trop Med Hyg 80: 268–278.
- 46. Kay BH, Boyd AM, Ryan PA, Hall RA (2007) Mosquito feeding patterns and natural infection of vertebrates with Ross River and Barmah Forest viruses in Brisbane, Australia. Am J Trop Med Hyg 76: 417–423.
- 47. Tempelis CH (1975) Host-feeding patterns of mosquitoes, with a review of advances in analysis of blood meals by serology. J Med Entomol 11: 635–653.
- 48. Kilpatrick AM, Daszak P, Jones MJ, Marra PP, Kramer LD (2006a) Host heterogeneity dominates West Nile virus transmission. Proc R Soc Lond B Biol Sci 273: 2327–2333.
- 49. Wheeler SS, Barker CM, Fang Y, Armijos MV, Carroll BD, et al. (2009) Differential impact of West Nile virus on California birds. Condor 111: 1–20.
- 50. Hamer GL, Chaves LF, Anderson TK, Kitron UD, Brawn JD, et al. (2011) Fine-scale variation in vector host use and force of infection drive localized patterns of West Nile virus transmission. Plos One 6: e23767.
- 51. Simpson JE, Hurtado PJ, Medlock J, Molaei G, Andreadis TG, et al. (2012) Vector host-feeding preferences drive transmission of multi-host pathogens: West Nile virus as a model system. Proc R Soc Lond B Biol Sci 279: 925–933.
- 52. Randolph SE (2010) Human activities predominate in determining changing incidence of tick-borne encephalitis in Europe. Euro Surveill 15: 24–31.
- 53. Bernabeu-Wittel M, Ruiz-Pérez M, del Toro MD, Aznar J, Muniain A, et al. (2007) West Nile virus past infections in the general population of Southern Spain. Enferm Infecc Microbiol Clin 25: 561–565.
- 54. Hubálek Z, Halouzka J (1999) West Nile fever - a reemerging mosquito-borne viral disease in Europe. Emerg Infect Dis 5: 643–650.
- 55. Kilpatrick MA, Fonseca D, Ebel GD, Reddy M, Kramer L (2008) Spatial and temporal variation in vector competence of Culex pipiens and Cx. restuans mosquitoes for West Nile virus. Am J Trop Med Hyg 83: 607–613.
- 56. McIntosh BM, Jupp PG, Santos ID, Meenehan GM (1976) Epidemics of West Nile and Sindbis viruses in South Africa with Culex (Culex) univitattus Theobald as vector. S Afr J Sci 72: 295–300.
- 57. Taylor RM, Work TH, Hurlbut HS, Rizk F (1956) A study of the ecology of West Nile virus in Egypt. Am J Trop Med Hyg 5: 579–620.
- 58. Jupp PG (2001) The ecology of West Nile virus in South Africa and the occurrence of outbreaks in humans. Ann N Y Acad Sci 951: 143–152.
- 59. Miller BR, Nasci RS, Godsey MS, Savage HM, Lutwama JJ, et al. (2000) First field evidence for natural vertical transmission of West Nile virus in Culex univittatus complex mosquitoes from Rift Valley province, Kenya. Am J Trop Med Hyg 62: 240–246.
- 60. Orshan L, Bin H, Schnur H, Kaufman A, Valinsky A, et al. (2008) Mosquito vectors of West Nile Fever in Israel. J Med Entomol 45: 939–947.
- 61. Almeida AP, Galao RP, Sousa CA, Novo MT, Parreira R, et al. (2008) Potential mosquito vectors of arboviruses in Portugal: species, distribution, abundance and West Nile infection. Trans R Soc Trop Med Hyg 102: 823–832.