Multiple lineages of R1 retrotransposable elements can coexist in the rDNA loci of Drosophila

Non-LTR R2 element evolutionary patterns: phylogenetic incongruences, rapid radiation and the maintenance of multiple lineages

Retrotransposons of the R2 superclade specifically insert within the 28S ribosomal gene. They have been isolated from a variety of metazoan genomes and were found vertically inherited even if their phylogeny does not always agree with that of the host species. This was explained with the diversification/extinction of paralogous lineages, being proved the absence of horizontal transfer. We here analyze the widest available collection of R2 sequences, either newly isolated from recently sequenced genomes or drawn from public databases, in a phylogenetic framework. Results are congruent with previous analyses, but new important issues emerge. First, the N-terminal end of the R2-B clade protein, so far unknown, presents a new zinc fingers configuration. Second, the phylogenetic pattern is consistent with an ancient, rapid radiation of R2 lineages: being the estimated time of R2 origin (850–600 Million years ago) placed just before the metazoan Cambrian explosion, the wide element diversity and the incongruence with the host phylogeny could be attributable to the sudden expansion of available niches represented by host’s 28S ribosomal genes. Finally, we detect instances of coexisting multiple R2 lineages showing a non-random phylogenetic pattern, strongly similar to that of the “library” model known for tandem repeats: a collection of R2s were present in the ancestral genome and then differentially activated/repressed in the derived species. Models for activation/repression as well as mechanisms for sequence maintenance are also discussed within this framework.

A population genetic model for the maintenance of R2 retrotransposons in rRNA gene loci

R2 retrotransposable elements exclusively insert into the tandemly repeated rRNA genes, the rDNA loci, of their animal hosts. R2 elements form stable long-term associations with their host, in which all individuals in a population contain many potentially active copies, but only a fraction of these individuals show active R2 retrotransposition. Previous studies have found that R2 RNA transcripts are processed from a 28S co-transcript and that the likelihood of R2-inserted units being transcribed is dependent upon their distribution within the rDNA locus. Here we analyze the rDNA locus and R2 elements from nearly 100 R2-active and R2-inactive individuals from natural populations of Drosophila simulans. Along with previous findings concerning the structure and expression of the rDNA loci, these data were incorporated into computer simulations to model the crossover events that give rise to the concerted evolution of the rRNA genes. The simulations that best reproduce the population data assume that only about 40 rDNA units out of the over 200 total units are actively transcribed and that these transcribed units are clustered in a single region of the locus. In the model, the host establishes this transcription domain at each generation in the region with the fewest R2 insertions. Only if the host cannot avoid R2 insertions within this 40-unit domain are R2 elements active in that generation. The simulations also require that most crossover events in the locus occur in the transcription domain in order to explain the empirical observation that R2 elements are seldom duplicated by crossover events. Thus the key to the long-term stability of R2 elements is the stochastic nature of the crossover events within the rDNA locus, and the inevitable expansions and contractions that introduce and remove R2-inserted units from the transcriptionally active domain.

Selfish transposable elements survive in eukaryotic genomes despite the elaborate mechanisms developed by the hosts to limit their activity. One accessible system that simplifies the complex interactions between element and host involves the R2 elements, which exclusively insert in the tandemly arranged rRNA genes. R2 exhibits remarkable stability in animal lineages even though each insertion inactivates one rRNA gene. Here we determine the size of the rDNA locus and R2 number in natural isolates of Drosophila simulans. Combined with previous data concerning the expression and regulation of R2, we develop a detailed population genetic model for rRNA gene and R2 evolution that duplicates all properties of the rRNA loci in natural populations. Critical components of the model are that only a contiguous 40 unit array of rRNA gene units are needed for transcription, that R2 elements are active only when present in this transcription domain, and that most of the crossovers in the rDNA loci occur in this domain. These results suggest that the key to the long-term survival of R2 is the redistribution of rDNA units in the locus brought about by the crossovers that maintain sequence identity in all rDNA units.

R2 and R2/R1 hybrid non-autonomous retrotransposons derived by internal deletions of full-length elements

BACKGROUND

