Horizontal transmission, vertical inactivation, and stochastic loss of mariner-like transposable elements

Horizontal transmission has been well documented as a major mechanism for the dissemination of mariner-like elements (MLEs) among species. Less well understood are mechanisms that limit vertical transmission of MLEs resulting in the "spotty" or discontinuous distribution observed in closely related species. In this article we present evidence that the genome of the common ancestor of the melanogaster species subgroup of Drosophila contained an MLE related to the mellifera (honey bee) subfamily. Horizontal transmission, approximately 3-10 MYA, is strongly suggested by the observation that the sequence of the MLE in Drosophila erecta is 97% identical in nucleotide sequence with that of an MLE in the cat flea, Ctenocephalides felis. The D. erecta MLE has a spotty distribution among species in the melanogaster subgroup. The element has a high copy number in D. erecta and D. orena, a moderate copy number in D. teissieri and D. yakuba, and was apparently lost ("stochastic loss") in the lineage leading to D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. In D. erecta, most copies are concentrated in the heterochromatin. Two copies from D. erecta, denoted De12 and De19, were cloned and sequenced, and they appear to be nonfunctional ("vertical inactivation"). It therefore appears that the predominant mode of MLE evolution is vertical inactivation and stochastic loss balanced against occasional reinvasion of lineages by horizontal transmission.      

Local changes in GC/AT substitution biases and in crossover frequencies on Drosophila chromosomes

I present here evidence of remarkable local changes in GC/AT substitution biases and in crossover frequencies on Drosophila chromosomes. The substitution pattern at 10 loci in the telomeric region of the X chromosome was studied for four species of the Drosophila melanogaster species subgroup. Drosophila orena and Drosophila erecta are clearly the most closely related species pair (the erecta complex) among the four species studied; however, the overall data at the 10 loci revealed a clear dichotomy in the silent substitution patterns between the AT-biased- substitution melanogaster and erecta lineages and the GC-biased-substitution yakuba and orena lineages, suggesting two or more independent changes in GC/AT substitution biases. More importantly, the results indicated a between- loci heterogeneity in GC/AT substitution bias in this small region independently in the yakuba and orena lineages. Indeed, silent substitutions in the orena lineage were significantly biased toward G and C at the consecutive yellow, lethal of scute, and asense loci, but they were significantly biased toward A and T at sta. The substitution bias toward G and C was centered in different areas in yakuba (significantly biased at EG:165H7.3, EG:171D11.2, and suppressor of sable). The similar silent substitution patterns in coding and noncoding regions, furthermore, suggested mutational biases as a cause of the substitution biases. On the other hand, previous study reveals that Drosophila yakuba has about 20-fold higher crossover frequencies in the telomeric region of the X chromosome than does D. melanogaster; this study revealed that the total genetic map length of the yakuba X chromosome was only about 1.5 times as large as that of melanogaster and that the map length of the X-telomeric y-sta region did not differ between Drosophila yakuba and D. erecta. Taken together, the data strongly suggested that an approximately 20- fold reduction in the X-telomeric crossover frequencies occurred in the ancestral population of D. melanogaster after the melanogaster-yakuba divergence but before the melanogaster-simulans divergence.     

Vertical Transmission of the Retrotransposable Elements R1 and R2 during the Evolution of the Drosophila Melanogaster Species Subgroup

R1 and R2 are non-long-terminal repeat retrotransposable elements that insert into specific sequences of insect 28S ribosomal RNA genes. These elements have been extensively described in Drosophila melanogaster. To determine whether these elements have been horizontally or vertically transmitted, we characterized R1 and R2 elements from the seven other members of the melanogaster species subgroup by genomic blotting and nucleotide sequencing. Each species was found to have homogeneous families of R1 and R2 elements with the exception of erecta and orena, which have no R2 elements. The DNA sequences of multiple R1 and R2 copies from each species indicated nucleotide divergence within each species averaged only 0.48% for R1 and 0.35% for R2, well below the level of divergence among the species. Most copies of R1 and R2 (40 of 47) sequenced from the seven species were potentially functional, as indicated by the absence of premature termination codons or translational frameshifts that would destroy the open reading frame of the element. The sequence relationships of both the R1 and R2 elements from the various members of the melanogaster subgroup closely followed that of the species phylogeny, suggesting that R1 and R2 have been stably maintained by vertical transmission since the origin of this species subgroup 17-20 million years ago. The remarkable stability of R1 and R2, compared to what has been suggested for transposable elements that insert at multiple locations in these same species, may be due to their unique specificity for sites in the rRNA gene locus. Under low copy number conditions, when it is essential for any mobile element to transpose, the insertion specificities of R1 and R2 ensure uniform developmentally regulated target sites that can be occupied with little or no detrimental effect on the host.     

