Anopheles larval ecology and physicochemical characterization of larval habitats in Dire Dawa: an area colonized by Anopheles stephensi in Eastern Ethiopia
Ephrem Abiy 1,2✉ Email
Teshome Degefa 1 Email
Meshesha Balkew 3 Email
Hailu Merga 4 Email
Denekew Zewdu 1
Alamayo Ayala 2 Email Email
Harrysone Atieli 5
Ming Chieh 5 Email Email
Guofa Zhuo 5 Email
Guiyan Yan 5 Email
Delenasaw Yewhalaw 1,2 Email
1 School of Medical Laboratory Sciences, Institute of Health Jimma University Jimma Ethiopia
2 Tropical and Infectious Diseases Research Center (TIDRC) Jimma University Jimma Ethiopia
3 PMI Evolve Project Former staff of Abt Global Addis Ababa Ethiopia
4 Department of Epidemiology, Institute of Health Jimma University Jimma Ethiopia
5
A
Program of Public Health University of California at Irvine USA
Ephrem Abiy1, 2*, Teshome Degefa1, Meshesha Balkew3, Hailu Merga4, Denekew Zewdu1, Alamayo Ayala2, Harrysone Atieli5, Ming Chieh5, Guofa Zhuo5, Guiyan Yan5 and Delenasaw Yewhalaw1,2
1School of Medical Laboratory Sciences, Institute of Health, Jimma University, Jimma, Ethiopia 2Tropical and Infectious Diseases Research Center (TIDRC), Jimma University, Jimma, Ethiopia
3 Former staff of Abt Global, PMI Evolve Project, Addis Ababa, Ethiopia
4Department of Epidemiology, Institute of Health, Jimma University, Jimma, Ethiopia
5Program of Public Health, University of California at Irvine, USA
*Correspondence: ephremabiy1@gmail.com
Email address:
EA: ephremabiy1@gmail.com
TD: teshedege@gmail.com
MB: meshesha_b@yahoo.com
HM: hailu.etea@ju.edu.et
DZ:danblacknote111@gmail.com
AA: alamayoayele22@gmail.com
HA: etemesi2012@yahoo.com
ZG: zhoug@hs.uci.edu
ML: mingchil@uci.edu
GY: guiyuny@hs.uci.edu
DY yewhalawd@gmail.com
Abstract
Background
Understanding mosquito larval ecology is essential for planning and implementations of vector control strategies. The biotic and abiotic factors affect larval occurrence, density, survival and morphogenesis of mosquitoes. Artificial containers are very suitable larval habitats for some species of Anopheles and Aedes mosquitoes in urban and peri-urban settings. Therefore, this study we identified, mapped and characterized larval habitats, estimated larval density and indices larval habitats. In addition, we determined the species composition of Anopheles mosquitoes and species evenness in urban, peri-urban and rural areas of Dire Dawa city adminstration
Methods
Larval habitats were surveyed and identified monthly for a period of 16 months from February 2023 to December 2024 in urban, peri-urban and rural areas of Dire Dawa. Mosquito larvae and pupae were collected and those larvae identified as Anopheles were fed on fish-food. Emerged adults were provided with 10% sucrose solution and kept under standard conditions in field insectary. Females Anopheles was identified morphologically and further species-specific PCR assay was employed to identify members of An. gambiae s.l. In addition, real time PCR assay was performed to identify An.stephensi and An.arabiensis. Water samples were taken from the larval habitats and the physico-chemical parameters were measured using HANNA Multi-parameter (H198194). Larval habitat diversity, larval abundance and distribution were assessed across the three ecological settings (urban, peri-urban and rural).
Results
A total of 23, 526 larvae and 1,808 pupae of Anopheles mosquitoes were collected from 909 man-made(uncovered cemented cisterns (Brick), plastic sheets, steel drums, Tire tracks, Canal ditch, plastic tanker/Barrel) and natural habitats (River edges, ponds, animal hoof prints and swaps) in urban, peri-urban areas of Dire Dawa. The highest mean larval density (51 larvae per dip) of Anopheles mosquitoes was recorded from peri-urban sites in uncovered water tanker (brick) followed by urban site in brick (46 larvae per dip). Anopheles larvae were not found in steel drum and plastic barrels in rural sites.
A
A total of 2,934 adult Anopheles mosquitoes were emerged from immatures collected from all sites, of which 75% (2194/2934) were An. stephensi, and 22% (636 / 2934) were An.arabiensis. The remaining 3.0% were An. pharoensis, An. coustani, An .amharicus, and An. pretoriensis. Anopheles stephensi, An. arabiensis and An. amharicus shared the same habitats across the three ecologies. Larval density was positively correlated with availability of brick making, proximity to houses, urban setting, presence of competitors /predators, vegetation cover, shade cover and substrate type. But larval presence was not correlated with presence/absence of intervention. Larva/pupa presence were positively correlated with pH (r = 0.264, p = 0.01) and water pressure (r = 0.21, p < 0.05).There was positive correlation among temperature, electrical conductivity (EC), total dissolved solids (TDS), salinity and dissolved oxygen (DO) and negative correlation among temperature with resistivity, pH with mvPH. Larval presence was positively correlated with water salinity and pH.
Conclusion
Anopheles stephensi was the predominant species found in the study area, followed by An.arabiensis, An.amharicu, An.pharoensis, An.coustani and An.pretoriensis. Uncovered water tankers (Bricks) were the most prolific artificial habitats in urban and peri-urban sites followed by plastic sheets while natural habitats such as hoof prints and river basins were the most efficient habitat types in rural and urban settings, respectively. Anopheles stephensi was found in natural habitats of Butuji and legehare rivers from urban site, and in rural sites from man-made habitat of plastic sheet.
Anopheles amharicus larvae was found in plastic sheet, an artificial habitat common in urban, peri-urban and rural areas.
We report here the occurrence of An.stephensi in rural areas, breeding in natural habitats, and co-existing with An. gambiae s.l complex.
Availability of brick making, shorter distance from living houses being in urban ecology were directly correlated with larval density. Habitat abundance and positivity of uncovered water tankers (cemented cisterns or bricks) in urban and pre-urban sites could indicate for feasibility and proper implementation of larval source management.
Key words:
Physico-chemical characterization
mosquito larval habitat
Anopheles stephensi
larval density
Ethiopia
A
Introduction
Anopheles larval habitat characterization refers to the process of identifying and describing the environmental conditions where Anopheles mosquitoes lay their eggs and where their larvae develop [1]. These habitats are typically shallow bodies of water, such as ponds, swamps, streams, rice fields, and containers, which can vary in size and water chemistry. The characterization includes factors like water depth, vegetation, temperature, pH, salinity, and the presence of organic matter, as these conditions influence the survival and development of the larvae [24]. Understanding these habitats is crucial for malaria control, as targeting larval habitats with appropriate interventions, such as biological or chemical control methods, can help reduce mosquito populations and limit the transmission of malaria.
Each malaria vector has its own larval habitat preferences for oviposition that might be linked with their genetics make up or environmental parameters that imposed behavioral variation [5].
The main malaria vector An.gambiae complex prefers to breed in clean water of temporary rain pools, hoof prints, natural pockets of river, saline water (some species).The An.funestus groups breed in emergent vegetation and year round river sides, breeds in rice fields, and irrigation schemes. But the invasive An.stephensi breeds in artificial small containers, jars, bricks, plastic sheets, barrels, iron tankers, steel drums, used tire, ditches sometimes even in polluted water and fountains [6, 7].
Mosquito breeding habitats have specific biotic and abiotic factors that make them suitable to oviposition and support development of mosquito life cycle or metamorphosis [1, 8]. Studies conducted in Ethiopia on Anopheles gambiae s.l. showed that the presence of predators, presence of other competing mosquito species, vegetation coverage and intervention, water temperature, electrical conductivity, pH, dissolved oxygen, salinity, total dissolved solution and turbidity are among the physicochemical factors that affect the density and distribution of mosquito larvae [9, 10]. Similar studies in India, Benin and Iran indicated the correlation between environmental factors and Anopheles stephensi and other Anopheles spss larval abundance [1113].But there are no studies that demonstrated the correlation between water physicochemical parameters and Anopheles stephensi larval abundance in Ethiopia.
Globally, LSM (larval source management) has a proven track record in malaria control and even elimination in some countries. Historical successes in Brazil, Egypt, and parts of Asia relied heavily on environmental management [14]. More recent studies show that larviciding using biological agents like Bacillus thuringiensis israelensis (Bti) can reduce larval and adult vector densities by up to 90% [15].In Ethiopia, larval source management is assumed to be a supplementary strategy besides the principal vector control interventions : use of insecticide treated nets and indoor residual spraying but there is inadequate larval identification and response to guide adoption of LSM approaches [16].
In Ethiopia, pilot studies have shown promising results with significant reductions in larval density and adult vector population reduction through regular Bti and Sumilarv application [17, 18] .Similarly, a study in Eastern part of the country on larval habitat suggested that LSM could be effective for focused control of An.stephensi [19]. In addition, LSM integrated with LLINs and IRS resulted in a notable reduction in malaria incidence in Tolay area of oromia region and in northern Cot^e Idvore’ ([20, 21]. In contrary, Ayana et al, 2025 reported that there was an increased malaria cases irrespective of LSM implementation and need to investigate factors that might affect the success of the intervention [22]. Therefore, this study was conducted to assess the physicochemical characteristics of larval habitats and relation with the ecology, season, larval density and distribution of larvae across the urban, peri-urban and rural areas of Dire Dawa.
Methods
Study area
The study was conducted in three ecological settings, of Dire Dawa City administration (515 Km East of Addis Ababa) which included urban, per-urban and rural settings. The urban setting considered for this study were Dire Dawa University, GTZ and Mermersa. The peri-urban sites were those around the industrial zone including villages of Melka-Jebdu and Jerba. The rural study sites were Aseliso, Hulahulul and Boren. The rural sites were agricultural sites with grass land and scanty vegetation coverage. (Fig. 1)[23]. During the study, larval source management was implemented jointly by the President’s Malaria Initiative (PMI) and the Ministry of Health of Ethiopia in the urban and peri-urban sites [24].
The study area has seasonal malaria transmission, with Plasmodium falciparum and P. vivax the main parasite species,and the vectors include Anopheles arabiensis and An. phareonsis. Recent studies reported the occurrence of An. stephensi as an additional vector [25, 26].
Fig. 1
Map of the study area in Dire Dawa City Administration, 2024.This map was created with ArcGIS Pro 3.5 by Esri.(www.esri.com).
Click here to Correct
Immature mosquito collections and larval habitat characterization
A
Larval surveys were conducted following standard WHO protocol and mosquito immatures were collected using larval dipper of 350ml capacity [27, 28]. Twenty dips were taken from each habitat by dipping at different sides of the habitat. Each larval instar and pupae were counted and recorded. Data were entered monthly using open data kit (ODK) immediately synced to the International Centers of Excellence for Malaria Research (ICEMR) database server to avoid data loss. Geographic coordinates, ecological setting (urban, peri-urban and rural), habitat type (Artificial and Natural) and picture of the habitat on spot, source of water (tap ,rain, or river), distance from household, land use (surrounding environment of the habitat), shade status/degree of exposure to sunlight, water flow, substrate type, distance from brick making, distance from nearest resting habitat and frequency of utilization, presence of canopy cover, presence of intervention, type of intervention and presence of predators were included among others. Mosquito immatures were classified as Anopheline, Aedine and Culicine based on morphological observations. The species composition of Anopheles mosquitoes were determined from adults reared from larvae or pupae.
Molecular species identification
The collected adult Anopheles were identified first to species using the morphological keys of Coetzee et al. 2020 and Glich et al. 1995 [29, 30]. Those identified as An. gamabiae s.l were further re-identified using species specific PCR. Real time PCR was also employed to confirm An. stephensi at the Genomics Laboratory of Jimma University Tropical Infectious and Diseases research center (TIDRC).
DNA extraction and amplification
A
DNA extraction was done using both the Chelex and Extracta protocols [31]. PCR amplification was carried out according to the methods of Scott et al. [32] using species-specific primers for An. arabiensis (AR: 5′-AAG TGT CCT TCT CCA TCC TA-3′) and An. amharicus, formerly Anopheles quadriannulatus B (QD: 5′-CAG ACC AAG AGA GAT GGT TAG TAT-3′)[33]. Anopheles gambiae (GA: 5′-CTG GTT TGG TCGGCA CGT TT-3′) and a universal primer (UN: 5′-GTGTGC CCC TTC CTC GATGT-3′). Then the amplicon was loaded on a 2% agarose gel stained with ethidium bromide and run for gel electrophoresis. Anopheles arabiensis from the Sekoru insectary colony and An. amharicus specimens from Arjo area were used as positive controls [34].
Physico-chemical analysis of larval habitats
Physicochemical parameters of each habitat were measured using HANNA (H198184) instrument of multi-parameters (Multiparameter Meter, Hanna Instruments Inc., Woonsocket, Rhode Island, 02895, USA) after proper calibration of the instrument with the manufacturer’s probe solutions that was provided together with multi-parameter. Measurements of acidity and alkalinity (pH/mV), oxidative reduction potential (ORP), electrical conductivity (EC), total dissolved solids (TDS), resistivity, salinity, seawater, dissolved oxygen (DO), atmospheric pressure and temperature were taken from 102 larval habitats.
Statistical Analysis
Data were entered in to SPSS version 27 and descriptive analysis was used to compare frequency, distribution, and magnitude of Anopheles larvae in the urban, peri-urban and rural interfaces. Multiple regression analysis was used to compare the effect of physicochemical variables of habitats on the Anopheles larva presence / absence, and larval density. Principal component analysis (PCA) was performed and Spearman’s correlation coefficient were used to dictate the correlation among water physicochemical parameters and larval presence and density. Cross tabs, Pearson’s Chi squared and Fisher’s exact tests were used to compare the association of presence/absence of predators, and presence/absence of intervention on larval presence and density. Larval density was calculated as number of larvae/pupae per 20 dips per habitat. Habitat index, container Index, and breteau Index were used to compare the larval density and abundance in each of the three ecological settings. Standard mathematical formulae [35] were used to calculate the larval indices as follows: 1) House Index (HI) was calculated by dividing the number of houses infested to the total number of houses inspected and multiplied by 100.
2) Container Index (CI) was calculated by dividing the number of containers positive for anopheles larvae to the total number of containers inspected and multiplied by 100. 3)Breteau Index (BI) was calculated by dividing the number of positive containers to total number of houses and multiplied by 100.
Results:
Larval habitat abundance and indices
A total of 909 habitats were surveyed in Dire Dawa, of which, 857 (94.3%) were artificial breeding habitats such as bricks/concrete water tankers, plastic sheets, plastic water tankers (barrels), steel drums, canal ditches and tire tracks while the other were natural habitats such as river edges, animal hoof prints, ponds and swamps (Table 1)( Fig. 2).
Fig. 2
Types of larval habitats in Dire Dawa City Administration
Click here to Correct
In Dire Dawa, uncovered water tanker/Cemented cistern or brick were the predominant habitat type 520/909 (57.2%) and followed by plastic sheet 221/909 (24.3%). Anopheles larvae were found in 297 (32.7%) of the houses surveyed and in 400(43.7%) of containers inspected in Dire Dawa (Table 1 and Table 2).
Table 1
Larval habitat abundance, proportion and larval presence positivity in urban, peri-urban and rural sites in Dire Dawa City Administration (2024)
Habitat type
n*(%)
positive for Anopheles(%)
 
