Skip to main content

Exposure of the mosquito vector Culex pipiens to the malaria parasite Plasmodium relictum: effect of infected blood intake on immune and antioxidant defences, fecundity and survival

Abstract

Background

The intake of a Plasmodium-infected blood meal may affect mosquito physiology and a series of trade-offs may occur, in particular between immune defences, reproduction and self-maintenance. We evaluated the cost of exposure to Plasmodium in the mosquito vector by investigating the effect of exposure on fecundity and survival and the implication of immune and antioxidant defences in mediating this cost.

Methods

We used the natural Culex pipiens-Plasmodium relictum association. We exposed female mosquitoes to increasing levels of parasites by allowing them to feed either on uninfected canaries, Serinus canaria, (unexposed mosquitoes) or on infected canaries with low (low exposure) or high (high exposure) parasitaemia. We recorded blood meal size, fecundity (laying probability and clutch size) and survival. We quantified the expression of genes involved in immune and antioxidant defences (nitric oxide synthase, NOS; superoxide dismutase, SOD; glucose-6-phosphate dehydrogenase, G6PDH).

Results

We found that the laying probability of exposed females decreased with increasing exposure to the parasite and with increasing SOD expression. Clutch size of exposed females was higher compared to unexposed ones for similar blood meal size and was positively correlated to the NOS expression. We found no effect of exposure on survival. After blood meal intake, SOD increased in the three groups, NOS increased in exposed females and G6PDH increased in highly exposed females only.

Conclusions

Our results illustrated a trade-off between fight against the parasite and reproduction and a cost of exposure which might be mediated by the investment in immune and/or antioxidant defences. They also showed that this trade-off could lead to opposed outcome, potentially depending on the vector physiological status. Finally, they highlighted that the ingestion of a Plasmodium-infected blood meal may affect mosquito life history traits in a complex way.

Background

Mosquito vectors are important protagonists in malaria epidemiology. They host the sexual reproduction of Plasmodium parasites responsible for malaria, leading to infectious sporozoite stages that reach the mosquito salivary glands, ready to be transmitted to a vertebrate host. Plasmodium infection is costly for mosquitoes, affecting (for instance) feeding behaviour [1, 2], resistance to nutritional stress [3], fecundity [4] and survival [5]. Yet, modifications of vector life history traits might occur due to the ingestion of a Plasmodium-infected blood meal only, even if the parasite does not develop.

Plasmodium-infected blood contains parasite stages as well as vertebrate host factors associated with infection that mosquitoes will be exposed to through their blood meal. Thus, the ingestion of a Plasmodium-infected blood meal may lead to physiological changes [612] that may be costly for the vector. First, exposure is followed by the activation of several immune effectors of the vector [13, 14] that may prevent Plasmodium development [1520], such as unspecific and highly reactive pro-oxidant compounds. Pro-oxidants are a broadly used cytotoxic defence in invertebrates [21] affecting Plasmodium development [14, 2224] but also harming the host which may suffer oxidative stress [25]. Mosquitoes may invest in antioxidant defences to protect themselves from oxidative damage [26, 27]. However parasites may also take advantage of the host antioxidant defences and a trade-off can occur in the mosquito between fight against the parasite and self-maintenance. For instance, it has been shown that mosquitoes experience local oxidative stress enhancement focused on the invasion site [26, 27] avoiding generalized oxidative damage [28]. Then, exposure is associated with vitellogenesis disruption, apoptosis and egg resorption in the ovaries [11, 2932] that can lead to delayed reproduction and/or decreased fecundity [33]. This highlights a second trade-off between fight against the parasite and reproduction [34]. This trade-off has been suggested to originate from the amino acid (arginine) requirement needed for both immune response [21] and egg production [35]. It may also originate from the antioxidant requirement needed for both self-maintenance during the fight against the parasite and reproduction. Indeed, antioxidants are important in vector fecundity as exemplified by the involvement of catalase and superoxide dismutase in egg production [36, 37]. Therefore, following the intake of a Plasmodium-infected blood meal, mosquitoes should suffer costs from the exposure even if the parasite does not develop.

These studies have enabled a better understanding of the mechanisms and pathways involved in the mosquito-Plasmodium interaction. Taken all together, they allow us to predict that the investment in immune and antioxidant defences would mediate the cost of exposure and would lead to differential resource allocation into fight against the parasite, self-maintenance or reproduction. We proposed to investigate the effect of exposure to malaria parasite Plasmodium relictum (pSGS1) on fitness components of the mosquito vector Culex pipiens, a natural host-parasite association. We experimentally exposed female mosquitoes to Plasmodium parasites by allowing them to feed on experimentally infected canaries, Serinus canaria. It has been shown that the exposure level is related to the number of parasite blood stages ingested during the blood meal [38]. We tested three exposure levels - no, low and high exposure - by using three groups of birds, respectively: uninfected birds, birds in the chronic phase of infection carrying low parasite intensity (i.e. parasitaemia) and birds in the acute phase of infection with high parasitaemia. We measured two fitness components of the female mosquitoes: the fecundity measured as the oviposition occurrence as well as the clutch size after the first gonotrophic cycle and the survival, measured as the time from blood meal to death. We also investigated if the effect of Plasmodium exposure was linked with investment in immune and antioxidant defences. We sacrificed part of the females at 5 days post-feeding (dpf, corresponding to Plasmodium oocyst development in the midgut basal side) and at 15 dpf (corresponding to Plasmodium sporozoite invasion of the salivary glands) to quantify gene expression. We measured the expression of the immune gene nitric oxide synthase (NOS) producing the pro-oxidant molecule nitric oxide (NO), a broadly used cytotoxic defence in invertebrates [21]. We measured the expression of the antioxidant defence superoxide dismutase (SOD) detoxifying the pro-oxidant molecule superoxide (O2 ), a pro-oxidant at the top of the pro-oxidative reaction chain. We also measured the expression of the glucose-6-phosphate dehydrogenase (G6PDH). It is an important enzyme implicated in the pentose phosphate pathway allowing the recycling of NADP+ to NADPH [39], which is involved in both antioxidant and immune pathways. In antioxidant pathways, NADPH can be used to produce the antioxidant glutathione and avoid oxidative damage. In immune pathways, it can be used by NOS to produce nitric oxide or by NADPH-oxidase (NOX) to produce superoxide pro-oxidants. This last mechanism, called respiratory burst, occurs in vertebrate macrophages and insect haemocytes during the defence against parasites [40, 41]. The pentose phosphate pathway is known to play a role in host-protozoa interaction [42], for instance a G6PDH mutation confers resistance to malaria in human populations [43]. We expected a trade-off between fight against the parasite and fecundity, with exposed females investing more in immune and antioxidant defences but laying less than unexposed females. We also expected the effect of exposure on investment in defences and on fecundity to increase with the exposure level.

