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

De Novo Transcriptome Analysis Reveals Potential Thermal Adaptation Mechanisms in the Cicada Hyalessa fuscata

Simple Summary

In metropolitan Seoul and its vicinity, cicadas of the species Hyalessa fuscata living in warmer areas could tolerate the heat better than those living in cooler areas, but genetic mechanisms involved in better heat tolerance remained unclear. In this study, we examined differences in gene expression of cicadas living in a warm urban area, a cool urban area and a suburban area in three experimental treatments: no heating, 10 min heating and heating until the cicadas lost their mobility. Cicadas from the warm urban area changed their gene expressions the most. Activated genes were mostly related to heat shock, energy metabolism, and detoxification. These results suggested that under heat stress, cicadas inhabiting warm areas could differentially express genes to increase their thermal tolerance.

Abstract

In metropolitan Seoul, populations of the cicada Hyalessa fuscata in hotter urban heat islands (“high UHIs”) exhibit higher thermal tolerance than those in cooler UHIs (“low UHIs”). We hypothesized that heat stress may activate the expression of genes that facilitate greater thermal tolerance in high-UHI cicadas than in those from cooler areas. Differences in the transcriptomes of adult female cicadas from high-UHI, low-UHI, and suburban areas were analyzed at the unheated level, after acute heat stress, and after heat torpor. No noticeable differences in unheated gene expression patterns were observed. After 10 min of acute heat stress, however, low-UHI and suburban cicadas expressed more heat shock protein genes than high-UHI counterparts. More specifically, remarkable changes in the gene expression of cicadas across areas were observed after heat torpor stimulus, as represented by a large number of up- and downregulated genes in the heat torpor groups compared with the 10 min acute heat stress and control groups. High-UHI cicadas expressed the most differentially expressed genes, followed by the low-UHI and suburban cicadas. There was a notable increase in the expression of heat shock, metabolism, and detoxification genes; meanwhile, immune-related, signal transduction, and protein turnover genes were downregulated in high-UHI cicadas versus the other cicada groups. These results suggested that under heat stress, cicadas inhabiting high-UHIs could rapidly express genes related to heat shock, energy metabolism, and detoxification to protect cells from stress-induced damage and to increase their thermal tolerance toward heat stress. The downregulation of apoptosis mechanisms in high-UHI cicadas suggested that there was less cellular damage, which likely contributed to their high tolerance of heat stress.

Amyrel, a novel glucose-forming α-amylase from Drosophila with 4-α-glucanotransferase activity by disproportionation and hydrolysis of maltooligosaccharides

The α-amylase paralogue Amyrel present in true flies (Diptera Muscomorpha) has been classified as a glycoside hydrolase in CAZy family GH13 on the basis of its primary structure. Here we report that, in fact, Amyrel is currently unique amongst Animals as it possesses both the hydrolytic α-amylase activity (EC 3.2.1.1) and a 4-α-glucanotransferase (EC 2.4.1.25) transglycosylation activity. Amyrel reacts specifically on α-(1-4) glycosidic bonds of starch and related polymers but produces a complex mixture of maltooligosaccharides, in sharp contrast with canonical animal α-amylases. With model maltooligosaccharides G2 (maltose) to G7, the Amyrel reaction starts by a disproportionation leading to Gn-1 and Gn + 1 products, which become themselves substrates for new disproportionation cycles. As a result, all detectable odd- and even-numbered maltooligosaccharides at least up to G12 were observed. However, hydrolysis of these products proceeds simultaneously, as shown by p-nitrophenyl-tagged oligosaccharides and microcalorimetry, and upon prolonged reaction, glucose is the major end product followed by maltose. The main structural determinant of these atypical activities was found to be a Gly-His-Gly-Ala deletion in the so-called flexible loop bordering the active site. Indeed, engineering this deletion in pig pancreatic and D. melanogaster α-amylases results in reaction patterns similar to those of Amyrel. It is proposed that this deletion provides more freedom to the substrate for subsites occupancy and allows a less constrained action pattern resulting in versatile activities at the active site.

New insights into the origin and evolution of α-amylase genes in green plants