Evidence for recent horizontal transfer of gypsy-homologous LTR-retrotransposon gtwin into Drosophila erecta followed by its amplification with multiple aberrations

Long terminal repeat (LTR) retrotransposon gtwin was initially discovered in silico, and then it was isolated as gypsy-homologous sequence from Drosophila melanogaster strain, G32. The presence of ORF3 suggests, that gtwin, like gypsy, may be an endogenous retrovirus, which can leave the cell and infect another one. Therefore, in this study we decided to investigate the distribution of gtwin in different species of the melanogaster subgroup in order to find out whether gtwin can be transferred horizontally as well as vertically. Gtwin was found in all 9 species of this subgroup, hence it seems to have inhabited the host genomes for a long time. In addition, we have shown that in the Drosophila erecta genome two gtwin families are present. The first one has 93% of identity to D. melanogaster element and is likely to be a descendant of gtwin that existed in Drosophila before the divergence of the melanogaster subgroup species. The other one has >99% of identity to D. melanogaster gtwin. The most reasonable explanation is that this element has been recently horizontally transferred between D. melanogaster and D. erecta. The number and variety of gtwin copies from the "infectious" family suggest that after the horizontal transfer into D. erecta genome, gtwin underwent amplification and aberrations, leading to the rise of its diverse variants.     

Reproductive Isolation and Morphogenetic Evolution in Drosophila Analyzed by Breakage of Ethological Barriers

This article reports the breaking of ethological barriers through the constitution of soma-germ line chimeras between species of the melanogaster subgroup of Drosophila, which are ethologically isolated. Female Drosophila yakuba and D. teissieri germ cells in a D. melanogaster ovary produced functional oocytes that, when fertilized by D. melanogaster sperm, gave rise to sterile yakuba-melanogaster and teissieri-melanogaster male and female hybrids. However, the erecta-melanogaster and orena-melanogaster hybrids were lethal, since female D. erecta and D. orena germ cells in a D. melanogaster ovary failed to form oocytes with the capacity to develop normally. This failure appears to be caused by an altered interaction between the melanogaster soma and the erecta and orena germ lines. Germ cells of D. teissieri and D. orena in a D. melanogaster testis produced motile sperm that was not stored in D. melanogaster females. This might be due to incompatibility between the teissieri and orena sperm and the melanogaster seminal fluid. A morphological analysis of the terminalia of yakuba-melanogaster and teissieri-melanogaster hybrids was performed. The effect on the terminalia of teissieri-melanogaster hybrids of a mutation in doublesex, a regulatory gene that controls the development of the terminalia, was also investigated.     

Local recombination and mutation effects on molecular evolution in Drosophila.

I studied the cause of the significant difference in the synonymous-substitution pattern found in the achaete-scute complex genes in two Drosophila lineages, higher codon bias in Drosophila yakuba, and lower bias in D. melanogaster. Besides these genes, the functionally unrelated yellow gene showed the same substitution pattern, suggesting a region-dependent phenomenon in the X-chromosome telomere. Because the numbers of A/T --> G/C substitutions were not significantly different from those of G/C --> A/T in the yellow noncoding regions of these species, a AT/GC mutational bias could not completely account for the synonymous-substitution biases. In contrast, we did find an approximately 14-fold difference in recombination rates in the X-chromosome telomere regions between the two species, suggesting that the reduction of recombination rates in this region resulted in the reduction of the efficacy of selection in D. melanogaster. In addition, the D. orena yellow showed a 5% increase in the G + C content at silent sites in the coding and noncoding regions since the divergence from D. erecta. This pattern was significantly different from those at the orena Adh and Amy loci. These results suggest that local changes in recombination rates and mutational pressures are contributing to the irregular synonymous-substitution patterns in Drosophila.     

Ancestral polymorphism and recent invasion of transposable elements in Drosophila species

