Molecular evolution of P transposable elements in the genus Drosophila. I. The saltans and willistoni species groups

A phylogenetic survey using the polymerase chain reaction (PCR) has identified four major P element subfamilies in the saltans and willistoni species groups of Drosophila. One subfamily, containing about half of the sequences studied, consists of elements that are very similar to the canonical (and active) P element from D. melanogaster. Within this subfamily, nucleotide sequence differentiation among different copies from the same species and among elements from different species is relatively low. This observation suggests that the canonical elements are relatively recent additions to the genome or, less likely, are evolving slowly relative to the other subfamilies. Elements belonging to the three noncanonical lineages are distinct from the canonical elements and from one another. Furthermore, there is considerably more sequence variation, on the average, within the noncanonical subfamilies compared to the canonical elements. Horizontal transfer and the coexistence of multiple, independently evolving element subfamilies in the same genome may explain the distribution of P elements in the saltans and willistoni species groups. Such explanations are not mutually exclusive, and each may be involved to varying degrees in the maintenance of P elements in natural populations of Drosophila.      

Evolution of P elements in natural populations of Drosophila willistoni and D. sturtevanti

To determine how population structure of the host species affects the spread of transposable elements and to assess the strength of selection acting on different structural regions, we sequenced P elements from strains of Drosophila willistoni and Drosophila sturtevanti sampled from across the distributions of these species. Elements from D. sturtevanti exhibited considerable sequence variation, and similarity among them was correlated to geographic distance between collection sites. By contrast, all D. willistoni elements sampled were essentially identical (pi < 0.2%) and exhibited patterns typical of a recent population expansion. While the canonical P elements sampled from D. sturtevanti appear to be long-time residents in that species, a rapid expansion of a very young canonical P-element lineage is suggested in D. willistoni, overcoming barriers such as large geographical distances and moderate levels of population subdivision. Between-species comparisons reveal selective constraints on P-element evolution, as indicated by significantly different substitution rates in noncoding, silent, and replacement sites. Most remarkably, in addition to replacement sites, selection pressure appears to be strong in the first and third introns and in the 3' and 5' flanking regions.     

Phylogeny of the subgenus sophophora (Diptera: drosophilidae) based on combined analysis of nuclear and mitochondrial sequences

Sequences from the nuclear (nu) alcohol dehydrogenase gene, the nu 28S ribosomal RNA locus, and the mitochondrial cytochrome oxidase II gene were used both individually and in combined analyses to infer the phylogeny of the subgenus Sophophora (Diptera: Drosophilidae). We used several optimality criteria, including maximum likelihood, maximum parsimony, and minimum evolution, to analyze these partitions to test the monophyly of the subgenus Sophophora and its four largest species groups, melanogaster, obscura, saltans, and willistoni. Our results suggest that the melanogaster and obscura species groups are each monophyletic and form a closely related clade. The Neotropical clade, containing the saltans and willistoni species groups, is also recovered, as previous studies have suggested. While the saltans species group is strongly supported as monophyletic, the results of several analyses indicate that the willistoni species group may be paraphyletic with respect to the saltans species group.     

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.     

Multiple events of horizontal transfer of the Minos transposable element between Drosophila species

In this study the Minos element was analyzed in 26 species of the repleta group and seven species of the saltans group of the genus Drosophila. The PCR and Southern blot analysis showed a wide occurrence of the Minos transposable element among species of the repleta and the saltans groups and also a low number of insertions in both genomes. Three different analyses, nucleotide divergence, historical associations, and comparisons between substitution rates (d(N) and d(S)) of Minos and Adh host gene sequences, suggest the occurrence of horizontal transfer between repleta and saltans species. These data reinforce and extend the Arca and Savakis [Genetica 108 (2000) 263] results and suggest five events of horizontal transfer to explain the present Minos distribution: between D. saltans and the ancestor of the mulleri and the mojavensis clusters; between D. hydei and the ancestor of the mulleri and the mojavensis clusters; between D. mojavensis and D. aldrichi; between D. buzzatii and D. serido; and between D. spenceri and D. emarginata. An alternative explanation would be that repeated events of horizontal transfer involving D. hydei, which is a cosmopolitan species that diverged from the others repleta species as long as 14Mya, could have spread Minos within the repleta group and to D. saltans. The data presented in this article support a model in which distribution of Minos transposon among Drosophila species is determined by horizontal transmission balanced by vertical inactivation and extinction.     

