Isolation and characterization of Drosophila melanogaster U2 small nuclear RNA genes

We describe here the organization of DNA sequences complementary to Drosophila melanogaster U2 small nuclear (sn) RNA. From a genomic library we isolated two recombinants containing two genes each. Genomic reconstruction experiments and Southern analysis revealed that D. melanogaster possesses only four to five U2 snRNA genes or very closely related sequences. The nucleotide sequence of one of the clones analysed shows 77% homology with rat U2 snRNA. A stretch of 12 nucleotides that has been implicated in heterogeneous nuclear RNA splicing is conserved between rat and Drosophila. The genomic organization of these genes is very similar in different melanogaster strains but diverges highly in different Drosophila species.     

A complete and a truncated U1 snRNA gene of Drosophila melanogaster are found as inverted repeats at region 82E of the polytene chromosomes.

A phage containing two sequences homologous to U1 snRNA was isolated from a Drosophila melanogaster genomic library, and identified with a previously cloned D. melanogaster U1 snRNA gene. DNA sequence analysis showed that complete and truncated U1 snRNA genes are present, both of which have base substitutions relative to U1 snRNA. These genes show conservation of 5' and 3' flanking regions relative to other U1 and U2 snRNA genes of Drosophila. Intramolecular renaturation experiments and electron microscope mapping demonstrates that the two U1 snRNA sequences are present as inverted repeats about 2.7kb apart, separated by a smaller pair of inverted repeats of an unrelated sequence. These U1 snRNA sequences were located by in situ hybridization at 82E, and related sequences were found at 21D and 95C on the polytene chromosome map. The results are discussed with reference to the origin and function of snRNAs.     

Evolutionary rearrangement of the amylase genomic regions between Drosophila melanogaster and Drosophila pseudoobscura

Two Drosophila pseudoobscura genomic clones have sequence similarity to the Drosophila melanogaster amylase region that maps to the 53CD region on the D. melanogaster cytogenetic map. The two clones with similarity to amylase map to sections 73A and 78C of the D. pseudoobscura third chromosome cytogenetic map. The complete sequences of both the 73A and 78C regions were compared to the D. melanogaster genome to determine if the coding region for amylase is present in both regions and to determine the evolutionary mechanism responsible for the observed distribution of the amylase gene or genes. The D. pseudoobscura 73A and 78C linkage groups are conserved with the D. melanogaster 41E and 53CD regions, respectively. The amylase gene, however, has not maintained its conserved linkage between the two species. These data indicate that amylase has moved via a transposition event in the D. melanogaster or D. pseudoobscura lineage. The predicted genes within the 73A and 78C regions show patterns of molecular evolution in synonymous and nonsynonymous sites that are consistent with previous studies of these two species.     

Drosophila melanogaster U1 snRNA genes

We have isolated and characterized a recombinant which contains a Drosophila melanogaster U1 small nuclear RNA (snRNA) gene colinear with the published snRNA sequence. Southern hybridizations of the fly genomic DNA, using as probe a plasmid containing only the coding region of the gene, shows that the fly contains at most three or four genes and very few related sequences for the small nuclear U1 RNA. These genes were localized by in situ hybridization at different chromosomal loci and show no spatial relationship to the U2 snRNA genes.     

Conservation of gene order, structure and sequence between three closely linked genes in Drosophila melanogaster and Drosophila virilis

In Drosophila melanogaster, the apparently unrelated genes anon-66Da, RpL14, and anon-66Db (from telomere to centromere) are located on a 5547 bp genomic fragment on chromosome arm 3L at cytological position 66D8. The three genes are tightly linked, and flanked by two relatively large genes with unknown function. We have taken a comparative genomic approach to investigate the evolutionary history of the three genes. To this end we isolated a Drosophila virilis 7.3 kb genomic fragment which is homologous to a 5.5 kb genomic region of D. melanogaster. Both fragments map to Muller's element D, namely to section 66D in D. melanogaster and to section 32E in D. virilis, and harbor the genes anon-66Da, RpL14, and anon-66Db. We demonstrate that the three genes exhibit a high conservation of gene topography in general and in detail. While most introns and intergenic regions reveal sequence divergences, there are, however, a number of interspersed conserved sequence motifs. In particular, two introns of the RpL14 gene contain a short, highly conserved 60 nt long sequence located at corresponding positions. This sequence represents a novel Drosophila small nucleolar RNA, which is homologous to human U49. Whereas DNA flanking the three genes shows no significant interspecies homologies, the 3'-flanking region in D. virilis contains sequences from the transposable element Penelope. The Penelope family of transposable elements has been shown to promote chromosomal rearrangements in the D. virilis species group. The presence of Penelope sequences in the D. virilis 7.3 kb genomic fragment may be indicative for a transposon-induced event of transposition which did not yet scramble the order of the three genes but led to the breakdown of sequence identity of the flanking DNA.     

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)     