R2 is a non-long terminal repeat (non-LTR) retrotransposable element that inserts site specifically into the 28S genes of the ribosomal (r)RNA gene loci. Encoded at the 5' end is a ribozyme that generates the precise 5' end by self-cleavage of a 28S gene cotranscript. Sequences at the 3' end are necessary for the R2 protein to bind RNA and initiate the target primed reverse transcription (TPRT) reaction. These minimal RNA requirements suggested that if recombination/DNA repair conjoined the 5' and 3' ends of R2, the result would be a non-autonomous element that could survive as long as autonomous R2 elements supplied the TPRT activity.

RESULTS

A PCR-based survey of 39 Drosophila species aided by genomic sequences from 12 of these species revealed two types of non-autonomous elements. We call these elements SIDEs (for 'Short Internally Deleted Elements'). The first consisted of a 5' ribozyme and a 3' end of an R2 element as predicted. Variation at the 5' junctions of the R2 SIDE copies was typical for R2 insertions suggesting their propagation by TPRT. The second class of SIDE contained sequences from R1 elements, another non-LTR retrotransposon that inserts into rRNA gene loci. These insertions had an R2 ribozyme immediately upstream of R1 3' end sequences. These hybrid SIDEs were inserted at the R1 site with 14 bp target site duplications typical of R1 insertions suggesting they used the R1 machinery for retrotransposition. Finally, the survey revealed examples of U12 small nuclear (sn)RNA and tRNA sequences at the 5' end of R2 elements suggesting the R2 reverse transcriptase can template jump from the R2 transcript to a second RNA during TPRT.

CONCLUSIONS

The R2 SIDE and R2/R1 hybrid SIDEs are rare examples of non-autonomous retrotransposons in the Drosophila genome. Associated non-autonomous elements and in vivo template jumps are two additional characteristics R2 shares with other non-LTR retrotransposons such as mammalian L1s. Analysis of the hybrid SIDEs provides supporting evidence that R1 elements, like R2 elements, recognize their 3' untranslated region (UTR) sequences and, thus, belong to the stringent class of non-LTR elements.

Wolbachia-mediated persistence of mtDNA from a potentially extinct species

Drosophila quinaria is polymorphic for infection with Wolbachia, a maternally transmitted endosymbiont. Wolbachia-infected individuals carry mtDNA that is only distantly related to the mtDNA of uninfected individuals, and the clade encompassing all mtDNA haplotypes within D. quinaria also includes the mtDNA of several other species of Drosophila. Nuclear gene variation reveals no difference between the Wolbachia-infected and uninfected individuals of D. quinaria, indicating that they all belong to the same interbreeding biological species. We suggest that the Wolbachia and the mtDNA with which it is associated were derived via interspecific hybridization and introgression. The sequences in the Wolbachia and the associated mtDNA are ≥6% divergent from those of any known Drosophila species. Thus, in spite of nearly complete species sampling, the sequences from which these mitochondria were derived remain unknown, raising the possibility that the donor species is extinct. The association between Wolbachia infection and mtDNA type within D. quinaria suggests that Wolbachia may be required for the continued persistence of the mtDNA from an otherwise extinct Drosophila species. We hypothesize that pathogen-protective effects conferred by Wolbachia operate in a negative frequency-dependent manner, thus bringing about a stable polymorphism for Wolbachia infection.

Histone and ribosomal RNA repetitive gene clusters of the boll weevil are linked in a tandem array

Histones are the major protein component of chromatin structure. The histone family is made up of a quintet of proteins, four core histones (H2A, H2B, H3 & H4) and the linker histones (H1). Spacers are found between the coding regions. Among insects this quintet of genes is usually clustered and the clusters are tandemly repeated. Ribosomal DNA contains a cluster of the rRNA sequences 18S, 5.8S and 28S. The rRNA genes are separated by the spacers ITS1, ITS2 and IGS. This cluster is also tandemly repeated. We found that the ribosomal RNA repeat unit of at least two species of Anthonomine weevils, Anthonomus grandis and Anthonomus texanus (Coleoptera: Curculionidae), is interspersed with a block containing the histone gene quintet. The histone genes are situated between the rRNA 18S and 28S genes in what is known as the intergenic spacer region (IGS). The complete reiterated Anthonomus grandis histone-ribosomal sequence is 16,248 bp.

