Molecular Evolution of P Transposable Elements in the Genus Drosophila. II. The obscura Species Group

A phylogenetic analysis of P transposable elements in the Drosophila obscura species group is described. Multiple P sequences from each of 10 species were obtained using PCR primers that flank a conserved region of exon 2 of the transposase gene. In general, the P element phylogeny is congruent with the species phylogeny, indicating that the dominant mode of transmission has been vertical, from generation to generation. One manifestation of this is the distinction of P elements from the Old World obscura and subobscura subgroups from those of the New World affinis subgroup. However, the overall distribution of elements within the obscura species group is not congruent with the phylogenetic relationships of the species themselves. There are at least four distinct subfamilies of P elements, which differ in sequence from each other by as much as 34%, and some individual species carry sequences belonging to different subfamilies. P sequences from D. bifasciata are particularly interesting. These sequences belong to two subfamilies and both are distinct from all other P elements identified in this survey. Several mechanisms are postulated to be involved in determining phylogenetic relationships among P elements in the obscura group. In addition to vertical transmission, these include retention of ancestral polymorphisms and horizontal transfer by an unknown mating-independent mechanism.     

Repeated horizontal transfer of P transposons between Scaptomyza pallida and Drosophila bifasciata

Two distinct P element subfamilies, designated M-type and O-type, reside in the genome of D. bifasciata. PCR-screening of 65 Drosophila species revealed that only D. bifasciata and its closest relative D. imaii possess O-type elements. Outside the genus, O-type elements were detected in Scaptomyza pallida. Restriction analyses show that the general structure of the O-type elements from S. pallida and D. bifasciata is the same. Sequence divergence turned out to be extremely low (0.43%). These results suggest that the O-type subfamily of D. bifasciata has been received by horizontal transfer from an external source, most probably from the genus Scaptomyza, as has been previously suspected for the M-type family. Since the sequence divergence between M-type elements from S. pallida and D. bifasciata is eighteenfold higher than that between O-type elements, two independent intergeneric transfer events have to be postulated. In order to re-examine the taxonomic status of S. pallida, a partial sequence (489 bp) of the Adh gene was analysed. The data clearly prove that S. pallida has to be placed far outside the D. obscura group.     

Is the evolutionary history of the O-type P element in the saltans and willistoni groups of Drosophila similar to that of the canonical P element?

We studied the occurrence of O-type P elements in at least one species of each subgroup of the saltans group, in order to better understand the phylogenetic relationships among the elements within the saltans group and with those of species belonging to the willistoni group. We found that the O-type subfamily has a patchy distribution within the saltans group (it does not occur in D. neocordata and D. emarginata), low sequence divergence among species of the saltans group as well as in relation to species of the willistoni group, a lower rate of synonymous substitution for coding sequences compared to Adh, and phylogenetic incongruities. These findings suggest that the evolutionary history of the O-type subfamily within the saltans and willistoni groups follows the same model proposed for the canonical subfamily of P elements, i.e., events of horizontal transfer between species of the saltans and willistoni groups.     

Phylogenetic Analysis of Mos1-Like Transposable Elements in the Drosophilidae

We have performed a phylogenetic analysis of 59 mariner elements in 14 Drosophilidae species that are related to the active Drosophila mauritiana Mos1 element. This includes 38 previously described sequences and 21 new sequences amplified by PCR from 10 species. Most of the elements detected are nonfunctional due to several frameshifts and deletions. They have been subdivided into four groups according to specific signatures in the nucleotidic and amino acid sequences. The mean nucleotide diversity is 4.8 ± 0.1% and reflects mainly the divergence of inactive elements over different periods. Although this probably gives rise to occasional homoplasies between distantly related taxa, the elements of each species remain grouped together. Horizontal transfer, reported previously between D. mauritiana and Zaprionus tuberculatus, can be extended to Z. verruca, while the Mos1-like element of Z. indianus belongs to another group. Interpretation of the phylogeny leads to a comparison of the influence of common ancestral sequences and putative horizontal transfers. Received: 31 May 1999 / Accepted: 28 June 1999     

Adh Nucleotide Variation in Drosophila willistoni: High Replacement Polymorphism in an Electrophoretically Monomorphic Protein

Drosophila willistoni was the subject of intensive allozyme studies and the locus coding for alcohol dehydrogenase (Adh) was found to be virtually monomorphic. DNA sequence analysis of 18 alleles throughout the distribution of the species has revealed six replacement polymorphisms. The ratio of replacement to silent polymorphisms is higher in D. willistoni than in any other Drosophila species studied for Adh nucleotide variation. Also in contrast to other species, the variation in introns and noncoding DNA is about the same as in the coding region. We speculate that both these differences indicate D. willistoni has historically had a small population size possibly related to Pleistocene refugia in the Neotropics. Received: 5 August 1996 / Accepted: 12 April 1997     

