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

Gut transcriptome of two bark beetle species stimulated with the same kairomones reveals molecular differences in detoxification pathways

Dendroctonus bark beetles are the most destructive agents in coniferous forests. These beetles come into contact with the toxic compounds of their host's chemical defenses throughout their life cycle, some of which are also used by the insects as kairomones to select their host trees during the colonization process. However, little is known about the molecular mechanisms by which the insects counteract the toxicity of these compounds. Here, two sibling species of bark beetles, D. valens and D. rhizophagus, were stimulated with vapors of a blend of their main kairomones (α-pinene, β-pinene and 3-carene), in order to compare the transcriptional response of their gut. A total of 48 180 unigenes were identified in D. valens and 43 704 in D. rhizophagus, in response to kairomones blend. The analysis of differential gene expression showed a transcriptional response in D. valens (739 unigenes, 0.58–10.36 Log2FC) related to digestive process and in D. rhizophagus (322 unigenes 0.87–13.08 Log2FC) related to xenobiotics metabolism. The expression profiles of detoxification genes mainly evidenced the up-regulation of COEs and GSTs in D. valens, and the up-regulation of P450s in D. rhizophagus. Results suggest that terpenes metabolism comes accompanied by an integral hormetic response, result of compensatory mechanisms, including the activation of other metabolic pathways, to ensure the supply of energy and the survival of organisms which is specific for each species, according to its life history and ecological strategy.

Genomic analysis of the carboxylesterase family in the salmon louse (Lepeophtheirus salmonis)

The pyrethroid deltamethrin and the macrocyclic lactone emamectin benzoate (EMB) are used to treat infestations of farmed salmon by parasitic salmon lice, Lepeophtheirus salmonis. While the efficacy of both compounds against Atlantic populations of the parasite has decreased as a result of the evolution of resistance, the molecular mechanisms of drug resistance in L. salmonis are currently not fully understood. The functionally diverse carboxylesterases (CaE) family includes members involved in pesticide resistance phenotypes of terrestrial arthropods. The present study had the objective to characterize the CaE family in L. salmonis and assess its role in drug resistance. L. salmonis CaE homologues were identified by homology searches in the parasite's transcriptome and genome. The transcript expression of CaEs predicted to be catalytically competent was studied using quantitative reverse-transcription PCR in drug susceptible and multi-resistant L. salmonis. The above strategy led to the identification of 21 CaEs genes/pseudogenes. Phylogenetic analyses assigned 13 CaEs to clades involved in neurodevelopmental signaling and cell adhesion, while three sequences were predicted to encode secreted enzymes. Ten CaEs were identified as being potentially catalytically competent. Transcript expression of acetylcholinesterase (ace1b) was significantly increased in multi-resistant lice compared to drug-susceptible L. salmonis, with transcript abundance further increased in preadult-II females following EMB exposure. In summary, results from the present study demonstrate that L. salmonis possesses fewer CaE gene family members than most arthropods characterized so far. Drug resistance in L. salmonis was associated with overexpression of ace1b.

Carboxylesterase genes in nitenpyram-resistant brown planthoppers, Nilaparvata lugens

Carboxylesterases (CarEs) represent one of the major detoxification enzyme families involved in insecticide resistance. However, the function of specific CarE genes in insecticide resistance is still unclear in the insect Nilaparvata lugens (Stål), a notorious rice crop pest in Asia. In this study, a total of 29 putative CarE genes in N. lugens were identified, and they were divided into seven clades; further, the β-esterase clade was significantly expanded. Tissue-specific expression analysis found that 17 CarE genes were abundantly distributed in the midgut and fat body, while 12 CarE genes were highly expressed in the head. The expression of most CarE genes was significantly induced in response to the challenge of nitenpyram, triflumezopyrim, chlorpyrifos, isoprocarb and etofenprox. Among these, the expression levels of NlCarE2, NlCarE4, NlCarE9, NlCarE17 and NlCarE24 were increased by each insecticide. Real-time quantitative polymerase chain reaction and RNA interference assays revealed the NlCarE1 gene to be a candidate gene mainly involved in nitenpyram resistance, while simultaneously silencing NlCarE1 and NlCarE19 produced a stronger effect than silencing either one individually, suggesting a cooperative relationship in resistance formation. These findings lay the foundation for further clarification of insecticide resistance mediated by CarE in N. lugens.

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.