Gene duplication is a source of genetic materials and evolutionary changes, and has been associated with gene family expansion. Functional divergence of duplicated genes is strongly directed by natural selections such as organism diversification and novel feature acquisition. We show that, plant α-amylase gene family (AMY) is comprised of six subfamilies (AMY1-AMY6) that fell into two ancient phylogenetic lineages (AMY3 and AMY4). Both AMY1 and AMY2 are grass-specific and share a single-copy ancestor, which is derived from grass AMY3 genes that have undergone massive tandem and whole-genome duplications during evolution. Ancestral features of AMY4 and AMY5/AMY6 genes have been retained among four green algal sequences (Chrein_08.g362450, Vocart_0021s0194, Dusali_0430s00012 and Monegl_16464), suggesting a gene duplication event following Chlorophyceae diversification. The observed horizontal gene transfers between plant and bacterial AMYs, and chromosomal locations of AMY3 and AMY4 genes in the most ancestral green body (C. reinhardtii), provide evidences for the monophyletic origin of plant AMYs. Despite subfamily-specific sequence divergence driven by natural selections, the active site and SBS1 are well-conserved across different AMY isoforms. The differentiated electrostatic potentials and hydrogen bands-forming residue polymorphisms, further imply variable digestive abilities for a broad substrates in particular tissues or subcellular localizations.

The Amylases of Insects

Alpha-amylases are major digestive enzymes that act in the first step of maltopolysaccharide digestion. In insects, these enzymes have long been studied for applied as well as purely scientific purposes. In many species, amylases are produced by multiple gene copies. Rare species are devoid of Amy gene. They are predominantly secreted in the midgut but salivary expression is also frequent, with extraoral activity. Enzymological parameters are quite variable among insects, with visible trends according to phylogeny: Coleopteran amylases have acidic optimum activity, whereas dipteran amylases have neutral preference and lepidopteran ones have clear alkaline preference. The enzyme structure shows interesting variations shaped by evolutionary convergences, such as the recurrent loss of a loop involved in substrate handling. Many works have focused on the action of plant amylase inhibitors on pest insect amylases, in the frame of crop protection by transgenesis. It appears that sensitivity or resistance to inhibitors is finely tuned and very specific and that amylases and their inhibitors have coevolved. The multicopy feature of insect amylases appears to allow tissue-specific or stage-specific regulation, but also to broaden enzymological abilities, such as pH range, and to overcome plant inhibitory defenses.

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.

LaaA, a novel high-active alkalophilic alpha-amylase from deep-sea bacterium Luteimonas abyssi XH031(T)

Alpha-amylase is a kind of broadly used industrial enzymes, most of which have been exploited from terrestrial organism. Comparatively, alpha-amylase from marine environment was largely undeveloped. In this study, a novel alkalophilic alpha-amylase with high activity, Luteimonas abyssi alpha-amylase (LaaA), was cloned from deep-sea bacterium L. abyssi XH031(T) and expressed in Escherichia coli BL21. The gene has a length of 1428bp and encodes 475 amino acids with a 35-residue signal peptide. The specific activity of LaaA reached 8881U/mg at the optimum pH 9.0, which is obvious higher than other reported alpha-amylase. This enzyme can remain active at pH levels ranging from 6.0 to 11.0 and temperatures below 45°C, retaining high activity even at low temperatures (almost 38% residual activity at 10°C). In addition, 1mM Na(+), K(+), and Mn(2+) enhanced the activity of LaaA. To investigate the function of potential active sites, R227G, D229K, E256Q/H, H327V and D328V mutants were generated, and the results suggested that Arg227, Asp229, Glu256 and Asp328 were total conserved and essential for the activity of alpha-amylase LaaA. This study shows that the alpha-amylase LaaA is an alkali-tolerant and high-active amylase with strong potential for use in detergent industry.

Evolution of behavioural and cellular defences against parasitoid wasps in the Drosophila melanogaster subgroup

