Concerted evolution of duplicated protein-coding genes in Drosophila

Genetic variation in recalcitrant repetitive regions of the Drosophila melanogaster genome

Many essential functions of organisms are encoded in highly repetitive genomic regions, including histones involved in DNA packaging, centromeres that are core components of chromosome segregation, ribosomal RNA comprising the protein translation machinery, telomeres that ensure chromosome integrity, piRNA clusters encoding host defenses against selfish elements, and virtually the entire Y chromosome. These regions, formed by highly similar tandem arrays, pose significant challenges for experimental and informatic study, impeding sequence-level descriptions essential for understanding genetic variation. Here, we report the assembly and variation analysis of such repetitive regions in Drosophila melanogaster, offering significant improvements to the existing community reference assembly. Our work successfully recovers previously elusive segments, including complete reconstructions of the histone locus and the pericentric heterochromatin of the X chromosome, spanning the Stellate locus to the distal flank of the rDNA cluster. To infer structural changes in these regions where alignments are often not practicable, we introduce landmark anchors based on unique variants that are putatively orthologous. These regions display considerable structural variation between different D. melanogaster strains, exhibiting differences in copy number and organization of homologous repeat units between haplotypes. In the histone cluster, although we observe minimal genetic exchange indicative of crossing over, the variation patterns suggest mechanisms such as unequal sister chromatid exchange. We also examine the prevalence and scale of concerted evolution in the histone and Stellate clusters and discuss the mechanisms underlying these observed patterns.

Insect α-Amylases and Their Application in Pest Management

Amylase is an indispensable hydrolase in insect growth and development. Its varied enzymatic parameters cause insects to have strong stress resistance. Amylase gene replication is a very common phenomenon in insects, and different copies of amylase genes enable changes in its location and function. In addition, the classification, structure, and interaction between insect amylase inhibitors and amylases have also invoked the attention of researchers. Some plant-derived amylase inhibitors have inhibitory activities against insect amylases and even mammalian amylases. In recent years, an increasing number of studies have clarified the effects of pesticides on the amylase activity of target and non-target pests, which provides a theoretical basis for exploring safe and efficient pesticides, while the exact lethal mechanisms and safety in field applications remain unclear. Here, we summarize the most recent advances in insect amylase studies, including its sequence and characteristics and the regulation of amylase inhibitors (α-AIs). Importantly, the application of amylases as the nanocide trigger, RNAi, or other kinds of pesticide targets will be discussed. A comprehensive foundation will be provided for applying insect amylases to the development of new-generation insect management tools and improving the specificity, stability, and safety of pesticides.

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.

The evolution of genomic GC content undergoes a rapid reversal within the genus Plasmodium

The genome of the malarial parasite Plasmodium falciparum is extremely AT rich. This bias toward a low GC content is a characteristic of several, but not all, species within the genus Plasmodium. We compared 4283 orthologous pairs of protein-coding sequences between Plasmodium falciparum and the less AT-biased Plasmodium vivax. Our results indicate that the common ancestor of these two species was also extremely AT rich. This means that, although there was a strong bias toward A+T during the early evolution of the ancestral Plasmodium lineage, there was a subsequent reversal of this trend during the more recent evolution of some species, such as P. vivax. Moreover, we show that not only is the P. vivax genome losing its AT richness, it is actually gaining a very significant degree of GC richness. This example illustrates the potential volatility of nucleotide content during the course of molecular evolution. Such reversible fluxes in nucleotide content within lineages could have important implications for phylogenetic reconstruction based on molecular sequence data.

The natural history of class I primate alcohol dehydrogenases includes gene duplication, gene loss, and gene conversion

BACKGROUND

Gene duplication is a source of molecular innovation throughout evolution. However, even with massive amounts of genome sequence data, correlating gene duplication with speciation and other events in natural history can be difficult. This is especially true in its most interesting cases, where rapid and multiple duplications are likely to reflect adaptation to rapidly changing environments and life styles. This may be so for Class I of alcohol dehydrogenases (ADH1s), where multiple duplications occurred in primate lineages in Old and New World monkeys (OWMs and NWMs) and hominoids.