Genomewide comparative analysis of the highly abundant transposable element DINE-1 suggests a recent transpositional burst in Drosophila yakuba

DINE-1 (Drosophila interspersed element) is the most abundant repetitive sequence in the Drosophila genome derived from transposable elements. It comprises >1% of the Drosophila melanogaster genome (DMG) and is believed to be a relic from an ancient transpositional burst that occurred approximately 5-10 MYA. We performed a genomewide comparison of the abundance, sequence variation, and chromosomal distribution of DINE-1 in D. melanogaster and D. yakuba. Unlike the highly diverged copies in the DMG (pairwise distance approximately 15%), DINE-1's in the Drosophila yakuba genome (DYG) have diverged by only 3.4%. Moreover, the chromosomal distribution of DINE-1 in the two species is very different, with a significant number of euchromatic insertions found only in D. yakuba. We propose that these different patterns are caused by a second transpositional burst of DINE-1's in the D. yakuba genome approximately 1.5 MYA. On the basis of the sequence of these recently transposed copies, we conclude that DINE-1 is likely to be a family of nonautomomous DNA transposons. Analysis of the chromosomal distribution of two age groups of DINE-1's in D. yakuba indicates that (1) there is a negative correlation between recombination rates and the density of DINE-1's and (2) younger copies are more evenly distributed in the chromosome arms, while older copies are mostly located near the centromere regions. Our results fit the predictions of a selection-transposition balance model. Our data on whole-genome comparison of a highly abundant TE among Drosophila sibling species demonstrate the unexpectedly dynamic nature of TE activity in different host genomes.     

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.     

P transposable element in Drosophila melanogaster: horizontal transfer]

The P transposable element family in Drosophila melanogaster is responsible for the syndrome of hybrid dysgenesis which includes chromosomal rearrangements, male recombination, high mutability and temperature sensitive agametic sterility (called gonadal dysgenesis sterility). P element activity is controlled by a complex regulation system, encoded by the elements themselves, which keeps their transposition rate low within the strain bearing P elements and limits copy number by genome. A second regulatory mechanism, which acts on the level of RNA processing, prevents P mobility to somatic cells. The oldest available strains, representing most major geographical regions of the world, exhibited no detectable hybridization to the P-element. In contrast, all recently collected natural populations that were tested carried P-element sequences. The available evidence is consistent with the hypothesis of a worldwide P-element invasion of D. melanogaster during the past 30 years. Timing and direction of the invasion are discussed. The lack of P-element in older strains of Drosophila melanogaster as well as in the species must closely related to Drosophila melanogaster, suggests that P entered the Drosophila melanogaster genome recently, probably by horizontal transfer from an other species. The analysis of P-element elsewhere in the genus Drosophila reveals that several more distantly related species carried transposable elements with sequences quite similar to P. The species with the best-matching P-element is D. willistoni. A P-element from this species was found to match all but one of the 2907 nucleotides of the Drosophila melanogaster P-element. The phylogenic distributions and the likely horizontal transfers of the two other Drosophila transposable elements are discussed.     

Mobile dispersed repeated DNA elements in the Drosophila genome