Methods

Experimental design

Experimental infection of canaries

We used three different groups of birds to expose female mosquitoes to Plasmodium parasites. The first bird group was constituted with four uninfected canaries (uninfected group). Female mosquitoes feeding on these birds were not exposed to parasites and composed the control (unexposed) group. The second bird group consisted of four canaries experimentally infected via an intraperitoneal injection of 200 μl of blood from wild great tits (Parus major) naturally infected by Plasmodium relictum mixed with phosphate buffer saline (PBS, 1:1) two years before the beginning of the experiment. These canaries passed the acute phase of infection (peak of parasitaemia) and were in the chronic phase of infection [44], they had then low parasitaemia (chronic group). Mosquitoes feeding on them were then assumed to be exposed to a low amount of Plasmodium blood stages (low exposure). The third bird group consisted of four canaries experimentally infected via an intraperitoneal injection of 75 μl of a blood mix prepared from the four canaries of the chronic group mixed with PBS (1:1). In order to ensure a high parasitaemia in this last bird group, experimental infections were performed 10 days before the first mosquito feeding session to coincide with the acute phase of infection (acute group). Mosquitoes feeding on these birds were assumed to be exposed to a high amount of Plasmodium blood stages (high exposure).

We performed eight mosquito feeding sessions, using one canary of each bird group per feeding session, and using each bird twice. Prior to the first feeding session and once all the feeding sessions were performed, we blood sampled the birds in order to determine bird haematocrit (i.e. the fraction of red blood cells in the total blood volume) and bird parasitaemia in infected ones. Canaries from the acute group had significantly higher parasitaemia than canaries from the chronic group (Wilcoxon rank test: W = 0, P = 0.029; Additional file 1: Figure S1a). Mean haematocrit did not differ between the three bird groups (ANOVA: F (2,9) = 1.08, P = 0.381, Additional file 1: Figure S1b).

Mosquito rearing and exposure to Plasmodium

We set up rainfall collecting containers (50 × 30 × 25 cm) in July 2014 at the forest of Dorigny (46°31′N, 6°34′E; alt. 380 m). Containers were initially filled up with water from Lake Geneva and baited with live yeast to favour visitation by gravid Culex pipiens females. From the 19th of August to the 10th of September, we collected Culex pipiens egg rafts in order to create independent clutch groups: eight groups of 25 ± 4 clutches (mean ± standard deviation). Clutches were allowed to hatch in the laboratory (24 °C, 65% relative humidity and a 12:12 h light–dark cycle) in individual containers filled up with 250 ml of mineral water. Larvae were then fed with commercial fish flakes (Tetra). The first emerging males of each clutch were morphologically determined to the species level and only clutches confirmed as Culex pipiens species were used. For each clutch group, seven days after the first emergence, adults were collected during 5 days, pooled in common cages (30 × 30 × 90 cm) and were provided with a 10% glucose solution. This allowed creating eight independent groups (each of them coming from one clutch group aforementioned and later called emergent group) of known age mosquitoes (age ± 2 days). Once an emergent group was created, females were kept 16 more days in cages with males for mating.

When females were 19 ± 2 days old, they were allowed to feed on one of the three different bird groups. For each emergent group, pools of 20 females were kept in cages (30 × 30 × 30 cm) 24 h before the feeding session and provided with water only. At the end of the 24 h, a first batch of three cages received one bird of each bird group (1 canary per cage). After 30 min, canaries were removed and placed in a second batch of three cages. This scenario was repeated to obtain as far as possible a total of 20 fed mosquitoes per bird over all the visited cages. Directly after each feeding session, fed females (presence of blood in the abdomen determined by eye) were placed separately in empty plastic vials (Sartsdet, 30 ml) provided with a 10% glucose solution during four days to collect haematin secretion for the estimation of the blood meal size. Haematin was stored at -80 °C until laboratory analyses and quantified via spectrophotometric measurement following previously described methods [45, 46]. Blood meal size was square root transformed to achieve normality.

At 5 days post-feeding (dpf), females were transferred into new vials containing 2.5 ml of water in order to allow them to oviposit and were fed with water only. Egg rafts were collected in order to measure clutch size as a proxy of fecundity. Females were then transferred into a third vial and provided with water only until death. We imposed a nutritional stress, a common stress factor in the wild, to the females because previous studies have shown that some effects of Plasmodium infection on vector life history traits might be stronger under stressful conditions [3]. All females were daily checked from the fifth dpf until death to record survival.

In order to quantify immune and antioxidant gene expression, ten females per exposure level and per feeding session were sacrificed at 5 and at 15 dpf (whenever possible five sacrificed females per dpf). To minimize variation between individuals, females were randomly chosen among females having already oviposited. Unexposed and exposed groups differed in the identity of the birds on which they fed, which can be an important confounding factor. We therefore compared these groups to a reference group. For this purpose, eight female mosquitoes of each emergent group were isolated before the feeding session and then sacrificed in order to assess the gene expression level prior to blood meal (reference baseline). Sacrificed and naturally dead females were stored at -80 °C until laboratory analyses. After mosquito dissection, wing size was measured as described in [47] and used as a proxy of body size.

