Mosquitoes cause significant human health issues. However, very few studies have attempted to examine the question of how abiotic and biotic factors affect the abundance of Culicidae in the larval habitat. The objective of this study was to analyze the influence of biotic and abiotic factors on the increase of the most common mosquito species in the Fez region (Central Morocco). Larvae mosquitoes were sampled by standard dipping technique in four different types of macrohabitats, between November 2015 and November 2016. Each mosquito specimen was morphologically identified by the Moroccan Culicidae key and the Brunhes key. The analysis was done using R analysis software. We collected a total of 772 mosquito larvae belonging to nine different species, five of which are considered of medical interest. Culex pipiens (Linnaeus, 1758), known as the major vector in the transmission of West Nile virus fever (WNV), was the most common species of all mosquito larvae collected. The results of Poisson regression analysis showed that factors such as the presence of green filamentous algae, vegetation cover, and debris were found to be positively significant in the distribution of the genus Culex. However, there was insufficient evidence to determine the parameters that are capable of estimating the abundance of Anopheles. The findings have also estimated that biotic and abiotic factors can lead to significant variation in the abundance of Culex perexiguus (Theobald, 1903), Culex theileri (Theobald, 1903), and Culex pipiens (Linnaeus, 1758). Identifying the priority parameters governing the proliferation of mosquitoes in the region of Fez can be one of the key elements for better vector control.

1. Introduction

Mosquitoes infect humans and other vertebrates. Their distribution is to a large extent influenced by climatic conditions and species habits across the globe [1]. Interactions between biotic and abiotic factors can conduct to significant variation in the reproductive landscape of insects [2]. They differ in their preference for the type, size, turbidity, algal cover and stability of the habitat [3]; these factors can lead to determining the density, size and vector competence [3]. Surveying biotic and abiotic factors for various mosquito fauna make it easier to monitor the potential modifications of larval habitats affected by rains, global climate change, and man-made activities [46].

In Morocco, many species of Culicidae family have been reported as important vectors of infectious diseases, including the presence of the major vectors of malaria (Anopheles labranchiae and Anopheles sergentii (Theobald, 1907)) [7, 8], the Aedes albopictus, responsible for the transmission of dengue, chikungunya, and Zika viruses [9], and Culex pipiens, the potential vector of both West Nile virus [10] and Rift Valley fever virus [11].

Antecedent studies on larval breeding sites have been a valuable source of information regarding the larval biology and ecology of different mosquito species. The search has provided that the factors salinity, pH, and water body permanence determine the occurrence and distribution of larval [12], and water depth and temperature influence larval propagation [13]. Moreover, the human factors like land-use and land-cover changes represent the principal key factor of the mosquito larvae occurrence [14].

Few studies have analyzed the influence of biotic and abiotic factors on the abundance of mosquito species in the Fez region (north central region of Morocco).

The acquisition of knowledge about the larval habitats and estimates of the biotic and abiotic factors associated with the distribution of mosquitos can contribute to an efficient control method.

2. Materials and Methods

2.1. Study Area

The study was carried out in the Fez region located in the northern central part of Morocco, characterized by mild and sunny Mediterranean climate and a large human population (1150131 in 2014) (Figure 1).

The study area was limited to four habitat types (dam, swamp, source, and El Oued), categorized into artificial and natural active larval habitats. The selection of habitats was based on their permanent nature and the presence of larvae identified in a preliminary survey.

2.1.1. Lgaâda Dam

The dam has a height of 30 m, with 2.90 Mm3 of retention capacity and 10205.21 Mm3 of total cumulated capacity. This is an artificial reservoir fed by runoff water which contributes to the supply of drinking water and allows for additional irrigation for crops during dry periods. It is located at 400 m above mean sea level, 34°01.155′N of latitude, and 004°57.213′W of longitude.

2.1.2. Awinat Elhajaj

This natural site is a shallow swamp tributary to channels of Lgaâda dam water. It is located at an altitude of 382 m, 34°01.226′N of latitude, and 004°57.105′W of longitude.

2.1.3. Douwar Lhandiya

It is an artificial source located at 410 m of altitude, 34°02.049′N of latitude, and 004°58.512′W of longitude. The station choice was made as it is a periurban area, close to the university campus Dhar El Mahraz of Fez, where foreign students with large numbers can be carriers of infectious diseases, typically imported from areas of endemicity.

2.1.4. Jnan EL Alami Oued

It is located downstream of the main water body crossing the city of Fez (Sebou Oued) at an altitude of 242 m, 34°04.421 N of latitude, and 004°57.710′W of longitude. It is one of the most polluted rivers in the region and contains pollutants from various sources such as sewage, fertilizers from agricultural fields, and industrial waste.

2.2. Larval Habitat Characterization

