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

Genomic and transcriptomic analyses in Drosophila suggest that the ecdysteroid kinase-like (EcKL) gene family encodes the 'detoxification-by-phosphorylation' enzymes of insects

Phosphorylation is a phase II detoxification reaction that, among animals, occurs near exclusively in insects, but the enzymes responsible have never been cloned or otherwise identified. We propose the hypothesis that members of the arthropod-specific ecdysteroid kinase-like (EcKL) gene family encode detoxicative kinases. To test this hypothesis, we annotated the EcKL gene family in 12 species of Drosophila and explored their evolution within the genus. Many ancestral EcKL clades are evolutionarily unstable and have experienced repeated gene gain and loss events, while others are conserved as single-copy orthologs. Leveraging multiple published gene expression datasets from D. melanogaster, and using the cytochrome P450s-a classical detoxification family-as a test case, we demonstrate relationships between xenobiotic induction, detoxification tissue-enriched expression and evolutionary instability in the EcKLs and the P450s. We devised a systematic method for identifying candidate detoxification genes in large gene families that is concordant with experimentally determined functions of P450 genes in D. melanogaster. Applying this method to the EcKLs suggested a significant proportion of these genes play roles in detoxification, and that the EcKLs may constitute a detoxification gene family in insects. Additionally, we estimate that between 11 and 16 uncharacterised D. melanogaster P450s are strong detoxification candidates. Lastly, we also found previously unreported genomic and transcriptomic variation in a number of EcKLs and P450s associated with toxic stress phenotypes using a targeted phenome-wide association study (PheWAS) approach in D. melanogaster, presenting multiple future avenues of research for detoxification genetics in this species.

Incongruent nuclear and mitochondrial genetic structure of new world screwworm fly populations due to positive selection of mutations associated with dimethyl- and diethyl-organophosphates resistance

Livestock production is an important economic activity in Brazil, which has been suffering significant losses due to the impact of parasites. The New World screwworm (NWS) fly, Cochliomyia hominivorax, is an ectoparasite and one of the most important myiasis-causing flies endemic to the Americas. The geographic distribution of NWS has been reduced after the implementation of the Sterile Insect Technique (SIT), being eradicated in North America and part of Central America. In South America, C. hominivorax is controlled by chemical insecticides, although indiscriminate use can cause selection of resistant individuals. Previous studies have associated the Gly137Asp and Trp251Leu mutations in the active site of carboxylesterase E3 to resistance of diethyl and dimethyl-organophosphates insecticides, respectively. Here, we have sequenced a fragment of the carboxylesterase E3 gene (ChαE7), comprising part of intron iII, exon eIII, intron iIII and part of exon eIV, and three mitochondrial gene sequences (CR, COI and COII), of NWS flies from 21 locations in South America. These markers were used for population structure analyses and the ChαE7 gene was also investigated to gain insight into the selective pressures that have shaped its evolution. Analysis of molecular variance (AMOVA) and pairwise FST analysis indicated an increased genetic structure between locations in the ChαE7 compared to the concatenated mitochondrial genes. Discriminant analysis of principal components (DAPC) and spatial analysis of molecular variance (SAMOVA) indicated different degrees of genetic structure for all markers, in agreement with the AMOVA results, but with low correlation to geographic data. The NWS fly is considered a panmitic species based on mitochondrial data, while it is structured into three groups considering the ChαE7 gene. A negative association between the two mutations related to organophosphate resistance and Fay & Wu’s H significant negative values for the exons, suggest that these mutations evolved under positive selection.

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.

Functional fat body proteomics and gene targeting reveal in vivo functions of Drosophila melanogaster α-Esterase-7

Carboxylesterases constitute a large enzyme family in insects, which is involved in diverse functions such as xenobiotic detoxification, lipid metabolism and reproduction. Phylogenetically, many insect carboxylesterases are represented by multienzyme clades, which are encoded by evolutionarily ancient gene clusters such as the α-Esterase cluster. Much in contrast to the vital importance attributed to carboxylesterases in general, the in vivo function of individual α-Esterase genes is largely unknown. This study employs a functional proteomics approach to identify esterolytic enzymes of the vinegar fly Drosophila melanogaster fat body. One of the fat body carboxylesterases, α-Esterase-7, was selected for mutational analysis by gene targeting to generate a deletion mutant fly. Phenotypic characterization of α-Esterase-7 null mutants and transgenic flies, which overexpress a chimeric α-Esterase-7:EGFP gene, reveals important functions of α-Esterase-7 in insecticide tolerance, lipid metabolism and lifespan control. The presented first deletion mutant of any α-Esterase in the model insect D. melanogaster generated by gene targeting not only provides experimental evidence for the endogenous functions of this gene family. It also offers an entry point for in vivo structure-function analyses of α-Esterase-7, which is of central importance for naturally occurring insecticide resistance in wild populations of various dipteran insect species.

