Esterase polymorphism in insecticide susceptible populations of the mosquito Culex pipiens

Parallel Evolution in Mosquito Vectors-A Duplicated Esterase Locus is Associated With Resistance to Pirimiphos-methyl in Anopheles gambiae

The primary control methods for the African malaria mosquito, Anopheles gambiae, are based on insecticidal interventions. Emerging resistance to these compounds is therefore of major concern to malaria control programs. The organophosphate (OP), pirimiphos-methyl, is a relatively new chemical in the vector control armory but is now widely used in indoor-residual spray campaigns. While generally effective, phenotypic resistance has developed in some areas in malaria vectors. Here, we used a population genomic approach to identify novel mechanisms of resistance to pirimiphos-methyl in A. gambiae s.l mosquitoes. In multiple populations, we found large and repeated signals of selection at a locus containing a cluster of detoxification enzymes, some of whose orthologs are known to confer resistance to OPs in Culex pipiens. Close examination revealed a pair of alpha-esterases, Coeae1f and Coeae2f, and a complex and diverse pattern of haplotypes under selection in A. gambiae, A. coluzzii and A. arabiensis. As in C. pipiens, copy number variants have arisen at this locus. We used diplotype clustering to examine whether these signals arise from parallel evolution or adaptive introgression. Using whole-genome sequenced phenotyped samples, we found that in West Africa, a copy number variant in A. gambiae is associated with resistance to pirimiphos-methyl. Overall, we demonstrate a striking example of contemporary parallel evolution which has important implications for malaria control programs.

Parallel evolution in mosquito vectors - a duplicated esterase locus is associated with resistance to pirimiphos-methyl in An. gambiae

The primary control methods for the African malaria mosquito, Anopheles gambiae, are based on insecticidal interventions. Emerging resistance to these compounds is therefore of major concern to malaria control programmes. The organophosphate, pirimiphos-methyl, is a relatively new chemical in the vector control armoury but is now widely used in indoor residual spray campaigns. Whilst generally effective, phenotypic resistance has developed in some areas in malaria vectors. Here, we used a population genomic approach to identify novel mechanisms of resistance to pirimiphos-methyl in Anopheles gambiae s.l mosquitoes. In multiple populations, we found large and repeated signals of selection at a locus containing a cluster of detoxification enzymes, some of whose orthologs are known to confer resistance to organophosphates in Culex pipiens. Close examination revealed a pair of alpha-esterases, Coeae1f and Coeae2f, and a complex and diverse pattern of haplotypes under selection in An. gambiae, An. coluzzii and An. arabiensis. As in Cx. pipiens, copy number variation seems to play a role in the evolution of insecticide resistance at this locus. We used diplotype clustering to examine whether these signals arise from parallel evolution or adaptive introgression. Using whole-genome sequenced phenotyped samples, we found that in West Africa, a copy number variant in Anopheles gambiae is associated with resistance to pirimiphos-methyl. Overall, we demonstrate a striking example of contemporary parallel evolution which has important implications for malaria control programmes.

Isoenzymes and protein polymorphism in Blaps polycresta and Trachyderma hispida (Forsskål, 1775) (Coleoptera: Tenebrionidae) as biomarkers for ceramic industrial pollution

