Horizontal transfer and selection in the evolution of P elements

The roles of selection and horizontal transfer in the evolution of the canonical subfamily of P: elements were studied in the saltans and willistoni species groups of the genus Drosophila (subgenus Sophophora). We estimate that the common ancestor of the canonical P: subfamily dates back 2-3 Myr at the most, despite the much older age (more than 40 Myr) of the P: family as a whole. The evolution of the canonical P: subfamily is characterized by weak selection at nonsynonymous sites. These sites have evolved at three quarters the rate of synonymous sites, in which no selective constraints were detected. Their recent horizontal transfer best explains the high degree of similarity among canonical P: elements from the saltans and willistoni species groups. These results are consistent with a model of P:-element evolution in which selective constraints are imposed at the time of horizontal transfer. Furthermore, it is estimated that the spread and diversification of the canonical subfamily involved a minimum of 11 horizontal transfer events among the 18 species surveyed within the past 3 Myr. The presence of multiple P: subfamilies in the saltans and willistoni species groups is likely to be the result of multiple invasions that have previously swept through these taxa in a succession of horizontal transfer events. These results suggest that horizontal transfer among eukaryotes might be more common than anticipated.     

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.     

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.     

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.     

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.     

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.     

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.     

Vege and Mar: two novel hAT MITE families from Drosophila willistoni

Two novel families of miniature inverted repeat transposable elements (MITEs), Vege and Mar, are described from Drosophila willistoni. Based on their structures, both element families are hypothesized to belong to the hAT superfamily of transposable elements. Both elements have perfect, inverted terminal repeats and 8-bp target site duplications and were found to have inserted within fixed copies of nonautonomous P elements. Vege is present in all studied D. willistoni populations and appears to have a relatively low copy number. Mar was identified in only a single D. willistoni population, and its copy number is presently unknown. Although MITEs occupy relatively large proportions of the genomes of a broad range of organisms, this may be their first unambiguous identification in any species of the genus Drosophila.     

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.     

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.      

Sexual Interactions Between Two Distantly Related Drosophila Species, D. melanogaster and D. willistoni (Diptera: Drosophilidae)

Molecular analysis suggests that the pomace fly Drosophila melanogaster acquired the P family of transposable elements from another Drosophila species, D. willistoni. Since the two species are distantly related, it has been assumed that transmission of P element DNA from D. willistoni to D. melanogaster was mediated by a vector. The possibility of an alternative mode of transmission was assessed by characterizing the sexual behaviors of D. willistoni males and females, then observing D. willistoni and D. melanogaster males and females to see whether males from one species interacted sexually with females from the other species in a laboratory setting. We observed that D. melanogaster males court D. willistoni females vigorously and, in some cases, stimulate the females to be receptive to copulation. However, D. willistoni males perform relatively little courtship in response to D. melanogaster females and do not attempt to copulate. Thus, it is unlikely that sexual interactions effected the transmission of P element DNA from D. willistoni to D. melanogaster in the flies' natural habitat.     

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.      

New Drosophila P-like elements and reclassification of Drosophila P-elements subfamilies

Genomic searches for P-like transposable elements were performed (1) in silico in the 12 available Drosophila genomes and (2) by PCR using degenerate primers in 21 Neotropical Drosophila species. In silico searches revealed P-like sequences only in Drosophila persimilis and Drosophila willistoni. Sixteen new P-like elements were obtained by PCR. These sequences were added to sequences of previously described P-like elements, and a phylogenetic analysis was performed. The subfamilies of P-elements described in the literature (Canonical, M, O, T, and K) were included in the reconstructed tree, and all were monophyletic. However, we suggest that some subfamilies can be enlarged, other subdivided, and some new subfamilies may be proposed, totalizing eleven subfamilies, most of which contain new P-like sequences. Our analyses support the monophyly of P-like elements in Drosophilidae. We suggest that, once these elements need host-specific factors to be mobilizable, the horizontal transfer (HT) of P-like elements may be inhibited among more distant taxa. Nevertheless, HT among Drosophilidae species appears to be a common phenomenon.     

Evolutionary dynamics of wAu-like Wolbachia variants in neotropical Drosophila spp

