Molecular evolution of duplicated amylase gene regions in Drosophila melanogaster: evidence of positive selection in the coding regions and selective constraints in the cis-regulatory regions

Structural and Functional Characterization of Drosophila melanogaster α-Amylase

Insects rely on carbohydrates such as starch and glycogen as an energy supply for growth of larvae and for longevity. In this sense α-amylases have essential roles under extreme conditions, e.g., during nutritional or temperature stress, thereby contributing to survival of the insect. This makes them interesting targets for combating insect pests. Drosophila melanogaster α-amylase, DMA, which belongs to the glycoside hydrolase family 13, sub family 15, has been studied from an evolutionary, biochemical, and structural point of view. Our studies revealed that the DMA enzyme is active over a broad temperature and pH range, which is in agreement with the fluctuating environmental changes with which the insect is confronted. Crystal structures disclosed a new nearly fully solvated metal ion, only coordinated to the protein via Gln263. This residue is only conserved in the subgroup of D. melanogaster and may thus contribute to the enzyme adaptive response to large temperature variations. Studies of the effect of plant inhibitors and the pseudo-tetrasaccharide inhibitor acarbose on DMA activity, allowed us to underline the important role of the so-called flexible loop on activity/inhibition, but also to suggest that the inhibition modes of the wheat inhibitors WI-1 and WI-3 on DMA, are likely different.

A single amino-acid substitution toggles chloride dependence of the alpha-amylase paralog amyrel in Drosophila melanogaster and Drosophila virilis species

In animals, most α-amylases are chloride-dependent enzymes. A chloride ion is required for allosteric activation and is coordinated by one asparagine and two arginine side chains. Whereas the asparagine and one arginine are strictly conserved, the main chloride binding arginine is replaced by a glutamine in some rare instances, resulting in the loss of chloride binding and activation. Amyrel is a distant paralogue of α-amylase in Diptera, which was not characterized biochemically to date. Amyrel shows both substitutions depending on the species. In Drosophila melanogaster, an arginine is present in the sequence but in Drosophila virilis, a glutamine occurs at this position. We have investigated basic enzymological parameters and the dependence to chloride of Amyrel of both species, produced in yeast, and in mutants substituting arginine to glutamine or glutamine to arginine. We found that the amylolytic activity of Amyrel is about thirty times weaker than the classical Drosophila α-amylase, and that the substitution of the arginine by a glutamine in D. melanogaster suppressed the chloride-dependence but was detrimental to activity. In contrast, changing the glutamine into an arginine rendered D. virilis Amyrel chloride-dependent, and interestingly, significantly increased its catalytic efficiency. These results show that the chloride ion is not mandatory for Amyrel but stimulates the reaction rate. The possible phylogenetic origin of the arginine/glutamine substitution is also discussed.

Enzymatic characterization of recombinant α-amylase in the Drosophila melanogaster species subgroup: is there an effect of specialization on digestive enzyme?

We performed a comparative study on the enzymological features of purified recombinant α-amylase of three species belonging to the Drosophila melanogaster species subgroup: D. melanogaster, D. erecta and D. sechellia. D. erecta and D. sechellia are specialist species, with host plant Pandanus candelabrum (Pandanaceae) and Morinda citrifolia (Rubiaceae), respectively. The temperature optima were around 57-60℃ for the three species. The pH optima were 7.2 for D. melanogaster, 8.2 for D. erecta and 8.5 for D. sechellia. The kcat and Km were also estimated for each species with different substrates. The specialist species D. erecta and D. sechellia display a higher affinity for starch than D. melanogaster. α-Amylase activity is higher on starch than on glycogen in all species. α-Amylases of D. erecta and D. sechellia have a higher activity on maltooligosaccharides (G6 and G7) than on starch, contrary to D. melanogaster. Such differences in the enzymological features between the species might reflect adaptation to different ecological niches and feeding habits.

Adaptation of Musca domestica L. field population to laboratory breeding causes transcriptional alterations

BACKGROUND

The housefly, Musca domestica, has developed resistance to most insecticides applied for its control. Expression of genes coding for detoxification enzymes play a role in the response of the housefly when encountered by a xenobiotic. The highest level of constitutive gene expression of nine P450 genes was previously found in a newly-collected susceptible field population in comparison to three insecticide-resistant laboratory strains and a laboratory reference strain.

RESULTS

