article: Characteristics of mosquito larval habitats are important in determining whether they can survive and successfully complete their developmental stages. Therefore, data on the ecological factors affecting mosquito density and abundance especially the physicochemical properties of water of their breeding sites, can possibly be helpful in implementing larval management programs. Mosquito larvae were collected using a standard 350 ml dipper from fixed habitats including: artificial pool, river edge, creek and etc, in 30 villages of 16 counties from May-December 2014. Water samples were collected during larval collection and temperature (°C), acidity (pH), turbidity (NTU), electrical conductivity (μS/cm), alkalinity (mg/l CaCO3), total hardness (mg/l), nitrate (mg/l), chloride (mg/l), phosphate (mg/l) and sulphate (mg/l) were measured using standard methods. Spearman correlation coefficient, Kruskal-Wallis test of nonparametric analysis, Chi-square (χ2) analysis, regression analysis and C8 interspecific correlation coefficient were used for data analysis. A total of 7,566 mosquito larvae belonging to 15 species representing three genera were collected from fixed larval breeding places. Culex pipiens was the dominant species except in four villages where An. maculipennis s. l. and Cx. torrentium were predominant. There was a significant positive correlation between the density of Cx. pipiens and electrical conductivity, alkalinity, total hardness and chloride, whereas no significant negative correlation was observed between physicochemical factors and larval density. The highest interspecific association of up to 0. 596 was observed between An. maculipennis s. l/An. pseudopictus followed by up to 0. 435 between An. maculipennis s. l/An. hyrcanus and An. hyrcanus/An. pseudopictus. The correlations observed between physicochemical factors and larval density, can possibly confirm the effect of these parameters on the breeding activities of mosquitoes, and may be indicative of the presence of certain mosquito fauna in a given region.
  Water characteristics of breeding places are important for oviposition and development of mosquitoes [1]. Different characteristics of the oviposition sites such as: vegetation, temperature, turbidity, pH, concentration of ammonia, nitrite and nitrate, sulphate, phosphate, chloride, calcium, and water hardness affect mosquito larval density [2–4]. Changing these factors in larval habitats may create conditions favorable or unfavorable for mosquito biology [5]. Temperature lower than 14–16°C and higher than 30°C reduces the rate of larval development of many species [2,6]. Larvae of most mosquito species are found in nature in pH 3. 3–10. 5 [2]. Studies in several micro and macrohabitats have revealed that distribution of Cx. pipiens, Cs. longiareolata, Cx. antennatus, Oc. caspius, Cx. vagans, Cx. decens, Cx. perexiguus, Cx. univittatus, An. multicolor are directly correlated with temperature, ammonia, nitrate, pH, dissolved oxygen and salinity [7,8]; An. culicifacies with temperature and dissolved oxygen [1,9]; Ae. albopictus with conductivity, total dissolved solids, nitrate, phosphate, sulphate, turbidity and salinity [10,11]; and An. varuna with calcium [1]. While larval densities of Cx. pipiens, Cx. perexiguus and Ae. albopictus are not directly affected by total nitrogen content, salinity and turbidity [8,10,11]. Dissolved nitrogen content can be a limiting factor on larval growth of the genus Aedes (formerly Ochlerotatus) by indirect effects on the trophic structure of tree-hole environments [12]. Although abundance of larvae can be associated with the soluble nitrate and phosphate levels in an area [13], it can equally be independent of those in another area because the nitrogen or phosphorus in water bodies can result in eutrophication and oxygen depletion, harmful algal blooms, toxic effects on fish and some aquatic organisms and overall reductions in aquatic biodiversity [4,14]. Therefore other factors can be effective in predicting larval density in different areas [15]. For example, the vegetation of mosquito larval habitats is considered as an important factor in the process of egg-laying and density of mosquito larvae [16]. Tall emergent aquatic plants can cover the surface and decrease mosquito larval density by acting as a barrier for egg-laying female. They may cause microbial growth and produce a high variety of predators [17]. Not much data currently exists regarding the physiochemical characteristics of mosquito larval habitats in Iran. Most of the available data is on anopheline mosquito larval habitats. The distribution of An. sacharovi is significantly associated with calcium bicarbonate, sodium sulphate and salinity in larval habitats in Ardebil Province, northwestern Iran[18]. Larval density of An. culicifacies s. l. , An. dthali, An. stephensi, An. superpictus, An. fluviatilis s. l. , An. turkhudi, An. moghulensis, and An. apoci is associated with temperature, EC, alkalinity, chloride and sulphate, total hardness, and dry residues in Bashagard and Rudan district, southern Iran [15,19]. A significant difference was observed between the density of An. culicifacies and phosphate, calcium and EC; An. turkhudi and pH, total hardness and nitrate; An. superpictus and total hardness and nitrate; An. stephensi and nitrate, An. multicolor and pH and sulphate in Iranshahr, southwestern Iran [20]. Temperature, pH, turbidity, EC, TDS, alkalinity, total hardness, calcium, chloride, fluoride, nitrite, nitrate, phosphate, sulphate were not significantly correlated with An. claviger, An. marteri, An. superpictus, An. turkhudi, Cx. arbieeni, Cx. hortensis, Cx. mimeticus, Cx. modestus, Cx. pipiens, Cx. territans, Cx. theileri, Cs. longiareolata, Cs. subochrea, Oc. caspius s. l. in Qom Province, central Iran [21] To date, there is no information on physicochemical characteristics of mosquito larval habitats and coefficient of interspecific association (C8) in Mazandaran Province and this is the first such study in the province. Differences in environmental and geographical characteristics of the North of Iran in comparison with other provinces of the country have caused a fertile environment for the development of mosquitoes. The area had a history of diseases transmitted by vectors, the most important of which has been Malaria. Mazandaran has a unique environment fit for migratory birds and thousands of different wild bird species that spend winter in numerous fresh water lakes and wetlands across the province. West Nile virus is circulating in the province between the migratory birds and humans mostly by Culex mosquitoes [22,23]. Also, the environment may allow the invasion and establishment of Aedes vectors of dengue and Zika based on the risk assessment of the latter, a huge national research activity is underway to verify the risk. Information on the ecological factors affecting mosquito larval biology such as the physicochemical properties of the water of the breeding places and interspecific associations are important in survival, spatio-temporal distribution [24], biodiversity, affinity and association indices of disease vectors [25–27]. The information may serve as the basis for designing and implementation of adequate vector control programs [28]. Despite the voluminous literature on the distribution of mosquito larvae and physicochemical factors, the data set seems to be inconclusive in leading to a prediction of the presence of larvae in different habitats. Therefore more studies and systematic reviews with proper generalization are required. The present study was conducted to document the relationship between physicochemical characteristics of mosquito larval habitats and the presence or absence of a given species, density or diversity and interspecific associations of larvae in different habitats in Mazandaran Province.
  This research has been approved by the Ethic committee of Mazandaran University of Medical Sciences under the code 1017. The study was carried out from May to December 2014 in Mazandaran Province located within the latitude of 35° 47′- 36° 35′ N and longitude of 50° 34′–54° 10′ E (Fig 1). The province has an area of 23756. 4 km2 and a population of 3,073,943 (2011 census) [29]. It is located between Golestan Province in the East, Gilan Province in the West, Caspian Sea in the North and Tehran and Semnan Provinces in the South. The inhabitants are mainly involved in agriculture, animal husbandry, production of farmed fish and tourism industry. The climate varies from mild and humid in the Caspian Sea shore to moderate and cold in mountainous regions. The minimum and maximum mean annual temperatures and rainfall ranges between 1. 2–29. 2°C and 0–755. 6 mm, respectively. Mosquito larvae were collected by standard 350 ml dipper from fixed oviposition sites (breeding places with permanent water during the sampling period) such as artificial pools, river edge, creek, marsh, large metal bucket, abandoned wells, water canals, pit with plastic floor (Fig 2), in 30 villages of 16 counties across the province (Fig 1). In each village, one fixed station was selected and visited for larval collection once a month. It is important to mention that one hundred staff members of the Mazandaran health centers were recruited and trained during two theoretical and practical workshops before undertaking sample collection. We met with county and village councils to seek cooperation of the villagers for assistance with sampling team. Collected larvae were conserved in lactophenol and were transferred to the Medical Entomology Laboratory at Faculty of Health, Mazandaran University of Medical Sciences. In the laboratory, microscope slides were prepared from each of the specimens using Berlese’s medium and larvae were identified by morphological characteristics according to appropriate keys [30,31]. Water samples were collected from the same 30 fixed larval habitats in which dip samples had been collected in designated villages and were kept in one liter appropriately labeled polyethylene bottles. The bottles were placed inside an ice box and were sent for analysis to the Faculty of Health. Temperature (°C), acidity (pH), turbidity (NTU), electrical conductivity (EC) (Micro-Simens/cm) (μS/cm), alkalinity (mg/l CaCO3), total hardness (mg/l), nitrate (mg/l), chloride (mg/l), phosphate (mg/l) and sulphate (mg/l) of the water were measured. Water temperature and pH were measured on site using a thermometer and pH probe (Eutech-Cyberscan PH5500), turbidity using a turbidimeter device (2100P Portable Turbidimeter at Hach), EC using a conductometer device (EC LYTIC-AQUA (CON200) before dipping. Alkalinity, total hardness and chloride were determined using direct titration techniques and nitrate, sulphate and phosphate were measured using a spectrophotometer device (Perkinelmer UV/Visible Lambda EZ201 and HACH). All analyses were conducted according to the standard methods used in Rice et al. [32]. The species data from each site was summed up and collectively reported in a graph. The means and standard deviations of physicochemical parameters of each breeding site was calculated using SPSS software version 19. The presumption for normality was tested using the Shapiro-Wilk test. The physicochemical parameters among sites were compared by Kruskal-Wallis test of nonparametric analysis. Chi-square (χ2) analysis was performed to determine whether there was any significant difference in distribution of the species in different counties. The Spearman correlation coefficient was used to examine the relation of the mosquito larval densities to the physicochemical factors adjusted by types of habitats. Regression analysis between larval densities and physicochemical factor were also performed to clarify the relationship and r-squared values were calculated. Hurlbert’s coefficient of interspecific association (C8) was used to measure the associations between co-occurring species using presence-absence data. Values of C8 range from -1 to +1 for negative and positive associations, respectively. Positive associations between species can probably show a common habitat preference or interspecific attraction, whereas negative associations may reveal different habitat preferences or interspecific repulsion. The formula used for calculating C8 is as follows: C8=ad−bc|ad−bc||Obsx2−Minx2Maxx2−Minx2| where a, b, c, and d are the values in four cells of a 2 × 2 contingency table; Obs χ2 is referred to the value of χ2 associated with the observed values of a, b, c and d; Max χ2 is referred to the value of χ2 when a is as large (if ad ≥bc) or as small (if ad < bc) as the marginal totals of the 2 × 2 table permit; Min χ2 is referred to the value of χ2 when the observed a differs from its expected value (â) by less than 1. 00 [33].
  A total of 7,566 mosquito larvae from three genera and 15 species were collected from eleven different types of fixed oviposition sites. Using random effect model, it is observed that habitats type affect the abundance of species (P < 0. 05). Culex pipiens (56. 22%), Cx. tritaeniorhynchus (10. 05%) and An. maculipennis s. l. (10. 50%), were the most abundant species observed in artificial pool. whereas, An. marteri (0. 06%), Cx. hortensis (0. 15%) and Cs. morsitans (0. 08%) were relatively uncommon in habitats such as concrete canal, artificial pool and tire; respectively (Fig 3). Using the coefficient of interspecific association (C8), significant positive associations was determined between An. maculipennis s. l/An. hyrcanus and An. hyrcanus/An. pseudopictus of up to 0. 435. This coefficient was up to 0. 596 between An. maculipennis s. l/An. pseudopictus (Table 1). There are significant differences in the distribution of the species in spatial scale (P < 0. 05). Culex pipiens had the largest distribution and was the dominant species in all villages except Tileno, Zavat, Asad Abad and Shah Mansur Mahale. The highest and the lowest number and percentage of Cx. pipiens were recorded in the villages of Firozabad (558,86. 11%) and Tileno (3,3. 26%), respectively. Culex tritaeniorhynchus was dispersed in all villages except Abbasa and Soleymanabad and showed the highest abundance in the village of Al Tappeh (95,24. 05%). Anopheles maculipennis s. l. were the dominant anopheline in the villages of Tileno, Zavat and Shah Mansur Mahale and showed the highest dispersion rate compared to the rest of the Anopheles. Culiseta morsitans and An. marteri showed the lowest distribution in the villages of Bishekola (0. 76%) and Chalmardi (0. 40%). Data regarding distribution of other species in spatial scale are shown in Fig 4. Means and standard deviations of the physicochemical factors from larval habitats in different villages were calculated for each species. Culiseta morsitans was found in breeding sites with higher pH up to 7. 44. It prefers lower temperature, alkalinity, nitrate, chloride and phosphate. Anopheles marteri prefers higher turbidity, alkalinity, total hardness, phosphate and sulphate, while An. claviger, Cx. mimeticus and Cx. perexiguus breed in habitats with higher temperature, EC and chloride, respectively. Culex mimeticus, Cx. hortensis, Cx. torrentium and An. claviger were collected in habitats with the lowest ranges of pH (6. 88±0. 40), turbidity (12. 50±0. 70), EC (69. 42±1051. 15), total hardness (192) and sulphate (14. 18), respectively (Table 2). Ranges of the physicochemical parameters of larval habitats of different species are presented in Table 2. Positive correlation exists between the larval abundance of Cx. pipiens and the physicochemical characteristics including EC, alkalinity, total hardness and chloride with Spearman rank correlations of 0. 575 (P<0. 001), 0. 617 (P<0. 001), 0. 495 (P<0. 005) and 0. 539 (P<0. 002), respectively. However, there was no significant negative correlation between physicochemical characteristics and larval abundance (Table 3). The regression relationship (R2) between Cx. pipiens and EC, alkalinity, total hardness and chloride were 0. 11,0. 20,0. 01 and 0. 06, respectively which is shown by Scatter plot (Fig 5).
  In conclusion, based on the findings of the present study, physicochemical factors of breeding sites including EC, alkalinity, total hardness and chloride may determine the distribution and abundance of Cx. pipiens in the area. Although seems important and expected to be observed, a significant positive correlation has not been detected for the rest of the species present in the study area. High interspecific association between pair of species An. maculipennis s. l/An. pseudopictus show that these species have common needs and adaptability for sympatry. These findings could be useful in comprehending the ecology of mosquito larvae that may be beneficial in designing and implementing larval control programs. This is the first study looking at the physicochemical parameters associated with larval abundance and diversity in Mazandaran Province. Further studies are needed to enable the use of this information in vector control programs with confidence. It should be noted that there is a limitation to this study and that is we measured the physicochemical parameters only once in the whole of the study areas and sampling locations whereas the density of mosquito larvae measured every month. The possibility of relation of larval habitats physicochemical parameters and larval density fluctuation might have well been established should the trend of changes of the physicochemical parameters be determined, this has further research implications.
summary: Determination of association between mosquito larval abundance and physicochemical factors in aquatic habitats in Mazandaran Province was the purpose of this study. In 30 villages of 16 counties in the province, Culex pipiens, Cx. tritaeniorhynchus and An. maculipennis s. l. were the most abundant species discovered in larval habitats. The pair of species An. maculipennis s. l/An. pseudopictus has the highest interspecific association in the area. Culex pipiens showed the greatest compatibility compared with other species with a range of densities from the lowest to the highest in different habitats containing varying levels of physicochemical factors. Measuring physicochemical factors of water in larval habitats indicated that values of chloride, electrical conductivity, alkalinity and total hardness were significantly associated with the larval abundance in different larval habitats. The electrical conductivity, alkalinity, total hardness and chloride had significant positive correlation with the density of Cx. pipiens, whereas no negative correlation was obtained between larval density and physicochemical factors. These differences in preferences and association with varying levels of physicochemical factors can have vector-borne diseases control implications, because habitat changes through the manipulation of these factors to the detriment of vector species can ultimately lead to reduced abundance of vector species. The data may also be useful to predict changes in the densities of mosquito larvae or their transition should natural or artificial environmental modifications lead to alteration in the physicochemical factors.
keywords: invertebrates
  medicine and health sciences
  ecology and environmental sciences
  chemical compounds
  phosphates
  salts
  animals
  nitrates
  developmental biology
  materials science
  habitats
  materials physics
  insect vectors
  infectious diseases
  turbidity
  sulfates
  life cycles
  chemistry
  chlorides
  disease vectors
  insects
  arthropoda
  physics
  mosquitoes
  biology and life sciences
  species interactions
  physical sciences
  larvae
  organisms
title: Correlation between mosquito larval density and their habitat physicochemical characteristics in Mazandaran Province, northern Iran