Habitat characteristics and spatial distribution of Anopheles mosquito larvae in malaria elimination settings in Dembiya District, Northwestern Ethiopia

doi:10.21203/rs.3.rs-830514/v1

Download PDF

Research Article

Habitat characteristics and spatial distribution of Anopheles mosquito larvae in malaria elimination settings in Dembiya District, Northwestern Ethiopia

https://doi.org/10.21203/rs.3.rs-830514/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

The persistence and productivity of larval habitat is a major factor that regulates adult Anopheles mosquitoes’ density and malaria transmission intensity. A study on Anopheles mosquitoes breeding habitat diversity, distribution, characteristics and larval density in different seasons across various habitat types is important to design effective larval control strategies. This study was aimed to investigate the Anopheles species composition, the productivity of larval breeding habitats, and their spatial distribution in selected localities of Dembiya District. A longitudinal study on characteristics and productivity of Anopheles larval breeding habitats was conducted from June 2018 to May 2019 in selected localities of Dembiya District. Anopheles larvae were collected using standard WHO dipper (350 ml capacity) and droppers depending on the size of the breeding habitats. Physicochemical characteristics of breeding habitats were measured and Anopheles mosquitoes were identified by using morphological keys and polymerase chain reaction (PCR). Logistic regression was used to assess the association of environmental factors with the presence or absence of Anopheles mosquitoes larvae. A total of 1,629 Anopheles larvae and 185 pupae were collected from both localities. These comprise Anopheles arabiensis, An. pharoensis, An. coustani, An. christyi, An. squamosus, An. demeilloni, An. danicalicus, and An. cinereus. The highest density of Anopheles larvae was collected from at a drying water canal (14.7 ± 3.5 larvae/dip) and the lowest larval density was recorded in rain water pools (0.2 ± 0.2 larvae/dip). The presence or absence of Anopheles larvae were significantly associated with physical characteristics of the breeding habitats such as turbidity (mid turbid) (AOR = 66.03; 95% CI: 2.01-2168.24, p = 0.019) and presence of grasses (AOR = 12.62; 95% CI: 1.29-122.78, p = 0.029). This study indicated that breeding sites persist and support Anopheles mosquito breeding activities and impact malaria control and elimination programs. Incorporating vector control strategies targeting Anopheles larvae as a part of malaria intervention strategies could enhance the malaria control and elimination program in the study area.

Anopheles larvae

breeding habitat

An. arabiensis

An. pharoensis

drainage canals

Dembiya

Mosquitoes are responsible for the transmission of different types of diseases such as malaria, yellow fever, dengue, and West Nile fever (Alonso et al. 2011; Braack et al. 2018). Malaria is a life-threatening infectious disease caused by protozoan parasites of the genus Plasmodium that are transmitted through the bites of infected female Anopheles mosquitoes. Globally, an estimate of 241 million malaria cases and 627,000 deaths were reported in 2020, of which 95% of malaria cases and highest proportion of deaths were reported from a WHO Africa region (WHO, 2021). In this region 80% of all malaria deaths are recorded in children’s under the age of 5. From the total malaria case across the globe 55% is recorded in sub-Sahara Africa (WHO, 2021). In Ethiopia, more than 60% of the population is at risk of malaria infection and 68% of the country’s landmass is malarious (FMoH, 2014). The diverse ecology and favorable environmental conditions of the country supported the rapid development of Anopheles mosquitoes and Plasmodium parasites (Taffese et al. 2018).

There are more than 400 different species of Anopheles mosquitoes worldwide, of which around 30 are known malaria vectors (WHO, 2018). In Ethiopia, An. arabiensis is the primary malaria vector, whereas An. funestus, An. pharoensis and An. nili are secondary vectors (Gillies and Coetzee, 1987). The most important stage of the Anopheles mosquito life cycle such as egg-laying, larval and pupal development, and adult emergence takes place in the aquatic environment (Oyewole et al. 2009). Different species of Anopheles mosquitoes have their own preferred aquatic habitats. For instance, An. arabiensis prefers to breed in temporary, small, sunlit, clear, and shallow fresh-water pools (Gimnig et al. 2001; Edillo et al. 2006; Himeidan and Rayah, 2008).

The abundance and distribution of potential breeding habitats of Anopheles mosquitoes determine the density of adult Anopheles mosquitoes and transmission of malaria parasites (Carter et al. 2000; Rejmánková et al. 2013). Several environmental characteristics like climate, physical, chemical, and biological conditions of the breeding habitats have effects on the development and survival of the Anopheles mosquito larvae (Mereta et al. 2013; Roux and Robert 2019). Anthropogenic factors like agricultural expansion, construction of dams and urbanization affects the diversity and distribution of Anopheles mosquito breeding habitat and larval development (De Silva and Marshall, 2012). Climate change associated with these anthropogenic effects favors rapid development of Anopheles mosquitoes and Plasmodium parasite in areas with antecedently low malaria transmission (Alonso et al. 2011).

In Ethiopia, malaria control program mainly depends on the management of clinical malaria cases or control efforts targeting adult mosquitoes by selective indoor residual spray (IRS) and insecticide-impregnated bed nets (LLINs) (FMoH, 2014). However, these strategies are challenged by the development of drug resistant Plasmodium parasites and insecticide resistant vector species (Messenger et al. 2017; Taffese et al. 2018; Loha et al. 2019). Considering this, interests are rising to use larval source management as a part of integrated vector management in the country (Gari and Lindtjørn, 2018; Asale et al. 2019). Nevertheless, our knowledge about mosquitoes breeding habitat diversity, distribution and characteristics is limited and insufficient to design effective malaria intervention strategies by larval control (Rejmánková et al., 2013).

Therefore, the aim of this study was to assess species composition, distribution, and ecology of Anopheles mosquito larvae and pupae in malaria endemic localities of Dembiya District, Northwestern Ethiopia. The result of this study will provide important information to design effective vector control strategies by larval source management.

Description of the study area