P elements in the saltans group of Drosophila: a new evaluation of their distribution and number of genomic insertion sites

Few are studies on P elements that have addressed the saltans group. These studies had shown that species from the cordata and elliptica subgroups were devoid of any discernible P homologous sequences, while species from the parasaltans, sturtevanti, and saltans subgroups all contain P element sequences. Our analyses showed the presence of one to 15 P element insertion sites in species of the saltans group, including Drosophila neocordata and Drosophila emarginata (cordata and elliptica subgroups, respectively). From these species, only those from the parasaltans, sturtevanti, and saltans subgroups harbor canonical P elements and, only those of the last two subgroups seem to harbor putative full-sized elements. Due to the low similarity of the sequences found in D. neocordata and D. emarginata to those earlier described, we suggest that these sequences might be rudimental P element derivatives that were present in the ancestral of the subgenus Sophophora.     

Molecular evolution of P transposable elements in the genus Drosophila. III. The melanogaster species group

Phylogenetic relationships were determined for 76 partial P-element sequences from 14 species of the melanogaster species group within the Drosophila subgenus Sophophora. These results are examined in the context of the phylogeny of the species from which the sequences were isolated. Sequences from the P-element family fall into distinct subfamilies, or clades, which are often characteristic for particular species subgroups. When examined locally among closely related species, the evolution of P elements is characterized by vertical transmission, whereby the P-element phylogeny traces the species phylogeny. On a broader scale, however, the P-element phylogeny is not congruent with the species phylogeny. One feature of P-element evolution in the melanogaster group is the presence of more than one P-element subfamily, differing by as much as 36%, in the genomes of some species. Thus, P elements from several individual species are not monophyletic, and a likely explanation for the incongruence between P-element and species phylogenies is provided by the comparison of paralogous sequences. In certain instances, horizontal transfer seems to be a valid alternative explanation for lack of congruence between species and P-element phylogenies. The canonical P-element subfamily, which represents the active, autonomous transposable element, is restricted to D. melanogaster. Thus, its origin clearly lies outside of the melanogaster species group, consistent with the earlier conclusion of recent horizontal transfer.      

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.     

Nonrecombining genes in a recombination environment: the Drosophila "dot" chromosome

Rate of recombination is a powerful variable affecting several aspects of molecular variation and evolution. A nonrecombining portion of the genome of most Drosophila species, the "dot" chromosome or F element, exhibits very low levels of variation and unusual codon usage. One lineage of Drosophila, the willistoni/saltans groups, has the F element fused to a normally recombining E element. Here, we present polymorphism data for genes on the F element in two Drosophila willistoni and one D. insularis populations, genes previously studied in D. melanogaster. The D. willistoni populations were known to be very low in inversion polymorphism, thus minimizing the recombination suppression effect of inversions. We first confirmed, by in situ hybridization, that D. insularis has the same E + F fusion as D. willistoni, implying this was a monophyletic event. A clear gradient in codon usage exists along the willistoni F element, from the centromere distally to the fusion with E; estimates of recombination rates parallel this gradient and also indicate D. insularis has greater recombination than D. willistoni. In contrast to D. melanogaster, genes on the F element exhibit moderate levels of nucleotide polymorphism not distinguishable from two genes elsewhere in the genome. Although some linkage disequilibrium (LD) was detected between polymorphic sites within genes (generally <500 bp apart), no long-range LD between F element loci exists in the two willistoni group species. In general, the distribution of allele frequencies of F element genes display the typical pattern of expectations of neutral variation at equilibrium. These results are consistent with the hypothesis that recombination allows the accumulation of nucleotide variation as well as allows selection to act on synonymous codon usage. It is estimated that the fusion occurred ～20 Mya and while the F element in the willistoni lineage has evolved "normal" levels and patterns of nucleotide variation, equilibrium may not have been reached for codon usage.     