METHODOLOGY/PRINCIPAL FINDINGS

To build a preferred model for the natural history of ADH1s, we determined the sequences of nine new ADH1 genes, finding for the first time multiple paralogs in various prosimians (lemurs, strepsirhines). Database mining then identified novel ADH1 paralogs in both macaque (an OWM) and marmoset (a NWM). These were used with the previously identified human paralogs to resolve controversies relating to dates of duplication and gene conversion in the ADH1 family. Central to these controversies are differences in the topologies of trees generated from exonic (coding) sequences and intronic sequences.

CONCLUSIONS/SIGNIFICANCE

We provide evidence that gene conversions are the primary source of difference, using molecular clock dating of duplications and analyses of microinsertions and deletions (micro-indels). The tree topology inferred from intron sequences appear to more correctly represent the natural history of ADH1s, with the ADH1 paralogs in platyrrhines (NWMs) and catarrhines (OWMs and hominoids) having arisen by duplications shortly predating the divergence of OWMs and NWMs. We also conclude that paralogs in lemurs arose independently. Finally, we identify errors in database interpretation as the source of controversies concerning gene conversion. These analyses provide a model for the natural history of ADH1s that posits four ADH1 paralogs in the ancestor of Catarrhine and Platyrrhine primates, followed by the loss of an ADH1 paralog in the human lineage.

Annelid Distal-less/Dlx duplications reveal varied post-duplication fates

Background

Dlx (Distal-less) genes have various developmental roles and are widespread throughout the animal kingdom, usually occurring as single copy genes in non-chordates and as multiple copies in most chordate genomes. While the genomic arrangement and function of these genes is well known in vertebrates and arthropods, information about Dlx genes in other organisms is scarce. We investigate the presence of Dlx genes in several annelid species and examine Dlx gene expression in the polychaete Pomatoceros lamarckii.

Results

Two Dlx genes are present in P. lamarckii, Capitella teleta and Helobdella robusta. The C. teleta Dlx genes are closely linked in an inverted tail-to-tail orientation, reminiscent of the arrangement of vertebrate Dlx pairs, and gene conversion appears to have had a role in their evolution. The H. robusta Dlx genes, however, are not on the same genomic scaffold and display divergent sequences, while, if the P. lamarckii genes are linked in a tail-to-tail orientation they are a minimum of 41 kilobases apart and show no sign of gene conversion. No expression in P. lamarckii appendage development has been observed, which conflicts with the supposed conserved role of these genes in animal appendage development. These Dlx duplications do not appear to be annelid-wide, as the polychaete Platynereis dumerilii likely possesses only one Dlx gene.

Conclusions

On the basis of the currently accepted annelid phylogeny, we hypothesise that one Dlx duplication occurred in the annelid lineage after the divergence of P. dumerilii from the other lineages and these duplicates then had varied evolutionary fates in different species. We also propose that the ancestral role of Dlx genes is not related to appendage development.

Drosophila genes that affect meiosis duration are among the meiosis related genes that are more often found duplicated

Using a phylogenetic approach, the examination of 33 meiosis/meiosis-related genes in 12 Drosophila species, revealed nine independent gene duplications, involving the genes cav, mre11, meiS332, polo and mtrm. Evidence is provided that at least eight out of the nine gene duplicates are functional. Therefore, the rate at which Drosophila meiosis/meiosis-related genes are duplicated and retained is estimated to be 0.0012 per gene per million years, a value that is similar to the average for all Drosophila genes. It should be noted that by using a phylogenetic approach the confounding effect of concerted evolution, that is known to lead to overestimation of the duplication and retention rate, is avoided. This is an important issue, since even in our moderate size sample, evidence for long-term concerted evolution (lasting for more than 30 million years) was found for the meiS332 gene pair in species of the Drosophila subgenus. Most striking, in contrast to theoretical expectations, is the finding that genes that encode proteins that must follow a close stoichiometric balance, such as polo, mtrm and meiS332 have been found duplicated. The duplicated genes may be examples of gene neofunctionalization. It is speculated that meiosis duration may be a trait that is under selection in Drosophila and that it has different optimal values in different species.

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.