Maintenance of multiple lineages of R1 and R2 retrotransposable elements in the ribosomal RNA gene loci of Nasonia

Sequencing reads from the Nasonia genome project were used to study the ribosomal RNA gene loci and the retrotransposons R1 and R2 that insert specifically into the 28S genes. Five highly divergent R1 and five highly divergent R2 families were identified in the three sequenced species, as well as a non-autonomous element that appears to use the retrotransposition machinery of R1. A duplication of the R1 target site within the spacer region of the rDNA units was also found to be extensively utilized by R1 elements. We document numerous instances where the R1 and R2 families appropriated parts of the retrotransposition machinery of other lineages and speculate that this enables rapid adaptation and the maintenance of multiple R1 and R2 families.

Origin of nascent lineages and the mechanisms used to prime second-strand DNA synthesis in the R1 and R2 retrotransposons of Drosophila

Comparative analysis of 12 Drosophila genomes reveals insights into the evolution and mechanism of integration of R1 and R2 retrotransposons.

Background

Most arthropods contain R1 and R2 retrotransposons that specifically insert into the 28S rRNA genes. Here, the sequencing reads from 12 Drosophila genomes have been used to address two questions concerning these elements. First, to what extent is the evolution of these elements subject to the concerted evolution process that is responsible for sequence homogeneity among the different copies of rRNA genes? Second, how precise are the target DNA cleavages and priming of DNA synthesis used by these elements?

Results

Most copies of R1 and R2 in each species were found to exhibit less than 0.2% sequence divergence. However, in many species evidence was obtained for the formation of distinct sublineages of elements, particularly in the case of R1. Analysis of the hundreds of R1 and R2 junctions with the 28S gene revealed that cleavage of the first DNA strand was precise both in location and the priming of reverse transcription. Cleavage of the second DNA strand was less precise within a species, differed between species, and gave rise to variable priming mechanisms for second strand synthesis.

Conclusions

These findings suggest that the high sequence identity amongst R1 and R2 copies is because all copies are relatively new. However, each active element generates its own independent lineage that can eventually populate the locus. Independent lineages occur more often with R1, possibly because these elements contain their own promoter. Finally, both R1 and R2 use imprecise, rapidly evolving mechanisms to cleave the second strand and prime second strand synthesis.

Role of recombination in the long-term retention of transposable elements in rRNA gene loci

Multiple theoretical studies have focused on the concerted evolution of the tandemly repeated rRNA genes of eukaryotes; however, these studies did not consider the transposable elements that interrupt the rRNA genes in many organisms. For example, in insects, R1 and R2 have been stable components of the rDNA locus for hundreds of millions of years, suggesting either that they have minimal effects on fitness or that they are unable to be eliminated. We constructed a simulation model of recombination and retrotransposition within the rDNA locus that addresses the population dynamics and fitness consequences associated with R1 and R2 insertions. The simulations suggest that even without R1 and R2 retrotransposition the frequent sister chromatid exchanges postulated from various empirical studies will, in combination with selection, generate rDNA loci that are much larger than those needed for transcription. These large loci enable the host to tolerate high levels of R1 and R2 insertions with little fitness consequences. Changes in retrotransposition rates are likely to be accommodated by adjustments in sister chromatid exchange (SCE) rate, rather than by direct selection on the number of uninserted rDNA units. These simulations suggest that the rDNA locus serves as an ideal niche for the long-term survival of transposable elements.

The Juan non-LTR retrotransposon in mosquitoes: genomic impact, vertical transmission and indications of recent and widespread activity

BACKGROUND

In contrast to DNA-mediated transposable elements (TEs), retrotransposons, particularly non-long terminal repeat retrotransposons (non-LTRs), are generally considered to have a much lower propensity towards horizontal transfer. Detailed studies on site-specific non-LTR families have demonstrated strict vertical transmission. More studies are needed with non-site-specific non-LTR families to determine whether strict vertical transmission is a phenomenon related to site specificity or a more general characteristic of all non-LTRs. Juan is a Jockey clade non-LTR retrotransposon first discovered in mosquitoes that is widely distributed in the mosquito family Culicidae. Being a non-site specific non-LTR, Juan offers an opportunity to further investigate the hypothesis that non-LTRs are genomic elements that are primarily vertically transmitted.

