Evolutionary genetics of Drosophila esterases

Comparative evolutionary analyses of eight whitefly Bemisia tabaci sensu lato genomes: cryptic species, agricultural pests and plant-virus vectors

BACKGROUND

The group of > 40 cryptic whitefly species called Bemisia tabaci sensu lato are amongst the world's worst agricultural pests and plant-virus vectors. Outbreaks of B. tabaci s.l. and the associated plant-virus diseases continue to contribute to global food insecurity and social instability, particularly in sub-Saharan Africa and Asia. Published B. tabaci s.l. genomes have limited use for studying African cassava B. tabaci SSA1 species, due to the high genetic divergences between them. Genomic annotations presented here were performed using the 'Ensembl gene annotation system', to ensure that comparative analyses and conclusions reflect biological differences, as opposed to arising from different methodologies underpinning transcript model identification.

RESULTS

We present here six new B. tabaci s.l. genomes from Africa and Asia, and two re-annotated previously published genomes, to provide evolutionary insights into these globally distributed pests. Genome sizes ranged between 616-658 Mb and exhibited some of the highest coverage of transposable elements reported within Arthropoda. Many fewer total protein coding genes (PCG) were recovered compared to the previously published B. tabaci s.l. genomes and structural annotations generated via the uniform methodology strongly supported a repertoire of between 12.8-13.2 × 10³ PCG. An integrative systematics approach incorporating phylogenomic analysis of nuclear and mitochondrial markers supported a monophyletic Aleyrodidae and the basal positioning of B. tabaci Uganda-1 to the sub-Saharan group of species. Reciprocal cross-mating data and the co-cladogenesis pattern of the primary obligate endosymbiont 'Candidatus Portiera aleyrodidarum' from 11 Bemisia genomes further supported the phylogenetic reconstruction to show that African cassava B. tabaci populations consist of just three biological species. We include comparative analyses of gene families related to detoxification, sugar metabolism, vector competency and evaluate the presence and function of horizontally transferred genes, essential for understanding the evolution and unique biology of constituent B. tabaci. s.l species.

CONCLUSIONS

These genomic resources have provided new and critical insights into the genetics underlying B. tabaci s.l. biology. They also provide a rich foundation for post-genomic research, including the selection of candidate gene-targets for innovative whitefly and virus-control strategies.

Derived esterase activity in Drosophila sechellia contributes to evolved octanoic acid resistance

The dietary specialist fruit fly Drosophila sechellia has evolved resistance to the secondary defence compounds produced by the fruit of its host plant, Morinda citrifolia. The primary chemicals that contribute to lethality of M. citrifolia are the medium-chain fatty acids octanoic acid (OA) and hexanoic acid. At least five genomic regions contribute to this adaptation in D. sechellia and whereas the fine-mapped major effect locus for OA resistance on chromosome 3R has been thoroughly analysed, the remaining four genomic regions that contribute to toxin resistance remain uncharacterized. To begin to identify the genetic basis of toxin resistance in this species, we removed the function of well-known detoxification gene families to determine whether they contribute to toxin resistance. Previous work found that evolution of cytochrome P450 enzymatic activity or expression is not responsible for the OA resistance in D. sechellia. Here, we tested the role of the two other major detoxification gene families in resistance to Morinda fruit toxins - glutathione-S-transferases and esterases - through the use of the pesticide synergists diethyl maleate and tribufos that inhibit the function of these gene families. This work suggests that one or more esterase(s) contribute to evolved OA resistance in D. sechellia.

The molecular basis for the neofunctionalization of the juvenile hormone esterase duplication in Drosophila

The Drosophila melanogaster enzymes juvenile hormone esterase (DmJHE) and its duplicate, DmJHEdup, present ideal examples for studying the structural changes involved in the neofunctionalization of enzyme duplicates. DmJHE is a hormone esterase with precise regulation and highly specific activity for its substrate, juvenile hormone. DmJHEdup is an odorant degrading esterase (ODE) responsible for processing various kairomones in antennae. Our phylogenetic analysis shows that the JHE lineage predates the hemi/holometabolan split and that several duplications of JHEs have been templates for the evolution of secreted β-esterases such as ODEs through the course of insect evolution. Our biochemical comparisons further show that DmJHE has sufficient substrate promiscuity and activity against odorant esters for a duplicate to evolve a general ODE function against a range of mid-long chain food esters, as is shown in DmJHEdup. This substrate range complements that of the only other general ODE known in this species, Esterase 6. Homology models of DmJHE and DmJHEdup enabled comparisons between each enzyme and the known structures of a lepidopteran JHE and Esterase 6. Both JHEs showed very similar active sites despite low sequence identity (30%). Both ODEs differed drastically from the JHEs and each other, explaining their complementary substrate ranges. A small number of amino acid changes are identified that may have been involved in the early stages of the neofunctionalization of DmJHEdup. Our results provide key insights into the process of neofunctionalization and the structural changes that can be involved.

Host use and host shifts in Drosophila

Over a thousand Drosophila species have radiated onto a wide range of feeding and breeding sites. These radiations involve adaptations for locating, accepting, and growing in hosts with highly differing characteristics. In a number of species, owing to the availability of sequenced genomes, particular steps in host specialization and genes that control them, are being identified. Many cases of specialization involve the ability to detoxify some component of the host. Examples include Drosophila sechellia and the octanoic acid in Morinda citrifolia, alpha-amanitin in mycophagous drosophilids, and the alkaloids in cactophilic species. Owing to the known ecologies of many species for which genomes exist, the Drosophila model system provides an unprecedented opportunity to simultaneously examine the genes underlying HOST LOCATION, HOST ACCEPTANCE and HOST USE, the types of selection acting upon them and any coevolutionary interactions among the genes underlying these steps.

Esterase mediated resistance in deltamethrin resistant reference tick colony of Rhipicephalus (Boophilus) microplus

Monitoring of acaricide resistance is considered as one of the important facets of integrated tick management. In an attempt of development of resistance monitoring indicators, in the present study two reference tick lines of Rhipicephalus (Boophilus) microplus maintained in the Entomology laboratory, Indian Veterinary Research Institute (IVRI), Izatnagar, India, were studied to determine the possible contributing factors involved in development of resistance to deltamethrin. Electrophoretic profiling of esterase enzymes detected high activities of EST-1 in reference resistant tick colony designated as IVRI-IV whereas it was not detectable in reference susceptible IVRI-I line of R. (B.) microplus. Esterases were further characterized as carboxylesterase or acetylcholinesterase based on inhibitor study using PMSF, eserine sulphate, malathion, TPP and copper sulphate. It was concluded that an acetylcholinesterase, EST-1, possibly plays an important role for development of deltamethrin resistance in IVRI-IV colony of R. (B.) microplus.