Genomic structure and evolution of the Pi2/9 locus in wild rice species

Rice blast, caused by the fungal pathogen Magnaporthe oryzae, is a devastating disease of rice worldwide. Among the 85 mapped resistance (R) genes against blast, 13 have been cloned and characterized. However, how these genes originated and how they evolved in the Oryza genus remains unclear. We previously cloned the rice blast R-genes Pi2, Pi9, and Piz-t, and analyzed their genomic structure and evolution in cultivated rice. In this study, we determined the genomic sequences of the Pi2/9 locus in four wild Oryza species representing three genomes (AA, BB and CC). The number of Pi2/9 family members in the four wild species ranges from two copies to 12 copies. Although these genes are conserved in structure and categorized into the same subfamily, sequence duplications and subsequent inversions or uneven crossing overs were observed, suggesting that the locus in different wild species has undergone dynamic changes. Positive selection was found in the leucine-rich repeat region of most members, especially in the largest clade where Pi9 is included. We also provide evidence that the Pi9 gene is more related to its homologues in the recurrent line and other rice cultivars than to those in its alleged donor species O. minuta, indicating a possible origin of the Pi9 gene from O. sativa. Comparative sequence analysis between the four wild Oryza species and the previously established reference sequences in cultivated rice species at the Pi2/9 locus has provided extensive and unique information on the genomic structure and evolution of a complex R-gene cluster in the Oryza genus.

Nonallelic gene conversion in the genus Drosophila

Nonallelic gene conversion has been proposed as a major force in homogenizing the sequences of paralogous genes. In this work, we investigate the extent and characteristics of gene conversion among gene families in nine species of the genus Drosophila. We carried out a genome-wide study of 2855 gene families (including 17,742 genes) and determined that conversion events involved 2628 genes. The proportion of converted genes ranged across species from 1 to 9% when paralogs of all ages were included. Although higher levels of gene conversion were found among young gene duplicates, at most 1-2% of the coding sequences of these duplicates were affected by conversion. Using a second approach relying on gene family size changes and gene-tree/species-tree reconciliation methods, we estimate that only 1-15% of gene trees are misled by gene conversion, depending on the lineage considered. Several features of paralogous genes correlate with gene conversion, such as intra-/interchromosomal location, level of nucleotide divergence, and GC content, although we found no definitive evidence for biased substitution patterns. After considering species-specific differences in the age and distance between paralogs, we found a highly significant difference in the amount of gene conversion among species. In particular, members of the melanogaster group showed the lowest proportion of converted genes. Our data therefore suggest underlying differences in the mechanistic basis of gene conversion among species.

Nucleotide substitution bias within the genus Drosophila affects the pattern of proteome evolution

The availability of complete genome sequences for 12 Drosophila species provides an unprecedented resource for large-scale studies of genome evolution. In this study, we looked for correlated shifts in the patterns of genome and proteome evolution within the genus Drosophila. Specifically, we asked if the nucleotide composition of the Drosophila willistoni genome--which is significantly less GC rich than the other 11 sequenced Drosophila genomes--is reflected in an altered pattern of amino acid substitutions in the encoded proteins. Our results show that this is indeed the case: There are large and highly significant asymmetries in the patterns of amino acid substitution between D. willistoni and Drosophila melanogaster, and they are in the direction predicted by the nucleotide biases. The implication of this result, combined with previous studies on long-term proteome evolution, is that substitutional biases at the DNA level can be a major factor in determining both the long-term and the short-term directions of proteome evolution.

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.

Identification of a biomarker for sleep drive in flies and humans