Brown marmorated stink bug, Halyomorpha halys (Stål), genome: putative underpinnings of polyphagy, insecticide resistance potential and biology of a top worldwide pest

BACKGROUND

Halyomorpha halys (Stål), the brown marmorated stink bug, is a highly invasive insect species due in part to its exceptionally high levels of polyphagy. This species is also a nuisance due to overwintering in human-made structures. It has caused significant agricultural losses in recent years along the Atlantic seaboard of North America and in continental Europe. Genomic resources will assist with determining the molecular basis for this species' feeding and habitat traits, defining potential targets for pest management strategies.

RESULTS

Analysis of the 1.15-Gb draft genome assembly has identified a wide variety of genetic elements underpinning the biological characteristics of this formidable pest species, encompassing the roles of sensory functions, digestion, immunity, detoxification and development, all of which likely support H. halys' capacity for invasiveness. Many of the genes identified herein have potential for biomolecular pesticide applications.

CONCLUSIONS

Availability of the H. halys genome sequence will be useful for the development of environmentally friendly biomolecular pesticides to be applied in concert with more traditional, synthetic chemical-based controls.

Detoxification Genes Differ Between Cactus-, Fruit-, and Flower-Feeding Drosophila

We use annotated genomes of 14 Drosophila species covering diverse host use phenotypes to test whether 4 gene families that often have detoxification functions are associated with host shifts among species. Bark, slime flux, flower, and generalist necrotic fruit-feeding species all have similar numbers of carboxyl/cholinesterase, glutathione S-transferase, cytochrome P450, and UDP-glucuronosyltransferase genes. However, species feeding on toxic Morinda citrifolia fruit and the fresh fruit-feeding Drosophila suzukii have about 30 and 60 more, respectively. ABC transporters show a different pattern, with the flower-feeding D. elegans and the generalist necrotic fruit and cactus feeder D. hydei having about 20 and >100 more than the other species, respectively. Surprisingly, despite the complex secondary chemistry we find that 3 cactophilic specialists in the mojavensis species cluster have variably fewer genes than any of the other species across all 4 families. We also find 82 positive selection events across the 4 families, with the terminal D. suzukii and M. citrifolia-feeding D. sechellia branches again having the highest number of such events in proportion to their respective branch lengths. Many of the genes involved in these host-use-specific gene number differences or positive selection events lie in specific clades of the gene families that have been recurrently associated with detoxification. Several genes are also found to be involved in multiple duplication and/or positive selection events across the species studied regardless of their host use phenotypes; the most frequently involved are the ABC transporter CG1718, which is not in a specific clade associated with detoxification, and the α-esterase gene cluster, which is.

Genome-Wide Analysis of Carboxylesterases (COEs) in the Whitefly, Bemisia tabaci (Gennadius)

The whitefly (Bemisia tabaci), an important invasive pest that causes severe damage to crops worldwide, has developed resistance to a variety of insecticides. Carboxylesterases (COEs) are important multifunctional enzymes involved in the growth, development, and xenobiotic metabolism of insects. However, systematic studies on the COEs of B. tabaci are scarce. Here, 42 putative COEs in different functional categories were identified in the Mediterranean species of B. tabaci (B. tabaci MED) based on a genome database and neighbor-joining phylogeny. The expression patterns of the COEs were affected by the development of B. tabaci. The expression levels of six COEs were positively correlated with the concentration of imidacloprid to which B. tabaci adults were exposed. The mortality of B. tabaci MED adults fed dsBTbe5 (67.5%) and dsBTjhe2 (58.4%) was significantly higher than the adults fed dsEGFP (41.1%) when treated with imidacloprid. Our results provide a basis for functional research on COEs in B. tabaci and provide new insight into the imidacloprid resistance of B. tabaci.