Brick concrete water tanker
520 (57.2)
187(36.0)
Drainage ditch/canal
10 (1.1)
2(20.0)
River edge
25 (2.8)
14 (56.0)
Tire Tracks
6 (0.7)
4 (66.6)
Foot Print
11(1.2)
9 (81.8)
Swamp
6 (0.6)
4 (66.7)
Pond
10 (1.2)
4 (40.0)
Steel Drum
25 (2.8)
5 (20.0)
covered drinking water tanker/Birca
60 (6.6)
8 (13.3)
Plastic sheet
221(24.3)
58 (26.2)
Other (small water containers, pipe, foundation, septic tank)
15 (1.7)
2 (13.3)
Total
909 (100)
297 (33)
As depicted below in Table 2, about 42, 34%, and 24% of larva habitats surveyed were from urban, peri-urban and rural settings respectively.
Nine hundred-fifteen containers were found in houses surveyed with high number of containers 350(38.3%) in urban, 316 (34.5%) in peri-urban and 249 (27.2%) in rural sites respectively.
Table 2
Larval indices and the distribution of habitats in urban, peri-urban and rural ecologies in Dire Dawa City Administration, Ethiopia
Site
Houses
surveyed
n (%)
# positive
houses
n (%)HI*
#containers
n (%)
#Positive
containers CI*
n (%)
Breteau
Index BI*
n (%)
Urban
385 (42.0)
133 (34.5)
350 (38.3)
175 (50)
175 (45.5)
Peri-urban
307 (34.0)
107 (34.8)
316 (34.5)
130 (41.1)
130 (42.3)
Rural
217 (24.0)
57 (26.3)
249 (27.2)
95 (38.2)
95 (43.8)
Total
909 (100)
297 (32.7)
915 (100)
400 (43.7)
400 (44)
N = number, *HI-House index, CI:Container index and BI: Breteau index
A total of 23,526 larvae and 1,808 pupae of Anopheles mosquitoes were collected from 909 habitats in urban, peri-uban, and rural areas of Dire Dawa. There was a difference in monthly larval density at each ecology and habitat type.
Paired-samples t-tests indicated that there is a significant difference in the density of Anopheles larvae across months ( t=-3.27, n = 34, p = .003) with a mean difference of − 6.98 larvae per dip (95% CI: −11.32 to − 2.63) and a medium effect size (Cohen’s dₚ = 0.55), indicating a moderate seasonal influence on larval abundance. Similarly, the paired-samples t-test showed that larval density differed significantly between ecological categories ( t = − 5.37, n = 34, p < .001) with a mean difference of − 11.52 larvae per dip (95% CI: −15.88 to − 7.16) and a large effect size (Cohen’s dₚ = 0.91). A stronger effect was observed across ecological settings, which suggests that ecological characteristics exert a greater influence on Anopheles larval density than seasonal variation alone.
In urban sites, where larval density was estimated to be constant throughout the year, the monthly Anopheles larval density was high in October and lowest in May in urban settings. In peri-urban sites, Anopheles larval density was highest in October and absent in May. In rural sites, larval density was peak in February and absent from June to October (Fig. 3), showing that there was a significance difference between months and ecology, where some ecological settings were absent in some months (t = 7.38, n = 34, p < 0.001), suggesting that temporal and ecological factors together influence larval dynamics.
Fig. 3
Mean Anopheles larval density by month and ecological setting in in Diredawa, Ethiopia (2024)
Click here to Correct
The mean larval density were highest in peri-urban site in artificial uncovered water tanker (51 larvae per dip), followed by natural ditch in urban site (46 larvae per dips), and lowest in rural and urban plastic sheet habitats, and lowest in plastic sheets in Urban sites (0 per dip/habitat).
Habitat type differs significantly across the season (t = 6.36, n = 77, p < 0.001) indicating seasonal variation affect habitat type. Anopheles larvae abundance and density were significantly different in different types of habitats (t = -4.8, n = 77, p < 0.001). There is a significant association between ecological settings (urban, peri-urban, and rural) and habitat types (t=-6.027, n = 77, p < 0.001). Seasonal variations, specifically dry and long rain seasons, strongly affect larval abundance and density (t = -4.95, n = 77, p < 0.001).
As construction of houses and production of bricks are common in urban and peri-urban areas of Dire Dawa, Anopheles larvae population density was high throughout the year in the artificial water storages, regardless of the seasonal variation in the urban and peri-urban ecologies. But in rural areas, Anopheles larval density increased following the short rain season (Fig. 4).
Fig. 4
Mean Anopheles larval density by season, habitat type and ecological setting in Dire Dawa City Administration, Ethiopia (2024)
Click here to Correct
Proximity of larval habitats to houses 80% (732 /909) were located in less than 50 meters from nearby living houses, 65% (589/909) of the habitats were at 0–50 meters away from brick making houses. Additionally, 56% (509/909) of the habitats surveyed were surrounded by construction site or land-use for brick making. Half of the total habitats substrate were made of concrete with 35% (320/909) covered with plastic sheet or lining and about 11% of the habitat substrate was mud. The results of logistic regression analysis indicated that man-made/artificial habitat and grass land of surrounding were significant predictors for the presence of Anopheles larvae/pupae (Table 3).
Table 3
Logistic regression analysis for the assessment of predictors of mosquito larval occurrence in urban, peri-urban and rural areas of Dire Dawa City Administration, Ethiopia
Variable
Estimate
ChiSquare
Prob > ChiSq
Lower 95%
Upper 95%
Eco setting[peri-urban]
0.249
1.34
0.2469
-0.172
0.675
Eco setting[rural]
-0.251
1.03
0.3103
-0.747
0.228
Habitat category[man made]
2.222
5.76
0.0164
0.823
4.805
Landuse2[construction]
0.199
0.19
0.6608
-0.638
1.148
Landuse2[farmland]
0.407
0.7
0.4023
-0.534
1.386
Landuse2[grassland]
1.504
4.45
0.0349
0.165
3.019
Landuse2[house]
-0.588
1.05
0.3058
-1.717
0.557
Landuse2[road]
-1.943
1.28
0.2587
-5.645
1.139
Vegetation Cover
-0.008
0.23
0.6321
-0.041
0.027
Shade coverage
0.001
0.01
0.942
-0.018
0.018
LSM
-0.525
0.88
0.3494
-1.693
0.540
Distance to house
-0.004
2.26
0.1325
-0.010
0.001
Distance to resting
0.000
0
0.988
-0.005
0.005
Distance to brick making
-0.003
1.69
0.1932
-0.009
0.001
Habitat length
-0.031
0.86
0.3531
-0.107
0.032
Both natural and man-made habitats were having same correlation with Anopheles and Culex larval presence and Aedes larval density. But rural sites were highly correlated with Aedes larval density than urban and peri-urban and sites. Conversely, An.stephensi larval presence were correlated with urban and peri-urban than rural sites. Similarly, the result of the PCA indicated positive correlation between land use (surrounding environment) and larval presence, habitat with in Grass-land environment, farm land, and availability of shrubs and in or near household were correlated with presence of Anopheles, Culicines and Aedeines larvae (Fig. 5 and Fig. 7).
Fig. 5
Variables grouped by PCA, Variables grouped together means they are correlated. PCA 1 (red bars) includes Anopheles and Culex larval presence, PCA 2 (green bars) includes Aedes larval density, PCA 3 (blue bars) includes An. stephensi presence
Click here to Correct
The results of principal component analysis (PCA) indicated that Larva/pupa presence were positively correlated with pH (r = 0.264, p = 0.01) and water pressure (r = 0.21, p < 0.05). Similarly, there was a significant positive correlation with temperature and conductivity (EC) (r = 0.31, p = 0.02), Total dissolved solids (TDS) ((, r = 0.38, p = 0.02) and, Salinity (r = 0.38, p = 0.02)(Fig. 6).
EC, TDS and salinity were positively correlated (r = 1, p = 0.00). Also Dissolved oxygen (DO) and water pressure (psi) showed positive correlation (r = 0.211, p = 0.040). But, temperature and resistivity had negative correlation (r = − 0.494, p = 0.00), pH and mvpH also had negative correlation (r= -0.993, p = 0.00). In addition, there was no significant correlation between temperature and Dissolved oxygen, Dissolved oxygen and Salinity (p > 0.05). Oxydative reduction potential (ORP) was the variable that had no-significant correlation with the dependent variable (Fig. 6).
Fig. 6
Results of Spearman’s correlation among predictors and occurrence of mosquito immatures.
Click here to Correct
The PCA result also showed that Anopheles larval/pupa density were positively correlated with shorter distance to living house and distance to resting places, availability of brick making, vegetation cover, shading, pH and salinity (Fig. 7).
Fig. 7
Correlation of Environmental and physicochemical parameters with presence and density of mosquito immatures.
Click here to Correct
The Pearson’s Chi squared test revealed a significance association between the predators and competitors presence on larval presence (ꭓ2=16.3, df:1, p < 0.01). Larvae were found more often in locations where predators were present (66.7%) compared to locations without predators (31.5%).The result indicated that larva presence and predator presence are not independent, and predator and competitors presence is significantly associated with higher likelihood of larva presence.
But, the Pearson’s chi squared test indicated that there is no significance association between presence of intervention (larviciding) and larval presence (ꭓ2 =2.66, df = 2, p = 0.264), the presence of intervention doesn’t affect larval presence (p > 0.05).
As depicted in Fig. 5, Clear water had positive correlation with presence of Anopheles stephensi. Moreover, water salinity, water temperature and water pH were among the strong predictors of the presence of An.stephensi larvae (Fig. 7)
The analysis of water salinity revealed that the majority of samples (88.3%) fell within the valid or normal salinity range, indicating predominantly fresh water conditions across the study area. Only 11.7% of the total 873 samples exhibited measurable salinity above 0 ppt. Among these, slightly saline water (2–3 ppt) was the most frequent, accounting for 5.6% of the samples, followed by low salinity levels between 0–1 ppt (2.4%). Higher salinity categories were much less common, with only 1.8% of samples in the 4–5 ppt range and less than 1% exceeding 6 ppt. Extremely saline samples (> 30 ppt) represented just 0.3% of the total (Table 5).
A
Table 4
Results of water salinity measurements of larval habitats in Dire Dawa City Administration, Ethiopia
Salinity (ppt)
N (%)
Category
0–1
21(2.4)
Very low salinity (freshwater).
2–3
49 (5.6)
Slightly saline (brackish tendency).
4–5
16 (18)
Moderately saline.
6–7
6 (0.7)
Noticeably saline.
8–9
2 (0.2)
Highly saline.
10–14
5 (0.6)
Very saline—approaching seawater levels.
> 30
3 (0.3)
Extremely saline (seawater or hyper-saline).
The water physicochemical analysis indicated that water temperature of larval habitat ranged from 19°C to 34°C. The water temperature for 80.4% (82/102) of larval habitats was from 20–27°C (room temperature), 16.6% (17/102) of the larval habitats water temperature was from 28–34°C, and for the remaining 2.9% (3/102) of habitats, the water temperature recorded were 19°C.
The water pH values measured and categorized in to three as acidic water (pH value < 7), neutral (pH value = 7) and alkaline (pH > 7).The result indicated that 86.3% (88/102) of habitats were alkaline, 12.7% (13/102) were neutral and 0.98% (1/102) of habitats were acidic.
Five types of substrate were considered (Gravel, Sand, Mud, Concrete and Plastic cloth) and the result of the study showed that 50% (52/102) of the habitats substrate were made of concrete and 25% (26/102) of the habitat substrate were made of plastic cloths/container. Fisher’s exact test indicated that there is a statistical association between larval presence and substrate type (ꭓ2=26.6, df = 6, p < 0.01).
Species diversity and evenness
A total of 2,934 adult Anopheles mosquitoes were emerged from immatures from all sites and An. stephensi was the predominant species 75% (2194/2934) followed by An.arabiensis 22% (636 / 2934), An. pharoensis 1.7% (50/2934), An. coustani 1.3% (38/2934), An.amharicus 0.3% (9/2934) and An. pretoriensis 0.2% (7/2934), respectively (Table 5).
Table 5
Anopheles mosquito species diversity, richness and evenness by habitat type in Dire Dawa City Administration, Ethiopia
Habitat
Species
Abundance
Proportion
(pi)
Lnpi
pi*Lnpi
Shannon
diversity index (H')
Species
eveness (J): H/LnS
Species
richness (S)
Artificial
An.stephensi
1873
0.91
-0.1
-1.
0.1
0.3
3
An.arabiensis
185
0.09
-2.4
-0.2
0.22
An.amharicus
8
0.00
-5.6
-0.1
0.02
An.pretoriensis
0
0
0
0
0
An coustani
0
0
0
0
0
An.pharoensis
0
0
0
0
0
Total
 
2066
1
0
0
0.33
Natural
An.stephensi
321
0.37
-1
-0.4
0.4
0.6
6
An.arabiensis
451
0.52
-1
-0.3
0.3
An.amharicus
1
0.00
-7
-0.01
0.01
An.pretoriensis
7
0.01
-5
-0.04
0.04
An.coustani
38
0.04
-3
-0.14
0.14
An.pharoensis
50
0.06
-3
-0.16
0.16
Total
 
868
1
0
0
1.06
Further paired samples t-test revealed significant differences in species composition between habitat categories (man-made and natural). Abundance differed significantly between habitats (t = − 2.22, n = 13, p = 0.004), with man-made habitats having higher Anopheles abundance. Significant habitat-related differences were also observed in Shannon diversity (t = 10.03, n = 13, p < 0.001), species evenness (t = 5.70, n = 13, p < 0.001), and species richness (t = 5.7, n = 13, p < 0.001), indicating greater diversity, more even species distribution, and higher species richness in natural habitat category.
An. stephensi, An.arabiensis and An.amharicus found to share same artificial habitats more specifically plastic sheets (in urban, peri-urban and rural), and in natural habitats more specifically a river Butuji river found in Dire Dawa city.
Molecular identification of An.gambiae s.l and An.stephensi
Of total 645 An.gambiae s.l sub-samples that were tested using PCR, 98.6% (n = 636) were found to be An.arabiensis, and 1.4% (n = 9) were An.amharicus. Similarly, real time PCR (qPCR) was performed to confirm morphologically identified An.stephensi specimens from An.gambiae s.l, and accordingly 200 sub-samples of each An.stephensi and An.arabiensis were proved distinct.
Discussion
Our study indicated that An.stephensi immature stages are abundant throughout the year regardless of the seasonal variation of rain, with peak larval densities in dry season in all the three ecological settings. Similar results were reported by Yared et al. 2023 [19] in Somali region, Jigjiga town, where larvae were present in dry season. This might be because both Dire Dawa city and Jigjiga town have similar in expansion of urbanization and construction which support the breeding of immature stages in dry season.
Moreover, larval habitats were found in less than 50m from living houses and brick making houses were located at distance ranging from 0-50m away from living houses. This finding has an important implication as far as supporting the life cycle of Anopheles and maximizing the probability of man vector contact. This finding could help to implement targeted larviciding, larval source management in habitats that harbor the immature stages throughout the year.
The other finding of our study is the co-existence of An.stephensi with An.gambiae s.l, sharing same both artificial and natural habitats and ecological settings. An.stephensi were known to be vector of malaria mainly in urban settings and in an artificial habitats [6, 11], but in our study An.stephensi were found in rural ecology in Dire Dawa, and in river edges of Legehare and Butuji rivers, the two known and permanent habitats in the city. This further remarks the expansion of An.stephensi to rural sites and adapting natural habitats to breed. Our finding coincides with a study done in South east Iran where An.stephensi were found in river edges with sandy substrates [12].
Unexpected result of our research in the study area is, the presence of An.amharicus which were reported in South Western, Central, Western, [33, 36] and most recently in North west part of the country, here in the Eastern part of Ethiopia ,we report the existence of An.amharicus in urban,peri-urban and rural ecologies of Dire Dawa for first time. Additionally, An.coustani, An.pretoriensis and An.phareonsis were found in the study area.
This finding recommends further investigation on the role of An.amharicus and the other vectors in malaria transmission in Dire Dawa and their general bionomics including its susceptibility to current insecticides in use are important.
Water physicochemical analysis indicted that larval distribution and abundance was not significantly different as the measures of the water physicochemical vary. Also a change in Ecology, season, habitat category doesn’t indicate a difference on presence and absence of larvae and density of larvae in man-made habitats. Our study align with the study conducted in Northern Iran which indicated no significance difference in physicochemical characteristics and larval density [37].
In contrast to our study, that indicated the analysis of water salinity revealed that the majority of samples (88.3%) fell within the valid or normal salinity range and the 86% fell in alkaline pH, indicating predominantly fresh water conditions across the study area, studies on correlation of Anopheles larval density with physicochemical characteristics in Benin and South east Iran indicated that there was high positive correlation of larval density with temperature [13], dissolved oxygen and salinity [12, 13].This difference might be because of the difference in the source of water in urban areas of Benin and South east Iran from our study area where almost 98% of source of water for larval habitats were ground water and chemical characters were in similar ranges.
In agreement with a study conducted by Mereta et al. 2013 [9], our study found that there was wide spread distribution of Anopheles larvae in small man-made aquatic habitats but in contrast to Mereta eta al. 2013 finding, in our study, Anopheles larva was abundant in both natural and artificial habitats.
Year- round presence of larvae in artificial habitats indicated larval density, distribution and abundance were not affected by season as far as water supply continued for the purpose of mainly for construction and brick making. Expansion of An.stephensi to rural sites could be due to the rapid urbanization and industrial expansion which is accompanied with construction of water storages. The An.stephensi samples found in rural sites (40km away from Dire Dawa city) in Aseliso was in plastic sheet habitat type implicated that where ever the artificial habitats are used they are suitable for the exotic species to breed in. This study indicated the strong association of some ecological settings with the density and abundance of An.stephensi in urban and peri-urban than rural sites this could be due to availability of enormous construction and brick manufacturing in urban and peri-urban settings.
Conclusion
An.stephensi was the predominant species found in the study area, followed by An.arabiensis, An.amharicu, An.pharoensis, An.coustani and An.pretoriensis. Uncovered water tankers were the most efficient habitat types in urban and peri-urban sites followed by plastic sheet habitat types while hoof prints and river basins under natural habitats were the most efficient habitat types in rural and urban sites respectively. An. stephensi was found in rivers of urban site and plastic sheets of rural sites co-existing with An.arabiensis in same sites. Both shared same ecological settings and habitat types. An.amharicus was found to breed in artificial habitats of the three ecological settings.
We report here the existence of An.amahricus for the first time in the study area, Eastern part of Ethiopia, expansion of An.stephensi in rural areas, breeding in natural habitats, and co-existing with An gambiae s.l .
Availability of brick making, shorter distance from living houses being in urban ecology were directly correlated with larval density. The result of our study revealed that there is significant correlation of ecological and seasonal variation in the density of Anopheles larvae. Habitat abundance and positivity of uncovered water tankers in urban and pre-urban sites could indicate for feasibility and proper implementation of larval source management.
Limitation of the study
The physicochemical tests were performed in a cross-sectional way but we believe to be done during main rain season to observe if seasonal variation affects the water chemistry.
A
Funding
This study obtained financial support from the National Institutes of Health (D43 TW001505, R01 A1050243 and U19 AI129326). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
A
Data Availability
The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.
A
A
Author Contribution
EA, MB, GY and DY conceived and designed the study; EA, HM, HA and TD were involved in data collection and field supervision. EA, TD and AA performed laboratory analysis. EA, ML, DZ and GZ performed the data analysis and EA drafted the manuscript. ML developed map of the study area and figures. TD, MB, DY and GY critically reviewed the manuscript. All the authors read and approved the final manuscript.
Ethics approval and consent to participate
Ethical clearance and amendment was obtained from the institutional review board (IRB) of the Institute of Health, Jimma University (Reference number) JUIH/IRB/236/25. Furthermore, permission was obtained from the Dire Dawa Health Bureau and health facilities.
Consent for publication
Not applicable
Competing interest
The authors declare that they have no competing interests
References
1.
WHO. Manual on practical entomology in Malaria. Part II. Geneva: World Health Organization; 1975.
2.
Wilson ML, Krogstad DJ, Arinaitwe E, Arevalo-Herrera M, Chery L, Ferreira MU, Ndiaye D, Mathanga DP, Eapen A. Urban malaria: understanding its epidemiology, ecology, and transmission across seven diverse ICEMR network sites. Am J Trop Med Hyg. 2015;93:110.
3.
Silver JB. Mosquito ecology: field sampling methods. springer science & business media; 2007.
4.
Githeko A, Mbogo C, Atieli F. Resting behaviour, ecology and genetics of malaria vectors in large scale agricultural areas of Western Kenya. Parassitologia. 1996;38:481–9.
5.
Jaenike J. Genetics of oviposition-site preference in Drosophila tripunctata. Heredity. 1987;59:363–9.
6.
Taylor R, Messenger LA, Abeku TA, Clarke SE, Yadav RS, Lines J. Invasive Anopheles stephensi in Africa: insights from Asia. Trends Parasitol. 2024;40:731–43.
7.
Balkew M, Mumba P, Yohannes G, Abiy E, Getachew D, Yared S, Worku A, Gebresilassie A, Tadesse FG, Gadisa E. An update on the distribution, bionomics, and insecticide susceptibility of Anopheles stephensi in Ethiopia, 2018–2020. Malar J. 2021;20:263.
8.
Service M. Medical entomology for students. 2000.
9.
Mereta ST, Yewhalaw D, Boets P, Ahmed A, Duchateau L, Speybroeck N, Vanwambeke SO, Legesse W, De Meester L, Goethals PL. Physico-chemical and biological characterization of anopheline mosquito larval habitats (Diptera: Culicidae): implications for malaria control. Parasit Vectors. 2013;6:320.
10.
Tsegaye A, Demissew A, Hawaria D, Abossie A, Getachew H, Habtamu K, Degefa T, Wang X, Lee M-C, Zhou G. Anopheles larval habitats seasonality and environmental factors affecting larval abundance and distribution in Arjo-Didessa sugar cane plantation, Ethiopia. Malar J. 2023;22:350.
11.
Thomas S, Ravishankaran S, Justin JA, Asokan A, Mathai MT, Valecha N, Thomas MB, Eapen A. Overhead tank is the potential breeding habitat of Anopheles stephensi in an urban transmission setting of Chennai, India. Malar J. 2016;15:274.
12.
Soleimani-Ahmadi M, Vatandoost H, Zare M. Characterization of larval habitats for anopheline mosquitoes in a malarious area under elimination program in the southeast of Iran. Asian Pac J Trop Biomed. 2014;4:S73–80.
13.
Hessou-Djossou D, Djègbè I, Ahadji-Dabla KM, Nonfodji OM, Tchigossou G, Djouaka R, Cornelie S, Djogbenou L, Akogbeto M, Chandre F. Diversity of larval habitats of Anopheles mosquitoes in urban areas of Benin and influence of their physicochemical and bacteriological characteristics on larval density. Parasit vectors. 2022;15:207.
14.
Keiser J, Singer BH, Utzinger J. Reducing the burden of malaria in different eco-epidemiological settings with environmental management: a systematic review. Lancet Infect Dis. 2005;5:695–708.
15.
Fillinger U, Kannady K, William G, Vanek MJ, Dongus S, Nyika D, Geissbühler Y, Chaki PP, Govella NJ, Mathenge EM. A tool box for operational mosquito larval control: preliminary results and early lessons from the Urban Malaria Control Programme in Dar es Salaam, Tanzania. Malar J. 2008;7:20.
16.
Ethiopia FMoH. Ethiopia malaria elimination strategic plan: 2021–2025. 2020.
17.
Teshome A, Erko B, Golassa L, Yohannes G, Irish SR, Zohdy S, Dugassa S. Laboratory-based efficacy evaluation of Bacillus thuringiensis var. israelensis and temephos larvicides against larvae of Anopheles stephensi in Ethiopia. Malar J. 2023;22:48.
18.
Yewhalaw D, Erena E, Degefa T, Kifle YG, Zemene E, Simma EA. Biological and residual activity of candidate larvicide formulation, SumiLarv 2MR, against an exotic invasive mosquito Anopheles stephensi Liston, 1901 (Diptera: Culicidae) in Ethiopia. Scie Rep. 2025;15:291.
19.
Yared S, Gebresilassie A, Aklilu E, Abdulahi E, Kirstein OD, Gonzalez-Olvera G, Che-Mendoza A, Bibiano-Marin W, Waymire E, Lines J. Building the vector in: construction practices and the invasion and persistence of Anopheles stephensi in Jigjiga, Ethiopia. Lancet Planet Health. 2023;7:e999–1005.
20.
Mutero CM, Okoyo C, Girma M, Mwangangi J, Kibe L, Ng’ang’a P, Kussa D, Diiro G, Affognon H, Mbogo CM. Evaluating the impact of larviciding with Bti and community education and mobilization as supplementary integrated vector management interventions for malaria control in Kenya and Ethiopia. Malar J. 2020;19:390.
21.
Tia J-PB, Tchicaya ES, Zahouli JZ, Ouattara AF, Vavassori L, Assamoi J-B, Small G, Koudou BG. Combined use of long-lasting insecticidal nets and Bacillus thuringiensis israelensis larviciding, a promising integrated approach against malaria transmission in northern Côte d'Ivoire. Malar J. 2024;23:168.
22.
Ayana GM, Jalilian A, Ashine T, Molla E, Hailemeskel E, Yemane DH, Yirgu H, Negash N, Teferi N, Teshome D. Larval source management in Ethiopia: modelling to assess its effectiveness in curbing malaria surge in dire Dawa and Batu Towns. Malar J. 2024;23:366.
23.
Esri. ArcGIS Pro Version 3.5 [Computer software].Redlands,CA:Esri. (2025). 2025.
24.
2023. TPVP: Ethiopia Final Entomological Report, April 2022-March 2023. 2023.
25.
Emiru T, Getachew D, Murphy M, Sedda L, Ejigu LA, Bulto MG, Byrne I, Demisse M, Abdo M, Chali W. Evidence for a role of Anopheles stephensi in the spread of drug-and diagnosis-resistant malaria in Africa. Nat Med. 2023;29:3203–11.
26.
Merga H, Degefa T, Birhanu Z, Abiy E, Lee M-C, Yan G, Yewhalaw D. Urban malaria and its determinants in Eastern Ethiopia: the role of Anopheles stephensi and urbanization. Malar J. 2024;23:303.
27.
WHO. Manual on practical entomology in Malaria. Part II. Methods and techniques. 1975.
28.
Organization WH. Malaria entomology and vector control. In Malaria entomology and vector control. 2013.
29.
Coetzee M. Key to the females of Afrotropical Anopheles mosquitoes (Diptera: Culicidae). Malar J. 2020;19:1–20.
30.
Glick JI. Illustrated key to the female Anopheles of southwestern Asia and Egypt (Diptera: Culicidae). Mosq Syst. 1992;24:125–53.
31.
Musapa M, Kumwenda T, Mkulama M, Chishimba S, Norris DE, Thuma PE, Mharakurwa S. A simple Chelex protocol for DNA extraction from Anopheles spp. J Vis Exp 2013:3281.
32.
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–9.
33.
Coetzee M, Hunt RH, Wilkerson R, Della Torre A, Coulibaly MB, Besansky NJ. Anopheles coluzzii and Anopheles amharicus, new members of the Anopheles gambiae complex. Zootaxa. 2013;3619:246–74.
34.
Demissew A, Hawaria D, Kibret S, Animut A, Tsegaye A, Lee M-C, Yan G, Yewhalaw D. Impact of sugarcane irrigation on malaria vector Anopheles mosquito fauna, abundance and seasonality in Arjo-Didessa, Ethiopia. Malar J. 2020;19:344.
35.
Silver JB. Mosquito ecology: field sampling methods. Springer; 2008.
36.
Demissew A, Animut A, Kibret S, Tsegaye A, Hawaria D, Degefa T, Getachew H, Lee M-C, Yan G, Yewhalaw D. Evidence of pyrethroid resistance in Anopheles amharicus and Anopheles arabiensis from Arjo-Didessa irrigation scheme, Ethiopia. PLoS ONE. 2022;17:e0261713.
37.
Nikookar SH, Fazeli-Dinan M, Azari-Hamidian S, Mousavinasab SN, Aarabi M, Ziapour SP, Esfandyari Y, Enayati A. Correlation between mosquito larval density and their habitat physicochemical characteristics in Mazandaran Province, northern Iran. PLoS Negl Trop Dis. 2017;11:e0005835.
Total words in MS: 5218
Total words in Title: 21
Total words in Abstract: 314
Total Keyword count: 5
Total Images in MS: 7
Total Tables in MS: 5
Total Reference count: 37