It is a common experience to sacrifice sleep to meet the demands of our 24-h society. Current estimates reveal that as a society, we sleep on average 2 h less than we did 40 years ago. This level of sleep restriction results in negative health outcomes and is sufficient to produce cognitive deficits and reduced attention and is associated with increased risk for traffic and occupational accidents. Unfortunately, there is no simple quantifiable marker that can detect an individual who is excessively sleepy before adverse outcomes become evident. To address this issue, we have developed a simple and effective strategy for identifying biomarkers of sleepiness by using genetic and pharmacological tools that dissociate sleep drive from wake time in the model organism Drosophila melanogaster. These studies have identified a biomarker, Amylase, that is highly correlated with sleep drive. More importantly, both salivary Amylase activity and mRNA levels are also responsive to extended waking in humans. These data indicate that the fly is relevant for human sleep research and represents a first step in developing an effective method for detecting sleepiness in vulnerable populations.

The evolution of hexapod engrailed-family genes: evidence for conservation and concerted evolution

Phylogenetic analyses imply that multiple engrailed-family gene duplications occurred during hexapod evolution, a view supported by previous reports of only a single engrailed-family gene in members of the grasshopper genus Schistocerca and in the beetle Tribolium castaneum. Here, we report the cloning of a second engrailed-family gene from Schistocerca gregaria and present evidence for two engrailed-family genes from four additional hexapod species. We also report the existence of a second engrailed-family gene in the Tribolium genome. We suggest that the engrailed and invected genes of Drosophila melanogaster have existed as a conserved gene cassette throughout holometabolous insect evolution. In total 11 phylogenetically diverse hexapod orders are now known to contain species that possess two engrailed-family paralogues, with in each case only one paralogue encoding the RS-motif, a characteristic feature of holometabolous insect invected proteins. We propose that the homeoboxes of hexapod engrailed-family paralogues are evolving in a concerted fashion, resulting in gene trees that overestimate the frequency of gene duplication. We present new phylogenetic analyses using non-homeodomain amino acid sequence that support this view. The S. gregaria engrailed-family paralogues provide strong evidence that concerted evolution might in part be explained by recurrent gene conversion. Finally, we hypothesize that the RS-motif is part of a serine-rich domain targeted for phosphorylation.

Basal relationships in the Drosophila melanogaster species group

The Drosophila melanogaster species group is a popular model for evolutionary studies due to its morphological and ecological diversity and its inclusion of the model species D. melanogaster. However, phylogenetic relationships among major lineages within this species group remain controversial. In this report, the phylogeny of 10 species representing each of the well-supported monophyletic clades in the melanogaster group was studied using the sequences of 14 loci that together comprise 9493 nucleotide positions. Combined Bayesian analysis using gene-specific substitution models produced a 100% credible set of two trees. In the strict consensus of these trees, the ananassae subgroup branches first in the melanogaster species group, followed by the montium subgroup. The remaining lineages form a monophyletic clade in which D. ficusphila and D. elegans branch first, followed by D. biarmipes, D. eugracilis, and the melanogaster subgroup. This strongly supported phylogeny resolves most basal relationships in the melanogaster species group, and provides a framework that can be extended in the future to encompass more species.

Inhibition of neurogenesis at the zebrafish midbrain-hindbrain boundary by the combined and dose-dependent activity of a new hairy/E(spl) gene pair

The intervening zone (IZ) is a pool of progenitor cells located at the midbrain-hindbrain boundary (MHB) and important for MHB maintenance, midbrain-hindbrain growth and the generation of midbrain-hindbrain neurons. Recently, we implicated the Hairy/E(spl) transcription factor Her5 in the formation of the medial (most basal) part of the IZ (MIZ) in zebrafish; the molecular bases for lateral IZ (LIZ) formation, however, remain unknown. We now demonstrate that her5 is physically linked to a new family member, him, displaying an identical MHB expression pattern. Using single and double knockdowns of him and her5, as well as a him+her5 deletion mutant background (b404), we demonstrate that Him and Her5 are equally necessary for MIZ formation, and that they act redundantly in LIZ formation in vivo. We show that these processes do not involve cross-regulation between Him and Her5 expression or activities, although Him and Her5 can heterodimerize with high affinity. Increasing the function of one factor when the other is depleted further shows that Him and Her5 are functionally interchangeable. Together, our results demonstrate that patterning and neurogenesis are integrated by the her5-him gene pair to maintain a progenitor pool at the embryonic MHB. We propose a molecular mechanism for this process where the global 'Him+Her5' activity inhibits ngn1 expression in a dose-dependent manner and through different sensitivity thresholds along the medio-lateral axis of the neural plate.