We hypothesised that the complete development of the parasite in the vector might be associated with particular physiological changes in the mosquito vector (different intensities or different mechanisms compared to an aborted development). As we were interested by the cost of exposure when the parasite does not complete its life cycle, females which developed sporozoites in the salivary glands were excluded from the analyses.

Laboratory measurements

In order to assess both infection status and investment in immune and antioxidant defences in a single individual, we dissected the mosquitoes and separated thorax and abdomen. This allowed us to extract DNA from the thoracic part for Plasmodium sporozoite detection in the salivary gland by PCR and to extract RNA from the abdominal part for gene expression quantification.

Plasmodium detection and quantification

Mosquito DNA was extracted from the thorax and bird DNA was extracted from the blood using the DNeasy blood and tissue extraction kit (Qiagen, Hilden, Germany), according to the manufacturer’s protocol. Plasmodium parasites were detected using PCR method previously developed by [48] and parasite quantification (in infected birds) was performed using quantitative real-time PCR as described in [49]. Parasitaemia was log-transformed to achieve normality.

Immune and antioxidant gene expression quantification

We measured the expression of one gene involved in mosquito immune response (nitric oxide synthase, NOS), one gene involved in antioxidant process (superoxide dismutase, SOD) and one gene involved in both immune and antioxidant pathways (glucose 6-phosphate dehydrogenase, G6PDH). Total RNA was extracted from the mosquito abdomen using Trizol reagent (Life Technologies, Invitrogen, Carlsbad, California, USA) and treated with Deoxyribonuclease I, Amplification Grade (DNase I, Amp Grade, Invitrogen, Carlsbad, California, USA), according to the manufacturers’ protocols. Total RNA was quantified using NanoDrop spectrophotometer (Thermo Fisher Scientific, Waltham, Massachusetts, USA). Retrotranscription of RNA to cDNA was performed from 50 ng of total RNA using PrimeScript RT reagent kit (Takara Bio, Kusatsu, Japan), following the manufacturer’s protocol. The qPCR reaction was performed using PowerUp SYBR Green Master Mix (Thermo Fisher Scientific) according to the manufacturer’s protocol and cDNA samples diluted to 1:10. The primer sequences used were: NOS, 5′-CGA GAA GGC CCA CAT CTA CG-3′ and 5′-CGA CAG CAT GTA CTT CTC CA-3′; SOD, 5′-GCA TTG CGA AAA CTT CCT TC-3′ and 5′-TGC CCA GAT CAT CAA TTT CA-3′; G6PDH, 5′-CGC GCA CGA GGA AAA GTA CG-3′ and 5′-GGT TTG CGG TCT TCC CAA CC-3′ and rpl19, 5′-CGC TTT GTT TGA TCG TGT GT-3′ and 5′-CCA ATC CAG GAG TGC TTT TG-3′ as reference gene [36, 50, 51]. The qPCRs were run in duplicates with an AB 7500 qPCR system (Applied Biosystems, Foster City, California, USA) as follows: 2 min at 50 °C, 10 min at 95 °C followed by 40 cycles: 15 s at 95 °C, 1 min at 60 °C. A melting curve was produced at the end of each run to control for amplification specificity. A standard curve was produced for each gene with 2-fold serial dilutions of cDNA. Expression of each gene of interest was calculated the same way that bird parasitaemia and as described in [49]. Briefly, for each measured gene, DNA concentration was calculated from a standard curve and the gene expression level was given by the ratio of the DNA concentration of the gene of interest on the DNA concentration of the reference gene. Expression level of each gene was log-transformed to achieve normality. One female, which was sacrificed at 5 dpf and which received high exposure, was discarded from the analyses because of particularly high gene expression values (outsider).

Statistical analyses

Statistical analyses were performed with R (version 3.1, [52]). We had two measurements per bird for parasitaemia (before the first feeding session and at the end of all the feeding sessions), then we calculated the mean values for each individual and used mean bird parasitaemia in the analyses.

We analysed blood meal size as a response variable in a linear mixed effects model (lmer function in lme4 package) including terms for body size, exposure level (unexposed - low exposure - high exposure) and the interaction between both. For the fecundity, we first analysed the laying probability as a response variable in a generalized mixed effects model (glmer function in lme4 package) including terms for blood meal size, exposure level and the interaction between both. Considering mosquitoes which oviposited only, we also analysed clutch size as a response variable in a lmer model including terms for body size, blood meal size, exposure level and the interactions between body size and exposure level and between blood meal size and exposure level. Considering all females which died naturally, we analysed survival as a response variable in mixed effects cox model (coxme function in coxme package) including terms for blood meal size, fecundity, exposure level and the interactions between blood meal size and exposure level and between fecundity and exposure level. For all models, we included terms for emergent group and canary identity as random factors. We also ran the models considering low and high exposure levels only, replacing exposure level by bird parasitaemia and removing canary identity from the random factors to investigate the effect of exposure level as a quantitative continuous variable. The detailed model structures are described in Additional file 2: Table S1.

For each gene measured, we first compared the gene expression level of each group of fed females (5 no - 5 low - 5 high - 15 no - 15 low - 15 high, dpf exposure) to the gene expression level of sacrificed unfed females with a one-way Anova and emergent group as a random factor (lme function in nlme package). Finally, in order to look at the effect of gene expression on fecundity, we considered exposed females sacrificed at 15 dpf and analysed laying probability and clutch size as response variables and included terms for expression levels of each gene and exposure level as well as the two-way interactions between each gene expression level and exposure level in lmer and glmer models (for clutch size and laying probability respectively) and considering emergent group and canary identity as random factors. The detailed model structures are described in Additional file 2: Table S1.

For each model, we performed backward model selection based on likelihood ratio tests until reaching the minimal adequate model [53]. After model selection, we performed contrast analyses to look at the effect of each term individually. The significant P-values given in the text come from the minimal adequate models and the non-significant P-values come from the likelihood ratio tests prior to the elimination of the non-significant term from the model.