We compared gene expression of five P450s by qPCR as well as global gene expression by RNAseq in the newly-acquired field population (845b) in generation F1, F13 and F29 to test how gene expression changes following laboratory adaption. Four (CYP6A1, CYP6A36, CYP6D3, CYP6G4) of five investigated P450 genes adapted to breeding by decreasing expression. CYP6D1 showed higher female expression in F29 than in F1. For males, about half of the genes accessed in the global gene expression were up-regulated in F13 and F29 in comparison with the F1 population. In females, 60% of the genes were up-regulated in F13 in comparison with F1, while 33% were up-regulated in F29. Forty potential P450 genes were identified. In most cases, P450 gene expression was decreased in F13 flies in comparison with F1. Gene expression then increased from F13 to F29 in males and decreased further in females.

CONCLUSION

The global gene expression changes massively during adaptation to laboratory breeding. In general, global expression decreased as a result of laboratory adaption in males, while female expression was not unidirectional. Expression of P450 genes was in general down-regulated as a result of laboratory adaption. Expression of hexamerin, coding for a storage protein was increased, while gene expression of genes coding for amylases decreased. This suggests a major impact of the surrounding environment on gene response to xenobiotics and genetic composition of housefly strains.

Both positive and negative selection pressures contribute to the polymorphism pattern of the duplicated human CYP21A2 gene

The human steroid 21-hydroxylase gene (CYP21A2) participates in cortisol and aldosterone biosynthesis, and resides together with its paralogous (duplicated) pseudogene in a multiallelic copy number variation (CNV), called RCCX CNV. Concerted evolution caused by non-allelic gene conversion has been described in great ape CYP21 genes, and the same conversion activity is responsible for a serious genetic disorder of CYP21A2, congenital adrenal hyperplasia (CAH). In the current study, 33 CYP21A2 haplotype variants encoding 6 protein variants were determined from a European population. CYP21A2 was shown to be one of the most diverse human genes (HHe=0.949), but the diversity of intron 2 was greater still. Contrary to previous findings, the evolution of intron 2 did not follow concerted evolution, although the remaining part of the gene did. Fixed sites (different fixed alleles of sites in human CYP21 paralogues) significantly accumulated in intron 2, indicating that the excess of fixed sites was connected to the lack of effective non-allelic conversion and concerted evolution. Furthermore, positive selection was presumably focused on intron 2, and possibly associated with the previous genetic features. However, the positive selection detected by several neutrality tests was discerned along the whole gene. In addition, the clear signature of negative selection was observed in the coding sequence. The maintenance of the CYP21 enzyme function is critical, and could lead to negative selection, whereas the presumed gene regulation altering steroid hormone levels via intron 2 might help fast adaptation, which broadly characterizes the genes of human CNVs responding to the environment.

Gene duplication and ectopic gene conversion in Drosophila

The evolutionary impact of gene duplication events has been a theme of Drosophila genetics dating back to the Morgan School. While considerable attention has been placed on the genetic novelties that duplicates are capable of introducing, and the role that positive selection plays in their early stages of duplicate evolution, much less attention has been given to the potential consequences of ectopic (non-allelic) gene conversion on these evolutionary processes. In this paper we consider the historical origins of ectopic gene conversion models and present a synthesis of the current Drosophila data in light of several primary questions in the field.

Short 5'-flanking regions of the Amy gene of Drosophila kikkawai affect amylase gene expression and respond to food environments

Evolution of the duplicated genes and regulation in gene expression is of great interest, especially in terms of adaptation. Molecular population genetic and evolutionary studies on the duplicated amylase genes of Drosophila species have suggested that their 5'-flanking (cis-regulatory) regions play an important role in evolution of these genes. For better understanding of evolution of the duplicated amylase genes and gene expression, we studied functional significance of the Amy1 gene of Drosophila kikkawai using in vitro deletion mutagenesis followed by P-element-mediated germline transformation. We found that a 1.6-kb of the 5'-flanking region can produce strikingly higher level of larval amylase activity on starch food compared with that on glucose food. We found two cis-regulatory elements, which increase larval amylase activity on starch food. We also found a larval cis-regulatory element, which responds to the food difference. This food-response element is necessary for the function of the element increasing larval activity on starch food. A 5-bp deletion in a putative GRE caused high amylase activity, indicating a cis-regulatory element decreasing amylase activity. These cis-regulatory elements identified in the 5'-flanking region could be the targets of natural selection.