ABSTRACT: BACKGROUND: During the evolutionary history of transposable elements, some processes, such as ancestral polymorphisms and horizontal transfer of sequences between species, can produce incongruences in phylogenies. We investigated the evolutionary history of the transposable elements Bari and 412 in the sequenced genomes of the Drosophila melanogaster group and in the sibling species D. melanogaster and D. simulans using traditional phylogenetic and network approaches. RESULTS: The maximum likelihood (ML) phylogenetic analyses revealed incongruences and unresolved relationships for both the Bari and 412 elements. The DNA transposon Bari within the D. ananassae genome is more closely related to the element of the melanogaster complex than to the sequence in D. erecta, which is inconsistent with the species phylogeny. Divergence analysis and the comparison of the rate of synonymous substitutions per synonymous site of the Bari and host gene sequences explain the incongruence as an ancestral polymorphism inherited stochastically by the derived species. Unresolved relationships were observed in the ML phylogeny of both elements involving D. melanogaster, D. simulans and D. sechellia. A network approach was used to attempt to resolve these relationships. The resulting tree suggests recent transfers of both elements between D. melanogaster and D. simulans. The divergence values of the elements between these species support this conclusion. CONCLUSIONS: We showed that an ancestral polymorphism and recent invasion of genomes due to introgression or horizontal transfer between species occurred during the evolutionary history of the Bari and 412 elements in the melanogaster group. These invasions likely occurred in Africa during the Pleistocene, before the worldwide expansion of D. melanogaster and D. simulans.     

Intra- and interspecies variation among Bari-1 elements of the melanogaster species group.

We have investigated the distribution of sequences homologous to Bari-1, a Tc1-like transposable element first identified in Drosophila melanogaster, in 87 species of the Drosophila genus. We have also isolated and sequenced Bari-1 homologues from D. simulans, D. mauritiana, and D. sechellia, the species constituting with D. melanogaster the melanogaster complex, and from D. diplacantha and D. erecta, two phylogenetically more distant species of the melanogaster group. Within the melanogaster complex the Bari-1 elements are extremely similar to each other, showing nucleotide identity values of at least 99.3%. In contrast, Bari-1-like elements from D. diplacantha and D. erecta are on average only 70% similar to D. melanogaster Bari-1 and are usually defective due to nucleotide deletions and/or insertions in the ORFs encoding their transposases. In D. erecta the defective copies are all located in the chromocenter and on chromosome 4. Surprisingly, while D. melanogaster Bari-1 elements possess 26-bp inverted terminal repeats, their D. diplacantha and D. erecta homologues possess long inverted terminal repeats similar to the terminal structures observed in the S elements of D. melanogaster and in several other Tc1-like elements of different organisms. This finding, together with the nucleotide and amino acid identity level between D. diplacantha and D. erecta elements and Bari-1 of D. melanogaster, suggests a common evolutionary origin and a rapid diversification of the termini of these Drosophila Tc1-like elements.     

A New Basal Subfamily of mariner Elements in Ceratitis rosa and Other Tephritid Flies

Several copies of highly related transposable elements, Crmar2, Almar1, and Asmar1, are described from the genomes of Ceratitis rosa, Anastrepha ludens, and A. suspensa, respectively. One copy from C. rosa, Crmar2.5, contains a full-length, uninterrupted ORF. All the other copies, from the three species contain a long deletion within the putative ORF. The consensus Crmar2 element has features typical of the mariner/Tc1 superfamily of transposable elements. In particular, the Crmar2 consensus encodes a D,D41D motif, a variant of the D,D34D catalytic domain of mariner elements. Phylogenetic analysis of the relationships of these three elements and other members of the mariner/Tc1 superfamily, based on their encoded amino acid sequences, suggests that they form a new basal subfamily of mariner elements, the rosa subfamily. BLAST analyses identified sequences from other diptera, including Drosophila melanogaster, which appear to be members of the rosa subfamily of mariner elements. Analyses of their molecular evolution suggests that Crmar2 entered the genome of C. rosa in the recent past, a consequence of horizontal transfer.     

Mitochondrial DNA evolution in themelanogaster species subgroup ofDrosophila

Detailed restriction maps (40 cleavage sites on average) of mitochondrial DNAs (mtDNAs) from the eight species of the melanogaster species subgroup of Drosophila were established. Comparison of the cleavage sites allowed us to build a phylogenetic tree based on the matrix of nucleotide distances and to select the most parsimonious network. The two methods led to similar results, which were compared with those in the literature obtained from nuclear characters. The three chromosomally homosequential species D. simulans, D. mauritiana, and D. sechellia are mitochondrially very related, but exhibit complex phylogenetic relationships. D. melanogaster is their closest relative, and the four species form a monophyletic group (the D. melanogaster complex), which is confirmed by the shared unusual length of their mt genomes (18-19 kb). The other four species of the subgroup (D. yakuba, D. teissieri, D. erecta, and D. orena) are characterized by a much shorter mt genome (16-16.5 kb). The monophyletic character of the D. yakuba complex, however, is questionable. Two species of this complex, D. yakuba and D. teissieri, are mitochondrially indistinguishable (at the level of our investigation) in spite of their noticeable allozymic and chromosomal divergence. Finally, mtDNA distances were compared with the nuclear-DNA distances thus far established. These sequences seem to evolve at rather similar rates, the mtDNA rate being barely double that of nuclear DNA.     