It may be intuitive to predict that host immune systems will evolve to counter a broad range of potential challenges through simultaneous investment in multiple defences. However, this would require diversion of resources from other traits, such as growth, survival and fecundity. Therefore, ecological immunology theory predicts that hosts will specialize in only a subset of possible defences. We tested this hypothesis through a comparative study of a cellular immune response and a putative behavioural defence used by eight fruit fly species against two parasitoid wasp species (one generalist and one specialist). Fly larvae can survive infection by melanotically encapsulating wasp eggs, and female flies can potentially reduce infection rates in their offspring by laying fewer eggs when wasps are present. The strengths of both defences varied significantly but were not negatively correlated across our chosen host species; thus, we found no evidence for a trade-off between behavioural and cellular immunity. Instead, cellular defences were significantly weaker against the generalist wasp, whereas behavioural defences were similar in strength against both wasps and positively correlated between wasps. We investigated the adaptive significance of wasp-induced oviposition reduction behaviour by testing whether wasp-exposed parents produce offspring with stronger cellular defences, but we found no support for this hypothesis. We further investigated the sensory basis of this behaviour by testing mutants deficient in either vision or olfaction, both of which failed to reduce their oviposition rates in the presence of wasps, suggesting that both senses are necessary for detecting and responding to wasps.

α-Amylase: an enzyme specificity found in various families of glycoside hydrolases

α-Amylase (EC 3.2.1.1) represents the best known amylolytic enzyme. It catalyzes the hydrolysis of α-1,4-glucosidic bonds in starch and related α-glucans. In general, the α-amylase is an enzyme with a broad substrate preference and product specificity. In the sequence-based classification system of all carbohydrate-active enzymes, it is one of the most frequently occurring glycoside hydrolases (GH). α-Amylase is the main representative of family GH13, but it is probably also present in the families GH57 and GH119, and possibly even in GH126. Family GH13, known generally as the main α-amylase family, forms clan GH-H together with families GH70 and GH77 that, however, contain no α-amylase. Within the family GH13, the α-amylase specificity is currently present in several subfamilies, such as GH13_1, 5, 6, 7, 15, 24, 27, 28, 36, 37, and, possibly in a few more that are not yet defined. The α-amylases classified in family GH13 employ a reaction mechanism giving retention of configuration, share 4-7 conserved sequence regions (CSRs) and catalytic machinery, and adopt the (β/α)8-barrel catalytic domain. Although the family GH57 α-amylases also employ the retaining reaction mechanism, they possess their own five CSRs and catalytic machinery, and adopt a (β/α)7-barrel fold. These family GH57 attributes are likely to be characteristic of α-amylases from the family GH119, too. With regard to family GH126, confirmation of the unambiguous presence of the α-amylase specificity may need more biochemical investigation because of an obvious, but unexpected, homology with inverting β-glucan-active hydrolases.

The female-specific doublesex isoform regulates pleiotropic transcription factors to pattern genital development in Drosophila

Regulatory networks driving morphogenesis of animal genitalia must integrate sexual identity and positional information. Although the genetic hierarchy that controls somatic sexual identity in the fly Drosophila melanogaster is well understood, there are very few cases in which the mechanism by which it controls tissue-specific gene activity is known. In flies, the sex-determination hierarchy terminates in the doublesex (dsx) gene, which produces sex-specific transcription factors via alternative splicing of its transcripts. To identify sex-specifically expressed genes downstream of dsx that drive the sexually dimorphic development of the genitalia, we performed genome-wide transcriptional profiling of dissected genital imaginal discs of each sex at three time points during early morphogenesis. Using a stringent statistical threshold, we identified 23 genes that have sex-differential transcript levels at all three time points, of which 13 encode transcription factors, a significant enrichment. We focus here on three sex-specifically expressed transcription factors encoded by lozenge (lz), Drop (Dr) and AP-2. We show that, in female genital discs, Dsx activates lz and represses Dr and AP-2. We further show that the regulation of Dr by Dsx mediates the previously identified expression of the fibroblast growth factor Branchless in male genital discs. The phenotypes we observe upon loss of lz or Dr function in genital discs explain the presence or absence of particular structures in dsx mutant flies and thereby clarify previously puzzling observations. Our time course of expression data also lays the foundation for elucidating the regulatory networks downstream of the sex-specifically deployed transcription factors.

The Drosophila NR4A nuclear receptor DHR38 regulates carbohydrate metabolism and glycogen storage