The expression levels of Esterase Isoenzyme and total soluble protein fractionation were studied in two coleopteran insects Blaps polycresta and Trachyderma hispida to evaluate the possible hazards from ceramic and plastic factories in the Khorshed Region, East of Alexandria, Egypt. Two insect collection sites were selected. The first site was the garden of the Faculty of Science, Moharram Bek, Alexandria University, which is considered a non-polluted site, and Khorshed district, considered as the polluted site. Percentages of heavy metals were estimated using SEM-X-ray microanalysis of soft tissues of both sexes of the two coleopteran insects. Esterase Isoenzymes were found to be overexpressed in B. Polycresta but not T. hispida. Female B. polycresta from the polluted site exhibited overexpression of the second and third loci. Furthermore, the females were found to be more affected than males, which only showed the overexpression of the second loci. T. hispida (females and males) collected from the reference site were found to have increased esterase activity compared with those sampled from the polluted site. The Snake-Skin™ Dialysis tubing technique, used for optimizing the protein extraction method, reflected the highest quantified proteins compared to other, traditional methods. SDS polyacrylamide gel electrophoresis (SDS-PAGE) analysis of the whole-body protein reflected definite variations between T. hispida and B. polycresta in fraction number and activity at the two sites. Varied expression levels for metallothionein (MT) heavy metal resistance proteins for B. polycresta and T. hispida were also detected in the study. Based on these results, we suggest that biochemical biomarkers could infer environmental hazards, B. polycresta and T. hispida are successful biomarkers for heavy metal pollution.

EVOLUTION OF RESISTANCE IN CULEX PIPIENS: ALLELE REPLACEMENT AND CHANGING ENVIRONMENT

Fixation of adaptive mutations in populations is often constrained by pleiotropic fitness costs. The evolutionary pathways that compensate such fitness disadvantages are either the occurrence of modifier genes or replacement of the adaptive allele by less costly ones. In this context, 23 years of evolution of insecticide resistance genes in the mosquito Culex pipiens from southern France are analyzed. The aim of this study is to answer the following points. Is there a fitness cost associated with these resistance genes in natural populations? Does evolution proceed through allele replacement or through selection of modifiers? And finally, how do environmental changes affect the evolution of resistance genes? Samples from the same transect, crossing the boundary between an insecticide-treated and a nontreated area, are analyzed. Clinal analyses indicate a variable fitness cost among the resistance genes and show that allele replacement has been the primary mechanism of resistance evolution in this area. It is also shown that replacement was probably due to environmental changes corresponding to modification in pesticide-treatment intensity.

Genetic Characterization of Esterase Activity Variants Associated with an Esterase Gene Amplification in a Strain of Culex pipiens from California

In the Culex pipiens complex, a common mechanism of insecticide resistance is amplification of esterase genes leading to overproduction of detoxifying esterase enzymes. A number of electrophoretic esterase alleles have been identified, and in field populations individuals with the same esterase electromorph can exhibit a wide range of esterase enzyme activities. We isolated and characterized esterase activity variants associated with the esterase B1 electromorph from a field strain. A mating scheme was used to isolate chromosomes with esterase genes from the strain into 45 families. Twenty-six of the families received esterase genes from the field strain that conferred elevated esterase activity. Mean esterase activities in these families ranged from 43 to 695 nmoles α-naphthyl acetate hydrolyzed/min/mg protein. Variance components indicated that genetic variance (i.e., genetic differences among families) accounted for 77% of the total variation in esterase activity. A comparison of mean esterase activities indicated that there were at least 11 different esterase activity variants contributing to the observed genetic variation in esterase activity among the 26 families. The relevance of these results to understanding the dynamics of amplified esterase genes in populations is discussed.

Isolation and spatial expression analysis of a partial esterase gene of Culex nigripalpus (Diptera: Culicidae)

Mosquitoes play an important role as vectors for a variety of pathogens. Insecticide use for control means results in selection of tolerant/resistant mosquitoes. An esterase cDNA was isolated from Culex nigripalpus Theobald (Diptera: Culicidae), and its spatial expression was characterized. GenBank tblastx and VectorBase searches with the translation product of the longest cDNA fragment (397 bp, CN Temsha est-1) of the Cx. nigripalpus esterase revealed similarity (95%) to the estalpha2 esterase gene product of Culex pipiens quinquefasciatus Say, suggesting it is a member of the esterase family. Transcripts encoding Cx. nigripalpus esterase were detected in midgut tissue, thoraces, and abdomens of female Cx. nigripalpus, findings that may support possible roles in feeding, digestion, or reproduction. Gene expression studies to determine the role of Cx. nigripalpus esterase in the mosquito are being conducted.