Sequence variation of alcohol dehydrogenase (Adh) paralogs in cactophilic Drosophila

This study focuses on the population genetics of alcohol dehydrogenase (Adh) in cactophilic Drosophila. Drosophila mojavensis and D. arizonae utilize cactus hosts, and each host contains a characteristic mixture of alcohol compounds. In these Drosophila species there are two functional Adh loci, an adult form (Adh-2) and a larval and ovarian form (Adh-1). Overall, the greater level of variation segregating in D. arizonae than in D. mojavensis suggests a larger population size for D. arizonae. There are markedly different patterns of variation between the paralogs across both species. A 16-bp intron haplotype segregates in both species at Adh-2, apparently the product of an ancient gene conversion event between the paralogs, which suggests that there is selection for the maintenance of the intron structure possibly for the maintenance of pre-mRNA structure. We observe a pattern of variation consistent with adaptive protein evolution in the D. mojavensis lineage at Adh-1, suggesting that the cactus host shift that occurred in the divergence of D. mojavensis from D. arizonae had an effect on the evolution of the larval expressed paralog. Contrary to previous work we estimate a recent time for both the divergence of D. mojavensis and D. arizonae (2.4 +/- 0.7 MY) and the age of the gene duplication (3.95 +/- 0.45 MY).

Codon bias differentiates between the duplicated amylase loci following gene duplication in Drosophila

We examined the pattern of synonymous substitutions in the duplicated Amylase (Amy) genes (called the Amy1- and Amy3-type genes, respectively) in the Drosophila montium species subgroup. The GC content at the third synonymous codon sites of the Amy1-type genes was higher than that of the Amy3-type genes, while the GC content in the 5'-flanking region was the same in both genes. This suggests that the difference in the GC content at third synonymous sites between the duplicated genes is not due to the temporal or regional changes in mutation bias. We inferred the direction of synonymous substitutions along branches of a phylogeny. In most lineages, there were more synonymous substitutions from G/C (G or C) to A/T (A or T) than from A/T to G/C. However, in one lineage leading to the Amy1-type genes, which is immediately after gene duplication but before speciation of the montium species, synonymous substitutions from A/T to G/C were predominant. According to a simple model of synonymous DNA evolution in which major codons are selectively advantageous within each codon family, we estimated the selection intensity for specific lineages in a phylogeny on the basis of inferred patterns of synonymous substitutions. Our result suggested that the difference in GC content at synonymous sites between the two Amy-type genes was due to the change of selection intensity immediately after gene duplication but before speciation of the montium species.

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.

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

In this study, we randomly sampled Drosophila melanogaster from Japanese and Kenyan natural populations. We sequenced duplicated (proximal and distal) Amy gene regions to test whether the patterns of polymorphism were consistent with neutral molecular evolution. F(st) between the two geographically distant populations, estimated from Amy gene regions, was 0.084, smaller than reported values for other loci, comparing African and Asian populations. Furthermore, little genetic differentiation was found at a microsatellite locus (DROYANETSB) in these samples (G'st = -0.018). The results of several tests (Tajima's, Fu and Li's, and Wall's tests) were not significantly different from neutrality. However, a significantly higher level of fixed replacement substitutions was detected by a modified McDonald and Kreitman test for both populations. This indicates that positive selection occurred during or immediately after the speciation of D. melanogaster. Sliding-window analysis showed that the proximal region 1, a part of the proximal 5' flanking region, was conserved between D. melanogaster and its sibling species, D. simulans. An HKA test was significant when the proximal region 1 was compared with the 5' flanking region of Alcohol dehydrogenase (Adh), indicating a severe selective constraint on the Amy proximal region 1. These results suggest that natural selection has played an important role in the molecular evolution of Amy gene regions in D. melanogaster.