Animals balance nutrient storage and mobilization to maintain metabolic homeostasis, a process that is disrupted in metabolic diseases like obesity and diabetes. Here, we show that DHR38, the single fly ortholog of the mammalian nuclear receptor 4A family of nuclear receptors, regulates glycogen storage during the larval stages of Drosophila melanogaster. DHR38 is expressed and active in the gut and body wall of larvae, and its expression levels change in response to nutritional status. DHR38 null mutants have normal levels of glucose, trehalose (the major circulating form of sugar), and triacylglycerol but display reduced levels of glycogen in the body wall muscles, which constitute the primary storage site for carbohydrates. Microarray analysis reveals that many metabolic genes are mis-regulated in DHR38 mutants. These include phosphoglucomutase, which is required for glycogen synthesis, and the two genes that encode the digestive enzyme amylase, accounting for the reduced amylase enzyme activity seen in DHR38 mutant larvae. These studies demonstrate that a critical role of nuclear receptor 4A receptors in carbohydrate metabolism has been conserved through evolution and that nutritional regulation of DHR38 expression maintains the proper uptake and storage of glycogen during the growing larval stage of development.

Digestive alpha-amylases of the flour moth Ephestia kuehniella--adaptation to alkaline environment and plant inhibitors

The digestive tract of lepidopteran insects is extremely alkaline. In the present work, molecular adaptation of amylolytic enzymes to this environment was investigated in the flour moth Ephestia kuehniella, an important stored-product pest. Three digestive alpha-amylases [Ephestia kuehniella alpha-amylase isoenzymes 1-3 (EkAmy1-3)] with an alkaline pH optimum were purified from larvae and biochemically characterized. These isoenzymes differ significantly in their sensitivity to alpha-amylase inhibitors of plant origin that are directed against herbivores as antifeedants. Such functional variability renders the amylolytic system less vulnerable to suppression by plant defensive molecules. Moreover, we found that expression of alpha-amylases is upregulated in larvae feeding on a diet enriched with an alpha-amylase inhibitor. The alpha-amylases are secreted into the larval midgut by an exocytotic mechanism, as revealed by immunogold microscopy. The cDNA sequence of EkAmy3 was determined, and a homology model of EkAmy3 was built in order to analyze the structural features responsible for adaptation to alkaline pH. First, the overall fold was found to be stabilized by remodeling of ion pairs. Second, molecular simulations supported by activity measurements showed that EkAmy3 does not bind a Cl(-), owing to an Arg-to-Gln mutation in a conserved binding site. The Cl(-)-binding residues are in contact with the catalytic residues, and this change might help to fine-tune the catalytic pK(a) values to an alkaline pH optimum. We conclude that lepidopteran alpha-amylases are evolutionarily adapted in terms of structure and expression dynamics for effective functioning in the digestive system.

Transcriptional signatures in response to wheat germ agglutinin and starvation in Drosophila melanogaster larval midgut

One function of plant lectins such as wheat germ agglutinin is to serve as defences against herbivorous insects. The midgut is one critical site affected by dietary lectins. We observed marked cellular, structural and gene expression changes in the midguts of Drosophila melanogaster third instar larvae that were fed wheat germ agglutinin. Some of these changes were similar to those observed in the midguts of starved D. melanogaster. Dietary wheat germ agglutinin caused shortening, branching, swelling, distortion and in some cases disintegration of the midgut microvilli. Starvation was accompanied primarily by shortening of the microvilli. Microarray analyses revealed that dietary wheat germ agglutinin evoked differential expression of 61 transcripts; seven of these were also differentially expressed in starved D. melanogaster. The differentially transcribed gene clusters in wheat germ agglutinin-fed larvae were associated with (1) cytoskeleton organization; (2) digestive enzymes; (3) detoxification reactions; and (4) energy metabolism. Four possible transcription factor binding motifs were associated with the differentially expressed genes. One of these exhibited substantial similarity to MyoD, a transcription factor binding motif associated with cellular structures in mammals. These results are consistent with the hypothesis that wheat germ agglutinin caused a starvation-like effect and structural changes of midgut cells of D. melanogaster third-instar larvae.

Grafting the molecular phylogenetic tree with morphological branches to reconstruct the evolutionary history of the genus Zaprionus (Diptera: Drosophilidae)