Company amplification of estalpha2/estbeta2 and correlation between esterase gene copy number and resistance to insecticides in the field Culex pipiens pallens strains collected from Beijing, China

Seven field strains of Culex pipiens pallens were evaluated in the summer of 2004 for their susceptibility to three insecticides. Larval bioassays were carried out in accordance with WHO mosquito susceptibility tests. Field samples showed resistance to dichlorvos and propoxur. Low levels of resistance to phoxim were detected compared with the laboratory-susceptible strain. Quantitative analysis showed that co-amplification of estalpha2 and estbeta2 occurred in the seven field Cx. pipiens pallens. However, the levels of transcription of estalpha2/estbeta2 genes in the field strains indicated a 15.9:1 ratio of estbeta32 cDNA over estalpha2 cDNA. Such observed differential transcription of the estalpha2/estbeta2 genes showed that the promoters regulating their expression are independent and of different strengths. Regression analysis of expression level of estalpha2/estbeta2 gene and LC50 value to insecticides showed a significant correlation between expression level ofestalpha2 gene and LC50 value to dichlorvos. However, correlation of estalpha2 gene and LC50 value to propoxur and phoxim, or estbeta2 gene and three kinds of insecticides, was not significant.

Genes conferring resistance to organophosphorus insecticides in Culex pipiens (Diptera: Culicidae) from Tunisia

In Tunisia, the mosquito Culex pipiens shows various organophosphate resistance alleles at Ester and ace-1 loci. The characterization and the distribution pattern of these alleles were studied among 20 populations sampled from north to center of Tunisia. At the Ester locus, Ester4, Ester5, and Ester(B12) were present. A new esterase characterized by the same electrophoretic migration as esterase A1 was identified: A13, encoded by Ester(A13) allele. At the ace-1 locus, the presence of the ace-1(R), ace-1(D), and F290V mutated alleles was also detected. A large heterogeneity in allelic frequencies at Ester and ace-1 loci was observed among samples, with a high significant genotypic differentiation considering both loci (F, = 0.077, P < 10(-5)), depicting variations of insecticide treatment intensity between areas. A comparison between populations collected in 1996 and 2005 showed an absence of significant resistance evolution. However, the high frequencies of resistance alleles in 2005 populations suggested that the selection pressures are still important in Tunisia. Strategies for resistance management are discussed in the context of the current knowledge of the Tunisian situation.

Characterization of novel esterases in insecticide-resistant mosquitoes

In the mosquito Culex pipiens complex (Diptera: Culicidae), the amplification of carboxylesterase genes is an important mechanism providing resistance to organophosphate insecticides. Various amplified alleles at the Ester locus have been identified over the world. In this study, two newly detected Ester alleles, Ester(B10) and Ester(11) (including associated Ester(A11) and Ester(B11)), coding for esterases B10 and A11-B11, respectively, are characterized qualitatively and quantitatively. A high molecular identity is observed both at the nucleotide level and at the deduced amino acid level among the known Ester alleles. Real-time quantitative PCR results suggest 2.5-fold amplification of the Ester(B10) allele, 36.5-fold amplification of the Ester(A11) allele, and 19.1-fold amplification of the Ester(B11) allele. The ca. 2-fold difference in amplification level between Ester(A11) and Ester(B11) may indicate a new model for the esterase amplification. Bioassays show that these two resistant Ester alleles only can confer moderate or low resistance to the tested organophosphate insecticides.

Transposable element polymorphism of Wolbachia in the mosquito Culex pipiens: evidence of genetic diversity, superinfection and recombination