Distribution and evolution of introns in drosophila amylase genes

While the two amylase genes of Drosophila melanogaster are intronless, the three genes of D. pseudoobscura harbor a short intron. This raises the question of the common structure of the Amy gene in Drosophila species. We have investigated the presence or absence of an intron in the amylase genes of 150 species of Drosophilids. Using polymerase chain reaction (PCR), we have amplified a region that surrounds the intron site reported in D. pseudoobscura and a few other species. The results revealed that most species contain an intron, with a variable size ranging from 50 to 750 bp, although the very majoritary size was around 60–80 bp. Several species belonging to different lineages were found to lack an intron. This loss of intervening sequence was likely due to evolutionarily independent and rather frequent events. Some other species had both types of genes: In the obscura group, and to a lesser extent in the ananassae subgroup, intronless copies had much diverged from intron-containing genes. Base composition of short introns was found to be variable and correlated with that of the surrounding exons, whereas long introns were all A-T rich. We have extended our study to non-Drosophilid insects. In species from other orders of Holometaboles, Lepidoptera and Hymenoptera, an intron was found at an identical position in the Amy gene, suggesting that the intron was ancestral. Received: 23 October 1995 / Accepted: 5 March 1996     

Characterization and Evolution of mariner Elements from Closely Related Species of Fruit Flies (Diptera: Tephritidae)

Mariner elements were amplified using the polymerase chain reaction from two species of tephritid flies, Ceratitis rosa and Trirhithrum coffeae. The sequences were ∼1.3 kb in length. None of these elements appeared to be functional, as in every case the open reading frame (ORF) was disrupted by the presence of frameshifts or stop codons. These elements, Crmar1 and Tcmar1, are very similar to the Ccmar1 element previously amplified from the closely related tephritid species C. capitata and are members of the mellifera subfamily of mariner elements. The phylogeny and pattern of divergence of these elements were examined in relation to the phylogeny of the host species. It is highly probable that the elements were present in the ancestral lineage prior to the divergence of the three species. The copy numbers of the elements within each species are very different, ranging from about 10 in T. coffeae to 5,000 in C. rosa. The possible mechanisms which determine the copy number of an element in the host genome are discussed. Received: 25 April 1997 / Accepted: 31 July 1997     

Accumulation of Transposable Elements in the Heterochromatin and on the Y Chromosome of Drosophila simulans and Drosophila melanogaster

The elements of the transposon families G, copia, mdg 1, 412, and gypsy that are located in the heterochromatin and on the Y chromosome have been identified by the Southern blotting technique in Drosophila simulans and D. melanogaster populations. Within species, the abundance of such elements differs between transposon families. Between species, the abundance in the heterochromatin and on the Y chromosome of the elements of the same family can differ greatly suggesting that differences within a species are unrelated to structural features of elements. By shedding some new light on the mechanism of accumulation of transposable elements in the heterochromatin, these data appear relevant to the understanding of the long-term interaction between transposable elements and the host genome. Received: 8 August 1997 / Accepted: 11 December 1997     

Identification of Multiple Gypsy LTR-Retrotransposon Lineages in Vertebrate Genomes

Gypsy LTR-retrotransposons have been identified in the genomes of many organisms, but only a small number of vertebrate examples have been reported to date. Here we show that members of this family are likely to be widespread in many vertebrate classes with the possible exceptions of mammals and birds. Phylogenetic analyses demonstrate that although there are several distinct lineages of vertebrate gypsy LTR-retrotransposons, the majority clusters into one monophyletic clade. Groups of fungal, plant, and insect elements were also observed, suggesting horizontal transfer between phyla may be infrequent. However, in contrast to this, there was little evidence to support sister relationships between elements derived from vertebrate and insect hosts. In fact, the majority of the vertebrate elements appeared to be most closely related to a group of gypsy LTR-retrotransposons present within fungi. This implies either that at least one horizontal transmission between these two phyla has occurred previously or that a gypsy LTR-retrotransposon lineage has been lost from insect taxa. Received: 22 December 1998 / Accepted: 6 April 1999     

Two distinct P element subfamilies in the genome of Drosophila bifasciata

