Organophosphate and carbamate insecticides irreversibly inhibit acetylcholinesterase causing death of insects. Resistance-modified acetylcholinesterases(AChEs) have been described in many insect species and sequencing of their genes allowed several point mutations to be described. However, their relative frequency and their cartography had not yet been addressed.
To analyze the most frequent mutations providing insecticide resistance in Drosophila melanogaster acetylcholinesterase, the Ace gene was cloned and sequenced in several strains harvested from different parts of the world. Sequence comparison revealed four widespread mutations, I161V, G265A, F330Y and G368A. We confirm here that mutations are found either isolated or in combination in the same protein and we show that most natural populations are heterogeneous, composed of a mixture of different alleles. In vitro expression of mutated proteins showed that combining mutations in the same protein has two consequences: it increases resistance level and provides a wide spectrum of resistance.
The presence of several alleles in natural populations, offering various resistance to carbamate and organophosphate compounds will complicate the establishment of resistance management programs.
Keywords:acetylcholinesterase; resistance; insecticide; mutation; population structure
Acetylcholinesterase (AChE; EC 188.8.131.52) is a key enzyme of the cholinergic system because it regulates the level of acetylcholine and terminates nerve impulses by catalyzing the hydrolysis of acetylcholine. Its inhibition causes death, so irreversible inhibitors have been developed as insecticides: organophosphates and carbamates. They have similar properties to acetylcholine but are hemisubstrates because they phosphorylate or carbamoylate the active-site serine leading to irreversible inhibition of the enzyme. This inhibition leads to an accumulation of acetylcholine in the synapses which in turn leaves the acetylcholine receptors permanently open, resulting in the death of the insect .
In 1961, Smissaert described the first case of AChE with a reduced sensitivity to pesticides . Since then, resistance-modified AChEs have been described in many insect species [3-7]. Sequencing of the gene encoding AChE in resistant strains showed that the modifications arose from point mutations (Table 1). Combinations of several point mutations in the same protein were found in several alleles where they induced higher levels of resistance [8,9]. Most mutations were identical in several species, suggesting that a low number of mutations can actually provide resistance. These findings presented a striking contrast to experiments of in vitro expression of mutagenesized AChE, which revealed that insecticide resistance should have genetic diversity . Besides qualitative modification of the enzyme, overproduction of AChE results in insecticide resistance as shown first experimentally by transforming Drosophila to increase the dose of enzyme and second by finding a positive correlation between the amount of AChE and resistance in natural populations [11,12].
Table 1. Point mutations in acetylcholinesterase involved in insecticide resistance.
The evolution of insecticide resistance in insects tends to be rapid because selection is strong, populations are large, and generation times are short. With the threat of insecticide resistance looming larger, it is absolutely necessary to investigate the molecular basis of AChE mutation distribution in natural populations from different parts of the world. In this paper, we report the most frequent point mutations of AChE genes associated with insecticide resistance and the different patterns of mutation combinations in natural populations of D. melanogaster.
Point mutations detected in D. melanogaster populations
To identify mutations involved in insecticide resistance, the Ace gene encoding AChE in 30 strains of D. melanogaster harvested throughout the world was sequenced. Sequences were compared to a reference, the sequence from a strain harvested before the utilisation of insecticide (Canton-S strain). Three clones were sequenced per strain, when variability was detected, three new clones were sequenced. We found 18 mutations resulting in a change of the amino-acid sequence of the protein. These mutations were verified by PASA to eliminate the possibility of a PCR artefact. Among these mutations some were found only once and their positions were far from the active site, so should not be responsible for insecticide resistance (V19M, F187L, E193V, N228D, E281G, N300S, M305T, S325P, E453V, G489D, K507D).
Two mutations (E72G and E81K) were detected in the strains STIC from Italy and NY from the U.S.A. respectively. They are located near the entrance of the active site gorge in a position favourable to affect insecticide sensitivity. Similarly, another mutation located at the rim of the gorge (F77S) has already been reported from the strain Saltillo, which affects the sensitivity to insecticide , but this mutation was not found again in the present screening. As these three mutations (E72G, E81K and F77S) were present only in one population, their effects were not studied in detail in this paper.
Four mutations, I161V, G265A, F330Y and G368A, which are located near the active site, were found in several natural populations (Table 2). These mutations have already been described in Drosophila and/or in other insects (Table 1). The reference sequence I161, G265, F330 and G368 (allele IGFG) most probably represents the wild type sequence since then same amino-acids were present in the other available AChE sequences from Brachycera (Musca domestica, Lucilia cuprina and Bractrocera oleae) and each mutation arose from a single mutation. At position 161, the mutation ATC to GTC changes the isoleucine to valine. In one strain (NY), we found another mutation (GTC to GCC) changing the valine to alanine (allele AGFG, Table 2). Compared to valine, this mutation provided higher resistance to some insecticides: for example, the alanine mutant exhibited 10- and 7-fold less sensitivity than the valine mutant to diazinon-oxon and monocrotophos respectively (data not shown). However, as this mutation has only been found once and is not widespread, we did not study it thoroughly. At position 368, two different codons (GCC or GCT) encode alanine suggesting either that the mutation of glycine to alanine originated from two independent events or that the C/T transition occurred in the GCC allele. The fact that we did not find any GGT alleles in natural populations favours the second hypothesis.
Table 2. Alleles found in natural populations of D. melanogaster.
Other mutations may be present in natural populations. Only 30 populations were analyzed, and some parts of the world were missing from the study or weakly represented. Furthermore the rearing of some populations in laboratory conditions for a long time without insecticide treatments may have caused the loss of certain mutations. However, I161V, G265A, F330Y and G368A seem to be the most frequent mutations providing insecticide resistance. This is consistent with the observation that these mutations were also found in other species (Table 1).
The four mutations can be amassed in one allele as shown in Table 2 (the allele IGFG represents the wild type allele and VAYA the allele with the four mutations). Various combinations were found in natural populations, including the allele with all four mutations. Although screening was not performed to estimate the proportion of each allele, we did, however, note that allele IAFG with a single mutation at position 265 and allele VAYG with a triple mutation (at positions 161, 265 and 368) were the most frequent (Table 3).
Table 3. Number of each allele found in natural populations of D. melanogaster.
Geographic distribution of point mutations I161V, G265A, F330Y and G368A
To tackle the geographic repartition of the four mutations, we checked their presence by PASA in several strains from several parts of the word (Table 4). At least ten flies were analyzed per strain. Some mutations were not detected by sequencing (Table 2) showing that the not all alleles in each population were detected. It appears that the four mutations are distributed all around the world. The worldwide distribution of these mutations might result from the ability of Drosophila to disperse either as larvae in fruits or as adults.
Table 4. Geographic distribution of point mutations at positions 161, 265, 330 and 368 determined by PASA.
Effect of mutations I161V, G265A, F330Y and G368A on insecticide resistance
Site-directed mutagenesis was used to study the effect of the four mutations on insecticide sensitivity. Mutated proteins were expressed in the baculovirus system and proteins were purified to homogeneity. The four single mutations, I161V, G265A, F330Y and G368A, affect insecticide sensitivity (Table 5) but show a strong disparity with regard to the 17 insecticides tested. For example resistance to coumaphos-oxon was obtained by mutation G265A (allele IAFG) and resistance to diazinon-oxon by mutations F330Y (allele IGYG).
Table 5. Effects of different patterns of mutation combinations on insecticide resistance represented as ki ratio (ki reference/ki mutant enzyme).
The combinations of mutations have two different types of effect. First, they provide higher resistance, as previously reported . For each insecticide, the allele which provides higher resistance corresponds to a combination of three or four point mutations (Table 5). However, this effect is not systematic, for example, the double mutant (allele VGYG) does not change the sensitivity to insecticides except for paraoxon, when compared with sensitivity of single mutants (alleles VGFG and IGYG). Second, combinations decrease the specificity of resistance: a single mutation provides resistance to 76% of insecticides (and sensitivity in 24%), double mutations provide resistance to 78% of insecticides, triple mutations provide resistance to 97% of insecticides and the co-occurrence of all four mutations provides resistance to all the insecticides tested.
Heterogeneity of populations and insecticide resistance
As shown by cDNA sequencing, several alleles coexist in most natural populations. For example, there were five alleles found in strain WC2 from USA and RIC from Costa Rica and four alleles in TORREON strain from Mexico and NDL strain from India. Because each allele manifested a specific resistance to certain insecticides, this allelic diversity in one natural population might enable the population to survive against a large variety of insecticide treatments. Correspondingly, subjecting field populations to a multiplicity of treatments might induce the emergence of multiple coexisting resistance alleles.
Practical consequences of the combination of mutations involved in insecticide resistance
In several cases, single mutations sensitised the enzyme to the insecticide (ki ratio was lower than 1.0, see Table 5). For instance, the mutation G368A rendered sensitivity to coumaphos oxon 4.0 times greater than that of the reference enzyme. The I161V/F330Y mutation caused 7.7 fold more sensitivity to monocrotophos than the reference enzyme.
One method to decrease the resistance of a field population is to apply an anti-resistance insecticide, which is an insecticide more active against the mutated protein than against the wild type protein. This strategy seems feasible when there is only one mutation in the population; for example, carbaryl can be used to decrease the frequency of the G368A mutant (allele IGFA). However, with the combination of multiple mutations, this strategy is less and less feasible, since the protein becomes resistant to all the insecticides. Furthermore, field populations are composed of a mixture of different alleles with different sensitivities to each insecticide, so treatment with one pesticide would eliminate one allele, but would select another one. For instance, treatment with carbaryl will decrease the frequency of the allele IGFA but will increase the frequency of the allele VAYG. So, the anti-resistance insecticide strategy does not seem to be efficient in field conditions.
Material and Methods
Strains of D. melanogaster were collected from several parts of the world (Table 2, 3, 4). They were maintained for several generations in laboratory conditions, without insecticide treatment. Canton-S strain, as a reference, has been cultured in the laboratory since its collection at the beginning of the 20th century in the USA.
RT-PCR, cloning, and sequencing
Total RNA was prepared using the RNeasy mini kit from Qiagen and mRNA was purified using oligodT-cellulose. First strand cDNA was obtained with random hexanucleotides as primers using the first-strand cDNA synthesis kit from Pharmacia Biotech. Ace cDNA were PCR amplified using primers hybridizing on the sequence encoding the signal peptide and on the sequence encoding the C-terminal hydrophobic peptide. As these two peptides were taken off during the maturation of the protein, the whole sequence of the mature protein was available. RT-PCR products were cloned into the P3T vector digested by Xcm1 and Sma1 .
Genomic DNA extraction and point mutations confirmed by PASA
Genomic DNA was isolated from each population of D. melanogaster by phenol/chloroform/isoamyl alcohol extraction. PCR amplification of specific alleles (PASA) was used to detect point mutations in natural populations of D. melanogaster. PASA is a modification of PCR that depends on a PCR oligonucleotide primer that precisely matches one of the alleles but mismatches with the other. When the mismatch occurs at the 3' end of the PCR primer, amplification is inefficient. Therefore, only amplification of the perfectly matched allele is obtained. It was performed using three oligonucleotides to exclude the possibility of PCR artefacts .
Site-directed mutagenesis and protein preparation
We generated point mutations using a standard PCR based strategy and verified the identities of individual clones through double stranded plasmid sequencing. DNA encoding soluble AChEs, reference type or mutated, of D. melanogaster were expressed with the baculovirus system . Secreted AChEs were purified to homogeneity and stabilized with 1 mg/ml BSA as previously reported . Amino-acid numbering follows that of the mature protein as in the structure published by Harel et al.  (http://www.rcsb.org/pdb webcite; PDB code: 1QO9) to facilitate the visualisation of the mutation on the structure. It corresponds to the numeration of the precursor used in a previous paper  minus the 38 amino-acids of the signal peptide.
Determination of enzyme sensitivity
The inhibition mechanism of AChE by organophosphate and carbamate compounds has been described by Aldridge .
with E = enzyme, PX = organophosphate or carbamate, X = leaving group. The insecticide phosphorylates or carbamoylates the active site serine of one enzyme molecule and the inhibition can be considered as irreversible during the first 30 min. Disappearance of free enzyme ([E]) follows second-order kinetics.
where [PXo] is the initial concentration of inhibitor, [E] is the free enzyme remaining at time t, [Eo] is the initial concentration of enzyme, t represents the time of incubation, and ki the biomolecular rate constant. [PXo] and t are known, [E]/ [Eo] is estimated by incubating the enzyme with the insecticide at 25°C in 25 mM phosphate buffer pH7. The variation of the remaining free enzyme [E]/ [Eo] with time was estimated by sampling aliquots at various times and recording the remaining activity ([A]/ [Ao]) with 1 mM acetylthiocholine since [A]/ [Ao] = [E]/ [Eo]. Kinetics studies were performed with at least three concentrations of insecticide. The values of ki were estimated by multiple non-linear regression with [PXo] and t as variables. Data were collected until standard deviation came below 10% of the ki value. The resistance level of the mutated enzyme was expressed as the ratio of ki (ki reference/ ki mutant enzyme)
PM sequenced alleles and verified the occurrence of mutations by PASA. MAS purified the recombinant proteins and estimated their sensitivity. AL performed in vitro mutagenesis and production of mutants. ZHT participated in coordination and DF conceived the study. All authors read and approved the final manuscript.
We thank Tom Wilson (Colorado State Univ.), Eric Bonnivard (Univ. Paris 6), Philip Daborn (Melbourne Univ.) for providing Drosophila strains. This research was supported by grants from the European Community (ACHEB, QLK3-CT-2000-00650 and SAFEGUARD, QLK3-CT-2000-000481) and from the Chinese-French Advanced Research Program (PRA BT01-01).
Science 1964, 143:129-131. PubMed Abstract
Pest Biochem Physiol 1994, 49:114-120. Publisher Full Text
Pestic Sci 1997, 49:90-96. Publisher Full Text
Devonshire AL, Byrne FJ, Moores GD, Williamson MS: Biochemical and molecular characterisation of insecticide-insensitive acetylcholinesterases. In In Structure and Function of Cholinesterases and Related Proteins. Edited by Doctor BP, Taylor P, Quinn DM, Rotundo RL, Gentry MK. Plenum Press, New York; 1998:491-496.
Pest Biochem Physiol 2000, 67:95-102. Publisher Full Text
Pest Biochem Physiol 2001, 70:100-107. Publisher Full Text
PCR Methods Appl 1992, 2:81-82. PubMed Abstract
Biotechniques 1992, 12:82-87. PubMed Abstract
Harel M, Kryger G, Rosenberry TL, Mallender WD, Lewis T, Fletcher RJ, Guss JM, Silman I, Sussman JL: Three-dimensional structures of Drosophila melanogaster acetylcholinesterase and of its complexes with two potent inhibitors.
Protein Sci 2000, 9:1063-1072. PubMed Abstract
Walsh SB, Dolden TA, Moores GD, Kristensen M, Lewis T, Devonshire AL, Williamson MS: Identification and characterization of mutations in housefly (Musca domestica) acetylcholinesterase involved in insecticide resistance.
Nabeshima T, Mori A, Kozaki T, Iwata Y, Hidoh O, Harada S, Kasai S, Severson DW, Konoa Y, Tomita T: An amino acid substitution attributable to insecticide-insensitivity of acetylcholinesterase in a Japanese encephalitis vector mosquito, Culex tritaeniorhynchus.