A molecular phylogeny for the drosophilid genus Zaprionus was inferred using a mitochondrial (CO-II) and a nuclear (Amyrel) gene using 22 available species. The combined molecular tree does not support the current classification, dubbed phylogenetic, based entirely upon a morphocline of forefemoral ornamentation. For species for which DNA was not available, phylogenetic positioning was only assigned using morphological characters. In order to avoid conflict between DNA and morphology in the combined analyses (supermatrix method), we developed a new method in which few morphological characters were sampled according to an a priori homoplasy assessment on the consensus molecular tree. At each internal node of the tree, a number of synapomorphies was determined, and species with no molecular sequences were grafted thereon. Analogously to tree vocabulary, we called our method 'morphological grafting'. New species groups and complexes were then defined in the light of our findings. Further, divergence times were estimated under a relaxed molecular clock, and historical biogeography was reconstructed under a maximum likelihood model. Zaprionus appears to be of recent origin in the Oriental region during the Late Miocene ( approximately 10 MYA), and colonization of Africa started shortly after ( approximately 7 MYA) via the maritime route of the Indian Ocean Islands. Most of the morphological and ecological diversification took place, later, in Western Africa during the Quaternary cyclic climatic changes. Furthermore, some species became recent invaders, with one, Zaprionus indianus, has successfully invaded South and North America during the last decade.

Origin and evolution of the Amyrel gene in the alpha-amylase multigene family of Diptera

Alpha-amylase genes often form multigene families in living organisms. In Diptera, a remote paralog, Amyrel, had been discovered in Drosophila, where this gene is currently used as a population and phylogenetic marker. The putative encoded protein has about 40% divergence with the classical amylases. We have searched the presence of the paralog in other families of Diptera to track its origin and understand its evolution. Amyrel was detected in a number of families of Muscomorpha (Brachycera-Cyclorrapha), suggesting an origin much older than previously thought. It has not been found elsewhere to date, and it is absent from the Anopheles gambiae genome. The intron-exon structures of the genes found so far suggest that the ancestral gene (before the duplication which gave rise to Amyrel) had two introns, and that subsequent, repeated and independent loss of one or both introns occurred in some Muscomorpha families. It seems that the Amyrel protein has experienced specific amino acid substitutions in regions generally well conserved in amylases, raising the possibility of peculiar, functional adaptations of this protein.

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.

Divergence between Drosophila santomea and allopatric or sympatric populations of D. yakuba using paralogous amylase genes and migration scenarios along the Cameroon volcanic line

We have used two paralogous genes (Amyrel and Amy) of the amylase multigene family to reconstruct the phylogeny of the nine Drosophila melanogaster subgroup sister species, including D. santomea, the newly discovered endemic from São Tomé island. The evolutionary divergence of these genes is of special interest as it is suspected to result from physiological evolution via gene duplication. This paper describes the relationship between the geographical origin of the various strains and the patterns of mating and phylogeny, focusing on the evolution of D. santomea and its relationship to other species and their niches. The Amyrel and Amy data indicate that, contrary to expectations, the sympatric insular D. yakuba population is less closely related to D. santomea than allopatric mainland ones, suggesting that the extant insular D. yakuba population on São Tomé results from a recent secondary colonization. Data for sympatric and allopatric D. yakuba suggest that D. santomea arose from a mainland D. yakuba parental stock when montane habitats of the Cameroon volcanic line extended to lower altitudes during colder and less humid periods. Despite their different modes of evolution and different functions, the Amyrel and Amy genes provide remarkably consistent topologies and hence reflect the same history, that of the species.

The population genetics of the origin and divergence of the Drosophila simulans complex species

The origins and divergence of Drosophila simulans and close relatives D. mauritiana and D. sechellia were examined using the patterns of DNA sequence variation found within and between species at 14 different genes. D. sechellia consistently revealed low levels of polymorphism, and genes from D. sechellia have accumulated mutations at a rate that is approximately 50% higher than the same genes from D. simulans. At synonymous sites, D. sechellia has experienced a significant excess of unpreferred codon substitutions. Together these observations suggest that D. sechellia has had a reduced effective population size for some time, and that it is accumulating slightly deleterious mutations as a result. D. simulans and D. mauritiana are both highly polymorphic and the two species share many polymorphisms, probably since the time of common ancestry. A simple isolation speciation model, with zero gene flow following incipient species separation, was fitted to both the simulans/mauritiana divergence and the simulans/sechellia divergence. In both cases the model fit the data quite well, and the analyses revealed little evidence of gene flow between the species. The exception is one gene copy at one locus in D. sechellia, which closely resembled other D. simulans sequences. The overall picture is of two allopatric speciation events that occurred quite near one another in time.