A longitudinal study on species composition, breeding habitat characteristics, and spatial distribution of Anopheles mosquito larvae and pupae was conducted from June 2018 to May 2019 in two localities (Guramba Bata and Arebiya) of Dembiya district in North Gondar Administrative Zone of Amhara Regional State, Northwestern Ethiopia (Fig. 1). The district is located at 12°39'59.99" N and 37°09'60.00" E. Kola Diba is the administrative center of the district, located 750 km north of Addis Ababa and 35 km southwest of Gondar. The district is bordered by Lake Tana in the south. Dembiya district has 45 localities with an estimated population of approximately 271,000. The majority of the population (91%) lives in rural areas, with most engaged in farming activities; the remaining 9% live in urban areas. The district has 49,528 rural households with 4.3 mean household sizes (CSA, 2007).

The elevation of Dembiya District ranges from 1500 m to 2600 m a.s.l. The agro-ecology of the District is mid-altitude (Woynadega) with a mean annual minimum temperature of 11 ⁰C and a maximum of 32 ⁰C, respectively. The mean annual rainfall ranges from 995 mm to 1175 mm. Land use data from the district agricultural bureau indicated that most of the land is considered as plain (87%), while mountainous, valleys, and wetland make up 5%, 4.8%, and 3.2%, respectively. Out of the total area of the District, 31% is cultivated land, 16% is none cultivable land, 5.6% forest and bush, 12.8% grazing, 8.1% is covered with water, 20.2% swamp and 4.3% is a residential area. The district receives bimodal rainfall, with the short rainy season from March to May and the main rainy season from June to September.

The major crops grown by small-holder farmers include teff (Eragrostis tef), maize (Zea mays), barley (Hordeum vulgare), sorghum (Sorghum bicolor), and finger millet (Eleusine coracana). Besides, legumes and pulses such as chickpeas (Cicer arietinum) and cowpeas (Vigna unguiculata) are also grown in the district. Farmers also grow pulse and spice crops like pepper (genus Capsicum), niger seed (Guizotia abyssinica), fenugreek (Trigonella foenum-graecum), black cumin (Nigella sativa), White cumin (Cuminum cyminum), and rice (Oryza sativa) on a limited amount of farmlands.