RESULTS

Systematic analysis of the ~1.3 Gbp Aedes aegypti (Ae. aegypti) genome sequence suggests that Juan-A is the only Juan-type non-LTR in Aedes aegypti. Juan-A is highly reiterated and comprises approximately 3% of the genome. Using minimum cutoffs of 90% length and 70% nucleotide (nt) identity, 663 copies were found by BLAST using the published Juan-A sequence as the query. All 663 copies are at least 95% identical to Juan-A, while 378 of these copies are 99% identical to Juan-A, indicating that the Juan-A family has been transposing recently in evolutionary history. Using the 0.34 Kb 5' UTR as the query, over 2000 copies were identified that may contain internal promoters, leading to questions on the genomic impact of Juan-A. Juan sequences were obtained by PCR, library screening, and database searches for 18 mosquito species of six genera including Aedes, Ochlerotatus, Psorophora, Culex, Deinocerites, and Wyeomyia. Comparison of host and Juan phylogenies shows overall congruence with few exceptions.

CONCLUSION

Juan-A is a major genomic component in Ae. aegypti and it has been retrotransposing recently in evolutionary history. There are also indications that Juan has been recently active in a wide range of mosquito species. Furthermore, our research demonstrates that a Jockey clade non-LTR without target site-specificity has been sustained by vertical transmission in the mosquito family. These results strengthen the argument that non-LTRs tend to be genomic elements capable of persistence by vertical descent over a long evolutionary time.

Finely orchestrated movements: evolution of the ribosomal RNA genes

Evolution of the tandemly repeated ribosomal RNA (rRNA) genes is intriguing because in each species all units within the array are highly uniform in sequence but that sequence differs between species. In this review we summarize the origins of the current models to explain this process of concerted evolution, emphasizing early studies of recombination in yeast and more recent studies in Drosophila and mammalian systems. These studies suggest that unequal crossover is the major driving force in the evolution of the rRNA genes with sister chromatid exchange occurring more often than exchange between homologs. Gene conversion is also believed to play a role; however, direct evidence for its involvement has not been obtained. Remarkably, concerted evolution is so well orchestrated that even transposable elements that insert into a large fraction of the rRNA genes appear to have little effect on the process. Finally, we summarize data that suggest that recombination in the rDNA locus of higher eukaryotes is sufficiently frequent to monitor changes within a few generations.

Chromatin structure and transcription of the R1- and R2-inserted rRNA genes of Drosophila melanogaster

About half of the rRNA gene units (rDNA units) of Drosophila melanogaster are inserted by the retrotransposable elements R1 and R2. Because transcripts to R1 and R2 were difficult to detect on blots and electron microscopic observations of rRNA synthesis suggested that only uninserted rDNA units were transcribed, it has long been postulated that inserted rDNA units are in a repressed (inactive) chromatin structure. Studies described here suggest that inserted and uninserted units are equally accessible to DNase I and micrococcal nuclease and contain similar levels of histone H3 and H4 acetylation and H3K9 methylation. These studies have low sensitivity, because psoralen cross-linking suggested few (estimated <10%) of the rDNA units of any type are transcriptionally active. Nuclear run-on experiments revealed that R1-inserted and R2-inserted units are activated for transcription at about 1/5 and 1/10, respectively, the rate of uninserted units. Most transcription complexes of the inserted units terminate within the elements, thus explaining why previous molecular and electron microscopic methods indicated inserted units are seldom transcribed. The accumulating data suggest that all units within small regions of the rDNA loci are activated for transcription, with most control over R1 and R2 activity involving steps downstream of transcription initiation.

The mode of evolution of aggregation pheromones in Drosophila species