Results

Mosquito exposure

At the end of all the feeding sessions, 156 female mosquitoes took a blood meal on uninfected canaries, 199 took a blood meal on canaries from the chronic group and 162 took a blood meal on canaries from the acute group. Details about sample sizes are given in Table 1.

Table 1 Sample sizes of female (♀) mosquitoes obtained at the end of all the feeding sessions

Blood meal size, fecundity and survival

Blood meal size increased with mosquito body size (Table 2). In exposed females, blood meal size slightly increased with bird parasitaemia (Table 2).

Table 2 Minimal adequate models for blood meal size, fecundity and survival

Laying probability increased with blood meal size (Table 2, Fig. 1) and tended to decrease with increasing exposure level (P = 0.064): females under high exposure had lower probability to lay than females under low exposure (P = 0.010) but there was a marginally significant difference between females under high exposure and unexposed females due to higher variation in the unexposed group (P = 0.066, Table 2, Fig. 2a). The effect of exposure level on laying probability was clearer when considering exposure level as a quantitative continuous variable and in exposed females, laying probability decreased with bird parasitaemia (Table 2, Fig. 2b).

Fig. 1
figure 1

Laying probability as a function of blood meal size. Fecundity measured as laying probability (in percent) of females that died naturally as a function of blood meal size from the minimal adequate model given in Table 2. Circle size represents mosquito sample size (n)

Fig. 2
figure 2

Laying probability as a function of exposure to Plasmodium parasites. Fecundity measured as laying probability (in percent) as a function of a exposure level in females that died naturally, and b bird parasitaemia (circle size represents mosquito sample size) in exposed females that died naturally from the minimal adequate model given in Table 2

Clutch size increased with body size (Table 2). There was also a significant interaction between blood meal size and exposure level to Plasmodium (Table 2): clutch size increased with blood meal size and for similar blood meal size exposed females laid more eggs compared to unexposed ones (Table 2, Fig. 3). This difference of blood-to-egg conversion increased with blood meal size (Fig. 3). There was no effect of bird parasitaemia on clutch size (Table 2).

Fig. 3
figure 3

Clutch size as a function of blood meal size and exposure level. Fecundity measured as clutch size in females that oviposited as a function of blood meal size and exposure level from the minimal adequate model given in Table 2

Only clutch size significantly explained survival and females laying more eggs died faster (Table 2).

Gene expression

Compared to unfed female mosquitoes, NOS level was significantly higher in exposed females at 5 dpf (low exposure: P = 0.021, high exposure: P = 0.027, Fig. 4a). NOS levels in the other groups were not significantly different from unfed females. All fed females had increased SOD level with a 10-fold increase at 5 dpf (no – low – high exposure at 5 dpf: all P < 0.001) and a 5-fold increase at 15 dpf (no - low exposure at 15 dpf: both P = 0.002, high exposure at 15 dpf: P < 0.001) compared to unfed females (Fig. 4b). Compared to unfed females, only highly exposed females at 5 dpf had a higher G6PDH level (Fig. 4c). None of the gene expression levels were correlated with body size, blood meal size or bird parasitaemia (in exposed mosquitoes), but SOD level was increasing with NOS and with G6PDH levels (Table 3).

Fig. 4
figure 4

Immune and antioxidant defences as a function of exposure level and days post- feeding. Gene expression levels of a nitric oxide synthase (NOS, log-transformed), b superoxide dismutase (SOD, log-transformed) and c glucose 6-phosphate dehydrogenase (G6PDH, log-transformed) in unfed female mosquitoes and in fed female mosquitoes: unexposed, with low exposure and with high exposure to Plasmodium parasites at 5 and 15 days post-feeding (dpf). Stars indicate significant statistical differences compared to unfed females at 0 dpf

Table 3 Pearson correlations

In exposed females sacrificed at 15 dpf, there was a marginal effect of the interaction between SOD and exposure level on laying probability. Laying probability decreased with SOD level and the decrease was more important in highly exposed females compared to females receiving a low exposure (Table 4, Fig. 5a). There was no effect of NOS or G6PDH levels (Table 4). Among females that laid eggs, clutch size increased with NOS level (Table 4, Fig. 5b) but there was no effect of SOD and G6PDH expressions or exposure level (Table 4).

Table 4 Minimal adequate models for fecundity as a function of gene expression
Fig. 5
figure 5

Fecundity and gene expression in exposed females sacrificed at 15 days post-feeding. a Laying probability as a function of superoxide dismutase level (SOD, log-transformed) and exposure level from the minimal adequate model given in Table 4. Circle size represents mosquito sample size. b Clutch size as a function of nitric oxide synthase level (NOS, log-transformed) from the minimal adequate model given in Table 4

Discussion

In this study, we investigated the effect of exposure to Plasmodium parasites on the mosquito vector Culex pipiens fecundity, survival and investment in immune and antioxidant defences. Firstly, we found a cost of exposure to Plasmodium parasites on the vector fecundity as laying probability decreased with increasing exposure level (measured as bird parasitaemia). Laying probability also decreased with increasing SOD expression in exposed females. This effect suggested a trade-off between fight against the parasite/self-maintenance and reproduction. Secondly, among females that laid eggs, exposed females had a higher blood-to-egg conversion compared to unexposed ones and their clutch size also increased with NOS expression. These results depicted a link between investment in immune defences and fecundity. Thirdly, only highly exposed females showed a higher expression of G6PDH, suggesting the involvement of that gene during a vector-Plasmodium interaction. Finally, we found no effect of exposure to Plasmodium on the vector survival.