Wolbachia is a group of maternally inherited endosymbiotic bacteria that infect and induce cytoplasmic incompatibility (CI) in a wide range of arthropods. In contrast to other species, the mosquito Culex pipiens displays an extremely high number of CI types suggesting differential infection by multiple Wolbachia strains. Attempts so far failed to detect Wolbachia polymorphism that might explain this high level of CI diversity found in C. pipiens populations. Here, we establish that Wolbachia infection is near to or at fixation in worldwide populations of the C. pipiens complex. Wolbachia polymorphism was addressed by sequence analysis of the Tr1 gene, a unique transposable element of the IS5 family, which allowed the identification of five C. pipiens Wolbachia strains, differing either by nucleotide substitution, presence or absence pattern, or insertion site. Sequence analysis also showed that recombination, transposition and superinfection occurred at very low frequencies. Analysis of the geographical distributions of each Wolbachia strain among C. pipiens populations indicated a strong worldwide differentiation independent from mosquito subspecies type, except in the UK. The availability of this polymorphic marker now opens the way to investigate evolution of Wolbachia populations and CI dynamics, in particular in regions where multiple crossing types coexist among C. pipiens populations.

Molecular comparisons of the Culex pipiens (L.) complex esterase gene amplicons

The amplification of carboxylesterase genes is a mechanism of organophosphate resistance in Culex mosquitoes. Amplified carboxylesterase genes from an insecticide resistant Culex pipiens strain collected in Cyprus were analysed and compared to other Culex amplified carboxylesterase alleles. A 12 kb section of genomic DNA containing two gene loci coding for carboxylesterase alleles A5 and B5 was cloned and sequenced. A comparison between this amplicon and one from a strain with co-amplified carboxylesterase alleles A2 and B2 revealed a number of differences. The intergenic spacer was 3.7 kb in length in the A5-B5 amplicon (2.7 kb in A2-B2) and contained putative Juan and transposable elements upstream of B5. A fragment of a gene with high homology to aldehyde oxidase was also present immediately downstream of A5. The comparison revealed no differences that would explain the successful spread of the A2-B2 amplicon worldwide whilst the A5-B5 amplicon is restricted to the Mediterranean.

Structural organization of the estalpha3(1) gene in a Colombian strain of Culex quinquefasciatus differs from that in Cuba

In Culex mosquitoes (Diptera: Culicidae), the most common mechanism for resistance to organophosphorus (OP) insecticides involves amplification of one or more esterases. Two esterase loci are often involved, with different allelic forms co-amplified. Estalpha3(1) is co-amplified with estbeta1 in a Colombian (COL) strain of Culex quinquefasciatus Say. These two alleles co-migrate on acrylamide gels, often leading to misscoring of the phenotype as elevation of a single estbeta enzyme. By sequencing COL genomic DNA, we determined the estalpha3(1) gene length is 1623 nucleotides. The open reading frame of estalpha3(1) encodes a 540 amino acid protein, as for estalpha2(1) in strain Pel RR from Sri Lanka. The intron/exon boundaries of estalpha3(1) are identical to those of estalpha2(1), suggesting that they are alleles of the same locus. The COL estalpha3(1) gene differs from estaalpha3(2) in strain MRES from Cuba, although they have equivalent electrophoretic mobility, showing that these two strains contain distinct resistance-associated amplicons. Twenty nucleotide differences were scored between the MRES partial 495 bp sequence and that in the COL strain, with two amino acid changes, demonstrating distinct estalpha enzymes. Our sequencing data show 95% identity between the three estalpha genes (each has six introns and seven exons) in OP-resistant Cx. quinquefasciatus. Amplified estalpha3(1) and estbeta1 are at least 10kb apart in temephos-selected COL and 2.7kb apart in Pel RR, whereas these non-amplified genes are only 1.7kb apart in the nonselected parental COL stock, as in Pel SS (susceptible Sri Lankan strain), demonstrating that this region of the genome is susceptible to expansion and contraction.

Tissue localization of overproduced esterases in the mosquito Culex pipiens (Diptera: Culicidae)

