Capacity of mosquitoes to transmit malaria depends on larval environment
Parasites & Vectors volume 7, Article number: 593 (2014)
Adult traits of holometabolous insects such as reproduction and survival can be shaped by conditions experienced during larval development. These “carry-over” effects influence not only individual life history and fitness, but can also impact interactions between insect hosts and parasites. Despite this, the implications of larval conditions for the transmission of human, wildlife and plant diseases that are vectored by insects remain poorly understood.
We used Anopheles stephensi mosquitoes and the rodent malaria, Plasmodium yoelii yoelii, to investigate whether quality of larval habitat influenced vectorial capacity of adult mosquitoes. Larvae were reared under two dietary conditions; one group received a diet commonly used for colony maintenance (0.3 mg/individual/day of Tetrafin fish food) while the other group received a reduced food diet (0.1 mg/individual/day). Upon emergence, adults were provided an infectious blood feed. We assessed the effects of diet on a range of larval and adult traits including larval development times and survival, number of emerging adults, adult body size and survival, gonotrophic cycle length, and mating success. We also estimated the effects of larval diet on parasite infection rates and growth kinetics within the adult mosquitoes.
Larval dietary regime affected larval survival and development, as well as size, reproductive success and survival of adult mosquitoes. Larval diet also affected the intensity of initial Plasmodium infection (oocyst stage) and parasite replication, but without differences in overall infection prevalence at either the oocyst or sporozoite stage.
Together, the combined effects led to a relative reduction in vectorial capacity (a measure of the transmission potential of a mosquito population) in the low food treatment of 70%. This study highlights the need to consider environmental variation at the larval stages to better understand transmission dynamics and control of vector-borne diseases.
Malaria is the most important vector-borne disease of humans worldwide, with approximately 219 million people infected annually, resulting in about 600,000 deaths per year . The transmission intensity of malaria is inextricably linked to the biology of the mosquito vectors and can be characterized using a summary metric known as the vectorial capacity (C). Vectorial capacity describes the rate at which future infections arise from a currently infected host (provided that all female mosquitoes become infected) and provides a measure of the transmission potential of a vector population ,. It is defined as:
where m is vector density (ratio of adult mosquitoes to humans), a is the daily probability of a human host being fed on by a vector, n is the extrinsic incubation period of the parasite, p is the daily probability of adult vector survival, and b is the proportion of mosquitoes with sporozoites disseminated in their salivary glands.
Any variation in environment that affects relevant aspects of vector biology could result in a change in transmission risk via effects on vectorial capacity -. Recent work shows that changes in temperature (both means and diurnal fluctuation) and rainfall events can have substantial effects on the transmission potential of malaria ,-. Other sources of environmental heterogeneity include differences in food resource availability, seasonality of habitats, and land use changes ,.
To date, many studies examining the effect of environment on mosquito biology and aspects of vectorial capacity have focused directly on the adult mosquitoes. This is logical, as it is only the adult female mosquitoes that transmit malaria and the frontline interventions used for control (such as insecticide treated bed nets, indoor insecticide sprays, screening, repellents etc.) primarily target the adult stage. Small changes in daily survival probability (p) and human biting rate (a), for example, can have very large effects on vectorial capacity, which explains in-part the effectiveness of insecticide treated nets as these act on both traits simultaneously -. Biting rate (determined by the duration of the gonotrophic cycle) has also been shown to be a major factor explaining variation in malaria incidence -. However, larval condition has become increasingly recognized as having an influence on adult mosquito life history traits -,-. Based on studies in other invertebrate systems, it is expected that variation in quality of larval habitats could feed through to impact adult life history, which in turn could affect transmission -. For terrestrial insects and other small invertebrates, estimates of body condition are often positively correlated with body size ,, and larger individuals often exhibit increased probability of survival, fecundity and ultimate overall fitness ,. Adult survival and vector density are key elements of vectorial capacity. Larval effects on adult body size might also be important in determining vector competence. Larger individuals could support more parasites due to greater availability of host resources ,. Alternatively, individuals with more reserves might be able to devote more energy toward immune defense ,.
Extensive studies on Aedes species have shown that larval environment can have considerable effects on life history traits important to transmission, such as development rate, adult longevity, and efficiency of egg development -. Larval environment has also been demonstrated to significantly shape vector competence in a variety of Aedes-arbovirus systems. For example, interactions between larval competition and density in Aedes aegypti and Aedes albopictus can significantly increase susceptibility to dengue virus and Sindbis virus ,. Nutritional stress in Aedes aegypti has also been shown to influence the interaction between humoral and cellular branches of the immune system, which could affect vector competence for a suite of pathogens .
A more limited number of studies on Anopheles vector spp. support the influence of larval condition on subsequent adult traits and vectorial capacity, but the patterns are not well understood ,-. Here, we use the Asian malaria vector, Anopheles stephensi (Liston), and a rodent model malaria, Plasmodium yoelii yoelii, to investigate whether nutritional quality of larval habitat affects vectorial capacity. We show trans-stadial impacts on a range of traits indicating potential for strong effects of larval rearing condition on subsequent transmission of malaria.
General experimental design
To manipulate larval environment quality, we fed An. stephensi larvae differing quantities of food throughout development. Larvae were collected from our standard lab colony at The Pennsylvania State University (this colony was initiated in 2011 with eggs from a longstanding colony maintained at Johns Hopkins University). Newly hatched (<24 h old) first instar larvae were transferred to 265 mL plastic cups containing 80 mL of distilled water at initial densities of 50 larvae per cup. Larvae were maintained on Tetrafin® fish food under standard insectary conditions (26°C ± 0.5°C, 80% humidity, and a 12:12 L:D photoperiod). We exposed larvae to one of two experimental food treatments: 1) a “high” food regime (0.3 mg/individual/day), which is consistent with our standard colony maintenance, and 2) a “low” food regime (0.1 mg/individual/day). These food treatments were selected based on a series of pilot studies. We maintained a constant food environment in each cup by counting larvae and pupae daily, replacing water (kept at a constant volume), and adjusting the amount of food allocated per individual. Each treatment group contained 48 replicate cups. The entire experiment was repeated in time, for a total of two experimental blocks (96 cups per treatment, in all). For mating success and gonotrophic cycle experiments (details below), each larval treatment group contained 15 cups in two experimental blocks (also repeated in time), totaling 30 cups per treatment.
Quantifying the effects of food treatment on mosquito life history traits
Larval survival and development time
To assess effects of experimental food treatment on larval survival, larvae and pupae in each cup were counted daily using an eyedropper. To estimate duration of larval development, and quantify adult mosquito production of each cup, the date of adult female emergence and number of females emerging were recorded.
Adult body size
To determine if larval food manipulation generated differences in adult female body size, we took one wing per mosquito (n = 100 per treatment, per experimental block) and mounted wings onto a glass microscope slide using clear nail lacquer. Wings were then measured using a micrometer eyepiece with a standard dissecting microscope. Measurements were taken from the tip of the wing (excluding fringe) to the distal end of the alula. Wing length is known to be positively correlated with body size in mosquitoes ,.
Gonotrophic cycle, mating success and daily adult survival
To determine effect of larval food treatment on gonotrophic cycle length, females that had emerged 3-5 days prior (and had been housed in mesh cages with males from their respective treatments) were fed to repletion on a Hemotek membrane feeder using pork sausage casing filled with human blood heated to 37°C. Blood-fed females (n = 50 per treatment in each experimental block, for a total of 100 per treatment) were then transferred to individual 50 mL plastic tubes covered with mesh, and provided with a cotton ball moistened with 10% glucose solution that was replenished daily. Each day, sugar moistened cotton balls were removed 4 hours prior to offering a blood meal (also via a Hemotek feeder) for 10 minutes each on four subsequent days. Tubes were monitored daily for eggs, which were counted for each individual female. Mortality was recorded until day 14 after the first blood meal.
Here, gonotrophic cycle is defined as the time (in days) from initial blood meal to laying of the first clutch of eggs, as few individuals laid twice throughout the experiment. Mating success was also assayed by dissection of spermathecae from females in both treatments that never laid eggs throughout the monitoring period. Presence of sperm, whether alive or dead, was considered a successful mating.
Adult survival was monitored for 14 days after the first blood meal, at which point adult females were 17-19 days old. While mosquitoes can live for many weeks under ideal lab conditions , this prolonged survival is difficult to interpret as it is generally accepted that few mosquitoes live beyond 2 weeks in field settings -. Moreover, 14 days after the blood meal is when female mosquitoes are potentially able to transmit malaria and before strong effects of senescence are expected (note the basic vectorial capacity equation is age-independent and assumes a constant rate of mortality in the absence of senescence).
Quantifying the effects of food treatment on vector competence
In order to assess the effects of larval food treatment on measures of vector competence, we first randomly allocated three-to-five day old individuals from each cup of both treatments across one of four replicate cages (80-100 females per cage). Females were allowed to feed on two anesthetized mice per cage (C57 females, > 6 weeks old) infected with Plasmodium yoelii yoelii (clone 17XNL) for thirty minutes. All vertebrate animal work was carried out by trained research technologists under Penn State University IACUC protocols specified in permit #27452. To ensure that mosquitoes received an infectious blood meal at the same adult age, first instar larvae for the low food treatment group were collected from the colony four days earlier than first instars to be allocated to the optimal food treatment group. This was done to adjust for the slower developmental time in the low food treatment group and allow for age-matched comparison of the groups at the same time after administration of the infectious blood meal. Individuals that did not feed were removed from the cage. Infectious feeds were performed at 26°C for optimum host seeking and probing behavior. Immediately following feeding, cages were transferred to a second incubator set at 24°C, as this is the thermal optimum for Plasmodium yoelii yoelii growth and replication. After infection, individuals were provided with cotton balls moistened with 10% glucose offered ad libitum and replaced daily. Daily survival of infected adults was monitored by counting and removing the dead individuals throughout the duration of the experiment.
Despite our high level of replication at the larval food level (48 cups per treatment per block), we chose to pool individuals from cups into cages (4 per treatment per block). Combining in this way means we cannot trace individual mosquitoes and their resulting infection dynamics back to a specific cup. However, we considered randomization across 4 treatment cages, with use of 2 mice per cage to account for natural mouse-to-mouse variation in gametocytemia, and complete replication through time, sufficient to detect any biologically relevant treatment effects. Putting individual mice on each cup would have increased the number of mice used to unethical levels.
Parasite prevalence and intensity
A subset of 25 individuals was taken from each replicate cage at seven days post-infection to assess infection prevalence and intensity. Midguts were dissected in 1X phosphate-buffered saline and examined under a light microscope. Presence or absence of oocysts and the number of established oocysts were quantified. Dissected midguts were saved individually in 1.5 mL tubes and stored in absolute ethanol at -80°C for future genomic DNA analysis (discussed below). At day 16-18 post-infection, 15 individuals were sampled from each replicate cage, and salivary glands were dissected out in 1X phosphate-buffered saline. Glands were examined under a light microscope and scored for presence or absence of sporozoites. To estimate vector competence (b, proportion of infectious bites on a susceptible host that lead to an infected host), we used the proportion of mosquitoes with sporozoites disseminated in the salivary glands. This is a standard approximation and assumes that if a mosquito has sporozoites in the salivary glands, it will likely transmit during feeding ,,.
To quantify how food treatment affected sporozoite production for each infected mosquito, we performed genomic DNA extraction and qPCR analysis for Plasmodium genomes in midguts saved from oocyst dissection. Plasmodium genomic DNA was extracted from midguts using the E.Z.N.A. MicroElute Genomic DNA kit (Omega Bio-Tek, as per the manufacturer’s protocol). DNA was eluted in 20 μL of molecular grade water, and the number of parasite genomes present in midguts was quantified using a previously developed qPCR assay . Briefly, reactions were run on an ABI Prism 7500 Sequence Detection System (TaqMan). Initial denaturation was 20 seconds at 95°C followed by 40 cycles of a three-second 95°C denaturation period and a 30-second 60°C period of annealing and extension. Primers and probes were designed to amplify DNA from several Plasmodium species. We constructed standard curves for P. yoelii genome detection by extracting DNA from a known number of infected mouse red blood cells using the BloodPrep kit (Applied Biosystems) on the ABI Prism 6100 Nucleic Acid Prep Station (as per the manufacturer’s protocol). Parasite production per oocyst was evaluated by dividing the total number of parasite genomes by the number of oocysts quantified for each midgut. We used both sporozoite production per midgut and per oocyst as measures of the efficiency of parasite replication.
All statistical analyses were carried out using IBM SPSS v.21 (Armonk, NY). For all analyses, full factorial models were reduced through backwards elimination of non-significant, higher order interactions, and henceforth non-significant higher order interactions are not reported in our discussion of the results or displayed in our model tables. All models were evaluated for goodness of fit by assessing model deviance per degrees of freedom, log likelihood values and residual plots. To assess significant pairwise comparisons, we used Bonferroni-adjusted post-hoc tests.
We used univariate general linear model (GLM) analysis to determine how food treatment affected adult body size and mating success. Larval and adult survival, day of emergence, number of adults emerged, fecundity and measures of vector competence (oocyst intensity and prevalence, sporozoite prevalence) were analyzed using generalized linear models (GZLM) so that non-normal error distributions could be used in the analysis. All distributions were chosen based on both best model fit and plots of raw data.
Normal distribution was assumed for analysis of day of adult female emergence and gonotrophic cycle. We assumed complimentary log-log, poisson, binomial, negative binomial, and gamma distributions in the analyses for mosquito survival, number of adults emerging and fecundity, oocyst and sporozoite prevalence, oocyst intensity, and number of sporozoites produced for each treatment group, respectively. For each dependent variable in our analyses, food treatment, cage replicate (for parasite traits), and experimental block were included as factors. We included oocyst intensity as a covariate in the model assessing treatment effects on sporozoite production. Across all models concerning characteristics of Plasmodium infection, replicate cage was nested within treatment to correct for the fact that mosquitoes in each cage received a blood meal from a different group of mice, and so were not related to one another across treatment groups.
Quantifying effects of larval food treatment in the context of the vectorial capacity equation
We calculated vectorial capacity using mean parameter estimates quantified from our empirical data (Table 1). In the current study we have no direct measure of the extrinsic incubation period (EIP) of the parasite (n), so we assumed the EIP for P. yoelii development at 24°C to be 12 days for both treatment groups based on previous research . For daily survival rates (p), we used the average rate over the entire 18-day monitoring period. Vector density (m) was estimated using the mean of total emerged females per replicate larval cup. We followed convention in using the reciprocal of the mean gonotrophic cycle length as a proxy for daily biting rate (a). The proportion of mosquitoes potentially infectious (b), was calculated using raw data means for sporozoite prevalence.
All values were calculated using means from our empirical data, except for extrinsic incubation period (EIP), which is assumed based on ideal conditions in previous work. Standard errors of means are displayed for vector density, biting rate, and proportion infectious (sporozoite prevalence).
The effects of food treatment on mosquito life history traits
Food treatment significantly affected larval survival (and henceforth the number of emerging adult females). Larvae reared on low food diets did not survive as well as those reared on high food diets (Figure 1a). However, this effect was strongly mediated by day (analyzed here as a covariate) and interactions among food treatment and experimental block and experimental block and day (Table 2, GZLM, interval-censored survival analysis). Although the trend of reduced larval survival and subsequent adult emergence in the low food group remained consistent across experimental blocks, the first block showed an initial (days 2-5 after hatching) sharp decline in the survival of the high food group, which did not replicate in the second experimental block resulting in the significant treatment x day x block interaction (Table 2). We did see significant replicate effects; however, upon plotting Kaplan-Meier survival curves for each of the 48 replicates per treatment per block, these effects were due to several cups that exhibited unusually high or low larval mortality (data not shown). However, the trends remain the same and when assessing the final productivity of each cup (total number of females), replicate is not a significant effect.
Low food environments produced, on average, significantly fewer total adult females per replicate cup than the high food environments (Figure 1c). In the low diet, mean emergence was 16.7 females (S.E. = 0.44), while high diet cups produced a mean of 19.7 (S.E. = 0.359) (p = <0.001, GZLM, Table 2). Adult females emerging from the high food treatment group had significantly larger wings than those emerging from the low food treatment group, with little overlap between the two groups (Figure 1b).
Length of gonotrophic cycle was significantly shaped by larval food treatment (Table 2), and females from low larval food diets had gonotrophic cycles that were on average one day longer than mosquitoes fed high food diets (Figure 2a). There was a significant effect of experimental block, with individuals in the second taking slightly longer to lay the first clutch of eggs, but this effect was consistent across treatments, resulting in no significant treatment x block interactions. Larval food treatment also significantly affected mating success (presence of sperm in spermatheca; Figure 2b) and fecundity. Females from the low food treatment were less likely to be mated than those from the high food treatment (mean proportion mated = 0.71 ± .45, 0.97 ± 0.18, respectively), and those in the low food treatment that did lay eggs, laid significantly fewer than high food females (mean size of first clutch = 28.67 ± 2.011, 64.14 ± 1.797, respectively; treatment p <0.001, GZLM, Poisson distribution).
Effects of food treatment on vector competence
Larval food treatment significantly affected adult infected female survival (p = 0.001, GZLM, interval-censored survival, Table 2). The difference was not significant until day five post-infection, when survival in the low food group fell sharply, then leveled off at day 10, with little mortality in either group for the remainder of the observation period (Figure 1c).
Despite differences in adult survivorship, food treatment had no significant effect on either oocyst or sporozoite prevalence (Table 3, Figure 3a). In contrast, food treatment significantly influenced oocyst intensity (Table 3, Figure 3b), and females from the low food group had an overall lower oocyst burden than those from the high larval food group. This effect was consistent across experimental blocks, despite the difference in infection levels between experiments.
Individuals fed a high food diet exhibited a greater number of parasite genomes per oocyst than mosquitoes fed a low food diet (mean of 1.69 × 105 and 2.2 × 105 genomes per midgut for low and high treatment, respectively) (Table 3). Unsurprisingly, the number of Plasmodium genomes per midgut increased with oocyst intensity. However, the slope of the positive relationship between oocyst intensity and number of Plasmodium genomes is steeper for individuals in the high food treatment (slope = 9.38E3, R2 = 0.787) than for the low food treatment (slope = 6.44E3, R2 = 0.747), as indicated by the significant interaction between food treatment and oocyst intensity (treatment x oocyst intensity, Table 3).
Here, the significance of experimental block is likely due to the fact that oocyst intensity in the second block (mean = 13.26) was lower than that of the first block (mean = 46.36). Though the trends replicate and remain significant, the magnitude of the relationship between intensity and genome count is much greater in the first block.
Effects of food treatment on vectorial capacity
For daily survival rates (p), we took the average rate over the entire 18-day monitoring period, which resulted in p = 0.973 and p = 0.983 for the low and high food groups, respectively. Mean adult female density (m) per replicate cup in the low diet was 16.7 females (S.E. = 0.44), while high diet cups produced a mean of 19.7 (S.E. = 0.359) (p = <0.001, GZLM, Table 2). Gonotrophic cycle for low food individuals was 4.29 days (the reciprocal of which gives an average biting rate per day of a = 0.233, S.E. = .009), while high food individuals had an average cycle length of 3.41 days (a = 0.293, S.E. = 0.006). Sporozoite prevalence (b) mean was 0.30 for both groups (S.E. = 0.042). These parameter values yielded mean vectorial capacities of C = 7.155 for mosquitoes from the low food environment and C = 22.462 for those from the high food environment.
In this study we demonstrate effects of variation in food abundance at the larval stage on both larval survival and a suite of adult mosquito life history traits that combine to determine vectorial capacity. The treatment effects equate to a c.70% relative reduction in malaria transmission potential for mosquitoes from the low food environment, or conversely, a relative increase in transmission potential of c.310% for mosquitoes from the high food environment.
The most obvious influence of larval food treatment was the impact on daily survival and prolonged development of larvae in poor nutritive environments. Longer development times and smaller proportion of mosquitoes successfully pupating and eclosing, has been reported in Anopheles gambiae, Aedes triseriatus and Anopheles darlingi,,. Further, although larval development is not a direct component of vectorial capacity, it will affect vector population growth rates and possibly vector density (depending on the nature of the density dependence ). Adult survival has also been shown to be affected by quality of larval habitat in other mosquitoes ,. One recent study suggested that quality of larval habitat had no effect on subsequent survival of adult An. stephensi and hence, no implications for transmission . However, this study monitored adult survival for 5 days only, which is much shorter than the incubation period of the parasite and so provides little insight into ultimate impacts on the number of adults potentially able to transmit the parasite (note that differences in survival between treatments in our study only emerged between days 5-9).
Larval food treatment also carried over to have a significant effect on adult body size. Several studies have shown positive correlation between body size and fitness in vector species. For example, in Aedes albopictus reared in both the laboratory and the field, pupal mass and wing length were consistently correlated positively with fecundity . In field-collected Aedes communis, wing length was a significant predictor of longevity when the adult females were food stressed . Studies on Anopheles gambiae show that reduction in resource availability through increasing larval density significantly decreases adult mass in males and females. Reduced body mass was demonstrated to be negatively correlated with both age at pupation and mating success, suggesting that both development time and body size have an influence on adult fitness -.
Larval food treatment also impacted duration of the gonotophic cycle, and hence the reciprocal estimate of daily biting rate. Changes in gonotrophic cycle length were relatively small (<1 day) but because biting rate enters into the vectorial capacity equation as a squared term (reflecting that one bite is required for a mosquito to acquire the parasite and at least one other required to pass it on), even small changes can contribute to differences in overall vectorial capacity ,,. The prolonged gonotrophic cycle length we observe in small females could be due to decreased nutritional reserves (as could the effects on fecundity). This reduction in reserves has been shown in Aedes aegypti females to prolong especially the first gonotrophic cycle, as more than 50% of the lipid resources needed to develop eggs are from larval stores . In nature, the timing of the initial blood meal could also be delayed if a smaller proportion of females engage in blood feeding, as seen in Anopheles darlingi emerging from poor larval habitats , or if smaller females have higher initial preference for sugar feeding as opposed to blood feeding . These latter effects would potentially further impact vectorial capacity.
Differences in fecundity and mating success, while they do not contribute directly to vectorial capacity, could further shape vector density. It has been demonstrated in several studies that fecundity is positively correlated with body size in several anopheline species, due to both available teneral reserves for ovarian development as well as larger blood meal size taken by larger females -. Possible explanations for impacts on mating success include depressed mate-seeking in small individuals of both sexes due to energetic demand ,. Anopheles gambiae males have been demonstrated to preferentially select larger females for mates , so low food females may have not presented as attractive mates for their conspecific males. Additionally, differences between treatments could affect the viability and number of sperm . In our studies, mosquitoes were restricted to mating within their respective treatment groups. These trends in mating success might change if mosquitoes were allowed to mate freely with individuals of differing sizes from diverse larval habitats .
The effects of larval diet on parasite development and overall vector competence were mixed. We observed differences in parasite intensity and genomic DNA replication, suggesting that larger mosquitoes might offer greater nutritional resources for parasite exploitation within the insect, but found no effect on parasite prevalence at either the oocyst or sporozoite stage. We used a rodent malaria as an established model  but there are recognized differences between typical infection intensities of rodent and human malarias , and the resultant mosquito immune responses ,. We also acknowledge that we are using a laboratory strain of An. stephensi that has known susceptibility to certain Plasmodium species and that natural sympatric vector-parasite pairings can differ in baseline competence . The functional importance of these differences is slightly unclear, but they complicate translation of results across systems and we do not suggest that our calculations of vectorial capacity can be applied directly to the field. That said, human malaria parasites do require mosquito resources for development - so there is no reason to think that host condition is unimportant. Several previous studies have demonstrated effects of larval density on adult traits including body size ,,, survival , blood intake, and mating competitiveness , suggesting our results to be robust. Comparative differences might be enhanced or diminished in the face of the complexity that exists in the field.
The effects we observe on parasite replication rate suggest the additional intriguing possibility that host condition could affect the extrinsic incubation period (EIP) of the parasite. We did not measure EIP in our study as we dissected salivary glands at a single time point during infection. However, differences in parasite replication rate due to the influence of temperature, for example, are known to affect length of the EIP, with important consequences for vectorial capacity , We are not aware of any studies exploring the effects of mosquito condition (larval through to adult) on malaria parasite development but nutritionally-dependent parasite replication has been observed in other invertebrate-parasite systems ,-. If EIP is affected by mosquito condition, this could add to the influence of larval habitat on malaria transmission.
Russell et al. (47) analysed a time series of An. gambiae abundance and condition (body size) data collected in a high transmission setting in central Tanzania. Their study revealed marked variation in adult body size across the year linked to rainfall-driven changes in availability (and likely quality) of larval breeding habitats. The density dependent feedbacks between mosquito condition and carrying capacity of larval habitats provided an elegant explanation for the observed adult mosquito population dynamics. Our data indicate that variation in quality of larval habitats could have important impacts on vectorial capacity beyond predicted effects on mosquito abundance.
Caveats regarding the laboratory nature of our experimental system notwithstanding, the magnitude of the effects we observe are non-trivial. A 70% reduction in vectorial capacity is of the same order as reported from well implemented control programs using insecticide sprays or insecticide treated nets -. Thus, variation in larval habitat quality could have a marked influence on the temporal and spatial dynamics of malaria transmission. Larval resource effects could also have important implications for control. For example, poorly implemented larval control, such as inefficient larviciding, could possibly alleviate density effects, increasing resource availability for the remaining larvae. This reduced larval competition might lead to inadvertent increases in transmission potential. Conversely, reduced nutrition via increases in larval competition could reduce transmission potential and might provide a partial explanation for the so-called ‘paddy paradox’, whereby malaria burden remains low in rice irrigation areas in spite of very high mosquito densities . The potential for such effects reinforces the need for better understanding the effects of environmental quality on larval and adult mosquito ecology .
World Malaria Report 2012. World Health Organization.
Garrett-Jones C:Prognosis for interruption of malaria transmission through assessment of the mosquito’s vectorial capacity. Nature. 1964, 204: 1173-1175. 10.1038/2041173a0.
Smith DL, McKenzie FE: Statics and dynamics of malaria infection in Anopheles mosquitoes. Malar J. 2004, 3: 13-10.1186/1475-2875-3-13.
Paaijmans KP, Blanford S, Chan BHK, Thomas MB: Warmer temperatures reduce the vectorial capacity of malaria mosquitoes. Biol Lett. 2012, 8: 465-468. 10.1098/rsbl.2011.1075.
Mordecai EA, Paaijmans KP, Johnson LR, Balzer C, Ben-Horin T, de Moor E, McNally A, Pawar S, Ryan SJ, Smith TC, Lafferty KD: Optimal temperature for malaria transmission is dramatically lower than previously predicted. Ecol Lett. 2013, 16: 22-30. 10.1111/ele.12015.
Afrane YA, Little TJ, Lawson BW, Githeko AK, Yan G:Deforestation and vectorial capacity of Anopheles gambiae Giles mosquitoes in malaria transmission, Kenya. Emerg Infect Dis. 2008, 14: 1533-1538. 10.3201/eid1410.070781.
Afrane YA, Githeko AK, Yan G:The ecology of Anopheles mosquitoes under climate change: case studies from the effects of deforestation in East African highlands. Ann N Y Acad Sci. 2012, 1249: 204-210. 10.1111/j.1749-6632.2011.06432.x.
Fontenille D, Simard F: Unravelling complexities in human malaria transmission dynamics in Africa through a comprehensive knowledge of vector populations. Comp Immunol Microbiol Infect Dis. 2004, 27: 357-375. 10.1016/j.cimid.2004.03.005.
Paaijmans KP, Read AF, Thomas MB: Understanding the link between malaria risk and climate. Proc Natl Acad Sci U S A. 2009, 106: 13844-13849. 10.1073/pnas.0903423106.
Paaijmans KP, Blanford S, Bell AS, Blanford JI, Read AF, Thomas MB: Influence of climate on malaria transmission depends on daily temperature variation. Proc Natl Acad Sci U S A. 2010, 107: 15135-15139. 10.1073/pnas.1006422107.
Blanford JI, Blanford S, Crane RG, Mann ME, Paaijmans KP, Schreiber KV, Thomas MB: Implications of temperature variation for malaria parasite development across Africa. Sci Rep. 2013, 3: 1300-10.1038/srep01300.
Afrane YA, Lawson BW, Githeko AK:Effects of microclimatic changes caused by land use and land cover on duration of gonotrophic cycles of Anopheles gambiae (Diptera: Culicidae) in Western Kenya Highlands. J Med Entomol. 2005, 42: 974-980. 10.1603/0022-2585(2005)042[0974:EOMCCB]2.0.CO;2.
Parham PE, Michael E: Modeling the effects of weather and climate change on malaria transmission. Environ Health Perspect. 2010, 118: 620-626. 10.1289/ehp.0901256.
Okech BA, Gouagna LC, Yan G, Githure JI, Beier JC: Larval habitats of Anopheles gambiae s.s. (Diptera: Culicidae) influences vector competence to Plasmodium falciparum parasites. Malar J. 2007, 6: 50-10.1186/1475-2875-6-50.
Lyimo EO, Takken W, Koella JC:Effect of rearing temperature and larval density on larval survival, age at pupation and adult size of Anopheles gambiae. Entomol Exp Appl. 1992, 63: 265-271. 10.1111/j.1570-7458.1992.tb01583.x.
Zakharova NF, Losev GI, Iakubovich VI: The effect of density and temperature on the larval populations of the malaria vector Anopheles sacharovi. Med Parazitol (Mosk). 1990, 1: 3-7.
Dye C: The analysis of parasite transmission by bloodsucking insects. Annu Rev Entomol. 1992, 37: 1-19. 10.1146/annurev.en.37.010192.000245.
Dye C: The epidemiological context of vector control. Trans R Soc Trop Med Hyg. 1994, 88: 147-149. 10.1016/0035-9203(94)90270-4.
Suwanchaichinda C, Paskewitz SM:Effects of larval nutrition, adult body size, and adult temperature on the ability of Anopheles gambiae (Diptera: Culicidae) to melanize Sephadex beads. J Med Entomol. 1998, 35 (2): 157-161.
Aboagye-Antwi F, Tripet F: Effects of larval growth condition and water availability on desiccation resistance and its physiological basis in adult Anopheles gambiae sensu stricto. Malar J. 2010, 9: 225-10.1186/1475-2875-9-225.
Afrane YA, Zhou G, Lawson BW, Githeko AK, Yan G:Life-table analysis of Anopheles arabiensis in western Kenya highlands: effects of land covers on larval and adult survivorship. Am J Trop Med Hyg. 2007, 77: 660-666.
Ellers J, Jervis M: Body size and the timing of egg production in parasitoid wasps. Oikos. 2003, 102: 164-172. 10.1034/j.1600-0706.2003.12285.x.
Pechenik JA, Wendt DE, Jarrett JN:Metamorphosis Is not a new beginning: larval experience influences juvenile performance. Bioscience. 2011, 48: 901-910. 10.2307/1313294.
De Block M, Stoks R:Short-term larval food stress and associated compensatory growth reduce adult immune function in a damselfly. Ecol Entomol. 2008, 33: 796-801. 10.1111/j.1365-2311.2007.00957.x.
Schulte-Hostedde AI, Zinner B, Millar JS, Graham J:Restitution of mass-size residuals: validating body condition indices. Ecology. 2013, 86: 155-163. 10.1890/04-0232.
Kingsolver JG, Huey RB:Size, temperature, and fitness: three rules. Evol Ecol Res. 2008, 10: 251-268.
Honek A:Intraspecific variation in body size and fecundity in insects: a general relationship. Oikos. 1993, 66: 483-492. 10.2307/3544943.
Kingsolver JG, Pfennig DW:Individual-level selection as a cause of Cope’s rule of phyletic size increase. Evolution (N Y). 2004, 58: 1608-1612.
Plaistow SJ, Troussard JP, Cézilly F: The effect of the acanthocephalan parasite Pomphorhynchus laevis on the lipid and glycogen content of its intermediate host Gammarus pulex. Int J Parasitol. 2001, 31: 346-351. 10.1016/S0020-7519(01)00115-1.
Pulkkinen K, Ebert D: Host starvation decreases parasite load and mean host size in experimental populations. Ecology. 2004, 85: 823-833. 10.1890/03-0185.
Moret Y, Schmid-Hempel P:Survival for immunity: the price of immune system activation for bumblebee workers. Science. 2000, 290 (80): 1166-1168. 10.1126/science.290.5494.1166.
Briegel H:Metabolic relationship between female body size, reserves, and fecundity of Aedes aegypti. J Insect Physiol. 1990, 36: 165-172. 10.1016/0022-1910(90)90118-Y.
Zhou G, Pennington JE, Wells MA:Utilization of pre-existing energy stores of female Aedes aegypti mosquitoes during the first gonotrophic cycle. Insect Biochem Mol Biol. 2004, 34: 919-925. 10.1016/j.ibmb.2004.05.009.
Levi T, Ben-Dov E, Shahi P, Borovsky D, Zaritsky A: Growth and development of Aedes aegypti larvae at limiting food concentrations. Acta Trop. 2014, 133: 42-44. 10.1016/j.actatropica.2014.02.001.
Alto BW, Lounibos LP, Mores CN, Reiskind MH: Larval competition alters susceptibility of adult Aedes mosquitoes to dengue infection. Proc Biol Sci. 2008, 275: 463-471. 10.1098/rspb.2007.1497.
Alto BW, Lounibos LP, Higgs S, Juliano SA: Larval competition differentially affects arbovirus infection in Aedes mosquitoes. Ecology. 2005, 86: 3279-3288. 10.1890/05-0209.
Telang A, Qayum AA, Parker A, Sacchetta BR, Byrnes GR:Larval nutritional stress affects vector immune traits in adult yellow fever mosquito Aedes aegypti (Stegomyia aegypti). Med Vet Entomol. 2012, 26: 271-281. 10.1111/j.1365-2915.2011.00993.x.
Paaijmans KP, Huijben S, Githeko AK, Takken W:Competitive interactions between larvae of the malaria mosquitoes Anopheles arabiensis and Anopheles gambiae under semi-field conditions in western Kenya. Acta Trop. 2009, 109: 124-130. 10.1016/j.actatropica.2008.07.010.
Takken W, Smallegange RC, Vigneau AJ, Johnston V, Brown M, Mordue-Luntz A, Billingsley PF: Larval nutrition differentially affects adult fitness and Plasmodium development in the malaria vectors Anopheles gambiae and Anopheles stephensi. Parasit Vectors. 2013, 6: 345-10.1186/1756-3305-6-345.
Araújo M, Gil LH, E-Silva A:Larval food quantity affects development time, survival and adult biological traits that influence the vectorial capacity of Anopheles darlingi under laboratory conditions. Malar J. 2012, 11: 261-10.1186/1475-2875-11-261.
Nasci R:Relationship of wing length to adult dry weight in several mosquito species (Diptera: Culicidae). J Med Entomol. 1990, 27: 716-719.
Lyimo EO, Takken W:Effects of adult body size on fecundity and the pre-gravid rate of Anopheles gambiae females in Tanzania. Med Vet Entomol. 1993, 7: 328-332. 10.1111/j.1365-2915.1993.tb00700.x.
Reisen WK, Aslamkhan M:A release-recapture experiment with the malaria vector, Anopheles stephensi Liston, with observations on dispersal, survivorship, population size, gonotrophic rhythm and mating behavior. Ann Trop Med Parasitol. 1979, 73 ((3)): 251-69.
Detinova T: Age structure of Insect Populations of Medical Importance. Annu Rev Entomol. 1966, 13: 247-
Clements A, Paterson GD: The analysis of mortality and survival rates in wild populations of mosquitoes. J Appl Ecol. 1981, 18: 373-399. 10.2307/2402401.
Smith DL, Battle KE, Hay SI, Barker CM, Scott TW, McKenzie FE:Ross, macdonald, and a theory for the dynamics and control of mosquito-transmitted pathogens. PLoS Pathog. 2012, 8: e1002588-10.1371/journal.ppat.1002588.
Bell AS, Blanford S, Jenkins N, Thomas MB, Read AF:Real-time quantitative PCR for analysis of candidate fungal biopesticides against malaria: technique validation and first applications. J Invertebr Pathol. 2009, 100: 160-168. 10.1016/j.jip.2009.01.006.
Hard JJ, Bradshaw WE, Malarkey DJ:Resource- and density-dependent development in tree-hole mosquitoes. Oikos. 1989, 54: 137-144. 10.2307/3565259.
Grieco JP, Rejmánková E, Achee NL, Klein CN, Andre R, Roberts D:Habitat suitability for three species of Anopheles mosquitoes: Larval growth and survival in reciprocal placement experiments. J Vector Ecol. 2007, 32: 176-187. 10.3376/1081-1710(2007)32[176:HSFTSO]2.0.CO;2.
Russell TL, Lwetoijera DW, Knols BGJ, Takken W, Killeen GF, Ferguson HM:Linking individual phenotype to density-dependent population growth: the influence of body size on the population dynamics of malaria vectors. Proc Biol Sci. 2011, 278: 3142-3151. 10.1098/rspb.2011.0153.
Hawley W:The effect of larval density on adult longevity of a mosquito, Aedes sierrensis: epidemiological consequences. J Anim Ecol. 2010, 54: 955-964. 10.2307/4389.
Ameneshewa B, Service M: The relationship between female body size and survival rate of the malaria vector Anopheles arabiensis in Ethiopia. Med Vet Entomol. 1996, 10: 170-172. 10.1111/j.1365-2915.1996.tb00724.x.
Armbruster P, Hutchinson RA:Pupal mass and wing length as indicators of fecundity in Aedes albopictus and Aedes geniculatus (Diptera: Culicidae). J Med Entomol. 2002, 39: 699-704. 10.1603/0022-2585-39.4.699.
Andersson HI:The effect of sugar meals and body size on fecundity and longevity of female Aedes communis(Diptera: Culicidae). Physiol Entomol. 1992, 17: 203-207. 10.1111/j.1365-3032.1992.tb01011.x.
Gimnig JE, Ombok M, Otieno S, Kaufman MG, Vulule JM, Walker ED:Density-dependent development of Anopheles gambiae (Diptera: Culicidae) larvae in artificial habitats. J Med Entomol. 2002, 39: 162-172. 10.1603/0022-2585-39.1.162.
Ng’habi KR, John B, Nkwengulila G, Knols BGJ, Killeen GF, Ferguson HM: Effect of larval crowding on mating competitiveness of Anopheles gambiae mosquitoes. Malar J. 2005, 4: 49-10.1186/1475-2875-4-49.
Muriu SM, Coulson T, Mbogo CM, Godfray HCJ:Larval density dependence in Anopheles gambiae s.s., the major African vector of malaria. J Anim Ecol. 2013, 82: 166-174. 10.1111/1365-2656.12002.
Dye C:Vectorial capacity: must we measure all its components?. Parasitol Today. 1986, 2: 203-209. 10.1016/0169-4758(86)90082-7.
Dye C, Hasibeder G:Population dynamics of mosquito-borne disease: effects of flies which bite some people more frequently than others. Trans R Soc Trop Med Hyg. 1986, 80: 69-77. 10.1016/0035-9203(86)90199-9.
Hancock RG, Foster WA:Larval and adult nutrition effects on blood/nectar choice of Culex nigripalpus mosquitoes. Med Vet Entomol. 1997, 11: 112-122. 10.1111/j.1365-2915.1997.tb00299.x.
Yaro AS, Touré AM, Guindo A, Coulibaly MB, Dao A, Diallo M, Traoré SF:Reproductive success in Anopheles arabiensis and the M and S molecular forms of Anopheles gambiae: do natural sporozoite infection and body size matter?. Acta Trop. 2012, 122: 87-93. 10.1016/j.actatropica.2011.12.005.
Fernandes L, Briegel H: Reproductive physiology of Anopheles gambiae and Anopheles atroparvus. J Vector Ecol. 2005, 30: 11-26.
Hogg JC, Hurd H:Malaria-induced reduction of fecundity during the first gonotrophic cycle of Anopheles stephensi mosquitoes. Med Vet Entomol. 1995, 9: 176-180. 10.1111/j.1365-2915.1995.tb00175.x.
Yuval B: Effect of body size on swarming behavior and mating success of male anopheles freeborni. J Insect Behav. 1993, 6: 333-342. 10.1007/BF01048114.
Yuval B:Mating systems of blood-feeding flies. Annu Rev Entomol. 2006, 51: 413-440. 10.1146/annurev.ento.51.110104.151058.
Okanda FM, Dao A, Njiru BN, Arija J, Akelo HA, Touré Y, Odulaja A, Beier JC, Githure JI, Yan G, Gouagna LC, Knols BGJ, Killeen GF:Behavioural determinants of gene flow in malaria vector populations: Anopheles gambiae males select large females as mates. Malar J. 2002, 1: 10-10.1186/1475-2875-1-10.
Ponlawat A, Harrington LC:Age and body size influence male sperm capacity of the dengue vector Aedes aegypti (Diptera: Culicidae). J Med Entomol. 2007, 44: 422-426. 10.1603/0022-2585(2007)44[422:AABSIM]2.0.CO;2.
Stoffolano JG, Gonzalez EY, Sanchez M, Kane J, Oquendo AL, Sakolsky G, Schafer P, Yin C, Vela K:Relationship between size and mating success in the blow fly Phormia regina (Diptera: Calliphoridae). Ann Entomol Soc Am. 2000, 93: 673-677. 10.1603/0013-8746(2000)093[0673:RBSAMS]2.0.CO;2.
Vaughan JA: Population dynamics of Plasmodium sporogony. Trends Parasitol. 2007, 23: 63-70. 10.1016/j.pt.2006.12.009.
Billingsley PF, Medley GF, Charlwood D, Sinden RE: Relationship between prevalence and intensity of Plasmodium falciparum infection in natural populations of Anopheles mosquitoes. Am J Trop Med Hyg. 1994, 51: 260-270.
Dong Y, Aguilar R, Xi Z, Warr E, Mongin E, Dimopoulos G: Anopheles gambiae immune responses to human and rodent Plasmodium parasite species. PLoS Pathog. 2006, 2: e52-10.1371/journal.ppat.0020052.
Garver LS, Dong Y, Dimopoulos G: Caspar controls resistance to Plasmodium falciparum in diverse anopheline species. PLoS Pathog. 2009, 5: e1000335-10.1371/journal.ppat.1000335.
Joy DA, Gonzalez-Ceron L, Carlton JM, Gueye A, Fay M, McCutchan TF, Su X:Local adaptation and vector-mediated population structure in Plasmodium vivax malaria. Mol Biol Evol. 2008, 25: 1245-1252. 10.1093/molbev/msn073.
Caragata EP, Rancès E, Hedges LM, Gofton AW, Johnson KN, O’Neill SL, McGraw EA: Dietary cholesterol modulates pathogen blocking by Wolbachia. PLoS Pathog. 2013, 9: e1003459-10.1371/journal.ppat.1003459.
Rivero A, Ferguson HM:The energetic budget of Anopheles stephensi infected with Plasmodium chabaudi: is energy depletion a mechanism for virulence?. Proc Biol Sci. 2003, 270: 1365-1371. 10.1098/rspb.2003.2389.
Preuss J, Jortzik E, Becker K:Glucose-6-phosphate metabolism in Plasmodium falciparum. IUBMB Life. 2012, 64: 603-611. 10.1002/iub.1047.
Rono MK, Whitten MMA, Oulad-Abdelghani M, Levashina EA, Marois E:The major yolk protein vitellogenin interferes with the anti-plasmodium response in the malaria mosquito Anopheles gambiae. PLoS Biol. 2010, 8: e1000434-10.1371/journal.pbio.1000434.
Takken W, Klowden MJ, Chambers GM:Effect of body size on host seeking and blood meal utilization in Anopheles gambiae sensu stricto (Diptera: Culicidae): the disadvantage of being small. J Med Entomol. 1998, 35: 639-645.
Ng’habi KR, Huho BJ, Nkwengulila G, Killeen GF, Knols BGJ, Ferguson HM:Sexual selection in mosquito swarms: may the best man lose?. Anim Behav. 2008, 76: 105-112. 10.1016/j.anbehav.2008.01.014.
Kollien AH, Schaub GA: The development of Trypanosoma cruzi in triatominae. Parasitol Today. 2000, 16: 381-387. 10.1016/S0169-4758(00)01724-5.
Krist AC, Jokela J, Wiehn J, Lively CM:Effects of host condition on susceptibility to infection, parasite developmental rate, and parasite transmission in a snail-trematode interaction. J Evol Biol. 2003, 17: 33-40. 10.1046/j.1420-9101.2003.00661.x.
Logan A, Ruiz-González MX, Brown MJF: The impact of host starvation on parasite development and population dynamics in an intestinal trypanosome parasite of bumble bees. Parasitology. 2005, 130: 637-642. 10.1017/S0031182005007304.
Garrett-Jones C, Grab B:The Assessment of Insecticidal Impact on the Malaria Mosquito ’s Vectorial Capacity, from Data on the Proportion of Parous Females. Bull World Health Organ. 1964, 31: 71-86.
Molineaux L, Shidrawi GR, Clarke J, Boulzaguet J, Ashkar T:Assessment of insecticidal impact on the malaria mosquito’s vectorial capacity, from data on the man-biting rate and age-composition. Bull World Health Organ. 1979, 57: 265-274.
Bayoh MN, Mathias DK, Odiere MR, Mutuku FM, Kamau L, Gimnig JE, Vulule JM, Hawley WA, Hamel MJ, Walker ED: Anopheles gambiae: historical population decline associated with regional distribution of insecticide-treated bed nets in western Nyanza Province, Kenya. Malar J. 2010, 9: 62-10.1186/1475-2875-9-62.
Chitnis N, Schapira A, Smith T, Steketee R:Comparing the effectiveness of malaria vector-control interventions through a mathematical model. Am J Trop Med Hyg. 2010, 83: 230-240. 10.4269/ajtmh.2010.09-0179.
Ijumba JN, Lindsay SW:Impact of irrigation on malaria in Africa: paddies paradox. Med Vet Entomol. 2001, 15: 1-11. 10.1046/j.1365-2915.2001.00279.x.
Ferguson HM, Dornhaus A, Beeche A, Borgemeister C, Gottlieb M, Mulla MS, Gimnig JE, Fish D, Killeen GF:Ecology: a prerequisite for malaria elimination and eradication. PLoS Med. 2010, 7: e1000303-10.1371/journal.pmed.1000303.
We thank members of the Thomas and Read laboratory groups for discussion, Andrew Read for statistical analysis advice, J. Teeple for insectary support, and M. Jones, S. Blanford and D. Sim for laboratory assistance. This research was funded by NIH-R21 (AI096036-01) and NSF-GRFP #DGE1255832. This project was also assisted by a grant from the Pennsylvania Department of Health using Tobacco Settlement Funds. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the National Science Foundation or the National Institutes of Health.
The authors declare that they have no competing interests.
Conceived and designed the experiments: LLM-J CCM MBT. Performed the experiments: LLM-J. Analyzed the data: LLM-J CCM. Wrote the manuscript: LLM-J CCM MBT. All authors read and approved the final submitted version of this manuscript.
Authors’ original submitted files for images
About this article
Cite this article
Moller-Jacobs, L.L., Murdock, C.C. & Thomas, M.B. Capacity of mosquitoes to transmit malaria depends on larval environment. Parasites Vectors 7, 593 (2014). https://doi.org/10.1186/s13071-014-0593-4
- Anopheles stephensi
- Disease ecology
- Food stress
- Host-parasite interactions
- Plasmodium yoelii yoelii
- Trans-stadial effects