Plasmodium parasites are known to affect vector fecundity with, for instance, egg resorption occurring in the mosquito ovaries during Plasmodium development [11, 2932]. Here, we found that laying probability increased with blood meal size: a particularly small blood meal being not sufficient to mature eggs [54]. This taken into account, we found that exposure to Plasmodium imposed a cost on the vector fecundity as laying probability decreased with increasing exposure to Plasmodium (measured as bird parasitaemia) as well as with increasing SOD expression. All fed females, regardless of their exposure level, expressed similarly SOD showing no specific investment in this antioxidant defence after Plasmodium exposure but simply a protective role after blood intake [55]. All exposed females, regardless of the low or high exposure, expressed more NOS. Although uninfected blood ingestion may trigger immune activation in the vector due to the presence of immunogenic factors in the vertebrate blood, a stronger activation as well as the activation of specific genes after a Plasmodium-infected blood meal have been described [610]. Our results suggested a higher oxidative stress in exposed females compared to unexposed ones (i.e. same amount of antioxidant, but more pro-oxidants in exposed females compared to unexposed ones). Finally, only highly exposed females were expressing more G6PDH at 5 dpf. It is the first time that this gene is shown to be involved during a vector-Plasmodium interaction. This result suggested that the activation of some physiological pathways might depend on the exposure level and occur above a certain threshold of exposure. G6PDH is an enzyme involved in the pentose phosphate pathway allowing the transformation of NADP+ in NADPH [39]. It is involved in both immune and antioxidant pathways through its implication in NO, O2 and glutathione production. Our data do not allow to discriminate between these processes. Yet, the laying probability of exposed females decreased with an increasing exposure level, an increasing SOD level and with the expression of G6PDH. All together, these results corroborate a trade-off between fight against parasites/self-maintenance and reproduction associated with an investment in immune and/or antioxidant defences that may mediate the trade-off.

Previous studies [11, 2932] have been looking at mosquito fecundity via dissection of the ovaries without allowing females to oviposit. Some studies measured fecundity as the clutch size after oviposition and observed a decreased clutch size with infected blood meal intake [33]. Here, exposed females that oviposited had a higher blood-to-egg conversion compared to unexposed ones. A positive effect of Plasmodium exposure on the vector fecundity has already been reported in [56] where the prevalence of gravid females (females producing eggs in the ovaries) was higher in the group taking an infected blood meal than in the group taking an uninfected blood meal. If this result would be due to different blood quality (e.g. amount of resources) between non-infected and infected blood, we could expect a fixed conversion difference between unexposed and exposed females. This was not the case here as the conversion difference increased with blood meal size. Alternatively, it has also been shown in Culex pipiens that parasitism influences reproductive life-history traits with for instance delayed hatching date of eggs laid after a Plasmodium-infected blood meal [33] or larvae pupating faster and reproducing earlier to compensate for decreased survival when infected by microsporidian parasites [57]. In our case, a higher blood-to-egg conversion might either compensate a potential impaired hatching success or, for similar hatching rate, enhance reproductive success in exposed females compared to unexposed ones. This could be part of a fight against parasites-reproduction trade-off. Indeed, the first gonotrophic cycle is important for female mosquitoes, because several risks increase with subsequent reproduction events: risk of death at each blood meal intake [5860], decreased fertility with age [37] and with infection [29, 61]. Therefore, maximizing the first reproduction when taking an infected blood meal may be beneficial. Finally, the positive link observed between clutch size and NOS expression could either reflect the condition status of the females: the better the female, the more eggs and the higher the immune response, or suggests that investing in immune defences could act as a signal for investing in egg maturation in a proportional way. Even though our data do not provide a terminal measurement of reproductive success (e.g. hatching rate) of the females and do not allow us to conclude about the underlying mechanism, these results showed a positive link between investment in immune defences and fecundity (measured as clutch size). This positive link contrasts with the former result showing a decreased fecundity (measured as laying probability) with an increased exposure to Plasmodium. This dual response could depend on physiological signals perceived by individuals. For instance, the amount of ingested parasites, the strength of a particular gene activation (or downregulation) or the engagement in a particular physiological pathway (e.g. G6PDH more expressed in highly exposed females only) may modify the vector physiological status. Such modifications could act as cues and contribute to determine if females will invest in self-maintenance or in reproduction. We could thus expect that the stronger the cues the female is receiving, the higher the investment in reproduction and the higher the probability of developing Plasmodium infection reaching the salivary glands. However, several factors may affect the probability of infection and lead to very contrasting results. For instance, if infection depends on individual quality, we would expect infected individuals to present a lower investment in immune and antioxidant defences and a lower reproductive success. In contrast, if infection depends on the number of ingested blood parasites only, we would expect no relationship between life history traits and infection. Supplementary data on Plasmodium-infected female mosquitoes could help understanding the involved underlying mechanisms.

While some of our present findings were in the direction of the predictions and supported previous studies, others were counterintuitive. The effect of Plasmodium parasites on vector life history traits will depend on intrinsic characteristics of the host-parasite association under study. Such characteristics are the interspecific interaction (what are the interacting species?), the nature of the association (i.e. natural versus artificial) and the time since the (host and parasite) populations have been isolated from the wild. Indeed, naturally occurring associations share an evolutionary history that has led the interacting species to express specific exploitation and defence strategies. These strategies will depend on the interacting species and might be absent in artificial associations. Furthermore, the isolation of individuals from the wild to controlled laboratory conditions for several generations alters the selective pressures imposed on them and may modify the expression of wild traits. In research conducted on vector-Plasmodium, studies may differ substantially in terms of these characteristics [62]. Here, we studied a natural (Culex pipiens-Plasmodium relictum) and wild (mosquitoes hatched from clutches collected on the field and Plasmodium strain was isolated since two years) vector-parasite association. The use of a natural vertebrate host (instead of Serinus canaria, an artificial host used here for practical reasons) would allow to study a fully natural host-vector-parasite association. As already stressed in the literature [5], there are important factors to consider that may explain the existence of conflicting results between studies.

Conclusions