Genetic coadaptation of the amylase gene system in Drosophila melanogaster: evidence for the selective advantage of the lowest AMY activity and of its epistatic genetic background

In natural populations of Drosophila melanogaster, an amylase isozyme with the lowest alpha-amylase activity (AMY(1,1)) is predominant. To evaluate the selective significance of AMY(1,1) and its regulatory factor(s), we examined selection experiments in laboratory populations on two distinct food environments. After 300 generations, AMY(1,1) became predominant (89%) in a glucose (a product of AMY)-rich environment, while an isozyme with higher alpha-amylase activity, AMY(1,6), became predominant (83%) in a starch (substrate)-rich environment. We found that the identical alleles of the amylase (Amy) gene, which encodes each of AMY(1,1) and AMY(1,6), were shared between the two populations in the different food environments, employing the nucleotide sequencing of the duplicated Amy genes. Nevertheless, AMY(1,6) homozygotes selected in the starch-rich environment had a twofold higher AMY enzyme activity than those selected in the glucose-rich environment, suggesting a coadaptation of the coding region and its regulatory factor(s) on the genetic background. Such a difference in AMY enzyme activity was not detected between AMY(1,1) homozygotes, suggesting that the effect of the genetic background is epistatic. Our results indicate that natural selection is working on the Amy gene system as a whole for flies to adapt to the various food environments of local populations.

The coalescent and infinite-site model of a small multigene family

The infinite-site model of a small multigene family with two duplicated genes is studied. The expectations of the amounts of nucleotide variation within and between two genes and linkage disequilibrium are obtained, and a coalescent-based method for simulating patterns of polymorphism in a small multigene family is developed. The pattern of DNA variation is much more complicated than that in a single-copy gene, which can be simulated by the standard coalescent. Using the coalescent simulation of duplicated genes, the applicability of statistical tests of neutrality to multigene families is considered.

Applications of selective neutrality tests to molecular ecology

This paper reviews how statistical tests of neutrality have been used to address questions in molecular ecology are reviewed. The work consists of four major parts: a brief review of the current status of the neutral theory; a review of several particularly interesting examples of how statistical tests of neutrality have led to insight into ecological problems; a brief discussion of the pitfalls of assuming a strictly neutral model if it is false; and a discussion of some of the opportunities and problems that molecular ecologists face when using neutrality tests to study natural selection.

A method for estimating the mutation, gene conversion and recombination parameters in small multigene families

A simple two-locus gene conversion model is considered to investigate the amounts of DNA variation and linkage disequilibrium in small multigene families. The exact solutions for the expectations and variances of the amounts of variation within and between two loci are obtained. It is shown that gene conversion increases the amount of variation within each locus and decreases the amount of variation between two loci. The expectation and variance of the amount of linkage disequilibrium are also obtained. Gene conversion generates positive linkage disequilibrium and the degree of linkage disequilibrium decreases as the recombination rate is increased. Using the theoretical results, a method for estimating the mutation, gene conversion, and recombination parameters is developed and applied to the data of the Amy multigene family in Drosophila melanogaster. The gene conversion rate is estimated to be approximately 60-165 times higher than the mutation rate for synonymous sites.

Nucleotide variation of the duplicated amylase genes in Drosophila kikkawai

We examined levels and patterns of the nucleotide polymorphism of the Amylase genes with a head-to-head duplication in Drosophila kikkawai. The levels of variation in D. kikkawai were comparable to those in Drosophila melanogaster. Tajima's test, Fu and Li's test, HKA test, and MK test did not show significant departure from neutrality. We found an excess of replacement changes in the within-locus class, representing polymorphism in one of the duplicated genes, compared with the between-locus class, representing polymorphism shared between the duplicated genes. Most replacement changes in the within-locus class were singletons. These results suggest that most replacement changes are deleterious. A contrasting evolutionary pattern, involving concerted evolution in the coding regions but differential evolution in the 5'-flanking regions, was observed. However, unlike the duplicated Amy genes of D. melanogaster, the coding regions of the duplicated genes in D. kikkawai tended to diverge. Using Ohta's model of the small multigene family, we found that recombination (interchromosomal equal crossing-over) rate was one order higher than gene conversion (unequal crossing-over) rate, resulting in a considerable but incomplete homogenization of the duplicated coding regions. Linkage disequilibria were found in the intron as well as within and around the regulatory cis-element sequences of one of the duplicated genes (Amy1). The possible causes of these linkage disequilibria were discussed.