Resistance gene complexes: evolution and utilization

More than 30 genes have been characterized from different plant species that provide resistance to a variety of different pathogen and pest species. The structures of most are consistent with a role in pathogen recognition and defense response signaling. Resistance genes are very abundant in plant genomes and most belong to tightly linked gene families. Evolution of R genes is driven by selection on allelic variation created by mutation and re-assorted by recombination between alleles and sometimes between different gene family members. Selection favors genes that can recognize pathogen avr gene products that are present in pathogen populations. Selection at linked gene families favors haplotypes with useful combinations of genes but a limited physiological cost to the plant. Future utilization of R genes will include transfer between related genera and identification or construction of genes that condition durable resistance to variable pathogens. Genes with durable resistance may interact with conserved pathogen elicitors or condition resistance responses that are independent of specific Avr gene interactions.

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.

Sequence and phylogenetic analysis of the SNF4/AMPK gamma subunit gene from Drosophila melanogaster

To optimize gene expression under different environmental conditions, many organisms have evolved systems which can quickly up- and down-regulate the activity of other genes. Recently, the SNF1 kinase complex from yeast and the AMP-activated protein kinase complex from mammals have been shown to represent homologous metabolic sensors that are key to regulating energy levels under times of metabolic stress. Using heterologous probing, we have cloned the Drosophila melanogaster homologue of SNF4, the noncatalytic effector subunit from this kinase complex. A sequence corresponding to the partial genomic sequence as well as the full-length cDNA was obtained, and shows that the D. melanogaster SNF4 is encoded in a 1944-bp cDNA representing a protein of 648 amino acids (aa). Southern analysis of Drosophila genomic DNA in concert with a survey of mammalian SNF4 ESTs indicates that in metazoans, SNF4 is a duplicated gene, and possibly even a larger gene family. We propose that one gene copy codes for a short (330 aa) protein, whereas the second locus codes for a longer version (<410 aa) that is extended at the carboxy terminus, as typified by the Drosophila homologue presented here. Phylogenetic analysis of yeast, invertebrate, and multiple mammalian isoforms of SNF4 shows that the gene duplication likely occurred early in the metazoan lineage, as the protein products of the different loci are relatively divergent. When the phylogeny was extended beyond the SNF4 gene family, SNF4 shares sequence similarity with other cystathionine-β-synthase domain-containing proteins, including IMP dehydrogenase and a variety of uncharacterized Methanococcus proteins.Key words: SNF4, AMPK gamma subunit, derepression, gene family, phylogeny.

Molecular analysis of the intergenic region of the duplicated Amy genes of Drosophila melanogaster and Drosophila teissieri

The intergenic regions between the duplicated amylase coding regions (Amy) of D. melanogaster and D. teissieri were sequenced. Their lengths in D. melanogaster and D. teissieri were 4,536 bp and 4,621 bp, respectively. Since homology between the upstream regions of the two duplicated genes was found up to 450 bp from the initiation codon of the Amy genes, the ancestral Amy coding region duplicated together with at least 450 bp of the 5'-flanking region as one unit. Comparison of the regions between the two species revealed that the level of divergence was very heterogeneous. Although the mean level of the nucleotide difference in this region was 0.107, no nucleotide substitution was found in four subregions whose sizes were more than 100 bp. Since the probability of these four subregions being completely conserved between D. melanogaster and D. teissieri was very low, these subregions were considered to have relatively important roles in evolution. Large insertions and deletions were not observed in this region but small ones were observed all over the region except for an about 1-kb subregion. This 1-kb region corresponded to an open reading frame encoding a protein which had some sequence identity with the proteins of the serine protease inhibitor superfamily (serpin). Since we could find a transcript of this gene and the synonymous substitution rate was higher than the replacement substitution rate, we suggest that this gene encodes an active serpin in Drosophila.

