The Evolution of Duplicate Glyceraldehyde-3-Phosphate Dehydrogenase Genes in Drosophila

In Drosophila melanogaster there are two genes which encode the enzyme glyceraldehyde-3-phosphate dehydrogenase (GAPDH), Gapdh-43E and Gapdh-13F. We have shown that Gapdh-43E codes for the GAPDH subunit with an apparently larger molecular weight while Gapdh-13F encodes the GAPDH subunit having an apparently smaller molecular weight. Immunoblots of sodium dodecyl sulfate gels were used to survey species from throughout the genus and results indicated that two classes of GAPDH subunits are present only in Drosophila species of the melanogaster and takahashi subgroups of the melanogaster group. Only the smaller subunit is found in species of the obscura group while all other species have only a large subunit. Drosophila hydei was analyzed at the DNA level as a representative species of the subgenus Drosophila. The genome of this species has a single Gapdh gene which is localized at a cytogenetic position likely to be homologous to Gapdh-43 E of D. melanogaster. Comparison of its sequence with the sequence of the D. melanogaster Gapdh genes indicates that the two genes of D. melanogaster are more similar to one another than either is to the gene from D. hydei. The Gapdh gene from D. hydei contains an intron following codon 29. Neither Gapdh gene of D. melanogaster has an intron within the coding region. Southern blots of genomic DNA were used to determine which species have duplicate Gapdh genomic sequences. Gene amplification was used to determine which species have a Gapdh gene that is interrupted by an intron. Species of the subgenus Drosophila have a single Gapdh gene with an intron. Species of the willistoni and saltans groups have a single Gapdh gene that does not contain an intron.(ABSTRACT TRUNCATED AT 250 WORDS)     

Evolution of the Adh locus in the Drosophila willistoni group: the loss of an intron, and shift in codon usage

We report here the DNA sequence of the alcohol dehydrogenase gene (Adh) cloned from Drosophila willistoni. The three major findings are as follows: (1) Relative to all other Adh genes known from Drosophila, D. willistoni Adh has the last intron precisely deleted; PCR directly from total genomic DNA indicates that the deletion exists in all members of the willistoni group but not in any other group, including the closely related saltans group. Otherwise the structure and predicted protein are very similar to those of other species. (2) There is a significant shift in codon usage, especially compared with that in D. melanogaster Adh. The most striking shift is from C to U in the wobble position (both third and first position). Unlike the codon-usage-bias pattern typical of highly biased genes in D. melanogaster, including Adh, D. willistoni has nearly 50% G + C in the third position. (3) The phylogenetic information provided by this new sequence is in agreement with almost all other molecular and morphological data, in placing the obscura group closer to the melanogaster group, with the willistoni group farther distant but still clearly within the subgenus Sophophora.      

Macro-evolution of the hairy enhancer in Drosophila species

It has been suggested that many of the changes in the developmental program might be in the cis-acting promoters and enhancer regions. Here I study the macro-evolutionary changes of an enhancer region for the early developmental gene hairy in Drosophila melanogaster, D. simulans, D. pseudoobscura, D. willistoni, D. nebulosa, D. hydei, and D. virilis. The enhancer region is characterized by small, highly conserved blocks interspersed among highly variable regions. Nevertheless, species phylogenies constructed by the enhancer sequences agree with the widely accepted phylogeny of these species. The evolution of the variable regions is consistent with a molecular clock, while the evolution of the conserved blocks is significantly different from a clock. In particular, the D. pseudoobscura lineage shows the highest degree of species-specific change consistent with changes in expression timing reported in an earlier study. It has been suggested that the variation in sequence length between highly conserved blocks may play a role in the coordination of regulatory processes, such as protein-protein interactions; thus, stabilizing selection has been suggested to act on the length variations. Here I develop a test for stabilizing selection on length variation and show that the hairy enhancer does not show statistically significant evidence for stabilizing selection. J. Exp. Zool. (Mol. Dev. Evol.) 291:175-185, 2001.     

Evidence for Horizontal Transmission of the P Transposable Element between Drosophila Species

Several studies have suggested that P elements have rapidly spread through natural populations of Drosophila melanogaster within the last four decades. This observation, together with the observation that P elements are absent in the other species of the melanogaster subgroup, has lead to the suggestion that P elements may have entered the D. melanogaster genome by horizontal transmission from some more distantly related species. In an effort to identify the potential donor in the horizontal transfer event, we have undertaken an extensive survey of the genus Drosophila using Southern blot analysis. The results showed that P-homologous sequences are essentially confined to the subgenus Sophophora. The strongest P hybridization occurs in species from the closely related willistoni group. A wild-derived strain of D. willistoni was subsequently selected for a more comprehensive molecular examination. As part of the analysis, a complete P element was cloned and sequenced from this line. Its nucleotide sequence was found to be identical to the D. melanogaster canonical P, with the exception of a single base substitution at position 32. When the cloned element was injected into D. melanogaster embryos, it was able to both promote transposition of a coinjected marked transposon and induce singed-weak mutability, thus demonstrating its ability to function as an autonomous element. The results of this study suggest that D. willistoni may have served as the donor species in the horizontal transfer of P elements to D. melanogaster.     