Drosophila mediopunctata P elements: a new example of horizontal transfer.

Sequences homologous to the P element of Drosophila melanogaster were previously identified in Drosophila mediopunctata, a member of the tripunctata group, subgenus Drosophila. We report here that the P element is present in about three to five copies in the D. mediopunctata genome. While one of the insertion sites appears to be fixed, others may be polymorphic, indicating relatively recent P element activity. Phylogenetic analysis revealed that the D. mediopunctata element belongs to the canonical subfamily of P elements and that divergence of the D. mediopunctata element from other members of this subfamily ranges from 2% to 5% at the nucleotide level. This is the first report of a canonical P element outside the subgenus Sophophora. Based primarily on the striking incongruence between P element and host species phylogenies, the presence of a canonical P element in D. mediopunctata is most likely explained by horizontal transfer between species.     

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.     

Recent horizontal transfer of a mariner transposable element among and between Diptera and Neuroptera

Transposable elements of the mariner family are widespread among insects and other invertebrates, and initial analyses of their relationships indicated frequent occurrence of horizontal transfers between hosts. A specific PCR assay was used to screen for additional members of the irritans subfamily of mariners in more than 400 arthropod species. Phylogenetic analysis of cloned PCR fragments indicated that relatively recent horizontal transfers had occurred into the lineages of a fruit fly Drosophila ananassae, the horn fly Haematobia irritans, the African malaria vector mosquito Anopheles gambiae, and a green lacewing Chrysoperla plorabunda. Genomic dot-blot analysis revealed that the copy number in these species varies widely, from about 17,000 copies in the horn fly to three copies in D. ananassae. Multiple copies were sequenced from genomic clones from each of these species and four others with related elements. These sequences confirmed the PCR results, revealing extremely similar elements in each of these four species (greater than 88% DNA and 95% amino acid identity). In particular, the consensus sequence of the transposase gene of the horn fly elements differs by just two base pairs out of 1,044 from that of the lacewing elements. The mosquito lineage has diverged from the other Diptera for over 200 Myr, and the neuropteran last shared a common ancestor with them more than 265 Myr ago, so this high similarity implies that these transposons recently transferred horizontally into each lineage. Their presence in only the closest relatives in at least the lacewing lineage supports this hypothesis. Such horizontal transfers provide an explanation for the evolutionary persistence and widespread distribution of mariner transposons. We propose that the ability to transfer horizontally to new hosts before extinction by mutation in the current host constitutes the primary selective constraint maintaining the sequence conservation of mariners and perhaps other DNA-mediated elements.      

Fluctuating Mutation Bias and the Evolution of Base Composition in Drosophila

The idea that the pattern of point mutation in Drosophila has remained constant during the evolution of the genus has recently been challenged. A study of the nucleotide composition focused on the Drosophila saltans group has evidenced unsuspected nucleotide composition differences among lineages. Compositional differences are associated with an accelerated rate of amino acid replacement in functionally less constrained regions. Here we reassess this issue from a different perspective. Adopting a maximum-likelihood estimation approach, we focus on the different predictions that mutation and selection make about the nonsynonymous-to-synonymous rate ratio. We investigate two gene regions, alcohol dehydrogenase (Adh) and xanthine dehydrogenase (Xdh), using a balanced data set that comprises representatives from the melangaster, obscura, saltans, and willistoni groups. We also consider representatives of the Hawaiian picture-winged group. These Hawaiian species are known to have experienced repeated bottlenecks and are included as a reference for comparison. Our results confirm patterns previously detected. The branch ancestral to the fast-evolving willistoni/saltans lineage, where most of the change in GC content has occurred, exhibits an excess of synonymous substitutions. The shift in mutation bias has affected the extent of the rate variation among sites in Xdh. Received: 4 May 1999 / Accepted: 26 July 1999     