Esterases of Varroa destructor (Acari: Varroidae), parasitic mite of the honeybee

Varroa destructor is an ectoparasite that causes serious damage to the population of the honeybee. Increasing resistance of the parasite to acaricides is related, among others, to metabolic adaptations of its esterases to facilitate decomposition of the chemicals used. Esterases are a large heterogeneous group of enzymes that metabolize a number of endogenous and exogenous substrates with ester binding. The aim of the present study was to determine the activity of esterases in the body extracts (BE) and excretion/secretion products (E/SP) of the mite. The enzymes contained in the E/SP should originate mainly from the salivary glands and the alimentary system and they may play a particularly important role in the first line of defence of the mite against acaricides. Activity of cholinesterases (ChEs) [acetylcholinesterase (AChE) and butyrylcholinesterase], carboxylesterases (CEs) and phosphatases [alkaline phosphatase (AP) and acid phosphatase (AcP)] was investigated. The activity of all the enzymes except AChE was higher in the E/SP than in the BE. ChEs from the BE and from the E/SP reacted differently on eserine, a ChE inhibitor. Eserine inhibited both enzymes from the BE, increased decomposition of acetylcholine, but did not influence hydrolysis of butyrylcholine by the E/SP. Activity of the CEs from the BE in relation to the esters of carboxylic acids can be presented in the following series: C10 > C12 > C14 > C8 > C2 > C4 = C16, while activity of the CEs from the E/SP was: C4 > C8 > C2 > C14 > C10 > C12 > C16. The inhibitor of CEs, triphenyl phosphate, reduced the activity of esterases C2–C8 and C14–C16; however, it acted in the opposite way to CEs C10 and C12. The activity of both phosphatases was higher in the E/SP than in the BE (AcP about twofold and AP about 2.6-fold); the activities of AP and AcP in the same material were similar. Given the role of esterases in resistance to pesticides, further studies are necessary to obtain complete biochemical characteristics of the enzymes currently present in V. destructor.

Computational methods of identification of pseudogenes based on functionality: entropy and GC content

Spectral entropy and GC content analyses reveal comprehensive structural features of DNA sequences. To illustrate the significance of these features, we analyze the β-esterase gene cluster, including the Est-6 gene and the ψEst-6 putative pseudogene, in seven species of the Drosophila melanogaster subgroup. The spectral entropies show distinctly lower structural ordering for ψEst-6 than for Est-6 in all species studied. However, entropy accumulation is not a completely random process for either gene and it shows to be nucleotide dependent. Furthermore, GC content in synonymous positions is uniformly higher in Est-6 than in ψEst-6, in agreement with the reduced GC content generally observed in pseudogenes and nonfunctional sequences. The observed differences in entropy and GC content reflect an evolutionary shift associated with the process of pseudogenization and subsequent functional divergence of ψEst-6 and Est-6 after the duplication event. The data obtained show the relevance and significance of entropy and GC content analyses for pseudogene identification and for the comparative study of gene-pseudogene evolution.

Heterologous expression and biochemical characterisation of fourteen esterases from Helicoverpa armigera

Esterases have recurrently been implicated in insecticide resistance in Helicoverpa armigera but little is known about the underlying molecular mechanisms. We used a baculovirus system to express 14 of 30 full-length esterase genes so far identified from midgut cDNA libraries of this species. All 14 produced esterase isozymes after native PAGE and the isozymes for seven of them migrated to two regions of the gel previously associated with both organophosphate and pyrethroid resistance in various strains. Thirteen of the enzymes obtained in sufficient yield for further analysis all showed tight binding to organophosphates and low but measurable organophosphate hydrolase activity. However there was no clear difference in activity between the isozymes from regions associated with resistance and those from elsewhere in the zymogram, or between eight of the isozymes from a phylogenetic clade previously associated with resistance in proteomic and quantitative rtPCR experiments and five others not so associated. By contrast, the enzymes differed markedly in their activities against nine pyrethroid isomers and the enzymes with highest activity for the most insecticidal isomers were from regions of the gel and, in some cases, the phylogeny that had previously been associated with pyrethroid resistance. Phospholipase treatment confirmed predictions from sequence analysis that three of the isozymes were GPI anchored. This unusual feature among carboxylesterases has previously been suggested to underpin an association that some authors have noted between esterases and resistance to the Cry1Ac toxin from Bacillus thuringiensis. However these three isozymes did not migrate to the zymogram region previously associated with Cry1Ac resistance.

Microarray analysis of gene regulations and potential association with acephate-resistance and fitness cost in Lygus lineolaris

The tarnished plant bug has become increasingly resistant to organophosphates in recent years. To better understand acephate resistance mechanisms, biological, biochemical, and molecular experiments were systematically conducted with susceptible (LLS) and acephate-selected (LLR) strains. Selection of a field population with acephate significantly increased resistance ratio to 5.9-fold, coupled with a significant increase of esterase activities by 2-fold. Microarray analysis of 6,688 genes revealed 329 up- and 333 down-regulated (≥2-fold) genes in LLR. Six esterase, three P450, and one glutathione S-transferase genes were significantly up-regulated, and no such genes were down-regulated in LLR. All vitellogenin and eggshell protein genes were significantly down-regulated in LLR. Thirteen protease genes were significantly down-regulated and only 3 were up-regulated in LLR. More than twice the number of catalysis genes and more than 3.6-fold of metabolic genes were up-regulated, respectively, as compared to those down-regulated with the same molecular and biological functions. The large portion of metabolic or catalysis genes with significant up-regulations indicated a substantial increase of metabolic detoxification in LLR. Significant increase of acephate resistance, increases of esterase activities and gene expressions, and variable esterase sequences between LLS and LLR consistently demonstrated a major esterase-mediated resistance in LLR, which was functionally provable by abolishing the resistance with esterase inhibitors. In addition, significant elevation of P450 gene expression and reduced susceptibility to imidacloprid in LLR indicated a concurrent resistance risk that may impact other classes of insecticides. This study demonstrated the first association of down-regulation of reproductive- and digestive-related genes with resistance to conventional insecticides, suggesting potential fitness costs associated with resistance development. This study shed new light on the understanding of the molecular basis of insecticide resistance, and the information is highly valuable for development of chemical control guidelines and tactics to minimize resistance and cross-resistance risks.