Interspecific Comparison of the Transformer Gene of Drosophila Reveals an Unusually High Degree of Evolutionary Divergence

The transformer (tra) gene of Drosophila melanogaster occupies an intermediate position in the regulatory pathway controlling all aspects of somatic sexual differentiation. The female-specific expression of this gene's function is regulated by the Sex lethal (Sxl) gene, through a mechanism involving sex-specific alternative splicing of tra pre-mRNA. The tra gene encodes a protein that is thought to act in conjunction with the transformer-2 (tra-2) gene product to control the sex-specific processing of doublesex (dsx) pre-mRNA. The bifunctional dsx gene carries out opposite functions in the two sexes, repressing female differentiation in males and repressing male differentiation in females. Here we report the results from an evolutionary approach to investigate tra regulation and function, by isolating the tra-homologous genes from selected Drosophila species, and then using the interspecific DNA sequence comparisons to help identify regions of functional significance. The tra-homologous genes from two Sophophoran subgenus species, Drosophila simulans and Drosophila erecta, and two Drosophila subgenus species, Drosophila hydei and Drosophila virilis, were cloned, sequenced and compared to the D. melanogaster tra gene. This comparison reveals an unusually high degree of evolutionary divergence among the tra coding sequences. These studies also highlight a highly conserved sequence within intron one that probably defines a cis-acting regulator of the sex-specific alternative splicing event.     

Constraints on Intron Evolution in the Gene Encoding the Myosin Alkali Light Chain in Drosophila

Interspecific comparisons of intron sequences reveal conserved blocks of invariant nucleotides and several other departures from the strictly neutral model of molecular evolution. To distinguish the past action of evolutionary forces in introns known to have regulatory information, we examined nucleotide sequence variation at 991 sites in a random sample of 16 Drosophila melanogaster alleles of the gene encoding the myosin alkali light chain (Mlc1). The Mlc1 gene of D. melanogaster encodes two MLC1 isoforms via developmentally regulated alternative pre-mRNA splicing. Analyses of these data reveal that introns 4 and 5, which flank the alternatively spliced exon 5, have reduced levels of both intraspecific polymorphism and interspecific divergence relative to intron 3. No polymorphism was observed in any of the exons examined in D. melanogaster. A genealogical analysis clearly demonstrates the occurrence of intragenic recombination in the ancestral history of Mlc1. Recombination events are estimated to be 13 times more likely than mutation events over the span of the sequenced region. Although there is little evidence for pairwise linkage disequilibrium in the Mlc1 region, higher order disequilibrium does seem to be present in the 5' half of the portion of the gene that was examined. Predictions of the folding free energy of the pre-mRNA reveal that sampled alleles have a significantly higher (less stable) free energy than do randomly permuted sequences. These results are consistent with the hypothesis that introns surrounding an alternatively spliced exon are subjected to additional constraints, perhaps due to specific aspects of secondary structure required for appropriate splicing of the pre-mRNA molecule.     

Intron length evolution in Drosophila

I present data on the evolution of intron lengths among 3 closely related Drosophila species, D. melanogaster, Drosophila simulans, and Drosophila yakuba. Using D. yakuba as an outgroup, I mapped insertion and deletion mutations in 148 introns (spanning approximately 30 kb) to the D. melanogaster and D. simulans lineages. Intron length evolution in the 2 sister species has been different: in D. melanogaster, X-linked introns have increased slightly in size, whereas autosomal ones have decreased slightly in size; in D. simulans, both X-linked and autosomal introns have decreased in size. To understand the possible evolutionary causes of these lineage- and chromosome-specific patterns of intron evolution, I studied insertion-deletion (indel) polymorphism and divergence in D. melanogaster. Small insertion mutations segregate at elevated frequencies and enjoy elevated probabilities of fixation, particularly on the X chromosome. In contrast, there is no detectable X chromosome effect on fixations in D. simulans. These findings suggest X chromosome-specific selection or biased gene conversion-gap repair favoring insertions in D. melanogaster but not in D. simulans. These chromosome- and lineage-specific patterns of indel substitution are not easily explained by existing general population genetic models of intron length evolution. Genomic data from D. melanogaster further suggest that the forces described here affect introns and intergenic regions similarly.     

