- Research
- Open access
- Published:
Impact of the intensity of infection in birds on Plasmodium development within Culex pipiens mosquitoes
Parasites & Vectors volume 18, Article number: 54 (2025)
Abstract
Background
In vector-borne diseases, invertebrate hosts are exposed to highly variable quantities of parasites during their blood meal. This heterogeneity may partly explain the overdispersed distribution of parasites within the vector population and the variability in the extrinsic incubation period (EIP) of the parasite. Indeed, the quantity of parasites ingested is often considered as a good predictor of the quantity of parasites that will develop within the vectors, as well as the speed at which they will develop (i.e. EIP). However, density-dependent processes can influence the relationship between parasite burden in the vertebrate host and in vectors, making this relationship unclear at times.
Methods
Here, we used an avian malaria system to investigate whether the proportion of red blood cells infected by sexual and/or asexual stages of Plasmodium relictum influences the intensity of infection and the EIP within vectors. For this purpose, we experimentally infected 12 birds in order to generate a range of infection intensity. More than 1000 mosquitoes took a blood meal on these hosts, and the development of Plasmodium within the vectors was followed for more than 20 days.
Results
Our study reveals a negative relationship between the intensity of infection in birds and the time until 10% of mosquitoes become infectious (EIP10). A period of only 4 days was sufficient to detect sporozoites in at least 10% of mosquitoes fed on the most infected hosts. However, the number of sporozoites did not vary significantly according to the vertebrate host intensity of infection, but was positively correlated to the oocyst burden (parasitic stage preceding the sporozoite stage).
Conclusions
While the quantity of ingested parasites had no impact on oocyst and sporozoite burden in infectious mosquitoes, the EIP10 was affected. Studies have demonstrated that small changes in the EIP can have a significant effect on the number of mosquitoes living long enough to transmit parasites. Here, we observed a difference of 4–6 days in the detection of the first sporozoites, depending on the intensity of infection of the bitten vertebrate host. Considering that a gonotrophic cycle lasts 3–4 days, the shortened EIP may have significant effects on Plasmodium transmission.
Graphical Abstract