Aggregation pheromones are used by fruit flies of the genus Drosophila to assemble on breeding substrates, where they feed, mate and oviposit communally. These pheromones consist of species-specific blends of chemicals. Here, using a phylogenetic framework, we examine how differences among species in these pheromone blends have evolved. Theoretical predictions, genetic evidence, and previous empirical analysis of bark beetle species, suggest that aggregation pheromones do not evolve gradually, but via major, saltational shifts in chemical composition. Using pheromone data for 28 species of Drosophila we show that, unlike with bark beetles, the distribution of chemical components among species is highly congruent with their phylogeny, with closely related species being more similar in their pheromone blends than are distantly related species. This pattern is also strong within the melanogaster species group, but less so within the virilis species group. Our analysis strongly suggests that the aggregation pheromones of Drosophila exhibit a gradual, not saltational, mode of evolution. We propose that these findings reflect the function of the pheromones in the ecology of Drosophila, which does not hinge on species specificity of aggregation pheromones as signals.

Competition between R1 and R2 transposable elements in the 28S rRNA genes of insects

R1 and R2 are non-LTR retrotransposons that insert in the 28S rRNA genes of arthropods. R1 elements insert into a site that is 74 bp downstream of the R2 insertion site, thus the presence of an R2 in the same 28S gene may inhibit the expression of R1. Consistent with such a suggestion, the R1 elements of Drosophila melanogaster have a strong bias against inserting into 28S genes already containing an R2 element. R2 elements, on the other hand, are only 2-3 fold inhibited from inserting into a 28S gene already containing an R1. D. melanogaster R1 elements are unusual in that they generate a 23-bp deletion of the target site upstream of the insertion. Using in vitro assays developed to study R2 integration, we show that the presence of R1 sequences 51 bp downstream of the R2 insertion site changes the nucleosomal structure that can be formed by the R2 target site. The R2 endonuclease is inhibited from cleaving these altered nucleosomes. We suggest that R1 elements have been selected to make this large deletion of the 28S gene to block the insertion of an upstream R2 element. These findings are consistent with the model that R1 and R2 are in competition for the limited number of insertion sites available within their host's genome.

APE-type non-LTR retrotransposons: determinants involved in target site recognition

Non-long terminal repeat (Non-LTR) retrotransposons represent a diverse and widely distributed group of transposable elements and an almost ubiquitous component of eukaryotic genomes that has a major impact on evolution. Their copy number can range from a few to several million and they often make up a significant fraction of the genomes. The members of the dominating subtype of non-LTR retrotransposons code for an endonuclease with homology to apurinic/apyrimidinic endonucleases (APE), and are thus termed APE-type non-LTR retrotransposons. In the last decade both the number of identified non-LTR retrotransposons and our knowledge of biology and evolution of APE-type non-LTR retrotransposons has increased tremendously.

Monitoring the mode and tempo of concerted evolution in the Drosophila melanogaster rDNA locus

Non-LTR retrotransposons R1 and R2 have persisted in rRNA gene loci (rDNA) since the origin of arthropods despite their continued elimination by the recombinational mechanisms of concerted evolution. This study evaluated the short-term evolutionary dynamics of the rDNA locus by measuring the divergence among replicate Drosophila melanogaster lines after 400 generations. The total number of rDNA units on the X chromosome of each line varied from 140 to 310, while the fraction of units inserted with R1 and R2 retrotransposons ranged from 37 to 65%. This level of variation is comparable to that found in natural population surveys. Variation in locus size and retrotransposon load was correlated with large changes in the number of uninserted and R1-inserted units, yet the numbers of R2-inserted units were relatively unchanged. Intergenic spacer (IGS) region length variants were also used to evaluate changes in the rDNA loci. All IGS length variants present in the lines showed significant increases and decreases of copy number. These studies, combined with previous data following specific R1 and R2 insertions in these lines, help to define the type and distribution, both within the locus and within the individual units, of recombinational events that give rise to the concerted evolution of the rDNA locus.

Decreased diversity but increased substitution rate in host mtDNA as a consequence of Wolbachia endosymbiont infection