Latrophilin mediates insecticides susceptibility and fecundity through two carboxylesterases, esterase4 and esterase6, in Tribolium castaneum

Latrophilin (LPH) is known as an adhesion G-protein-coupled receptor which involved in multiple physiological processes in organisms. Previous studies showed that lph not only involved the susceptibility to anticholinesterase insecticides but also affected fecundity in Tribolium castaneum. However, its regulatory mechanisms in these biological processes are still not clear. Here, we identified two potential downstream carboxylesterase (cce) genes of Tclph, esterase4 and esterase6, and further characterized their interactions with Tclph. After treatment of T. castaneum larvae with carbofuran or dichlorvos insecticides, the transcript levels of Tcest4 and Tcest6 were significantly induced from 12 to 72 h. RNAi against Tcest4 or Tcest6 led to the higher mortality compared with the controls after the insecticides treatment, suggesting that these two genes play a vital role in detoxification of insecticides in T. castaneum. Furthermore, with insecticides exposure to Tclph knockdown beetles, the expression of Tcest4 was upregulated but Tcest6 was downregulated, indicating that beetles existed a compensatory response against the insecticides. Additionally, RNAi of Tcest6 resulted in 43% reductions in female egg laying and completely inhibited egg hatching, which showed the similar phenotype as that of Tclph knockdown. These results indicated that Tclph affected fecundity by positively regulating Tcest6 expression. Our findings will provide a new insight into the molecular mechanisms of Tclph involved in physiological functions in T. castaneum.

Carboxylesterase genes in pyrethroid resistant house flies, Musca domestica

Carboxylesterases are one of the major enzyme families involved in the detoxification of pyrethroids. Up-regulation of carboxylesterase genes is thought to be a major component of insecticide resistant mechanisms in insects. Based on the house fly transcriptome and genome database, a total of 39 carboxylesterase genes of different functional clades have been identified in house flies. In this study, eleven of these genes were found to be significantly overexpressed in the resistant ALHF house fly strain compared with susceptible aabys and wild-type CS strains. Eight up-regulated carboxylesterase genes with their expression levels were further induced to a higher level in response to permethrin treatments, indicating that constitutive and inductive overexpression of carboxylesterases are co-responsible for the enhanced detoxification of insecticides. Spatial expression studies revealed these up-regulated genes to be abundantly distributed in fat bodies and genetically mapped on autosome 2 or 3 of house flies, and their expression could be regulated by factors on autosome 1, 2 and 5. Taken together, these results demonstrate that multiple carboxylesterase genes are co-upregulated in resistant house flies, providing further evidence for their involvement in the detoxification of insecticides and development of insecticide resistance.

Esterase Profile in Drosophila mercatorum pararepleta (Diptera; Drosophilidae), a Non-cactophilic Species of the repleta Group: Development Patterns and Aspects of Genetic Variability