The genome of Drosophila bifasciata harbours two distinct subfamilies of P-homologous sequences, designated M-type and O-type elements based on similarities to P element sequences from other species. Both subfamilies have some general features in common: they are of similar length (M-type: 2935 bp, O-type: 2986 bp), are flanked by direct repeats of 8 by (the presumptive target sequence), contain terminal inverted repeats, and have a coding region consisting of four exons. The splice sites are at homologous positions and the exons have the coding capacity for proteins of 753 amino acids (M-type) and 757 amino acids (O-type). It seems likely that both types of element represent functional transposons. The nucleotide divergence of the two P element subfamilies is high (31%). The main structural difference is observed in the terminal inverted repeats. Whereas the termini of M-type elements consist of 31 by inverted repeats, the inverted repeats of the O-type elements are interrupted by non-complementary stretches of DNA, 12 by at the 5 end and 14 by at the 3 end. This peculiarity is shared by all members of the O-type subfamily. Comparison with other P element sequences indicates incongruities between the phylogenies of the species and the P transposons. M-type and O-type elements apparently have no common origin in the D. bifasciata lineage. The M-type sequence seems to be most closely related to the P element from Scaptomyza pallida and thus could be considered as a more recent invader of the D. bifasciata gene pool. The origin of the O-type elements cannot be unequivocally deduced from the present data. The sequence comparison also provides new insights into conserved domains with possible implications for the function of P transposons.     

Two distinct P element subfamilies in the genome of Drosophila bifasciata

The genome of Drosophila bifasciata harbours two distinct subfamilies of P-homologous sequences, designated M-type and O-type elements based on similarities to P element sequences from other species. Both subfamilies have some general features in common: they are of similar length (M-type: 2935 bp, O-type: 2986 bp), are flanked by direct repeats of 8 by (the presumptive target sequence), contain terminal inverted repeats, and have a coding region consisting of four exons. The splice sites are at homologous positions and the exons have the coding capacity for proteins of 753 amino acids (M-type) and 757 amino acids (O-type). It seems likely that both types of element represent functional transposons. The nucleotide divergence of the two P element subfamilies is high (31%). The main structural difference is observed in the terminal inverted repeats. Whereas the termini of M-type elements consist of 31 by inverted repeats, the inverted repeats of the O-type elements are interrupted by non-complementary stretches of DNA, 12 by at the 5′ end and 14 by at the 3′ end. This peculiarity is shared by all members of the O-type subfamily. Comparison with other P element sequences indicates incongruities between the phylogenies of the species and the P transposons. M-type and O-type elements apparently have no common origin in the D. bifasciata lineage. The M-type sequence seems to be most closely related to the P element from Scaptomyza pallida and thus could be considered as a more recent invader of the D. bifasciata gene pool. The origin of the O-type elements cannot be unequivocally deduced from the present data. The sequence comparison also provides new insights into conserved domains with possible implications for the function of P transposons.     

Codon Usage by Transposable Elements and Their Host Genes in Five Species

We compared the codon usage of sequences of transposable elements (TEs) with that of host genes from the species Drosophila melanogaster, Arabidopsis thaliana, Caenorhabditis elegans, Saccharomyces cerevisiae, and Homo sapiens. Factorial correspondence analysis showed that, regardless of the base composition of the genome, the TEs differed from the genes of their host species by their AT-richness. In all species, the percentage of A + T on the third codon position of the TEs was higher than that on the first codon position and lower than that in the noncoding DNA of the genomes. This indicates that the codon choice is not simply the outcome of mutational bias but is also subject to selection constraints. A tendency toward higher A + T on the third position than on the first position was also found in the host genes of A. thaliana, C. elegans, and S. cerevisiae but not in those of D. melanogaster and H. sapiens. This strongly suggests that the AT choice is a host-independent characteristic common to all TEs. The codon usage of TEs generally appeared to be different from the mean of the host genes. In the AT-rich genomes of Arabidopsis thaliana, Caenorhabditis elegans, and Saccharomyces cerevisiae, the codon usage bias of TEs was similar to that of weakly expressed genes. In the GC-rich genome of D. melanogaster, however, the bias in codon usage of the TEs clearly differed from that of weakly expressed genes. These findings suggest that selection acts on TEs and that TEs may display specific behavior within the host genomes. Received: 2 May 2001 / Accepted: 29 October 2001     

Phylogenetic Position of the Subgenus Lordiphosa of the Genus Drosophila (Diptera: Drosophilidae) Inferred from Alcohol Dehydrogenase (Adh) Gene Sequences