Genetic variability in esterases and the insecticide resistance in brazilian strains of Oryzaephilus mercator and Oryzaephilus surinamensis (Coleoptera: Silvanidae)

Oryzaephilus mercator and O. surinamensis are stored grains and processed food pests, the latter being responsible for major economical losses. Polyacrylamide gel electrophoresis was used to analyse esterase patterns during insect development. Seven esterases, three cholinesterases, two carboxylesterases and two acetylesterases, were identified in O. mercator, one of which was proper to adults. Five esterases, of which four were cholinesterases, occurred in O. surinamensis. Strains of O. mercator and O. surinamensis were also exposed to malathion and chlorpyrifos-methyl. According to the LC50 estimates, OmLC-M and OmLA strains of O. mercator and OsLB strain of O. surinamensis were the most resistant to both insecticides. However, higher sensitivity to malathion and chlorpyrifos-methyl has also been verified in some of its esterases. Cholinesterases OmEST-1 and OsEST-5 seem to be involved in this resistance. These results suggest that organophosphate tolerance may be related to genetic variability in esterase isoenzymes.

Esterase patterns in species of the triatome bug Rhodnius

The aim of the present study was to analyse esterase patterns in three triatomine species of Rhodnius genus. Four loci, Est 1, Est 2, Est 3 and Est 4, were found. The corresponding enzymes were characterized as carboxylesterases (E.C. 3.1.1.1) or cholinesterases (E.C. 3.1.1.8) based on inhibitory experiments, using eserine sulphate, malathion, mercury chloride, p-chloromercuribenzoate (pCMB) and iodoacetamide. Low genetic variability was observed: Est 1, Est 2 and Est 3 were monomorphic in Rhodnius domesticus, Rhodnius robustus and Rhodnius neivai, whereas locus Est 4 was polymorphic in the first two species. The UPGMA analysis based on esterase genotypic frequencies indicated greater similarity between R. domesticus and R. robustus when compared with R. neivai. The present study expands our knowledge about genetic variability among triatomines and accords with the hypothesis that R. domesticus is a species derived from R. robustus.

Gene identification and proteomic analysis of the esterases of the cotton bollworm, Helicoverpa armigera

Some of the resistance of Helicoverpa armigera to conventional insecticides such as organophosphates and synthetic pyrethroids appears to be due to metabolic detoxification by carboxylesterases. To investigate the H. armigera carboxyl/cholinesterases, we created a data set of 39 putative paralogous H. armigera carboxyl/cholinesterase sequences from cDNA libraries and other sources. Phylogenetic analysis revealed a close relationship between these sequences and 70 carboxyl/cholinesterases from the recently sequenced genome of the silkworm, Bombyx mori, including several conserved clades of non-catalytic proteins. A juvenile hormone esterase candidate from H. armigera was identified, and B. mori orthologues were proposed for 31% of the sequences examined, however low similarity was found between lepidopteran sequences and esterases previously associated with insecticide resistance from other insect orders. A proteomic analysis of larval esterases then enabled us to match seven of the H. armigera carboxyl/cholinesterase sequences to specific esterase isozymes. All identified sequences were predicted to encode catalytically active carboxylesterases, including six proteins with N-terminal signal peptides and N-glycans, with two also containing C-terminal signals for glycosylphosphatidylinositol anchor attachment. Five of these sequences were matched to zones of activity on native PAGE at relative mobility values previously associated with insecticide resistance in this species.

Birth and death of genes and functions in the beta-esterase cluster of Drosophila

Here we analyze the molecular evolution of the β-esterase gene cluster in the Drosophila genus using the recently released genome sequences of 12 Drosophila species. Molecular evolution in this small cluster is noteworthy because it contains contrasting examples of the types and stages of loss of gene function. Specifically, missing orthologs, pseudogenes, and null alleles are all inferred. Phylogenetic analyses also suggest a minimum of 9 gene gain–loss events; however, the exact number and age of these events is confounded by interparalog recombination. A previous enigma, in which allozyme loci were mapped to β-esterase genes that lacked catalytically essential amino acids, was resolved through the identification of neighbouring genes that contain the canonical catalytic residues and thus presumably encode the mapped allozymes. The originally identified genes are evolving with selective constraint, suggesting that they have a “noncatalytic” function. Curiously, 3 of the 4 paralogous β-esterase genes in the D. ananassae genome sequence have single inactivating (frame-shift or nonsense) mutations. To determine whether these putatively inactivating mutations were fixed, we sequenced other D. ananassae alleles of these four loci. We did not find any of the 3 inactivating mutations of the sequenced strain in 12 other strains; however, other inactivating mutations were observed in the same 3 genes. This is reminiscent of the high frequency of null alleles observed in one of the β-esterase genes (Est7/EstP) of D. melanogaster.

Variability of esterase patterns in adult flies of the saltans species group of Drosophila (subgenus Sophophora)

Esterases are known for their involvement in several physiological processes and high degree of polymorphism, in many organisms. Such polymorphism has been used to characterize species and species groups and to study genetic changes occurred in their evolutionary history. In the present study, the esterase patterns of 19 strains from 10 species representative of the five subgroups of the saltans species group were analyzed using polyacrylamide gel electrophoresis and alpha- and beta- naphthyl acetates as substrates. Fifty-one esterase bands were detected and classified as 31 alpha-esterases, 18 beta-esterases and two alpha/beta-esterases. On the basis of the inhibition patterns using Malathion and eserine sulfate, 34 bands were classified as carboxylesterases, 14 as acethylesterases and three as cholinesterases. Ten gene loci were tentatively established on the basis of data on band position in the gel, substrate preference and inhibition pattern. Twenty bands were species-specific, the remaining being shared by species from the same or different subgroups. Bands detected exclusively in males and bands with a different frequency or degree of expression between sexes were also detected. In the gels prepared for analysis of gene expression in the body parts (head, thorax and abdomen), the degree of expression of the beta-esterases was higher in the thorax, while the alpha-esterases were expressed predominantly in the abdomen and thorax. A global view of the data available at present on the esterases of the species from the saltans group and their degree of polymorphism are presented, as well as the possibility of using some beta-esterases, because of their characteristics in the gels, as markers for species identification.

Molecular population genetics of the alpha-esterase5 gene locus in original and colonized populations of Drosophila buzzatii and its sibling Drosophila koepferae