Luciana Paes de Barros Machado, Natalia Silva Alves, Jaqueline de Oliveira Prestes, Gabriela Ronchi Salomón, Daiane Biegai, Thais Wouk, and Rogério Pincela Mateus (2017) Esterases are a diversified group of isozymes that performs several metabolic functions in Drosophila. In the D. repleta group, this class of enzymes was well described in cactophilic species, existing a lack of studies considering substrate speci city and life cycle expression in the non-cactophilic species. The larvae of cactophilic species of the D. repleta group develop in rotting cacti cladodes, but adults are generalists. Thus, different patterns expression can be found for esterases throughout development. In this work we analyzed esterase pro le and substrate speci city during development, and genetic variability aspects in D. mercatorum pararepleta, a non-cactophilic and generalist species of D. repleta group that was understudied hitherto. Samples of 3rd (F3) and 104th (F104) generations of three D. mercatorum pararepleta strains, obtained after collections in xerophytic enclaves of southeastern Brazil (ITI and SER in São Paulo state and RIP in Paraná state), and of D33 strain (obtained from Cristalina-GO, Midwest of Brazil, and established in the laboratory in 1987) were analyzed. Eight esterase loci, EST-1 to EST-8, were detected. EST-1 and EST-2 were adult exclusive. Only EST-3 and EST-8 were monomorphic; all the others presented between two (EST-6) and six (EST-7) alleles. EST-7 was the only dimeric locus and also the only one that showed to be a preferably β-esterase regarding affinity to α- and β-naphthyl acetates as substrates. The other seven loci were divided into three classes: α-esterase exclusive (EST-2); preferably α-esterase (EST-3, EST-4, EST-5 and EST-8); and α/β-esterase (EST-1 and EST-6). The EST-3, EST-5 and EST-6 loci were not detected in all samples, suggesting that they could have become pseudogenes due to the mutation accumulation after the gene duplication. The allele frequency of EST-7 locus, which showed the highest number of alleles, in adults of D33 and SER-F3 evidenced a higher variability and diversity in the oldest strain (six alleles, Ho = 0.46) than in the youngest ( five alleles, Ho = 0.26). Moreover, the analysis of SER-F104 revealed that this locus became monomorphic. The higher variability in the strain established in the laboratory at least two decades ago, together with the allele fixation in the SER-F104, indicate that the SER strain probably suffered a more severe action of founder effect/bottleneck when it was established in the laboratory and, therefore, even if the maintenance afterwards was performed using a high number of individuals, it did not assured the conservation of the existing genetic variability.

Natural genomic amplification of cholinesterase genes in animals

Tight control of the concentration of acetylcholine at cholinergic synapses requires precise regulation of the number and state of the acetylcholine receptors, and of the synthesis and degradation of the neurotransmitter. In particular, the cholinesterase activity has to be controlled exquisitely. In the genome of the first experimental models used (man, mouse, zebrafish and drosophila), there are only one or two genes coding for cholinesterases, whereas there are more genes for their closest relatives the carboxylesterases. Natural amplification of cholinesterase genes was first found to occur in some cancer cells and in insect species subjected to evolutionary pressure by insecticides. Analysis of the complete genome sequences of numerous representatives of the various metazoan phyla show that moderate amplification of cholinesterase genes is not uncommon in molluscs, echinoderms, hemichordates, prochordates or lepidosauria. Amplification of acetylcholinesterase genes is also a feature of parasitic nematodes or ticks. In these parasites, over-production of cholinesterase-like proteins in secreted products and the saliva are presumed to have effector roles related to host infection. These amplification events raise questions about the role of the amplified gene products, and the adaptation processes necessary to preserve efficient cholinergic transmission. This is an article for the special issue XVth International Symposium on Cholinergic Mechanisms.

Identification of carboxylesterase genes and their expression profiles in the Colorado potato beetle Leptinotarsa decemlineata treated with fipronil and cyhalothrin

Based on the Leptinotarsa decemlineata transcriptome dataset and the GenBank sequences, 70 novel carboxylesterases and 2 acetylcholinesterases were found. The 72 members belong to a multifunctional carboxylesterase/cholinesterase superfamily (CCE). A phylogenetic tree including the 72 LdCCEs and the CCEs from Tribolium castaneum, Drosophila melanogaster and Apis mellifera revealed that all CCEs fell into three main phylogenetic groups: dietary/detoxification, hormone/semiochemical processing, and neurodevelopmental classes. Numbers of L. decemlineata CCEs in the three classes were 52, 12 and 8, respectively. The dietary/detoxification class includes two clades: coleopteran xenobiotic metabolizing and α-esterase type CCEs. CCEs in the two clades have independently expanded in L. decemlineata. The hormone/semiochemical processing class has three clades: integument CCEs, β- and pheromone CCEs and juvenile hormone CCEs. Integument CCEs in L. decemlineata have also expanded. The neurodevelopmental CCEs are implicated the most ancient class, containing acetylcholinesterase, neuroligin, neurotactin, glutactin, gliotactin and others. Among the 70 novel CCE genes, KM220566, KM220530, KM220576, KM220527 and KM220541 were fipronil-inducible, and KM220578, KM220566, KM220542, KM220564, KM220561, KM220554, KM220527, KM220538 and KM220541 were cyhalothrin-inducible. They were the candidates involving in insecticide detoxification. Moreover, our results also provided a platform to understand the functions and evolution of L. decemlineata CCE genes.