A substantial fraction of insects and other terrestrial arthropods are infected with parasitic, maternally transmitted endosymbiotic bacteria that manipulate host reproduction. In addition to imposing direct selection on the host to resist these effects, endosymbionts may also have indirect effects on the evolution of the mtDNA with which they are cotransmitted. Patterns of mtDNA diversity and evolution were examined in Drosophila recens, which is infected with the endosymbiont Wolbachia, and its uninfected sister species D. subquinaria. The level of mitochondrial, but not nuclear, DNA diversity is much lower in D. recens than in D. subquinaria, consistent with the hypothesized diversity-purging effects of an evolutionarily recent Wolbachia sweep. The d(N)/d(S) ratio in mtDNA is significantly greater in D. recens, suggesting that Muller's ratchet has brought about an increased rate of substitution of slightly deleterious mutations. The data also reveal elevated rates of synonymous substitutions in D. recens, suggesting that these sites may experience weak selection. These findings show that maternally transmitted endosymbionts can severely depress levels of mtDNA diversity within an infected host species, while accelerating the rate of divergence among mtDNA lineages in different species.

Characterization of active R2 retrotransposition in the rDNA locus of Drosophila simulans

The rRNA gene (rDNA) loci of all arthropod lineages contain non-LTR retrotransposable elements that have evolved to specifically insert into the 28S rRNA genes. Extensive in vitro experiments have been conducted to investigate the mechanism of R2 retrotransposition but little is known of the insertion frequency or cellular factors that might regulate R2 activity. In this article, isofemale lines obtained from a population of Drosophila simulans were surveyed for recent R2 insertions. Within most lines, all individuals showed the same collection of R2 insertions, providing no evidence for recent R2 activity. However, in a few of the isofemale lines, virtually all individuals differed in their R2 insertion profiles. The descendants of individual pairs of flies from these "active lines" rapidly accumulated new insertions. The frequent insertion of new R2 elements was associated with the elimination of old R2 elements from the rDNA locus. The existence of lines in which R2 retrotransposes frequently and lines in which the elements appear dormant suggests that cellular mechanisms that can regulate the activity of R2 exist. Retrotransposition activity was correlated with the number of full-length R2 elements but not with the size of the rDNA locus or the number of uninserted units.

Evolution of the Transposable Element Pokey in the Ribosomal DNA of Species in the Subgenus Daphnia (Crustacea: Cladocera)

Pokey is a member of the piggyBac (previously called the TTAA-specific) family of transposons and inserts into a conserved region of the large subunit ribosomal RNA gene. This location is a "hot spot" for insertional activity, as it is known to contain other arthropod transposable elements. However, Pokey is unique in that it is the first DNA transposon yet known to insert into this region. All other insertions are class I non-LTR retrotransposons. This study surveyed variation in Pokey elements through phylogenetic analysis of the 3' ends of Pokey elements from ribosomal DNA (rDNA) in species from the nominate subgenus of the genus Daphnia (Crustacea: Cladocera). The results suggest that Pokey has been stably, vertically inherited within rDNA over long periods of evolutionary time. No evidence was found to support horizontal transfer, which commonly occurs in other DNA transposons, such as P and mariner. Furthermore, Pokey has diverged into sublineages that have persisted across speciation events in some groups. In addition, a new highly divergent paralogous Pokey element was discovered in the rDNA of one species.

R1 and R2 retrotransposition and deletion in the rDNA loci on the X and Y chromosomes of Drosophila melanogaster

The non-LTR retrotransposons R1 and R2 insert into the 28S rRNA genes of arthropods. Comparisons among Drosophila lineages have shown that these elements are vertically inherited, while studies within species have indicated a rapid turnover of individual copies (elimination of old copies and the insertion of new copies). To better understand the turnover of R1 and R2, 200 retrotranspositions and nearly 100 eliminations have been scored in the Harwich mutation-accumulation lines of Drosophila melanogaster. Because the rDNA arrays in D. melanogaster are present on the X and Y chromosomes and no exchanges were detected in these lines, it was possible to show that R1 retrotranspositions occur predominantly in the male germ line, while R2 retrotranspositions were more evenly divided between the germ lines of both sexes. The rate of elimination of elements from the Y rDNA array was twice that of the X rDNA array with both chromosomal loci containing regions where the rate of elimination was on average eight times higher. Most R1 and R2 eliminations appear to occur by large intrachromosomal events (i.e., loop-out events) that involve multiple rDNA units. These findings are interpreted in light of the known abundance of R1 and R2 elements in the X and Y rDNA loci of D. melanogaster.

