- Open Access
Does bird metabolic rate influence mosquito feeding preference?
Parasites & Vectorsvolume 11, Article number: 110 (2018)
Host selection by mosquitoes plays a central role in the transmission of vector-borne infectious diseases. Although interspecific variation in mosquito attraction has often been reported, the mechanisms underlying intraspecific differences in hosts’ attractiveness to mosquitoes are still poorly known. Metabolic rate is related to several physiological parameters used as location cues by mosquitoes, and so potentially affect host-vector contact rates. Therefore, individual hosts with higher metabolic rates should be more attractive to host-seeking mosquitoes. Here, we experimentally investigated the role of bird metabolic rate in the feeding preferences of Culex pipiens (Linnaeus), a widespread mosquito vector of many pathogens affecting human and wildlife health.
Passer domesticus (Linnaeus) pairs containing one bird treated with 2,4-dinitrophenol (DNP) and the other injected with phosphate-buffered saline solution (PBS) (i.e. control) were simultaneously exposed overnight to mosquitoes. The treatment did not affect the proportion of mosquitoes biting on each individual. However, mosquito feeding preference was negatively associated with bird resting metabolic rate but positively with bird body mass. These two variables explained up to 62.76% of the variations in mosquito feeding preference.
The relationships between mosquito feeding preferences and individual host characteristics could be explained by enhanced anti-mosquito behaviour associated with higher metabolic rates. The potential role of cues emitted by hosts is also discussed. Thus, individuals with high metabolism may actively avoid being bitten by mosquitoes, despite releasing more attractant cues. Since metabolic rates can be related to individual differences in personality and life history traits, differences in mosquitoes’ feeding preferences may be related to intraspecific differences in exposure to vector-borne pathogens.
Mosquitoes (Diptera: Culicidae) are responsible for the transmission of multiple vector-borne pathogens that cause diseases such as malaria, West Nile fever and yellow fever . Host selection by mosquitoes is recognized as a key factor affecting pathogen amplification and transmission risk since it drives host-vector contact rates [2, 3]. Differential mosquito biting preferences have been reported at host interspecific level [3,4,5,6,7,8], but also among individuals within species [9,10,11]. Recent studies demonstrated that host characteristics such as body size, age, sex, reproductive status, breath constituents, and health status are proximate factors driving variations in individual host attractiveness to mosquitoes [10,11,12,13,14,15]. These studies hypothesized that mosquito feeding preferences may be driven by host cues including CO2 plumes, odours, sweat and movements, which are often associated with host metabolism. Understanding the individual characteristics underlying these asymmetries in mosquito blood-feeding is of great importance as variability in contact rates could result in heterogeneous individual risk for pathogen transmission [16, 17]. However, the mechanisms underlying host intraspecific differences in mosquito attraction are still poorly understood .
Mosquitoes use visual, thermal and chemical cues to detect their hosts . Metabolic rates of host animals are related to the activity and physiology of an individual . This in turn is directly linked to the emission of CO2, heat and humidity , which may enhance mosquito attraction . On the other hand, metabolic rate may also be associated with defensive behaviour, with individuals of higher metabolism being more restless and hence more difficult for mosquitoes to bite. In addition, body mass (BM) is usually a positive correlate of metabolic rate in many organisms . In birds, metabolic rate is positively associated with BM as supported by a study on 231 species , but a negative association was also reported within small-sized hummingbirds . Owing to the potential link with metabolic rate, BM may also affect the emission of multiple host-seeking cues as well as the defensive behaviour, which may affect mosquito attraction. A positive link between BM and blood-sucking insect attraction has been reported by a high number of authors working with groups of arthropod vectors such as mosquitoes [25, 26], biting midges Culicoides  and blackflies [28, 29]. However, the specific role of avian BM in mosquito attraction at intraspecific level remains poorly studied. Finally, a potential bias in studies of mosquito preference is that much research to date has focused on mosquito attractants and whether or not one particular individual will be chosen as a host before another one. However, under natural conditions individuals exposed to mosquito bites are in many cases surrounded by other potential hosts. Therefore, the likelihood of being bitten may not only depend on these attractant factors but also on the host composition, that is, whether or not they are surrounded by more attractive and/or susceptible counterparts (i.e. the infection intensity by blood parasites) .
Despite the potential importance of host metabolic rate in mosquito feeding preferences, to the best of our knowledge no study has yet experimentally tested the link between host metabolic rate and mosquito host selection. To do that, we assigned house sparrows Passer domesticus (Linnaeus) to two experimental treatments: birds injected with 2,4-dinitrophenol (DNP) or treated as controls, and exposed them to the bites of the mosquito Culex pipiens (Linnaeus) in pairs. Culex pipiens is widespread and acts as the main bridge vector for a number of pathogens affecting humans and wildlife including house sparrows (e.g. Haemoproteus , Plasmodium , West Nile virus and Saint Louis encephalitis virus, reviewed in ). Our experimental approach simulates a situation of host selection where birds with different physiological conditions grouped together during the breeding season or at roosts, and mosquitoes are active in seeking hosts.
Here we tested the hypothesis that mosquitoes preferentially bite individuals with higher metabolic rates, as they may release more attractant cues for host-seeking mosquitoes. Alternatively, they will be bitten less than their counterparts with lower metabolic rates, as other factors related to bird metabolism, such as anti-mosquito behaviour, may overrule the effect of attractant cues on determining mosquito bites.
Mosquito larvae were collected from Cañada de los Pájaros (Seville, Spain) in summer 2014 and were reared in plastic trays containing water in climatic chambers. Larvae were supplied with shrimp food (Mikrozell 20 ml/22 g; Dohse Aquaristik GmbH & Co. KG, D-53501, Gelsdorf, Germany). Mosquitoes were kept at 27 (±1) °C and 65–70% relative humidity (RH) under a photoperiod of 12:12 h (Light:Dark). One to 5 days after emergence, adult mosquitoes were anaesthetized with diethyl ether  and then sexed and identified to species level  on chilled Petri dishes using a stereomicroscope (Nikon SMZ-645, Tokyo, Japan). Female Cx. pipiens were retained and placed in insect-rearing cages (BugDorm-43030F, 32.5 × 32.5 × 32.5 cm; MegaView Science, Taiwan) in the same chamber conditions as above with ad libitum access to 1% sugar solution until 10–19 days-old. Mosquitoes were deprived of the sugar solution 24 h prior to the experiment and maintained with water until 12 h before the experiment.
Bird sampling and maintenance
Thirty juvenile house sparrows were trapped in Huelva province (southern Spain) in July 2014 using mist nets. We chose wild house sparrows as vertebrate hosts because this species is a natural reservoir for multiple vector-borne pathogens [34,35,36] and has been reported to be one of the preferred hosts of several mosquito species including Cx. pipiens [7, 37]. Yearling birds were individually marked with metal rings and weighed with a digital scale (Pesola-MS500, Pesola©, Switzerland). A small blood sample from each bird was taken using jugular venipuncture for future molecular analyses (see below). Subsequently, birds were transported to the Animal Experimentation Unit at the Estación Biológica de Doñana (EBD-CSIC) and kept in pairs in cages (58.5 × 25 × 36 cm) in a vector-free room at 22 ± 1 °C and a 12 L:12D photocycle. Water and food (mixed grain) were provided ad libitum. Two to 5 days after finishing the experiment, birds were released at the site of capture.
Measurements of bird resting metabolic rates and body mass
The 30 house sparrows used in this study consisted of 21 males (10 DNP and 11 control) and 9 females (5 DNP and 4 control). Of the 15 pairs, 9 contained a male and a female bird and 6 pairs included 2 male birds.
The RMR of each bird was measured as the minimum oxygen consumption under post-absorptive digestive conditions during its resting cycle [38, 39]. RMR was measured during a 12 h period from 20:00 h to 08:00 h using an open-circuit respirometer (Sable Systems International, Las Vegas, NV, USA). Oxygen consumption (ml O2/min) was estimated as the lowest value of the averages of 10 min runs . Birds’ BM was recorded before RMR measurements were taken.
Bird treatments and blood-feeding assays
The night following the RMR measurements, half of the birds (n = 15) were randomly injected subcutaneously with 0.2 mg of DNP diluted in 0.04 ml of phosphate-buffered saline solution (PBS) (DNP group), while the remaining birds (n = 15) were injected with the same volume of PBS (control group). DNP is an artificial decoupler of oxidative phosphorylation  and, acting as a protonophore, facilitates the leak of the protons that build up the force to drive ATP synthesis and results in poor connection between oxidation and phosphorylation. This induces an increase in the metabolic rate (i.e. oxygen consumption) to compensate for mitochondrial inefficiency and to meet energy demands . Immediately after injection, a pair of birds consisting of a DNP and a control bird was exposed to 10–19 days-old unfed female mosquitoes in the dark for 12 h from 20:00 h to 08:00 h (activity peak of Cx. pipiens, see [43, 44]). In all, 15 trials over 3 nights were conducted. In each trial, a birdcage (38.5 × 26 × 5.5 cm) containing a pair of birds was exposed to an average of 190 (range: 181–198) unfed Cx. pipiens females in insect-rearing tents (BugDorm-3120, white, 60 × 60 × 60 cm). Mosquitoes and birds were allowed to move and come into contact without any restrictions, as mosquitoes were able to freely enter the birdcages. At the end of each trial, blood-engorged mosquitoes were aspirated from inside tents, counted and stored at -20 °C.
The RMR could not be measured immediately following DNP and PBS injection because the mosquito exposure trials were taking place. Thus, in order to assess the effect of DNP injection on bird RMR, we captured eight additional house sparrows; four of them were injected with DNP and the other four with PBS as in the previous experiment. Immediately after the injection, the RMR of these individuals was recorded during 12 h using the same approach reported above.
We isolated genomic DNA from blood samples taken from birds using the DNA Kit Maxwell® 16LEV (Promega, Madison, WI, USA) . Birds were molecularly sexed and their Plasmodium, Haemoproteus and Leucocytozoon infection status were determined . To reduce any potential effect of host infection status on mosquito host selection (see [12, 15]), only birds without detectable infection by these parasites were included in the experimental procedure.
Thirty engorged mosquitoes were randomly selected after each trial to determine the origin of their blood meals; the only exception was one trial that produced only nine engorged mosquitoes, which were all analysed. Mosquito abdomens were separated from the head-thorax using sterile pipette tips and Petri dishes on an ice surface. Genomic DNA of the blood meal was isolated using the HotSHOT procedure (see [47, 48]). DNA samples were stored at -20 °C until PCR amplification analyses.
The sex of birds bitten by each mosquito was determined from the blood meal [49, 50]. We used the primer pair P2 (5′-TCT GCA TCG CTA AAT CCT TT-3′) and P8 (5′-CTC CCA AGG ATG AGR AAY TG-3′) that targets the sex-related chromo-helicase-DNA-binding gene (CHD). PCR amplification was carried out in a total volume of 25 μl in thermal cyclers (BIO-RAD T100, Hercules, USA; and Agilent Sure Cycler 8800, Santa Clara, USA). The reaction conditions and cycle temperatures are described in . Positive amplifications were visualized in 3% agarose gels. This procedure was used to partially identify the origin of the mosquito blood meals. In particular, for the nine pairs containing a male and a female bird, the blood meals with one-band amplification were identified as male-derived blood meals. Blood meals providing two bands of amplification were identified as the blood meals taken from female birds exclusively or as the mixed blood meals taken from both male and female birds.
These two-band samples and those from six bird pairs including two males were processed using eight different primer pairs to target different microsatellite fragments of the genotyped birds (see Additional file 1: Table S1 ). Microsatellite amplifications were conducted with a total volume of 20 μl for each sample containing 2 μl of extracted DNA sample, 2 μl of PCR buffer (10×), 0.6 μl of MgCl2 (50 mM), 0.16 μl of dNTPs (25 mM), 0.1 μl of Taq, 13.54 μl of H2O and 0.8 μl of primer for two DNA strands, respectively. Positive amplifications were visualized in 3% agarose gels to identify homozygous (one band) and heterozygous (two bands) individuals for each microsatellite and compared between birds from the same trial pair.
For pairs composed of two males, we selected pairs of microsatellite primers having mutually exclusive amplification patterns for each bird of the pair. This procedure allows birds to be identified and reduce the cost of sequencing. Samples with one-band amplification for either of the pair of primers were identified as blood meals from either one of the pair of birds, while two-banded amplifications for the two microsatellites were identified as mixed-blood meals. For those cases where birds showed a similar amplification pattern, we sequenced four different microsatellites (Pdo A08, B01, D09 and F09; see also Additional file 1: Table S1) from bird blood samples and mosquito blood meals using the 3130xl ABI Genetic Analyzer (Applied Biosystems, Foster City, USA). Alleles were scored using GENEMAPPER v.3.7 (Applied Biosystems). The origin of the remaining samples was resolved by comparing the size of alleles amplified by multiple primer pairs. Consequently, for each trial we obtained the number of mosquitoes that had bitten each individual and the number of mosquitoes that had bitten both birds. To assess the reliability of the assignment of blood origin, both the sex determination and microsatellite genotyping were run in duplicate for 52 and 12 samples, respectively. No inconsistent results were found.
The RMR in the DNP and control groups prior to the experimental injection was compared with Generalized Linear Models (GLMs) with normally distributed errors. The same procedure was also used to test for differences in RMR between control and DNP birds immediately after injection. Generalized Linear Mixed Models (GLMMs) with binomial error and logit link function were used to test the potential differences in the number of fed mosquitoes in relation to total introduced mosquitoes (included as the binomial denominator) between heterogeneous (containing one male and one female bird) and homogenous (containing two males) trials. We also used GLMMs to test the relationship between bird RMR and mosquito feeding preference in the trials. In this case, the number of engorged mosquitoes on a particular bird was analysed as a binomial variable with the total number of engorged mosquitoes as the binomial denominator. This response variable ‘feeding preference’ was incorporated into models codified as the number of mosquitoes that fed on one bird with respect to the number of mosquitoes that fed on the other bird within a pair (without mixed blood meals) using the cbind function. Before fitting models, bird RMR was logistic-transformed to attain normality . In the models, BM and RMR were introduced as covariates; bird sex, treatment, the interaction between sex and treatment and between treatment and RMR were incorporated as explanatory factors; and bird identity was included as a random factor in order to cope with the overdispersion found in models with a count response . In addition, bird pair was also included as a random factor as some pairs were composed of birds of different gender and so direct comparisons between pair members could not be conducted without controlling for confounding variables such as bird sex and BM. The multi-collinearity of explanatory variables was first assessed by calculating the generalized variance inflation factors (gVIFs) and, as these gVIF values were < 4 for the two continuous variables, both were incorporated into further analyses . Model selection was based on the second order Akaike’s information criteria (AICc). Delta AICc (ΔAICc) was calculated as the difference in AICc between the model with the lowest AICc and other models.
First, we fitted a global model containing all the predictors using the lme4 package v.1.1 . We standardized input variables before model analysis using the arm package v.1.8 . We then derived a set of sub-models (including the null model, which contained only the intercept) from the global model by using the dredge function implemented in the MuMIn package v.1.15 . A ‘top model set’ was created by selecting those models with a difference of ΔAICc < 2. If more than one model was selected in the top model set, we performed a model-averaging approach to summarize the results using the MuMIn package v.1.15 . Finally, as a measure of goodness-of-fit for mixed models, we calculated the explained variance (conditional R2) for each of the selected top models . All analyses were carried out in R software v.3.2.5 .
The percentage of fed mosquitoes in relation to total introduced mosquitoes in each trial varied from 5.46 to 40.6% (mean ± SE = 27.2 ± 2.51%; see also Additional file 2: Table S2). The number of fed mosquitoes between heterogeneous (containing one male and one female bird) and homogenous (containing two males) trials did not differ significantly (estimate ± SE = 0.358 ± 0.282, z = 1.142, P = 0.25). The blood meal origin of 429 mosquitoes was identified to the individual level. The mean (± SE) number of mosquitoes with a blood meal derived from a single individual was 12.63 ± 1.74 (range: 0–30; Additional file 2: Table S2). An average of 3.33 (range: 0–7) mosquitoes contained mixed blood meals per trial. In subsequent analyses we only present the results excluding these mixed blood meals as results including this data were qualitatively the same (data not shown).
Prior to the treatment, the resting metabolic rate (RMR) did not differ significantly between birds assigned to the DNP and control groups (estimate ± SE = -0.001 ± 0.017, t(29) = -0.027, P = 0.98). Two top models (Table 1) were selected according to the AICc criterion. The explained variance (conditional R2) was 62.76% (model with BM and RMR (logistic-transformed) as explanatory variables) and 54.84% (model with only BM). Neither of these models included the experimental treatment and the treatment had no significant effect when added to the model with the lowest AICc (estimates ± SE = -0.086 ± 0.766, z = -0.112, P = 0.91). The averaged estimates indicated that feeding preference was positively associated with BM but negatively correlated to RMR (Table 2). The relative importance of BM and RMR was 1.00 and 0.63, respectively. The 95% confidence intervals for the parameter estimates did not include zero, indicating that these two predictors significantly influenced feeding preference (Table 2, Fig. 1). These results were not an artefact caused by collinearity between both variables since the correlation coefficient between RMR and BM was very low and not significant (r = 0.223, P = 0.235), and the results did not change qualitatively when using the residuals of RMR against BM as a predictor instead of RMR (results not shown).
In order to identify causes explaining the non-significant effect of the experimental treatment, we measured the RMR of eight birds immediately after the injection with DNP (n = 4) or treated as controls (n = 4). The RMR of birds did not statistically differ between the two experimental groups (estimate ± SE = -0.073 ± 0.245, t(7) = -0.297, P = 0.78).
In this study we tested the relationship between host metabolic rate and the feeding preference of Cx. pipiens mosquitoes. We found that these mosquitoes preferred to feed on birds with higher BM but lower RMR. To date this is the first evidence that host metabolic rate does affect mosquito feeding preference. As our focal species Cx. pipiens is an important vector for multiple infectious diseases, identifying factors affecting the biting preferences of this mosquito species may throw some light on the epidemiology of these pathogens.
The positive association between feeding preference and BM suggests that larger individuals may release more cues that facilitate their detection and location by mosquitoes along multiple pathways such as vision, motion and odour. On the other hand, host breath may contain allomonal properties that potentially reduce mosquito attraction, as shown in human hosts with a wind-tunnel experiment by Mukanaba et al. . In our case, however, birds and mosquitoes were kept in close proximity but under complete darkness and in a windless environment. Hence the role of visual cues is probably limited here and olfactory cues as well as heat and humidity from two birds may mix up and so may not serve for host discrimination but may rather indicate the specific body parts as biting sites for mosquitoes [21, 60]. In addition, rather than the number of attracted mosquitoes, our study measured mosquito host preference in term of the number of blood-fed mosquitoes, which may be subject to the influence of a host defensive response. In this context, motion (including anti-mosquito behaviour) rather than vision, odour or breath could be the most important factor determining mosquito feeding patterns. Nonetheless, we cannot rule out the possibility that the latter cues could also affect, at least in part, our results. Smaller individuals tend to move more frequently than larger ones , and avian defensive behaviour can greatly affect mosquito feeding success [62,63,64,65]. This could explain why birds with higher BM were bitten more frequently by mosquitoes. Another non-mutually exclusive explanation could be that larger individuals are an easier prey for mosquitoes since, when compared to smaller or more active (higher RMR) individuals, they may offer larger biting surfaces and be less proficient at avoiding bites. Our study adds to the large body of evidence showing a positive correlation between host body size/mass and the attraction of different insect including mosquitoes [25, 26], biting midges  and blackflies [28, 29] and highlights the importance of host size on mosquito blood-feeding at the intraspecific level.
Contrary to our prediction, birds with lower RMR suffered more mosquito bites than individuals with higher RMR. Higher metabolic rate is expected to be associated with the increased emission of the cues used by host-seeking mosquitoes . Mosquito blood-feeding is a complex behaviour that includes different phases, from appetitive behaviour to a consummatory reaction and the cessation of feeding . In our study, birds were placed close (within 60 cm) to mosquitoes, and initially the heat and humidity released by hosts may have been used as clues for the detection by mosquitoes [21, 67]. However, after approaching their hosts, the success of blood-feeding is largely determined by bird behaviour since mosquitoes avoid those individuals/species that are more active at the time of biting . In our study, birds were able to move freely during the exposure to mosquitoes and hence may have performed anti-mosquito behaviour to protect themselves from bites. The defensive behaviour displayed by birds against mosquitoes (i.e. foot stomping, head and wing movement, tail shaking ) may reduce the ability of mosquitoes to complete a blood meal [63,64,65] but are also energetically costly . Animals with high RMR may be more active, aggressive, explorative and bold, while their low RMR counterparts may be calmer and shyer and have the tendency to avoid novel situations . In this context, although birds with higher metabolic rates might have attracted more mosquitoes due to a greater emission of host location cues, the final feeding number of mosquitoes on this type of birds could be lower than birds with lower metabolic rates that may perform less intense anti-mosquito behaviour. Potential differences in anti-mosquito behaviours between bird classes, could explain results from  who found that adult pigeons (Columba livia) attracted more mosquitoes than juvenile individuals, but mosquito feeding success (proportion fed) was greater on juveniles than adult birds, Therefore, intensive movements powered by higher RMR could explain why mosquitoes bite preferably birds with lower RMR. Nonetheless, we did not perform direct observation of host anti-mosquito behaviour during night exposure to mosquitoes, and so the effect of host defensiveness on mosquito feeding preference calls for further research.
Despite the initial differences in RMR, we did not find any significant effect of the experimental treatment on the mosquito feeding preference. The DNP administration did not significantly affect the RMR of birds, as no differences in RMR were found between DNP-treated and control birds during the following 12 h after injection. Previous studies have recorded an increase in metabolic rates as a result of mild mitochondrial uncoupling by DNP administration in species including invertebrates , amphibians , birds  and mammals . However, the efficacy of DNP in most of the cases was very short in time . DNP can be quickly eliminated from the organism, and, for example, within 24 h up to 98% of DNP have been eliminated in ducks and rabbits , and the metabolic rate returned to normal values a day after injection . If the efficacy of treatment did not significantly affect the RMR of birds or/and the efficacy lasted only during the period before mosquitoes were able to bite birds, this could explain, at least in part, why no significant effects were found in this study. In addition, the efficacy of DNP on RMR of animals may depend on the dose , route of administration and experimental conditions, such as temperature . The subcutaneous administration of DNP to birds in our study may have resulted in a slower release of drug into the blood than oral DNP administration, delaying the impact on bird metabolic rate. Even so, it is possible that the dose injected to birds was not enough to modify the bird RMR, as an oral dose of 5 mg/l of DNP was reportedly insufficient to noticeably affect the metabolic rate of zebra finches . Further studies are necessary in order to identify the effective dose of DNP to modify the RMR of wild house sparrows without increasing mortality or producing long-term damages in bird health.
Our study represents mosquitoes’ choice between different hosts emitting different physiological and behavioural cues. Although a bird with higher metabolism may release more cues that attract mosquitoes, the final outcome of mosquito feeding patterns may also depend on the surrounding hosts, which could reduce the individual risk of being attacked [79, 80], i.e. the per capita bird exposure to infected mosquitoes may be less given the encounter-dilution effect .
Hosts’ metabolic rates and body mass may influence mosquito feeding preference at intraspecific host level. As metabolism is closely related to individual differences in personality , behaviour  and life history traits , these findings may have important implications for individual exposure to mosquito bites and consequently for the amount of exposure to vector-borne diseases.
Akaike information criterion with a correction for finite sample sizes
Analysis of variance
Chromo helicase DNA binding gene
- CO2 :
Generalized linear mixed models
Generalized variance-inflation factors
Phosphate buffered saline solution
Polymerase chain reaction
Resting metabolic rate
Becker N, Petrić D, Boase C, Lane J, Zgomba M, Dahl C, et al. Mosquitoes and their control. 2nd ed. Heidelberg, Dordrecht, London, New York: Springer; 2010.
Kilpatrick AM, Daszak P, Jones MJ, Marra PP, Kramer LD. Host heterogeneity dominates West Nile virus transmission. Proc R Soc B. 2006;273(1599):2327–33.
Muñoz J, Ruiz S, Soriguer R, Alcaide M, Viana DS, Roiz D, et al. Feeding patterns of potential West Nile virus vectors in south-west Spain. PLoS One. 2012;7(6):1–9.
Gingrich JB, Williams GM. Host-feeding patterns of suspected West Nile virus mosquito vectors in Delaware, 2001–2002. J Am Mosq Control Assoc. 2005;21(2):194–200.
Molaei G, Andreadis TG, Armstrong PM, Bueno R, Dennett JA, Real SV, et al. Host feeding pattern of Culex quinquefasciatus (Diptera: Culicidae) and its role in transmission of West Nile virus in Harris County, Texas. Am J Trop Med Hyg. 2007;77(1):73–81.
Farajollahi A, Fonseca DM, Kramer LD, Kilpatrick AM. “Bird biting” mosquitoes and human disease: a review of the role of Culex pipiens complex mosquitoes in epidemiology. Infect Genet Evol. 2011;11(7):1577–85.
Martínez-de la Puente J, Muñoz J, Capelli G, Montarsi F, Soriguer R, Arnoldi D, et al. Avian malaria parasites in the last supper: identifying encounters between parasites and the invasive Asian mosquito tiger and native mosquito species in Italy. Malar J. 2015;14(1):32.
Rizzoli A, Bolzoni L, Chadwick E, Capelli G, Montarsi F, Grisenti M, et al. Understanding West Nile virus ecology in Europe: Culex pipiens host feeding preference in a hotspot of virus emergence. Parasit Vectors. 2015;8(1):213.
Kelly DW. Why are some people bitten more than others? Trends Parasitol. 2001;17(12):578–81.
Gervasi SS, Burkett-Cadena N, Burgan SC, Schrey AW, Hassan HK, Unnasch TR, et al. Host stress hormones alter vector feeding preferences, success, and productivity. Proc R Soc B. 2016;283(1836):20161278.
Mukabana WR, Takken W, Killeen GF, Knols BG. Allomonal effect of breath contributes to differential attractiveness of humans to the African malaria vector Anopheles gambiae. Malar J. 2004;3(1):1.
Cornet S, Nicot A, Rivero A, Gandon S. Malaria infection increases bird attractiveness to uninfected mosquitoes. Ecol Lett. 2013;16(3):323–9.
Fernandez-Grandon GM, Gezan SA, Armour JA, Pickett JA, Logan JG. Heritability of attractiveness to mosquitoes. PLoS One. 2015;0(4):e0122716.
Busula A, Takken W, Boer J, Mukabana W, Verhulst N. Variation in host preferences of malaria mosquitoes is mediated by skin bacterial volatiles. Med Vet Entomol. 2017;31(3):320–6.
Yan J, Martínez-de la Puente J, Gangoso L, Gutiérrez-López R, Soriguer R, Figuerola J. Avian malaria infection intensity influences mosquito feeding patterns. Int J Parasitol. 2017; https://doi.org/10.1016/j.ijpara.2017.09.005.
Paull SH, Song S, McClure KM, Sackett LC, Kilpatrick AM, Johnson PT. From superspreaders to disease hotspots: linking transmission across hosts and space. Front Ecol Environ. 2012;10(2):75–82.
VanderWaal KL, Ezenwa VO. Heterogeneity in pathogen transmission: mechanisms and methodology. Funct Ecol. 2016;30(10):1606–22.
Takken W, Verhulst NO. Host preferences of blood-feeding mosquitoes. Annu Rev Entomol. 2013;58:433–53.
Nilsson JÅ. Metabolic consequences of hard work. Proc R Soc Lond B. 2002;269(1501):1735–9.
Blaxter K. Energy metabolism in animals and man. Cambridge: CUP Archive; 1989.
van Breugel F, Riffell J, Fairhall A, Dickinson MH. Mosquitoes use vision to associate odor plumes with thermal targets. Curr Biol. 2015;25(16):2123–9.
Gillooly JF, Brown JH, West GB, Savage VM, Charnov EL. Effects of size and temperature on metabolic rate. Science. 2001;293(5538):2248–51.
McKechnie AE, Freckleton RP, Jetz W. Phenotypic plasticity in the scaling of avian basal metabolic rate. Proc R Soc B. 2006;73(1589):931–7.
Opazo JC, Soto-Gamboa M, Fernández MJ. Cell size and basal metabolic rate in hummingbirds. Rev Chil Hist Nat. 2005;78(2):261–5.
Port G, Boreham P, Bryan JH. The relationship of host size to feeding by mosquitoes of the Anopheles gambiae Giles complex (Diptera: Culicidae). Bull Entomol Res. 1980;70(1):133–44.
Estep LK, McClure CJ, Burkett-Cadena ND, Hassan HK, Unnasch TR, Hill GE. Developing models for the forage ratios of Culiseta melanura and Culex erraticus using species characteristics for avian hosts. J Med Entomol. 2012;49(2):378–87.
Martínez-de la Puente J, Merino S, Tomás G, Moreno J, Morales J, Lobato E, et al. Factors affecting Culicoides species composition and abundance in avian nests. Parasitology. 2009;136(9):1033–41.
Malmqvist B, Strasevicius D, Hellgren O, Adler PH, Bensch S. Vertebrate host specificity of wild-caught blackflies revealed by mitochondrial DNA in blood. Proc R Soc Lond B. 2004;271(Suppl. 4):S152–5.
Martínez-de la Puente J, Merino S, Lobato E, Rivero-de Aguilar J, del Cerro S, Ruiz-de-Castañeda R, et al. Nest-climatic factors affect the abundance of biting flies and their effects on nestling condition. Acta Oecol. 2010;36(6):543–7.
Gutiérrez-López R, Martínez-de la Puente J, Gangoso L, Yan J, Soriguer RC, Figuerola J. Do mosquitoes transmit the avian malaria-like parasite Haemoproteus? An experimental test of vector competence using mosquito saliva. Parasit Vectors. 2016;9:609.
Ferraguti M, Martínez-de la Puente J, Muñoz J, Roiz D, Ruiz S, Soriguer R, et al. Avian Plasmodium in Culex and Ochlerotatus mosquitoes from southern Spain: effects of season and host-feeding source on parasite dynamics. PLoS One. 2013;8(6):e66237.
Lipnick RL. Narcosis induced by ether and chloroform. In: Lipnick RL, editor. Studies of narcosis. Dordrecht: Springer; 1991. p. 93–107.
Schaffner F, Angel G, Geoffroy B, Hervy J, Rhaiem A, Brunhes J. The mosquitoes of Europe, an identification and training programme. Montpellier: IRD; 2001.
Komar N, Panella NA, Burns JE, Dusza SW, Mascarenhas TM, Talbot TO. Serologic evidence for West Nile virus infection in birds in the New York City vicinity during an outbreak in 1999. Emerg Infect Dis. 2001;7(4):621.
Arrigo NC, Adams AP, Watts DM, Newman PC, Weaver SC. Cotton rats and house sparrows as hosts for North and South American strains of eastern equine encephalitis virus. Emerg Infect Dis. 2010;16(9):1373–80.
LaPointe DA, Atkinson CT, Samuel MD. Ecology and conservation biology of avian malaria. Ann N Y Acad Sci. 2012;1249(1):211–26.
Martínez-de la Puente J, Ferraguti M, Ruiz S, Roiz D, Soriguer RC, Figuerola J. Culex pipiens forms and urbanization: effects on blood feeding sources and transmission of avian Plasmodium. Malar J. 2016;15(1):589.
McNab BK. On the utility of uniformity in the definition of basal rate of metabolism. Physiol Zool. 1997;70(6):718–20.
Rodríguez A, Broggi J, Alcaide M, Negro JJ, Figuerola J. Determinants and short-term physiological consequences of PHA immune response in lesser kestrel nestlings. J Exp Zool Part A. 2014;321(7):376–86.
Hill RW. Determination of oxygen consumption by use of the paramagnetic oxygen analyzer. J Appl Physiol. 1972;33(2):261–3.
Williams GC. Natural selection, the costs of reproduction, and a refinement of Lack's principle. Am Nat. 1966;100(916):687–90.
Nicholls DG, Ferguson S. Bioenergetics. 4th ed. Cambridge: Academic Press; 2013.
Chiba Y, Kubota M, Nakamura Y. Differential effects of temperature upon evening and morning peaks in the circadian activity of mosquitoes, Culex pipiens pallens and C. pipiens molestus. Biol Rhythm Res. 1982;13(1):55–60.
Anderson JF, Main AJ, Ferrandino FJ, Andreadis TG. Nocturnal activity of mosquitoes (Diptera: Culicidae) in a West Nile virus focus in Connecticut. J Med Entomol. 2007;44(6):1102–8.
Gutiérrez-López R, Martínez-de la Puente J, Gangoso L, Soriguer RC, Figuerola J. Comparison of manual and semi-automatic DNA extraction protocols for the barcoding characterization of hematophagous louse flies (Diptera: Hippoboscidae). J Vector Ecol. 2015;40(1):11–5.
Hellgren O, Waldenström J, Bensch S. A new PCR assay for simultaneous studies of Leucocytozoon, Plasmodium, and Haemoproteus from avian blood. J Parasitol. 2004;90(4):797–802.
Alcaide M, Rico C, Ruiz S, Soriguer R, Muñoz J, Figuerola J. Disentangling vector-borne transmission networks: a universal DNA barcoding method to identify vertebrate hosts from arthropod bloodmeals. PLoS One. 2009;4(9):e7092.
Martínez-de la Puente J, Ruiz S, Soriguer R, Figuerola J. Effect of blood meal digestion and DNA extraction protocol on the success of blood meal source determination in the malaria vector Anopheles atroparvus. Malar J. 2013;12(1):109–10.
Ellegren H. First gene on the avian W chromosome (CHD) provides a tag for universal sexing of non-ratite birds. Proc Biol Sci. 1996;263(1377):1635–41.
Griffiths R, Double MC, Orr K, Dawson RJG. A DNA test to sex most birds. Mol Ecol. 1998;7(8):1071–5.
Garnier S, Durand P, Arnathau C, Risterucci AM, Esparza-Salas R, Cellier-Holzem E, et al. New polymorphic microsatellite loci in the house sparrow, Passer domesticus. Mol Ecol Resour. 2009;9(3):1063–5.
Maxwell A. The logistic transformation in the analysis of paired-comparison data. Br J Math Stat Psychol. 1974;27(1):62–71.
Harrison XA. Using observation-level random effects to model overdispersion in count data in ecology and evolution. PeerJ. 2014;2:e616.
O’brien RM. A caution regarding rules of thumb for variance inflation factors. Qual Quan. 2007;41(5):673–90.
Bates D, Maechler M, Bolker B, Walker S. lme4: Linear mixed-effects models using Eigen and S4. R package version. 2014;1(7):1–23.
Gelman A, Su Y-S, Yajima M, Hill J, Pittau MG, Kerman J, et al. arm: Data analysis using regression and multilevel/hierarchical models (R package, version 9.01). 2009.
Bartoń K. MuMIn: Multi-model inference. R package version 1.10.0. Retrieved May 14, 2014 from http://cran.r-project.org/package=MuMIn.
Nakagawa S, Schielzeth H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol Evol. 2013;4(2):133–42.
R Core Development Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2016.
Cardé RT. Multi-cue integration: how female mosquitoes locate a human host. Curr Biol. 2015;25(18):R793–5.
Mooring MS, Benjamin JE, Harte CR, Herzog NB. Testing the interspecific body size principle in ungulates: the smaller they come, the harder they groom. Anim Behav. 2000;60(1):35–45.
Day JF, Edman JD. Mosquito engorgement on normally defensive hosts depends on host activity patterns. J Med Entomol. 1984;21(6):732–40.
Darbro JM, Harrington LC. Avian defensive behavior and blood-feeding success of the West Nile vector mosquito, Culex pipiens. Behav Ecol. 2007;18(4):750–7.
Klowden MJ, Lea AO. Effect of defensive host behavior on the blood meal size and feeding success of natural populations of mosquitoes (Diptera: Culicidae). J Med Entomol. 1979;15(5–6):514–7.
Edman JD, Scott TW. Host defensive behaviour and the feeding success of mosquitoes. Int J Trop Insect Sci. 1987;8(5–6):617–22.
Browne SM, Bennett GF. Response of mosquitoes (Diptera: Culicidae) to visual stimuli. J Med Entomol. 1981;18(6):505–21.
Biro PA, Stamps JA. Do consistent individual differences in metabolic rate promote consistent individual differences in behavior? Trends Ecol Evol. 2010;25(11):653–9.
Careau V, Thomas D, Humphries M, Réale D. Energy metabolism and animal personality. Oikos. 2008;117(5):641–53.
Blackmore JS, Dow RP. Differential feeding of Culex tarsalis on nestling and adult birds. Mosq News. 1958;18(1):15–7.
Padalko V. Uncoupler of oxidative phosphorylation prolongs the lifespan of Drosophila. Biochem Mosc. 2005;70(9):986–9.
Salin K, Luquet E, Rey B, Roussel D, Voituron Y. Alteration of mitochondrial efficiency affects oxidative balance, development and growth in frog (Rana temporaria) tadpoles. J Exp Biol. 2012;215(5):863–9.
Stier A, Bize P, Roussel D, Schull Q, Massemin S, Criscuolo F. Mitochondrial uncoupling as a regulator of life-history trajectories in birds: an experimental study in the zebra finch. J Exp Biol. 2014;217(19):3579–89.
Caldeira da Silva CC, Cerqueira FM, Barbosa LF, Medeiros MH, Kowaltowski AJ. Mild mitochondrial uncoupling in mice affects energy metabolism, redox balance and longevity. Aging Cell. 2008;7(4):552–60.
Harper J, Dickinson K, Brand M. Mitochondrial uncoupling as a target for drug development for the treatment of obesity. Obes Rev. 2001;2(4):255–65.
Gehring P, Buerge J. The cataractogenic activity of 2, 4-dinitrophenol in ducks and rabbits. Toxicol Appl Pharmacol. 1969;14(3):475–86.
Dominguez S, Menkel J, Fairbrother A, Williams B, Tanner R. The effect of 2, 4-dinitrophenol on the metabolic rate of bobwhite quail. Toxicol Appl Pharmacol. 1993;123(2):226–33.
Cassuto Y. Metabolic adaptations to chronic heat exposure in the golden hamster. Am J Phys. 1968;214(5):1147–51.
Calder WA. Gaseous metabolism and water relations of the zebra finch, Taeniopygia castanotis. Physiol Zool. 1964;37(4):400–13.
Cresswell W. Flocking is an effective anti-predation strategy in redshanks, Tringa totanus. Anim Behav. 1994;47(2):433–42.
Janousek WM, Marra PP, Kilpatrick AM. Avian roosting behavior influences vector-host interactions for West Nile virus hosts. Parasit Vectors. 2014;7:399.
Krebs BL, Anderson TK, Goldberg TL, Hamer GL, Kitron UD, Newman CM, et al. Host group formation decreases exposure to vector-borne disease: a field experiment in a ‘hotspot’of West Nile virus transmission. Proc R Soc B. 2014;281(1796):20141586.
Martina Ferraguti, Alberto Pastoriza, Esmeralda Pérez and Isabel Martín helped during the field and/or laboratory work. We are also very grateful to Plácido and Maribel for access to collect mosquito larvae in Cañada de los Pájaros. Three reviewers provided valuable comments in a previous version of the manuscript.
This work was supported by the Spanish Ministry of Science and Innovation (CGL2012-30759) and European Regional Development Fund (CGL2015-65055-P). JY was supported by the State Scholarship Fund from China Scholarship Council, JB and JMP by Juan de la Cierva contracts, LG by a contract from the Excellence Projects of the Junta de Andalucía (RNM-6400) and RGL by an FPI grant.
Availability of data and materials
All data supporting the conclusions of this article are included within the article and its additional files.
All experimental procedures were approved by the CSIC Ethics Committee and Animal Health authorities and complied with Spanish legislation (CEBA-EBD-12-40).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Primers used in this study for genotyping of house sparrows. *Primers used in DNA sequencing; +Primer labelled with FAM, VIC, NED or PET; adapted from Garnier et al. . (DOC 35 kb)
Table S2. Bird characteristics, experimental settings and mosquito blood-feeding data. Abbreviations: ID, bird identity; M/F, male/female; C/DNP, control/2,4-dinitrophenol; BM, body mass (g); RMR, resting metabolic rate (ml O2 min-1); TIM, total introduced mosquitoes (count); TFM, total fed mosquitoes (count); PTFM, percentage of total fed mosquitoes (%); TAM, total analysed mosquitoes (exclusively non-mixed blood meals, count); FM, fed mosquitoes (exclusively non-mixed blood meals, count); FP, feeding preference (ratio). (DOC 70 kb)
About this article
- Mosquito blood-feeding
- Resting metabolic rate
- Host attractiveness
- Host body mass
- Host defensive behaviour
- Intraspecific difference