The weather and the characteristics of the macrohabitats were recorded, including water use (livestock drink, industry, washing, and irrigation), average depth (measured in three various points of each pool), water flow, slope (zero slope, median slope, and strong slope), color and origin of water (source, river, dam, dam pipeline, and wastewater discharges), proximity to mosquito breeding habitat, surface debris, presence of vegetation, and filamentous algae. In addition, pH, temperature, dissolved oxygen, and total dissolved solids were taken using Consort Multiparameter Analyser C561. The salinity was taken at the same time as electrical conductivity using water quality instrument YSI scientific, and turbidity was recorded by Lovibond ® turbidimeter Turbidirect.

2.3. Mosquito Sampling

Larval sampling was carried out every two weeks using a dipping method for a period of one year from November 2015 to November 2016. The sampling was done in about ten ladle dives in several levels. The specimens were preserved in ethanol (70%).

The typical aquatic habitats sampled are shown in Figure 2. Mosquito larvae collections of each site were labeled with a sampling date, stored in glass flacons, and transferred to the Laboratory of Biotechnology and Preservation of Natural Resources. All mosquito samples were morphologically identified by microscopy at the species level, using the Moroccan key of identification of Culicidae [15] and Mediterranean Africa mosquito identification software [16].

2.4. Statistical Analysis

The statistical analysis has been used to assess the relationships between biotic and abiotic factors and abundance of the most common species especially the species of medical interest. R Software (version 3.4.3 (2017)) was used for data processing. Data exploration and modeling were used to detect the interactions between variables explained (number of species) and different explanatory variables (biotic and abiotic factors). The “number of each species” is the dependent variable, whereas “presence of vegetation,” “water used for washing and water for livestock,” “water used for industry,” “water used for irrigation,” “debris,” “color,” “weather,” “slope,” “proximity of industrial establishments,” “proximity of construction,” and “presence of green filamentous algae” are all nominal independent variables. The quantitative explanatory variables were taken into account according to nine physicochemical parameters such as water temperature (°C), pH, salinity (‰), dissolved oxygen (PPM), electrical conductivity (µs/cm), flow rate m³/s, total dissolved solids (g/l), turbidity (NTU), and water depth (cm).

The following models were used:(1)the logistic regression for low number species like Anopheles maculipennis s.l., Anopheles sergentii (Theobald, 1907), Culiseta longiareolata (Macquart, 1838) and Uranotaenia unguiculata (Edwards, 1913),(2)the Poisson regression for more abundant mosquito species, especially Culex theileri (Theobald, 1903), Culex perexiguus (Theobald, 1903), and Culex pipiens (Linnaeus, 1758).

The “vif” function in the car package was used to select independent explanatory variables and testing collinearity, while “glm” and “step” functions are used to select the best models.

3. Results

A total of 772 specimens belonging to nine species were collected. Among these species three species of anophelines were identified, including An. maculipennis s.l., An. sergentii (Theobald, 1907) and An. Cinereus. The Culex was more prevalent, wherein Culex pipiens (Linnaeus, 1758) was more common than Culex perexiguus (Theobald, 1903) and Culex theileri (Theobald, 1903), and Culex hortensis was represented by single species. The genera of Culiseta and Uranotaenia were represented by one species (Culiseta longiareolata (Macquart, 1838) and Uranotaenia unguiculata (Edwards, 1913)). The highest species richness (seven species) was observed in Douwar Lhandiya source and Awinat Elhajaj. The most abundant populations were collected, respectively, in Douwar Lhandiya source (448 specimens) and Awinat Elhajaj (217 specimens). The low abundant populations were collected in Lgaâda dam with 448 specimens and Jnan El Alami Oued with 28 specimens. The species common at all four sites was Cx. perexiguus (Theobald, 1903) (Figure 3).

The results of the logistic regression showed that the species Anopheles maculipennis s.l., Anopheles sergentii (Theobald, 1907), Culiseta longiareolata (Macquart, 1838), and Uranotaenia unguiculata (Edwards, 1913), which are less abundant in the present study, had no statistically significant relationship between their abundance and the biotic and abiotic factors.

The results of Poisson regression analysis (Table 1) showed that the species Cx. perexiguus (Theobald, 1903) has been positively correlated with site proximity to industrial establishments and the presence of debris and green filamentous algae and significantly increased with increasing depth and conductivity but negatively correlated with water flow and dissolved oxygen. Cx. theileri (Theobald, 1903) larval density was observed to negatively correlate with the total dissolved solids and proximity of the sites to building construction and positively correlate with dissolved oxygen, pH, and salinity, The most important observation was that, in relation to the presence of green filamentous algae and vegetation cover, they were statistically positively significant. The result shows negative, highly significant associations between the distribution of Cx. pipiens (Linnaeus, 1758) and the water used for washing and watering livestock, conductivity, and dissolved oxygen and slightly significant associations with pH, total dissolved solids, and water flow. However, the presence of emergent plants, the green and incolor water, and the water used for irrigation were slightly positively associated. In addition, a highly positively and significant relationship was found in terms of zero slope, salinity, and water depth.