One of the study localities, Guramba Bata (12⁰21’57.75’’N and 37⁰20'25.31'' E, altitude 1,795–1,820 m.a.s.l.), has a seasonal river “Ahya gedel” or “Nededo” which forms intermittent mosquito breeding water bodies until the end of December. Guramba Bata has one health post and one health center, 1113 households with 6008 inhabitants (2974 males and 3034 females) in 2017/18 (District Health Office report) (Fig. 1).

The second study locality, Arebiya (12⁰20'26.59''N and 37⁰22'16.04'' E) has “Megech” river serves as a water source during a dry season and flows into Lake Tana. This locality has 1976 households and a total of 8632 inhabitants (4298 males and 4384 females) in 2017/18. Arebiya has only one health post (District Health Office report) (Fig. 1).

Collection of Anopheles mosquitoes larvae and pupae

The larvae and pupae of Anopheles mosquitoes were collected from different breeding sites such as drainage canals, artificial pit shelters, hoof prints, and rain pools (Fig. 2). The larvae collection was performed for two consecutive days per month for 12 months. Before larval sampling, habitats were inspected for the presence or absence of Anopheles larvae and pupae. When Anopheles larvae were present, repeated dips were taken using standard WHO dippers (350 ml capacity) depending on the size of each larval habitat (WHO 1975). For small habitats like hoof prints, samples taken with droppers from several sites were pooled to get appropriate larval sampling (Soleimani-Ahmadi et al. 2014). Sampling was performed in the morning (09:00–12:00) or in the afternoon (14:00–17:00) for about 30 minutes at each larval breeding habitat.

Species Identification

Anopheles mosquito larvae and pupae were separated from Culicines using a hand lens (10X) based on gross morphological characteristics (Verrone 1962; Gillies and Coetzee, 1987). The 3rd and 4th instar Anopheles larvae, and Anopheles pupae collected from each type of breeding habitat, were transferred to separately labeled mosquito breeders (BioQuip, Dimensions: 7 − 3/4″ (195 mm) high x 3–5/8″ (92 mm) diameter) and reared in the field laboratory based on WHO guidelines (WHO, 1975). The adult Anopheles that emerged from field collected larvae and pupae were used for species identification based on morphological characteristics (Verrone 1962; Gillies and Coetzee 1987). Morphologically identified adult An. gambiae sensu lato was individually preserved in Eppendorf tubes with silica gel and cotton for further molecular identification.

Species identification using rDNA–polymerase chain reaction (PCR)

Preserved adult An. gambiae s.l specimens from field collected and reared larvae and pupae were further identified to sibling species using a ribosomal DNA polymerase chain reaction (PCR) by including the primers for An. gambiae s.s, An. arabiensis, An. quadriannulatus and An. amharicus (Scott et al. 1993).

Physico-chemical characteristics of larval habitats

Physical characteristics of the breeding habitats (positive or negative for Anopheles larvae or pupae), including habitat type, water depth, turbidity, vegetation, presence of algae, bottom substrate, habitat stability, lotic or lentic water, sunlight intensity, and distance from nearby houses were measured and recorded (Minakawa et al. 1999).

Chemical characteristics of the breeding habitats (positive or negative for Anopheles larvae or pupae) such as temperature, pH, and conductivity were measured at the field with HANNA© HI 98130 Combo pH & EC tester (Hanna Instruments Inc., Kehl am Rhein, Germany) with the probes placed 2 to 3 cm below the water surface.

Data analysis

The density of Anopheles larvae and pupae was expressed as the total number of Anopheles larvae per total number of dips taken. After checking for normality, all dependent variables were log₁₀(x + 1) transformed, and subjected to statistical analysis. Since the data were found to be normally distributed after transformation, parametric testes such as one-way analysis of variance (ANOVA) and student t-test (for independent variables with two categories) were used to analyze differences in mean larval densities among breeding habitat types and other environmental variables. When significant differences were observed in one-way ANOVA, means were separated using Tukey’s HSD (Tukey’s Honestly Significant Difference) test at α = 0.05. Multiple logistic regression analysis was used to detect the best predictor environmental variables associated with the presence or absence of Anopheles mosquito larvae (Sattler et al. 2005). The percentages of species composition of Anopheles mosquitoes collected from each breeding habitat was calculated (number of Anopheles mosquitoes species * 100/ total number of species identified). Pearson correlation analysis was used to assess the relationship between larval densities and chemical characteristics such as pH, temperature, and conductivity. The data were analyzed using SPSS version 20 (Armonk, NY: IBM Corp), p ≤ 0.05 were considered as significant.

Breeding habitat types and abundance of Anopheles larvae and pupae

During the one-year study period, a total of 108 potential larval habitats (60.2% from Arebiya and 39.8% from Guramba Bata) were assessed for the presence of Anopheles larvae. From the total sampled potential habitats only 41 were positive for Anopheles larvae and pupae (Table 1). From which, the predominantly encountered larval habitats were rain pools (17.6% (n = 19)), and river pool (17.6% (n = 19)). More than 80% of the larval breeding habitats were recorded during the long rainy season. Whereas during the dry season, the distribution of Anopheles larvae and pupae were restricted to water pools at riversides, pits dug for plastering a house, and temporary habitats around the hand pump water well.

The highest mean densities of Anopheles larvae were recorded from breeding sites in drainage canals (14.7 ± 3.5 larvae/dip), followed by abandoned burrow pits dug for plastering houses (8.8 ± 3.1 larvae/dip), swamps (3.8 ± 1.2 larvae/dip), and hoof prints (3 ± 1.2 larvae/dip). In addition, Anopheles larvae were also recorded from puddles (2.7 ± 2.7 larvae/dip), water pools at riversides (2.0 ± 0.9 larvae/dip), and tire tracks (0.4 ± 0.4 larvae/dip) (Table 1). The difference in mean densities of Anopheles larvae and pupae among habitats was statistically significant (F_{8, 99}= 9.85, p = 0.000) and (F_{8, 99}= 3.46, p = 0.001) respectively (Table 1).

Table 1

Distribution of *Anopheles* mosquito larvae and pupae in different breeding habitats in the two localities of Dembiya District, Northwestern Ethiopia.
Habitat type	No. of breeding sites surveyed (%)	No. of dips	No. of larvae collected	Mean (larvae/dips) ± se	No. of pupae collected	Mean (pupae/dips) ± se
Riversides	19 (17.6 )	99	279	2.0 ± 0.9^ab	48	0.3 ± 0.1^ab
Burrow pits	14 (12.96)	70	421	8.8 ± 3.1^bc	49	0.99 ± 0.4^ab
Drainage canals	12 (11)	48	566	14.7 ± 3.5^c	59	1.6 ± 0.6^b
Tire tracks	8 (7.4)	24	9	0.4 ± 0.4^a	3	0.1 ± 0.1^a
Hoof prints	7 (6.5)	19	63	3.0 ± 1.2^ab	5	0.2 ± 0.2^ab
Swamps	14 (12.96)	56	205	3.8 ± 1.2^ab	9	0.2 ± 0.1^ab
Rain pools	19 (17.6 )	69	30	0.2 ± 0.2^a	0	0.0 ± 0.0^a
Puddles Streams	7 (6.5) 8 (7.4)	33 30	56 0	2.7 ± 2.7^a 0.0 ± 0.0^a	12 0	0.6 ± 0.6^ab 0.0 ± 0.0^a
Total	108(100%)	448	1629		185

The mean larval and pupal densities with different letter designations in a column are significantly different with Tukey’s HSD post hoc analysis at α = 0.05.

Species composition and monthly distribution of Anopheles larvae and pupae

The species composition of Anopheles mosquitoes identified from the two study sites during the study period is presented in Table 2. A total of 1,629 Anopheles larvae and 185 pupae were collected from the two localities. From the total collected immature stages of Anopheles mosquitoes, 52.3% (852) larvae and 65.9 (122) pupae were from Arebiya and 47.7% (777) larvae and 34.1% (63) pupae were from Guramba Bata. The difference in mean Anopheles larval and pupal density between the two study sites were not statistically significant (t (106) = -0.454, p = 0.651) and (t (106) = 0.70, p = 0.485) respectively.

From the total collected Anopheles larvae and pupae, 835 females and 788 males have successfully emerged into adults. The rest 191 larvae and pupae were not able to emerge to adult. Therefore, only 835 female Anopheles mosquitoes were subjected to species identification based on morphological features. All An. gambiae s.l samples used for species identification using PCR were found to be An. arabiensis. Eight species of Anopheles such as An. arabiensis, An. pharoensis, An. coustani, An. christyi, An. squamosus, An. demeilloni, An. danicalicus, and An. cinereus were identified (Table 2). From the total identified Anopheles species, An. arabiensis (59.2%) was dominant followed by An. pharoensis (35.3%) and An. coustani (2.99%). The least common species were An. dancalicus and An. cinereus (Table 2).

The monthly distribution of Anopheles larvae and pupae showed that the highest density of Anopheles larvae in Arebiya was collected in June, September, October, and May (Fig. 3a). In the meantime, less density of Anopheles larvae was recorded during July and August, which is corresponding to pick monthly rainfall, which creates unstable larval breeding habitat due to over flooding. Similarly, in Guramba Bata high larval density was recorded starting from June, August, September, and October (Fig. 3b). However, the density of Anopheles larval density sharply declines in both study areas during a dry season (January, February, and March) (Fig. 3a & b).

Table 2

Species composition and abundance of *Anopheles* mosquitoes collected from the two localities of Dembiya District, Northwestern Ethiopia (Values in parenthesis are percentages).
Species	Arebiya	Guramba Bata	Total
An. arabiensis	297	197	494 (59.2)
An. phareonsis	113	182	295 (35.3)
An. christyi	6	2	8 (1)
An. squamosus	3	0	3 (0.4)
An. coustani	9	16	25 (2.99)
An. demeilloni	2	4	6 (0.7)
An. dancalicus	0	1	1 (0.1)
An. cinereus	0	3	3 (0.4)
Total	430 (51.5)	405 (48.5)	835 (100)

Species specific monthly distribution of Anopheles mosquitoes

The species-specific monthly distribution of the dominant Anopheles mosquitoes is indicated in Fig. 4a & b. The highest density of An. arabiensis in Arebiya was recorded during June (16.5%), September (9.1%), and May (21.6%) (Fig. 4a). Similarly, in Guramba Bata the highest density of An. arabiensis was recorded during June (12.8%), July (10.4%), and September (12.3%), but its number sharply declined after the end of the long rainy season (Fig. 4b). The number of An. pharoensis reached its peak around the end of the main rainy season in the two study areas (Fig. 4a & b).

Anopheles species specific breeding habitat types

The associations of species-specific Anopheles larvae in different breeding habitat types are presented in Table 3. An. arabiensis larvae were collected from drainage canals, burrow pits, water pools at riversides, tire tracks, hoof prints, and puddles (Table 3). A high number of An. arabiensis was collected near the edge of a small temporary and permanent habitat with still and mid turbid water, grass, and had full sunlight access. An. pharoensis was collected from a wide range of permanent habitats such as water pools at riversides, swamps, and grassy burrow pits (Table 3). The breeding habitats of An. pharoensis were relatively turbid, full of vegetation, and had partial sunlight access. An. coustani was collected from permanent habitats such as riverside water pools, burrow pits, and swamps (Table 3). The most common breeding sites of this species were usually permanent habitats with relatively turbid water, vegetation, and partial sunlight access.

Table 3

Species-specific spatial distribution of *Anopheles* mosquitoes in the two localities of Dembiya District, Northwestern Ethiopia (Values in parentheses are the percentages).
	Anopheles mosquito larval breeding habitat types
Species	Riverside	Burrow pit	Drainage canal	Tire track	Hoof print	Swamp	Rain pool	Puddles	Total
An. arabiensis	78 (15.8)	148 (29.9)	174 (35.2)	5 (1)	26 (5.3)	26 (5.3)	11 (2.2)	26 (5.3)	494 (59.2)
An. pharoensis	83 (28.1)	69 (23.4)	64 (21.7)	-	-	79 (26.8)	-	-	295 (35.3)
An. coustani	9 (36)	5 (20)	4 (16)	-	-	7 (28)	-	-	25 (2.9)
An. christyi	1 (12.5)	5 (62.5)	1 (12.5)	-	-	1 (12.5)	-	-	8 (0.96)
An. demeilloni	-	2 (33.3)	4 (66.7)	-	-	-	-	-	6 (0.7)
An. squamosus	-	2 (66.7)	-	-	-	1 (33.3)	-	-	3 (0.4)
An. cinereus	3 (100)	-	-	-	-	-	-	-	3 (0.4)
An. dancalicus	-	-	1 (100)	-	-	-	-	-	1 (0.1)
Total	174 (20.8)	231 (27.7)	248 (29.7)	5 (0.6)	26 (3.1)	114 (13.7)	11 (1.3)	26 (3.1)	835 (100)

Association of Anopheles larval densities with physical characteristics of larval breeding habitat types

The associations of mean larval densities with physical variables of the breeding habitat types are presented in Table 4. The densities of Anopheles larvae were significantly associated with shallow depth (≤ 0.5 m), medium turbidity, availability of grasses, muddy bottom substrates, distance to the nearest house (≤ 100 m), and presence of algae (p ≤ 0.05) (Table 4).

Table 4

Association of physical characteristics of breeding habitats with *Anopheles* mosquito larval density, Dembiya District, Northwestern Ethiopia.
Physical factors	Variables
		Densities of Anopheles larvae
		Mean ± se	p	F
Average depth	≤ 0.5 m	5.2 ± 1.01	0.040	2.144
Average depth	≥ 0.5 m	1.8 ± 0.57
Turbidity	Low	0.63 ± 0.45	0.000	19.899
	Med	6.96 ± 1.24
	High	0.55 ± 0.37
Vegetation	No vegetation	0.69 ± 0.43	0.000	-4.89
Vegetation	Grass	7.44 ± 1.29
Bottom substrate	Muddy	6.74 ± 1.08	0.000	8.231
	Clay	0.44 ± 0.34
	Sandy	3.21 ± 3.21
Water current	Stagnant	4.69 ± 0.96	0.180	-1.35
Water current	Flowing	2.34 ± 0.92
Light intensity	Full sunlight	4.28 ± 0.81	0.687	0.404
Light intensity	Shaded	3.57 ± 1.84
Distance	≤ 100m	7.04 ± 1.22	0.000	7.76
	100 -200m	0.53 ± 0.53
	201-300m	0.00 ± 0.00
	301–400m	0.17 ± 0.17
Algae	Present	10.16 ± 2.16	0.000	5.66
Algae	Absent	1.74 ± 0.40
Water persistency	Temporary	2.6 ± 0.70	0.99	1.66
Water persistency	Permanent	5.15 ± 1.18

A bivariate analysis showed that average depth (COR = 2.54; 95% CI: 0.98–6.6, p ≤ 0.05) turbidity (COR = 28.4; 95% CI: 3.61-224.27, p ≤ 0.05), vegetation (COR = 17.82; 95% CI: 5.62–56.54, p ≤ 0.05), distance to the nearest house (COR = 32.7; 95% CI: 4.15-257.27, p ≤ 0.05) and presence of algae (COR = 4.92; 95% CI: 1.99–12.11, p ≤ 0.05) were significantly associated with the presence or absence of Anopheles larvae (Table 5).

Table 5

A bivariate and multivariate analysis of association of physical factors of breeding habitats with presence or absence of *Anopheles* mosquito larvae, Dembia District, Northwestern Ethiopia.
Factors	Variables	OR (95% CI)		p-value
		COR	AOR
Average depth	≤ 0.5	2.5 (0.98–6.60)*	0.4 (0.04–4.32)	0.450
Average depth	≥ 0.5	1	1
Turbidity	Low	3.6 (0.35–37.36)	22.5 (0.66-774.51)	0.084
	Med	28.4 (3.61-224.27)*	66 (2.01-2168.24)	0.019
	High	1	1
Vegetation	No vegetation	1	1
Vegetation	Grass	17.8 (5.62–56.54) *	12.6 (1.29-122.78)	0.029
Bottom substrate	Muddy	21.7 (2.64-177.76)	8.6 (0.24-308.82)	0.238
	Clay	0.00 (0.00)	0.00 (0.00)	0.997
	Sandy	1	1
Water current	Stagnant	0.36 (0.12–1.05)	0.1 (0.01–1.29)	0.077
Water current	Flowing	1	1
Light intensity	Full sunlight	2.1 (0.76–5.79)	2.21 (0.25–19.92)	0.479
Light intensity	Shaded	1	1
Distance	≤ 100m	32.7 (4.15-257.27)*	27.99 (0.65-1215.17)	0.083
	100 -200m	1.39 (0.08–23.71)	0.99 (0.01-132.14)	1.00
	201-300m	0.00	0.00
	301–400m	1
Algae	Present	4.92 (1.99–12.11)*	6.94 (0.78–62.09)	0.083
Algae	Absent	1
Water persistency	Temporary	1.15 (0.51–2.61)	1.51 (0.23–9.99)	0.67
Water persistency	Permanent	1
Surface debris	Low	1.81 (0.54–6.08)	7.90 (0.53-117.17)	0.133
	Medium	1.11 (0.34–3.68)	0.723 (0.07–7.89)	0.791
	High	1	1
* indicates statistically significant values at p = 0.05

The final model for the parameters associated with presence or absence of Anopheles larvae showed that turbidity (AOR = 66.03; 95% CI: 2.01-2168.24, p = 0.019) and vegetation cover (AOR = 12.62; 95% CI: 1.29-122.78, p = 0.029) were key physical factors which determine the densities of Anopheles larvae (Table 5).

Correlation of chemical characteristics of larval habitat with Anopheles density

The total density of Anopheles larvae was positively correlated with temperature (r = 0.331 and p = 0.013) and pH (r = 0.697 and p = 0.00). However, the density Anopheles larvae was negatively correlated with conductivity (r= -0.321, p = 0.016) (Table 6).

Table 6

Correlation of chemical characteristics of breeding habitat with *Anopheles* mosquito larval density in Dembiya District, Northwestern Ethiopia.
Chemical factors	Larval density
Temperature (°C)	0.331^*
pH	0.697^**
Conductivity (µS/m)	-0.321^*
**. Correlation is significant at the 0.01 level
*. Correlation is significant at the 0.05 level

This study showed that the spatiotemporal distribution and species composition of Anopheles mosquito larval density were greatly affected by the physicochemical characteristics of the breeding habitat and the rainfall pattern of the localities. Describing larval habitat characteristics in terms of environmental attributes and identifying relationships between breeding habitat and larval density is important to develop novel methods of vector control by targeting the aquatic stage of Anopheles mosquitoes in areas with high vector intervention strategies.

During this study rain pools, river pools, burrow pit, swamp, drainage canal, tire track, hoof print, puddle, and stream were the encountered larval breeding habitats, of which rain pools and river pools were dominantly observed breeding habitats during the rainy season. Rain during the rainy season produces many rain pools and river edges, which are potential sites for larval development. In concurrent with this study, Anopheles mosquitoes prefer to 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 (Shililu et al. 2003; Yohannes et al. 2005; Omlin et al. 2007). In addition, similar habitat types were recorded in previous studies elsewhere in Ethiopia (Mereta et al. 2013) and Kenya (Imbahale et al. 2011).

The distribution of Anopheles larvae during the dry season was limited to water pools at riversides, pits dug for plastering a house, drainage canals, and water pools around hand pump water well. This is because during the dry season the water will be confined to temporary habitats such as riversides, burrow pits, drainage canals, and swamps which are important for the reproduction of Anopheles mosquitoes. Similarly, studies showed that the number and size of Anopheles breeding habitat are reduced during the dry season (Animut and Negash 2018). This restricted distribution of Anopheles mosquito breeding habitat during the dry season makes them more vulnerable for larval management. Larval management during a dry season is less costly than management during a wet season in areas where the dry season is associated with limited breeding habitats (WHO, 2013). Dry season larval management will hamper the exponential reproduction rate of Anopheles larvae at the end of the long rainy season hence it will limit malaria transmission (Animut and Negash 2018). This approach could be effective to reduce malaria transmission associated with insecticide resistant malaria vectors and outdoor host seeking mosquitoes, as it manages the immature stage (egg, larvae, and pupae) confined in a small aquatic environment (Killeen et al. 2002).

The result of this study indicated that the density of Anopheles larvae varies among habitat types, where a significantly highest density was recorded from water pools at drainage canal and the grassy edge of burrow pits. This larval density variation in the different habitats could be explained by the spatiotemporal differences in food resource and predation pressure in different habitats and the complex interaction between physicochemical factors such as water turbidity, depth, temperature, salinity, and dissolved oxygen (Mala and Irungu 2011; Kipyab et al. 2015; Roux and Robert 2019).

Anopheles mosquito species such as An. arabiensis, An. pharoensis, An. coustani, An. christyi, An. squamosus, An. demeilloni, An. danicalicus, and An. cinereus were identified from the two study areas. This result is in line with finding from the west Gojam zone (Animut and Negash 2018), south-central Ethiopia (Animut et al. 2012), and Central Ethiopia (Kenea et al. 2011). The highest density of Anopheles larvae was recorded during the rainy seasons. The reason is that rainfall and humidity strongly affect the availability of larval habitat, Anopheles species, and distribution (Imbahale et al. 2011).

An. arabiensis was the dominant species identified in the study area, similar to report from southwest Ethiopia (Getachew et al. 2020), Addis Zemen, South Gondar, Ethiopia (Kindu et al. 2018), South-central Ethiopia (Animut et al. 2012), and Western Kenya (Kweka et al. 2012). The highest density of An. arabiensis was recorded in June and September when the rainy season starts and retreats. Similar observations have shown that populations of An. arabiensis usually increase as the rains withdraw (Getachew et al. 2020). The result of this study indicated that An. arabiensis was dominantly distributed in small temporary habitats with still and less turbid water and grasses and full sunlight access. The reason for this could be the presence of less larval predation pressure in small temporary habitats than permanent habitats and temporary habitats with full sunlight access provide warmer water, which resulted in a high algal density (source of food for larvae) and rapid development of larvae to pupae (Gimnig et al. 2002). In concurrent with this study, An. arabiensis was reported to breed in small, temporary habitats with algae such as footprints, rain pools, puddles, tire tracks, and garden wells (Mattah et al. 2017). These species were also identified from temporary habitats such as swamps, irrigation canals, sand pools, canal leakage pools, water harvesting pools, and brick-making pits in central Ethiopia (Kenea et al. 2011).

Anopheles pharoensis was dominantly distributed in aquatic habitats with turbid, full of vegetation, and partial sunlight access. Different reports also supported this observation, claiming that the density of An. pharoensis was higher in aquatic habitats with floating vegetation and with shady conditions (Teklu et al. 2010).

Physical characteristics of the Anopheles mosquito breeding habitat such as low turbidity and the presence of grass were the limiting factors that determine the presence or absence of Anopheles mosquito larvae in this study. In addition, the density of Anopheles mosquito larvae was positively correlated with temperature and pH and negatively correlated with water conductivity. Conductivity measures the amount of inorganic matter and ions in water, therefore as turbidity of the water increases due to flooding conductivity also increases proportionally, which in turn affects the development of the larval population (Edillo et al. 2006). Low water turbidity and full sunlight access increase the water temperature and hence leads to a rapid larval development (Paaijmans et al. 2010). This result coincides with previous works conducted in the highlands of Ethiopia (Dejenie et al. 2011), Kenya (Minakawa et al. 1999), and Tanzania (Emidi et al. 2017).

Anopheles mosquitoes such as An. arabiensis, An. pharoensis, An. coustani, An. christyi, An. squamosus, An. demeilloni, An. danicalicus, and An. cinereus were identified during this study, of which An. arabiensis was the dominant species. Breeding habitats such as drainage canals and burrow pits served as a potential reproduction site of Anopheles mosquitoes in the study area. The distribution of Anopheles mosquito larvae is greatly affected by physicochemical characteristics of the breeding habitats such as water turbidity, vegetation cover, pH, temperature, conductivity, and season. A malaria control and elimination program in the study area should incorporate larval management strategies by using larvicides and source reduction as a part of IVM (community-based habitat modification during a dry season when the number of breeding habitat is limited). Further study on identification of the natural predators of Anopheles mosquitoes larvae/pupae in the study area is recommended.

Confidence Limit

Degree of freedom

COD

Crud odd ratio

AOD

Adjusted odd ratio

Standard error

PCR

Polymerase chain reaction

Ethics approval

Not applicable

Consent to participate

Not applicable

Consent for publication

Not applicable

Availability of data and materials

The data sets supporting the conclusions of this article are provided in the manuscript.

Conflicts of interest/Competing interests

The authors declare that there is no conflict of interest

Funding

This study was financed by Addis Ababa University.

Authors’ contributions

MT, HT, YW and SD designed the study. HT, YW and SD supervised and MT conducted the field work, rearing and identification experiments. MT conducted the statistical analyses. MT developed first draft, HT, YW and SD revised the manuscript. All authors read and approved the final manuscript.

Acknowledgements

The authors want to acknowledge Woldia University for providing a study leave and Addis Ababa University for providing financial support from Thematic Research Fund. We also express our appreciation to Aklilu lemma institute of pathobiology, Addis Ababa University, staff members for their technical assistance.

Alonso D, Bouma MJ, Pascual M (2011). Epidemic malaria and warmer temperatures in recent decades in an East African highland. Proc. Royal Soc 7;278(1712):1661-9. https://doi.org/10.1098/rspb.2010.2020.

Animut A, Gebre-Michael T, Balkew M, Lindtjørn B (2012). Abundance and dynamics of anopheline larvae in a highland malarious area of south-central Ethiopia. Parasit Vectors 1;5(1):117. https://doi.org/10.1186/1756-3305-5-117.

Animut A, Negash Y (2018). Dry season occurrence of Anopheles mosquitoes and implications in Jabi Tehnan District, West Gojjam Zone, Ethiopia. Malar. J 17(1):1-7. https://doi.org/10.1186/s12936-018-2599-4.

Asale A, Kussa D, Girma M, Mbogo C, Mutero CM (2019). Community based integrated vector management for malaria control: lessons from three years’ experience (2016–2018) in Botor-Tolay district, southwestern Ethiopia. BMC Public Health 1;19(1):1318. https://doi.org/10.1186/s12889-019-7606-3.

Braack L, de Almeida AP, Cornel AJ, Swanepoel R, De Jager C (2018). Mosquito-borne arboviruses of African origin: review of key viruses and vectors. Parasit Vectors 11(1):29. https://doi.org/10.1186/s13071-017-2559-9.

Carter R, Mendis KN, Robert SD (2000). Spatial targeting of interventions against malaria. Bull. World Health Organ 78(12): 246–50.

CSA (2007). Population and Housing Census of Ethiopia 2007: Statistical Report. Federal Democratic Republic of Ethiopia, Central Statistical Agency, Addis Ababa, Ethiopia. https://open.africa/dataset/d34d8ca9-97d3-4fb3-9249-47e2d3e3d322/resource/1af0860e-2fa6-4711-b9a7-53ddb04e721f/download/ethiopia-summary-and-statistical-report-of-the-2007-population-and-housing-census.pdf. Accessed 10 April 2021.

De Silva PM, Marshall JM (2012). Factors contributing to urban malaria transmission in sub-Saharan Africa: a systematic review. J. Trop. Med 18 ; 2012. https://doi.org/10.1155/2012/819563.

Dejenie T, Yohannes M, Assmelash T (2011). Characterization of mosquito breeding sites in and in the vicinity of Tigray microdams. Ethiop. J. Health Sci 21(1):57-66. https://doi.org/10.4314/ejhs.v21i1.69045.

Edillo, F. E., Tripét, F., Touré, Y. T., Lanzaro, G. C., Dolo, G., and Taylor, C. E (2006). Water quality and immatures of the M and S forms of Anopheles gambiae ss and An. arabiensis in a Malian village. Malar. J 5(1); 1-10. https://doi.org/10.1186/1475-2875-5-35.

Emidi B, Kisinza WN, Mmbando BP, Malima R, Mosha FW (2017). Effect of physicochemical parameters on Anopheles and Culex mosquito larvae abundance in different breeding sites in a rural setting of Muheza, Tanzania. Parasit Vectors 10(1):1-2. https://doi.org/10.1186/s13071-017-2238-x.

FMoH, 2014 . An epidemiological profile of malaria in Ethiopia. Ethiopia: Addis Ababa. https://webarchive.lshtm.ac.uk/www.linkmalaria.org/sites/link/files/content/country/profiles/Ethiopia%20Epi%20Report%20(240314).pdf. Accessed 10 April 2020.

Gari T, Lindtjørn B (2018). Reshaping the vector control strategy for malaria elimination in Ethiopia in the context of current evidence and new tools: opportunities and challenges. Malar. J 17(1):1-8. https://doi.org/10.1186/s12936-018-2607-8.

Getachew D, Balkew M, Tekie H (2020). Anopheles larval species composition and characterization of breeding habitats in two localities in the Ghibe River Basin, southwestern Ethiopia. Malar. J 19(1):1-3. https://doi.org/10.1186/s12936-020-3145-8.

Gillies, M.T. and Coetzee, M (1987). A supplement to the Anophelinae of Africa South of the Sahara. Publ S Afr Inst Med Res. 55: 1-143.

Gimnig, J. E., Ombok, M., Kamau, L., and Hawley, W. A. (2001). Characteristics of larval anopheline (Diptera: Culicidae) habitats in Western Kenya. J. Med. Entomol 38(2), 282-288. https://doi.org/10.1603/0022-2585-38.2.282.

Gimnig, J.E.,Ombok,M.,Otieno, S.,Kaufman,M.G.,Vulule, J.M. and Walker, E.D (2002). Density-dependent development of Anopheles gambiae (Diptera: Culicidae) larvae in artificial habitats. J. Med. Entomol. 39, 162–172.

Himeidan, Y. E., and Rayah, E. E. A. (2008). Role of some environmental factors on the breeding activity of Anopheles arabiensis in New Halfa town, eastern Sudan. EMHJ 14 (2), 252-259. https://apps.who.int/iris/handle/10665/117432.

Imbahale, SS, Paaijmans KP, Mukabana WR, Van Lammeren R, Githeko AK, Takken W. A. (2011). Longitudinal study on Anopheles mosquito larval abundance in distinct geographical and environmental settings in western Kenya. Malar. J 10(1):81. https://doi.org/10.1186/1475-2875-10-81.

Kenea, O, Balkew M, Gebre-Michael T (2011). 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 1; 48(2):85.

Killeen GF, Fillinger U, Knols BG (2002). Advantages of larval control for African malaria vectors: low mobility and behavioural responsiveness of immature mosquito stages allow high effective coverage. Malar. J 1(1):1-7. https://doi.org/10.1186/1475-2875-1-8.

Kindu M, Aklilu E, Balkew M, Gebre-Michael T (2018). Study on the species composition and ecology of Anophelines in Addis Zemen, South Gondar, Ethiopia. Parasit Vectors 1;11(1):215. https://doi.org/10.1186/s13071-018-2701-3.

Kipyab PC, Khaemba BM, Mwangangi JM, Mbogo CM (2015). The physicochemical and environmental factors affecting the distribution of Anopheles merus along the Kenyan coast. Parasit Vectors 1;8(1):221. https://doi.org/10.1186/s13071-015-0819-0.

Kweka EJ, Zhou G, Munga S, Lee MC, Atieli HE, Nyindo M, Githeko AK, Yan G (2012). Anopheline larval habitats seasonality and species distribution: a prerequisite for effective targeted larval habitats control programmes. PloS One 7(12). https://doi.org/10.1371/journal.pone.0052084.

Loha E, Deressa W, Gari T, Balkew M, Kenea O, Solomon T, Hailu A, Robberstad B, Assegid M, Overgaard HJ, Lindtjørn B (2019). Long-lasting insecticidal nets and indoor residual spraying may not be sufficient to eliminate malaria in a low malaria incidence area: results from a cluster randomized controlled trial in Ethiopia. Malar. J 1;18(1):141. https://doi.org/10.1186/s12936-019-2775-1.

Mala AO, Irungu LW (2011). Factors influencing differential larval habitat productivity of Anopheles gambiae complex mosquitoes in a western Kenyan village. J. Vector Borne Dis 1;48(1):52.

Mattah PA, Futagbi G, Amekudzi LK, Mattah MM, de Souza DK, Kartey-Attipoe WD, Bimi L, Wilson MD (2017). Diversity in breeding sites and distribution of Anopheles mosquitoes in selected urban areas of southern Ghana. Parasit Vectors 10(1):1-5. https://doi.org/10.1186/s13071-016-1941-3.

Mereta ST, Yewhalaw D, Boets P, Ahmed A, Duchateau L, Speybroeck N, Vanwambeke SO, Legesse W, De Meester L, Goethals PL (2013). Physico-chemical and biological characterization of anopheline mosquito larval habitats (Diptera: Culicidae): implications for malaria control. Parasit Vectors 1;6(1):320. https://doi.org/10.1186/1756-3305-6-320.

Messenger, L.A., Shililu, J., Irish, S.R., Anshebo, G.Y., Tesfaye, A.G., Ye-Ebiyo, Y., Chibsa, S., Dengela, D., Dissanayake, G., Kebede, E. and Zemene, E (2017). Insecticide resistance in Anopheles arabiensis from Ethiopia (2012–2016): a nationwide study for insecticide resistance monitoring. Malar. J 16(1): 1-14. https://doi.org/10.1186/s12936-017-2115-2.

Minakawa N, Mutero CM, Githure JI, Beier JC, Yan G (1999). Spatial distribution and habitat characterization of anopheline mosquito larvae in western Kenya. Am. J. Trop. Med. Hyg 1;61(6):1010-6. https://doi.org/10.4269/ajtmh.1999.61.1010.

Omlin FX, Carlson JC, Ogbunugafor CB, Hassanali A. (2007). Anopheles gambiae exploits the tree hole ecosystem in western Kenya: a new urban malaria risk?. Am J Trop Med Hyg 77:264–269. https://doi.org/10.4269/ajtmh.2007.77.264.

Oyewole IO, Momoh OO, Anyasor GN, Ogunnowo AA, Ibidapo CA, Oduola OA, Obansa JB, Awolola TS (2009). Physico-chemical characteristics of Anopheles breeding sites: Impact on fecundity and progeny development. AJEST 3(12).

Paaijmans KP, Imbahale SS, Thomas MB, Takken W (2010). Relevant microclimate for determining the development rate of malaria mosquitoes and possible implications of climate change. Malar. J 9(1):1-8. https://doi.org/10.1186/1475-2875-9-196.

Rejmánková E, Grieco J, Achee N, Roberts DR (2013). Ecology of larval habitats. Anopheles mosquitoes new insights into Malaria vectors. Intech 24:397-446. http://dx.doi.org/10.5772/5529.

Roux O, Robert V (2019). Larval predation in malaria vectors and its potential implication in malaria transmission: an overlooked ecosystem service? Parasit Vectors 1;12(1):217. https://doi.org/10.1186/s13071-019-3479-7.

Sattler, M. A., Mtasiwa, D., Kiama, M., Premji, Z., Tanner, M., Killeen, G. F., and Lengeler, C (2005). Habitat characterization and spatial distribution of Anopheles sp. mosquito larvae in Dar es Salaam (Tanzania) during an extended dry period. Malar. J 4(1): 1-15.

Scott JA, Brogdon WG, Collins FH (1993). Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am. J. Trop. Med. Hyg 1;49(4):520-9. https://doi.org/10.4269/ajtmh.1993.49.520.

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 (2003). Larval habitat diversity and ecology of anopheline larvae in Eritrea. J Med Entomol 40:921–929. https://doi.org/10.1603/0022-2585-40.6.921.

Soleimani-Ahmadi, M., Vatandoost, H. and Zare, M (2014). Characterization of larval habitats for Anopheline mosquitoes in a malarious area under elimination program in the southeast of Iran. Asian Pac. J. Trop. Biomed 4:S73-S80. https://doi.org/10.12980/APJTB.4.2014C899.

Taffese, H.S., Hemming-Schroeder, E., Koepfli, C., Tesfaye, G., Lee, M.C., Kazura, J., Yan, G.Y. and Zhou, G.F (2018). Malaria epidemiology and interventions in Ethiopia from 2001 to 2016. Infect. Dis. Poverty 7(1): 1-9. https://doi.org/10.1186/s40249-018-0487-3.

Tarekegn, M., Tekie, H., Dugassa, S. and Wolde-Hawariat, Y (2021). Malaria prevalence and associated risk factors in Dembiya district, North-western Ethiopia. Malar. J 20(1): 1-11. https://doi.org/10.1186/s12936-021-03906-9.

Teklu BM, Tekie H, McCartney M, Kibret S (2010). The effect of physical water quality and water level changes on the occurrence and density of Anopheles mosquito larvae around the shoreline of the Koka reservoir, central Ethiopia. Hydrol Earth Syst Sci 1;14(12):2595. 10.5194/hess-14-2595-2010.

Verrone GA (1962). Outline for the Determination of Malarial Mosquitoes in Ethiopia. Part II. Anopheline Larvae. Mosquito News. 22(4):394-401.

WHO (1975). Manual on practical entomology in Malaria. Part II. Methods and techniques. WHO Division of Malaria and other Parasitic Diseases. https://apps.who.int/iris/handle/10665/ 42481. Accessed 6 Feb2020

WHO (2013). Larval Source Management: A supplementary measure for malaria vector control: an operational manual. Geneva: World Health Organization. https://www.who.int/malaria/publications/atoz/larval_source_management_2pager_eng.pdf. Accessed 6 Feb 2020.

WHO (2018). World malaria report. Geneva; 2018 https://www.mmv.org/sites/default/files/uploads/docs/publications/WorldMalariaReport2018.pdf. Accessed 10 April 2020.

WHO (2021). World malaria report. Geneva; 2021. https://www.mmv.org/sites/default/files/uploads/docs/publications/World_Malaria_Report_2021.pdf. Accessed 4 May, 2022.

Yohannes M, Haile M, Ghebreyesus TA, Witten KH, Getachew A, Byass P, Lindsay SW (2005). Can source reduction of mosquito larval habitat reduce malaria transmission in Tigray, Ethiopia?. TM & IH 10(12):1274-85. https://doi.org/10.1111/j.1365-3156.2005.01512.x.

Download PDF

Editor assigned by journal
06 May, 2022
First submitted to journal
05 May, 2022

You are reading this latest preprint version

Habitat characteristics and spatial distribution of Anopheles mosquito larvae in malaria elimination settings in Dembiya District, Northwestern Ethiopia

Status:

Version 1

Abstract

Figures

Introduction

Methods

Description of the study area

Physico-chemical characteristics of larval habitats

Data analysis

Results

Discussion

Conclusions

Abbreviations

Declarations

References

Status:

Version 1