We have investigated the tissue distribution of overproduced esterases A (A1 and A2) and B (B1 and B2) in strains of Culex pipiens L. by immunocytochemistry. S-LAB mosquitoes, lacking overproduced esterases, were used as reference. Tissues showing a strong specific reaction (fluorescence) were observed with anti-esterase A1 antiserum in S54 (with A1) and BOUAKE (with A2) strains, and with anti-esterase B1 antiserum in TEM-R and EDIT (with B1) and BOUAKE (with B2) strains. Overproduction of esterases A and B was tissue-specific. The most constant pattern for the two types of esterases was their overproduction in the alimentary canal and Malpighian tubes, although fluorescence varied in intensity depending on strains and developmental stages. There was no difference in the tissue distribution of esterases Al and A2. In contrast, esterases B pattern was highly variable among strains. Differences between TEM-R and EDIT were explained by the different overall overproduction and number of copies of the amplified gene (10-fold higher in TEM-R). The most striking difference in esterase B1 and B2 tissue localization concerned the nervous system where neurons were intenisely fluorescent in TEM-R and EDIT (B1), but not in BOUAKE (B2). All esterase B positive tissues in TEM-R contained large quantities of esterase B1 mRNA (in situ hybridization), indicating that at least part of the protein revealed by immunochemistry was produced in the tissues where it was observed. Our results are discussed in terms of the protection that the different esterases can confer during exposition to organophosphorous insecticides.

Molecular and kinetic evidence for allelic variants of esterase Estbeta1 in the mosquito Culex quinquefasciatus

Elevated esterase Estbeta1 was purified from larvae of newly isolated strains of the mosquito Culex quinquefasciatus from Colombia (COL) and Trinidad (TRI) with resistance to organophosphate (OP) insecticides. Insecticide interactions were compared with those of elevated Estbeta1(2) from the OP-resistant Habana strain and the non-elevated Estbeta1(3) from the susceptible PelSS strain. On the basis of insecticide binding efficiency, all elevated Estbeta1 esterases were readily distinguishable. Differences between the EcoRI restriction fragment patterns of the amplified estbeta1 gene in COL and TRI strains compared with each other, and between amplified estbeta1(1), estbeta1(2) and the non-amplified estbeta1(3), suggest differences in their nucleotide sequence. Considering their variable insecticide binding efficiencies, these genetic differences would imply that, in contrast to estalpha2 and estbeta2, amplification of estbeta1 has occurred several times independently. Generally, the elevated Estbeta1s were more reactive with insecticides than the non-elevated Estbeta1(3). This supports the hypothesis that the elevated esterase-based mechanism confers resistance through amplification of alleles coding for esterases which have a greater specificity for the insecticides they sequester than the esterases coded by their non-amplified counterparts.

An overview of the evolution of overproduced esterases in the mosquito Culex pipiens

Insecticide resistance genes have developed in a wide variety of insects in response to heavy chemical application. Few of these examples of adaptation in response to rapid environmental change have been studied both at the population level and at the gene level. One of these is the evolution of the overproduced esterases that are involved in resistance to organophosphate insecticides in the mosquito Culex pipiens. At the gene level, two genetic mechanisms are involved in esterase overproduction, namely gene amplification and gene regulation. At the population level, the co-occurrence of the same amplified allele in distinct geographic areas is best explained by the importance of passive transportation at the worldwide scale. The long-term monitoring of a population of mosquitoes in southern France has enabled a detailed study to be made of the evolution of resistance genes on a local scale, and has shown that a resistance gene with a lower cost has replaced a former resistance allele with a higher cost.

Resistance to insecticides in insect vectors of disease: est alpha 3, a novel amplified esterase associated with amplified est beta 1 from insecticide resistant strains of the mosquito Culex quinquesfasciatus