Several studies have suggested that esterase-2 (EST-2) may be the target of natural selection in the cactophilic fly Drosophila buzzatii. In this work, we analyzed nucleotide variation in a fragment of alpha-esterase5 (alphaE5), the gene encoding EST-2, in original (Argentinian) and colonized (Australian) populations of D. buzzatii and in its sibling D. koepferae. Estimates of nucleotide heterozygosity in D. buzzatii were similar in Australia and Argentina, although we detected a loss of singletons in colonized populations, suggesting a moderate founder effect. Interspecific comparisons revealed that D. buzzatii was more polymorphic for nonsynonymous variation, whereas D. koepferae was more variable for synonymous and noncoding sites. The two major chromosomal arrangements (2st and 2j) in D. buzzatii displayed similar levels of nucleotide variation, whereas 2jz3 was monomorphic. The sequenced region allowed the discrimination of a greater number of EST-2 protein variants in the Australian sample than in the Argentinean sample. In D. koepferae, nucleotide variation in alphaE5 does not depart from neutral expectations, although tests of population structure were significant for silent variation. In contrast, D. buzzatii has probably undergone a recent population expansion in its South American range. In addition, the McDonald and Kreitman test revealed an excess of nonsynonymous polymorphism in both original and colonized populations of this species.

Further study on the esterase patterns of sibling species in the Drosophila saltans subgroup (saltans group): intraspecific and interspecific variations in the development

Twenty of the 32 esterase bands previously detected in the adults of D. prosaltans, D. saltans and D. austrosaltans were found in larvae and pupae studied in this work. The results showed that, in addition to expressing the highest number of esterase bands, the adult stage of the three species exhibited the highest degree of expression (amount of synthesis) for most of the bands. Differences between larval and pupal stages were detected in the degree of expression (amount of synthesis) of the bands and in the frequency of samples expressing them. The frequencies of expression of the bands corresponding to genes in loci 1-3 were greater in pupae than in larvae while the frequencies of expression of the bands corresponding to genes in loci 4-9 were predominantly expressed in larvae or were equal in both developmental stages. Like the adults, larvae, pupae and empty pupal cases (which were also studied in this work) showed specific esterases. Taken together, the observations showed that, in the species studied, every developmental stage is characterized by specific bands and by specific frequency and degree of expression of the bands shared with other stages.

Characterization of Esterases in a Brazilian Population of Zaprionus Indianus (Diptera: Drosophilidae)

The aim of this study was to characterize esterases in Zaprionus indianus, a drosophilid recently introduced into Brazil. A further aim was study the variation of activity of esterases in the presence of inhibitors and their expression according to sex, sexual activity and age of individual flies. Polymorphisms were detected in two esterase loci (Est-2 and Est-3) and monomorphisms in four others (Est-1, Est-4, Est-5 and Est-6). Biochemical tests using alpha- and beta-naphthyl acetate and the inhibitors malathion, eserine sulphate and PMSF allowed us to classify EST-2 and EST-5 as beta-esterases, both carboxyl-esterases, and EST-1, EST-3, EST-4 and EST-6 as alpha-esterases. EST-1 and EST-3 were classified as carboxyl-esterases and EST-4 and EST-6 as cholinesterases. EST-5 activity was more pronounced in males and EST-2 was restricted to them or to recently copulated females. EST-4, rarely detected, was not characterized. Based on their biochemical characteristics possible roles for these enzymes are suggested.

Genetic variation along time in a Brazilian population of Aedes aegypti (Diptera: Culicidae), detected by changes in the esterase patterns

Aedes aegypti from the Brazilian cities of São José do Rio Preto (SJ) and Goiânia (GO) were analyzed as to their esterase patterns and the results were compared with data obtained about 5 years before for SJ population. Esterase bands not detected in the previous study were now observed in mosquitoes from both SJ and GO populations, being the last considered a population resistant to insecticides. Other similarities between SJ and GO populations in this study, and some differences in comparison with the previous data on SJ were observed, involving, in addition to changes in band type, changes in frequency of mosquitoes expressing them and differential gene activation during development. As it is generally true for genetic features, changes in the esterase patterns are expected to be the result of factors such as selection by environmental conditions and genetic drift. In the present case, continuous use of insecticides aiming mosquito population size control in SJ by sanitary authorities could be involved in the observed changes. Changed esterases were classified as carboxylesterases and cholinesterases, which are enzymes already shown to take part in the development of resistance in several organisms. In addition, data obtained in the elapsed time by authorities responsible for the mosquito control has shown increasing insecticide resistance of SJ population mosquitoes parallel to increase in the total amount of esterases, reinforcing the mentioned possibility.

Entropy and GC Content in the beta-esterase gene cluster of the Drosophila melanogaster subgroup

We perform spectral entropy and GC content analyses in the beta-esterase gene cluster, including the Est-6 gene and the psiEst-6 putative pseudogene, in seven species of the Drosophila melanogaster species subgroup. psiEst-6 combines features of functional and nonfunctional genes. The spectral entropies show distinctly lower structural ordering for psiEst-6 than for Est-6 in all species studied. Our observations agree with previous results for D. melanogaster and provide additional support to our hypothesis that after the duplication event Est-6 retained the esterase-coding function and its role during copulation, while psiEst-6 lost that function but now operates in conjunction with Est-6 as an intergene. Entropy accumulation is not a completely random process for either gene. Structural entropy is nucleotide dependent. The relative normalized deviations for structural entropy are higher for G than for C nucleotides. The entropy values are similar for Est-6 and psiEst-6 in the case of A and T but are lower for Est-6 in the case of G and C. The GC content in synonymous positions is uniformly higher in Est-6 than in psiEst-6, which agrees with the reduced GC content generally observed in pseudogenes and nonfunctional sequences. The observed differences in entropy and GC content reflect an evolutionary shift associated with the process of pseudogenization and subsequent functional divergence of psiEst-6 and Est-6 after the duplication event.

Enhanced esterase gene expression and activity in a malathion-resistant strain of the tarnished plant bug, Lygus lineolaris