Wolbachia bacteria are common intracellular symbionts of arthropods and have been extensively studied in Drosophila. Most research focuses on two Old Word hosts, Drosophila melanogaster and Drosophila simulans, and does not take into account that some of the Wolbachia associations in these species may have evolved only after their fast global expansion and after the exposure to Wolbachia of previously isolated habitats. Here we looked at Wolbachia of Neotropical Drosophila species. Seventy-one lines of 16 Neotropical Drosophila species sampled in different regions and at different time points were analyzed. Wolbachia is absent in lines of Drosophila willistoni collected before the 1970s, but more recent samples are infected with a strain designated wWil. Wolbachia is absent in all other species of the willistoni group. Polymorphic wWil-related strains were detected in some saltans group species, with D. septentriosaltans being coinfected with at least four variants. Based on wsp and ftsZ sequence data, wWil of D. willistoni is identical to wAu, a strain isolated from D. simulans, but can be discriminated when using a polymorphic minisatellite marker. In contrast to wAu, which infects both germ line and somatic tissues of D. simulans, wWil is found exclusively in the primordial germ line cells of D. willistoni embryos. We report on a pool of closely related Wolbachia strains in Neotropical Drosophila species as a potential source for the wAu strain in D. simulans. Possible evolutionary scenarios reconstructing the infection history of wAu-like Wolbachia in Neotropical Drosophila species and the Old World species D. simulans are discussed.     

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.     

A simulation of P element horizontal transfer in Drosophila

Experimental data suggest that the P transposable element has invaded the Drosophila melanogaster genome after a horizontal transfer from the phylogenetically distant species Drosophila willistoni. The differences between P element phylogeny and that of the Drosophila genus could in part be explained by horizontal transfers. In vivo experiments show that P elements are able to transpose in the genomes of other Drosophila species. This suggests that horizontal transmission of P elements could have taken place in many species of this genus. The regulation, transposition, and deleterious effects of the P element in D. melanogaster were formalized and integrated in a global model to produce a simulation program that simulates a P element invasion. The simulations show that our knowledge of the P element in D. melanogaster can explain its behavior in the Drosophila genus. The equilibrium state of the invaded population of a new species depends on its ability to repair damage caused by P element activity. If repair is efficient, the equilibrium state tends to be of the P type state, in which case the element could subsequently invade other populations of the species. Conversely, the equilibrium state is of the M′ type state when the ability to repair damage is low. The invasion of the P element into other populations of this new species can then only occur by genetic drift and it is likely to be lost. The success of a P element invasion into a new species thus greatly depends on its ability to produce dysgenic crosses. This revised version was published online in August 2006 with corrections to the Cover Date.     

Replicon fusions promoted by the inverted repeats of Tn5. The right repeat is an insertion sequence

Members of three repetitive sequence families were isolated from recombinant λ-genome libraries, and were used to investigate sequence relationships within these families. Studies presented elsewhere show that members of all three repeat sequence families are transcribed tissue-specifically. The thermal stability of intrafamilial heteroduplexes was measured, and the extent of colinearity between related sequences was determined by restriction mapping, heteroduplex visualization, gel blot hybridization, and direct sequencing. One large and very divergent family, named 2108, was shown to consist of an assemblage of many small repeat sequence subfamilies. Each subfamily includes <40 members which are not contiguous in the genome but are very closely related colinear sequence elements several thousand nucleotides in length. The different 2108 subfamilies share only small sequence subelements, which in each subfamily occur in a different linear order and are surrounded by different sequences. A second divergent family consisting of short repetitive sequences, the 2109 family, includes many small internally homologous subfamilies as well. A third family, 2034, displays little internal sequence divergence and no apparent subfamily structure. The repeat sequence subfamilies may be biologically significant units of repetition. Thus specific 2108 subfamilies were shown to be evolutionary conserved to a remarkable degree. Highly homologous 2108 sequences were found shared among sea urchin species which diverged almost 200 million years ago, although only about 10% of the single copy DNA sequences of these species are now homologous enough to crossreact.     

Gene conversion as a secondary mechanism of short interspersed element (SINE) evolution.

The Alu repetitive family of short interspersed elements (SINEs) in primates can be subdivided into distinct subfamilies by specific diagnostic nucleotide changes. The older subfamilies are generally very abundant, while the younger subfamilies have fewer copies. Some of the youngest Alu elements are absent in the orthologous loci of nonhuman primates, indicative of recent retroposition events, the primary mode of SINE evolution. PCR analysis of one young Alu subfamily (Sb2) member found in the low-density lipoprotein receptor gene apparently revealed the presence of this element in the green monkey, orangutan, gorilla, and chimpanzee genomes, as well as the human genome. However, sequence analysis of these genomes revealed a highly mutated, older, primate-specific Alu element was present at this position in the nonhuman primates. Comparison of the flanking DNA sequences upstream of this Alu insertion corresponded to evolution expected for standard primate phylogeny, but comparison of the Alu repeat sequences revealed that the human element departed from this phylogeny. The change in the human sequence apparently occurred by a gene conversion event only within the Alu element itself, converting it from one of the oldest to one of the youngest Alu subfamilies. Although gene conversions of Alu elements are clearly very rare, this finding shows that such events can occur and contribute to specific cases of SINE subfamily evolution.