Changes in the recombinational environment affect divergence in the yellow gene of Drosophila

The complete coding region of the yellow (y) gene was sequenced in different Drosophila species. In the species of the melanogaster subgroup (D. melanogaster, D. simulans, D. mauritiana, D. yakuba, and D. erecta), this gene is located at the tip of the X chromosome in a region with a strong reduction in recombination rate. In contrast, in D. ananassae (included in the ananassae subgroup of the melanogaster group) and in the obscura group species (D. subobscura, D. madeirensis, D. guanche, and D. pseudoobscura), the y gene is located in regions with normal recombination rates. As predicted by the hitchhiking and background selection models, this change in the recombinational environment affected synonymous divergence in the y-gene-coding region. Estimates of the number of synonymous substitutions per site were much lower between the obscura group species and D. ananassae than between the species of the obscura group and the melanogaster subgroup. In fact, a highly significant increase in the rate of synonymous substitution was detected in all lineages leading to the species of the melanogaster subgroup relative to the D. ananassae lineage. This increase can be explained by a higher fixation rate of mutations from preferred to unpreferred codons (slightly deleterious mutations). The lower codon bias detected in all species of the melanogaster subgroup relative to D. ananassae (or to the obscura group species) would be consistent with this proposal. Therefore, at least in Drosophila, changes in the recombination rate in different lineages might cause deviations of the molecular-clock hypothesis and contribute to the overdispersion of the rate of synonymous substitution. In contrast, the change in the recombinational environment of the y gene has no detectable effect on the rate of amino acid replacement in the Yellow protein.     

Evidence for interspecific transfer of the transposable element mariner betweenDrosophila andZaprionus

Summary The transposable element mariner occurs widely in themelanogaster species group ofDrosophila. However, in drosophilids outside of themelanogaster species group, sequences showing strong DNA hybridization with mariner are found only in the genusZaprionus. the mariner sequence obtained fromZaprionus tuberculatus is 97% identical with that fromDrosophila mauritiana, a member of themelanogaster species subgroup, whereas a mariner sequence isolated fromDrosophila tsacasi is only 92% identical with that fromD. mauritiana. BecauseD. tsacasi is much more closely related toD. mauritiana than isZaprionus, the presence of mariner inZaprionus may result from horizontal transfer. In order to confirm lack of a close phylogenetic relationship between the genusZaprionus and themelanogaster species group, we compared the alcohol dehydrogenase (Adh) sequences among these species. The results show that the coding region of Adh is only 82% identical betweenZ. tuberculatus andD. mauritiana, as compared with 90% identical betweenD. tsacasi andD. mauritiana. Furthermore, the mariner gene phylogeny obtained by maximum likelihood and maximum parsimony analyses is discordant with the species phylogeny estimated by using the Adh genes. The only inconsistency in the mariner gene phylogeny is in the placement of theZaprionus mariner sequence, which clusters with mariner fromDrosophila teissieri andDrosophila yakuba in themelanogaster species subgroup. These results strongly suggest horizontal transfer.     

Molecular Evolution of the Duplicated Amy Locus in the Drosophila Melanogaster Species Subgroup: Concerted Evolution Only in the Coding Region and an Excess of Nonsynonymous Substitutions in Speciation

From the analysis of restriction maps of the Amy region in eight sibling species belonging to the Drosophila melanogaster species subgroup, we herein show that the patterns of duplication of the Amy gene are almost the same in all species. This indicates that duplication occurred before speciation within this species subgroup. From the nucleotide sequence data, we show a strong within-species similarity between the duplicated loci in the Amy coding region. This is in contrast to a strong similarity in the 5' and 3' flanking regions within each locus (proximal or distal) throughout the species subgroup. This means that concerted evolution occurred only in the Amy coding region and that differentiated evolution between the duplication occurred in the flanking regions. Moreover, when comparing the species, we also found a significant excess of nonsynonymous substitutions. In particular, all the fixed substitutions specific to D. erecta were found to be nonsynonymous. We thus conclude that adaptive protein evolution occurred in the lineage of D. erecta that is a ``specialist' species for host plants and probably also occurs in the process of speciation in general.     

Identical satellite DNA sequences in sibling species of Drosophila