The molecular evolution of cytochrome P450 genes within and between drosophila species

We map 114 gene gains and 74 gene losses in the P450 gene family across the phylogeny of 12 Drosophila species by examining the congruence of gene trees and species trees. Although the number of P450 genes varies from 74 to 94 in the species examined, we infer that there were at least 77 P450 genes in the ancestral Drosophila genome. One of the most striking observations in the data set is the elevated loss of P450 genes in the Drosophila sechellia lineage. The gain and loss events are not evenly distributed among the P450 genes-with 30 genes showing no gene gains or losses whereas others show as many as 20 copy number changes among the species examined. The P450 gene clades showing the fewest number of gene gain and loss events tend to be those evolving with the most purifying selection acting on the protein sequences, although there are exceptions, such as the rapid rate of amino acid replacement observed in the single copy phantom (Cyp306a1) gene. Within D. melanogaster, we observe gene copy number polymorphism in ten P450 genes including multiple cases of interparalog chimeras. Nonallelic homologous recombination (NAHR) has been associated with deleterious mutations in humans, but here we provide a second possible example of an NAHR event in insect P450s being adaptive. Specifically, we find that a polymorphic Cyp12a4/Cyp12a5 chimera correlates with resistance to an insecticide. Although we observe such interparalog exchange in our within-species data sets, we have little evidence of it between species, raising the possibility that such events may occur more frequently than appreciated but are masked by subsequent sequence change.

Purification and characterization of a specific late-larval esterase from two species of the Drosophila repleta group: contributions to understand its evolution

Background

After duplication, one copy of an original gene can become redundant and decay toward a pseudogene status or functionally diverge. Here, we performed the purification and biochemical characterization of EST-4 (a late larval β-esterase) from two Drosophila repleta group species, Drosophila mulleri and Drosophila arizonae, in order to establish comparative parameters between these enzymes in these species and to contribute to better understand their evolution.

Results

In D. mulleri, EST-4 had an optimal activity in temperatures ranging from 40° to 45°C and at pH 7.5, maintaining stability in alkaline pH (8.0 to 10.0). It was classified as serine esterase as its activity was inhibited by PMSF. No ion negatively modulated EST-4 activity, and iron had the most positive modulating effect. In D. arizonae, it showed similar optimum temperature (40°C), pH (8.0), and was also classified as a serine esterase, but the enzymatic stability was maintained in an acidic pH (5.5 to 6.5). Fe+2 had the opposite effect found in D. mulleri, that is, negative modulation. Al+3 almost totally inhibited the EST-4 activity, and Na+ and Cu+2 had a positive modulation effect. Kinetic studies, using ρ-nitrophenyl acetate as substrate, showed that EST-4 from D. mulleri had higher affinity, while in D. arizonae, it showed higher V max and catalytic efficiency in optimal reaction conditions.

Conclusions

EST-4 from D. mulleri and D. arizonae are very closely related and still maintain several similar features; however, they show some degree of differentiation. Considering that EST-4 from D. mulleri has more conspicuous gel mobility difference among all EST-4 studied so far and a lower catalytic efficiency was observed here, we proposed that after duplication, this new copy of the original gene became redundant and started to decay toward a pseudogene status in this species, which probably is not occurring in D. arizonae.

