- Open Access
Predation efficiency of Anopheles gambiae larvae by aquatic predators in western Kenya highlands
Parasites & Vectorsvolume 4, Article number: 128 (2011)
The current status of insecticide resistance in mosquitoes and the effects of insecticides on non-target insect species have raised the need for alternative control methods for malaria vectors. Predation has been suggested as one of the important regulation mechanisms for malaria vectors in long-lasting aquatic habitats, but the predation efficiency of the potential predators is largely unknown in the highlands of western Kenya. In the current study, we examined the predation efficiency of five predators on Anopheles gambiae s.s larvae in 24 hour and semi- field evaluations.
Predators were collected from natural habitats and starved for 12 hours prior to starting experiments. Preliminary experiments were conducted to ascertain the larval stage most predated by each predator species. When each larval instar was subjected to predation, third instar larvae were predated at the highest rate. Third instar larvae of An. gambiae were introduced into artificial habitats with and without refugia at various larval densities. The numbers of surviving larvae were counted after 24 hours in 24. In semi-field experiments, the larvae were counted daily until they were all either consumed or had developed to the pupal stage. Polymerase chain reaction was used to confirm the presence of An. gambiae DNA in predator guts.
Experiments found that habitat type (P < 0.0001) and predator species (P < 0.0001) had a significant impact on the predation rate in the 24 hour evaluations. In semi-field experiments, predator species (P < 0.0001) and habitat type (P < 0.0001) were significant factors in both the daily survival and the overall developmental time of larvae. Pupation rates took significantly longer in habitats with refugia. An. gambiae DNA was found in at least three out of ten midguts for all predator species. Gambusia affins was the most efficient, being three times more efficient than tadpoles.
These experiments provide insight into the efficiency of specific natural predators against mosquito larvae. These naturally occurring predators may be useful in biocontrol strategies for aquatic stage An. gambiae mosquitoes. Further investigations should be done in complex natural habitats for these predators.
Mosquitoes of the Anopheles gambiae complex contain the most efficient vector species for malaria transmission in sub-Saharan Africa . Controlling An. gambiae s.l. (hereafter referred to as An. gambiae) populations is a priority in reducing malaria incidence in both endemic and epidemic areas. In recent years, there has been evidence of increasing resistance of An. gambiae to pyrethroids and dichloro-diphenyl-trichloroethane (DDT) which are used in bed net treatment and indoor residual spraying respectively [2, 3].
There are concerns about the effect of chemicals on non-target organisms, including beneficial and non-beneficial insects [4, 5] fish  and other aquatic mosquito predators . The use of chemicals for malaria vector control may lead to high mortalities of predators in aquatic larval habitats and a subsequent increase in mosquito larval habitat productivity . There is renewed interest in establishing sustainable alternative control methods to complement existing vector control tools using biological resources. Naturally occurring predators have been shown to be a significant ecological factor in the regulation of An. gambiae larvae population [7–10]. Blaustein and Chase  found that predator and larvae associations are likely to reduce the mosquito populations and thus could be an effective management tool for their control. Predators such as notonectids , belostomatids , dytiscid beetles [8, 11], crustaceans , copepods , Odonata [14, 15], wolf spiders (Araneae: Lycosidae)  and amphibians  have been shown to be potential biological control agents against mosquito species in various habitats such as agricultural drainages, rice fields and small water bodies. These habitats are some of the predominant larval habitats in the region for An. gambiae. The concept of importing non-indigenous species for biological control has come under discussion because of the potential undesirable effects of predation, parasitism, and competition on non-target native fauna and flora [18, 19]. However, among the predators described above, notonectids, belostomatids, odonates, Gambusia affins and tadpoles have been found to coexist in natural mosquito breeding habitats in the western Kenya highlands and other parts of the world [7, 16]. These predators have wide range of prey [9, 20], and they are likely to regulate the abundance of larval mosquitoes that share the same habitats . Predator-larvae interactions have been found to be one of the most important factors in the mortality of mosquito larvae in natural habitats .
The present study aimed to evaluate the predation rate and efficacy of five main mosquito larvae predators found in natural habitats against An. gambiae larvae in different habitat types in highlands of western Kenya. Backswimmers (Hemiptera: Notonectidae), tadpoles (Anura: Myobatrachidae), belostomatids (Hemiptera: Belostomatidae), dragon fly nymphs (Odonata: Anisoptera) and Gambusia affins were evaluated in both 24 hour evaluation and semi-field experiments. Evaluation of larval predator efficiency has implications for developing and establishing biological control programs against An. gambiae and other pest mosquito larvae.
Materials and methods
Study area description
24 hour evaluations (0.16256N; 34.74408E) and semi-field (0.16825N; 34.71632E) experiments were conducted in the western Kenyan highlands. Western Kenya is considered to be highly prone to malaria epidemics. Swamp reclamation for agriculture and deforestation for timber and firewood in this area have increased the density of potential breeding habitats . An. gambiae s.s and An. funestus are the primary and secondary malaria vectors in this study area, respectively .
Predator collections were made from natural mosquito breeding habitats in Iguhu village in the western Kenya highlands. Predators were sampled using a standard dipper (350 ml), transferred to basins and transported to the field insectary in Iguhu. They were subsequently placed in individual basins to avoid possible cannibalism as was found to occur in previous study . Predators were introduced in habitats after a 12 hour starvation period for both 24 hour evaluations and semi-field experiments.
Preliminary predation assessment
Preliminary experiments were set-up to ascertain which larval instar was most highly to be preyed upon by each predator. Larvae were introduced in semi-natural habitats at 7:00 Hours and counted after 12 and 24 hours. Each larval instar was tested with each predator species. The potential aquatic predators species found not feeding on any larval stage such as water beetles (Coleoptera: Hydrophilidae) were not considered for semi-field evaluation.
24 hour evaluation and Semi-field experimental designs
These experiments were conducted in two habitat types, with and without refugia, in semi-natural environments. Semi-natural habitats without refugia consisted of two kilograms of soil, 2500 mls of rain water and mosquito larvae (Figure 1A). The habitats with refugia were made up of contents similar to habitats without refugia except that they contained stones and grasses which mimicked hiding structures found in natural habitats for larvae against predators (Figure 1B). The third instar larval was used in the 24 hour evaluation and semi-field experiments. Larval densities of 20, 30 and 40 larvae per basin were used to test the effects of larval density on predator efficiency.
24 hour evaluations were started in either the morning (08:00 h) or evening (18:00 h). The numbers of larvae surviving in each habitat over 24 hours were recorded. Each experiment in each density and habitat type was replicated 10 times.
In the semi-field evaluations (controlled experiments set up in field sites where natural habitats are found), experiments were set-up at 08:00 h and larvae were counted once daily until they were all either consumed or had developed to the pupal stage. Experiments with each predator for each larval density and habitat type were replicated 10 times.
Predator midgut analysis
To verify predation of larvae, ten predators of each species were randomly selected for dissection of the midguts after the experiment. Predators were preserved in 96% ethanol immediately after being taken from experimental basins to stop further digestion of larval DNA . DNA was extracted from the midgut of ten individuals of each of the five predator species tested, and a PCR reaction was performed with An. gambiae s.s. specific primers as described by Scott et al. . Five μl of DNA extract was amplified in a 25 μl PCR-mix containing 1X taq Buffer (Qiagen, Valencia CA, U.S.A.), 2 mM of MgCl2, 0.2 mM of each dNTP, 0.5 ng/μl of primer UN [5'-GTG TGC CCC TTC CTC GAT GT-3'], 0.25 ng/μl of primer GA [5'-CTG GTT TGG TCG GCA CGT TT-3'], 0.73 ng/μl of primer AR [5'-AAG TGT CCT TCT CCA TCC TA-3'], 1 ng/μl primer QD [5'-CAG ACC AAG ATG GTT AGT AT-3'], 0.5 ng/μl primer ME [5'-TGA CCA ACC CAC TCC CTT GA-3'] and 0.05 U/μl HotstartTaq polymerase (Qiagen, Valencia CA, U.S.A.). The PCR reaction was carried out with an initial step of 10 min at 94°C followed by 30 cycles, each consisting of 5 min denaturation at 94°C, 30 s annealing at 50°C and 30 s extension at 72°C; the final cycle products were extended for 10 min at 72°C. Fragments were run through an ethidium bromide 2% agarose gel and photographed under ultraviolet light illumination.
Data analyses were done using PWAS statistics program, version 18, (SPSS Inc., Chicago) for windows and statistica version 6.0. The predation efficiency is defined as mortality rate of larval tested. Daily larval survival among predator species in 24 hour evaluation experiments was compared using chi-square tests, and predation efficiency between predator species and other factors were compared using multivariate analysis of variance (MANOVA) using the Tukey-Kramer HSD test. The comparisons of the larval survival proportions in 12 hours and 24 hours experiments (i.e. morning and evening set-ups) were compared using chi-squared with the adjusted proportion of surviving larvae. In the semi-field experiments, the influence of predator species, habitat type and prey density on larvae survivals were analyzed using MANOVA and differences were compared using the Tukey-Kramer HSD test. The daily survival rates comparisons were computed by the use of one way analysis of variance (ANOVA). The Tukey-Kramer HSD test was used for the analysis of predator species contribution towards developmental time reduction for An. gambiae larvae.
Preliminary predation assessment
Preliminary predation assessments confirmed five of the eleven predators evaluated to be feeding on larvae of An. gambiae Notably, all of the confirmed predators demonstrated highest predation efficiency feeding on third instar larvae. Third instars were used in the subsequent 24 hour evaluation and semi-field experiments.
24 hour evaluation of predation efficiency
The predation efficiency of the five predator species varied significantly in both morning (χ2 = 33.06; d.f. = 12, P < 0.001) and evening (χ2 = 40.54; d.f. = 12, P < 0.0001) experimental set-ups. More larvae were consumed during the night phase hours (Figure 2). Multifactorial analysis of variance (MANOVA) showed that predator species, habitat type and prey density were significant factors affecting in predation efficiency of evaluated predators (Table 1). The chi-squared test for 12 and 24 hour predation differences found no significant differences for each predator in adjusted larval survival proportions (Figure 2).
In the semi-field experiments, there were significant variations among the five predators relative to control in An. gambiae s.s. larvae daily survival rate (Table 2) and pupation rate (Table 3) reduction relative to controls (Pair-wise comparison using Tukey-Kramer HSD quartile value q* = 2.87, P < 0.05).
The variation of factors influencing pupation rates using MANOVA results showed that predator species, prey densities, habitat type, interactions between predator species and prey density, predator species and habitat type had significant influence on pupation rate (Table 4). The prey density × habitat type × predator species and the prey density × habitat type interactions had no significant influence on pupation rates (Table 4). Daily survival rates (predation rates), habitat type and predator species had a significant influence on survival rate reduction (Table 5). The prey density interactions between predator species and prey density, between predator species and habitat type, between prey density and habitat type and among predator, and between prey density and habitat type had no significant influence on the predation rates (Table 5).
In the semi-field experiments, An. gambiae DNA was found in all five of the predators evaluated, confirming actual ingestion of larvae by each species. Of the ten replicates for each predator tadpoles, backswimmers, belestomatids, dragon fly nymphs and Gambusia affins; 30%, 80%, 90%, 90% and 100% of the samples were positive for An. gambiae DNA, respectively (Figure 3).
The results of this study have demonstrated availability of potential biological resources for controlling malaria vectors in the western Kenya highlands. In the 24 hour evaluation experiments, all evaluated predators were shown to be more efficient nocturnal predators. Gambusia affins was most efficient while tadpoles were the least efficient predators among all. Habitat type (with and without refugia) had a significant effect on predation in both 24 hour evaluation and in semi-field experiments, which suggests that habitat refugia may be a significant factor in increased larval survival in habitats with efficient predators. In 24 hour evaluation experiments, the set-up time (i.e. morning or evening) had no effect on overall predation of each predator after 24 hours of observation. In semi-field experiments, larval density did not affect the predation rate, which suggests that predator effectiveness will not be hindered by this factor in the long term. The predation rate of predator species and the survival rates of An. gambiae larvae in habitats with and without refugia were similar to results found by other studies in Kenya and elsewhere [8, 9, 11, 25]. In our experiments, all predators were shown to consume intermediate size prey (third instars larvae) at the highest rate. This may be due to a body size capture and handling trade-off. Similar findings have been reported in Cx. quinquefasciatus and Ae. albopictus predation experiments [26–28].
Differences in larvae body size and shapes of prey are known to influence predators capture and prey strategies [29, 30]. Mosquito developmental stages produce drastic changes in body shape; from a linear first instar larva, through stages increasing in size and finally to a smaller, round pupa. The latter stage was generally the least or not at all preferred by predators, or similarly, the least vulnerable to predation. Behavioral analysis indicated slightly lower capture success and greater handling times for mosquito pupae than first to fourth instar larvae [29, 30].
When the predation rates of the five predators were examined with respect to prey density, the rate of consumption varied among predator species. Daily and overall survival rates in semi-field experiments varied with predator species and not larval density and habitat type. The number of larvae consumed remained high in both habitat types when predators were introduced for evaluation. This reflects the combined effects of searching ability and consumption of An. gambiae larvae by the predator species on a temporal and spatial scale . Of the predators evaluated, Gambusia affins and backswimmers were most efficient in predation. The predators evaluated are known to feed on other mosquito species and they have been reported to coexist in several aquatic habitat types that are readily found in the study region [10, 11, 32]. The effects of these aquatic predators on daily survivorship and pupation rates of An. gambiae s.s larvae in western Kenya are being reported here for the first time. In recent years, culicine and Aedes larval population regulation by dytiscid beetles have been noted in different parts of the world , and this study has demonstrated their efficacy in predation of An. gambiae larvae.
The efficacy of existing malaria vector control methods in several parts of Africa has been reduced due to insecticide resistance selection pressure among insecticide classes used for indoor residual spray and bed net treatments [34, 35]. This has led to a revival of interest in the use of locally available biological resources for sustainable, cost effective control of aquatic mosquito stages. Previous studies in India , Kenya  and Australia  have shown the use of aquatic predators to be effective in reducing malaria vector populations and disease incidence. The integration of aquatic predators in broad scale malaria vector control campaigns may lead to more effective control programs [10, 38–40].
Several studies have demonstrated strong top-down regulation of mosquito larvae by aquatic predators . In the current study, variation in predation rate was mostly associated with the presence or absence of refugia and predator species. Predator candidates for the biocontrol of An. gambiae larvae would ideally be able to increase their population size in the absence of An. gambiae by relying on alternative prey [40, 42]. Given that larval An. gambiae habitats in western Kenya often contain a suite of controphic species, future studies should address competitive advantages of co-occurring species. The positive effect of refugia on larval survivorship also suggests that species specific ability to avoid predation likely exist. While predators clearly have the ability to regulate larval populations, An. gambiae are typically associated with ephemeral pools where predators may not be abundant or present at all. In these scenarios, larval populations are likely regulated by hydroperiod  and/or controphic and intraspecific exploitative competition . In habitats with larval predators, top-down and bottom-up processes are likely important as a joint determinants of community structure.
The larval and predator species composition and abundance in natural habitats are influenced by the ecological characteristics of habitat [44, 45]. From the viewpoint of biological control, the aquatic predators should have a wide range of adaptability in the habitats apart from the predation of target mosquito larvae. Therefore, it may be beneficial to advocate the use of these predators in the wider community as a biological control tool against the aquatic stages of An. gambiae in western Kenya.
Four of five predators evaluated in the 24 hour evaluation and semi-field experiments were able to consume a significant number of mosquito larvae and reduce survival and pupation rates considerably in both habitats, with and without refugia. Our results suggest that, the efficiency of a predator depends on detectability of prey in habitats. The predation risk was shown to be body size (larval instar) dependent. The evaluated predators may play an important role in larval population regulation and thereby impart a positive effect on malaria vector reduction in western Kenya.
Coetzee M, Craig M, le Sueur D: Distribution of African malaria mosquitoes belonging to the Anopheles gambiae complex. Parasitol Today. 2000, 16: 74-77. 10.1016/S0169-4758(99)01563-X.
Hargreaves K, Hunt RH, Brooke BD, Mthembu J, Weeto MM, Awolola TS, Coetzee M: Anopheles arabiensis and An. quadriannulatus resistance to DDT in South Africa. Med Vet Entomol. 2003, 17: 417-422. 10.1111/j.1365-2915.2003.00460.x.
N'Guessan R, Boko P, Odjo A, Chabi J, Akogbeto M, Rowland M: Control of pyrethroid and DDT-resistant Anopheles gambiae by application of indoor residual spraying or mosquito nets treated with a long-lasting organophosphate insecticide, chlorpyrifos-methyl. Malar J. 2010, 9: 44-10.1186/1475-2875-9-44.
Bernadou A, Demares F, Couret-Fauvel T, Sandoz JC, Gauthier M: Effect of fipronil on side-specific antennal tactile learning in the honeybee. J Insect Physiol. 2009, 55: 1099-1106. 10.1016/j.jinsphys.2009.08.019.
Rao JV, Kavitha P: In vitro effects of chlorpyrifos on the acetylcholinesterase activity of euryhaline fish, Oreochromis mossambicus. Z Naturforsch C. 2010, 65: 303-306.
Service MW: Mortalities of the immature stages of species B of the Anopheles gambiae complex in Kenya: comparison between rice fields and temporary pools, identification of predators, and effects of insecticidal spraying. J Med Entomol. 1977, 13: 535-545.
Blaustein L, Chase JM: Interactions between mosquito larvae and species that share the same trophic level. Annu Rev Entomol. 2007, 52: 489-507. 10.1146/annurev.ento.52.110405.091431.
Chandra G, Mandal SK, Ghosh AK, Das D, Banerjee SS, Chakraborty S: Biocontrol of larval mosquitoes by Acilius sulcatus (Coleoptera: Dytiscidae). BMC Infect Dis. 2008, 8: 138-10.1186/1471-2334-8-138.
Carlson JC, Dyer LA, Omlin FX, Beier JC: Diversity cascades and malaria vectors. J Med Entomol. 2009, 46: 460-464. 10.1603/033.046.0307.
Matias JR, Adrias AQ: The use of annual killifish in the biocontrol of the aquatic stages of mosquitoes in temporary bodies of fresh water; a potential new tool in vector control. Parasit Vectors. 2010, 3: 46-10.1186/1756-3305-3-46.
Ohba SY, Kawada H, Dida GO, Juma D, Sonye G, Minakawa N, M T: Predators of Anopheles gambiae sensu lato (Diptera: Culicidae) larvae in wetlands, western Kenya: confirmation by polymerase chain reaction method. J Med Entomol. 2010, 47: 783-787. 10.1603/ME10051.
Su T, Mulla MS: Spatial occurrence and hatch of field eggs of the tadpole shrimp Triops newberryi (Notostraca: Triopsidae), a potential biological control agent of immature mosquitoes. J Vector Ecol. 2002, 27: 128-137.
Marten GG, Reid JW: Cyclopoid copepods. J Am Mosq Control Assoc. 2007, 23 (2 Suppl): 65-92.
Mandal SK, Ghosh A, Bhattacharjee I, Chandra G: Biocontrol efficiency of odonate nymphs against larvae of the mosquito, Culex quinquefasciatus Say, 1823. Acta Trop. 2008, 106: 109-114. 10.1016/j.actatropica.2008.02.002.
Manna B, Aditya G, Banerjee S: Vulnerability of the mosquito larvae to the guppies (Poecilia reticulata) in the presence of alternative preys. J Vector Borne Dis. 2008, 45: 200-206.
Futami K, Sonye G, Akweywa P, Kaneko S, Minakawa N: Diving behavior in Anopheles gambiae (Diptera: Culicidae): avoidance of a predacious wolf spider (Araneae: Lycosidae) in relation to life stage and water depth. J Med Entomol. 2008, 45: 1050-1056. 10.1603/0022-2585(2008)45[1050:DBIAGD]2.0.CO;2.
Minakawa N, Munga S, Atieli F, Mushinzimana E, Zhou G, Githeko AK, Yan G: Spatial distribution of anopheline larval habitats in Western Kenyan highlands: effects of land cover types and topography. Am J Trop Med Hyg. 2005, 73: 157-165.
Ewel JJ, O'Dowd DJ, Bergelson J, Curtis CD, Carla MD, Luis DG, Doria RG, Hobbs RJ, Alan H, Keith RH: Deliberate introductions of species: research needs. Bioscience. 1999, 49: 619-630. 10.2307/1313438.
Simberloff D, Stiling P: Risks of species introduced for biological control. Biol Conservation. 1996, 78: 85-192.
Raghavendra K, Sharma P, Dash AP: Biological control of mosquito populations through frogs: opportunities & constrains. Indian J Med Res. 2008, 128: 22-25.
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.
Sih A: Foraging strategies and the avoidance of predation by an aquatic insect, Notonecta hoffmanni. Ecology. 1982, 63: 786-796. 10.2307/1936799.
Schielke E, Costantini C, Carchini G, Sagnon N, Powell J, Caccone A: Development of a molecular assay to detect predation on Anopheles gambiae complex larval stages. Am J Trop Med Hyg. 2007, 77: 464-466.
Scott JA, Brogdon WG, Collins FH: Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am J Trop Med Hyg. 1993, 49: 520-529.
Chandra G, Bhattacharjee I, Chatterjee SN, Ghosh A: Mosquito control by larvivorous fish. Indian J Med Res. 2008, 127: 13-27.
Allan JD: Effects of predation: predation. Science. 1987, 238: 93-94. 10.1126/science.238.4823.93.
Kesavaraju B, Alto BW, Lounibos LP, Juliano SA: Behavioural responses of larval container mosquitoes to a size-selective predator. Ecol Entomol. 2007, 32: 262-272. 10.1111/j.1365-2311.2006.00846.x.
Kesavaraju B, Juliano SA: Behavioral Responses of Aedes albopictus to a Predator Are Correlated with Size-Dependent Risk of Predation. Ann Entomol Soc Am. 2008, 101: 1150-1153. 10.1603/0013-8746-101.6.1150.
Kerfoot WC: The divergence of adjacent populations. Ecology. 1975, 56: 1298-1313. 10.2307/1934698.
Kerfoot WC: Implications of copepod predation. Limnol Oceanog. 1977, 22: 316-325. 10.4319/lo.1977.22.2.0316.
Mänd T, Tammaru T, Mappes J: Size dependent predation risk in cryptic and conspicuous insects. Evol Ecol. 2007, 21: 485-498. 10.1007/s10682-006-9130-z.
Munga S, Minakawa N, Zhou G, Barrack OO, Githeko AK, Yan G: Effects of larval competitors and predators on oviposition site selection of Anopheles gambiae sensu stricto. J Med Entomol. 2006, 43: 221-224. 10.1603/0022-2585(2006)043[0221:EOLCAP]2.0.CO;2.
Mogi M: Insects and other invertebrate predators. J Am Mosq Control Assoc. 2007, 23: 93-109. 10.2987/8756-971X(2007)23[93:IAOIP]2.0.CO;2.
Yewhalaw D, Bortel WV, Denis L, Coosemans M, Duchateau L, Speybroeck N: First evidence of high knockdown resistance frequency in Anopheles arabiensis (Diptera: Culicidae) from Ethiopia. Am J Trop Med Hyg. 2010, 83: 122-125. 10.4269/ajtmh.2010.09-0738.
Yewhalaw D, Wassie F, Steurbaut W, Spanoghe P, Van Bortel W, Denis L, Tessema DA, Getachew Y, Coosemans M, Duchateau L, Speybroeck N: Multiple insecticide resistance: an impediment to insecticide-based malaria vector control program. PLoS One. 2011, 6: e16066-10.1371/journal.pone.0016066.
Howard AF, Zhou G, Omlin FX: Malaria mosquito control using edible fish in western Kenya: preliminary findings of a controlled study. BMC Public Health. 2007, 7: 199-10.1186/1471-2458-7-199.
Brown MD, Kay BH, Hendrikz JK: Evaluation of Australian Mesocyclops (Cyclopoida: Cyclopidae) for mosquito control. J Med Entomol. 1991, 28 (5): 618-623.
Terlouw DJ, Morgah K, Wolkon A, Dare A, Dorkenoo A, Eliades MJ, Vanden Eng J, Sodahlon YK, ter Kuile FO, Hawley WA: Impact of mass distribution of free long-lasting insecticidal nets on childhood malaria morbidity: the Togo National Integrated Child Health Campaign. Malar J. 2010, 9: 199-10.1186/1475-2875-9-199.
Corbel V, Henry MC: Prevention and control of malaria and sleeping sickness in Africa: Where are we and where are we going?. Parasit Vectors. 2011, 4: 37-10.1186/1756-3305-4-37.
Reichard M, Watters BR, Wildekamp RH, Sonnenberg R, Nagy B, Polacik M, Valdesalici S, Cellerino A, Cooper BJ, Hengstler H, Rosenstock J, Sainthouse I: Potential negative impacts and low effectiveness in the use of African annual killifish in the biocontrol of aquatic mosquito larvae in temporary water bodies. Parasit Vectors. 2010, 3: 89-10.1186/1756-3305-3-89.
Louca V, Lucas MC, Green C, Majambere S, Fillinger U, Lindsay SW: Role of fish as predators of mosquito larvae on the floodplain of the Gambia River. J Med Entomol. 2009, 46: 546-556. 10.1603/033.046.0320.
Frank SD, Shrewsbury PM, Denno RF: Effects of alternative food on cannibalism and herbivore suppression by carabid larvae. Ecological Entomology. 2010, 35: 61-68. 10.1111/j.1365-2311.2009.01156.x.
Wellborn GA, Skelly DK, Werner EE: Mechanisms creating community structure across a freshwater habitat gradient. Ann rev ecology and systematics. 1996, 27: 337-363. 10.1146/annurev.ecolsys.27.1.337.
Aditya G, Pramanik MK, Saha GK: Larval habitats and species composition of mosquitoes in Darjeeling Himalayas, India. J Vector Borne Dis. 2006, 43: 7-15.
Mwangangi JM, Muturi EJ, Shililu J, Muriu SM, Jacob B, Kabiru EW, Mbogo CM, Githure J, Novak R: Contribution of different aquatic habitats to adult Anopheles arabiensis and Culex quinquefasciatus (Diptera: Culicidae) production in a rice agroecosystem in Mwea, Kenya. J Vector Ecol. 2008, 33: 129-138. 10.3376/1081-1710(2008)33[129:CODAHT]2.0.CO;2.
Authors wish to thank field assistants for assisting during the semi-field experiments. We extend our appreciations to Amruta Dixit for constructive comments she gave during write up of this manuscript. This study was funded by the National Institute of Health (NIH) USA (D43 TW01505 and R01 AI-50243 to GY). This article is published with the permission from the Director of Kenya Medical Research Institute.
The authors declare that they have no competing interests.
EJK conceived, designed and performed experiments. GY and AKG coordinated and supervised the study progress. EJK and ZG did data analysis and interpretation. EJK wrote this manuscript. GY, TMG, ZG, AKG, YA and MN edited the manuscript. All authors have read the manuscript and approved for submission.