In this study, we showed a dual response after the ingestion of a Plasmodium-infected blood meal in the mosquito vector. Following exposure, an allocation trade-off occurred between self-maintenance and reproduction, which tended to lead either to investing in defences but not reproducing or to reproducing more efficiently under infection risk. We also showed an involvement of immune and/or antioxidant defences that may act as mediators. These results illustrate a complex effect of exposure to Plasmodium on vector life history traits. The understanding of the mechanisms underlying the variation between individuals, but also between studies, is particularly important for malaria epidemiology.

Abbreviations

cDNA:

Complementary deoxyribonucleic acid

DNA:

Deoxyribonucleic acid

dpf:

Days post-feeding

G6PDH:

Glucose-6-phosphate dehydrogenase

NADP:

Nicotinamide adénine dinucléotide phosphate

NO:

Nitric oxide

NOS:

Nitric oxide synthase

NOX:

NADPH-oxidase

O2 :

Superoxide

PBS:

Phosphate buffer saline

PCR:

Polymerase chain reaction

qPCR:

Quantitative polymerase chain reaction

RNA:

Ribonucleic acid

SOD:

Superoxide dismutase

References

  1. Koella JC, Rieu L, Paul REL. Stage-specific manipulation of a mosquito’s host-seeking behavior by the malaria parasite Plasmodium gallinaceum. Behav Ecol. 2002;13:816–20.

    Article  Google Scholar 

  2. Ferguson HM, Read AF. Mosquito appetite for blood is stimulated by Plasmodium chabaudi infections in themselves and their vertebrate hosts. Malar J. 2004;3:1–8.

    Article  Google Scholar 

  3. Lalubin F, Delédevant A, Glaizot O, Christe P. Natural malaria infection reduces starvation resistance of nutritionally stressed mosquitoes. J Anim Ecol. 2014;83:850–7.

    Article  PubMed  Google Scholar 

  4. Rossignol PA, Ribeiro JM, Spielman A. Increased biting rate and reduced fertility in sporozoite-infected mosquitoes. Am J Trop Med Hyg. 1986;35:277–9.

    CAS  PubMed  Google Scholar 

  5. Ferguson HM, Read AF. Why is the effect of malaria parasites on mosquito survival still unresolved? Trends Parasitol. 2002;18:256–61.

    Article  PubMed  Google Scholar 

  6. Dimopoulos G, Richman A, Müller H-M, Kafatos FC. Molecular immune responses of the mosquito Anopheles gambiae to bacteria and malaria parasites. Proc Natl Acad Sci USA. 1997;94:11508–13.

  7. Dimopoulos G, Seeley D, Wolf A, Kafatos FC. Malaria infection of the mosquito Anopheles gambiae activates immune-responsive genes during critical transition stages of the parasite life cycle. EMBO J. 1998;17:6115–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Oduol F, Xu J, Niare O, Natarajan R, Vernick KD. Genes identified by an expression screen of the vector mosquito Anopheles gambiae display differential molecular immune response to malaria parasites and bacteria. Proc Natl Acad Sci USA. 2000;97:11397–402.

  9. Dimopoulos G, Christophides GK, Meister S, Schultz J, White KP, Barillas-Mury C, Kafatos FC. Genome expression analysis of Anopheles gambiae: responses to injury, bacterial challenge, and malaria infection. Proc Natl Acad Sci USA. 2002;99:8814–9.

  10. Dong Y, Manfredini F, Dimopoulos G. Implication of the mosquito midgut microbiota in the defense against malaria parasites. PLoS Pathog. 2009;5:e1000423.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Gray EM, Bradley TJ. Malarial infection in Aedes aegypti: effects on feeding, fecundity and metabolic rate. Parasitology. 2006;132:169–76.

    Article  CAS  PubMed  Google Scholar 

  12. Zhao YO, Kurscheid S, Zhang Y, Liu L, Zhang L, Loeliger K, Fikrig E. Enhanced survival of Plasmodium-infected mosquitoes during starvation. PLoS One. 2012;7:e40556.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. de Almeida OG, Lieberman J, Barillas-mury C. Epithelial nitration by a peroxidase/NOX5 system mediates mosquito antiplasmodial immunity. Science. 2012;335:856–9.

    Article  Google Scholar 

  14. Gonçalves RLS, Oliveira JHM, Oliveira GA, Andersen JF, Oliveira MF, Oliveira PL, Barillas-Mury C. Mitochondrial reactive oxygen species modulate mosquito susceptibility to Plasmodium infection. PLoS One. 2012;7:e41083.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Saraiva RG, Kang S, Simões ML, Angleró-Rodríguez YI, Dimopoulos G. Mosquito gut antiparasitic and antiviral immunity. Dev Comp Immunol. 2016;in press.

  16. Blandin S, Levashina EA. Mosquito immune responses against malaria parasites. Curr Opin Immunol. 2004;16:16–20.

    Article  CAS  PubMed  Google Scholar 

  17. Michel K, Kafatos FC. Mosquito immunity against Plasmodium. Insect Biochem Mol Biol. 2005;35:677–89.

    Article  CAS  PubMed  Google Scholar 

  18. Dimopoulos G. Insect immunity and its implication in mosquito - malaria. Cell Microbiol. 2003;5:3–14.

  19. Whitten MMA, Shiao SH, Levashina EA. Mosquito midguts and malaria: cell biology, compartmentalization and immunology. Parasite Immunol. 2006;28:121–30.

    Article  CAS  PubMed  Google Scholar 

  20. Crompton PD, Moebius J, Portugal S, Waisberg M, Hart G, Garver LS, et al. Malaria immunity in man and mosquito: insights into unsolved mysteries of a deadly infectious disease. Annu Rev Immunol. 2014;32:157–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Rivero A. Nitric oxide: an antiparasitic molecule of invertebrates. Trends Parasitol. 2006;22:219–25.

    Article  CAS  PubMed  Google Scholar 

  22. Luckhart S, Vodovotz Y, Cui L, Rosenberg R. The mosquito Anopheles stephensi limits malaria parasite development with inducible synthesis of nitric oxide. Proc Natl Acad Sci USA. 1998;95:5700–5.

  23. Kumar S, Christophides GK, Cantera R, Charles B, Han YS, Meister S, et al. The role of reactive oxygen species on Plasmodium melanotic encapsulation in Anopheles gambiae. Proc Natl Acad Sci USA. 2003;100:14139–44.

  24. Kumar S, Gupta L, Yeon SH, Barillas-Mury C. Inducible peroxidases mediate nitration of Anopheles midgut cells undergoing apoptosis in response to Plasmodium invasion. J Biol Chem. 2004;279:53475–82.

    Article  CAS  PubMed  Google Scholar 

  25. Sadd BM, Siva-Jothy MT. Self-harm caused by an insect’s innate immunity. Proc R Soc B Biol Sci. 2006;273:2571–4.

    Article  Google Scholar 

  26. Bahia AC, Oliveira JHM, Kubota MS, Araújo HRC, Lima JBP, Ríos-Velásquez CM, et al. The role of reactive oxygen species in Anopheles aquasalis response to Plasmodium vivax infection. PLoS One. 2013;8:e57014.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Molina-Cruz A, DeJong RJ, Charles B, Gupta L, Kumar S, Jaramillo-Gutierrez G, Barillas-Mury C. Reactive oxygen species modulate Anopheles gambiae immunity against bacteria and Plasmodium. J Biol Chem. 2008;283:3217–23.

    Article  CAS  PubMed  Google Scholar 

  28. Kumar S, Barillas-Mury C. Ookinete-induced midgut peroxidases detonate the time bomb in anopheline mosquitoes. Insect Biochem Mol Biol. 2005;35:721–7.

    Article  CAS  PubMed  Google Scholar 

  29. Hogg JC, Carwardine S, Hurd H. The effect of Plasmodium yoelii nigeriensis infection on ovarian protein accumulation by Anopheles stephensi. Parasitol Res. 1997;83:374–9.

    Article  CAS  PubMed  Google Scholar 

  30. Hopwood JA, Ahmed AM, Polwart A, Williams GT, Hurd H. Malaria-induced apoptosis in mosquito ovaries: a mechanism to control vector egg production. J Exp Biol. 2001;204:2773–80.

    CAS  PubMed  Google Scholar 

  31. Ahmed AM, Hurd H. Immune stimulation and malaria infection impose reproductive costs in Anopheles gambiae via follicular apoptosis. Microbes Infect. 2006;8:308–15.

    Article  CAS  PubMed  Google Scholar 

  32. Hogg JC, Hurd H. The effects of natural Plasmodium falciparum infection on the fecundity and mortality of Anopheles gambiae s. l. in north east Tanzania. Parasitology. 1997;114:325–31.

    Article  PubMed  Google Scholar 

  33. Vézilier J, Nicot A, Gandon S, Rivero A. Plasmodium infection decreases fecundity and increases survival of mosquitoes. Proc R Soc B Biol Sci. 2012;279:4033–41.

    Article  Google Scholar 

  34. Schwenke RA, Lazzaro BP, Wolfner MF. Reproduction-immunity trade-offs in insects. Annu Rev Entomol. 2016;61:239–56.

    Article  CAS  PubMed  Google Scholar 

  35. Uchida K. Balanced amino acid composition essential for infusion-induced egg development in the mosquito (Culex pipiens pallens). J Insect Physiol. 1993;39:615–21.

    Article  CAS  Google Scholar 

  36. Sim C, Denlinger DL. Catalase and superoxide dismutase-2 enhance survival and protect ovaries during overwintering diapause in the mosquito Culex pipiens. J Insect Physiol. 2011;57:628–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. DeJong RJ, Miller LM, Molina-Cruz A, Gupta L, Kumar S, Barillas-Mury C. Reactive oxygen species detoxification by catalase is a major determinant of fecundity in the mosquito Anopheles gambiae. Proc Natl Acad Sci USA. 2007;104:2121–6.

  38. Vézilier J, Nicot A, Gandon S, Rivero A. Insecticide resistance and malaria transmission: infection rate and oocyst burden in Culex pipiens mosquitoes infected with Plasmodium relictum. Malar J. 2010;9:379.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Stanton RC. Glucose-6-phosphate dehydrogenase, NADPH, and cell survival. IUBMB Life. 2012;64:362–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Lavine MD, Strand MR. Insect hemocytes and their role in immunity. Insect Biochem Mol Biol. 2002;32:1295–309.

    Article  CAS  PubMed  Google Scholar 

  41. Strand MR. The insect cellular immune response. Insect Sci. 2008;15:1–14.

    Article  CAS  Google Scholar 

  42. Barrett MP. The pentose phosphate pathway and parasitic protozoa. Parasitol Today. 1997;13:11–6.

    Article  CAS  PubMed  Google Scholar 

  43. Ruwende C, Hill A. Glucose-6-phosphate dehydrogenase deficiency and malaria. J Mol Med. 1998;76:581–8.

    Article  CAS  PubMed  Google Scholar 

  44. Valkiunas G. Avian Malaria Parasites and Other Haemosporidia. Boca Ratón: CRC Press; 2005.

    Google Scholar 

  45. Rivero A, Ferguson HM. The energetic budget of Anopheles stephensi infected with Plasmodium chabaudi: is energy depletion a mechanism for virulence? Proc R Soc B Biol Sci. 2003;270:1365–71.

    Article  CAS  Google Scholar 

  46. Briegel H. Determination of uric acid and hematin in a single sample of excreta from blood-fed insects. Experientia. 1980;36:1428.

    Article  CAS  Google Scholar 

  47. Mpho M, Holloway GJ, Callaghan A. Fluctuating wing asymmetry and larval density stress in Culex quinquefasciatus (Diptera: Culicidae). Bull Entomol Res. 2000;90:279–83.

    Article  CAS  PubMed  Google Scholar 

  48. Waldenström J, Bensch S, Hasselquist D, Östman Ö. A new nested polymerase chain reaction method very efficient in detecting Plasmodium and Haemoproteus infections from avian blood. J Parasitol. 2004;90:191–4.

    Article  PubMed  Google Scholar 

  49. Jenkins T, Delhaye J, Christe P. Testing local adaptation in a natural great tit-malaria system: an experimental approach. PLoS One. 2015;10:e0141391.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Vézilier J, Nicot A, Lorgeril J, Gandon S, Rivero A. The impact of insecticide resistance on Culex pipiens immunity. Evol Appl. 2013;6:497–509.

    Article  PubMed  Google Scholar 

  51. Sim C, Denlinger DL. Transcription profiling and regulation of fat metabolism genes in diapausing adults of the mosquito Culex pipiens. Physiol Genomics. 2009;39:202–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Development Core Team. R: A Language and Environment for Statistical Computing. 2011. R Foundation for Statistical Computing.

    Google Scholar 

  53. Crawley MJ. The R, Book. Chichester: Wiley; 2007.

    Book  Google Scholar 

  54. Clements AN, Boocock MR. Ovarian development in mosquitoes: stages of growth and arrest, and follicular resorption. Physiol Entomol. 1984;9:1–8.

    Article  Google Scholar 

  55. Graça-Souza AV, Maya-Monteiro C, Paiva-Silva GO, Braz GRC, Paes MC, Sorgine MHF, et al. Adaptations against heme toxicity in blood-feeding arthropods. Insect Biochem Mol Biol. 2006;36:322–35.

    Article  PubMed  Google Scholar 

  56. Ferguson HM, Gouagna LC, Obare P, Read AF, Babiker H, Githure J, Beier JC. The presence of Plasmodium falciparum gametocytes in human blood increases the gravidity of Anopheles gambiae mosquitoes. Am J Trop Med Hyg. 2005;73:312–20.

    PubMed  Google Scholar 

  57. Agnew P, Bedhomme S, Haussy C, Michalakis Y. Age and size at maturity of the mosquito Culex pipiens infected by the microsporidian parasite Vavraia culicis. Proc R Soc B Biol Sci. 1999;266:947–52.

    Article  Google Scholar 

  58. Roitberg BD, Mondor EB, Tyerman JGA. Pouncing spider, flying mosquito: blood acquisition increases predation risk in mosquitoes. Behav Ecol. 2003;14:736–40.

    Article  Google Scholar 

  59. Walker ED, Edman JD. The influence of host defensive behavior on mosquito (Diptera: Culicidae) biting persistence. J Med Entomol. 1985;22:370–2.

    Article  CAS  PubMed  Google Scholar 

  60. Reid JN, Hoffmeister TS, Hoi AG, Roitberg BD. Bite or flight: the response of mosquitoes to disturbance while feeding on a defensive host. Entomol Exp Appl. 2014;153:240–5.

    Article  Google Scholar 

  61. Hogg JC, Hurd H. Plasmodium yoelii nigeriensis: the effect of high and low intensity of infection upon the egg production and bloodmeal size of Anopheles stephensi during three gonotrophic cycles. Parasitology. 1995;111:555–62.

    Article  PubMed  Google Scholar 

  62. Tripet F. Ecological immunology of mosquito-malaria interactions: of non-natural versus natural model systems and their inferences. Parasitology. 2009;136:1935–42.

    Article  CAS  PubMed  Google Scholar 

  63. Nakagawa S, Schielzeth H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol Evol. 2013;4:133–42.

    Article  Google Scholar 

