An evolutionary model for the duplication and divergence of esterase genes in Drosophila

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.

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.

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.

Positive and negative selection in the beta-esterase gene cluster of the Drosophila melanogaster subgroup

We examine the pattern of molecular evolution of the beta-esterase gene cluster, including the Est-6 and psiEst-6 genes, in eight species of the Drosophila melanogaster subgroup. Using maximum likelihood estimates of nonsynonymous/synonymous rate ratios, we show that the majority of Est-6 sites evolves under strong (48% of sites) or moderate (50% of sites) negative selection and a minority of sites (1.5%) is under significant positive selection. Est-6 sites likely to be under positive selection are associated with increased intraspecific variability. One positively selected site is responsible for the EST-6 F/S allozyme polymorphism; the same site is responsible for the EST-6 functional divergence between species of the melanogaster subgroup. For psiEst-6 83.7% sites evolve under negative selection, 16% sites evolve neutrally, and 0.3% sites are under positive selection. The positively selected sites of psiEst-6 are located at the beginning and at the end of the gene, where there is reduced divergence between D. melanogaster and D. simulans; these regions of psiEst-6 could be involved in regulation or some other function. Branch-site-specific analysis shows that the evolution of the melanogaster subgroup underwent episodic positive selection. Collating the present data with previous results for the beta-esterase genes, we propose that positive and negative selection are involved in a complex relationship that may be typical of the divergence of duplicate genes as one or both duplicates evolve a new function.

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.

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.

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.

The role of gene conversion in determining sequence variation and divergence in the Est-5 gene family in Drosophila pseudoobscura

Nucleotide sequences of eight Est-5A and Est-5C genes corresponding to previously sequenced Est-5B genes in Drosophila pseudoobscura were determined to compare patterns of polymorphism and divergence among members of this small gene family. The three esterase genes were also sequenced from D. persimilis and D. miranda for interspecific comparisons. The data provide evidence that gene conversion between loci contributes to polymorphism and to the homogenization of the Est5 genes. For Est-5B, which encodes one of the most highly polymorphic proteins in Drosophila, 12% of the segregating amino acid variants appear to have been introduced via gene conversion from other members of the gene family. Interlocus gene conversion can also explain high sequence similarity, especially at synonymous sites, between Est-5B and Est-5A. Tests of neutrality using interspecific comparisons show that levels of polymorphism conform to neutral expectations at each Est-5 locus. However, McDonald-Kreitman tests based on intraspecific gene comparisons indicate that positive selection on amino acids has accompanied Est-5 gene duplication and divergence in D. pseudoobscura.

Divergent and conserved features in the spatial expression of the Drosophila pseudoobscura esterase-5B gene and the esterase-6 gene of Drosophila melanogaster

The regulatory regions of homologous genes encoding esterase 6 (Est-6) of Drosophila melanogaster and esterase 5B (Est-5B) of Drosophila pseudoobscura show very little similarity. We have undertaken a comparative study of the pattern of expression directed by the Est-5B and Est-6 5'-flanking DNA to attempt to reveal conserved elements regulating tissue-specific expression in adults. Esterase regulatory sequences were linked to a lacZ reporter gene and transformed into D. melanogaster embryos. Est-5B, 5' upstream elements, give rise to a beta-galactosidase expression pattern that coincides with the wild-type expression of Est-5B in D. pseudoobscura. The expression patterns of the Est-5B/lacZ construct are different from those of a fusion gene containing the upstream region of Est-6. Common sites of expression for both kinds of constructs are the third segment of antenna, the maxillary palps, and salivary glands. In vitro deletion mutagenesis has shown that the two genes have a different organization of regulatory elements controlling expression in both the third segment of antenna and maxillary palps. The results suggest that the conservation of the expression pattern in genes that evolved from a common ancestor may not be accompanied by preservation of the corresponding cis-regulatory elements.

Esterase polymorphism in insecticide susceptible populations of the mosquito Culex pipiens

Gene amplification involving a particular haplotype has been found at the esterase B locus of mosquitoes from various countries. This similarity has been explained by a unique amplification event followed by migration and selection by organophosphate (OP) insecticides. This assumes that the polymorphism of non-amplified esterase haplotypes is so large that the chance of independent amplification in two distinct populations is negligible. In order to test this assumption, three susceptible populations from northern Europe were sampled and analysed for esterase and haplotype polymorphism. At the protein level, 18 and 16 alleles were found for esterase A and B respectively in one French population (n = 74), and 16 and 14 in an English one (n = 50). At the DNA level, 24 alleles at the esterase B locus were detected in a sample of 72 mosquitoes from one population, with the use of only one restriction enzyme (EcoR V). Restriction maps of two non-amplified haplotypes randomly sampled from a single breeding site in Belgium were built with six restriction enzymes. 60% of all restriction sites were different among the two maps. The huge polymorphism found in northern Europe requires specific explanations for its stability, but it considerably strengthens the hypothesis of migration of amplified haplotypes.