The results of the rare observations of An. maculipennis s.l., An. sergentii (Theobald, 1907), Cs. longiareolata (Macquart, 1838), Ur. unguiculata (Edwards, 1913), An. Cinereus (Theobald, 1901), and Cx. hortensis (Ficalbi, 1889) were not explained because we do not have an explicit form of the relationship between variables explained and a set of explanatory variables.

4. Discussion

As shown in the statistical analyses, both abiotic and biotic parameters such as component environments and physicochemical parameters strongly affected the abundance of Culicidae, which is already mentioned in literature [3].

Our study showed that the abundance of Cx. pipiens (Linnaeus, 1758) larvae within the breeding site was not correlated with water pH which is in agreement with Rydzanicz et al. (2016) [17]. These authors suggest a strong correlation between Cx. pipiens (Linnaeus, 1758) larval abundance and the physical and hydrological characteristics of the aquatic systems, like electrical conductivity. By contrast, this finding was not positively correlated in our data. Other recent studies by Amara Korba et al. (2016) showed that pH, temperature, and dissolved oxygen exhibited comparable values without any relation with larval densities of Cx. pipiens (Linnaeus, 1758) [18], and there were no marked differences in the presence of this species and water flow, distance to the nearest house, artificial or natural larval breeding sites, soluble solid, and vegetation [19]. These results were not in complete agreement with our result. However, other results from another country were in agreement with our study’s findings [20]. It was confirmed that members of Cx. pipiens (Linnaeus, 1758) were potentially affected by water flow, water depth, salinity, soluble solid, and vegetation in artificial or natural larval breeding sites [20].

Larvae of Cx. pipiens (Linnaeus, 1758) are able to breed in a great variety of habitats. However, a highly polluted breeding site can inhibit their development [18]. Overall, our findings confirm that polluted water characterized by a yellow color affected the abundance of this species in larval habitat.

Other important factors were strongly associated with the presence of Cx. perexiguus (Theobald, 1903) larvae and Cx. theileri (Theobald, 1903). Both of them were increasing in the presence of algae [3]. Our results confirmed that species of mosquitoes (Cx. perexiguus (Theobald, 1903), Cx. theileri (Theobald, 1903)) are positively affected by green filamentous algae [3]. In the literature, it is confirmed that mosquito species were generally more present in natural areas than in urban and rural landscapes [2123]. In our work, the data obtained demonstrate that these two species were less common in the sites near to residential buildings.

Distribution and abundance of Cx. perexiguus (Theobald, 1903) were recognized to show a close relationship with industrial activities and correlate with slow-running water environments, debris, electrical conductivity, and water depth. The last parameter was confirmed to be correlated with this species in another study [24].

Considering our results, water environments with vegetation, dissolved oxygen, pH, salinity, and algae could imply faster larvae development of Cx. theileri (Theobald, 1903). Similar findings have been reported in previous studies conducted in Iran, which reported that the larvae of Cx. theleri were found with vegetation in natural and permanent habitats [25].

Considering the low number and the rare observations of the larvae (Cs. longiareolata (Macquart, 1838), Ur. unguiculata (Edwards, 1913), Cx. hortensis (Ficalbi, 1889), and Anopheles), the relationship between these species and abiotic and biotic parameters remain unclear. However, the exclusive presence of Cs. longiareolata (Macquart, 1838) in Douwar Lhandiya’s source confirms that this species is found only in fresh water pools [26]. All of these species were present in low-water habitats, which confirm that all species are more likely to occur in shallower water [27].

5. Conclusion

Our data provide some important findings into the larval habitat. In conclusion, our findings suggest that species of the genus Culex differ in the factors potentially affecting their presence. The preferred habitats of Cx. perexiguus (Theobald, 1903) were generally characterized by proximity to industrial establishments and the presence of debris and green filamentous algae. The population of this species was increased with increasing depth and high conductivity.

The species of Cx. theileri (Theobald, 1903) are frequently correlated with dissolved oxygen, pH and salinity, and they were plentiful in the presence of green filamentous algae and vegetation cover.

The most important factors that greatly affect the abundance of Cx. pipiens (Linnaeus, 1758) in the region were zero slope, salinity, and water depth.

Overall, our findings provide some results about the influence of biotic and abiotic parameters on the distribution of mosquito species. The results of this study could be useful for the authorities in entomological surveillance to efficiently develop important vector controls.

Further work must be done in order to better characterize the dynamics of these mosquitoes in this region.

Data Availability

The data used in this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


The authors are grateful to the members of the Laboratory of Medical Entomology, National Hygiene Institute, Rabat.