Identification and experimental validation of splicing regulatory elements in Drosophila melanogaster reveals functionally conserved splicing enhancers in metazoans

RNA sequence elements involved in the regulation of pre-mRNA splicing have previously been identified in vertebrate genomes by computational methods. Here, we apply such approaches to predict splicing regulatory elements in Drosophila melanogaster and compare them with elements previously found in the human, mouse, and pufferfish genomes. We identified 99 putative exonic splicing enhancers (ESEs) and 231 putative intronic splicing enhancers (ISEs) enriched near weak 5' and 3' splice sites of constitutively spliced introns, distinguishing between those found near short and long introns. We found that a significant proportion (58%) of fly enhancer sequences were previously reported in at least one of the vertebrates. Furthermore, 20% of putative fly ESEs were previously identified as ESEs in human, mouse, and pufferfish; while only two fly ISEs, CTCTCT and TTATAA, were identified as ISEs in all three vertebrate species. Several putative enhancer sequences are similar to characterized binding-site motifs for Drosophila and mammalian splicing regulators. To provide additional evidence for the function of putative ISEs, we separately identified 298 intronic hexamers significantly enriched within sequences phylogenetically conserved among 15 insect species. We found that 73 putative ISEs were among those enriched in conserved regions of the D. melanogaster genome. The functions of nine enhancer sequences were verified in a heterologous splicing reporter, demonstrating that these sequences are sufficient to enhance splicing in vivo. Taken together, these data identify a set of predicted positive-acting splicing regulatory motifs in the Drosophila genome and reveal regulatory sequences that are present in distant metazoan genomes.     

Nomadic enhancers: tissue-specific cis-regulatory elements of yellow have divergent genomic positions among Drosophila species

cis-regulatory DNA sequences known as enhancers control gene expression in space and time. They are central to metazoan development and are often responsible for changes in gene regulation that contribute to phenotypic evolution. Here, we examine the sequence, function, and genomic location of enhancers controlling tissue- and cell-type specific expression of the yellow gene in six Drosophila species. yellow is required for the production of dark pigment, and its expression has evolved largely in concert with divergent pigment patterns. Using Drosophila melanogaster as a transgenic host, we examined the expression of reporter genes in which either 5' intergenic or intronic sequences of yellow from each species controlled the expression of Green Fluorescent Protein. Surprisingly, we found that sequences controlling expression in the wing veins, as well as sequences controlling expression in epidermal cells of the abdomen, thorax, and wing, were located in different genomic regions in different species. By contrast, sequences controlling expression in bristle-associated cells were located in the intron of all species. Differences in the precise pattern of spatial expression within the developing epidermis of D. melanogaster transformants usually correlated with adult pigmentation in the species from which the cis-regulatory sequences were derived, which is consistent with cis-regulatory evolution affecting yellow expression playing a central role in Drosophila pigmentation divergence. Sequence comparisons among species favored a model in which sequential nucleotide substitutions were responsible for the observed changes in cis-regulatory architecture. Taken together, these data demonstrate frequent changes in yellow cis-regulatory architecture among Drosophila species. Similar analyses of other genes, combining in vivo functional tests of enhancer activity with in silico comparative genomics, are needed to determine whether the pattern of regulatory evolution we observed for yellow is characteristic of genes with rapidly evolving expression patterns.     

Convergent intron loss of MRP1 in Drosophila and mosquito species

Previous study revealed that the MRP1 gene ortholog DMRP1/CG6214 of Drosophila melanogaster contains 12 exons in the coding region. In the current study, the genes of DMRP1/CG6214 from D. melanogaster and Drosophila virilis were compared, and the result indicated that D. virilis had an extra intron located in exon 2, implying that intron loss or gain might have occurred at this locus. To track the evolution of the extra intron (Intron Z), orthologous nucleotide sequences of 37 arthropod species were cloned or annotated. Based on phylogenetic analysis, we found that Intron Z should present in the common ancestor of arthropod species, more than 420 Ma. In addition, we found that Sophophora subgenus species and mosquito (Culex pipiens) lost Intron Z independently, showing evolutionary convergence.     

Interspecific and intraspecific comparisons of the period locus in the Drosophila willistoni sibling species