The evolution of simple satellite DNAs was examined by DNA-DNA hybridization of ten Drosophila melanogaster satellite sequences to DNAs of the sibling species, Drosophila simulans and Drosophila erecta. Seven of these repeat types are present in tandem arrays in D. simulans and each of the ten sequences is repeated in D. erecta. In thermal melts, six of the seven satellite sequences in D. simulans and seven of the ten sequences in D. erecta melted within 1 deg.C of the corresponding values in D. melanogaster. The remaining sequences melted within 3 deg.C of the homologous hybrids. Therefore, there is little or no alteration in those satellite sequences held in common, despite a period of about ten million years since the divergence of D. melanogaster and D. simulans from a common ancestor. Simple satellite sequences appear to be more highly conserved than coding regions of the genome, on a per nucleotide basis. Since multiple copies of three satellite sequences could not be detected in D. simulans yet are present in D. erecta, a species more distantly related to D. melanogaster than is D. simulans, these sequences show discontinuities in evolution. There were major quantitative variations between species, showing that satellite DNAs are prone to massive amplification or diminution events over timespans as short as those separating sibling species. In D. melanogaster, these sequences amount to 21% of the genome but only 5% in D. simulans and 0.4% in D. erecta. There was a general trend of lower abundance with evolutionary distance for most satellites, suggesting that the amounts of different satellite sequences do not vary independently during evolution.     

Comparative evolutionary analysis of rDNA ITS regions in Drosophila

The internal transcribed spacer (ITS) of the ribosomal DNA is generally considered to be under low functional constraint, and it is therefore often treated as a typical nonfunctional spacer sequence. We have analyzed the ITS regions of five species from the Drosophila melanogaster subgroup, two Drosophila species from outside this group (D. pseudoobscura and D. virilis), as well as from the more distantly related dipteran fly Musca domestica. The sequence comparisons show a distinctive conservation/divergence pattern, indicating that some regions are more conserved than others. Moreover, secondary-structure calculations indicate several conserved structural elements within the ITS regions. On the other hand, a statistical test that allows us to estimate the fraction of sites that are not under selective constraint suggests that more than half of the spacer is apparently free to diverge and evolves with a rate that is close to the neutral rate of sequence evolution in Drosophila. The ITS sequences can be used to derive a molecular phylogeny for the species under study. We find that the ITS tree is largely in line with the so-far-known phylogeny of this group of species, with one difference. The species most distant within the D. melanogaster subgroup is D. yakuba, rather than D. orena, as is normally assumed.      

Structure and Evolution of the Adh Genes of Drosophila Mojavensis

The nucleotide sequence of the Adh region of Drosophila mojavensis has been completed and the region found to contain a pseudogene, Adh-2 and Adh-1 arranged in that order. Comparison of the sequence divergence of these genes to one another and to the Adh region of Drosophila mulleri and other species has allowed the development of a model for the evolution of the duplication of the Adh genes. There have been two major events. An initial duplication of an Adh gene whose dual promoter structure was similar to Drosophila melanogaster, resulted in a species with two Adh genes, one of which may have had only a proximal promoter. A second duplication of this gene generated an Adh region containing three genes. It is proposed that one of these is the ancestral gene having dual promoters, while the other two possess only proximal promoters. Subsequent events have resulted in both a change in the regulation of Adh-2 such that it is expressed as if it had a "distal" type promoter and the mutational inactivation of the most upstream gene resulting in the creation of a pseudogene. The sequence of the D. mojavensis Adh region has also revealed the presence of an element which is composed of juxtaposed inverted imperfectly repeated elements. There is a surprising and not fully explainable strong similarity of the nucleotide sequence of the 5' flanking region of the pseudogene in D. mojavensis and D. mulleri.     

Sequence Analysis of Active Mariner Elements in Natural Populations of Drosophila Simulans

Active and inactive mariner elements from natural and laboratory populations of Drosophila simulans were isolated and sequenced in order to assess their nucleotide variability and to compare them with previously isolated mariner elements from the sibling species Drosophila mauritiana and Drosophila sechellia. The active elements of D. simulans are very similar among themselves (average 99.7% nucleotide identity), suggesting that the level of mariner expression in different natural populations is largely determined by position effects, dosage effects and perhaps other factors. Furthermore, the D. simulans elements exhibit nucleotide identities of 98% or greater when compared with mariner elements from the sibling species. Parsimony analysis of mariner elements places active elements from the three species into separate groups and suggests that D. simulans is the species from which mariner elements in D. mauritiana and D. sechellia are most likely derived. This result strongly suggests that the ancestral form of mariner among these species was an active element. The two inactive mariner elements sequenced from D. simulans are very similar to the inactive peach element from D. mauritiana. The similarity may result from introgression between D. simulans and D. mauritiana or from selective constraints imposed by regulatory effects of inactive elements.