The molecular and cytogenetic organizations of 19 nonhomologous dispersed repeated sequence families were studied in 15 different laboratory strains of Drosophila melanogaster. Elements from each of the families appear to undergo transposition within the Drosophila genome, because there were striking differences in both the number and chromosomal locations of these elements between strains. A significant fraction (greater than 1%) of Drosophila DNA therefore has an unstable genomic organization. Each middle repetitive family exhibited similar variations in the chromosomal distribution of elements between the strains. Although the movements of these elements are not limited to a small number of genomic sites, there are chromosomal regions where elements from the different dispersed repeated DNA families appear to be clustered. The locations of such preferred integration sites are different in each of the D. melanogaster strains examined.     

Inferring the history of interchromosomal gene transposition in Drosophila using n-dimensional parsimony

Gene transposition puts a new gene copy in a novel genomic environment. Moreover, genes moving between the autosomes and the X chromosome experience change in several evolutionary parameters. Previous studies of gene transposition have not utilized the phylogenetic framework that becomes possible with the availability of whole genomes from multiple species. Here we used parsimonious reconstruction on the genomic distribution of gene families to analyze interchromosomal gene transposition in Drosophila. We identified 782 genes that have moved chromosomes within the phylogeny of 10 Drosophila species, including 87 gene families with multiple independent movements on different branches of the phylogeny. Using this large catalog of transposed genes, we detected accelerated sequence evolution in duplicated genes that transposed when compared to the parental copy at the original locus. We also observed a more refined picture of the biased movement of genes from the X chromosome to the autosomes. The bias of X-to-autosome movement was significantly stronger for RNA-based movements than for DNA-based movements, and among DNA-based movements there was an excess of genes moving onto the X chromosome as well. Genes involved in female-specific functions moved onto the X chromosome while genes with male-specific functions moved off the X. There was a significant overrepresentation of proteins involving chromosomal function among transposed genes, suggesting that genetic conflict between sexes and among chromosomes may be a driving force behind gene transposition in Drosophila.     

Characterizing gene movements between chromosomes in Drosophila

Genes occasionally change their location in the genome through inter-chromosomal duplication and loss. These changes happen as mistakes during recombination or through retrotransposition. In Han and Hahn 2011,(1) we surveyed the genomes of ten Drosophila species, to identify and characterize the gene transposition events in the history of these species. In the paper, we showed that the rate of gene transposition in Drosophila is higher than previously appreciated. To understand the process of gene transposition, we examined the sequences, locations, and functions of the transposed genes. Based on the elevated rate of sequence evolution in transposed genes and the frequent movements near the centromeres and telomeres, we could not reject the hypothesis that these are mutations fixed through relaxed selection. But, by examining the functions of transposed genes more carefully, we found that genes with male-specific functions and genes with female-specific functions move in opposite directions involving the X chromosome. We also found an over-representation of chromosome related functions among the transposed genes. These observations suggest the possibility of particular selection pressures contributing to gene transpositions in Drosophila.     

Molecular Evolution between Drosophila Melanogaster and D. Simulans: Reduced Codon Bias, Faster Rates of Amino Acid Substitution, and Larger Proteins in D. Melanogaster

Both natural selection and mutational biases contribute to variation in codon usage bias within Drosophila species. This study addresses the cause of codon bias differences between the sibling species, Drosophila melanogaster and D. simulans. Under a model of mutation-selection-drift, variation in mutational processes between species predicts greater base composition differences in neutrally evolving regions than in highly biased genes. Variation in selection intensity, however, predicts larger base composition differences in highly biased loci. Greater differences in the G+C content of 34 coding regions than 46 intron sequences between D. melanogaster and D. simulans suggest that D. melanogaster has undergone a reduction in selection intensity for codon bias. Computer simulations suggest at least a fivefold reduction in N(e)s at silent sites in this lineage. Other classes of molecular change show lineage effects between these species. Rates of amino acid substitution are higher in the D. melanogaster lineage than in D. simulans in 14 genes for which outgroup sequences are available. Surprisingly, protein sizes are larger in D. melanogaster than in D. simulans in the 34 genes compared between the two species. A substantial fraction of silent, replacement, and insertion/deletion mutations in coding regions may be weakly selected in Drosophila.     

The actin loci in the genus Drosophila: establishment of chromosomal homologies among five nearctic species of the Drosophila obscura group by in situ hybridization

