- Open Access
Abundance and dynamics of anopheline larvae in a highland malarious area of south-central Ethiopia
Parasites & Vectors volume 5, Article number: 117 (2012)
Malaria is a public health problem in Ethiopia, and increasingly so in highland areas, possibly because of global warming. This study describes the distribution, breeding habitat and monthly dynamics of anopheline larvae in Butajira, a highland area in south-central Ethiopia.
A study of the abundance and dynamics of Anopheles larvae was undertaken at different sites and altitudes in Butajira from July 2008 to June 2010. The sites included Hobe (1817 m.a.s.l), Dirama (1995m.a.s.l.) and Wurib (2196m.a.s.l.). Potential anopheline larval habitats were surveyed once per month in each village. The recorded characteristics of the habitats included habitat type, pH, surface debris, emergent plants, algae, substrate, turbidity, temperature, length, width, depth, distance to the nearest house and anophelines. The Spearman correlation coefficient and Mann–Whitney U test were used to calculate the degree of association between the density of anopheline species and key environmental factors.
Among the different types of habitat surveyed, the Odamo, Akamuja and Assas streams and Beko swamp were positive for anopheline larvae. A total of 3,957 third and fourth instar larvae were collected from the three localities, and they represented ten species of anophelines. These were: Anopheles cinereus (32.5%), An. arabiensis (31.4%), An. chrysti (23%), An. demeilloni (12.2%), An. pretoriensis (0.6%), An. azaniae (0.1%), An. rufipes(0.1%), An. sergentii (0.06%), An. garnhami (0.06%) and An. pharoensis (0.03%). The density of anopheline larvae was highest during the dry months. An. arabiensis was widely distributed, and its density decreased from the lowest elevation in Hobe to the highest in Wurib. The density of An. arabiensis larvae was correlated positively with larval habitat temperature (r = 0.33, p < 0.05) and negatively with depth of larval habitat (r = −0.56, p < 0.05).
Ten species of anophelines were identified, including two known vectors of malaria (An. arabiensis and An. pharoensis), along streams in Butajira. Larvae of An. arabiensis were found in streams at 2200m.a.s.l. This possible expansion of the malaria vector to highland areas indicates an increasing risk of malaria because a large proportion of the Ethiopian population live above this altitude.
Malaria is the leading cause of mortality and morbidity in Ethiopia in areas up to 2500 metres above sea level (m.a.s.l) , although cases have been reported up to 3000m.a.s.l. in some areas . About 70% of the population is estimated to be at risk of infection every year . Transmission of the disease is unstable in many highland areas of the country, where the population has low immunity, and these regions experience malaria epidemics . Anopheles arabiensis, a member of the An. gambiae complex, is the main vector of malaria in the country , while An. pharoensis An. funestus and An. nili represent secondary vectors [6, 7].
The transmission of malaria in high altitude areas of Ethiopia [1, 2, 8] might possibly be due to global warming , land use practices [10, 11] and ecological changes [12, 13] that could favour the breeding and survival of vectors. Warmer weather and increased water temperature enhance malaria transmission in the highlands by shortening the development time from eggs to adult mosquitoes , increasing the number of human blood meals taken by adults, increasing the frequency of egg laying and increasing the survival rate of adult mosquitoes [14, 15]. Increased warmth also shortens the sporogonic cycle of the parasite in the vector, which results in increased intensity of malaria transmission [16, 17]. The continuation of global climate change could therefore allow malaria to expand into the highlands of east Africa , threatening the lives of millions of people.
The existing malaria intervention strategy, which includes indoor residual insecticide spraying, nets treated with long-lasting insecticide, and case management, has been reducing the impact of the disease in Ethiopia. Nevertheless, spread of insecticide-resistant vectors [19, 20] and drug-resistant malaria parasites [21, 22] may result in disease outbreaks. Therefore, control of larvae, which has so far been given little attention, should be reintroduced and implemented together with the existing strategy. Larval control will result in the reduction of the adult mosquito population, subsequently limiting malaria transmission . However, current knowledge of the distribution and dynamics of the aquatic stages of mosquitoes is not adequate. Anopheles mosquitoes breed at the edges of rivers and streams, in temporary rain pools, ponds, dams, drainage ditches, burrow pits, rice fields, swamp margins, roadside puddles and in tree holes close to human dwellings [23–25]. However, mosquitoes differ in their preference for the type, size, turbidity, algal cover and stability of the habitat [26–28]; these factors determine the density, size and disease transmission competence of vectors [25, 29]. Although malaria has become an important health problem in the south-central highland area of Butajira [30, 31], information on the dynamics of the immature stages of the vectors is scarce. The aim of this study was to describe the species distribution and seasonal dynamics of anopheline larvae in the south-central highland area of Butajira. Such information is important in order to implement effective interventions and establish an early warning system for the disease in this country.
The study was undertaken in the Butajira area in the south-central highlands of the Southern Nations and Nationalities Regional State of Ethiopia, which is located 135 km south of Addis Ababa (Figure 1).
For larval sampling, three study sites were selected. These included Hobe (1817 metres above sea level), Dirama (1995m.a.s.l.) and Wurib (2196 m.a.s.l.). The sites are villages close to the Odamo, Akamuja and Assas streams, respectively. They were selected by Health Extension Workers on the basis of habitat availability, accessibility and malaria case reports. Rainfall data for the area (July 2008 to March 2010) was obtained from the National Meteorological Agency of Ethiopia from the only station in Butajira town, which is located 5–20 km from the study areas. On the basis of the previous thirty years of meteorological data from the area (National Meteorological Agency of Ethiopia), the average monthly rainfall is 94.6 mm and the relative humidity is 60.8%, while the average maximum and minimum temperatures are 25.5°C and 11.5°C, respectively. Peak rainfall occurs between July and August, while the lowest level occurs in November and December, with little rain between March and May.
Larval sampling and processing
Sampling for anopheline larvae was undertaken once a month from July 2008 to June 2010. Streams, water wells, small rain pools, pools in hoof- or foot-prints and false banana (Ensete ventricular) leaf axils were surveyed for the presence of larvae, and collections were made by applying a standard sampling procedure [27, 32–34]. Three to ten samples were taken with a soup ladle (350 ml capacity) from each breeding habitat, depending on the size of the habitat and the availability of larvae. In streams, dipping was performed at the edges and stream beds for a distance of 600 to 1600 m, depending on presence of larvae. Along the streams the average distance between two consecutive larval sampling points was 100 m.
Larvae were sorted into culicines and anophelines. All anopheline larvae sampled from each sampling point were identified as 1st, 2nd, 3rd or 4th larval instars, and the corresponding counts were recorded after transferring the larvae from the sampling dipper to white enamel trays. All culicine larvae and the 1st and 2nd anopheline larval instars were discarded. All late anopheline instars (3rd and 4th) were preserved in 70% alcohol after being killed in hot water (ca. 60°C) . In the laboratory, the larvae were mounted in gum-chloral mountant on slides and the species identified on the basis of morphology under a microscope . Furthermore, about 10% of the larval species identified morphologically by the first author (AA) were selected randomly and subjected to reidentification and confirmation by one of the senior and more experienced co-authors (MB). Larvae that were members of the An. gambiae complex were inferred from the results of species-specific PCR conducted on the adults collected from the same study sites (manuscript under preparation). After identification of the late instars, the density of the most common species was expressed as the number of larvae per 100 dips.
Characterization of larval habitat
The types of larval habitats and their characteristics, such as speed, length, width, depth, pH, turbidity, trees nearby (shade), distance to the nearest inhabited house, availability of emergent plants and substrate types were described by technicians and the first author (AA). The flow speed of aquatic habitat was described visually as fast flowing, slowly flowing or stagnant (not flowing). Habitat length, width, depth and distance to the nearest house was measured using measuring tape; shade was recorded as present or absent by observing terrestrial vegetation and/or trees and their branches near the breeding habitat. Emergent plants included both aquatic and immersed terrestrial vegetation . Turbidity was measured by placing a water sample in a clean glass test tube and holding it against a white background; it was classified into four levels: clear, low, medium and high . Substrate type was classified as muddy or sandy. The pH of the water was measured using a portable pH meter, and the water temperature was measured using a minitherm HI 8753 (Romania) digital thermometer.
The data were entered and analysed using SPSS version 16.0 statistical software (SPSS Inc., Chicago, IL). Monthly dynamics of the density of the major anopheline species and the corresponding monthly rainfall data are presented in line charts. The association of the density of the major species with habitat characteristics such as temperature, depth and pH was analysed using the Spearman correlation coefficient, while associations with substrate type (muddy or sandy), turbidity (low or medium), surface debris (present or absent), and surface algae (present or absent) were analysed using the Mann–Whitney U test. The extended Mantel–Haenszel chi-square test for linear trend was used to investigate the trends in major anopheline density at the Hope, Dirama and Wurib sites. The Kappa value was calculated to study the agreement between the researchers in the identification of larval species.
Potential anopheline breeding habitats surveyed from July 2008 to June 2010 in Hobe, Dirama and Wurib villages are presented in Table 1. Among the different types of habitat surveyed, three streams (Odamo, Akamuja and Assas) and one swamp (Beko) were found to harbour anopheline larvae. No anopheline larvae were found in water wells, false banana axils, hoof-prints and most temporary rain pools.
During the study period, 9532 immature anopheline larvae were collected, of which 3171 (33.3%) were 1st instars, 2414 (25.3%) were 2nd instars, 2266 (23.8%) were 3rd instars and 1681 (17.6%) were 4th instars. Among the total sampled, 2302 were from Odamo stream, 37 from a rain pool in Hobe, 1961 from Assas stream, 2294 from Beko swamp, 2925 from Akamuja Stream and 13 from a foot-print in Dirama village.
Of 3947 late (3rd and 4th) instar Anopheles larvae, 3100 (78.5%) were identified to species level (Table 2). The remaining 847 (21.5%) were either lost or could not be identified because of damage to larval parts during processing, or were not mounted on slides for identification. Ten percent (n = 305) of the morphologically identified larvae were selected randomly and subjected to re-identification by a second researcher. There was very good agreement (Kappa = 0.89, p < 0.01) between the researchers in the morphological identification of the anopheline larvae to species level. Anopheles cinereus was the dominant species (32.5%), followed by An. gambiae s.l. (= An. arabiensis in the present work) (31.4%), An. chrysti (23%) and An. demeilloni (12.2%).
Larvae of An. arabiensis were found in the four main breeding sites, with the highest density in Hobe (lowest elevation area) and the lowest density in Wurib (highest elevation area). Larval density declined significantly from the lowland to the highland areas (chi-square for linear trend = 1794, p < 0.01). On the other hand, the density of An. cinereus, An. chrysti and An. demeilloni increased from Hobe to Wurib. The six other species, An. pretoriensis, An. rufipes, An. sergentii, An. azaniae, An. garnhami and An. Pharoensis, were rare; the first five were sampled from the high altitude village while the last species was obtained from the intermediate altitude.
Figure 2 shows the seasonal density of the four common Anopheles species, expressed as the number of larvae per 100 ladle dips, and the corresponding monthly rainfall of the area. An. arabiensis larvae were predominant in Hobe, with high density from December 2008 to April 2009. This was the dry season, when the monthly rainfall was below 40 mm. The density of Anopheles larvae was generally lowest during July and August, corresponding to the highest amount of monthly rainfall. The density of An. demeilloni, An. chrysti and An. cinereus larvae showed similar trends. Among the three villages, Wurib had diverse species of anopheline larvae.
Beko swamp showed the presence of anopheline larvae most frequently (during 16 surveys), followed by Akamuja stream (11 surveys), among the 24 larval surveys (Table 3). The highest average water temperature was recorded along Odamo stream (26°C) and the lowest along Assas stream and in Beko Swamp (23°C). On average, a 1600 m stretch of the Akamuja stream was surveyed once each month for the presence of anophelines, and the shortest habitat distance surveyed, 600 m, was along the Beko swamp. Beko was the deepest permanent breeding habitat (5.3 ± 1.5 cm) and had the closest human inhabitants (20 m).
The density of An. arabiensis late instars increased significantly with increasing habitat temperature (r = 0.33, p < 0.01) and also with decreasing depth of habitat (r = −0.56, p < 0.05) (Table 4). Analysis using the Mann–Whitney U test revealed significantly higher larval density in sandy habitats (z = −3.648, p < 0.01) when compared with habitats with muddy substrate. The density of An. demeilloni was negatively associated with habitat temperature (r = −0.387, p < 0.05). An. arabiensis, An. chrysti, An. cinereus and An. demeilloni were not significantly associated with habitat characteristics such as pH, turbidity, surface debris and surface algae in any of the streams. These habitats supported larval development at their shallow edges, where the speed of flow was low, and on their beds in small and stagnant pools. No emergent vegetation was available along the three streams, but was present in Beko swamp. There was no canopy cover along the anopheline-positive habitats, except for some scattered trees, with no measurable shade on the breeding habitats. All the land close to the breeding habitats was cultivated by farmers.
Ten anopheline species were identified in Butajira. The predominant species was An. Arabiensis, which is the main vector of malaria in the country [5, 6]. Its density decreased from Hobe at the lowest elevation (about 1800 m.a.s.l.) to Wurib at 2200 m.a.s.l. Two of these species (An. gambiae s.l, presumably An. arabiensis and An. chrysti) have been reported from neighbouring villages at about the same altitude . This shows that malaria transmission in the area [30, 31] decreases with increasing altitude. Malaria-related mortality in the area was reported previously to follow a similar altitudinal trend .
We found that An. arabiensis breeds at 2196 m.a.s.l., which is above the altitude reported previously from Kenyan highlands [38, 39]. This suggests that malaria vectors are breeding in highland areas, and global warming  could be one explanation for the expansion of An. arabiensis and An. pharoensis in the Butajira highlands. The 1958 malaria epidemic that affected most highland areas, including at 2600 m.a.s.l. , and a recent report of malaria prevalence of 3.2% at an altitudinal range of 2500 to 3000 m.a.s.l.  could be attributed partly to the expansion of the vectors into areas of higher elevation. Expansion of these vectors to highland areas is a serious threat because most of the Ethiopian population lives in the highlands.
The study revealed that the edges and beds of streams serve as anopheline breeding habitats in the Butajira area during months with low precipitation, as reported previously in the central Rift Valley of Ethiopia  and Western Kenya . Streambed pools were also productive breeding habitats of An. arabiensis during low rainfall seasons in Eritrea . Similar findings have also been documented in other areas of East Africa [32, 42]. Streams can produce large vector populations during dry seasons, and hence larval management that targets streambed pools and stream edges may bring substantial reduction in vector density, and subsequently the incidence of malaria [23, 43], in the south-central highlands of Ethiopia. The absence of larvae along the streams during the rainy months could result from increased stream flow, which carries away immature stages of mosquitoes from their breeding points. Heavy rainfall could also kill larvae directly .
Larvae could not be investigated in temporary water collections formed during rain, in water wells, or in Ensete leaf axils. The absence of larvae from most of the temporary collections of surface water could be due to rapid infiltration of the rain water into the soil and high evaporation. Many permanent water wells did not support anopheline larvae, except culicines. This could be due to their depth, which ranges from 15 to 20 m from the surface, and their water volume, which prevents the entry of direct sunlight and could in turn lower habitat temperature and reduce the availability of the food necessary for larval development. Although temporary habitats may dry out or be flushed out before immature anophelines complete their development , they are unpredictable in occurrence and may make a small contribution to overall adult productivity . In addition, their contribution to vector breeding should not be ignored  because some may support anopheline breeding. Given that we were not able to perform weekly or daily sampling of larvae, we might have missed some potential and temporary breeding habitats between the monthly surveys, and this could have biased our results. We recommend that future studies should be carried out at frequent intervals to produce more detailed information on the dynamics of anopheline larvae.
The anopheline breeding points were shallow edges and beds of streams that were sunlit, slow flowing or stagnant, with or without debris and surface algae. Similar habitat types were reported from the Ethiopian Rift Valley  and Eritrea [25, 42]. The larval density of An. arabiensis increased with increasing habitat temperature and decreasing habitat depth. The occurrence of An. arabiensis larvae in Beko Swamp is an indication of its adaptation to habitats with emergent grass and its expansion to higher elevations, which results in an increased risk of highland malaria. Variability in the pH, turbidity, surface debris and surface algae of the streams did not affect the density of An. arabiensis An. chrysti An. cinereus and An. demeilloni larvae significantly. An. arabiensis is adapted to diverse habitats [25, 27]. The density of An. chrysti An. demeilloni and An. cinereus was not significantly correlated with habitat temperature and depth, which indicates that these anophelines can breed at a greater range of depths and temperatures than An. arabiensis. The lower density of the vector (An. arabiensis) in the Beko and Assas habitats of Wurib village may have been due to the relatively low temperature in the area, which may affect its breeding negatively. However, this area supported more of other anopheline species for much of the study period, when compared with the other three permanent breeding habitats. This may be because the grass present in this habitat might have prevented the loss of immature forms in running water or by the direct splashing of rainfall, and the grass might have served as a resting site for newly emerging and gravid anopheline mosquitoes .
This study has revealed that An. arabiensis breeds on the edges and beds of streams in south-central Ethiopia at elevations up to 2200 m.a.s.l. during the dry months. This observation underlines the importance of streams as breeding habitats of An. arabiensis during dry periods. The edges and pools of streams may be important for maintenance of the Anopheles population and for small-scale transmission of malaria during dry seasons. Hence, policy makers and organizations involved in malaria control activities need to consider options for the management of larvae that target streams during dry seasons. This strategy may reduce An. arabiensis density, and thus reduce the risk of malaria transmission [23, 32, 42]. However, streams might not be the only breeding habitats for anophelines in the area, and hence weekly surveys of all the available habitats and habitat chemistry need to be performed to design a comprehensive and effective larval control strategy.
MoH Federal Democratic Ethiopia: National Five-Year Strategic Plan for Malaria Prevention and Control in Ethiopia. 2006, Addis Ababa, Ethiopia
Graves PM, Richards FO, Ngondi J, Emerson PM, Shargie EB, Endeshaw T, Ceccato P, Ejigsemahu Y, Mosher AW, Hailemariam A, Zerihun M, Teferi T, Ayele B, Mesele A, Yohannes G, Tilahun A, Gebre T: Individual, household and environmental risk factors for malaria infection in Amhara, Oromia and SNNP regions of Ethiopia. Trans R Soc Trop Med Hyg. 2009, 103: 1211-1220. 10.1016/j.trstmh.2008.11.016.
Yeshiwondim AK, Gopal S, Hailemariam AT, Dengela DO, Patel HP: Spatial analysis of malaria incidence at the village level in areas with unstable transmission in Ethiopia. Int J Health geogr. 2009, 8: 5-10.1186/1476-072X-8-5.
Abeku TA, van Oortmarssen GJ, Borsboom G, de Vlas SJ, Habbema JD: Spatial and temporal variations of malaria epidemic risk in Ethiopia: factors involved and implications. Acta Trop. 2003, 87: 331-340. 10.1016/S0001-706X(03)00123-2.
White GB, Tessfaye F, Boreham PFL, Lemma G: Malaria vector capacity ofAnopheles arabiensisandAn. quadriannulatusin Ethiopia: chromosomal interpretation after 6 years storage of field preparations. Trans R Soc Trop Med Hyg. 1980, 74: 683-684.
Krafsur ES, Armstrong JC: An integrated view of entomological and parasitological observations on falciparum malaria in Gambela, Western Ethiopian Lowlands. Trans R Soc Trop Med Hyg. 1978, 72: 348-356. 10.1016/0035-9203(78)90125-6.
Balkew M, Elhassen I, Ibrahim M, Gebre-Michael T, Engers H: Very high DDT-resistant population ofAnopheles pharoensisTheobald (Diptera: Culicidae) from Gorgora, northern Ethiopia. Parasite. 2006, 13: 327-329.
Fontaine RE, Najjar AE, Prince JS: The 1958 malaria epidemic in Ethiopia. Am J Trop Med Hyg. 1961, 10: 795-803.
Costello A, Abbas M, Allen A, Ball S, Bell S, Bellamy R, Friel S, Groce N, Johnson A, Kett M, Lee M, Levy C, Maslin M, McCoy D, McGuire B, Montgomery H, Papier D, Pagel C, Patel J, Antonio J, de Oliveira P, Redclift N, Rees H, Rogger D, Scott J, Stephenson J, Twigg J, Wolff J, Patterson C: Managing the health effects of climate change: Lancet and University College London Institute for Global Health Commission. Lancet. 2009, 373: 1693-1733. 10.1016/S0140-6736(09)60935-1.
Kibret S, Alemu Y, Boelee E, Tekie H, Alemu D, Petros B: The impact of a small-scale irrigation scheme on malaria transmission in Ziway area, Central Ethiopia. Trop Med Int Health. 2010, 15: 41-50.
Shanks GD, Hay SI, Omumbo JA, Snow RW: Malaria in Kenya’s western highlands. Emerg Infec Dis. 2005, 11: 1425-1432.
Ghebreyesus TA, Haile M, Witten KH, Getachew A, Yohannes AM, Yohannes M, Teklehaimanot HD, Lindsay SW, Byass P: Incidence of malaria among children living near dams in northern Ethiopia: community based incidence survey. BMJ. 1999, 319: 663-666. 10.1136/bmj.319.7211.663.
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.
Bayoh MN, Lindsay SW: Effect of temperature on the development of the aquatic stages ofAnopheles gambiaesensu stricto (Diptera: Culicidae). Bull Entomol Res. 2003, 93: 375-381.
Lindsay SW, Martens WJ: Malaria in the African highlands: past, present and future. Bull World Health Organ. 1998, 76: 33-45.
Lindsay SW, Birley MH: Climate change and malaria transmission. Ann Trop Med Parasitol. 1996, 90: 573-588.
Ye Y, Louis VR, Simboro S, Sauerborn R: Effect of meteorological factors on clinical malaria risk among children: an assessment using village-based meteorological stations and community-based parasitological survey. BMC Public Health. 2007, 7: 101-10.1186/1471-2458-7-101.
Pascual M, Ahumada JA, Chaves LF, Rodo X, Bouma M: Malaria resurgence in the East African highlands: temperature trends revisited. Proc Natl Acad Sci U S A. 2006, 103: 5829-5834. 10.1073/pnas.0508929103.
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.
Balkew M, Ibrahim M, Koekemoer LL, Brooke BD, Engers H, Aseffa A, Gebre-Michael T, Elhassen I: Insecticide resistance inAnopheles arabiensis(Diptera: Culicidae) from villages in central, northern and south west Ethiopia and detection of kdr mutation. Parasit Vectors. 2010, 3: 40-10.1186/1756-3305-3-40.
Ketema T, Getahun K, Bacha K: Therapeutic efficacy of chloroquine for treatment ofPlasmodium vivaxmalaria cases in Halaba district. South Ethiopia. Parasit Vectors. 2011, 4: 46-
Schunk M, Kumma WP, Miranda IB, Osman ME, Roewer S, Alano A, Loscher T, Bienzle U, Mockenhaupt FP: High prevalence of drug-resistance mutations inPlasmodium falciparumandPlasmodium vivaxin southern Ethiopia. Malar J. 2006, 5: 54-10.1186/1475-2875-5-54.
Yohannes M, Haile M, Ghebreyesus TA, Witten KH, Getachew A, Byass P, Lindsay SW: Can source reduction of mosquito larval habitat reduce malaria transmission in Tigray, Ethiopia?. Trop Med Int Health. 2005, 10: 1274-1285. 10.1111/j.1365-3156.2005.01512.x.
Omlin FX, Carlson JC, Ogbunugafor CB, Hassanali A: Anopheles gambiaeexploits the treehole ecosystem in western Kenya: a new urban malaria risk?. Am J Trop Med Hyg. 2007, 77: 264-269.
Shililu J, Ghebremeskel T, Seulu F, Mengistu S, Fekadu H, Zerom M, Ghebregziabiher A, Sintasath D, Bretas G, Mbogo C, Mbogo C, Githure J, Brantly E, Novak R, Beier JC: Larval habitat diversity and ecology of anopheline larvae in Eritrea. J Med Entomol. 2003, 40: 921-929. 10.1603/0022-2585-40.6.921.
Minakawa N, Sonye G, Yan G: Relationships between occurrence of Anopheles gambiae s.l. (Diptera: Culicidae) and size and stability of larval habitats. J Med Entomol. 2005, 42: 295-300. 10.1603/0022-2585(2005)042[0295:RBOOAG]2.0.CO;2.
Minakawa N, Mutero CM, Githure JI, Beier JC, Yan G: Spatial distribution and habitat characterization of anopheline mosquito larvae in Western Kenya. Am J Trop Med Hyg. 1999, 61: 1010-1016.
Gimnig JE, Ombok M, Kamau L, Hawley WA: Characteristics of larval anopheline (Diptera: Culicidae) habitats in Western Kenya. J Med Entomol. 2001, 38: 282-288. 10.1603/0022-2585-38.2.282.
Gimnig JE, Ombok M, Otieno S, Kaufman MG, Vulule JM, Walker ED: Density-dependent development ofAnopheles gambiae(Diptera: Culicidae) larvae in artificial habitats. J Med Entomol. 2002, 39: 162-172. 10.1603/0022-2585-39.1.162.
Deressa W, Ali A, Enqusellassie F: Self-treatment of malaria in rural communities, Butajira, southern Ethiopia. Bull World Health Organ. 2003, 81: 261-268.
Deressa W, Fantahun M, Ali A: Malaria-related mortality based on verbal autopsy in an area of low endemicity in a predominantly rural population in Ethiopia. Malar J. 2007, 6: 128-10.1186/1475-2875-6-128.
Minakawa N, Githure JI, Beier JC, Yan G: Anopheline mosquito survival strategies during the dry period in western Kenya. J Med Entomol. 2001, 38: 388-392. 10.1603/0022-2585-38.3.388.
Majambere S, Fillinger U, Sayer DR, Green C, Lindsay SW: Spatial distribution of mosquito larvae and the potential for targeted larval control in The Gambia. Am J Trop Med Hyg. 2008, 79: 19-27.
WHO: Manual on practical entomology in malaria. Part II. Methods and Techniques. Division of Malaria and Other Parasitic Diseases Geneva. 1975
WHO: Malaria entomology and vector control: Part 1. Learner’s Guide, Trial Edition. HIV/AIDS, Tuberculosis and Malaria, Roll Back Malaria. 2003
Gillies MT, Coetzee M: A supplement to the Anophelinae of Africa south of the Sahara (Afrotropical Region). Publications of the South African Institute for Medical Research. 1987, 55: 1-143.
Tesfaye S, Belyhun Y, Teklu T, Mengesha T, Petros B: Malaria prevalence pattern observed in the highland fringe of Butajira. Southern Ethiopia: a longitudinal study from parasitological and entomological survey. Malar J. 2011, 10: 153-
Chen H, Githeko AK, Zhou G, Githure JI, Yan G: New records ofAnopheles arabiensisbreeding on the Mount Kenya highlands indicate indigenous malaria transmission. Malar J. 2006, 5: 17-10.1186/1475-2875-5-17.
Kweka EJ, Zhou G, Lee MC, Gilbreath TM, Mosha F, Munga S, Githeko AK, Yan G: Evaluation of two methods of estimating larval habitat productivity in western Kenya highlands. Parasit Vectors. 2011, 4: 110-10.1186/1756-3305-4-110.
Kenea O, Balkew M, Gebre-Michael T: Environmental factors associated with larval habitats of anopheline mosquitoes (Diptera: Culicidae) in irrigation and major drainage areas in the middle course of the Rift Valley, central Ethiopia. J Vector Borne Dis. 2011, 48: 85-92.
Imbahale SS, Paaijmans KP, Mukabana WR, van Lammeren R, Githeko AK, Takken W: A longitudinal study on Anophelesmosquito larval abundance in distinct geographical and environmental settings in western Kenya. Malar J. 2011, 10: 81-10.1186/1475-2875-10-81.
Shililu J, Mbogo C, Ghebremeskel T, Githure J, Novak R: Mosquito larval habitats in a semiarid ecosystem in Eritrea: impact of larval habitat management onAnopheles arabiensispopulation. Am J Trop Med Hyg. 2007, 76: 103-110.
Mutuku FM, Bayoh MN, Gimnig JE, Vulule JM, Kamau L, Walker ED, Kabiru E, Hawley WA: Pupal habitat productivity ofAnopheles gambiaecomplex mosquitoes in a rural village in western Kenya. Am J Trop Med Hyg. 2006, 74: 54-61.
Paaijmans KP, Wandago MO, Githeko AK, Takken W: Unexpected high losses ofAnopheles gambiaelarvae due to rainfall. PLoS One. 2007, 2: e1146-10.1371/journal.pone.0001146.
We thank Prof. Getachew Tilahun (ex-director of the Aklilu Lemma Institute of Pathobiology) for his encouragement and facilitating field trips during the surveys. We also thank Nega Nigussie, Yohannes Negash, Wossen Sisay and Zerihun Tesfaye for their technical help both in the field and in the laboratory. Rainfall and temperature data of the Butajira area was obtained from The National Meteorology Agency of Ethiopia. This study was financially supported by NUFU (Project No: NUFUPRO-2007/10121).
The authors declare that they have no competing interests.
AA designed the study, collected data in the field, carried out the data analysis and wrote the first draft of the manuscript. TGM participated in the study design, interpretation of the results and editing of the manuscript. MB participated in the conception of the study, in the study design and editing of the manuscript. BL conceived the idea for the study and took part in the study design, data entry and analysis, data interpretation and editing the manuscript. All authors have read and approved the final manuscript.