- Open Access
Evaluating larval mosquito resource partitioning in western Kenya using stable isotopes of carbon and nitrogen
Parasites & Vectorsvolume 6, Article number: 353 (2013)
In sub-Saharan Africa, malaria, transmitted by the Anopheles mosquito, remains one of the foremost public health concerns. Anopheles gambiae, the primary malaria vector in sub-Saharan Africa, is typically associated with ephemeral, sunlit habitats; however, An. gambiae larvae often share these habitats with other anophelines along with other disease-transmitting and benign mosquito species. Resource limitations within habitats can constrain larval density and development, and this drives competitive interactions among and between species.
We used naturally occurring stable isotope ratios of carbon and nitrogen to identify resource partitioning among co-occurring larval species in microcosms and natural habitats in western Kenya. We used two and three source mixing models to estimate resource utilization (i.e. bacteria, algae, organic matter) by larvae.
Laboratory experiments revealed larval δ13C and δ15N composition to reflect the food sources they were reared on. Resource partitioning was demonstrated between An. gambiae and Culex quinquefasciatus larvae sharing the same microcosms. Differences in larval δ13C and δ15N content was also evident in natural habitats, and Anopheles species were consistently more enriched in δ13C when compared to culicine larvae.
These observations demonstrate inter-specific resource partitioning between Cx. quinquefasciatus and An. gambiae larvae in natural habitats in western Kenya. This information may be translated into opportunities for targeted larval control efforts by limiting specific larval food resources, or through bio-control utilizing competitors at the same trophic level.
Larval density and inter-/intra-specific competition have a significant effect on the development and growth rates of mosquito larvae in their aquatic habitats [1–3]. Higher larval densities and competition between An. gambiae s.s. and An. arabiensis result in extended development times, and several studies indicate that this may be a direct result of increased competition for microbial food resources [3–5]. This is particularly noteworthy for An. gambiae s.l. larvae, because they typically inhabit temporary, sunlit pools and increased development times may exceed the duration of aquatic habitats [6–8]. Although larval competition for food resources likely plays an important role in mosquito habitat productivity, relatively few studies have examined the actual partitioning of resources between and among species .
Extensive work has been carried out in characterizing larval habitats of malaria vectors. An. gambiae larvae have been found in a suite of habitat types including abandoned goldmines, brick making pits, ponds, agricultural drainage ditches and puddles [10, 11]. When habitat size falls below ~1 m3, habitat stability and pupal occurrence decreases rapidly ; however, predation of larvae is typically lower and water temperatures are often higher in smaller habitats [12, 13]. Minakawa found a significant correlation between habitat availability and the density of indoor resting mosquitoes in homes .
Stable isotopes represent a valuable yet not fully realized technology for mosquito ecology. Stable isotopes are convenient, naturally occurring tracers that are safe and non-radioactive. There are few examples of stable isotopes being used in vector biology studies; however, Hood-Nowotny and Knols illustrated a wide range of potential applications of stable isotopes in entomological studies and the falling costs and rapidly improving technology of stable isotope analysis . Hood-Nowotny and others demonstrated the use of 13C enrichment as a population marker in a sterile insect technique context . Recently, stable isotopes were used to identify the detrital resource base of Aedes triseriatus in tree holes . Although natural abundance stable isotope studies have been used for several decades to clarify trophic relationships in aquatic ecosystems [17–19], few have specifically looked at resource partitioning among co-occurring species of mosquito larvae. Stable isotope ratios of carbon and nitrogen in invertebrate larvae reflect the larval diet , thus differences in larval diet should be reflected in larval isotope signatures. Winters and Yee were able to demonstrate a higher proportion of δ15N in Aedes albopictus relative to Culex restuans, which was indicative of a stronger association of Ae. albopictus larvae with detritus as a food resource .
Stable isotopes may also provide a means to examine the relationship between larvae and their potential predators. We have previously determined dragon fly nymphs to be one of the most efficient larval predators in the region and in microcosm studies, An. gambiae DNA was detected in dragonfly nymph guts .
In this study, we sought to use naturally occurring stable isotopes to examine inter- and intra-specific resource partitioning among larval mosquitoes. We examined inter- and intra-specific resource partitioning by measuring the carbon and nitrogen stable isotope ratios of larvae sharing the same habitats in western Kenya. We hypothesized that stable isotope ratios of carbon and nitrogen would be able to identify resource partitioning among larval species, and that Anopheles sp. would have isotopic signatures indicative of strong competition for resources, while differences in culicine and anopheline feeding strategies would demonstrate significant resource partitioning among the two genera. Dragon fly nymphs were also collected to identify stable isotopic signatures indicative of food choice. We also conducted a laboratory control study to assess trophic shift, or the difference in δ13C and δ15N composition of larvae compared to their food resources caused by larval metabolism . Finally, to compliment our laboratory and field experiments, we conducted microcosm experiments to demonstrate resource partitioning between genera in relatively simple environments, lacking the complexity of natural habitats.
Anopheles gambiae G3 strain mosquitoes, obtained from Malaria Research and Reference Reagent Resource Center (http://www.mr4.org), were used in the laboratory experiments. Culex quinquefaciatus larvae reared on Tetramin fish food (Spectrum Brands, Inc., Madison, WI, USA) were obtained from the Anthony James laboratory at UCI. Mosquitoes were reared in a walk-in insectary regulated at a relative humidity of 75%, temperature of 26°C, and 12:12 hour light–dark cycle. Mosquito larvae were reared in trays of 100 first-instar larvae per liter of water. Cx. quinquefaciatus and An. gambiae larvae were reared from first instar on fish food in order to determine the trophic shift associated with larval metabolism. This also served to identify differences in isotopic discrimination between the two species. An. gambiae larvae were reared on two different food types, fish food or brewer’s yeast, to test the utility of stable isotopes in determining what food resources were consumed. Third instar larvae were collected into paper cups and approximately 30 minutes were allowed for gut clearance before larvae were rinsed three times with distilled water. Larvae (1 mg/sample) and food samples (2.6 mg/sample) were oven dried overnight at 50°C prior to sample weighing and processing. Three samples were analyzed for each sample type, and for larval samples up to four larvae were pooled to obtain adequate sample weight for analysis.
To examine resource partitioning in a more controlled setting relative to natural habitats, we established an array of microcosms at the Kenya Medical Research Institute (KEMRI) in Kisian, Kenya. The microcosms were constructed using 60 cm diameter washtubs and approximately four liters of local top-soil were added to each. Soil was well mixed prior to addition to keep soil as similar as possible across microcosms. Rainwater was collected from roof runoff in Kisian, Kenya and approximately six liters were added to each microcosm. Microcosms were placed in full sun, screened with a light mesh to prevent oviposition by wild mosquitoes and left for five days to allow for microbial colonization prior to addition of larvae as done previously . Colonies of Cx. quinquefasciatus and Kisian strain An. gambiae s.s. were used for microcosm experiments. The colonies were maintained at the insectary in Kisian. One hundred first instar larvae of each species were added to each microcosm. Water levels were maintained as needed by the addition of rainwater. Third and fourth instar larvae were collected, held for 30 minutes for gut clearance and rinsed three times with distilled water. Detritus and filamentous algae samples were also collected. Filamentous algae were rinsed thoroughly with distilled water prior to desiccation. Samples were dried for 48 hours in a desiccator, and stored with desiccant for transport to the United States further processing. Three replicates were taken for analysis for each sample type, and for larval samples desiccated larvae were ground and up to four larvae were pooled per sample to obtain adequate weight for analysis.
Natural habitat sample collection
Sampling of natural habitats was conducted in Iguhu, Kenya (1,500 m above sea level) in the western Kenya highlands between October 2009 and August 2010. Four selected natural habitats consisted of abandoned goldmines and agricultural drainage ditches. Two habitat types were utilized in order to determine if the degree of resource partitioning was consistent between habitats. Habitats were chosen based on the co-occurrence of two or more mosquito species, at least one of which was An. gambiae. Sample collections were made from each habitat on a single day. Third and fourth instar larvae were collected using a standard dipper and transferred to the lab on ice in water from their respective habitats. Larvae were allowed several hours for gut clearance. Detritus and filamentous algae samples were collected and prepared using the same methods as the microcosm experiments. In 2011, allochthonous C3 plant material cow dung and dragon fly nymphs were collected in addition to larvae. Dragonfly nymphs were held for several hours to allow for gut clearance, and cleaned with distilled water. Samples were prepared for processing following the same methods as those in the microcosm experiments.
Stable isotope analysis
Samples of insects (0.75 mg/sample), detritus (15 mg/sample), algae (1.5 mg/sample), and cow dung (1.5 mg/sample) were analyzed for 13C and 15 N isotopes at the University of California, Davis Stable Isotope Facility. Detritus samples were fumed for 12 hours with concentrated HCl to remove inorganic carbon. Samples were analyzed using a PDZ Europa ANCA-GSL elemental analyzer interfaced to a PDZ Europa 20–20 isotope ratio mass spectrometer (Sercon Ltd., Cheshire, UK). Samples were combusted at 1000°C in a reactor packed with chromium oxide and silvered cobaltous/cobaltic oxide. Following combustion, oxides were removed in a reduction reactor (reduced copper at 650°C). The helium carrier then flowed through a magnesium perchlorate water trap. N2 and CO2 were separated on a Carbosieve GC column (65°C, 65 mL/min) before entering the IRMS. During analysis, samples were interspersed with several replicates of at least two different laboratory standards. These laboratory standards, which were selected to be compositionally similar to the samples being analyzed, have been previously calibrated against NIST Standard Reference Materials (IAEA-N1, IAEA-N2, IAEA-N3, USGS-40, and USGS-41). Each sample’s preliminary isotope ratio was measured relative to the reference gases analyzed with each sample. These preliminary values were finalized by correcting the values for the entire batch based on the known values of the included laboratory standards, which results in corrected δ13C and δ15N values. δ13C values are typically depleted in the heavier isotope (negative) while δ15N values are enriched (positive) relative to international standards . Samples from goldmine three and the microcosm experiments were obtained in 2011 and were analyzed using similar, previously published protocols at the University of California, Irvine Stable Isotope Facility .
Data analyses were carried out using JMP Statistical Discovery Software statistics program (SAS Institute, Cary, NC). For laboratory and microcosm experiments, differences in δ13C and δ15N between Cx. quinquefaciatus and An. gambiae were compared using the paired student’s t-test. Welch’s t-test was used to compare δ13C and δ15N of mosquito species with samples that had unequal variance. Separate tests were used for 13C and 15N. We used analysis of variance (ANOVA) to compare δ13C and δ15N signatures of the co-occurring mosquito species within each habitat for the natural habitat experiments. To examine the relative contribution of potential food sources to target species, we used two and three source mixing models made available by Phillips and Gregg (IsoError, version 1.04; http://www.epa.gov/wed/pages/models/stableIsotopes/isotopes.htm) . Two and three source mixing models allow for estimation of the proportional contributions of food sources to the consumer, and they take into account the mean, standard deviations and sample number for both δ13C and δ15N values.
We conducted experiments in the laboratory to determine whether carbon and nitrogen stable isotope ratios were able to indicate food resource differences and similarities between Cx. quinquefasciatus and An. gambiae larvae. Carbon and nitrogen isotopic signatures of Cx. quinquefasciatus and An. gambiae larvae reared from first instar on Tetramin fish food were not significantly different (δ13C, t = 0.66, df = 4, P = 0.54; δ15N, t = 0.09, df = 4, P = 0.93) (Figure 1). Both larval species were slightly enriched in 13C when compared to the dietary source; however, the carbon enrichment was negligible relative to isotopic discrimination found in microcosm and natural habitat experiments. When An. gambiae first instar larvae from the same laboratory colony were reared on yeast or fish food alone, the carbon and nitrogen isotopic signatures reflected the food source (An. gambiae reared on fish food and yeast; Welch’s t-test; δ13C, t = 104.20, df = 3.0, P < 0.0001; δ15N, t = 19.97, df = 2.2, P = 0.0015) (Figure 1).
In microcosms with equal proportions of Cx. quinquefasciatus and An. gambiae we demonstrated significant differences in both δ13C (t = 2.57, df = 5, P = 0.0499) and δ15N (t = 3.33, df = 5, P = 0.0209) larval body tissue (Figure 2). A two source mixing model suggested that detrital and algal food resources contributed to mosquito tissue of both species. Elemental analysis of An. gambiae body content showed that 85.3% of and 38.1% of body nitrogen was attributed to algal resources, compared to 59.6% and 57.2% in Cx. quinquefasciatus, respectively (Table 1). While Cx. quinquefasciatus appeared to forage about equally between the algal and detrital, the isotopic signatures for An. gambiae was skewed (enriched in 15 N and depleted in 13C), suggesting another, un-analyzed food source may have been contributing to the An. gambiae diet (Figure 2).
Natural habitat sample collections
In all four habitats Culex sp. were consistently depleted in δ13C relative to An. gambiae. However, δ15N did not clearly demonstrate partitioning of resources (Table 2). In all habitats with multiple anopheline species, the anophelines generally clustered together which may reflect the differences in morphology and foraging strategies between the two genera (Figure 3, Table 2). In 2011, we sampled cow dung from near the habitat edge as a third potential food source. A three source, dual isotope mixing model was not able to resolve the relationships of larvae to the potential food sources, suggesting the need for a more comprehensive sampling of potential food sources. Dragonfly nymphs were sampled from the same habitat. In the present study, a mixing model using isotopes of carbon and nitrogen suggested that An. funestus (53%) and An. gambiae (47%) make up the majority of contribution to dragonfly nymphs when compared to Culex sp. (0.1%), (Figure 4).
Laboratory experiments to determine the degree of trophic shift associated with larval metabolism demonstrated low levels of fractionation. Larvae of both An. gambiae and Cx. quinquefasciatus closely reflected their laboratory chow (fish food) diet. Further, when An. gambiae from the same colony were reared on brewer’s yeast, their isotopic composition was indicative of this alternative diet. In the subsequent field experiments, larval genera and species showed variable isotopic compositions consistent with laboratory results, which is indicative of resource partitioning. Kaufman also found that larval carbon and nitrogen isotopes reflected their food source ; however, higher fractionation in larvae was observed when compared to the current study. Stable isotope values of δ15N in Ae. albopictus and Cx. restuans have suggested differences in assimilation or utilization efficiency of nitrogen . In our laboratory studies the low level of trophic shift associated with An. gambiae and Cx. quinquefasciatus feeding allowed us to directly interpret differences in larval isotopic composition as evidence of resource partitioning in the subsequent microcosm and natural habitat experiments.
The partitioning of resources in the microcosms was striking, as evidenced by the relative differences between larval An. gambiae and Cx. quinquefasciatus δ13C and δ15N signatures in microcosms. This is indicative of differential feeding on and/or assimilation of the resources available in microcosms. We hypothesized that partitioning would be less apparent in microcosms when compared to natural habitats. Bacteria, algae and protozoans make up a large proportion of the larval diet and can vary widely in different habitats depending on habitat duration, soil type and allochthonous inputs [26–28]. Due to limitations of food source diversity (i.e. little to no allochthonous plant material, debris, and a limited period for microbial growth), opportunities for larvae to forage on distinct resources are limited. Although the resources in microcosms are presumably less diverse than those found in natural habitats, we found more pronounced differences in microcosms. In addition to a highly diverse menu for larvae in natural habitats, difficulties in identifying resource partitioning in natural habitats using stable isotope analysis alone is also likely due to the long term mixing and nutrient turnover these habitats experience . The resource portioning we found between these commonly co-occurring larvae may offer opportunities for species-specific targeted control. Shading habitats or using algaecides may have a negative impact on larvae associated with emergent algae and photosynthetic microbes, such as An. gambiae. For example, in western Kenya, shading agricultural drainage ditches with Napier grass, a regularly cultivated cow fodder, reduced larval populations of An. gambiae s.l. by 75-80% . As observed, the shading did result in lower temperatures which can affect larval development ; however, it is likely that changes in the aquatic microbial community structure in response to shading may have played a role in the developmental success of larvae, and perhaps oviposition via volatile cues to the gravid female. Some of these interventions may not fully eliminate the target larval species, but in some cases, population reduction may lead to competitive exclusion of the target species.
Dragonfly nymphs had an isotopic composition indicative of a primarily anopheline diet. Although this result should be interpreted somewhat carefully since only three food sources were taken into account (Cx. quinquefasciatus, An. gambiae and An. funestus), studies suggest that culicine larvae may have a competitive advantage over anophelines when dragon fly nymphs are present in the habitat [2, 31]. Our group reported a 70% reduction of An. gambiae s.s. larvae exposed to dragonfly nymphs, and used PCR to confirm the presence of An. gambiae DNA in the nymph gut content . There may be other significant contributors to the nymph diet, but these resources would presumably occupy the same trophic level of the anophelines. Microcosm experiments with known predators and prey may help to evaluate the role of tropic structure in determining mosquito species success and habitat productivity.
In addition to more comprehensive food source sampling, stable isotopes may be useful in studying other aspects of vector ecology. Stable isotope ratios of oxygen and hydrogen are significantly affected by processes of evaporation and infiltration of ground water and rainwater, which will vary based on habitat size [14, 32]. Thus, naturally occurring stable isotopes may potentially be used to address other problems in mosquito ecology, such as tracking adult dispersal from larval habitats. Variation in larval development conditions results in variable adult mosquito size and fitness [3, 33], and determining the aquatic source of “successful” adult mosquitoes will provide useful insights for targeted vector control programs. Commercially available enriched compounds can also be used for such studies in the laboratory and the field. For example, Hamer et al. recently used 15 N-labeled potassium nitrate and 13C-labeled glucose to monitor dispersal of naturally breeding Culex pipiens mosquitoes .
As resource utilization and partitioning among larval mosquitoes can play a major role in determining mosquito species composition and survival in aquatic habitats, many studies have addressed the larval mosquito diet and identified differences in food preference largely based on gut content analyses [35, 36]. The advantage of using a stable isotope method to examine resource partitioning is that carbon and nitrogen isotopes are assimilated into the larval body and contribute to the larval body mass whereas gut content analysis does not necessarily indicate assimilation.
Understanding the trophic relationships in larval mosquito habitats has broad implications for malaria vector management. Habitats with co-occurring culicine and anopheline larvae are common throughout western Kenya, and competitive interactions among co-occurring mosquito larval species play a role in regulating adult populations [1, 2]. The current study demonstrated a high degree of resource partitioning between Cx. quinquefasciatus and An. gambiae larvae. This was most evident in the microcosm studies conducted in Kenya. Further, resource partitioning among culicines and anophelines in natural habitats was also demonstrated. This information can be translated into potential opportunities for targeted larval control efforts by limiting specific larval food resources or through bio-control, utilizing competitors at the same trophic level. Given that culicine and anopheline larvae have been shown to partition the resources available in their habitats, it may be possible to target key food resources as an indirect method of species-specific malaria control.
No specific permits were required for the described field studies. These locations were not protected land, and the field studies did not involve endangered or protected species.
Kirby M, Lindsay S: Effect of temperature and inter-specific competition on the development and survival of Anopheles gambiae sensu stricto and An. arabiensis larvae. Acta Trop. 2009, 109 (2): 118-123. 10.1016/j.actatropica.2008.09.025.
Kweka EJ, Zhou G, Beilhe LB, Dixit A, Afrane Y, Gilbreath TM, Munga S, Nyindo M, Githeko AK, Yan G: Effects of co-habitation between Anopheles gambiae ss and Culex quinquefasciatus aquatic stages on life history traits. Parasit Vectors. 2012, 5: 33-10.1186/1756-3305-5-33.
Gimnig J, Ombok M, Otieno S, Kaufman M, Vulule J, Walker E: Density-dependent development of Anopheles gambiae (Diptera: Culicidae) larvae in artificial habitats. J Med Entomol. 2002, 39 (1): 162-172. 10.1603/0022-2585-39.1.162.
Kaufman M, Walker E, Smith T, Merritt R, Klug M: Effects of larval mosquitoes (Aedes triseriatus) and stemflow on microbial community dynamics in container habitats. Appl Environ Microbiol. 1999, 65 (6): 2661-
Walker E, Lawson D, Merritt R, Morgan W, Klug M: Nutrient dynamics, bacterial populations, and mosquito productivity in tree hole ecosystems and microcosms. Ecology. 1991, 72 (5): 1529-1546. 10.2307/1940953.
Minakawa N, Sonye G, Yan G: Relationships between occurrence of Anopheles gambiae sl (Diptera: Culicidae) and size and stability of larval habitats. J Med Entomol. 2005, 42 (3): 295-300. 10.1603/0022-2585(2005)042[0295:RBOOAG]2.0.CO;2.
Bayoh M, Lindsay S: Temperature related duration of aquatic stages of the Afrotropical malaria vector mosquito Anopheles gambiae in the laboratory. Med Vet Entomol. 2004, 18 (2): 174-179. 10.1111/j.0269-283X.2004.00495.x.
Bayoh M, Lindsay S: Effect of temperature on the development of the aquatic stages of Anopheles gambiae sensu stricto (Diptera: Culicidae). Bull Entomol Res. 2003, 93 (05): 375-381.
Winters AE, Yee DA: Variation in performance of two co-occurring mosquito species across diverse resource environments: insights from nutrient and stable isotope analyses. Ecol Entomol. 2012, 37 (1): 56-64. 10.1111/j.1365-2311.2011.01337.x.
Gimnig J, Ombok M, Kamau L, Hawley W: Characteristics of larval anopheline (Diptera: Culicidae) habitats in Western Kenya. J Med Entomol. 2001, 38 (2): 282-288. 10.1603/0022-2585-38.2.282.
Minakawa N, Mutero C, Githure J, Beier J, Yan G: Spatial distribution and habitat characterization of anopheline mosquito larvae in Western Kenya. Am J Trop Med Hyg. 1999, 61 (6): 1010-
Service M: Mortalities of the immature stages of species B of the Anopheles gambiae complex in Kenya: comparison between rice fields and temporary ponds, identification of predators and effects of insecticidal spraying. J Med Entomol. 1977, 13: 535-545.
Sunahara T, Ishizaka K, Mogi M: Habitat size: a factor determining the opportunity for encounters between mosquito larvae and aquatic predators. J Vector Ecol. 2002, 27: 8-20.
Hood Nowotny R, Knols B: Stable isotope methods in biological and ecological studies of arthropods. Entomol Exp Appl. 2007, 124 (1): 3-16. 10.1111/j.1570-7458.2007.00572.x.
Hood-Nowotny R, Mayr L, Knols BGJ: Use of carbon-13 as a population marker for Anopheles arabiensis in a sterile insect technique (SIT) context. Malar J. 2006, 5 (1): 6-10.1186/1475-2875-5-6.
Kaufman MG, Pelz-Stelinski KS, Yee DA, Juliano SA, Ostrom PH, Walker ED: Stable isotope analysis reveals detrital resource base sources of the tree hole mosquito, Aedes triseriatus. Ecol Entomol. 2010, 35 (5): 586-593. 10.1111/j.1365-2311.2010.01217.x.
Edward D, Newton J, Gilburn A: Investigating dietary preferences in two competing dipterans, Coelopa frigida and Coelopa pilipes, using stable isotope ratios of carbon and nitrogen. Entomol Exp Appl. 2008, 127 (3): 169-175. 10.1111/j.1570-7458.2008.00692.x.
McNeely C, Finlay J, Power M: Grazer traits, competition, and carbon sources to a headwater-stream food web. Ecology. 2007, 88 (2): 391-401. 10.1890/0012-9658(2007)88[391:GTCACS]2.0.CO;2.
Brito E, Moulton T, De Souza M, Bunn S: Stable isotope analysis indicates microalgae as the predominant food source of fauna in a coastal forest stream, south east Brazil. Austral Ecol. 2006, 31 (5): 623-633. 10.1111/j.1442-9993.2006.01610.x.
Post D: Using stable isotopes to estimate trophic position: models, methods, and assumptions. Ecology. 2002, 83 (3): 703-718. 10.1890/0012-9658(2002)083[0703:USITET]2.0.CO;2.
Kweka E, Zhou G, Gilbreath T, Afrane Y, Nyindo M, Githeko A, Yan G: Predation efficiency of Anopheles gambiae larvae by aquatic predators in western Kenya highlands. Parasit Vectors. 2011, 4 (1): 128-10.1186/1756-3305-4-128.
McCutchan JH, Lewis WM, Kendall C, McGrath CC: Variation in trophic shift for stable isotope ratios of carbon, nitrogen, and sulfur. Oikos. 2003, 102 (2): 378-390. 10.1034/j.1600-0706.2003.12098.x.
Wang Y, Gilbreath TM, Kukutla P, Yan G, Xu J: Dynamic Gut Microbiome across Life History of the Malaria Mosquito Anopheles gambiae in Kenya. PLoS One. 2011, 6 (9): e24767-10.1371/journal.pone.0024767.
Wang W, Pataki DE: Drivers of spatial variability in urban plant and soil isotopic composition in the Los Angeles basin. Plant Soil. 2011, 350: 323-338.
Phillips DL, Gregg JW: Uncertainty in source partitioning using stable isotopes. Oecologia. 2001, 127 (2): 171-179. 10.1007/s004420000578.
Merritt R, Craig D, Wotton R, Walker E: Feeding behavior of aquatic insects: case studies on black fly and mosquito larvae. Invertebr Biol. 1996, 115 (3): 206-217. 10.2307/3226931.
Merritt RW, Dadd RH, Walker ED: Feeding behavior, natural food, and nutritional relationships of larval mosquitoes. Annu Rev Entomol. 1992, 37 (1): 349-374. 10.1146/annurev.en.37.010192.002025.
Kaufman M, Wanja E, Maknojia S, Bayoh M, Vulule J, Walker E: Importance of algal biomass to growth and development of Anopheles gambiae larvae. J Med Entomol. 2006, 43 (4): 669-676. 10.1603/0022-2585(2006)43[669:IOABTG]2.0.CO;2.
Wamae P, Githeko A, Menya D, Takken W: Shading by Napier Grass Reduces Malaria Vector Larvae in Natural Habitats in Western Kenya Highlands. Eco Health. 2010, 7: 485-497.
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.
Kumar R, Ramakrishna Rao T: Predation on mosquito larvae by Mesocyclops thermocyclopoides (Copepoda: Cyclopoida) in the presence of alternate prey. Int Rev Hydrobiol. 2003, 88 (6): 570-581. 10.1002/iroh.200310631.
Peterson B, Fry B: Stable isotopes in ecosystem studies. Annu Rev Ecol Syst. 1987, 18: 293-320. 10.1146/annurev.es.18.110187.001453.
Afrane Y, Zhou G, Lawson B, Githeko A, Yan G: Effects of microclimatic changes caused by deforestation on the survivorship and reproductive fitness of Anopheles gambiae in western Kenya highlands. Am J Trop Med Hyg. 2006, 74 (5): 772-
Hamer GL, Donovan DJ, Hood-Nowotny R, Kaufman MG, Goldberg TL, Walker ED: Evaluation of a stable isotope method to mark naturally-breeding larval mosquitoes for adult dispersal studies. J Med Entomol. 2012, 49 (1): 61-70. 10.1603/ME11076.
Garros C, Ngugi N, Githeko A, Tuno N, Yan G: Gut content identification of larvae of the Anopheles gambiae complex in western Kenya using a barcoding approach. Mol Ecol Resour. 2008, 8 (3): 512-518. 10.1111/j.1471-8286.2007.02013.x.
Walker ED, Olds EJ, Merritt RW: Gut content analysis of mosquito larvae (Diptera: Culicidae) using DAPI stain and epifluorescence microscopy. J Med Entomol. 1988, 25 (6): 551-554.
We thank the members of the Yan Laboratory, Prof. Kathleen Treseder and Prof. Timothy Bradley at UC, Irvine for helpful discussion and expertise. This paper is published with the permission of the Director KEMRI.
All authors declare they have no competing interest in this study.
TMG, GY and AKG conceived and designed experiments. TMG, YAA and EJK implemented and supervised the study. TMG did data analysis and interpretation. TMG wrote the manuscript. TMG, EJK and GY revised the manuscript. All authors approved the final version for submission.