Molecular phylogeny of the Drosophila melanogaster species subgroup

Although molecular and phenotypic evolution have been studied extensively in Drosophila melanogaster and its close relatives, phylogenetic relationships within the D. melanogaster species subgroup remain unresolved. In particular, recent molecular studies have not converged on the branching orders of the D. yakuba-D. teissieri and D. erecta-D. orena species pairs relative to the D. melanogaster-D. simulans-D. mauritiana-D. sechellia species complex. Here, we reconstruct the phylogeny of the melanogaster species subgroup using DNA sequence data from four nuclear genes. We have employed "vectorette PCR" to obtain sequence data for orthologous regions of the Alcohol dehydrogenase (Adh), Alcohol dehydrogenase related (Adhr), Glucose dehydrogenase (Gld), and rosy (ry) genes (totaling 7164 bp) from six melanogaster subgroup species (D. melanogaster, D. simulans, D. teissieri, D. yakuba, D. erecta, and D. orena) and three species from subgroups outside the melanogaster species subgroup [D. eugracilis (eugracilis subgroup), D. mimetica (suzukii subgroup), and D. lutescens (takahashii subgroup)]. Relationships within the D. simulans complex are not addressed. Phylogenetic analyses employing maximum parsimony, neighbor-joining, and maximum likelihood methods strongly support a D. yakuba-D. teissieri and D. erecta-D. orena clade within the melanogaster species subgroup. D. eugracilis is grouped closer to the melanogaster subgroup than a D. mimetica-D. lutescens clade. This tree topology is supported by reconstructions employing simple (single parameter) and more complex (nonreversible) substitution models.

Complex evolution of gypsy in Drosophilid species

In an endeavor to contribute to the comprehension of the evolution of transposable elements (TEs) in the genome of host species, we investigated the phylogenetic relationships of sequences homologous to the retrotransposon gypsy of Drosophila melanogaster in 19 species of Drosophila, in Scaptodrosophila latifasciaeformis, and in Zaprionus indianus. This phylogenetic study was based on approximately 500 base pairs of the env gene. Our analyses showed considerable discrepancy between the phylogeny of gypsy elements and the relationship of their host species, and they allow us to infer a complex evolutionary pattern that could include ancestral polymorphism, vertical transmission, and several cases of horizontal transmission.

Evolution of multiple components of virulence in Drosophila-nematode associations

Virulence is of central importance in host-parasite interactions, yet little is known about how it changes over extended evolutionary periods. In this study, all four species in the testacea species group of Drosophila were experimentally infected with sympatric and allopatric nematodes in the Howardula aoronymphium species complex, and the effect of parasite infection on three components of host fitness was determined. The Drosophila species show striking differences in their responses to infection, with reductions reaching 80% in adult lifespan, 100% in female fertility, and 90% in male fertility. Female sterility appears to be determined by the host; species that are sterilized by their local nematodes are also sterilized by the other allopatric nematodes in the H. aoronymphium complex. Host species that are not sterilized by their local parasite are not sterilized by other nematodes in the complex. In contrast, reductions in host adult lifespan and male fertility depend on both the host and the parasite. Whereas all nematodes reduced the survival of their local host species equally (about 40-45%), survival of two host species was drastically reduced (about 80%) when infected with an allopatric parasite. Thus, virulence is evolutionarily labile in associations between Drosophila testacea group species and their Howardula parasites. The data suggest that changes in the sterility component of virulence are due primarily to host evolution, whereas changes in the host mortality component are due in large part to parasite evolution.

Transcription of endogenous and exogenous R2 elements in the rRNA gene locus of Drosophila melanogaster