Nucleotide sequence of a novel arylesterase gene from Vibro mimicus and characterization of the enzyme expressed in Escherichia coli

A gene coding for an arylesterase of Vibrio mimicus was cloned. Sequence determination reveals that the esterase gene has an open reading frame of 600 nucleotides which encodes a protein of M(r) 22,300. The deduced amino acid sequence contain a pentapeptide GDSLS (residues 27-31), which was also found in the phospholipid-cholesterol acyltransferase from Aeromonas hydrophila. Substitution of Ser-29 by alanine or cysteine in the cloned gene abolished the esterase activity in the tributyrin plate assay. On the other hand, the activity was not lost when Ser-31 was changed to alanine. The cloned gene was expressed in Escherichia coli, and the protein purified by a four-step procedure. The purified protein migrated on SDS/PAGE as a single band with an apparent M(r) of 22,100. This enzyme favoured the hydrolysis of several arylesters and was classified as an arylesterase (EC 3.1.1.2). N-Terminal analysis showed that Ser-20 was the first amino acid of the mature secreted protein, suggesting that the N-terminal 19 hydrophobic amino acids served as a signal peptide.

Sequential gel electrophoretic analysis of esterase-2 in two populations of Drosophila buzzatii

Sequential electrophoresis, using three different buffer systems on cellulose acetate gels, was used to characterize the allelic variation for esterase-2 in two populations of D. buzzatii in Australia that are separated by 550 km. Twenty-five alleles were detected, of which nine were unique to one population, eight unique to the other, and only eight were common to both populations. Allele frequencies within each population were significantly different between the two major chromosome sequences (standard and j inversion), and for each chromosome sequence allele frequencies were significantly different between populations. Observed allelic frequency distributions were not significantly different from those predicted for selective neutrality using the homozygosity test statistic. However, estimates of the effective sizes of the populations derived from their observed differentiation, together with the history of the species in Australia, provide support for some form of balancing selection affecting at least some of the alleles.

Evolutionary genetics of Drosophila esterases

Over 30 carboxylester hydrolases have been identified in D. melanogaster. Most are classified as acetyl, carboxyl or cholinesterases. Sequence similarities among most of the carboxyl and all the cholinesterases so far characterised from D. melanogaster and other eukaryotes justify recognition of a carboxyl/cholinesterase multigene family. This family shows minimal sequence similarities with other esterases but crystallographic data for a few non-drosophilid enzymes show that the family shares a distinctive overall structure with some other carboxyl and aryl esterases, so they are all put in one superfamily of/beta hydrolases. Fifteen esterase genes have been mapped in D. melanogaster and twelve are clustered at two chromosomal sites. The constitution of each cluster varies across Drosophila species but two carboxyl esterases in one cluster are sufficiently conserved that their homologues can be identified among enzymes conferring insecticide resistance in other Diptera. Sequence differences between two other esterases, the EST6 carboxyl esterase and acetylcholinesterase, have been interpreted against the consensus super-secondary structure for the carboxyl/cholinesterase multigene family; their sequence differences are widely dispersed across the structure and include substantial divergence in substrate binding sites and the active site gorge. This also applies when EST6 is compared across species where differences in its expression indicate a difference in function. However, comparisons within and among species where EST6 expression is conserved show that many aspects of the predicted super-secondary structure are tightly conserved. Two notable exceptions are a pair of polymorphisms in the substrate binding site of the enzyme in D. melanogaster. These polymorphisms are associated with differences in substrate interactions in vitro and demographic data indicate that the alternative forms are not selectively equivalent in vivo.

Evolutionary genetics of the Drosophila alcohol dehydrogenase gene-enzyme system

Evolutionary genetics embodies a broad research area that ranges from the DNA level to studies of genetic aspects in populations. In all cases the purpose is to determine the impact of genetic variation on evolutionary change. The broad range of evolutionary genetics requires the involvement of a diverse group of researchers: molecular biologists, (population) geneticists, biochemists, physiologists, ecologists, ethologists and theorists, each of which has its own insights and interests. For example, biochemists are often not concerned with the physiological function of a protein (with respect to pH, substrates, temperature, etc.), while ecologists, in turn, are often not interested in the biochemical-physiological aspects underlying the traits they study. This review deals with several evolutionary aspects of the Drosophila alcohol dehydrogenase gene-enzyme system, and includes my own personal viewpoints. I have tried to condense and integrate the current knowledge in this field as it has developed since the comprehensive review by van Delden (1982). Details on specific issues may be gained from Sofer and Martin (1987), Sullivan, Atkinson and Starmer (1990); Chambers (1988, 1991); Geer, Miller and Heinstra (1991); and Winberg and McKinley-McKee (1992).