Download references

Acknowledgements

We are thankful to Anne-Lyse Ducrest for her advice on laboratory protocols and Laura Bonny for her laboratory assistance.

Funding

The project was funded by the Swiss National Science Foundation (grant 31003A-138187 and 31003A-159600).

Availability of data and materials

The datasets supporting the main conclusions of this study are included within the manuscript and its additional files. Individual data are available from the corresponding author on reasonable request.

Authors’ contributions

Conceived and designed the experiment: JD PC OG. Performed the experiment: JD CA. Analysed the data: JD. Wrote the manuscript: JD PC OG. All authors read and approved the final version of the manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval

This study was conducted under license according to the Swiss animal legislation (authorization number 1730.1).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Jessica Delhaye.

Additional files

Additional file 1: Figure S1.

Bird parasitaemia and haematocrit. a Parasitaemia (log-transformed, arbitrary unit) before the first feeding session, at the end of all the feeding sessions and mean value for canaries in the chronic group (1 to 4) and in the acute group (5 to 8). Mean values significantly differed between the two groups (Wilcoxon rank test: W = 0, P = 0.029). b Haematocrit (fraction of red blood cells in the total blood volume) before the first feeding session, at the end of all the feeding sessions and mean value for canaries in the chronic group (1 to 4), in the acute group (5 to 8) and in the uninfected group (9 to 12). Mean values did not significantly differ between the three groups (F (2,9) = 1.08, P = 0.381). (PDF 183 kb)

Additional file 2: Table S1.

Detailed model structures and initial full models prior to model selection. For each response variable (blood meal size, clutch size, survival and laying probability), details are given about which females (♀) have been considered, the model type as well as the fixed and the random factors considered. a for each response variable, the same models were also ran considering females under low and high exposure only, replacing exposure level by bird parasitaemia and removing canary identity from the random factors. b Abbreviationsdpf: days post-feeding, NOS: nitric oxide synthase, SOD: superoxide dismutase, G6PDH: glucose-6-phosphate dehydrogenase. (DOCX 13 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Delhaye, J., Aletti, C., Glaizot, O. et al. Exposure of the mosquito vector Culex pipiens to the malaria parasite Plasmodium relictum: effect of infected blood intake on immune and antioxidant defences, fecundity and survival. Parasites Vectors 9, 616 (2016). https://doi.org/10.1186/s13071-016-1905-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-016-1905-7

Keywords