Extensive use of insecticides on cotton in the mid-South has prompted resistance development in the tarnished plant bug, Lygus lineolaris (Palisot de Beauvois). A field population of tarnished plant bugs in Mississippi with 11-fold higher resistance to malathion was used to examine how gene regulation conferred resistance to this organophosphate insecticide. In laboratory bioassays, synergism by the esterase inhibitors S,S,S,-tributylphosphorotrithioate (DEF) and triphenylphosphate (TPP) effectively abolished resistance and increased malathion toxicity by more than 80%. Esterase activities were compared in vitro between malathion susceptible and resistant (selected) strains. More than 6-, 3- and 10-fold higher activities were obtained with the resistant strain using alpha-naphthyl acetate, beta-naphthyl acetate, and p-nitrophenyl acetate, respectively. Up to 95% and 89% of the esterase activity in the susceptible and resistant strains, respectively, was inhibited by 1 mM DEF. Inhibition of esterase activity up to 75% and 85% in the susceptible and resistant strains, respectively, was obtained with 0.03 mM TPP. Esterase activities in field populations increased by up to 5.4-fold during the fall season. The increase was synchronized with movement of the insect into cotton where exposure to pesticides occurred. Esterase cDNA was cloned and sequenced from both malathion susceptible and resistant strains. The 1818-nucleotide cDNA contained a 1710-bp open reading frame coding a 570 amino acid protein which was similar to many insect esterases conferring organophosphate resistance. No amino acid substitution was observed between susceptible and resistant strains, indicating that esterase gene mutation was not involved in resistance development in the resistant strain in Mississippi. Further examination of esterase gene expression levels using quantitative RT-PCR revealed that the resistant strain had a 5.1-fold higher level of esterase mRNA than the susceptible strain. The results of this study indicated that up-regulation of the esterase gene appeared to be related to the development of resistance in the tarnished plant bug.

Hydrolysis of organophosphorus insecticides by in vitro modified carboxylesterase E3 from Lucilia cuprina

Resistance of the blowfly, Lucilia cuprina, to organophosphorus (OP) insecticides is due to mutations in LcalphaE7, the gene encoding carboxylesterase E3, that enhance the enzyme's ability to hydrolyse insecticides. Two mutations occur naturally, G137D in the oxyanion hole of the esterase, and W251L in the acyl binding pocket. Previous in vitro mutagenesis and expression of these modifications to the cloned gene have confirmed their functional significance. G137D enhances hydrolysis of diethyl and dimethyl phosphates by 55- and 33-fold, respectively. W251L increases dimethyl phosphate hydrolysis similarly, but only 10-fold for the diethyl homolog; unlike G137D however, it also retains ability to hydrolyse carboxylesters in the leaving group of malathion (malathion carboxylesterase, MCE), conferring strong resistance to this compound. In the present work, we substituted these and nearby amino acids by others expected to affect the efficiency of the enzyme. Changing G137 to glutamate or histidine was less effective than aspartate in improving OP hydrolase activity and like G137D, it diminished MCE activity, primarily through increases in Km. Various substitutions of W251 to other smaller residues had a broadly similar effect to W251L on OP hydrolase and MCE activities, but at least two were quantitatively better in kinetic parameters relating to malathion resistance. One, W251G, which occurs naturally in a malathion resistant hymenopterous parasitoid, improved MCE activity more than 20-fold. Mutations at other sites near the bottom of the catalytic cleft generally diminished OP hydrolase and MCE activities but one, F309L, also yielded some improvements in OP hydrolase activities. The results are discussed in relation to likely steric effects on enzyme-substrate interactions and future evolution of this gene.

Pseudogenes: are they "junk" or functional DNA?

Pseudogenes have been defined as nonfunctional sequences of genomic DNA originally derived from functional genes. It is therefore assumed that all pseudogene mutations are selectively neutral and have equal probability to become fixed in the population. Rather, pseudogenes that have been suitably investigated often exhibit functional roles, such as gene expression, gene regulation, generation of genetic (antibody, antigenic, and other) diversity. Pseudogenes are involved in gene conversion or recombination with functional genes. Pseudogenes exhibit evolutionary conservation of gene sequence, reduced nucleotide variability, excess synonymous over nonsynonymous nucleotide polymorphism, and other features that are expected in genes or DNA sequences that have functional roles. We first review the Drosophila literature and then extend the discussion to the various functional features identified in the pseudogenes of other organisms. A pseudogene that has arisen by duplication or retroposition may, at first, not be subject to natural selection if the source gene remains functional. Mutant alleles that incorporate new functions may, nevertheless, be favored by natural selection and will have enhanced probability of becoming fixed in the population. We agree with the proposal that pseudogenes be considered as potogenes, i.e., DNA sequences with a potentiality for becoming new genes.

Molecular population genetics of theβ-esterase gene cluster ofDrosophila melanogaster

We have investigated nucleotide polymorphism at the beta-esterase gene cluster including the Est-6 gene and psiEst-6 putative pseudogene in four samples of Drosophila melanogaster derived from natural populations of southern Africa (Zimbabwe), Europe (Spain), North America (USA: California), and South America (Venezuela). A complex haplotype structure is revealed in both Est-6 and psiEst-6. Total nucleotide diversity is twice in psiEst-6 as in Est-6; diversity is higher in the African sample than in the non-African ones. Strong linkage disequilibrium occurs within the beta-esterase gene cluster in non-African samples, but not in the African one. Intragenic gene conversion events are detected within Est-6 and, to a much greater extent, within psiEst-6; intergenic gene conversion events are rare. Tests of neutrality with recombination are significant for the beta-esterase gene cluster in the non-African samples but not significant in the African one. We suggest that the demographic history (bottleneck and admixture of genetically differentiated populations) is the major factor shaping the pattern of nucleotide polymorphism in the beta-esterase gene cluster. However there are some 'footprints' of directional and balancing selection shaping specific distribution of nucleotide polymorphism within the cluster. Intergenic epistatic selection between Est-6 and psiEst-6 may play an important role in the evolution of the beta-esterase gene cluster preserving the putative pseudogene from degenerative destruction and reflecting possible functional interaction between the functional gene and the putative pseudogene. Est-6 and psiEst-6 may represent an indivisible intergenic complex ('intergene') in which each single component (Est-6 or psiEst-6) cannot separately carry out the full functional role.

Nucleotide Variation of theEst-6Gene Region in Natural Populations ofDrosophila melanogaster

We have investigated nucleotide polymorphism in the Est-6 gene region in four samples of Drosophila melanogaster derived from natural populations of East Africa (Zimbabwe), Europe (Spain), North America (California), and South America (Venezuela). There are two divergent sequence types in the North and South American samples, which are not perfectly (North America) or not at all (South America) associated with the Est-6 allozyme variation. Less pronounced or no sequence dimorphism occurs in the European and African samples, respectively. The level of nucleotide diversity is highest in the African sample, lower (and similar to each other) in the samples from Europe and North America, and lowest in the sample from South America. The extent of linkage disequilibrium is low in Africa (1.23% significant associations), but much higher in non-African populations (22.59, 21.45, and 37.68% in Europe, North America, and South America, respectively). Tests of neutrality with recombination are significant in non-African samples but not significant in the African sample. We propose that demographic history (bottleneck and admixture of genetically different populations) is the major factor shaping the nucleotide patterns in the Est-6 gene region. However, positive selection modifies the pattern: balanced selection creates elevated levels of nucleotide variation around functionally important (target) polymorphic sites (RsaI–/RsaI+ in the promoter region and F/S in the coding region) in both African and non-African samples; and directional selection, acting during the geographic expansion phase of D. melanogaster, creates an excess of very similar sequences (RsaI– and S allelic lineages, in the promoter and coding regions, respectively) in the non-African samples.

