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).
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.
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).
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).
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).
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).
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.