The period (per) locus has received much attention in molecular evolution studies because it is one of the best studied "behavioral genes" and because it offers insight into the evolution of repetitive sequences. We studied most of the coding region of per in Drosophila willistoni and confirmed previously observed patterns of conservation and divergence among distantly related species. Five regions are so highly diverged that they cannot be aligned, whereas a region encompassing the PAS domain is very conserved. Structural and nucleotide polymorphism patterns in the willistoni group are not the same as those observed in previously studied species. We sequenced the region homologous to the highly polymorphic threonine-glycine repeat of D. melanogaster in multiple strains of D. willistoni, as well as in other members of willistoni group, and found an unusual amount of conservation in this region. However, the next nonconserved region downstream in the sequence is quite variable and polymorphic for the number of repeated glycines. The glycine codon usage is significantly different in this glycine repeat as compared to other parts of the gene. We were able to plot the directionality of change in the glycine repeat region onto a phylogeny and find that the addition of glycines is the general trend with the diversification of the willistoni group.      

Distribution of transposable elements in neotropical species of Drosophila

The phylogenetic distribution of transposable families, P, gypsy, hobo, I, and mariner has been analyzed in 33 species of 11 groups of neotropical Drosophila and a Drosophilidae species Zygotrica vittimaculosa, using squash blot and dot blot. Genomic DNA of almost all neotropical species tested hybridized with gypsy probe and some species showed a particularly strong hybridization signal, as D. gaucha, D. virilis, and species of flavopilosa group. The hobo element was restricted to melanogaster group and some strains of D. willistoni. Only D. simulans DNA showed hybridization to mariner probe in all species tested and D. simulans and D. melanogaster showed hybridization with I element probe. P element homologous sequence was present in D. melanogaster and all species and strains of the willistoni and saltans groups tested. The presence of at least one P-homologous sequence was detected in Drosophila mediopunctata. This one was the only P-bearing species of all six tested from the tripunctata group. Four different pairs of primers homologous to segments of the canonical sequence of D. melanogaster's P were used to amplify specific sequences from D. mediopunctata DNA, showing the occurrence of seemingly well-conserved P-homologous sequences. This revised version was published online in July 2006 with corrections to the Cover Date.     

Introns and their flanking sequences of Bombyx mori rDNA.

We obtained two different clones (16 kb and 13 kb) of B. mori rDNA with intron sequence within the 28S-rRNA coding region. The sequence surrounding the intron was found to be highly conserved as indicated in several eukaryotes (Tetrahymena, Drosophila and Xenopus). The 28S rRNA-coding sequence of 16 kb and 13 kb clone was interrupted at precisely the same sites as those where the D. melanogaster rDNA interrupted by the type I and type II intron, respectively. The intron sequences of B. mori were different from those of D. melanogaster. In 16 kb clone, the intron was flanked by 14 bp duplication of the junction sequence, which was also present once within the 28S rRNA-coding region of rDNA without intron. This 14 bp sequence was identical with those surrounding the introns of Dipteran rDNAs.     

Molecular Evolution of the Metallothionein Gene Mtn in the Melanogaster Species Group: Results from Drosophila Ananassae

Three distinctly different alleles of the metallothionein gene Mtn have been identified in natural Drosophila melanogaster populations: Mtn(.3), Mtn(1), and Dp(Mtn(1)), where the latter designates a tandem duplication of Mtn(1). In Drosophila simulans, only Mtn(.3)-type alleles have been found. It has been suggested that Mtn(.3) is the ancestral allele and demonstrated that a presumed two-step transition from Mtn(.3) to Mtn(1) to Dp(Mtn(1)) is accompanied by an approximate 5-fold increase in RNA levels. We analyzed the evolutionary genetics of the Mtn locus of Drosophila ananassae, a distant relative of D. melanogaster and D. simulans within the melanogaster species group. The Mtn gene of D. ananassae is most similar to Mtn(.3). (i) it is identical with Mtn(.3) at the amino acid level, but differs from Mtn(1) in its terminal codon; (ii) its 3'' UTR contains a characteristic extra DNA segment of about 50 bp which is present in Mtn(.3), but lacking in Mtn(1); (iii) duplications of Mtn were not found in a worldwide sample of 110 wild D. ananassae chromosomes. However, the intron of the Mtn gene in D. ananassae is only 69 bp long, whereas the length of the Mtn(.3) and Mtn(1) introns is 265 bp; and it lacks a polypyrimidine stretch upstream of the 3'' splice site in contrast to the much greater pyrimidine-richness found in the Mtn(.3) and Mtn(1) introns. A short intron (67 bp) was also identified in a D. pseudoobscura Mtn allele, suggesting that the short intron is the ancestral form and that the transition from the short to the long intron occurred within the melanogaster species group. We discuss the significance of this observation with regard to the recently proposed classification of D. melanogaster introns into two groups: short introns (<90 bp) which tend to lack polypyrimidine stretches, and longer ones which have strong 3'' splice signals similar to mammalian introns. A database search revealed that this length dimorphism is an evolutionarily conserved feature of Drosophila introns; transitions from one size class to the other appear to be rare between closely related species (e.g., within the melanogaster subgroup).