Is esterase-P encoded by a cryptic pseudogene in Drosophila melanogaster?

We have amplified and sequenced the gene encoding Esterase-P (Est-P) in 10 strains of Drosophila melanogaster. Three premature termination codons occur in the coding region of the gene in two strains. This observation, together with other indirect evidence, leads us to propose that Est-P may be a pseudogene in D. melanogaster. Est-P would be a "cryptic" pseudogene, in the sense that it retains intact the coding sequence (without stop codons and other alterations usually observed in pseudogenes) in most D. melanogaster strains. We conjecture that the beta-esterase cluster may consist in other Drosophila species of functional and nonfunctional genes. We also conjecture that the rarity of detected pseudogenes in Drosophila may be due to the difficulty of discovering them, because most of them are cryptic.

Precise limitation of concerted evolution to ORFs in mosquito Hsp82 genes

Two Hsp82 genes were isolated from the malaria vector Anopheles albimanus in a single lambda phage clone. The two genes are in a head-to-head arrangement separated by approx. 0.9 kbp. Northern hybridizations and 5' RACE demonstrate that both genes are transcribed, have moderate levels of constitutive transcription, and are also heat-inducible with maximum transcript accumulation occurring after 40 degrees C heat shocks. Both genes have typical heat-shock promoters and conserved intron boundaries in the untranslated leaders. The open reading frames are 99.6% identical differing in only nine silent nucleotide positions in the 2166 bp ORFs. However, precisely outside the ORFs, the flanking DNA of the two genes shows no evidence of common derivation. The high degree of identity between the two ORFs appears to be a result of gene conversion occurring by a process similar to that previously suspected in the A. albimanus Hsp70 genes and several D. melanogaster genes arranged as palindromes. This process probably involves a stem-loop intermediate and is restricted in extent by flanking sequence divergence. These Hsp82 genes clearly demonstrate the extreme precision with which gene conversion can lead to protein-coding-region homogeneity yet allow flanking DNA divergence.

Evolutionary relationships and sequence variation of alpha-amylase variants encoded by duplicated genes in the Amy locus of Drosophila melanogaster

To infer the genealogical relationships of alpha-amylase electromorphs of Drosophila melanogaster, we determined the nucleotide sequences of a collection of electromorphs sampled throughout the world. On average there were 1.0 amino acid substitutions between identical electromorphs and 3.9 between different electromorphs, respectively. We found that the evolution of AMY1 through AMY6 electromorphs occurred by sequential accumulation of single amino acid substitutions each causing one charge difference. The nucleotide diversities at synonymous sites within Amy1,Amy2,Amy3,Amy4 and Amy6 were 0.0321, 0.0000, 0.0355, 0.0059 and 0.0030, respectively. We also obtained evidence of genetic exchanges, such as intrachromosomal recombination, interchromosomal recombination or gene conversion, between the two duplicated Amy genes as well as among the alleles.

Molecular evolution of the duplicated Amy locus in the Drosophila melanogaster species subgroup: concerted evolution only in the coding region and an excess of nonsynonymous substitutions in speciation

From the analysis of restriction maps of the Amy region in eight sibling species belonging to the Drosophila melanogaster species subgroup, we herein show that the patterns of duplication of the Amy gene are almost the same in all species. This indicates that duplication occurred before speciation within this species subgroup. From the nucleotide sequence data, we show a strong within-species similarity between the duplicated loci in the Amy coding region. This is in contrast to a strong similarity in the 5' and 3' flanking regions within each locus (proximal or distal) throughout the species subgroup. This means that concerted evolution occurred only in the Amy coding region and that differentiated evolution between the duplication occurred in the flanking regions. Moreover, when comparing the species, we also found a significant excess of nonsynonymous substitutions. In particular, all the fixed substitutions specific to D. erecta were found to be nonsynonymous. We thus conclude that adaptive protein evolution occurred in the lineage of D. erecta that is a "specialist" species for host plants and probably also occurs in the process of speciation in general.