Developmental expression and gene/enzyme identifications in the alpha esterase gene cluster of Drosophila melanogaster

Here we show how the 10 genes of the alpha esterase cluster of Drosophila melanogaster have diverged substantially in their expression profiles. Together with previously described sequence divergence this suggests substantial functional diversification. By peptide mass fingerprinting and in vitro gene expression we have also shown that two of the genes encode the isozymes EST9 (formerly ESTC) and EST23. EST9 is the major 'alpha staining' esterase in zymograms of gut tissues in feeding stages while orthologues of EST23 confer resistance to organophosphorus insecticides in other higher Diptera. The results for EST9 and EST23 concur with previous suggestions that the products of the alpha esterase cluster function in digestion and detoxification of xenobiotic esters. However, many of the other genes in the cluster show developmental or tissue-specific expression that seems inconsistent with such roles. Furthermore, there is generally poor correspondence between the mRNA expression patterns of the remaining eight genes and isozymes previously characterized by standard techniques of electrophoresis and staining, suggesting that the alpha cluster might only account for a small minority of the esterase isozyme profile.

DNA polymorphism in the beta-Esterase gene cluster of Drosophila melanogaster

We have analyzed nucleotide polymorphism within a 5.3-kb region encompassing the functional Est-6 gene and the psiEst-6 putative pseudogene in 28 strains of Drosophila melanogaster and one of D. simulans. Two divergent sequence types were detected, which are not perfectly associated with Est-6 allozyme variation. The level of variation (pi) is very close in the 5'-flanking region (0.0059) and Est-6 gene (0.0057), but significantly higher in the intergenic region (0.0141) and putative pseudogene (0.0122). The variation in the 3'-flanking region is intermediate (0.0083). These observations may reflect different levels of purifying selection in the different regions. Strong linkage disequilibrium occurs within the region studied, with the largest values revealed in the putative pseudogene and 3'-flanking region. Moreover, recombination is restricted within psiEst-6. Gene conversion is detected both within and (to a lesser extent) between Est-6 and psiEst-6. The data indicate that psiEst-6 exhibits some characteristics that are typical of nonfunctional genes, while other characteristics are typically attributed to functional genes; the same situation has been observed in other pseudogenes (including Drosophila). The results of structural entropy analysis demonstrate higher structural ordering in Est-6 than in psiEst-6, in accordance with expectations if psiEst-6 is indeed a pseudogene. Taking into account that the function of psiEst-6 is not known (but could exist) and following the terminology of J. Brosius and S. J. Gould, we suggest that the term "potogene" may be appropriate for psiEst-6, indicating that it is a potential gene that may have acquired some distinctive but unknown function.

Genetic polymorphism at two linked loci, Sod and Est-6, in Drosophila melanogaster

We have examined the patterns of polymorphism at two linked loci, Sod and Est-6, separated by nearly 1000 kb on the left arm of chromosome 3 of Drosophila melanogaster. The evidence suggests that natural selection has been involved in shaping the polymorphisms. At the Sod locus, a fairly strong (s>0.01) selective sweep, started >or=2600 years ago, increased the frequency of a rare haplotype, F(A), to about 50% frequency in populations of Europe, Asia, and the Americas. More recently, an F(A) allele mutated to an S allele, which has increased to frequencies 5-15% in populations of Europe, Asia and North America. All S alleles are identical (or very nearly) in sequence and differ by one nucleotide substitution (which accounts for the F-->S electrophoretic difference) from F(A) alleles. At the Est-6 locus, the evidence indicates both directional and balancing selection impacting separately the promoter and the coding regions of the gene, with linkage disequilibrium occurring within each region. Some linkage disequilibrium also exists between the two genes.

African dung beetle Onthophagus gazella Fabricius (Coleoptera: Scarabaeidae) esterase isozymes

African beetles Onthophagus gazella from both sexes were analyzed by electrophoresis for an investigation of esterase isozymes using alpha-naphthyl propionate and methylumbelliferyl propionate as substrates. Only one of the esterases (Est. 6) reacted with one of the substrates (alpha-naphthyl propionate). Six areas of activity were found, two of them being polymorphic (Est. 3 and Est. 4). For presence of Est. 3, 337 individuals were analyzed, including descendants of 32 controlled crossings: two alleles were identified, whose frequencies are Est. 3A = 0.447 and Est. 3B = 0.553. The population is in equilibrium for this locus (qui-square = 4.18; 0.2 > P > 0.1). For Est. 4, 338 individuals, descendants of 32 controlled crossings, were analysed. In this case, three alleles were identified whose frequencies are: Est. 4A = 0.277; Est. 4B = 0.661; and Est. 4C = 0.062. The population is not in equilibrium for this locus (qui-square = 40.259; p < 0.001). Two esterases were detected only in the pupal stage and another one in larvae. Of the 23 loci analyzed in these insects up to now, 3 are polymorphic (13%), which indicates very low variability in the population here studied.

Preliminary characterisation of esterase and platelet-activating factor (PAF)-acetylhydrolase activities from cat flea (Ctenocephalides felis) salivary glands

Naphthyl esterase and platelet-activating factor (PAF)-acetylhydrolase activities were detected in the salivary glands of the cat flea, Ctenocephalides felis. Salivary naphthyl esterase activity is disgorged during exploratory probing. Whole extracts of salivary glands contain esterase activity against the short-chain naphthyl esters alpha-naphthyl acetate (approximately 210pmol/min/gland pair; 10.0micromol/min/mg specific activity; K(m) approximately 59microM) and beta-naphthyl acetate (approximately 110pmol/min/gland pair; 5.2micromol/min/mg specific activity; K(m) approximately 132microM). Salivary gland extracts have PAF-acetylhydrolase activity (approximately 5pmol/min/gland pair; 0.24micromol/min/mg specific activity) but do not have detectable acetylcholinesterase activity. Native-PAGE and IEF resolve three and six salivary gland naphthyl esterase bands, respectively, and both patterns are different from carcass esterases. Salivary gland naphthyl esterase activity binds reversibly to Concanavalin A, and enzymatic deglycosylation with glycopeptidase F produced a new, fast-migrating salivary gland naphthyl esterase band on Native-PAGE. Renaturation of esterase activity after SDS-PAGE gave approximately 56kDa, approximately 57kDa and approximately 58kDa naphthyl-esterase-positive bands. On gel filtration naphthyl esterase and PAF-acetylhydrolase activities co-elute as a single peak with an apparent molecular weight of approximately 59kDa. This partially purified pool of enzyme had esterase activity against a series of short-chain alpha- and beta-naphthyl esters. The heterogeneity of salivary gland esterases, their relationship to PAF-acetylhydrolase, and the possible physiological functions of salivary gland PAF-acetylhydrolase activity are discussed.