R2 retrotransposons insert into the rRNA-encoding units (rDNA units) that form the nucleoli of insects. We have utilized an R2 integration system in Drosophila melanogaster to study transcription of foreign sequences integrated into the R2 target site of the 28S rRNA genes. The exogenous sequences were cotranscribed at dramatically different levels which closely paralleled the level of transcription of the endogenous R1 and R2 elements. Transcription levels were inversely correlated with the number of uninserted rDNA units, variation in this number having been brought about by the R2 integration system itself. Females with as few as 20 uninserted rDNA units per X chromosome had expression levels of endogenous and exogenous insertion sequences that were 2 orders of magnitude higher than lines that contained over 80 uninserted rDNA units per chromosome. R2 insertions only 167 bp in length exhibited this range of transcriptional regulation. Analysis of transcript levels in males suggested R2 insertions on the Y chromosome are not down-regulated to the same extent as insertions on the X chromosome. These results suggest that transcription of the rDNA units can be tightly regulated, but this regulation gradually breaks down as the cell approaches the minimum number of uninserted genes needed for survival.

Rates of R1 and R2 retrotransposition and elimination from the rDNA locus of Drosophila melanogaster

R1 and R2 elements are non-LTR retrotransposons that insert specifically into the 28S rRNA genes of arthropods. The process of concerted evolution of the rDNA locus should give rise to rapid turnover of these mobile elements compared to elements that insert at sites throughout a genome. To estimate the rate of R1 and R2 turnover we have examined the insertion of new elements and elimination of old elements in the Harwich mutation accumulation lines of Drosophila melanogaster, a set of inbred lines maintained for >350 generations. Nearly 300 new insertion and elimination events were observed in the 19 Harwich lines. The retrotransposition rate for R1 was 18 times higher than the retrotransposition rate for R2. Both rates were within the range previously found for retrotransposons that insert outside the rDNA loci in D. melanogaster. The elimination rates of R1 and R2 from the rDNA locus were similar to each other but over two orders of magnitude higher than that found for other retrotransposons. The high rates of R1 and R2 elimination from the rDNA locus confirm that these elements must maintain relatively high rates of retrotransposition to ensure their continued presence in this locus.

Coevolution of the telomeric retrotransposons across Drosophila species

As in other eukaryotes, telomeres in Drosophila melanogaster are composed of long arrays of repeated DNA sequences. Remarkably, in D. melanogaster these repeats are produced, not by telomerase, but by successive transpositions of two telomere-specific retrotransposons, HeT-A and TART. These are the only transposable elements known to be completely dedicated to a role in chromosomes, a finding that provides an opportunity for investigating questions about the evolution of telomeres, telomerase, and the transposable elements themselves. Recent studies of D. yakuba revealed the presence of HeT-A elements with precisely the same unusual characteristics as HeT-A(mel) although they had only 55% nucleotide sequence identity. We now report that the second element, TART, is also a telomere component in D. yakuba; thus, these two elements have been evolving together since before the separation of the melanogaster and yakuba species complexes. Like HeT-A(yak), TART(yak) is undergoing concerted sequence evolution, yet they retain the unusual features TART(mel) shares with HeT-A(mel). There are at least two subfamilies of TART(yak) with significantly different sequence and expression. Surprisingly, one subfamily of TART(yak) has >95% sequence identity with a subfamily of TART(mel) and shows similar transcription patterns. As in D. melanogaster, other retrotransposons are excluded from the D. yakuba terminal arrays studied to date.

Dynamics of R1 and R2 elements in the rDNA locus of Drosophila simulans

The mobile elements R1 and R2 insert specifically into the rRNA gene locus (rDNA locus) of arthropods, a locus known to undergo concerted evolution, the recombinational processes that preserve the sequence homogeneity of all repeats. To monitor how rapidly individual R1 and R2 insertions are turned over in the rDNA locus by these processes, we have taken advantage of the many 5' truncation variants that are generated during the target-primed reverse transcription mechanism used by these non-LTR retrotransposons for their integration. A simple PCR assay was designed to reveal the pattern of the 5' variants present in the rDNA loci of individual X chromosomes in a population of Drosophila simulans. Each rDNA locus in this population was found to have a large, unique collection of 5' variants. Each variant was present at low copy number, usually one copy per chromosome, and was seldom distributed to other chromosomes in the population. The failure of these variants to spread to other units in the same rDNA locus suggests a strong recombinational bias against R1 and R2 that results in the individual copies of these elements being rapidly lost from the rDNA locus. This bias suggests a significantly higher frequency of R1 and R2 retrotransposition than we have previously suggested.