- Open Access
Correlation between carboxylesterase alleles and insecticide resistance in Culex pipiens complex from China
Parasites & Vectorsvolume 4, Article number: 236 (2011)
In China, large amounts of chemical insecticides are applied in fields or indoors every year, directly or indirectly bringing selection pressure on vector mosquitoes. Culex pipiens complex has evolved to be resistant to all types of chemical insecticides, especially organophosphates, through carboxylesterases. Six resistant carboxylesterase alleles (Ester) were recorded previously and sometimes co-existed in one field population, representing a complex situation for the evolution of Ester genes.
In order to explore the evolutionary scenario, we analyzed the data from an historical record in 2003 and a recent investigation on five Culex pipiens pallens populations sampled from north China in 2010. Insecticide bioassays showed that these five populations had high resistance to pyrethroids, medium resistance to organophosphates, and low resistance to carbamates. Six types of Ester alleles, EsterB1, Ester2, Ester8, Ester9, EsterB10, and Ester11 were identified, and the overall pattern of their frequencies in geographic distribution was consistent with the report seven years prior to this study. Statistical correlation analysis indicated that Ester8 and Ester9 positively correlated with resistance to four insecticides, and EsterB10 to one insecticide. The occurrences of these three alleles were positively correlated, while the occurrence of EsterB1 was negatively correlated with Ester8, indicating an allelic competition.
Our analysis suggests that one insecticide can select multiple Ester alleles and one Ester allele can work on multiple insecticides. The evolutionary scenario of carboxylesterases under insecticide selection is possibly "one to many".
Mosquitoes, due to their special behavior, physiology and close relationship with humans, act as ideal transmitters of a wide variety of human disease agents, including filariasis, Japanese encephalitis, West Nile virus, dengue and malaria . Chemical insecticides have been extensively used since the 1940s to control mosquito vectors, and four major categories of insecticides have been sequentially applied historically: organochlorines, organophosphates, carbamates and pyrethroids [2–4]. In China, insecticide production is over 3.9 × 108 kg and application is around 2.5 × 108 kg annually since 1997; directly or indirectly bringing heavy selection pressure on mosquitoes . As one of seven major species of vector mosquitoes in China, the Culex pipiens complex has evolved to be resistant to all types of chemical insecticides, except carbamates, in many regions .
Three gene families of detoxification enzymes are involved in metabolic resistance to chemical insecticides: carboxylesterases, P450 monooxygenases, and glutathione S-transferases . Carboxylesterases often mediate resistance to organophosphates, carbamates, and to a lesser extent, pyrethroids . In mosquitoes, non-specific carboxylesterases sequester, rather than rapidly metabolize, the pesticide molecules . Two gene loci, Est-2 and Est-3, encode non-specific esterases through gene amplification or up-regulation in the C. pipiens complex [8, 9]. These loci are always in complete linkage disequilibrium and referred to as the Ester superlocus. To date, twelve Ester alleles conferring OP resistance have been identified at the Ester locus (the corresponding overproduced esterases are named in parentheses): Ester1 (A1), Ester2 (A2-B2), Ester4 (A4-B4), Ester5 (A5-B5), Ester8 (A8-B8), Ester9 (A9-B9), Ester11 (A11-B11), EsterB1 (B1), EsterB6 (B6), EsterB7 (B7), EsterB10 (B10) and EsterB12 (B12) [5, 9–11].
An unusual diversity of Ester alleles was observed in field populations from China, where EsterB1, Ester2, EsterB6, EsterB7, Ester8, Ester9, EsterB10 and Ester11 were reported, and as many as six Ester alleles co-existed in one population . This situation is atypical, and may represent a complex situation for the evolution of insecticide resistance genes in China. Accumulation and analysis of monitoring data will be helpful to understand the evolution of insecticide resistance genes in field mosquito populations.
In recent years, we monitored the resistance levels and Ester phenotype frequencies of C. pipiens complex field populations to organophosphates, carbamates, and pyrethroids at either the large scale or locally in China . In order to explore the evolutionary scenario of Ester resistance alleles, we analyzed the data from five Culex pipiens pallens populations sampled in north China, and from an historical record from 2003 . We investigated the resistance levels of these five populations towards three organophosphate, two pyrethroid, and two carbamate insecticides. Esterase phenotype frequencies of Ester locus were also quantified. Lastly, we combined our data with that of a large scale investigation from 2003 , and statistically analyzed the relationship between resistance level and Ester phenotype frequency for each insecticide, and among phenotype frequencies of various Ester alleles.
Mosquito samples and strains
Five field populations of C. p. pallens were collected as egg-rafts, larvae, or pupae in north China from July to August of 2010 (Table 1). They were raised for several generations in the laboratory at 25 ± 1° and a photoperiod of 14L : 10D. Several standard strains were used as references in bioassays or starch gel electrophoresis: S-LAB, an insecticide-susceptible strain without any known resistance genes ; SB1 strain, homozygous for EsterB1, displaying overproduced esterase B1; SA2 strain, homozygous for Ester2, displaying overproduced esterases A2-B2 ; MAO2 strain, homozygous for Ester8, displaying overproduced esterases A8-B8 ; LING strain, homozygous for Ester9, displaying overproduced esterases A9-B9 ; KARA2 strain, homozygous for EsterB10, displaying overproduced esterase B10 ; and WU strain, homozygous for Ester11, displaying overproduced esterases A11-B11 .
Resistance characteristics of larvae were determined by bioassays on fourth-instar larvae, following the methods described in Raymond and Marquine (1994) . The test insecticides were three organophosphates, dichlorvos, fenitrothion, malathion, two pyrethroids, deltamethrin, permethrin, and two carbamates, propoxur, BPMC (2-sec-Butylphenyl methyl carbamate). Four or five doses and four replicates (20 larvae per replicate) per dose were performed with each insecticide. The mortalities of larvae were recorded after 24 hours of treatment. Larvae from S-LAB strain were tested at the same time as susceptible controls. Based on Finney (1971)  data were analyzed by the log-probit program of Raymond (1993) , which provides LCs, slope for each mortality line, parallelism between different mortality lines and resistance ratios with 95% confidence intervals.
Starch gel electrophoresis
Ester allele phenotypes of individual adult mosquitoes were revealed using starch gel electrophoresis (TME 7.4 buffer system) as described by Pasteur et al. (1981; 1988) [20, 21]. Mosquitoes from six standard strains (SB1, SA2, MAO2, LING, KARA2 and WU) were used as markers in electrophoresis to indicate the esterase phenotypes of field collected mosquitoes.
Statistical analysis of correlation
The data of 5 populations from this study and 20 populations from the investigation in 2003  were combined for correlation analysis between resistance levels of insecticides (dichlorvos, parathion, chlorpyrifos, propoxur and BPMC) and Ester phenotype frequencies and among phenotype frequencies of various Ester alleles using SPSS 13.0.
Insecticide resistance status of field populations
Five field populations of C. p. pallens, collected in north China (Table 1), were tested for resistance levels towards frequently applied organophosphate, pyrethroid, and carbamate insecticides with larvae bioassays. The tests were finished within three generations (F1-F3) after the mosquitoes were brought back and raised in the lab. Bioassay results indicate that the five populations showed highest resistance levels towards pyrethroids, followed by organophosphates, and last to carbamates (Table 2). For the two pyrethroid insecticides, resistance to permethrin was higher than to deltamethrin. Most of the populations had a resistance ratio (RR) higher than 10 to permethrin, with the highest reaching 68 folds in the population TAA, while lower than 10 to deltamethrin. For the three organophosphate insecticides, several populations showed medium resistance (RR between 10 and 20) to dichlorvos, while for fenitrothion and malathion only lower than 4-fold resistance was observed. All the populations had low resistance (RR lower than 10) to the two carbamate insecticides, propoxur and BPMC.
Populations from Henan (WHE and YUY) and Shandong (TAA) provinces were more resistant to insecticides than those from Beijing (BJI) and Liaoning (LPU) province. BJI and LPU showed low levels of resistance to all tested insecticides. Populations WHE, YUY, and TAA had medium resistance to dichlorvos, and high resistance (RR higher than 20) to permethrin.
Identification of esterase phenotypes
In order to detect the types of high active carboxylesterases, at least thirty individual adult mosquitoes of each population were applied to starch gel electrophoresis. It has been reported that these high active carboxylesterases are expressed in both adults and larvae [22, 23]. A total of 200 adult mosquitoes were checked for the five populations (Table 3). Six types of overproduced esterases were identified in this investigation, i.e. B1, A2-B2, A8-B8, A9-B9, B10, and A11-B11. In the populations from Henan province (WHE and YUY), all six types of esterases were detected, and four esterases co-existed in one population, WHE. In the other three populations, only B1 and A2-B2 were detected (Figure 1). B1 was the predominant esterase, with frequencies over 70 % in most populations. Although the frequencies of A2-B2 were not as high as B1, it still had a distribution as broad as B1. Esterases A8-B8, A9-B9, B10 and A11-B11 were minor, with limited distribution and low frequencies.
Correlation analyses between insecticide resistance levels (dichlorvos, parathion, chlorpyrifos, propoxur and BPMC) and Ester phenotype frequencies showed that the frequencies of esterases A8-B8 and A9-B9 positively correlated to the resistance levels of dichlorvos, parathion, chlorpyrifos and propoxur. Esterase B10 only positively correlated to dichlorvos resistance, even though the correlation coefficients (R2) were not high (Figure 2). None of the six types of esterases displayed a linear correlation to the BPMC resistance in mosquitoes. Frequencies of esterases B1, A2-B2, and A11-B11 did not linearly correlate to any resistance of the five insecticides. For each insecticide, except for BPMC, multiple types of esterases had a positive relationship with its resistance in mosquitoes.
The correlation among phenotype frequencies of various Ester alleles was also checked. The results demonstrate that there were positive correlations between phenotype frequencies of A8-B8 and A9-B9, A8-B8 and B10, A9-B9 and B10, and a negative correlation between A8-B8 and B1 (Figure 3). The existence of A2-B2 or A11-B11 had no linear relationship with other esterases.
As a periodic investigation, we surveyed the resistance levels of five field populations of C. p. pallens towards seven insecticides commonly used in mosquito control in China. High and medium resistance appeared to pyrethroids and organophosphates, in some populations, respectively, and only low resistance was observed for carbamates in the five populations. Compared with previous investigations in 2003  and 2006 , there is no remarkable change in the organophosphate and carbamate resistance status in the 2010 survey. The resistance to dichlorvos is still leading in organophosphate resistance. Low resistance to the two carbamate insecticides implies that there is still no target resistance mediated by G119S mutation in the acetylcholinesterase 1. On the other hand, the pyrethroid resistance seems elevated, especially the resistance to permethrin. Most of the populations in this survey had the permethrin resistance ratio higher than 10, including the Beijing population BJI, while the other four Beijing populations surveyed in 2006 were susceptible or only showed resistance below 2 fold towards permethrin and tetramethrin . This is the aftermath of an increase of pyrethroid insecticides applied in recent years in China.
In this study we observed six highly active carboxylesterases in these northern populations. The overall pattern of their frequencies in geographic distribution is consistent with the reports seven years prior to this study [5, 12]. Henan province still has impressive diversity of Ester alleles. Shandong and Liaoning provinces also only have EsterB1 and Ester2. This means the Ester alleles did not extend their geographic range during the seven years in these areas. In the Beijing population, where four Ester alleles were recorded previously, only EsterB1 and Ester2 were observed. One possible reason for this difference could be that our sample size in the Beijing region was not large enough. More attention should be given to this region in future.
Among the six types of Ester alleles, Ester8 and Ester9 have linear correlation with resistance to four organophosphate or carbamate insecticides and EsterB10 to one organophosphate insecticide in field populations. These three alleles are adaptive to the selection pressure endemic to China, and the occurrences of these alleles are positively correlated at a certain coefficient. In contrast, two ubiquitously distributed alleles around the world, EsterB1 and Ester2, do not show linear correlation with resistance to any one of the five insecticides, implying that they play a general role in organophosphate and carbamate resistance. Although the endemic allele Ester11 did not linearly correlate with resistance to the five insecticides, this can not rule out the possible correlation of Ester11 with resistance to malathion, fenitrothion, or other insecticides used in China.
Our results suggest that one insecticide can select multiple Ester alleles and one Ester allele can work on multiple insecticides. So the evolutionary scenario of carboxylesterases under insecticide selection in the field is most likely "one to many", not "one to one". The relationship of EsterB1 and Ester8 is very interesting. Their occurrence seems to contradict each other. This allelic competition was illustrated by the situation in southern France, where Ester1 had been replaced by Ester4 over a 10 year period , and then followed by the local occurrence of Ester2 . It is likely that, in the future, as the result of the allelic competition, EsterB1 will be eliminated from regions predominated by Ester8, such as the Guangdong province, and Ester8 will be prevented from invading, or eliminated from, regions predominated by EsterB1 such as Beijing, and the Shandong province. Future work is needed to identify the parameters driving the competition between these two alleles.
Our analysis suggests that one insecticide can select multiple Ester alleles and one Ester allele can work on multiple insecticides. The evolutionary scenario of carboxylesterases under insecticide selection is possibly "one to many". This study will shed light on the understanding of the evolution of insecticide resistance genes in field populations and guide the management of insecticide resistance.
Hemingway J, Ranson H: Insecticide resistance in insect vectors of human disease. Annu Rev Entomol. 2000, 45: 371-391.
Liu N, Xu Q, Zhu F, Zhang L: Pyrethroid resistance in mosquitoes. Insect Sci. 2006, 13: 159-166.
Kamgang B, Marcombe S, Chandre F, Nchoutpouen E, Nwane P, Etang J, Corbel V, Paupy C: Insecticide susceptibility of Aedes aegypti and Aedes albopictus in Central Africa. Parasit Vectors. 2011, 4: 79-
Lima EP, Paiva MH, de Araújo AP, da Silva EV, da Silva UM, de Oliveira LN, Santana AE, Barbosa CN, de Paiva Neto CC, Goulart MO, Wilding CS, Ayres CF, de Melo Santos MA: Insecticide resistance in Aedes aegypti populations from Ceará, Brazil. Parasit Vectors. 2011, 4: 5-
Cui F, Tan Y, Qiao CL: Filariasis vector in China: Insecticide resistance and population structure of mosquito Culex pipiens complex. Pest Manag Sci. 2007, 63: 453-458.
Hemingway J, Hawkes NJ, McCarroll L, Ranson H: The molecular basis of insecticide resistance in mosquitoes. Insect Biochem Mol Biol. 2004, 34: 653-665.
Kadous AA, Ghiasuddin SM, Matsumura F, Scott JG, Tanaka K: Difference in the picrotoxinin receptor between the cyclodiene-resistant and susceptible strains of the German cockroach. Pestic Biochem Physiol. 1983, 19: 157-166.
Rooker S, Guillemaud T, Berge J, Pasteur N, Raymond M: Coamplification of esterase A and B genes as a single unit in Culex pipiens mosquitoes. Heredity. 1996, 77: 555-561.
Raymond M, Chevillon C, Guillemaud T, Lenormand T, Pasteur N: An overview of the evolution of overproduced esterases in the mosquito Culex pipiens. Philos Trans R Soc Lond B Biol Sci. 1998, 353: 1707-1711.
Raymond M, Berticat C, Weill M, Pasteur N, Chevillon C: Insecticide resistance in the mosquito Culex pipiens: what have we learned about adaptation?. Genetica. 2001, 112-113: 287-296.
Cheikh RB, Berticat C, Berthomieu A, Pasteur N, Cheikh HB, Weill M: Characterization of a novel high-activity esterase in Tunisian populations of the mosquito Culex pipiens. J Econ Entomol. 2008, 101: 484-491.
Cui F, Lin LF, Qiao CL, Xu Y, Marquine M, Weill M, Raymond M: Insecticide resistance in Chinese populations of the Culex pipiens complex through esterase overproduction. Entomol Exp Appl. 2006, 120: 211-220.
Yan S, Wu Z, Cui F, Zhao Q, Qiao C: Dynamics of esterase alleles in Culex pipiens complex mosquitoes in Beijing. J Econ Entomol. 2008, 101: 1897-1902.
Georghiou GP, Metcalf RL, Gidden FE: Carbamate-resistance in mosquitoes: Selection of Culex pipiens fatigans Wiedemann. (= C. quinquefasciatus Say) for resistance to Baygon. Bull World Health Organ. 1966, 35: 691-708.
Berticat C, Rousset F, Raymond M, Berthomieu A, Weill M: High Wolbachia density in insecticide-resistant mosquitoes. Proc Biol Sci. 2002, 269: 1413-1416.
Weill M, Marquine M, Berthomieu A, Dubois MP, Bernard C, Qiao CL, Raymond M: Identification and characterization of novel organophosphate detoxifying esterase alleles in the Guangzhou area of China. J Am Mosq Control Assoc. 2001, 17: 238-244.
Raymond M, Marquine M: Evolution of insecticide resistance in Culex pipiens populations: The Corsican paradox. J Evol Biol. 1994, 7: 315-337.
Finney DJ: Probit Analysis. 1971, Cambridge: Cambridge University Press
Raymond M: PROBIT CNRS-UMII. Licence L93019, Avenix, 24680 St. Georges d'Orques, France. 1993
Pasteur N, Iseki A, Georghiou GP: Genetic and biochemical studies of the highly active esterases A' and B associated with organophosphate resistance in mosquitoes of the Culex pipiens complex. Biochem Genet. 1981, 19: 909-919.
Pasteur N, Pasteur G, Bonhomme F, Britton-Davidian J: Practical Isozyme Genetics. 1988, Chichester: Ellis Horwood Press
Callaghan A, Boiroux V, Raymond M, Pasteur N: Prevention of changes in the electrophoretic mobility of overproduced esterases from organophosphate-resistant mosquitoes of the Culex pipiens complex. Med Vet Entomol. 1994, 8: 391-394.
Pasteur N, Nancé E, Bons N: Tissue Localization of Overproduced Esterases in the mosquito Culex pipiens (Diptera: Culicidae). J Med Entomol. 2001, 38: 791-801.
Guillemaud T, Lenormand T, Bourguet D, Chevillon C, Pasteur N, Raymond M: Evolution of resistance in Culex pipiens: allele replacement and changing environment. Evolution. 1998, 52: 443-453.
Labbé P, Lenormand T, Raymond M: On the worldwide spread of one insecticide resistance gene: historical chance, genetic drift or competitive advantage. J Evol Biol. 2005, 18: 1471-1484.
We are very grateful to Montana Atwater for providing constructive comments on the manuscript. This work is supported by the grant from the Major State Basic Research Development Program of China (973 Program) (No. 2012CB114102) and the Important National Science & Technology Specific Project (No. 2012ZX10004219-003).
The authors declare that they have no competing interests.
YL and HZ performed the study, analyzed the data and drafted the manuscript. CQ supervised the study and helped draft the manuscript. XL and FC conceived and designed the study, and helped draft the manuscript. All authors approved the final version of the manuscript.