Identification of a point mutation in an esterase gene in different populations of the southern cattle tick, Boophilus microplus

Two esterase cDNA sequences were obtained from susceptible and organophosphorus resistant strains of Boophilus microplus. Both sequences have a high degree of homology to carboxylesterase B. One gene has identical sequences in both strains and the other showed two point mutations. One mutation produces an amino acid substitution when the amino acid sequence is deduced, this mutation was detected in six different populations susceptible and resistant to insecticides, but a pyrethroid resistant strain was the only one that showed only the mutant allele. Identification of this mutation and the strong signal detected in southern blot with this strain, suggest that esterases are contributing to detoxification of pyrethroid compounds, as a resistant mechanism in Mexican strains of the southern cattle tick.

Carboxyl/cholinesterases: a case study of the evolution of a successful multigene family

The evolution of organismal diversity among the Metazoa is dependent on the proliferation of genes and diversification of functions in multigene families. Here we analyse these processes for one highly successful family, the carboxyl/cholinesterases. One key to the expansion of the functional niche of this group of enzymes is associated with versatile substrate binding and catalytic machinery. Qualitatively new functions can be obtained by substitution of one or a very few amino acids. This crudely adapted new functionality is then refined rapidly by a pulse of change elsewhere in the molecule; in one case about 13% amino acid divergence occurred in 5-10 million years. Furthermore, we postulate that the versatility of the substrate binding motifs underpins the recruitment of several family members to additional noncatalytic signal transduction functions.

Molecular evolution of two linked genes, Est-6 and Sod, in Drosophila melanogaster

We have obtained 15 sequences of Est-6 from a natural population of Drosophila melanogaster to test whether linkage disequilibrium exists between Est-6 and the closely linked Sod, and whether natural selection may be involved. An early experiment with allozymes had shown linkage disequilibrium between these two loci, while none was detected between other gene pairs. The Sod sequences for the same 15 haplotypes were obtained previously. The two genes exhibit similar levels of nucleotide polymorphism, but the patterns are different. In Est-6, there are nine amino acid replacement polymorphisms, one of which accounts for the S-F allozyme polymorphism. In Sod, there is only one replacement polymorphism, which corresponds to the S-F allozyme polymorphism. The transversion/transition ratio is more than five times larger in Sod than in Est-6. At the nucleotide level, the S and F alleles of Est-6 make up two allele families that are quite different from each other, while there is relatively little variation within each of them. There are also two families of alleles in Sod, one consisting of a subset of F alleles, and the other consisting of another subset of F alleles, designed F(A), plus all the S alleles. The Sod F(A) and S alleles are completely or nearly identical in nucleotide sequence, except for the replacement mutation that accounts for the allozyme difference. The two allele families have independent evolutionary histories in the two genes. There are traces of statistically significant linkage disequilibrium between the two genes that, we suggest, may have arisen as a consequence of selection favoring one particular sequence at each locus.

The same amino acid substitution in orthologous esterases confers organophosphate resistance on the house fly and a blowfly

Organophosphate (OP) insecticide resistance in certain strains of Musca domestica is associated with reduction in the carboxylesterase activity of a particular esterase isozyme. This has been attributed to a 'mutant ali-esterase hypothesis', which invokes a structural mutation to an ali-esterase resulting in the loss of its carboxylesterase activity but acquisition of OP hydrolase activity. It has been shown that the mutation in Lucilia cuprina is a Gly137-->Asp substitution in the active site of an esterase encoded by the Lc alpha E7 gene (Newcomb, R.D., Campbell, P.M., Ollis, D.L., Cheah, E., Russell, R.J., Oakeshott, J.G., 1997. A single amino acid substitution converts a carboxylesterase to an organophosphate hydrolase and confers insecticide resistance on a blowfly. Proc. Natl. Acad. Sci. USA 94, 7464-7468). We now report the cloning and characterisation of the orthologous M. domestica Md alpha E7 gene, including the sequencing of cDNAs from the OP resistant Rutgers and OP susceptible sbo and WHO strains. The Md alpha E7 gene has the same intron structure as Lc alpha E7 and encodes a protein with 76% amino acid identity to Lc alpha E7. Comparisons between susceptible and resistance alleles show resistance in M. domestica is associated with the same Gly137-->Asp mutation as in L. cuprina. Bacterial expression of the Rutgers allele shows its product has OP hydrolase activity. The data indicate identical catalytic mechanisms have evolved in orthologous Md alpha E7 and Lc alpha E7 molecules to endow diazinon-type resistance on the two species of higher Diptera.

Differential mRNA expression levels and gene sequences of a putative carboxylesterase-like enzyme from two strains of the parasitoid Anisopteromalus calandrae (Hymenoptera: Pteromalidae)

Carboxylesterase-like enzyme cDNAs have been cloned and sequenced from malathion-resistant and susceptible strains of the parasitoid Anisopteromalus calandrae (Howard) (Hymenoptera: Pteromalidae). The cDNAs consist of 1963 nucleotides including a 35 bp untranslated 5'-end, a 1596 bp open reading frame, and a 332 bp untranslated 3'-end. The open reading frame encodes 532 amino acid residues. The predicted protein sequence from these cDNAs includes 2 potential N-glycosylation sites, a carboxylesterase type-B serine active site FGGDSENVTIFGESAG, and conserved residues Ser187, Glu317, and His432 to function as the catalytic triad. The predicted carboxylesterase-like enzyme sequence is most similar to that of the carboxylesterase from the peach-potato aphid, Myzus persicae with 45% sequence identity. Alignment of the parasitoid carboxylesterase-like enzyme cDNAs revealed that there are two nucleotide differences in the open reading frame between the parasitoid strains, including a silent mutation and a point mutation that presumably causes a gene product difference. A nucleotide thymine at position 658 in the susceptible strain cDNA is replaced by a guanine in the resistant strain cDNA. This substitution leads to an amino acid change from tryptophan (Trp220) in the susceptible strain to glycine (Gly220) in the resistant strain. This substitution is genetically linked to resistance but it is not known how or if this amino acid substitution affects detoxification of malathion. Northern blot analyses demonstrated that expression level of the carboxylesterase-like enzyme mRNA in adult A. calandrae is approximately 30-fold higher in the resistant strain relative to that in the susceptible strain. Southern analysis indicated that Pst I or Eco RI restriction sites are different in the two strains. Both a modified gene structure and an increase in expression of carboxylesterase may be responsible for the high level of resistance found in this beneficial wasp.