Overexpressed esterases in a fenvalerate resistant strain of the cotton bollworm, Helicoverpa armigera

Enhanced detoxification is the major mechanism responsible for pyrethroid resistance in Chinese populations of Helicoverpa armigera. Previous work has shown that enhanced oxidation contributes to resistance in the fenvalerate-selected Chinese strain, YGF. The current study provides evidence that enhanced hydrolysis by esterase isozymes also contributes to the resistance in this strain. The average esterase activity of third instar YGF larvae was 1.9-fold compared with that of a susceptible SCD strain. Much of this difference was attributed to isozymes at two zones which hydrolysed the model carboxylester substrate 1-naphthyl acetate and also a 1-naphthyl analogue of fenvalerate. A preparation enriched for enzymes migrating to one of these zones from YGF was shown to hydrolyse fenvalerate with a specific activity of about 2.9 nmol/min/mg. This material was also matched by mass spectrometry with four putative carboxylesterase genes, all of which clustered within a phylogenetic clade of secreted midgut esterases. Quantitative PCR on these four genes showed several-fold greater expression in tissues of YGF compared to SCD but no differences was found in the number of copies of the genes between the strains.

Genomic analysis of carboxyl/cholinesterase genes in the silkworm Bombyx mori

Background

Carboxyl/cholinesterases (CCEs) have pivotal roles in dietary detoxification, pheromone or hormone degradation and neurodevelopment. The recent completion of genome projects in various insect species has led to the identification of multiple CCEs with unknown functions. Here, we analyzed the phylogeny, expression and genomic distribution of 69 putative CCEs in the silkworm, Bombyx mori (Lepidoptera: Bombycidae).

Results

A phylogenetic tree of CCEs in B. mori and other lepidopteran species was constructed. The expression pattern of each B. mori CCE was also investigated by a search of an expressed sequence tag (EST) database, and the relationship between phylogeny and expression was analyzed. A large number of B. mori CCEs were identified from a midgut EST library. CCEs expressed in the midgut formed a cluster in the phylogenetic tree that included not only B. mori genes but also those of other lepidopteran species. The silkworm, and possibly also other lepidopteran species, has a large number of CCEs, and this might be a consequence of the large cluster of midgut CCEs. Investigation of intron-exon organization in B. mori CCEs revealed that their positions and splicing site phases were strongly conserved. Several B. mori CCEs, including juvenile hormone esterase, not only showed clustering in the phylogenetic tree but were also closely located on silkworm chromosomes. We investigated the phylogeny and microsynteny of neuroligins in detail, among many CCEs. Interestingly, we found the evolution of this gene appeared not to be conserved between B. mori and other insect orders.

Conclusions

We analyzed 69 putative CCEs from B. mori. Comparison of these CCEs with other lepidopteran CCEs indicated that they had conserved expression and function in this insect order. The analyses showed that CCEs were unevenly distributed across the genome of B. mori and suggested that neuroligins may have a distinct evolutionary history from other insect order. It is possible that such an uneven genomic distribution and a unique neuroligin evolution are shared with other lepidopteran insects. Our genomic analysis has provided novel information on the CCEs of the silkworm, which will be of value to understanding the biology, physiology and evolution of insect CCEs.

Metabolic enzymes associated with xenobiotic and chemosensory responses in Nasonia vitripennis

The numbers of glutathione S-transferase, cytochrome P450 and esterase genes in the genome of the hymenopteran parasitoid Nasonia vitripennis are about twice those found in the genome of another hymenopteran, the honeybee Apis mellifera. Some of the difference is associated with clades of these families implicated in xenobiotic resistance in other insects and some is in clades implicated in hormone and pheromone metabolism. The data support the hypothesis that the eusocial behaviour of the honeybee and the concomitant homeostasis of the nest environment may obviate the need for as many gene/enzyme systems associated with xenobiotic metabolism as are found in other species, including N. vitripennis, that are thought to encounter a wider range of potentially toxic xenobiotics in their diet and habitat.

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.