Vector control programmes in many countries face the dual problems of parasite drug resistance and insecticide resistance in the insect vectors of the disease. Here we report for the first time a new esterase-based insecticide resistance mechanism in the filariasis vector Culex quinquefasciatus. The field collected COL strain of C. quinquefasciatus from Columbia was heterogeneous for organophosphorus insecticide resistance. On native polyacrylamide gels it had an elevated beta-naphthyl acetate specific esterase with the same Rf as that for the Est beta 1s involved in insecticide resistance in other strains of this mosquito species. After five generations of temephos insecticide selection, both the esterase specific activity with p-nitrophenyl acetate and the temephos LC50 values were increased, suggesting that elevation of esterase activity was the underlying mechanism of resistance. Western blots with antisera raised to Est alpha 2(1) and Est beta 2(1) from C. quinquefasciatus indicated that the COL strain had an elevated Est alpha 3 enzyme which co-migrated on native gels with Est beta 1. Southern blots indicated that an est alpha 3 gene was amplified in the COL strain and a Cuban mosquito strain (MRes), although the restriction digest patterns of the est beta 1 genes in these two strains are different. In contrast, the Californian TEMR strain, with the amplified est beta 1(1) gene, had no associated elevated Est alpha. Restriction digest patterns for COL and TEMR DNA suggest that they contain an identical est beta 1(1) gene, but our data suggest that the est alpha 3 gene occurs on the same amplicon as an est beta 1 gene although the genes are probably > 10 kb apart. Hence, either the COL strain has two est beta 1 genes or the est beta 1(1) amplicon in TEMR has been disrupted at some stage during the long colonisation of this strain and the amplified est alpha has been lost.

Contrasting levels of variability between cytoplasmic genomes and incompatibility types in the mosquito Culex pipiens

Reproductive incompatibilities called cytoplasmic incompatibilities are known to affect a large number of arthropod species and are mediated by Wolbachia, a maternally transmitted microorganism. The crossing relationships between strains of potential hosts define their incompatibility types and it is generally assumed that differences between strains of Wolbachia induce different crossing types. Among all the described host species, the mosquito, Culex pipiens, displays the greatest variability of cytoplasmic incompatibility crossing types. We analysed mitochondrial and bacterial DNA variability in Culex pipiens in order to investigate some possible causes of incompatibility crossing type variability. We sequenced fragments of the ftsZ gene, and the A + T-rich control region of the mtDNA. We also sequenced the second subunit of the mitochondrial cytochrome oxidase (COII) gene, in Culex pipiens and a closely related species, C. torrentium, in order to verify the usefulness of the A + T-rich region for the present purposes. No variability was found in the Wolbachia ftsZ gene fragment, and very limited variation of the mitochondrial marker whatever the compatibility type or the origin of the host. A low variability was found in the A + T-rich region and comparison of divergence of the A + T-rich region and COII gene between C. pipiens and C. torrentium did not reveal any special constraints affecting this region. In contrast to observations in other host species, variability of incompatibility crossing types is not due to multiple infections by distantly related Wolbachia strains.

Testing the unique amplification event and the worldwide migration hypothesis of insecticide resistance genes with sequence data

In the mosquito Culex pipiens, over-production of esterases is a common insecticide resistance mechanism. Different alleles at the esterase A and B loci are known to occur in natural populations, as shown by enzyme electrophoresis and RFLP studies on the esterase B locus. Here we analyse the variability of the esterase A locus at the nucleotide level in mosquitoes possessing or lacking over-produced esterase A. A surprisingly high level of nucleotide polymorphism is found in coding and noncoding regions, and the extent of polymorphism detected is higher than that previously described for the most polymorphic loci in Drosophila. We also show that eight strains from different localities (Africa, America, Asia) that possess the over-produced esterase A2 share the same nucleotide sequence at the esterase A locus, strongly supporting the evolutionary scenario of a unique event of amplification of this gene followed by dispersal through migration rather than the occurrence of multiple independent mutational events of the gene.