cDNA cloning, baculovirus-expression and kinetic properties of the esterase, E3, involved in organophosphorus resistance in Lucilia cuprina

Resistance to organophosphorus insecticides (OPs) in the sheep blowfly, Lucilia cuprina, is associated with a non-staining phenotype of the carboxylesterase isozyme, E3 (E.C. 3.1.1.1). Here, we show that a member of alpha-esterase multigene family, Lc alpha E7, encodes E3. An Lc alpha E7 cDNA has been isolated from an OP-susceptible strain and expressed in a baculovirus. The expressed product is the same as E3 in its electrophoretic mobility and preference for alpha-over beta-naphthyl acetate as substrate. Its preference (kcat/K(m)) for a range of carboxylester substrates is alpha-naphthyl butyrate > alpha-naphthyl propionate > alpha-naphthyl acetate > methylthiobutyrate > p-nitrophenyl acetate. The enzyme is potently inhibited by OPs (ki [paraoxon] = 6.3 +/- 1.4 x 10(7)/M/min, ki [chlorfenvinphos] = 5.9 +/- 0.6 x 10(7)/M/min) and exhibits a high turnover of methylthiobutyrate (1009/s), consistent with its proposed homology to the ali-esterase that is thought to mutate to confer OP resistance in Musca domestica. E3 shares 64% amino acid identity with its Drosophila melanogaster homologue, Dm alpha E7, and is also closely related to other esterases involved in OP resistance such as the B1 esterase of Culex pipiens (38%) and E4 of Myzus persicae (30%).

Duplication and divergence of the genes of the alpha-esterase cluster of Drosophila melanogaster

The alpha-esterase cluster of D. melanogaster contains 11 esterase genes dispersed over 60 kb. Embedded in the cluster are two unrelated open reading frames that have sequence similarity with genes encoding ubiquitin-conjugating enzyme and tropomyosin. The esterase amino acid sequences show 37-66% identity with one another and all but one have all the motifs characteristic of functional members of the carboxyl/cholinesterase multigene family. The exception has several frameshift mutations and appears to be a pseudogene. Patterns of amino acid differences among cluster members in relation to generic models of carboxyl/cholinesterase protein structure are broadly similar to those among other carboxyl/cholinesterases sequenced to date. However the alpha-esterases differ from most other members of the family in: their lack of a signal peptide; the lack of conservation in cysteines involved in disulfide bridges; and in four indels, two of which occur in or adjacent to regions that align with proposed substrate-binding sites of other carboxyl/cholinesterases. Phylogenetic analyses clearly identify three simple gene duplication events within the cluster. The most recent event involved the pseudogene which is located in an intron of another esterase gene. However, relative rate tests suggest that the pseudogene remained functional after the duplication event and has become inactive relatively recently. The distribution of indels also suggests a deeper node in the gene phylogeny that separates six genes at the two ends of the cluster from a block of five in the middle.

Isolation of alpha cluster esterase genes associated with organophosphate resistance in Lucilia cuprina

PCR primers designed from the alpha-esterase gene cluster of Drosophila melanogaster have been used to isolate fragments from eight esterase genes in the Australian sheep blowfly, Lucilia cuprina. Phylogenetic analysis suggests that three are homologues of the alpha E7, alpha E8 and alpha E9 genes of the alpha-esterase cluster of D. melanogaster. A further three are also probably alpha-esterases, whereas the remaining two more closely resemble beta-esterases. Transcripts for five of the alpha-esterase genes were detected by PCR in adult midgut, consistent with a role for these enzymes in digestion and/or detoxification. Based on the tissue distribution of these transcripts, Lc alpha E7 may possibly encode the esterase, E3, which is involved in organophosphate resistance.

Excessive homoplasy in an evolutionarily constrained protein

The evolution of monomorphic proteins among closely related species has not been examined in detail. To investigate this phenomenon, the glycerol-3-phosphate dehydrogenase (Gpdh) locus was sequence in a broad range of Drosophila species. Although purifying selection to remove amino acid variation is the dominant force in the evolution of Gpdh, some replacements have occurred. The sequences were compared in the context of the phylogeny of the genus, revealing a high proportion of amino acid parallelism and reversal (homoplasy) at four sites. The level of homoplasy is significantly greater than that seen in other proteins for which multiple sequences are available, showing that Gpdh is strongly constrained by both the number of amino acid differences and the types of changes allowed. These four sites evolve at a much higher rate than do the other variable positions in the protein, accounting for half of the interspecific amino acid replacements. However, unlike typical hypervariable sites, where multiple changes to several different amino acids are seen, evolutionary 'flip-flopping' between two amino acid states defines this new class of hypervariable site.

Molecular cloning of an alpha-esterase gene cluster on chromosome 3r of Drosophila melanogaster

All or part of the alpha-esterase gene cluster in Drosophila melanogaster has been isolated by screening a YAC clone that spans cytological region 84D3-10 with consensus carboxyl/cholinesterase oligonucleotides. The cluster encompasses 11 putative esterase genes within 65 kb of genomic DNA and is one of the largest clusters of related protein-coding genes yet reported in Drosophila. The cluster must include the gene encoding the major alpha-esterase isozyme, EST9, which has previously been mapped to 84D3-5. It probably also includes the genes encoding the EST23, MCE and ALI esterases that have previously been mapped to 84D3-E2. The latter three are homologs of genes involved in organophosphate insecticide resistance in the sheep blowfly, Lucilia cuprina and the housefly, Musca domestica. Sequencing of one of the putative esterase genes in the Drosophila cluster, alpha E1, shows that it would encode features characteristic of an active carboxyl/cholinesterase, including the so-called catalytic triad, the nucleophilic elbow and oxyanion hole. It also shows that the closest relative of alpha E1 amongst previously published esterase sequences is ESTB1, which confers organophosphate resistance in Culex mosquitoes. We argue that we have cloned the D. melanogaster version of a major cluster of esterase genes which have variously mutated to confer organophosphate resistance in diverse Diptera.