Functional conservation of a glucose-repressible amylase gene promoter from Drosophila virilis in Drosophila melanogaster

Previous studies have demonstrated that the expression of the alpha-amylase gene is repressed by dietary glucose in Drosophila melanogaster. Here, we show that the alpha-amylase gene of a distantly related species, D. virilis, is also subject to glucose repression. Moreover, the cloned amylase gene of D. virilis is shown to be glucose repressible when it is transiently expressed in D. melanogaster larvae. This cross-species, functional conservation is mediated by a 330-bp promoter region of the D. virilis amylase gene. These results indicate that the promoter elements required for glucose repression are conserved between distantly related Drosophila species. A sequence comparison between the amylase genes of D. virilis and D. melanogaster shows that the promoter sequences diverge to a much greater degree than the coding sequences. The amylase promoters of the two species do, however, share small clusters of sequence similarity, suggesting that these conserved cis-acting elements are sufficient to control the glucose-regulated expression of the amylase gene in the genus Drosophila.

The salivary glands of the vector mosquito, Aedes aegypti, express a novel member of the amylase gene family

Several cDNA clones with similarity to alpha-amylases have been characterized from a library made from adult female salivary gland RNA isolated from the vector mosquito, Aedes aegypti. The corresponding gene, designated Amylase I (Amy I), is expressed specifically in the proximal-lateral lobes of the adult female salivary gland, a pattern overlapping that of another gene, Mal I, involved in carbohydrate metabolism. The deduced amino acid sequence of Amy I indicates that this gene encodes a protein, approximate M(r) = 81,500, that appears to be a novel member of the amylase gene family. The mosquito protein contains a putative signal peptide for secretion and several consensus sites for asparagine-linked glycosylation. The Amy I protein shows significant similarity to invertebrate and vertebrate amylases including the conservation of four reactive and substrate binding sites. However, the amino-terminal region of the Amy-I protein is unique to the mosquito. Similarity with the Drosophila melanogaster protein is evident only after the first 260 amino acids in the mosquito sequence. The identification of this gene and its expression pattern adds to the observed relationship between spatial-specific gene expression in the female salivary glands and the specific feeding mode of the adult mosquito.

Two genes encoding 1-aminocyclopropane-1-carboxylate synthase in zucchini (Cucurbita pepo) are clustered and similar but differentially regulated

A 17-kilobase (kb) region of the zucchini (Cucurbita pepo) genome has been sequenced and contains two genes, CP-ACC1A and CP-ACC1B, encoding 1-aminocyclopropane-1-carboxylate synthase (ACC synthase; S-adenosyl-L-methionine methylthioadenosine-lyase, EC 4.4.1.14). The genes are transcribed convergently and are separated by a 5.7-kb intergenic region. Their coding regions are interrupted by four introns located in identical positions. While the DNA identity in their coding regions is 97%, their 5' and 3' flanking regions are highly divergent. Transcription of CP-ACCIA is rapidly induced by wounding in fruit and etiolated hypocotyls and by indoleacetic acid (IAA)/benzyladenine/LiCl only in fruit tissue. Conditions that induce CP-ACC1B expression have not been found. Protein synthesis inhibition derepresses the expression of CP-ACC1A and other unidentified ACC synthase genes, suggesting that they may be under negative control. The amino acid sequences deduced from the nucleotide sequences of the genes are 493 and 494 residues long with 95% identity. The most notable feature of the amino acid sequence is the presence of 11 of the 12 invariant amino acid residues involved in the binding of the substrate and pyridoxal-5'-phosphate in various aminotransferases. We conclude that ACC synthase is encoded by a multigene family of which certain members are differentially induced by auxin in a tissue-specific manner. Furthermore, ACC synthase, a pyridoxal-containing enzyme, may have an evolutionary relationship with the superfamily of aminotransferases.