J Arthropod-Borne Dis, December 2019, 13(4): 378–390 SH Nikookar et al.: First Report of … 378 http://jad.tums.ac.ir Published Online: December 31, 2019 Original Article First Report of Biochemical Mechanisms of Insecticide Resistance in the Field Population of Culex pipiens (Diptera: Culicidae) From Sari, Mazandaran, North of Iran Seyed Hassan Nikookar1; Mahmoud Fazeli-Dinan2; Seyyed Payman Ziapour3; Fatemeh Ghorbani4; Yaser Salim-Abadi5; Hassan Vatandoost4,6; Ahmad Ali Hanafi-Bojd4; *Ahmadali Enayati1 1Department of Medical Entomology and Vector Control, School of Public Health and Health Sciences Research Center, Addiction Institute, Mazandaran University of Medical Sciences, Sari, Iran 2Department of Medical Entomology and Vector Control, School of Public Health and Health Sciences Research Center, Mazandaran University of Medical Science, Sari, Iran 3Department of Parasitology, Zoonosis Research Center, Pasteur Institute of Iran, Amol, Iran 4Department of Medical Entomology and Vector Control, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran 5Department of Health Services and Health Promotion, School of Health, Rafsanjan University of Medical Sciences, Rafsanjan, Iran 6Department of Chemical Pollutants and Pesticides, Institute for Environmental Research, Tehran University of Medical Sciences, Tehran, Iran (Received 03 June 2019; accepted 29 Dec 2019) Abstract Background: Culex pipiens play an important role in transmission of infectious diseases. Vector control by chemical pesticides, leads inevitably to resistance development. Understanding the underlying resistance mechanisms can help improve the control programmes and insecticide resistance management. Methods: The total contents of cytochrome p450s and the activities of glutathione S-transferases, alpha- and beta- esterases and inhibition rates of acetylcholine esterase (by propoxur) were measured in the field population of Cx. pipiens collected from Sari County, North of Iran, in 2016 and the results were compared with those of the laboratory susceptible strain according to the biochemical assay methods of WHO for adult mosquitoes. Independent sample t-test was used to compare the mean values of enzyme activities/contents between filed and laboratory susceptible popula- tions. Results: The enzyme ratio of cytochrome p450s, alpha- and beta-esterases in the field population was 2.07, 3.72 and 1.36 respectively when compared with the results of the laboratory population. Although not statistically significant, the mean GSTs activities in the field population was marginally less than the laboratory population (ER=0.92). Ace- tylcholinesterase was insensitive to propoxur in 62.82% of the individuals of the tested field population. There was a significant difference (P< 0.05) between all values of the activities/contents of the enzyme in the field population except for GSTs compared with the laboratory susceptible strain. The highest enzyme activity was related to alpha esterase. Conclusion: The present study showed a range of metabolic mechanisms, comprising p450s and esterases combined with target site insensitivity of AChE, contributing to organophosphate, carbamate and pyrethroid resistance in the field population of Cx. pipiens. Keywords: Culex pipiens; Insecticide resistance; Enzyme; Iran Introduction Mosquitoes are the most important arthro- pods of medical importance transmitting vari- ous diseases to humans and animals. Among mosquitoes, the genus Culex, especially, Cx. pipiens complex are known vectors of several human pathogens including West Nile virus, *Corresponding author: Dr Ahmad Ali Enayati, E-mail: ahmadali_enayati@yahoo.com http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2019, 13(4): 378–390 SH Nikookar et al.: First Report of … 379 http://jad.tums.ac.ir Published Online: December 31, 2019 St. Louis encephalitis virus, Rift Valley fever virus, Japanese encephalitis virus, Wuchereria bancrofti (1-3) as well as zoonotic pathogens including dog heartworm (4) and avian malaria (5). Mosquito-borne diseases immensely affect public health and cause problems in terms of economy and development all over the world (6). Culex pipiens is widely distributed in al- most all continents throughout the world (7, 8). Lately, Cx. pipiens has been arrived on New- foundland’s Avalon Peninsula, Canada and southern Sweden (9, 10). It shows a vast geo- graphical distribution in most parts of Iran (11- 13). The distribution, dominance and high abun- dance of Cx. pipiens have been reported from northern Iran (14-18). The use of chemical in- secticide is one of the most widely practiced strategies for mosquito-borne diseases control (19). Given the important role of pesticides in mosquito control programms, the development of resistance to insecticides is of great concern for diseases control (20, 21), recently led to many outbreaks of mosquito-borne diseases (22). Resistance to insecticides in Cx. pipiens has been reported from many parts of the world (23, 24). In Iran, approximately 14,000 tons of active ingredient agricultural pesticides were imported or produced (25). Population and environmen- tal studies have alarmingly revealed high levels of pesticide residues in the environment (25), it could possibly be a cause of concern for devel- opment of resistance to pesticide in mosquito populations that have a water-related life cycle. Culex pipiens was resistant to DDT, propoxur, lambda-cyhalothrin, cyfluthrin and deltamethrin (26–29), tolerant to deltamethrin (27, 28) and susceptible to malathion (28). In a recent study, however, revealed relatively high resistance to propoxur, malathion, fenitro- thion and dieldrin in the field populations of Cx. pipiens from Sari County, North of Iran (30). The two basic mechanisms of insecticide resistance in mosquitoes are metabolic resistance and target site insensitivity. Metabolic resistance is mainly caused by elevated levels of enzyme activities that lead to detoxification or seques- tration of insecticides before they reach to their target site (31). Glutathione S-transferases (GSTs), cytochrome p450s and esterase’s (ESTs) are three main enzyme groups responsible for metabolic resistance to the main classes of in- secticide applied against insects of public health importance (32, 33). Elevated level of active- ties of these enzymes was shown to be respon- sible for resistance to pesticides in a variety of mosquitoes including Anopheles stephensi (34, 35), Aedes aegypti (36, 37) and An. culicifacies (38). Biochemical studies showed the involve- ment of cytochrome p450s in the metabolic resistance to four major insecticide classes in An. funestus, An. stephensi and An. darlingi (39–41). Increased levels of esterases have been associated with resistance to organophosphates, carbamates and pyrethroids in Cx. quinquefas- ciatus and Helicoverpa armigera (42). Involve- ment of glutathione S-transferases has been ob- served in resistance to organophosphates, or- ganochlorines and pyrethroid insecticides in An. subpictus, Cx. pipiens and Rhipicephalus bur- sa, Musca domestica, Drosophila melanogaster and Nilaparvata lugens (43–45). Glutathione S- trasferases are principally responsible for re- sistance to organochlorine insecticides and their involvement in resistance to other groups of in- secticides is secondary (44). There is also evi- dence for involvement of the insensitive ace- tylcholinesterase (AChE) in resistance to organ- ophosphates and carbamates in Cx. pipiens, Ae. aegypti and Ae. albopictus (46, 47). Esterases showed more important role in the detoxification of organophosphate insecticides used against Cx. pipiens in Egypt (48). DDT, organophosphate and pyrethroid resistance was characterized in Cx. pipiens from Turkey and cytochrome p450s and esterases were involved in resistance (21). Overproduced esterases and insensitive AChE were involved in enzymatic resistance of the field populations of Cx. pipiens from Grand Tunis area, Northeast Tunisia (49). Other resistance mechanisms, including target- site resistance mutations, can be involved in in- secticide resistance in Cx. pipiens populations http://jad.tums.ac.ir/ https://www.google.com/url?sa=t&rct=j&q=&esrc=s&source=web&cd=4&cad=rja&uact=8&ved=0ahUKEwjHuOSe3sbUAhWGcRQKHTznCd0QFgg2MAM&url=http%3A%2F%2Fwww.mosquito.org%2Fmosquito-borne-diseases&usg=AFQjCNHPKnitXjPoWEVYwXT4ivRjHS8JxA J Arthropod-Borne Dis, December 2019, 13(4): 378–390 SH Nikookar et al.: First Report of … 380 http://jad.tums.ac.ir Published Online: December 31, 2019 from the Mediterranean region (50, 51). The kdr mutation (L1014F) has been reported in association with resistance to lambda-cyhalo- thrin from Cx. pipiens in Morocco region (51). The L1014S mutation has been observed in Cx. pipiens pallens from Japan and China. The L1014C and V1016G mutations have only been reported in Cx. pipiens form molestus from Chi- na and Saudi Arabia, respectively (23). The kdr mutation along with elevated levels of p450s and GSTs, are associated with pyrethroid resistance in the populations of Cx. pipiens from China (52). Different amplified esterases (encoded by the Ester super locus) and two substitutions on the AChE1 (encoded by the ace- 1gene) including the G119S and the F290V substitutions have been reported to confer re- sistance to a great variety of organophosphate and carbamate insecticides in Cx. pipiens (50). The G119S mutation is common in Cx. pipiens populations in Europe, Africa and China but F290V mutation does not have a large disper- sal in Cx. pipiens populations outside the Med- iterranean region (53). Culex pipiens is a dominant species in Iran (11), especially Mazandaran Province (14, 15, 18). It is an important vector of the West Nile Virus in the world (1). Evidence of the West Nile virus circulating in the north of the coun- try (54) and the presence of wetlands for mi- grating birds (reservoir hosts), is a cause for concern for the spread of the virus in the prov- ince. In a recent study in Mazandaran Province, resistance to several insecticides has been doc- umented in Cx. pipiens (30), however, there is no report on the mechanisms of insecticide re- sistance in this species from Iran. Therefore, determination of the underlying insecticide re- sistance mechanisms of Cx. pipiens are im- portant for better insecticide resistance man- agement and implementation of effective con- trol programs if need be. Materials and Methods Study areas and sample collection The present study was undertaken in Sari Township, the capital of Mazandaran Province, northern Iran, in 2016 with geographical co- ordinates of 36◦33′48.70˝ N and 53◦03′36.35 E. It is situated between the Caspian Sea in the north and the Alborz mountain ranges in the south, Miandorud Township in the east and Qaemshahr Township in the west. The study ar- ea has a population of approximately 296,417 according to the 2011 census with an average annual temperature of 15 °C and precipitation of 789.2mm. It included two distinct geograph- ical areas of mountainous and plain/littoral. Lar- vae collection was carried out using the dip- ping sampling method from artificial habitats (drainage channels) in the village of Qajar Kheil located in a plain area of Sari Township, in summer 2015. Ecologically, the sampling site had features including stagnant water, with depth of 50cm, expanse of 5m, muddy floor, without plant and shadow–sun conditions. Water tem- perature and average pH were 20 °C and 7.12, respectively. The main occupations of the res- idents of the sampling area are agriculture, hor- ticulture, animal husbandry, fishing and hand- crafts Diazinon, chlorpyrifos, imidacloprid, cy- permethrin and glyphosate were the most com- monly used pesticides by farmers in agriculture and animal husbandry in plain areas of Sari Township especially Ghajar kheil. Rearing and preparation of samples Larvae were reared in standard insectary conditions 28±2 °C and 70±5% relative humid- ity in a 12: 12 dark: light photoperiod. Live un- fed 2–3d old adult females were put in the -80 °C freezer for conducting biochemical assays. A susceptible laboratory strain of Cx. pipiens provided by the Department of Medical Ento- mology and Vector Control, School of Public Health, Tehran University of Medical Scienc- es, Iran was used for comparison. Biochemical assays Activity levels of glutathione S-transferases (GSTs), alpha- and beta- esterases, and inhibi- tion rates of acetylcholine esterase (by propoxur), contents of cytochrome p450s and total pro- http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2019, 13(4): 378–390 SH Nikookar et al.: First Report of … 381 http://jad.tums.ac.ir Published Online: December 31, 2019 tein for adult mosquito specimens were meas- ured according to (45, 55, 56) with minor mod- ifications. Buffer solutions including 0.625M potassium phosphate pH 7.2, 0.025M sodium acetate pH 5, 1% triton sodium phosphate pH 7.8, 0.1M sodium phosphate pH 7, 0.02M so- dium phosphate pH 7.2, 5% sodium dodecyl sulphate (SDS) diluted in 0.1M sodium pho- sphate pH 7, 0.1M sodium phosphate 6.5 and 50mM sodium phosphate pH 7.4 were prepared fresh and used for a maximum of two weeks. The rest of the solutions and reagents were made fresh during the tests being done. Forty-five deep-frozen female adult mosquitoes (45 spec- imens of the field and 45 susceptible) were individually put in wells of a 96-well micro- titer plate (Maxwell®, China), initially, 100μl of cold dwater was added to each well and the specimens were homogenized by a steel pes- tle on ice. Then another 250μl of cold dwater was added and the plate was centrifuged at 1100g (3000rpm) in a refrigerated centrifuge (Beckman Coulter®, Inc., California, USA) at 4 °C for 20min. The supernatants were used as enzyme source for biochemical assays. Absorbance levels were measured spec- trophotometrically with a microplate reader (ELX808 Ultra Microplate Reader BIO- TEK® Inc., California, USA) activated by KC-Junior software, at specific wavelengths for each enzyme. The mean absorbance was calculated based on data for the two replicate wells per mosquito. Three blank replicates were considered in all enzyme assays that in- clude all reagents and working solutions re- lated to activities of each enzyme except the enzyme source (instead of which distilled wa- ter was used). Cytochrome p450s assay In a fresh microtiter plate, reaction mix- ture in each well consisting of 20 microliters of the homogenate in duplicate, 80 microliters of 0.625M potassium phosphate buffer pH 7.2, 200 microliters of 3,3′,5,5′ tetramethylbenzidine (TMBZ) solution (0.01g TMBZ dissolved in 5ml methanol plus 15ml of 0.25M sodium acetate buffer pH 5.0) and 25 microliters of 3% hydrogen peroxide were mixed. Then, plates were kept at room temperature for 2h and absorbance was read at 450nm as an end- point. The enzyme contents were described as equivalent units of cytochrome (EUC) p450s /mg protein amended for the known heme content of cytochrome C and p450s using a standard curve of purified cytochrome C. Acetylcholinesterase (AChE) assay Twenty-five microliters of homogenates were put in each well in duplicates followed by adding 145l of Triton phosphate buffer (1% Triton X-100 in 0.1M phosphate buffer pH 7.8) to each replicate to solubilize AChE. Ten l of DTNB solution (0.01M dithiobis-2- nitrobenzoic acid in 0.1M phosphate buffer pH 7.0) and 25l of the substrate ASCHI (0.01M acetylthiocholine iodide) were added to one replicate to initiate the reaction. The latter solution was replaced by 25l of the substrate ASCHI containing 0.2% of the in- hibitor propoxur (0.1M) for the second test replicates. The kinetics of the enzyme reac- tion was monitored continuously at 405nm for 5 min. The percentages of inhibition of AChE activity by propoxur in the test were calculat- ed relative to the uninhibited wells. The assay conditions were predetermined so that indi- viduals without an insensitive AChE-based resistance mechanism had >60% inhibition of the AChE activity. General esterase assay The alpha- and beta-naphthyl acetate sub- strates were used to measure the general es- terases activity. Reaction mixtures included 20 microliters of the homogenate in duplicate (for each of the alpha- and beta esterases) in adjacent microtiter plate wells (assigned alpha and beta) and 200 microliters of alpha- and be- ta-naphthyl acetate solution (120l of 30mM alpha- or beta-naphthyl acetate dissolved in 12ml 0.02M phosphate buffer pH 7.2) respec- http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2019, 13(4): 378–390 SH Nikookar et al.: First Report of … 382 http://jad.tums.ac.ir Published Online: December 31, 2019 tively. After 30min incubation period at room temperature, 50 microliters of fast blue solu- tion (0.023g fast blue dissolved in 2.25ml distilled water and 5.25ml of 5% SDS 0.1M sodium phosphate buffer pH 7) was added to each well. After another incubation period at room temperature for 5min, the absorbance was measured at 570nm as an endpoint. The re- sulting optical densities (OD) were converted to product concentration using standard curves of ODs for known concentrations of the prod- ucts alpha- and beta-naphthol, respectively. The enzyme activity was stated as µM of product formed/min/mg protein. GSTs assay The reaction mixture contains 10µl of the homogenate in duplicate plus 200µl of work- ing solution (10mM reduced glutathione dis- solved in 0.1M phosphate buffer pH 6.5 and 3mM CDNB originally dissolved in methanol). The absorbance was kinetically read at 340nm for 5min. The enzyme activity was reported as mM of conjugate produced/min/mg protein using the extinction coefficient of CDNB cor- rected for the path length of the solution in the microtiter plate well. Protein assay Protein content of mosquito homogenates was measured using Bradford method by add- ing 300 microliters of Bio-Rad reagent (pre- pared as 1: 4 dilution in ddH2O) to 10 micro- liters of supernatant in duplicates. The absorb- ance was measured at 570nm after the mixture was incubated for 5min at room temperature. Absorbance was converted into protein con- centration using a bovine serum albumin stand- ard curve obtained with the same method and reagents. Analysis of biochemical assays data Raw data of readings from the plate reader were transformed into enzyme activities/ con- tents using the Microsoft Excel version 10 and the equations for each enzyme group. Independ- ent-sample t-test using SPSS version 19 soft- ware (IBM, USA) was employed to compare the mean values of enzyme activities or con- tents between filed and laboratory susceptible populations. Enzyme ratios (ER) were com- puted by dividing the mean activities or con- tent of different enzyme groups of the field strain with those of the laboratory susceptible strain. The P<0.05 was considered as statisti- cally significant. The percentage inhibition of AChE by propoxur compared with non-inhib- ited reaction was computed and the threshold of >60% was considered as AChE insensitive to propoxur. Results Our previous study on Cx. pipiens revealed that this species is highly resistant to all four major classes of insecticides including orga- nochlorine, organophosphate, carbamate and pyrethroids, but the resistance level was lower to pyrethroids compared to other insecticides. The highest and the least mortality rates were produced by etofenprox and propoxur with a rate of 76.47% and 4.25%, respectively. The mean and standard deviation of cyto- chrome p450s contents in the field and labor- atory susceptible populations were 49×10-6± 299×10-7 and 282×10-7±235×10-7EUC cyto- chrome p450s/mg protein, respectively. The ratio of cytochrome p450s in the field popula- tion was 2.07 when compared with that of the susceptible laboratory population (Table 1). The differences between the contents of p450 in the field and laboratory susceptible popula- tions were statistically significant (P= 0.04). The mean GST activity in the field popu- lation was marginally less than that in the la- boratory susceptible population (0.119vs 0.129 mM/min/mg protein) with a ratio of 0.92 (Ta- ble 1). However, the difference of GSTs ac- tivity in the field and laboratory susceptible populations was statistically not significant (P> 0.05) at 5% level. The mean activity of alpha- and beta-naph- thyl acetate was 878×10-6 and 1×10-3μM/min/ http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2019, 13(4): 378–390 SH Nikookar et al.: First Report of … 383 http://jad.tums.ac.ir Published Online: December 31, 2019 mg protein in the field population and 236× 10-6 and 7×10-4μM/min/mg protein in the la- boratory susceptible population, respectively. The enzyme ratios in the field population were 3.72 and 1.36 for alpha- and beta-esterases when compared with those of the laboratory susceptible population (Table 1). There was a significant difference between the field and laboratory susceptible populations regarding the activity of alpha- (P< 0.001) and beta-es- terases (P= 0.001). Acetylcholinesterase biochemical assay showed that 62.82 % of the tested field pop- ulation had insensitive AChE to propoxur as the inhibition of AChE was less than the threshold of 60% (Table 2). Scatter diagram shows variations in enzyme levels in each individual in field and laborato- ry populations (Fig. 1). Fig. 1. Scatter diagram of distribution of enzyme activities/contents levels in each individual in the field and labora- tory susceptible populations of Cx. pipiens from Sari, Iran in 2016 http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2019, 13(4): 378–390 SH Nikookar et al.: First Report of … 384 http://jad.tums.ac.ir Published Online: December 31, 2019 Table 1. Mean enzyme activity/content and enzyme ratio (ER) of MFO, GST, Alpha and Beta esterases in the field and susceptible laboratory populations of Cx. pipiens, summer 2016 Enzymes Population Mean ER MFO Field Population 49 × 10-6 2.07 laboratory susceptible Population 282 × 10-7 1 GST Field Population 0.119 0.92 laboratory susceptible Population 0.129 1 Alpha Field Population 878 × 10-6 3.72 laboratory susceptible Population 236 × 10-6 1 Beta Field Population 1 × 10-3 1.36 laboratory susceptible Population 7 × 10-4 1 Table 2. Inhibition of AChE activity by propoxur in the filed population of Cx. pipiens collected in summer 2016 No of tested Cx. pipiens Inhibition rates of AChE No of tested Cx. pipiens Inhibition rates of AChE No of tested Cx. pipiens Inhibition rates of AChE No of tested Cx. pipiens Inhibition rates of AChE 1 74.92 21 65.8 41 36.5* 61 45.8* 2 52.3* 22 25.4* 42 86.7 62 60.4* 3 71.3 23 63.1 43 52.8* 63 53.1* 4 59.4* 24 41* 44 54* 64 54.9* 5 54.8* 25 13.4* 45 60.5* 65 71.2 6 45.7* 26 68.8 46 17.3* 66 57* 7 52.8* 27 36.6* 47 67 67 56.8* 8 38.1* 28 57.7* 48 59.9* 68 56.2* 9 73.1 29 58.5* 49 86.5 69 60.3* 10 76.76 30 38.5* 50 60.5* 70 82.5 11 68.4 31 78.1 51 65.3 71 65.6 12 73.4 32 49.6* 52 24.4* 72 51.9* 13 52.5* 33 37.2* 53 43.7* 73 65.2 14 63.5 34 76 54 50.3* 74 43.1* 15 32.3* 35 46.8* 55 21.7* 75 21.4* 16 76.3 36 64.7 56 59.2* 76 84.2 17 68.7 37 67.4 57 33.4* 77 16.8* 18 48.6* 38 59.5* 58 35.2* 78 64 19 84.3 39 75.9 59 10.8* 20 66.4 40 2.6* 60 45.7* *The inhabitation rate is lower than the threshold of 60% http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2019, 13(4): 378–390 SH Nikookar et al.: First Report of … 385 http://jad.tums.ac.ir Published Online: December 31, 2019 Discussion Culex pipiens is considered to be a com- petent vector of pathogens in the world (1). This species is one of the most common mos- quito species in Iran (11) and has a widespread distribution in the North of the country in the Caspian Sea region (14, 15). In a recent study in the area Cx. pipiens showed resistance to organochlorine, organo- phosphorus, carbamates and pyrethroid insec- ticides (30). Development of insecticide re- sistance can be caused by direct selection pres- sure from insecticides used in vector control programmes and or indirectly by the use of ag- rochemicals (57). In this follow up study, ele- vated enzyme activities of alpha- and beta- es- terases, cytochrome P450s and insensitive AChE in the filed population of Cx. pipiens confirmed the occurrence of resistance and also elucidat- ed its underlying resistance mechanisms. The results of this study indicated that dif- ferent enzyme groups may play a role in insec- ticide resistance in the field population of Cx. pipiens. Alpha- and beta-esterases and cyto- chrome p450s levels were significantly higher in the filed population of Cx. pipiens compared to the laboratory susceptible population. Involve- ment of esterases and cytochrome p450s in DDT, malathion and pyrethroid resistance was described in populations of Cx. pipiens from Turkey (21). Zayed et al. showed an associa- tion between elevated levels of esterases ac- tivity in Cx. pipiens larvae with organophos- phorus resistance (43) which is in accordance with the findings of our research. Enhancements of cytochrome p450s in py- rethroid resistance were reported in Cx. pipiens quinquefasciatus (58), Cx. pipiens pallens (59), An. stephensi (32), Ae. aegypti and Ae. albopic- tus (60, 61). General esterases were mostly re- ported to be involved in resistance to organo- phosphates and carbamates, elevated levels of these enzymes have also been associated with resistance to organophosphates in Cx. pipiens quinquefasciatus (58), Anopheles stephensi (35) and cross resistance to pyrethroids (34). In the present study, significant increase in alpha-esterases activity was discovered in the field population compared with other enzyme groups relative to the laboratory susceptible population. This higher level of enzyme activ- ity can possibly suggest the significance of de- velopment of organophosphorus resistance in the field populations of Cx. pipiens. The de- velopment of resistance to insecticide in the study area could possibly be due to the heavy use of pesticides. Mazandaran is the leading province in terms of pesticide use based on its toxic load in the Iran (62). In recent years farm- ers have converted fields not previously culti- vated into horticulture, and this is probably one reason for the high use of pesticides for pest control in the area Pesticides applied in agri- culture can pollute streams running off farms and this could be the indirect path of how these insecticides may end up in larval habitats of mosquitoes. Consequently high usage of insec- ticides in horticulture and agriculture practices leads to emergence of resistance to insecticides in mosquito populations in many areas (63). Similar observations have also been made in populations of Cx. quinquefasciatus from Malaysia and India (58, 64) and Ae. aegypti and Ae. albopictus from Thailand (47) associ- ated with elevated alpha-esterases levels. In this study the difference in the activities of glutathione S-transferases in susceptible la- boratory and field populations were not statis- tically significant from each other. This is ob- served in similar studies in Iran (29). GSTs are mainly involved in organochlorine insecticides resistance (44). Whereas p450s and esterases are contributing to resistance to most groups of insecticides. As the population tested in the current study was highly resistant to many dif- ferent organophosphorous (OP) and carbamates insecticides, this is quite natural, according to the economics of resistance, that the insects over- produce those enzymes with higher role in de- http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2019, 13(4): 378–390 SH Nikookar et al.: First Report of … 386 http://jad.tums.ac.ir Published Online: December 31, 2019 toxifying more insecticide groups. Acetylcholinesterase (AChE) is a key en- zyme in the nervous system, hydrolyzing ace- tylcholine neurotransmitter and terminating neu- ral impulses, and are the target for both organ- ophosphates and carbamate insecticides, two main classes of pesticides applied for pest and vector control in agriculture and public health (22). In the present study, inhibition rates of AChE by propoxur showed the frequency of insensitive AChE gene in the filed population of Cx. pipiens to be higher than the threshold of <60% which confirmed its involvement in the resistance to organophosphates and carba- mates. High insecticide resistance caused by in- sensitivity or reduced sensitivity in response to OPs and carbamate insecticides was observed in An. gambiae, An. albimanus, Cx. vishnui, Cx. pipiens, and Cx. quinquefasciatus (22). High in- sensitivity of acetylcholinesterase can be main- ly due to mutations in the ace-1gene (22). Thus, molecular studies are recommended to verify the presence of these mutations in ace-1 gene in the field populations of Cx. pipiens in Ma- zandaran Province. Conclusion The findings of this research provide the first evidence on the involvement of metabol- ic mechanisms including alpha- and beta-es- terases, mixed-function oxidases and acetyl- cholinesterase in insecticide resistance in the field population of Cx. pipiens from Sari, Ma- zandaran Province, Iran. There is evidence of detection of WNV in the mosquito populations and records of human infection with the virus in North and North West of Iran. The findings of the present study provide valuable information on the underlying resistance mechanisms of this species which is pivotal in implementing suc- cessful control program of this important vector. Acknowledgements This study was financially supported by the National Institute for the Medical Research De- velopment (NIMAD) grant number 958707. The authors declare that there is no conflict of interest. References 1. Weitzel T, Jawień P, Rydzanicz K, Lonc E, Becker N (2015) Culex pipiens sl and Culex torrentium (Culicidae) in Wrocław area (Poland): occurrence and breeding site preferences of mosquito vectors. Parasi- tol Res. 114: 289–295 2. Turell MJ (2012) Members of the Culex pipiens complex as vectors of viruses. J Am Mosq Control Assoc. 28: 123–126. 3. de Wispelaere M, Despres P, Choumet V (2017) European Aedes albopictus and Culex pipiens Are Competent Vectors for Japanese Encephalitis Virus. PLoS Negl Trop Dis. 11: e0005294. 4. Bravo-Barriga D, Parreira R, Almeida AP, Calado M, Blanco-Ciudad J, Serrano- Aguilera FJ, Pérez-Martín JE, Sánchez- Peinado J, Pinto J, Reina D, Frontera E (2016) Culex pipiens as a potential vector for transmission of Dirofilaria immitis and other unclassified Filarioidea in South- west Spain.Vet Parasitol. 223: 173–180. 5. Ventim R, Ramos JA, Osorio H, Lopes RJ, Perez-Tris J, Mendes L (2012) Avian ma- laria infections in western European mos- quitoes. Parasitol Res. 111: 637–645. 6. Shi L, Hu H, Ma K, Zhou D, Yu J, Zhong D, Fang F, Chang X, Hu Sh, Zou F (2015) Development of resistance to pyrethroid in Culex pipiens pallens population un- der different insecticide selection pres- sures. PloS Negl Trop Dis. 9: e0003928. 7. Gao Q, Xiong C, Su F, Cao H, Zhou J, Jiang Q (2016) Structure, Spatial and Tem- poral Distribution of the Culex pipiens Complex in Shanghai, China. Int J Envi- ron Res Public Health. 13: 1150. 8. Vinogradova E, Shaikevich E, Ivanitsky A (2007) A study of the distribution of the http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2019, 13(4): 378–390 SH Nikookar et al.: First Report of … 387 http://jad.tums.ac.ir Published Online: December 31, 2019 Culex pipiens complex (Insecta: Diptera: Culicidae) mosquitoes in the European part of Russia by molecular methods of iden- tification. Comp Cytogenet. 1: 129–138. 9. Hesson JC, Schäfer M, Lundström JO (2016) First report on human-biting Culex pipiens in Sweden. Parasit Vectors. 9: 632. 10. Chaulk AC, Carson KP, Whitney HG, Fon- seca DM, Chapman TW (2016) The ar- rival of the northern house mosquito Culex pipiens (Diptera: Culicidae) on New- foundland's Avalon Peninsula. J Med Entomol. 53: 1364–1369. 11. Zaim M (1987) The distribution and larval habitat characteristics of Iranian Culic- inae. Journal of the American Mosquito Control Association. 3: 568–573. 12. Moosa-Kazemi S, Vatandoost H, Nikookar H, Fathian M (2009) Culicinae (Diptera: culicidae) mosquitoes in Chabahar Coun- ty, Sistan and Baluchistan Province, south- eastern iran. Iran J Arthropod Borne Dis. 3: 29–35. 13. Dehghan H, Sadraei J, Moosa-Kazemi SH, Abolghasemi E, Solimani H, Jaffari No- doshan A, Najafi MH (2016) A pictorial key for Culex pipiens Complex (Diptera: Culicidae) In Iran. J Arthropod Borne Dis. 10: 291–302. 14. Nikookar S, Fazeli-Dinan M, Azari-Ha- midian S, Mousavinasab S, Arabi M, Zia- pour S, Shojaee J, Enayati AA (2017) Species composition and abundance of mosquito larvae in relation with their hab- itat characteristics in Mazandaran Prov- ince, northern Iran. Bull Entomol Res. 107: 598–610. 15. Nikookar SH, Azari-Hamidian S, Fazeli- Dinan M, Nasab SNM, Aarabi M, Zia- pour SP, Ziapour P, Enayati AA (2016) Species composition, co-occurrence, asso- ciation and affinity indices of mosquito larvae (Diptera: Culicidae) in Mazanda- ran Province, northern Iran. Acta Trop. 157: 20–29. 16. Nikookar SH, Moosa-Kazemi SH, Yaghoobi- Ershadi MR, Vatandoost H, Oshaghi MA, Ataei A, Anjamrooz M (2015) Fauna and larval habitat characteristics of mosqui- toes in Neka County, northern Iran. J Ar- thropod Borne Dis. 9: 253–266. 17. Nikookar SH, Moosa-Kazemi SH, Oshaghi MA, Vatandoost H, Yaghoobi-Ershadi MR, Enayati AA, Motevali-Haghi F, Zia- pour SP, Fazeli-Dinan M (2015) Biodi- versity of culicid mosquitoes in rural Neka township of Mazandaran Province, north- ern Iran. J Vector Borne Dis. 52: 63–72. 18. Nikookar SH, Moosa-Kazemi SH, Oshaghi M, Yaghoobi-Ershadi M, Vatandoost H, Kianinasab A (2010) Species composi- tion and diversity of mosquitoes in neka County, Mazandaran Province, northern iran. J Arthropod-Borne Dis. 4: 26–34. 19. McCarroll L, Hemingway J (2002) Can in- secticide resistance status affect parasite transmission in mosquitoes? Insect Bi- ochem Mol Biol. 32: 1345–1351. 20. Labbe P, Berticat C, Berthomieu A, Unal S, Bernard C, Weill M, Lenormand T (2007) Forty years of erratic insecticide resistance evolution in the mosquito Cu- lex pipiens. PLoS Genet. 3(11): e205. 21. Akıner MM, Ekşi E (2015) Evaluation of insecticide resistance and biochemical mechanisms of Culex pipiens L. in four localities of east and middle mediterra- nean basin in Turkey. Int J Mosq Res. 2: 39–44. 22. Liu N (2015) Insecticide resistance in mos- quitoes: impact, mechanisms, and research directions. Annu Rev Entomol. 60: 537– 559. 23. Scott JG, Yoshimizu MH, Kasai S. Pyre- throid resistance in Culex pipiens mos- quitoes (2015) Pestic Biochem Physiol. 120: 68–76. 24. Zhang H, Meng F, Qiao C, Cui F (2012) Identification of resistant carboxylester- ase alleles in Culex pipiens complex via PCR-RFLP. Parasit Vectors. 5: 209. 25. Zaim M, Mousavi SB, Baghestani MA, Aitio http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2019, 13(4): 378–390 SH Nikookar et al.: First Report of … 388 http://jad.tums.ac.ir Published Online: December 31, 2019 A (2017) An assessment of agricultural pesticide use in Iran, 2012–2014. J En- viron Health Sci Eng. 15: 1–8. 26. Naseri-Karimi N, Vatandoost H, Bagheri M, Chavshin AR (2015) Susceptibility status of Culex pipiens against deltame- thrin and DDT, Urmia County, West Azerbaijan Province, northwestern Iran. Asian Pac J Trop Dis. 5: S77–S79. 27. Ataie A, Moosa-Kazemi SH, Vatandoost H, Yaghoobi-Ershadi MR, Bakhshi H, Anjomruz M (2015) Assessing the sus- ceptibility status of mosquitoes (Diptera: Culicidae) in a dirofilariasis focus, north- western Iran. J Arthropod Borne Dis. 9 (1): 7–21. 28. Fathian M, Vatandoost H, Moosa-Kazemi SH, Raeisi A, Yaghoobi-Ershadi MR, Oshaghi MA, Sedaghat MM (2015) Sus- ceptibility of Culicidae mosquitoes to some insecticides recommended by WHO in a malaria endemic area of southeastern Iran. J Arthropod Borne Dis. 9(1): 22–34. 29. Salim-Abadi Y, Oshaghi MA, Enayati AA, Abai MR, Vatandoost H, Eshraghian MR, Mirhendi H, Hanafi-Bojd AA, Gorouhi MA, Rafi F (2016) High Insecticides Re- sistance in Culex pipiens (Diptera: Cu- licidae) from Tehran, Capital of Iran. J Arthropod Borne Dis. 10(4): 483–492. 30. Ghorbani F, Vatandoost H, Hanafi-Bojd AA, Abai MR, Nikookar H, Enayati AA (2018) High resistance of vector of West Nile virus, Culex pipiens Linnaeus (Dip- tera: Culicidae) to different insecticides recommended by WHO in northern Iran. J Arthropod Borne Dis. 12: 24–30. 31. Brown FV, Logan RAE, Wilding CS (2019) Carbamate resistance in a UK population of the halophilic mosquito Ochlerotatus detritus implicates selection by agricultural usage of insecticide. Int J Pest Manag. 65: 284–292. 32. Safi NHZ, Ahmadi AA, Nahzat S, Ziapour SP, Nikookar SH, Fazeli-Dinan M, Enayati AA, Hemingway J (2017) Evidence of metabolic mechanisms playing a role in multiple insecticides resistance in Anoph- eles stephensi populations from Afghan- istan. Malar J. 16: 100. 33. Ziapour SP, Kheiri S, Fazeli-Dinan M, Sahraei-Rostami F, Mohammadpour RA, Aarabi M, Nikookar SH, Sarafrazi M, Asgarian F, Enayati A (2017) Pyrethroid resistance in Iranian field populations of Rhipicephalus (Boophilus) annulatus. Pes- tic Biochem Physiol. 136: 70–79 34. Enayati AA, Vatandoost H, Ladonni H, Townson H, Hemingway J (2003) Molec- ular evidence for a kdr‐like pyrethroid re- sistance mechanism in the malaria vec- tor mosquito Anopheles stephensi. Med Vet Entomol. 17: 138–144. 35. Soltani A, Vatandoost H, Oshaghi MA, Maleki-Ravasan N, Enayati AA, Asgar- ian F (2015) Resistance mechanisms of Anopheles stephensi (Diptera: Culicidae) to Temephos. J Arthropod Borne Dis. 9: 71–83. 36. Marcombe S, Mathieu RB, Pocquet N, Riaz MA, Poupardin R, Sélior S, Darriet F, Reynaud S, Yébakima A, Corbel V (2012) Insecticide resistance in the dengue vector Aedes aegypti from Martinique: distribu- tion, mechanisms and relations with en- vironmental factors. PLoS One. 7: e30989. 37. Morou E, Dowd AJ, Rajatileka S, Steven A, Hemingway J, Ranson H, Paine M, Vontas J (2010) A simple colorimetric as- say for specific detection of glutathione-s transferase activity associated with DDT resistance in mosquitoes. PLoS Negl Trop Dis. 4: e808. 38. Gunasekaran K, Muthukumaravel S, Sahu S, Vijayakumar T, Jambulingam P (2011) Glutathione S Transferase activity in In- dian vectors of malaria: a defense mech- anism against DDT. J Med Entomol. 48: 561–569. 39. Enayati AA, Ladonni H (2006) Biochemi- cal assay baseline data of permethrin re- sistance in Anopheles stephensi (Diptera: http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2019, 13(4): 378–390 SH Nikookar et al.: First Report of … 389 http://jad.tums.ac.ir Published Online: December 31, 2019 Culicidae) from Iran. Pakistan J Biol Sci. 9: 1265–1270. 40. Brooke B, Kloke G, Hunt R, Koekemoer L, Tem E, Taylor M, Small G, Heming- way J, Coetzee M (2001) Bioassay and biochemical analyses of insecticide re- sistance in southern African Anopheles funestus (Diptera: Culicidae). Bull Ento- mol Res. 91(4): 265–272. 41. Fonseca-González I, Quiñones ML, McAl- lister J, Brogdon WG (2009) Mixed-func- tion oxidases and esterases associated with cross-r esistance between DDT and lambda-cyhalothrin in Anopheles darlingi Root 1926 populations from Colombia. Mem Inst Oswaldo Cruz. 104: 18–26. 42. Achaleke J, Martin T, Ghogomu RT, Vais- sayre M, Brévault T (2009) Esterase‐ mediated resistance to pyrethroids in field populations of Helicoverpa armigera (Lep- idoptera: Noctuidae) from Central Africa. Pest Manag Sci. 65: 1147–1154. 43. Zayed ABB, Szumlas DE, Hanafi HA, Fryauff DJ, Mostafa AA, Allam KM, Brogdon WG (2006) Use of bioassay and microplate assay to detect and measure insecticide resistance in field populations of Culex pipiens from filariasis endemic areas of Egypt. J Am Mosq Control Assoc. 22(3): 473–482. 44. Che-Mendoza A, Penilla RP, Rodríguez DA (2009) Insecticide resistance and gluta- thione S-transferases in mosquitoes: a review. Afr J Biotechnol. 8: 1386–1397. 45. Enayati AA, Asgarian F, Amouei A, Sharif M, Mortazavi H, Boujhmehrani H, Hem- ingway J (2010) Pyrethroid insecticide re- sistance in Rhipicephalus bursa (Acari, Ixodidae). Pestic Biochem Physiol. 97: 243–248. 46. Bourguet D, Capela R, Raymond M (1996) An insensitive acetylcholinesterase in Cu- lex pipiens (Diptera: Culicidae) from Por- tugal. J Econ Entomol. 89: 1060–1066. 47. Pethuan S, Jirakanjanakit N, Saengtharatip S, Chareonviriyaphap T, Kaewpa D, Rongnoparut P (2007) Biochemical stud- ies of insecticide resistance in Aedes (Steg- omyia) aegypti and Aedes (Stegomyia) al- bopictus (Diptera: Culicidae) in Thailand. Trop Biomed. 24: 7–15. 48. Emtithal AE-S, Thanaa AE-B (2012) Ef- ficacy of some insecticides on field pop- ulations of Culex pipiens (Linnaeus) from Egypt. J Basic Appl Zool. 65: 62–73. 49. Tabbabi A, Daaboub J, Laamari A, Cheikh R, Cheikh H (2017) Pirimiphos-Methyl resistance status of field populations of Culex pipiens (Diptera: Culicidae) from Grand Tunis area, northeast Tunisia. He- reditary Genet. 6(1): 175. 50. Alout H, Labbe P, Berthomieu A, Ma- koundou P, Fort P, Pasteur N, Weill M (2016) High chlorpyrifos resistance in Culex pipiens mosquitoes: strong syner- gy between resistance genes. Heredity (Edinb). 116: 224–231. 51. Bkhache M, Tmimi FZ, Charafeddine O, Faraj C, Failloux AB, Sarih Mh (2016) First report of L1014F-kdr mutation in Culex pipiens complex from Morocco. Parasit Vectors. 9: 644. 52. Xu W, Liu S, Zhang Y, Gao J, Yang M, Liu X, Tao L (2017) Cypermethrin re- sistance conferred by increased target insensitivity and metabolic detoxification in Culex pipiens pallens Coq. Pestic Bi- ochem Physiol. 142: 77–82. 53. Johnson BJ, Fonseca DM (2016) Insecti- cide resistance alleles in wetland and res- idential populations of the West Nile virus vector Culex pipiens in New Jersey. Pest Manag Sci. 72: 481–488. 54. Bagheri M, Terenius O, Oshaghi MA, Mota- zakker M, Asgari S, Dabiri F, Vatandoost H, Mohammadi Bavani M, Chavshin AR (2015) West Nile virus in mosquitoes of Iranian wetlands. Vector Borne Zoono- tic Dis. 15: 750–754. 55. Hemingway J (1998) Field and Laborato- ry Manual for the Mechanistic Detection http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2019, 13(4): 378–390 SH Nikookar et al.: First Report of … 390 http://jad.tums.ac.ir Published Online: December 31, 2019 of Insecticide Resistance in Insects. World Health Organization, Geneva. 56. Penilla PR, Rodriguez AD, Hemingway J, Torres JL, Arredondo-Jimenez Ji, Rodri- guez MH (1998) Resistance management strategies in malaria vector mosquito con- trol. Baseline data for a large‐scale field trial against Anopheles albimanus in Mex- ico. Med Vet Entomol. 12: 217–233. 57. Reid MC, McKenzie FE (2016) The con- tribution of agricultural insecticide use to increasing insecticide resistance in Afri- can malaria vectors. Malar J. 15: 107. 58. Low VL, Chen CD, Lee HL, Tan TK, Chen CF, Leong CS, Lim YAL, Lim PE, Norma-Rashid Y, Sofian-Azirun M (2013) Enzymatic characterization of insecticide resistance mechanisms in field popula- tions of Malaysian Culex quinquefascia- tus Say (Diptera: Culicidae). PloS one. 8: e79928. 59. Shen B, Dong HQ, Tian H-S, Ma L, Li XL, Wu GL, Zhu CL (2003) Cytochrome p450 genes expressed in the deltamethrin- susceptible and-resistant strains of Culex pipiens pallens. Pestic Biochem Physiol. 75: 19–26. 60. Wan-Norafikah O, Nazni WA, Lee HL, Zainol-Ariffin P, Sofian-Azirun M (2010) Permethrin resistance in Aedes aegypti (Linnaeus) collected from Kuala Lum- pur, Malaysia. J Asia Pac Entomol. 13: 175–182. 61. Wan-Norafikah O, Nazni WA, Lee HL, Zainol-Ariffin P, Sofian-Azirun M (2013) Susceptibility of Aedes albopictus Skuse (Diptera: Culicidae) to permethrin in Kua- la Lumpiur, Malaysia. Asian Biomed. 7: 51–62. 62. Vatandoost H, Zaim M, Hanafi-Bojd AA, Mousavi B, Nikpour F (2015) Supply of chemical pesticides in agricultural in Iran (2015–2016). institute for environmental research, Tehran university of Medical Sciences. 63. Yadouleton AWM, Asidi A, Djouaka RF, Braïma J, Agossou CD, Akogbeto MC (2009) Development of vegetable farm- ing: a cause of the emergence of insecti- cide resistance in populations of Anoph- eles gambiae in urban areas of Benin. Malar J. 8: 103. 64. Sarkar M, Bhattacharyya I, Borkotoki A, Goswami D, Rabha B, Baruah I, Srivasta- va R (2009) Insecticide resistance and detoxifying enzyme activity in the prin- cipal bancroftian filariasis vector, Culex quinquefasciatus, in northeastern India. Med Vet Entomol. 23: 122–131. http://jad.tums.ac.ir/