Selective constraints on P-element evolution.

P elements, like mariners, inhabit eukaryotic genomes and transpose via a DNA intermediate. Mutant and wild-type elements in the same genome should be transposed with equal probability by trans-acting transposase, and so no selection should counteract the accumulation of inactivating mutations in transposase genes. Thus, copies of mariner elements diverge within a host species under no selection (Robertson and Lampe 1995). It is unknown whether or not this pattern holds for P elements, which are unrelated to mariner elements but share the same life history. Publicly available P-element sequences were analyzed for evidence of conservative selection for the function of P-element-encoded proteins. Results were compared to predictions derived from several hypotheses that could explain selection, or the lack of it. P-element protein-coding sequences do evolve under conservative selection but apparently because of more than one selective force. Of the four exons in the P-element transposase, the first three (exons 0, 1, and 2) can be translated alone into a repressor of transposition, while the last (exon 3) is only expressed as part of the full-length transposase and probably serves a transposition-specific role. As full-length P-element copies diverge from each other within a host population, selection maintains exons 0-2 but apparently not exon 3. The selection acting on exons 0-2 may act at the host level for repression of transposition (since host level selection does act on orthologous truncated elements that contain only exons 0-2). Evidence of selection on exon 3 is only found in comparisons of more diverged elements from different species, suggesting that selection for transposition acts primarily at horizontal transfer events. Thus, horizontal transfer events may be the sole source of the selection that is crucial to the maintenance of autonomous P elements in the face of mutation (as suggested by Robertson and Lampe 1995). The predictions derived here suggest a strategy for collecting sequence data that could definitively answer these questions.     

An unusual evolutionary behaviour of a sea urchin histone gene cluster

DNA sequences of cloned histone coding sequences and spacers of sea urchin species that diverged long ago in evolution were compared. The highly repeated H4 and H3 genes active during early embryogenesis had evolved (in their silent sites) at a rate (0.5-0.6% base changes/Myr) similar to single-copy protein-coding genes and nearly as fast as spacer DNA (0.7% base changes/Myr) and unique DNA. Thus, evolution in the major histone genes conforms to a universal evolutionary clock based on the rate of base sequence change. By contrast, the H4 and H3 coding sequences and a non-transcribed spacer of the DNA clone h19 of Psammechinus miliaris show an exceptionally low rate of sequence evolution only 1/100 to 1/200 that predicted from the clock hypothesis. According to the classical model of gene inheritance, the h19 DNA sequences in the Psammechinus genome require unusual conservation mechanisms by selection at the level of the gene and spacer sequences. An alternative explanation could be recent horizontal gene transfer of a histone gene cluster from the very distantly related Strongylocentrotus dröbachiensis to the P. miliaris genome.     

Widespread occurence of mariner transposons in coastal crabs

Mariner-like elements (MLEs) are ubiquitous DNA mobile elements found in almost all eukaryote genomes. Nevertheless most of the known copies are inactive and the question of the genome invasion by MLEs remains largely hypothetical. We have previously reported the presence of highly homologous copies of MLEs in the genome of phylogenetically distant crustacea living in the same hydrothermal environment suggesting the possibility of horizontal transfer. In order to further support the hypothesis that horizontal transmission of MLEs might occur between crustacean sympatric species, we described here 85 MLE sequences found in the genome of a large spectrum of coastal crab species. The number of the MLEs copies in genomes was variable. Half of these MLEs fit with the irritans subfamily of MLEs whereas the second half grouped in a new subfamily called marmoratus. In addition, a molecular phylogeny of crabs was established by using the 16S information. The comparison between 16S and MLEs based trees reveals their incongruence, and suggests either the existence of horizontal transfer events between phylogenetically distant species, or an ancestral MLE polymorphism followed by different evolution and stochastic loss.