Selective sweeps at the organophosphorus insecticide resistance locus, Rop-1, have affected variation across and beyond the α-esterase gene cluster in the Australian sheep blowfly, Lucilia cuprina

A major theoretical consequence of selection at a locus is the genetic hitchhiking of linked sites (selective sweep). The extent of hitchhiking around a gene is related to the strength of selection and the rate of recombination, with its impact diminishing with distance from the selected site. At the Rop-1 locus of the sheep blowfly, Lucilia cuprina, polymorphisms at two different sites within the LcαE7 gene encode forms of the protein that confer organophosphorus insecticide resistance. To assess the impact of selection at these two sites on variation around LcαE7, we sequenced regions within six other genes along chromosome IV across isogenic (IV) strains of L. cuprina. High levels of linkage disequilibrium, characterized by low haplotype number (K) and diversity (H), and significant R(2) values were observed for two genes, LcαE1 and LcαE10, both members of the same α-esterase gene cluster as LcαE7. A significant R(2) value was also observed for a gene predicted to be the next closest to LcαE7, AL03, but not for any of the other genes, LcRpL13a, Lcdsx, or LcAce. Skews in the site frequency spectra toward high-frequency variants were significant for LcαE1 (Fay and Wu's H = -2.91), LcαE10 (H = -1.85), and Lcdsx (H = -2.00). Since the selective sweeps, two forms of likely returning variation were observed, including variation in microsatellites in an intron of LcαE10 and a recombination event between LcαE7 and LcαE10. These data suggest that two incomplete soft sweeps have occurred at LcαE7 that have significantly affected variation across, and beyond, the α-esterase gene cluster of L. cuprina. The speed and impact of these selective sweeps on surrounding genomic variation and the ability of L. cuprina to respond to future environmental challenges are discussed.

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.

Gene family evolution across 12 Drosophila genomes

Comparison of whole genomes has revealed large and frequent changes in the size of gene families. These changes occur because of high rates of both gene gain (via duplication) and loss (via deletion or pseudogenization), as well as the evolution of entirely new genes. Here we use the genomes of 12 fully sequenced Drosophila species to study the gain and loss of genes at unprecedented resolution. We find large numbers of both gains and losses, with over 40% of all gene families differing in size among the Drosophila. Approximately 17 genes are estimated to be duplicated and fixed in a genome every million years, a rate on par with that previously found in both yeast and mammals. We find many instances of extreme expansions or contractions in the size of gene families, including the expansion of several sex- and spermatogenesis-related families in D. melanogaster that also evolve under positive selection at the nucleotide level. Newly evolved gene families in our dataset are associated with a class of testes-expressed genes known to have evolved de novo in a number of cases. Gene family comparisons also allow us to identify a number of annotated D. melanogaster genes that are unlikely to encode functional proteins, as well as to identify dozens of previously unannotated D. melanogaster genes with conserved homologs in the other Drosophila. Taken together, our results demonstrate that the apparent stasis in total gene number among species has masked rapid turnover in individual gene gain and loss. It is likely that this genomic revolving door has played a large role in shaping the morphological, physiological, and metabolic differences among species.

Though comparative genome sequencing has revealed vast similarities in the total number of genes contained within closely related species, this similarity hides enormous complexities in the identity and number of constituent proteins. Species can differ in their complement of genes through both gene duplication and loss. Here we investigated the gain and loss of genes from the genomes of 12 fully sequenced Drosophila (fruit flies). We find high rates of gain and loss in all species and estimate that approximately one new gene is gained or lost every 60,000 years. We also find several hundred cases of extremely rapid gene turnover, with dozens of genes gained or lost in only a few million years. The highest turnover in gene number occurs in genes involved in sex and reproduction. Taken together, our results demonstrate that the apparent stasis in total gene number among species has masked rapid turnover in individual gene gain and loss. It is likely that this evolutionary revolving door has played a large role in shaping the morphological, physiological, and metabolic differences among species.

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.

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.

Exploring the evolutionary history of the alcohol dehydrogenase gene (Adh) duplication in species of the family tephritidae

In the olive fruit fly Bactrocera oleae and the med fly Ceratitis capitata previous studies have shown the existence of two Adh genes in each species. This observation, in combination with the former finding that various Drosophila species of virilis and repleta group encode two isozymes of ADH which are the result of a gene duplication, challenged us to address a scenario dealing with the evolutionary history of the Adh gene duplication in Tephritidae. In our lab we proceeded to the cloning and sequence analysis of Adh genes from more tephritid species, a prerequisite for further study of this issue. Here we show that phylogenetic trees produced from either the nucleotide or the amino acid sequences of 14 tephritid Adh genes consisted of two main clusters, with Adh sequences of the same "type" grouping together (i.e., Adh1 sequences form a cluster and Adh2 sequences form a second one), as expected if there was one duplication event before speciation within the family Tephritidae. We used the amount of divergence between the two isozymic forms of Adh of the species carrying both Adh1 and Adh2 genes to obtain an estimate of the age of the duplication event. Interestingly, our data again support the hypothesis that the duplication of an ancestral Adh single gene in the family Tephritidae occurred before the emergence of the genera Bactrocera and Ceratitis, thus suggesting that Adh duplication was based on a prespeciation rather than a postspeciation event that might have involved two independent duplication events, one in each of the two genera.

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.

Expression of enhancers is altered in Drosophila melanogaster hybrids

The molecular foundations of evolution are difficult to trace because most protein sequences are virtually identical in closely related species. The largest fraction of sequence within the genome, however, is composed of noncoding sequences where regulatory elements locate to various sites. It has been suggested that changes in the activity of these elements may trigger evolutionary change. In Drosophila, the enhancer trap procedure identifies regulatory sequences in the genome after the insertion of a P-element-based construct. We generated new insertions and characterized their expression domains in the adult eye and larval imaginal disks using the white and LacZ reporter genes. Lines with robust expression patterns in D. melanogaster were analyzed in hybrids to test the conservation of regulatory mechanisms between species. Most of the enhancers used in this study modified their expression in hybrids with the mating species D. mauritiana and D. simulans. Expression changes resulted either in gain or loss of expression and were cell-type or hybrid-genome specific. Further characterization of a limited number of enhancers in D. melanogaster showed that expression domains could adapt to changes in cell number during development but not after the completion of cell proliferation. Also, expression of some enhancers appeared to be sensitive to heterochromatin from the Y but not the X chromosome. Taken together, these results demonstrate the high sensitivity of regulatory mechanisms of gene expression as a prime source of evolutionary change and suggest quantitative changes in available transcription factors as one of the mechanisms involved.

Overproduction of a P450 that metabolizes diazinon is linked to a loss-of-function in the chromosome 2 ali-esterase (MdalphaE7) gene in resistant house flies

Up-regulation of detoxifying enzymes in insecticide-resistant strains of the house fly is a common mechanism for metabolic resistance. However, the molecular basis of this increased insecticide metabolism is not well understood. In the multiresistant Rutgers strain, several cytochromes P450 and glutathione S-transferases are constitutively overexpressed at the transcriptional level. Overexpression is the result of trans-regulation, and a regulatory gene has been located on chromosome 2. A Gly137 to Asp point mutation in alphaE7 esterase gene, leading to the loss of carboxylesterase activity, has been associated with organophosphate resistance in the house fly and the sheep blowfly. We show here that purified recombinant CYP6A1 is able to detoxify diazinon with a high efficiency. We also show that either the Gly137 to Asp point mutation in alphaE7 esterase gene or a deletion at this locus confer resistance and overproduction of the CYP6A1 protein. Based on these findings, we propose it is the absence of the wild-type Gly137 allele of the alphaE7 gene that releases the transcriptional repression of genes coding for detoxification enzymes such as CYP6A1, thereby leading to metabolic resistance to diazinon.

Molecular and biochemical survey of acaricide resistance mechanisms in larvae from Mexican strains of the southern cattle tick, Boophilus microplus

We examined the larvae of several organophosphate and pyrethroid-resistant Mexican strains of Boophilus microplus using biochemical and molecular tests to investigate the mechanisms conferring acaricide resistance. The electrophoretic profiles of esterase activity in protein extracts from coumaphos and permethrin-resistant strains compared to the susceptible strain revealed distinct differences, which inhibitor studies attributed to carboxylesterases. Esterase hydrolysis assays showed significant enhancement of both total and permethrin hydrolysis in one pyrethroid-resistant strain, with no enhancement in two other strains with very high resistance to pyrethroids. Sequence analysis of sodium channel mRNA fragments in all pyrethroid-resistant strains determined that they did not possess the classic kdr and super-kdr mutations known to confer pyrethroid resistance in several insect species. Using reverse transcriptase-polymerase chain reaction (RT-PCR) with degenerate primers designed from conserved regions of insect esterase amino acid sequences, a B. microplus larval cDNA fragment was isolated whose deduced amino acid sequence was significantly similar to esterases from a wide range of species. In Northern blot RNA analysis the cDNA hybridized to a 2.1 kb mRNA that was abundant in all resistant strains except one, in which a very low abundance could provide a marker for the mechanism conferring resistance in this strain.

The evolution of an alpha-esterase pseudogene inactivated in the Drosophila melanogaster lineage

Previous analyses of the alpha-esterase cluster of Drosophila melanogaster revealed 10 active genes and the DmalphaE4a-Psi pseudogene. Here, we reconstruct the evolution of the pseudogene from the sequences of 12 alleles from widely scattered D. melanogaster populations and single alleles from Drosophila simulans and Drosophila yakuba. All of the DmalphaE4a-Psi alleles contain numerous inactivating mutations, suggesting that pseudogene alleles are fixed in natural populations. Several lines of evidence also suggest that DmalphaE4a is now evolving without selective constraint in the D. melanogaster lineage. There are three polymorphic indels which result in frameshifts; a key nucleotide of the intron splice acceptor is polymorphic; the neutral mutation parameter is the same for replacement and silent sites; one of the nonsilent polymorphisms results in a stop codon; only 1 of the 13 replacement polymorphisms is biochemically conservative; residues that are conserved among active esterases have different states in DmalphaE4a-Psi; and there are about half as many transitional polymorphisms as transversional ones. In contrast, the D. simulans and D. yakuba orthologs DsalphaE4a and DyalphaE4a do not have the inactivating mutations of DmalphaE4a-Psi and appear to be evolving under the purifying selection typical of protein- encoding genes. For instance, there have been more substitutions in the introns than in the exons, and more in silent sites than in replacement sites. Furthermore, most of the amino acid substitutions that have occurred between DyalphaE4a and DsalphaE4a are located in sites that typically vary among active alpha-esterases rather than those that are usually conserved. We argue that the original alphaE4a gene had a function which it has lost since the divergence of the D. melanogaster and D. simulans lineages.

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.

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.

Amyrel, a paralogous gene of the amylase gene family in Drosophila melanogaster and the Sophophora subgenus

We describe a gene from Drosophila melanogaster related to the alpha-amylase gene Amy. This gene, which exists as a single copy, was named Amyrel. It is strikingly divergent from Amy because the amino acid divergence is 40%. The coding sequence is interrupted by a short intron at position 655, which is unusual in amylase genes. Amyrel has also been cloned in Drosophila ananassae, Drosophila pseudoobscura, and Drosophila subobscura and is likely to be present throughout the Sophophora subgenus, but, to our knowledge, it has not been detected outside. Unexpectedly, there is a strong conservation of 5' and 3' flanking regions between Amyrel genes from different species, which is not the case for Amy and which suggests that selection acts on these regions. In contrast to the Amy genes, Amyrel is transcribed in larvae of D. melanogaster but not in adults. However, the protein has not been detected yet. Amyrel evolves about twice as fast as Amy in the several species studied. We suggest that this gene could result from a duplication of Amy followed by accelerated and selected divergence toward a new adaptation.

Evidence that the E4 and FE4 esterase genes responsible for insecticide resistance in the aphid Myzus persicae (Sulzer) are part of a gene family

The amplification of genes encoding the esterases E4 and FE4 is a widespread mechanism of insecticide resistance in the peach-potato aphid, Myzus persicae (Sulzer). We present evidence that in susceptible aphids the two genes are adjacent to each other in a head-to-tail arrangement with E4 upstream of FE4 and with approx. 19 kb of intervening sequence. There are also at least two other closely related sequences which might come from other members of an esterase gene family, in line with reports of other insect gene families encoding detoxifying enzymes. The close identity between E4 and FE4 genes indicates a recent duplication and divergence. The subsequent amplifications giving multiple copies of either E4 or FE4 must have involved two separate events, each probably occurring once and then being selected by insecticide exposure and spread by migration. The cloning of sequences upstream of the FE4 gene suggest, by comparison with E4, that the two genes are regulated in different ways. FE4 has sequences corresponding to a conventional promoter (TATA box and CAP site) that are not present in E4; on the other hand, FE4 lacks the CpG island present 5' of E4 genes that may control expression through changes in DNA methylation. The differences are likely to have occurred by the duplication event that gave rise to E4 and FE4 leading to different 5' sequences.

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%).