We analyzed the phylogenetic relationship between the species of Lordiphosa and other Drosophilidae using alcohol dehydrogenase (Adh) gene sequences. The phylogenetic trees consistently show that the four species Drosophila kurokawai, D. collinella, D. stackelbergi, and D. clarofinis, which include three species groups of Lordiphosa, form a monophyletic clade. This clade is placed as a sister group to the willistoni and saltans groups of Sophophora. On the other hand, three species of Lordiphosa, D. tenuicauda, D. pseudotenuicauda, and D. acutissima, all of which belong to the tenuicauda group, are not shown to be related to the major Lordiphosa lineage. In the phylogenetic trees, these species are included into the clade comprised of Drosophila and Hirtodrosophila, although it remains uncertain whether the tenuicauda group is a monophyletic group or not. These results indicate that Lordiphosa is polyphyletic and that most of the members of the subgenus have a close relationship to the neotropical groups of Sophophora. The above conclusion is compatible with the hypothesis of Okada (Mushi [1963] 37:79–100) and Lastovka and Máca (Acta Ent Bohemoslov [1978] 75:404–420) that Lordiphosa is most closely related to Sophophora; in contrast, our results contradict the hypothesis of Grimaldi (Bull Am Mus Nat Hist [1990] 197:1–139) that Lordiphosa is a sister group to the genus Scaptomyza. Received: 12 May 1999 / Accepted: 14 April 2000     

Determination of the evolutionary relationships in Rattus sensu lato (Rodentia : Muridae) using L1 (LINE-1) amplification events

We determined ∼215 bp of DNA sequence from the 3′-untranslated region (UTR) of 240 cloned L1 (LINE-1) elements isolated from 22 species of Rattus sensu lato and Rattus sensu stricto murine rodents. The sequences were sorted into different L1 subfamilies, and oligonucleotides cognate to them were hybridized to genomic DNA of various taxa. From the distribution of the L1 subfamilies in the various species, we inferred the partial phylogeny of Rattus sensu lato. The four Maxomys species comprise a well-defined clade separate from a monophyletic cluster that contains the two Leopoldamys and four Niviventer species. The Niviventer/Leopoldamys clade, in turn, shares a node with the clade that contains Berylmys, Sundamys, Bandicota, and Rattus sensu stricto. The evolutionary relationships that we deduced agree with and significantly extend the phylogeny of Rattus sensu lato established by other molecular criteria. Furthermore, the L1 amplification events scored here produced a unique phylogenetic tree, that is, in no case did a character (a given L1 amplification event) appear on more than one branch. The lack of homoplasy found in this study supports the robustness of L1 amplification events as phylogenetic markers for the study of mammalian evolution. Received: 8 November 1996 / Accepted: 11 April 1997     

Tirant: A New Retrotransposon-Like Element in Drosophila melanogaster

In this paper we report a new retrotransposon-like element of Drosophila melanogaster called Tirant. This sequence is moderately repeated in the genome of this species and it has been found to be widely dispersed throughout its distribution area. From Southern blot and in situ analyses, this sequence appears to be mobile in D. melanogaster, since its chromosome location and the hybridization patterns vary among the different strains analyzed. In this way, partial sequencing of Tirant ends suggests that it is a retrotransposon, since it is flanked by two LTRs. The presence of sequences homologous to Tirant has been also investigated in 28 species of the genus Drosophila by means of Southern analyses. These sequences were only detected in species from melanogaster and obscura groups. These data suggest that ancestral sequences of Tirant appeared after the Sophophora radiation and before the divergence of those groups. Received: 1 January 1995 / Accepted: 20 August 1995     

EVOLUTIONARY IMPLICATIONS OF DNA DIVERGENCE IN THE DROSOPHILA OBSCURA GROUP

Using DNA–DNA hybridization, we have determined the degree of single-copy DNA (scDNA) divergence among eight species of the Drosophila obscura group. These include Old World and New World species as well as members of two subgroups. Contrary to classical systematics, members of the affinis subgroup are more closely related to American members of the obscura subgroup than are Old World species. The Old World species are not a monophyletic group. The degree of scDNA divergence among species is not necessarily correlated with morphology, chromosomal divergence, or ability to form hybrids. A unique pattern of hybrid formation was found: species separated by a ΔT_m of 6.5°C can form hybrids whereas species separated by a ΔT_m of 2.5°C cannot. As with other groups of Drosophila, the obscura group has discrete parts of the genome evolving at very different rates. The slow evolving fraction of the nuclear genome is evolving at about the same rate as mitochondrial DNA. The additional scDNA divergence accompanying the step from partial reproductive isolation (between North American pseudoobscura and the isolated Bogotà population) to full isolation is very small. The resolution of the technique was challenged by five closely related taxa with a maximum ΔT_m of 2.5°C separating them; the taxa were unambiguously resolved and the “correct” phylogeny recovered. Finally, there is some indication that scDNA in the obscura group may be evolving considerably slower than in the melanogaster subgroup.     