Evidence for a high ancestral GC content in Drosophila

Study of the nucleotide composition in Drosophila, focusing on the saltans and willistoni groups, has revealed unanticipated differences in nucleotide composition among lineages. Compositional differences are associated with an accelerated rate of nucleotide substitution in functionally less constrained regions. These observations have been set forth against the extended opinion that the pattern of point mutation has remained constant during the evolution of the genus. A crucial assumption has been that the most recent common ancestor of the subgenus Sophophora had an elevated GC content. Until now, this assumption has been supported by indirect arguments, consisting of extrapolations from closely related outgroups and limited by the robustness of mathematical descriptions concerning the extensive nucleotide composition differences among sequences. The present study seeks to test the assumption of a high ancestral GC content using realistic representations of the nucleotide substitution process to account for potential biases induced by the heterogeneous GC content of the taxa. The analysis of eight nuclear genes unambiguously corroborates that the common ancestor of Sophophora had an elevated GC content.

Evolution of nucleotide substitutions and gene regulation in the amylase multigenes in Drosophila kikkawai and its sibling species

In order to determine evolutionary changes in gene regulation and the nucleotide substitution pattern in a multigene family, the amylase multigenes were characterized in Drosophila kikkawai and its sibling species. The nucleotide substitution pattern was investigated. Drosophila kikkawai has four amylase genes. The Amy1 and Amy2 genes are a head-to-head duplication in the middle of the B arm of the second chromosome, while the Amy3 and Amy4 genes are a tail-to-tail duplication near the centromere of the same chromosome. In the sibling species of D. kikkawai (Drosophila bocki, Drosophila leontia, and Drosophila lini), sequencing of the Amy1, Amy2, Amy3, and Amy4 genes revealed that the Amy1 and Amy2 gene group diverged from Amy3 and Amy4 after duplication. In the Amy1 and Amy2 genes, the divergent evolution occurred in the flanking regions; in contrast, the coding regions have evolved in concerted fashion. The electrophoretic pattern of AMY isozymes was also examined. In D. kikkawai and its siblings, two or three electrophoretically different isozymes are encoded by the Amy1 and Amy2 genes (S isozyme) and by the Amy3 and Amy4 genes (F (M) isozymes). The S and F (M) isozymes show different patterns of band intensity when larvae and flies were fed in different media. Amy1 and Amy2, which encode the S isozyme, are more strikingly regulated than Amy3 and Amy4, which encode the F (M) isozyme. The GC content and codon usage bias were higher for the Amy1 and Amy2 genes than for the Amy3 and Amy4 genes. Although the ratio of synonymous and replacement substitutions within the Amy1 and Amy2 gene group was not significantly different from that within the Amy3 and Amy4 gene group, the synonymous substitution rate in the lineage of Amy1 and Amy2 was lower than that of Amy3 and Amy4. In conclusion, after the first duplication but before speciation of four species, the synonymous substitution rate between the two lineages and the electrophoretic pattern of the isozymes encoded by them changed, although we do not know whether there was any evolutionary relationship between the two.

The Amylase gene cluster on the evolving sex chromosomes of Drosophila miranda

On the basis of chromosomal homology, the Amylase gene cluster in Drosophila miranda must be located on the secondary sex chromosome pair, neo-X (X2) and neo-Y, but is autosomally inherited in all other Drosophila species. Genetic evidence indicates no active amylase on the neo-Y chromosome and the X2-chromosomal locus already shows dosage compensation. Several lines of evidence strongly suggest that the Amy gene cluster has been lost already from the evolving neo-Y chromosome. This finding shows that a relatively new neo-Y chromosome can start to lose genes and hence gradually lose homology with the neo-X. The X2-chromosomal Amy1 is intact and Amy2 contains a complete coding sequence, but has a deletion in the 3'-flanking region. Amy3 is structurally eroded and hampered by missing regulatory motifs. Functional analysis of the X2-chromosomal Amy1 and Amy2 regions from D. miranda in transgenic D. melanogaster flies reveals ectopic AMY1 expression. AMY1 shows the same electrophoretic mobility as the single amylase band in D. miranda, while ectopic AMY2 expression is characterized by a different mobility. Therefore, only the Amy1 gene of the resident Amy cluster remains functional and hence Amy1 is the dosage compensated gene.