The actin genes of five nearctic species of the Drosophila obscura group were mapped by in situ hybridization, using the 5C actin gene of D. melanogaster as a probe. In all species but D. azteca eight actin loci were observed variously dispersed over all five (A- E) chromosomal elements. In D. azteca ten actin hybridization sites were found; four of which most probably originated by duplications or by transposition events. Although the five nearctic species differ from all other Drosophila species of the D. obscura group so far studied in the number of loci as well as in the chromosomal distribution and location of the actin loci, the uniformity of the main pattern with six actin loci throughout the genus Drosophila reinforces the hypothesis that the chromosomal elements have maintained their essential identities during the course of evolution. Our findings are in accordance with the conclusion that the nearctic D. obscura species have differentiated from a common ancestor of the palearctic species and that they belong to two distinct subgroups, the pseudoobscura and the affinis subgroups.     

Molecular phylogeny and divergence times of drosophilid species

The phylogenetic relationships and divergence times of 39 drosophilid species were studied by using the coding region of the Adh gene. Four genera--Scaptodrosophila, Zaprionus, Drosophila, and Scaptomyza (from Hawaii)--and three Drosophila subgenera--Drosophila, Engiscaptomyza, and Sophophora--were included. After conducting statistical analyses of the nucleotide sequences of the Adh, Adhr (Adh-related gene), and nuclear rRNA genes and a 905-bp segment of mitochondrial DNA, we used Scaptodrosophila as the outgroup. The phylogenetic tree obtained showed that the first major division of drosophilid species occurs between subgenus Sophophora (genus Drosophila) and the group including subgenera Drosophila and Engiscaptomyza plus the genera Zaprionus and Scaptomyza. Subgenus Sophophora is then divided into D. willistoni and the clade of D. obscura and D. melanogaster species groups. In the other major drosophilid group, Zaprionus first separates from the other species, and then D. immigrans leaves the remaining group of species. This remaining group then splits into the D. repleta group and the Hawaiian drosophilid cluster (Hawaiian Drosophila, Engiscaptomyza, and Scaptomyza). Engiscaptomyza and Scaptomyza are tightly clustered. Each of the D. repleta, D. obscura, and D. melanogaster groups is monophyletic. The splitting of subgenera Drosophila and Sophophora apparently occurred about 40 Mya, whereas the D. repleta group and the Hawaiian drosophilid cluster separated about 32 Mya. By contrast, the splitting of Engiscaptomyza and Scaptomyza occurred only about 11 Mya, suggesting that Scaptomyza experienced a rapid morphological evolution. The D. obscura and D. melanogaster groups apparently diverged about 25 Mya. Many of the D. repleta group species studied here have two functional Adh genes (Adh-1 and Adh-2), and these duplicated genes can be explained by two duplication events.      

Chromosomal evolution of elements B and C in the Sophophora subgenus of Drosophila: evolutionary rate and polymorphism

The locations of 77 markers along the chromosomal elements B (41 markers) and C (36 markers) of Drosophila subobscura, D. pseudoobscura, and D. melanogaster were obtained by in situ hybridization on polytene chromosomes. In comparisons between D. subobscura and D. pseudoobscura, 10 conserved segments (accounting for 32% of the chromosomal length) were detected on element B and eight (17% of the chromosomal length) on element C. The fixation rate of paracentric inversions inferred by a maximum likelihood approach differs significantly between elements. Muller's element C (0.17 breakpoints/Mb/million years) is evolving two times faster than element B (0.08 breakpoints/Mb/million years). This difference in the evolutionary rate is paralleled by differences in the extent of chromosomal polymorphism in the corresponding lineages. Element C is highly polymorphic in D. subobscura, D. pseudoobscura, and in other obscura group species such as D. obscura and D. athabasca. In contrast, the level of polymorphism in element B is much lower in these species. The fixation rates of paracentric inversions estimated in the present study between species of the Sophophora subgenus are the highest estimates so far reported in the genus for the autosomes. At the subgenus level, there is also a parallelism between the high fixation rate and the classical observation that the species of the Sophophora subgenus tend to be more polymorphic than the species of the Drosophila subgenus. Therefore, the detected relationship between level of polymorphism and evolutionary rate might be a general characteristic of chromosomal evolution in the genus Drosophila.