Horizontal transmission versus vertical inheritance of P elements in Drosophila and Scaptomyza: has the M-type subfamily spread from East Asia?

Polymerase chain reaction (PCR) screening for P elements was carried out in 77 species with a primer set highly specific for the M-type subfamily. In the course of this search M-type elements were detected in 29 species: In the melanogaster (subgroups montium and rhopaloa ) and obscura species groups of Drosophila (25 out of 71 species examined), and in the genus Scaptomyza (four out of six species). M-type elements are present in all species of the montium subgroup investigated so far (21), but occur only sporadically in other groups. Within the montium subgroup 20 species possess only incomplete copies, only one species ( D. lacteicornis ) harbours apparently full-sized elements. In contrast, outside the montium subgroup almost all species with M-type elements carry full-sized copies suggesting transpositional activity, at least in the recent past. The interior section of the full-sized M-type element of D. lacteicornis was partially sequenced (936 bp). In addition, the complete sequences of four internally deleted M-type elements of D. lacteicornis, D. rufa, D. quadraria , and D. triauraria were determined. Sequence comparisons (including sequence data from previous investigations) revealed striking discrepancies between P element phylogeny and the cladogenesis of their host species. Among several possible pathways for interspecific transfers of M-type elements, we favour the hypothesis assuming the invasion of Scaptomyza as well as the obscura group species of Drosophila via independent transmission routes originating from Asian species of the montium subgroup of Drosophila . The logical geographic scenario for these events would be East-Asia.     

Sex-specific methylation in <Emphasis Type="Italic">Drosophila</Emphasis>: an investigation of the <Emphasis Type="Italic">Sophophora</Emphasis> subgenus

Epigenetic phenomena have been widely characterized in the genomes of vertebrates and DNA methylation is a key mechanism of epigenetic regulation. The DNA methylation systems of invertebrates and vertebrates show several notable differences. However, the evolutionary implications of those differences only recently began to be revealed. Our study investigated the recurrence of sex-specific methylation, as previously described for the species Drosophila willistoni, in other species of the Sophophora subgenus that present close evolutionary relationship. The MSRE and Southern blot techniques were used to analyze rDNA of some species of the willistoni, melanogaster, saltans and obscura groups of Drosophila and the results suggested that differential DNA methylation between sexes only occurs in Drosophila tropicalis and D. insularis, two sibling species of the willistoni subgroup. However, only using the MSRE technique we could detect sex-specific patterns of DNA methylation in all species of willistoni subgroup. These results indicate that DNA methylation may present important differences, even between closely related species, shedding new light on this Neotropical species complex.     

Evolution of the Larval Cuticle Proteins Coded by the Secondary Sex Chromosome Pair: X2 and Neo-Y of Drosophila miranda: I. Comparison at the DNA Sequence Level

The larval cuticle protein genes (Lcps) represent a multigene family located at the right arm of the metacentric autosome 2 (2R) in Drosophila melanogaster. Due to a chromosome fusion the Lcp locus of Drosophila miranda is situated on a pair of secondary sex chromosomes, the X2 and neo-Y chromosome. Comparing the DNA sequences from D. miranda and D. melanogaster organization and the gene arrangement of Lcp1–Lcp4 are similar, although the intergene distances vary considerably. The greatest difference between Lcp1 and Lcp2 is due to the occurrence of a pseudogene in D. melanogaster which is not present in D. miranda. Thus the cluster of the four Lcp genes existed already before the separation of the melanogaster and obscura group. Intraspecific homogenizations of different cluster units must have occurred repeatedly between the Lcp1/Lcp2 and Lcp3/Lcp4 sequence types. The most obvious example is exon 2 of the Lcp3 gene in D. miranda, which has been substituted by the corresponding section of the Lcp4 gene rather recently. The homogenization must have occurred before the translocation which generated the neo-Y chromosome. Lcp3 of D. melanogaster has therefore no orthologous partner in D. miranda. Rearrangements in the promoter regions of the D. miranda Lcp genes have generated new, potentially functional CAAT-box motifs. Since three of the Lcp alleles on the neo-Y are not expressed and Lcp3 is expressed only at a reduced level, it is suggestive to speculate that the rearrangements might be involved as cis-regulatory elements in the up-regulation of the X2-chromosomal Lcp alleles, in Drosophila an essential process for dosage compensation. The Lcp genes on the neo-Y chromosome have accumulated more base substitutions than the corresponding alleles on the X2. Received: 27 December 1995 / Accepted: 30 April 1996