Background
The rate at which a population grows or declines is influenced not only by environmental conditions but also by the number of individuals already present in the population. This phenomenon, studied in population ecology, is known as a density-dependent process [1,2,3,4]. Density-dependent mechanisms are considered of central importance in stabilizing the population growth of both free-living [5] and parasitic organisms [6, 7]. Parasites are especially subject to density-dependent processes, as they often have complex life cycles involving bottlenecks upon transmission followed by exponential growth [2, 7,8,9]. Both positive and negative density-dependent processes are common in host–parasite associations, and both may influence the parasite transmission rate [10, 11] and resilience of infection and control interventions [12, 13]. Parasites are typically overdispersed within host populations (i.e. are aggregated, and the variance is considerably greater than the average [14, 15]), and the strength of density-dependent processes therefore differs greatly among hosts. Understanding the causes of parasite overdispersion is an important step towards assessing the impact of density-dependent processes on parasite transmission dynamics [16, 17].
Parasite overdispersion has been widely described in several vector-borne diseases [18,19,20,21]. However, most studies have focused on the causes and consequences of aggregate parasite distribution in vertebrate host populations, leaving aside what happens within invertebrate host populations (i.e. the vector). Yet, for many vector-borne diseases, key traits determining the transmission dynamics of the parasite (e.g. lifespan, feeding behaviour of vectors) may depend on the intensity of parasite infection in the vector [22,23,24,25,26,27]. Fluctuations in parasite density within vector populations may be due to the genetic background or age of the invertebrate hosts, which influence their susceptibility to parasite infection [18, 28, 29]. Environmental factors may also be involved, such as temperature, which affects the parasite’s ability to complete its life cycle within the invertebrate host [30, 31]. The number of parasites ingested by the vector during their infestation (usually a blood meal) can also generate infection heterogeneity (e.g. [32, 33]). Indeed, parasite numbers vary between vertebrate hosts and also spatio-temporally within individuals [34, 35]. Although a positive correlation is expected between the quantity of parasites present in the vertebrate host blood and the density of transmissible stages of the parasite found in vectors (after the incubation period), the relationship is not always obvious. This is particularly true for Plasmodium parasites, the aetiological agents of one of the most studied vector-borne diseases, malaria. The relationship between the proportion of red blood cells of the vertebrate host infected by gametocytes (i.e. the transmissible stage) and the parasite burden in the mosquito seems to be specific to Plasmodium–mosquito associations. For the human malaria parasite Plasmodium falciparum and both Anopheles coluzzii and Anopheles gambiae mosquitoes, the parasite burden within the mosquito midgut (i.e. the oocyst stage) increases with the density of gametocytes [36,37,38,39]. On the other hand, this relationship is either negative or unclear in the case of rodent and avian malaria, respectively [37, 40].
To date, the majority of the studies trying to elucidate the relationship between vertebrate host gametocytaemia (i.e. percentage of red blood cells infected by a mature gametocyte) and the parasite burden in mosquitoes only consider the oocyst stage, which is the simplest and least expensive measure of the mosquito infection. Even though oocyst burden seems to be positively correlated with the likelihood of transmitting sporozoites to a new vertebrate host [41], the extrinsic incubation period (EIP) of the parasite should be considered. EIP refers to the time required for the parasite to replicate and disseminate to the mosquito salivary glands after ingestion of the infected blood meal. It determines the time at which the vector will be infectious and able to transmit the parasite to a new vertebrate host: the shorter the EIP, the greater the opportunity for the mosquito to transmit the parasite [42].
The aim of this study was to assess the role of the vertebrate host parasitaemia (i.e. the percentage of red blood cells infected by either sexual or asexual stages of Plasmodium) and gametocytaemia on the within-vector parasite infection dynamics. We investigated whether the parasitaemia and gametocytaemia in the vertebrate host determine the density of oocysts in the mosquito midgut and/or the EIP of Plasmodium. Specifically, we studied the time for 10%, 50% and 90% of the infected mosquitoes to become infectious (denoted by EIP10, EIP50 and EIP90, respectively). We predicted that mosquitoes biting highly infected hosts would have higher oocyst burden and be infectious earlier than mosquitoes feeding on hosts with low parasitaemia and gametocytaemia. To test this prediction, we used a natural avian malaria system [40, 43], with domestic canaries (Serinus canaria) as the vertebrate host, a haemosporidian parasite recently isolated from the wild (Plasmodium relictum), and a natural vector of this parasite originally collected from a natural population of mosquitoes and reared under laboratory conditions thereafter (Culex pipiens [44]). Avian malaria is a relevant model for studying the ecology and evolution of host–parasite interactions, both in the field and in the laboratory. Indeed, in addition to the fact that avian malaria is widespread in birds and mosquitoes throughout the world, this biological system is amenable to laboratory experiments at every stage in the parasite’s life cycle (see [40, 43] for review).
Methods
Biological material
The experiment was carried out in May 2021 using a P. relictum (lineage SGS1) strain isolated from an infected house sparrow (Passer domesticus) captured in December 2020 on the campus of the University of Lausanne, Switzerland (46°31′25.607″N, 6°34′40.714″E). This strain was maintained through six regular passages across our stock canaries (S. canaria) using intraperitoneal (i.p.) injections until the beginning of the experiment [40].
The Cx. pipiens mosquito population used in this experiment was initiated from wild clutches collected in Lausanne in August 2017 and maintained in the insectary thereafter. Mosquitoes were reared using standard protocols [45]. On the day prior to a mosquito blood meal, around 1000 7–10-day-old females were randomly captured from different emergence cages and placed inside experimental cages (90 females per cage, 40 × 40 × 40 cm). During this time, females had access to water but were deprived of sugar solution, to maximize the biting rate.
Experimental design
Infection of vertebrate hosts
Twelve domestic canaries (S. canaria) were infected with P. relictum (lineage SGS1) by i.p. from our stock canaries kept in the lab since January 2021, i.e. for 5 months prior to the experiment. One bird died on day 5 post-infection and was therefore removed from the analyses. The dynamics of infection were monitored for each bird every 2 days from day 5 to day 20 post-infection by measuring the parasitaemia and gametocytaemia by microscopic examination using Giemsa-stained blood smears. Blood samples (3–5 µl) were taken from the medial metatarsal vein. The number of red blood cells infected by sexual and/or asexual Plasmodium stages were counted per 3000–4000 erythrocytes in randomly chosen fields on the blood smears [46].
Mosquito exposure to P. relictum
On day 12 after vertebrate host infection, corresponding to the peak of infection for our parasite strain, birds were exposed individually to 90 uninfected mated female mosquitoes for 3 h (6–9 p.m.). At the end of the experiment, the number of blood-fed and unfed mosquitoes was counted, and only fully blood-fed females were kept individually in plastic tubes under standard laboratory conditions (25 °C, 70% relative humidity [RH]) with 10% glucose ad libitum. A total of 36 fully blood-fed mosquitoes per bird were used as samples for dissections (see Table S1 for details).
Mosquito dissections
Every 2 days starting from day 4 to day 20 post-blood meal, 4 ± 1 mosquitoes per bird were randomly sampled in order to monitor the dynamics of both oocyst formation and sporozoite production. Each mosquito was dissected to (i) count the number of oocysts in its midgut with the aid of a binocular microscope and (ii) quantify the sporozoites in its head/thorax using real-time quantitative polymerase chain reaction (qPCR; see molecular analyses below). The appearance of sporozoites in the head/thorax homogenate is strongly correlated with the appearance of sporozoites in the salivary glands [47]. The amount of haematin (a by-product of the decomposition of haemoglobin) deposited in the individual resting tubes was used to estimate the size of the blood meal. Quantification of the haematin was carried out by absorbance as described in previous papers [48]).
Molecular analyses
To study the presence/absence and the density of Plasmodium sporozoites in the mosquito head/thorax, qPCR was conducted. Beads were added to each head/thorax sample, after which the sample was flash-frozen in liquid nitrogen and then ground. Next, DNA was extracted using the Qiagen DNeasy Blood & Tissue Kit following the manufacturer’s instructions, except that the samples were incubated overnight. For each individual (i.e. mosquito), two qPCRs were carried out: one targeting the mitochondrial DNA (mtDNA) cytochrome b (cytb) gene of Plasmodium (primers L4050Plasmo 5′-GCTTTATGTATTGTATTTATAC-3′, H4121Plasmo 5′-GACTTAAAAGATTTGGATAG-3′, probe TexasRed-CYTB-BHQ2 5′-CCTTTAGGGTATGATACAGC-3′) and the other targeting the CQ11 gene of Cx. pipiens mosquitoes (primers 1725-F 5′-GCGGCCAAATATTGAGACTT-3′, 1726-R 5′-CGTCCTCAAACATCCAGACA-3′, probe FAM-CQ11-BHQ1 5′-GGAACATGTTGAGCTTCGGK-3′ [49]. All samples were run in triplicate (QuantStudio 6 and 7 Pro Real-Time PCR Systems), and samples with a cycle threshold (Ct) value higher than 35 for the parasite were considered uninfected. Relative quantification (RQ) values were calculated to assess the parasite prevalence and can be interpreted as the fold amount of the target gene (Plasmodium CYTB) with respect to the amount of the reference gene (Cx. pipiens CQ11), calculated as 2−(CtCYTBPlasmodium −CtCQ11Cx. pipiens).
Statistical analyses
Analyses were carried out using R statistical software (v. 4.2.1). The relationship between the birds’ gametocytaemia and parasitaemia was studied using a linear regression. We verified the model’s assumptions by plotting residuals using the simulateResiduals function from the DHARMa [Residual Diagnostics for Hierarchical (Multi-Level/Mixed) Regression Models] package [50]. This checking procedure was applied to all the models presented below. Because of the high degree of collinearity between parasitaemia and gametocytaemia, these two variables were never fitted as response variables in the same model. The influence of parasitaemia or gametocytaemia on the day post-blood meal at which the maximum burden of oocysts or sporozoites was reached, as well as the day at which 10%, 50% or 90% of the blood-fed mosquitoes in a batch were detected as sporozoite carriers (i.e. EIP10, EIP50 and EIP90) was analysed using linear models. The average blood meal size of each mosquito batch was also added as an explanatory variable in the models. The influence of parasitaemia or gametocytaemia on mosquitoes’ blood meal size, oocyst burden and sporozoite density was analysed by fitting bird as a random factor into the models using lmer or glmer.nb (package: lme4 [51]) according to whether the errors were normally (haematin quantity, sporozoite density) or negative-binomially distributed (oocyst burden). Blood meal size and oocyst burden were added as explanatory variables in the models when the statistical unit was a mosquito and when the response variable might be influenced by these traits (see Table S2). The average blood meal size of the batch of mosquitoes was added as an explanatory variable in all models where the statistical unit was the batch of mosquitoes and when the response variable might be influenced by these traits (see Table S2). Maximal models, including all higher-order interactions, were simplified by sequentially eliminating non-significant terms and interactions to establish a minimal model [52]. The significance of the explanatory variables was established using a likelihood ratio test which is approximately distributed as a Chi-square distribution or an F test [53]. The significant Chi-square or F values given in the text are for the minimal model, whereas non-significant values correspond to those obtained before the deletion of the variable from the model. The different statistical models built to analyse the data are described in Table S2 in the electronic supplementary material.
Results
Parasitaemia and gametocytaemia of vertebrate host
This experiment was carried out to investigate whether parasitaemia and gametocytaemia in the vertebrate host on the day of the infectious blood meal underpinned the parasite development within the invertebrate vectors. Eleven vertebrate hosts out of the 12 that were experimentally infected survived until day 12, allowing a range of parasitaemia and gametocytaemia values to be generated. The parasitaemia on the day of the mosquito blood meal (i.e. day 12 post-infection) followed an aggregated distribution (variance-to-mean ratio = 6.89, mean ± standard deviation [SD]: 3.47 ± 4.89). Three hosts had a parasitaemia value of less than 1%, four individuals had parasitaemia ranging between 1% and 2%, and four individuals had a parasitaemia value of over 3%, one of which presented a parasitaemia value of over 17% (Fig. 1A). The percentage of red blood cells infected with mature gametocytes was on average low (mean ± SD: 0.0015 ± 0.002) relative to what has been observed in other studies using the same Plasmodium strain. As reported in a previous study (e.g. [40, 54]), we observed a strong positive correlation between parasitaemia and gametocytaemia (model 1: F = 27.68, P < 0.001, adjusted R-squared: 0.724, Fig. 1A). As the conclusions of the subsequent analyses are the same whether parasitaemia or gametocytaemia was fitted as explanatory variable, we focus hereafter solely on the influence of parasitaemia. Parasitaemia measurement was indeed much less prone to slight error, since it only involved counting the number of infected red blood cells. In the case of gametocytaemia measurement, it was necessary to differentiate the parasite stages and their maturation levels, which at times can be somewhat subjective (at present, there are no molecular tools for identifying mature gametocytes). All analyses incorporating gametocytaemia instead of parasitaemia fitted as explanatory variables are presented in Appendix 1.
A Relationship between vertebrate host parasitaemia and gametocytaemia. Eleven individuals were included in this analysis. B Temporal dynamics of Plasmodium development in mosquitoes. Average oocyst burden (green, left axis) and sporozoite counts (salmon, right axis) at each dissection day. Green and salmon shadows represent standard error. The left axis represents the average number of oocysts counted per female. The right axis represents the amount of sporozoites quantified by qPCR
Temporal dynamics of Plasmodium development in mosquitoes
The temporal dynamics of oocyst production inside the mosquito midgut were consistent between mosquito batches fed on the 11 infected birds (Fig. 1B, Fig. S1). Oocysts were detected in all batches of mosquitoes from the first day of dissection, i.e. the fourth day after the blood meal. The peak oocyst burden was reached between days 6 and 8 (Fig. 1B, Fig. S1). Subsequently, a swift reduction in the oocyst count was noted, such that by day 14 post-blood meal, only a residual burden ranging from 1% to 20% of the peak oocysts persisted. The first sporozoites appear on the head/thorax homogenate mainly after day 6 (Fig. 1B). However, for two mosquito batches, sporozoites were detected in some females on day 4 post-blood meal (Fig. S1). The peak concentration of sporozoites was reached between days 10 and 12 post-blood meal for half of the batches of mosquitoes and between days 16 and 18 for the other half (Fig. 1B, Fig. S1).
Effect of vertebrate host parasitaemia on oocyst production
The size of the blood meal of mosquitoes was not influenced by the parasitaemia of the bird on which they fed (mean ± standard error [SE], 20.07 µg ± 0.055, model 2: χ2 = 0.256, P = 0.613). Since all mosquitoes already had oocysts in their midgut when we started the dissections (day 4 post-blood meal), we were unable to study the relationship between parasitaemia and the first appearance of oocysts. We showed that parasitaemia of vertebrate hosts had no influence on the day post-blood meal on which the oocyst peak was reached (model 3: χ2 = 0.456, P = 0.518, Fig. 2A) or on the density of oocysts reached during this peak (model 4: χ2 = 2.334, P = 0.127, Fig. 2B). While the blood meal size influenced the oocyst burden (the larger the blood meal, the greater the number of oocysts, model 4: χ2 = 12.65, P = 0.003), it did not influence the day post-blood meal on which the oocyst peak was reached (model 3: χ2 = 0.216, P = 0.655).
Relationship between vertebrate host parasitaemia and development of oocysts in the midgut of mosquitoes. A Parasitaemia did not significantly influence the number of days post-blood meal at which the maximum oocyst burden was reached. Each green dot corresponds to the batch of mosquitoes fed on a vertebrate host where the maximum oocyst burden was reached. B Parasitaemia of the bitten bird did not significantly influence the oocyst burden measured on the day when the oocyst burden was maximal in a focal mosquito batch. Each point represents the number of oocysts in the midgut of a dissected mosquito on the day when the mosquito batch associated with a bird, and therefore with a specific parasitaemia, reached maximum oocyst load
Effect of vertebrate host parasitaemia on sporozoite production
The time required for sporozoites to be detected in the head/thorax homogenate of 10% of infected mosquitoes was influenced by the parasitaemia of the vertebrate hosts on which the mosquitoes fed (EIP10, model 5: F = 9.82, P= 0.016, Fig. 3A). Batches of mosquitoes that had bitten the most infected vertebrate hosts showed sporozoites as early as day 4 post-blood meal, which was up to 4–6 days earlier than females that had bitten the least infected vertebrate hosts (Fig. 3A). Fitting the quadratic term (parasitaemia2) slightly improved the model fit (model 5: χ2 = 6.45, P = 0.039), suggesting that the day on which the first sporozoites were detected was a decelerating polynomial function of vertebrate host parasitaemia. The average blood meal size of each mosquito batch did not influence EIP10 (model 5: χ2 = 2.18, P = 0.183). The effect of bird parasitaemia on the number of days after the blood meal required for 50% and 90% of the mosquitoes to become infectious was however no longer significant (EIP50, model 6: χ2 = 2.65, P = 0.148, EIP90, model 7: χ2 = 0.81, P = 0.397). The average blood meal size of each mosquito batch also had no effect (EIP50, model 6: χ2 = 1.93, P = 0. = 207, EIP90, model 7: χ2 = 0.39, P = 0.551). As previously reported, peak sporozoite concentrations in mosquitoes were reached between day 10 and 18 post-blood meal, but the time to reach the peak was not influenced by vertebrate host parasitaemia (model 8: F = 0.581, P = 0.468) or by blood meal size (model 8: F = 1.31, P = 0.727). The maximum sporozoite density reached during the peak was also not influenced by the parasitaemia of the vertebrate hosts (model 9: χ2 = 0.05, P = 0.815) or by the blood meal size (model 9: χ2 = 2.67 P = 0.102) but was influenced by the oocyst load reached during the oocyst peak (model 9: χ2 = 7.575, P = 0.006; Fig. 3B), where the higher the oocyst burden, the higher the sporozoite load.
A Significant influence of parasitaemia on the day post-blood meal at which the EIP10 was reached (i.e. at least 10% of infected mosquitoes carried sporozoites). Each salmon dot corresponds to a batch of mosquitoes fed on a vertebrate host. B Positive relationship between average oocyst burden and average sporozoite load quantified in each batch of mosquitoes on the day of maximum oocyst and sporozoite load, respectively
Discussion
The main finding presented in this study reveals a negative relationship between the intensity of malaria infection in the vertebrate host and the delay required after biting for the first mosquito vectors to become infectious. We showed that the time required for 10% of infected mosquitoes to become infectious (EIP10) was reached by the batches of vectors fed on the most infected vertebrate hosts as early as day 4 after the blood meal. However, taking a blood meal on highly infected vertebrate hosts did not affect the day after the blood meal required for 50% and 90% of mosquitoes to become infectious or the intensity of parasite infection within the mosquito. The number of oocysts and the concentration of sporozoites were not significantly affected by the vertebrate host parasitaemia or gametocytaemia (see Appendix 1 for gametocytaemia).
In several natural systems, such as P. falciparum and An. coluzzii, P. falciparum and An. gambiae, Plasmodium vivax and Anopheles dirus, or Plasmodium mexicanum and Lutzomyia vexator, it has been found that the oocyst burden in the invertebrate vector is positively correlated to the gametocyte density [32, 37, 39, 55]. This relationship, however, was fairly weak in many studies, and there was considerable variability in the number of oocysts associated with a given parasitaemia or gametocytaemia value (e.g. [55,56,57,58]). Density-dependent processes are at play in the regulation of the parasite journey within the vector, with major bottlenecks occurring during the life-stage transition from gametocyte to oocyst [7, 25, 59, 60]. Here, we observe that the maximum oocyst burden reached during peak oocyst development was not significantly influenced by parasitaemia or gametocytaemia (see Appendix 1 for gametocytaemia). The number of vertebrate hosts included in the study and the widely dispersed distribution of gametocytaemia (the majority of birds had very low gametocytaemia levels, variance-to-mean ratio = 38) could explain why the relationship between these two parameters was not statistically significant. Similar conclusions were however reported in a previous study carried out with the same biological system on 40 birds with a much wider range of gametocytaemia (i.e. 0.1–9% [40]). Therefore, our results and those published previously point in the same direction and strongly suggest that in the natural P. relictum/Cx. pipiens system, oocyst burden is not strongly influenced by gametocyte density [33, 40]. This could be explained by the fact that not all gametocytes counted in blood smears were yet infectious (i.e. ready to produce micro- or macro-gametes in the mosquito’s midgut). The mechanisms underlying gametocyte infectivity remain poorly understood. Although we know that gametocytes go through several stages of development before reaching the mature stage [61], we do not know whether the mature stage is systematically infectious.
Parasitaemia and gametocytaemia in birds also had no influence on the maximum sporozoite density measured in the head/thorax of Cx. pipiens. However, we observed a positive relationship between maximum sporozoite density and the maximum average oocyst burden. Although the number of sporozoites produced per oocyst can vary from a few dozen to several thousand depending on the parasite/mosquito pair [62], a positive correlation between these two parasite stages is commonly reported [47, 55, 63,64,65]. The strength of the relationship between oocyst burden, sporozoite density and the rate of transmission of Plasmodium to vertebrate hosts is however difficult to assess due to the methodological complexity associated with the quantitative monitoring of the different parasite stages within a single mosquito. Nonetheless, a recent study conducted on both human and rodent malaria systems has clearly demonstrated a significant positive correlation between mosquito parasite load (oocyst and sporozoite stages) and the quantity of sporozoites injected during probing (i.e. inoculum size [65]).
In addition to inoculum size, the time lapse between mosquito infection (i.e. ingestion of gametocytes during the blood meal) and the moment when they become infectious (i.e. EIP) also has a major influence on the transmission dynamics of Plasmodium. Small changes in EIP can have a large effect on the number of mosquitoes living long enough to transmit parasites and subsequently on the number of hosts a vector might infect over its lifespan (e.g. [29, 66]). For the majority of Plasmodium species, EIP is on average between 8 and 16 days (see Table 1 in [67] and e.g. [68,69,70,71] for studies carried out on P. relictum). Several studies have shown that EIP can be highly heterogeneous between mosquitoes due to differences in nutrition (e.g. [66, 72,73,74]) and environmental conditions (e.g. temperature, see Fig. 1 in [42]); however, considerable variation remains even after accounting for these factors. Our results suggest that variation in EIP could be explained in part by the quantity of parasites ingested by the mosquitoes. Indeed, we detected the first sporozoites in batches of mosquitoes fed on the most infected hosts, on average 4–6 days earlier than in batches of mosquitoes fed on hosts characterized by lower parasitaemia and gametocytaemia. Interestingly, a previous study carried out on the same biological system, but with only three infected vertebrate hosts, highlighted a similar trend (see Fig. S1 in [75]). However, the effect of host parasitaemia on the evolution of the prevalence of infectious mosquitoes over time gradually diminished. We note that parasitaemia no longer had a significant effect on the EIP50 or EIP90. This result could be explained by the distribution of the quantity of parasites ingested by each mosquito within a batch of mosquitoes fed on the same host. Indeed, within a batch, the quantity of parasites ingested by each vector can vary greatly (e.g. Fig. 2 in [76]). The mosquitoes which ingested a greater quantity of parasites could correspond to those which we identified as being infectious as early as days 4–6 post-blood meal within the batches of mosquitoes fed on the most infected birds. To explore this question, a method would need to be developed to estimate the quantity of parasites ingested from the quantity of haematin excreted and the proportion of infected red blood cells in the host’s blood. To our knowledge, the only existing method of approximating the quantity of parasites ingested involves killing the mosquitoes before they digest the blood meal (e.g. [34]).
Several studies have suggested that the higher the parasite load in the mosquito midgut, the lower their survival rate (e.g. [23, 77]). In a situation of overcrowding, an adaptive strategy for the parasite to ensure transmission could therefore accelerate its development in order to colonize the mosquito’s salivary glands as quickly as possible, so that it can be transmitted to a new host before the vector dies. This pattern may be explained by density-dependent processes which may operate during the first steps of Plasmodium development within the mosquito midgut [9, 25, 59]. A high gametocyte density in the blood meal may lead to competition between parasites at different steps of early sporogony (e.g. fecundation, resources, colonization of the midgut wall). This competition could accelerate the development of all parasites and/or favour fast-growing parasite variants, both leading to early development of the first oocysts and, consequently, the rapid appearance of the first sporozoites. A theoretical study supports this hypothesis, but their results also suggest that above a certain density, the number of parasites in the mosquito midgut is so high that it could have a negative effect on the EIP (e.g. by negatively impacting oocyst survival rates or delaying oocyst burst, for example [27]).
Conclusions
In summary, our study has demonstrated for the first time that the intensity of infection within the vertebrate host on which a mosquito feeds influences the EIP of Plasmodium. The first sporozoites were indeed observed in the batches of mosquitoes fed on the most infected hosts. It should be noted, however, that in this study we only used acutely infected hosts, and therefore with high parasitaemia. Indeed, as this was a pioneering study, we wanted to generate the widest possible range of infection intensities to detect a signal that might have been weak. In the field, the vast majority of vertebrate hosts are in the chronic stage of infection [78,79,80], and it would then be relevant to determine whether the results we obtained are translatable to milder parasitaemia. It should also be mentioned that in this study, to control as many parameters as possible (e.g. genetic background of the mosquitoes, health of the vertebrate hosts, parasite strains), we used a colony of mosquitoes maintained in the laboratory since 2017, domestic canaries from a professional breeder and a Plasmodium strain isolated from the field approximately 8 months prior to the experiment. It will now be relevant to compare our results with those that could be obtained using a more natural system. Finally, our study showed that oocysts may be detected much earlier than the literature would suggest (e.g. [25, 75]). A future study on the subject should begin monitoring earlier after the blood meal to determine the influence of the vertebrate host parasitaemia on the appearance of the first oocysts.
Data availability
All data supporting the conclusions of this paper will be available on the Dryad Digital Repository.
Abbreviations
- EIP:
-
Extrinsic incubation period
- RH:
-
Relative humidity
References
Fowler CW. A review of density dependence in populations of large mammals. In: Genoways HH, editor. Current mammalogy. Boston: Springer, US; 1987. p. 401–41. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/978-1-4757-9909-5_10.
Stiling PD. The Frequency of density dependence in insect host-parasitoid systems. Ecology. 1987;68:844–56.
Chen L, Swenson NG, Ji N, Mi X, Ren H, Guo L, et al. Differential soil fungus accumulation and density dependence of trees in a subtropical forest. Science. 2019;366:124–8.
Albery GF, Newman C, Ross JB, MacDonald DW, Bansal S, Buesching C. Negative density-dependent parasitism in a group-living carnivore. Proc Royal Soc B Biol Sci. 2020;287:20202655.
May RM. Biological populations obeying difference equations: stable points, stable cycles, and chaos. J Theor Biol. 1975;51:511–24.
Anderson RM, May RM. Regulation and stability of host-parasite population interactions: I regulatory processes. J Anim Ecol. 1978;47:219.
Churcher TS, Ferguson NM, Basáñez M-G. Density dependence and overdispersion in the transmission of helminth parasites. Parasitology. 2005;131:121–32.
Hughes WO, Petersen KS, Ugelvig LV, Pedersen D, Thomsen L, Poulsen M, et al. Density-dependence and within-host competition in a semelparous parasite of leaf-cutting ants. BMC Evol Biol. 2004;4:45.
Churcher TS, Dawes EJ, Sinden RE, Christophides GK, Koella JC, Basáñez M-G. Population biology of malaria within the mosquito: density-dependent processes and potential implications for transmission-blocking interventions. Malar J. 2010;9:311.
Michael E, Snow LC, Bockarie MJ. Ecological meta-analysis of density-dependent processes in the transmission of lymphatic filariasis: survival of infected vectors. Me. 2009;46:873–80.
Pollitt Laura C, Churcher TS, Dawes EJ, Khan SM, Sajid M, Basáñez M-G, et al. Costs of crowding for the transmission of malaria parasites. Evol Appl. 2013;6:617–29.
Dietz K. Density-dependence in parasite transmission dynamics. Parasitol Today. 1988;4:91–7.
Hamley JID, Milton P, Walker M, Basáñez M-G. Modelling exposure heterogeneity and density dependence in onchocerciasis using a novel individual-based transmission model, EPIONCHO-IBM: Implications for elimination and data needs. PLoS Negl Trop Dis. 2019;13:e0007557.
Poulin R. Explaining variability in parasite aggregation levels among host samples. Parasitology. 2013;140:541–6.
Poulin R. Are there general laws in parasite ecology? Parasitology. 2007;134:763–76.
Luong LT, Vigliotti BA, Hudson PJ. Strong density-dependent competition and acquired immunity constrain parasite establishment: Implications for parasite aggregation. Int J Parasitol. 2011;41:505–11.
Vazquez-Prokopec GM, Perkins TA, Waller LA, Lloyd AL, Reiner RC, Scott TW, et al. Coupled heterogeneities and their impact on parasite transmission and control. Trends Parasitol. 2016;32:356–67.
Medley GF, Sinden RE, Fleck S, Billingsley PF, Tirawanchap N, Rodriguez MH. Heterogeneity in patterns of malarial oocyst infections in the mosquito vector. Parasitology. 1993;106:441–9.
Woolhouse MEJ, Dye C, Etard J-F, Smith T, Charlwood JD, Garnett GP, et al. Heterogeneities in the transmission of infectious agents: implications for the design of control programs. Proc Natl Acad Sci. 1997;94:338–42.
Cooper L, Kang SY, Bisanzio D, Maxwell K, Rodriguez-Barraquer I, Greenhouse B, et al. Pareto rules for malaria super-spreaders and super-spreading. Nat Commun. 2019;10:3939.
Martin LB, Addison B, Bean AGD, Buchanan KL, Crino OL, Eastwood JR, et al. Extreme competence: keystone hosts of infections. Trends Ecol Evol. 2019;34:303–14.
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. 2018;48:257–64.
Gutiérrez-López R, Martínez-de la Puente J, Gangoso L, Yan J, Soriguer R, Figuerola J. Experimental reduction of host Plasmodium infection load affects mosquito survival. Sci Rep. 2019;9:8782.
Basáñez MG, Townson H, Williams JR, Frontado H, Villamizar NJ, Anderson RM. Density-dependent processes in the transmission of human onchocerciasis: relationship between microfilarial intake and mortality of the simuliid vector. Parasitology. 1996;113:331–55.
Dawes EJ, Zhuang S, Sinden RE, Basáñez M-G. The temporal dynamics of Plasmodium density through the sporogonic cycle within Anopheles mosquitoes. Trans R Soc Trop Med Hyg. 2009;103:1197–8.
Courtenay O, Peters NC, Rogers ME, Bern C. Combining epidemiology with basic biology of sand flies, parasites, and hosts to inform leishmaniasis transmission dynamics and control. PLoS Pathog. 2017;13:e1006571.
Childs LM, Prosper OF. The impact of within-vector parasite development on the extrinsic incubation period. R Soc open sci. 2020;7:192173.
Pigeault R, Nicot A, Gandon S, Rivero A. Mosquito age and avian malaria infection. Malar J. 2015;14:383.
Somé BM, Guissou E, Da DF, Richard Q, Choisy M, Yameogo KB, et al. Mosquito ageing modulates the development, virulence and transmission potential of pathogens. Proc Royal Soc B Biol Sci. 2024;291:20232097.
Blanford JI, Blanford S, Crane RG, Mann ME, Paaijmans KP, Schreiber KV, et al. Implications of temperature variation for malaria parasite development across Africa. Sci Rep. 2013;3:1300.
Hlavacova J, Votypka J, Volf P. The effect of temperature on leishmania (Kinetoplastida: Trypanosomatidae) development in sand flies. J Med Entomol. 2013;50:955–8.
Schall JJ. Transmission success of the malaria parasite Plasmodium mexicanum into its vector: role of gametocyte density and sex ratio. Parasitology. 2000;121(6):575–80. http://www.journals.cambridge.org/abstract_S0031182000006818.
Enriquez GF, Bua J, Orozco MM, Macchiaverna NP, Otegui JAA, Argibay HD, et al. Over-dispersed Trypanosoma cruzi parasite load in sylvatic and domestic mammals and humans from Northeastern Argentina. Parasites Vectors. 2022;15:37.
Isaïa J, Rivero A, Glaizot O, Christe P, Pigeault R. Last-come, best served? Mosquito biting order and Plasmodium transmission. Proc Royal Soc B Biol Sci. 2020;287:20202615.
Pigeault R, Isaïa J, Yerbanga RS, Dabiré KR, Ouédraogo J-B, Cohuet A, et al. Different distribution of malaria parasite in left and right extremities of vertebrate hosts translates into differences in parasite transmission. Sci Rep. 2020;10:10183.
Bousema T, Drakeley C. Epidemiology and infectivity of Plasmodium falciparum and Plasmodium vivax gametocytes in relation to malaria control and elimination. Clin Microbiol Rev. 2011;24:377–410.
Da DF, Churcher TS, Yerbanga RS, Yaméogo B, Sangaré I, Ouedraogo JB, et al. Experimental study of the relationship between Plasmodium gametocyte density and infection success in mosquitoes; implications for the evaluation of malaria transmission-reducing interventions. Exp Parasitol. 2015;149:74–83.
Bradley J, Stone W, Da DF, Morlais I, Dicko A, Cohuet A, et al. Predicting the likelihood and intensity of mosquito infection from sex specific Plasmodium falciparum gametocyte density. eLife. 2018;7:e34463.
Guissou E, Da DF, de Hien DFS, Yameogo KB, Yerbanga SR, Ouédraogo GA, et al. Intervention reducing malaria parasite load in vector mosquitoes: No impact on Plasmodium falciparum extrinsic incubation period and the survival of Anopheles gambiae. PLOS Pathog. 2023;19:e1011084.
Pigeault R, Vézilier J, Cornet S, Zélé F, Nicot A, Perret P, et al. Avian malaria: a new lease of life for an old experimental model to study the evolutionary ecology of Plasmodium. Philos Trans Royal Soc B Biol Sci. 2015;370:20140300.
Miura K, Swihart BJ, Deng B, Zhou L, Pham TP, Diouf A, et al. Strong concordance between percent inhibition in oocyst and sporozoite intensities in a Plasmodium falciparum standard membrane-feeding assay. Parasit Vectors. 2019;12:206.
Ohm JR, Baldini F, Barreaux P, Lefevre T, Lynch PA, Suh E, et al. Rethinking the extrinsic incubation period of malaria parasites. Parasites Vectors. 2018;11(1):178. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13071-018-2761-4.
Rivero A, Gandon S. Evolutionary Ecology of Avian Malaria: Past to Present. Trends Parasitol. 2018;34:712–26.
Glaizot O, Fumagalli L, Iritano K, Lalubin F, Rooyen JV, Christe P. High prevalence and lineage diversity of Avian malaria in wild populations of Great Tits (Parus major) and mosquitoes (Culex pipiens). PLoS ONE. 2012;7:e34964.
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.
Valkiunas G. Avian Malaria Parasites and other Haemosporidia. CRC Press; 2004. https://www.taylorfrancis.com/books/9780203643792
Zollner GE, Ponsa N, Garman GW, Poudel S, Bell JA, Sattabongkot J, et al. Population dynamics of sporogony for Plasmodium vivax parasites from western Thailand developing within three species of colonized Anopheles mosquitoes. Malar J. 2006;5:68.
Rivero A, Ferguson HM. The energetic budget of Anopheles stephensi infected with Plasmodium chabaudi: is energy depletion a mechanism for virulence? Proc Royal Soc Lond Series B Biol Sci. 2003;270:1365–71.
Rudolf M, Czajka C, Börstler J, Melaun C, Jöst H, von Thien H, et al. First nationwide surveillance of Culex pipiens complex and Culex torrentium mosquitoes demonstrated the presence of Culex pipiens biotype pipiens/molestus hybrids in Germany. PLoS ONE. 2013;8:e71832.
Hartig F, Lohse L, Leite M de S. DHAR Ma: Residual diagnostics for hierarchical (multi-level/mixed) regression models. R Packag version 020. 2018. https://cran.r-project.org/web/packages/DHARMa/index.html.
Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.
Crawley MJ. The R book. Hoboken: John Wiley & Sons; 2012.
Bolker BM. Ecological models and data in R. Princeton: Princeton University Press; 2008.
Bantuchai S, Imad H, Nguitragool W. Plasmodium vivax gametocytes and transmission. Parasitol Int. 2022;87:102497.
Sattabongkot J, Maneechai N, Rosenberg R. Plasmodium vivax: gametocyte infectivity of naturally infected Thai adults. Parasitology. 1991;102:27–31.
Ponnudurai T, Lensen AHW, Gemert GJAV, Bensink MPE, Bolmer M, Meuwissen JHET. Infectivity of cultured Plasmodium falciparum gametocytes to mosquitoes. Parasitology. 1989;98:165–73.
Robert V, Read AF, Essong J, Tchuinkam T, Mulder B, Verhave J-P, et al. Effect of gametocyte sex ratio on infectivity of Plasmodium falciparum to Anopheles gambiae. Trans R Soc Trop Med Hyg. 1996;90:621–4.
Balabaskaran Nina P, Mohanty AK, Ballav S, Vernekar S, Bhinge S, D’souza M, et al. Dynamics of Plasmodium vivax sporogony in wild Anopheles stephensi in a malaria-endemic region of Western India. Malar J. 2017;16:284.
Sinden RE, Dawes EJ, Alavi Y, Waldock J, Finney O, Mendoza J, et al. Progression of Plasmodium berghei through Anopheles stephensi is density-dependent. PLoS Pathog. 2007;3:e195.
Poudel SS, Newman RA, Vaughan JA. Rodent Plasmodium population dynamics of early sporogony within Anopheles stephensi mosquitoes. Para. 2008;94:999–1008.
Gautret P, Motard A. Periodic infectivity of Plasmodium gametocytes to the vector. A review. Parasite. 1999;6:103–11.
Vaughan JA. Population dynamics of Plasmodium sporogony. Trends Parasitol. 2007;23:63–70.
Vaughan JA, Noden BH, Beier JC. Population dynamics of Plasmodium falciparum sporogony in laboratory-infected Anopheles gambiae. J Parasitol. 1992;78:716–24.
Gamage-Mendis AC, Rajakaruna J, Weerasinghe S, Mendis C, Carter R, Mendis KN. Infectivity of Plasmodium vivax and P. falciparum to Anopheles tessellatus; relationship between oocyst and sporozoite development. Trans Royal Soc Trop Med Hygiene. 1993;87:3–6.
Kanatani S, Stiffler D, Bousema T, Yenokyan G, Sinnis P. Revisiting the Plasmodium sporozoite inoculum and elucidating the efficiency with which malaria parasites progress through the mosquito. Nat Commun. 2024;15:748.
Shaw WR, Holmdahl IE, Itoe MA, Werling K, Marquette M, Paton DG, et al. Multiple blood feeding in mosquitoes shortens the Plasmodium falciparum incubation period and increases malaria transmission potential. PLoS Pathog. 2020;16:e1009131.
Stopard IJ, Churcher TS, Lambert B. Estimating the extrinsic incubation period of malaria using a mechanistic model of sporogony. PLoS Comput Biol. 2021;17:e1008658.
Žiegytė R, Bernotienė R, Bukauskaitė D, Palinauskas V, Iezhova T, Valkiūnas G. Complete sporogony of Plasmodium relictum (lineages pSGS1 and pGRW11) in mosquito Culex pipiens pipiens form molestus, with implications to avian malaria epidemiology. J Parasitol. 2014;100:878–82.
Kazlauskienė R, Bernotienė R, Palinauskas V, Iezhova TA, Valkiūnas G. Plasmodium relictum (lineages pSGS1 and pGRW11): complete synchronous sporogony in mosquitoes Culex pipiens pipiens. Exp Parasitol. 2013;133:454–61.
LaPointe DA, Goff ML, Atkinson CT. Thermal constraints to the sporogonic development and altitudinal distribution of avian Malaria Plasmodium relictum in Hawai’i. Para. 2010;96:318–24.
Palinauskas V, Žiegytė R, Iezhova TA, Ilgūnas M, Bernotienė R, Valkiūnas G. Description, molecular characterisation, diagnostics and life cycle of Plasmodium elongatum (lineage pERIRUB01), the virulent avian malaria parasite. Int J Parasitol. 2016;46:697–707.
Shapiro LLM, Murdock CC, Jacobs GR, Thomas RJ, Thomas MB. Larval food quantity affects the capacity of adult mosquitoes to transmit human malaria. Proc R Soc B. 2016;283:20160298.
Hien DFDS, Dabiré KR, Roche B, Diabaté A, Yerbanga RS, Cohuet A, et al. Plant-mediated effects on mosquito capacity to transmit human malaria. PLOS Pathogens. 2016;12:e1005773.
Pathak AK, Shiau JC, Freitas RCS, Kyle DE. Blood meals from ‘dead-end’ vertebrate hosts enhance transmission potential of malaria-infected mosquitoes. One Health. 2023;17:100582.
Sekar V, Rivero A, Pigeault R, Gandon S, Drews A, Ahren D, et al. Gene regulation of the avian malaria parasite Plasmodium relictum, during the different stages within the mosquito vector. Genomics. 2021;113:2327–37.
Talman AM, Ouologuem DTD, Love K, Howick VM, Mulamba C, Haidara A, et al. Uptake of Plasmodium falciparum gametocytes during mosquito bloodmeal by direct and membrane feeding. Front Microbiol. 2020;11:246. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2020.00246/full.
Dawes EJ, Churcher TS, Zhuang S, Sinden RE, Basáñez M-G. Anopheles mortality is both age- and Plasmodium-density dependent: implications for malaria transmission. Malar J. 2009;8:228.
van Rooyen J, Lalubin F, Glaizot O, Christe P. Altitudinal variation in haemosporidian parasite distribution in great tit populations. Parasites Vectors. 2013;6:139.
Loiseau C. Avian Plasmodium parasitaemia as an indicator of reproduction investment. PCIInfections. 2024;1:100194.
Asghar M, Westerdahl H, Zehtindjiev P, Ilieva M, Hasselquist D, Bensch S. Primary peak and chronic malaria infection levels are correlated in experimentally infected great reed warblers. Parasitology. 2012;139:1246–52.
Funding
Open access funding provided by University of Lausanne. This project was funded by the Swiss National Science Foundation (SNSF), grant 31003A_179378, to P.C.
Author information
Authors and Affiliations
Contributions
J.I., R.P., O.G. and P.C. conceived the study and elaborated the experimental design. J.I., M.B. and S.M. performed the experiments. M.B. and J.W. performed the molecular analyses. J.I. and R.P. analysed the data. J.I. wrote the first draft of the manuscript, and all authors contributed substantially to revision.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
This study was approved by the Ethical Committee of the Vaud Canton veterinary authority, authorization number 1730.4.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
13071_2024_6652_MOESM1_ESM.pdf
Additional file 1: Fig. S1 Temporal dynamics of Plasmodium development in batches of mosquitoes for each bird bitten/parasitaemia. Average oocyst burden (green, left axis) and sporozoite counts (salmon, right axis) in mosquitoes for each bird (parasitaemia) at each dissection day. Green and salmon shadows represent standard error. The left axis represents the average number of oocysts counted per female. The right axis represents the amount of sporozoites quantified by qPCR.
13071_2024_6652_MOESM2_ESM.pdf
Additional file 2: Table. S1 Blood meal rate per bird. BF = the number of fully blood-fed mosquitoes at the end of the experiment. NBF = the number of unfed mosquitoes at the end of the experiment. TOTAL = the number of uninfected mosquitoes placed in each cage at the beginning of the experiment. BM_RATE = the blood meal rate, which is the proportion of blood-fed mosquitoes for each bird.
13071_2024_6652_MOESM3_ESM.pdf
Additional file 3: Table. S2 Description of statistical models used in the study. N = sample size. “Maximal Model” includes the complete set of explanatory variables. “Minimal model” is the model containing only the significant variables and their interactions. Avg_Haematin corresponds to the average blood meal size of each mosquito batch. Square brackets indicate variables fitted as random factors. Curly brackets indicate the error structure used (n: normal errors, b: binomial errors). The response variable was not transformed unless otherwise stated.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Isaïa, J., Baur, M., Wassef, J. et al. Impact of the intensity of infection in birds on Plasmodium development within Culex pipiens mosquitoes. Parasites Vectors 18, 54 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13071-024-06652-4
Received:
Accepted:
Published:
DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13071-024-06652-4