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

Insights into ribosomal DNA dominance and magnification through characterization of isogenic deletion alleles

The major loci for the large primary ribosomal RNA (rRNA) genes (35S rRNAs) exist as hundreds to thousands of tandem repeats in all organisms and dozens to hundreds in Drosophila. The highly repetitive nature of the ribosomal DNA (rDNA) makes it intrinsically unstable, and many conditions arise from the reduction in or magnification of copy number, but the conditions under which it does so remain unknown. By targeted DNA damage to the rDNA of the Y chromosome, we created and investigated a series of rDNA alleles. We found that complete loss of rDNA leads to lethality after the completion of embryogenesis, blocking larval molting and metamorphosis. We find that the resident retrotransposons-R1 and R2-are regulated by active rDNA such that reduction in copy number derepresses these elements. Their expression is highest during the early first instar, when loss of rDNA is lethal. Regulation of R1 and R2 may be related to their structural arrangement within the rDNA, as we find they are clustered in the flanks of the nucleolus organizing region (NOR; the cytological appearance of the rDNA). We assessed the complex nucleolar dominance relationship between X- and Y-linked rDNA using a histone H3.3-GFP reporter construct and incorporation at the NOR and found that dominance is controlled by rDNA copy number as at high multiplicity the Y-linked array is dominant, but at low multiplicity the X-linked array becomes derepressed. Finally, we found that multiple conditions that disrupt nucleolar dominance lead to increased rDNA magnification, suggesting that the phenomena of dominance and magnification are related, and a single mechanism may underlie and unify these two longstanding observations in Drosophila.

The Dynamic Interplay Between Ribosomal DNA and Transposable Elements: A Perspective From Genomics and Cytogenetics

Although both are salient features of genomes, at first glance ribosomal DNAs and transposable elements are genetic elements with not much in common: whereas ribosomal DNAs are mainly viewed as housekeeping genes that uphold all prime genome functions, transposable elements are generally portrayed as selfish and disruptive. These opposing characteristics are also mirrored in other attributes: organization in tandem (ribosomal DNAs) versus organization in a dispersed manner (transposable elements); evolution in a concerted manner (ribosomal DNAs) versus evolution by diversification (transposable elements); and activity that prolongs genomic stability (ribosomal DNAs) versus activity that shortens it (transposable elements). Re-visiting relevant instances in which ribosomal DNA-transposable element interactions have been reported, we note that both repeat types share at least four structural and functional hallmarks: (1) they are repetitive DNAs that shape genomes in evolutionary timescales, (2) they exchange structural motifs and can enter co-evolution processes, (3) they are tightly controlled genomic stress sensors playing key roles in senescence/aging, and (4) they share common epigenetic marks such as DNA methylation and histone modification. Here, we give an overview of the structural, functional, and evolutionary characteristics of both ribosomal DNAs and transposable elements, discuss their roles and interactions, and highlight trends and future directions as we move forward in understanding ribosomal DNA-transposable element associations.

Harnessing eukaryotic retroelement proteins for transgene insertion into human safe-harbor loci

Current approaches for inserting autonomous transgenes into the genome, such as CRISPR-Cas9 or virus-based strategies, have limitations including low efficiency and high risk of untargeted genome mutagenesis. Here, we describe precise RNA-mediated insertion of transgenes (PRINT), an approach for site-specifically primed reverse transcription that directs transgene synthesis directly into the genome at a multicopy safe-harbor locus. PRINT uses delivery of two in vitro transcribed RNAs: messenger RNA encoding avian R2 retroelement-protein and template RNA encoding a transgene of length validated up to 4 kb. The R2 protein coordinately recognizes the target site, nicks one strand at a precise location and primes complementary DNA synthesis for stable transgene insertion. With a cultured human primary cell line, over 50% of cells can gain several 2 kb transgenes, of which more than 50% are full-length. PRINT advantages include no extragenomic DNA, limiting risk of deleterious mutagenesis and innate immune responses, and the relatively low cost, rapid production and scalability of RNA-only delivery.

R2 and Non-Site-Specific R2-Like Retrotransposons of the German Cockroach, Blattella germanica

The structural and functional organization of the ribosomal RNA gene cluster and the full-length R2 non-LTR retrotransposon (integrated into a specific site of 28S ribosomal RNA genes) of the German cockroach, Blattella germanica, is described. A partial sequence of the R2 retrotransposon of the cockroach Rhyparobia maderae is also analyzed. The analysis of previously published next-generation sequencing data from the B. germanica genome reveals a new type of retrotransposon closely related to R2 retrotransposons but with a random distribution in the genome. Phylogenetic analysis reveals that these newly described retrotransposons form a separate clade. It is shown that proteins corresponding to the open reading frames of newly described retrotransposons exhibit unequal structural domains. Within these retrotransposons, a recombination event is described. New mechanism of transposition activity is discussed. The essential structural features of R2 retrotransposons are conserved in cockroaches and are typical of previously described R2 retrotransposons. However, the investigation of the number and frequency of 5′-truncated R2 retrotransposon insertion variants in eight B. germanica populations suggests recent mobile element activity. It is shown that the pattern of 5′-truncated R2 retrotransposon copies can be an informative molecular genetic marker for revealing genetic distances between insect populations.

Obligatory parthenogenesis and TE load: Bacillus stick insects and the R2 non-LTR retrotransposon

Transposable elements (TEs) are selfish genetic elements whose self-replication is contrasted by the host genome. In this context, host reproductive strategies are predicted to impact on both TEs load and activity. The presence and insertion distribution of the non-LTR retrotransposon R2 was here studied in populations of the strictly bisexual Bacillus grandii maretimi and of the obligatory parthenogenetic Bacillus atticus atticus. Furthermore, data were also obtained from the offspring of selected B. a. atticus females. At the population level, the gonochoric B. g. maretimi showed a significantly higher R2 load than the obligatory parthenogenetic B. a. atticus. The comparison with bisexual and unisexual Bacillus rossius populations showed that their values were higher than those recorded for B. a. atticus and similar, or even higher, than those of B. g. maretimi. Consistently, an R2 load reduction is scored in B. a. atticus offspring even if with a great variance. On the whole, data here produced indicate that in the obligatory unisexual B. a. atticus R2 is active and that mechanisms of molecular turnover are effective. Furthermore, progeny analyses show that, at variance of the facultative parthenogenetic B. rossius, the R2 activity is held at a lower rate. Modeling parental-offspring inheritance, suggests that in B. a. atticus recombination plays a major role in eliminating insertions rather than selection, as previously suggested for unisexual B. rossius progeny, even if in both cases a high variance is observed. In addition to this, mechanisms of R2 silencing or chances of clonal selection cannot be ruled out.

Hybridogenesis and a potential case of R2 non-LTR retrotransposon horizontal transmission in Bacillus stick insects (Insecta Phasmida)

Horizontal transfer (HT) is an event in which the genetic material is transferred from one species to another, even if distantly related, and it has been demonstrated as a possible essential part of the lifecycle of transposable elements (TEs). However, previous studies on the non-LTR R2 retrotransposon, a metazoan-wide distributed element, indicated its vertical transmission since the Radiata-Bilateria split. Here we present the first possible instances of R2 HT in stick insects of the genus Bacillus (Phasmida). Six R2 elements were characterized in the strictly bisexual subspecies B. grandii grandii, B. grandii benazzii and B. grandii maretimi and in the obligatory parthenogenetic taxon B. atticus. These elements were compared with those previously retrieved in the facultative parthenogenetic species B. rossius. Phylogenetic inconsistencies between element and host taxa, and age versus divergence analyses agree and support at least two HT events. These HT events can be explained by taking into consideration the complex Bacillus reproductive biology, which includes also hybridogenesis, gynogenesis and androgenesis. Through these non-canonical reproductive modes, R2 elements may have been transferred between Bacillus genomes. Our data suggest, therefore, a possible role of hybridization for TEs survival and the consequent reshaping of involved genomes.

Integration, Regulation, and Long-Term Stability of R2 Retrotransposons

R2 elements are sequence specific non-LTR retrotransposons that exclusively insert in the 28S rRNA genes of animals. R2s encode an endonuclease that cleaves the insertion site and a reverse transcriptase that uses the cleaved DNA to prime reverse transcription of the R2 transcript, a process termed target primed reverse transcription. Additional unusual properties of the reverse transcriptase as well as DNA and RNA binding domains of the R2 encoded protein have been characterized. R2 expression is through co-transcription with the 28S gene and self-cleavage by a ribozyme encoded at the R2 5' end. Studies in laboratory stocks and natural populations of Drosophila suggest that R2 expression is tied to the distribution of R2-inserted units within the rDNA locus. Most individuals have no R2 expression because only a small fraction of their rRNA genes need to be active, and a contiguous region of the locus free of R2 insertions can be selected for activation. However, if the R2-free region is not large enough to produce sufficient rRNA, flanking units - including those inserted with R2 - must be activated. Finally, R2 copies rapidly turnover within the rDNA locus, yet R2 has been vertically maintained in animal lineages for hundreds of millions of years. The key to this stability is R2's ability to remain dormant in rDNA units outside the transcribed regions for generations until the stochastic nature of the crossovers that drive the concerted evolution of the rDNA locus inevitably reshuffle the inserted and uninserted units, resulting in transcription of the R2-inserted units.

Dead element replicating: degenerate R2 element replication and rDNA genomic turnover in the Bacillus rossius stick insect (Insecta: Phasmida)

R2 is an extensively investigated non-LTR retrotransposon that specifically inserts into the 28S rRNA gene sequences of a wide range of metazoans, disrupting its functionality. During R2 integration, first strand synthesis can be incomplete so that 5’ end deleted copies are occasionally inserted. While active R2 copies repopulate the locus by retrotransposing, the non-functional truncated elements should frequently be eliminated by molecular drive processes leading to the concerted evolution of the rDNA array(s). Although, multiple R2 lineages have been discovered in the genome of many animals, the rDNA of the stick insect Bacillus rossius exhibits a peculiar situation: it harbors both a canonical, functional R2 element (R2Br^fun) as well as a full-length but degenerate element (R2Br^deg). An intensive sequencing survey in the present study reveals that all truncated variants in stick insects are present in multiple copies suggesting they were duplicated by unequal recombination. Sequencing results also demonstrate that all R2Br^deg copies are full-length, i. e. they have no associated 5' end deletions, and functional assays indicate they have lost the active ribozyme necessary for R2 RNA maturation. Although it cannot be completely ruled out, it seems unlikely that the degenerate elements replicate via reverse transcription, exploiting the R2Br^fun element enzymatic machinery, but rather via genomic amplification of inserted 28S by unequal recombination. That inactive copies (both R2Br^deg or 5'-truncated elements) are not eliminated in a short term in stick insects contrasts with findings for the Drosophila R2, suggesting a widely different management of rDNA loci and a lower efficiency of the molecular drive while achieving the concerted evolution.

Evolutionary dynamics of R2 retroelement and insertion inheritance in the genome of bisexual and parthenogenetic Bacillus rossius populations (Insecta Phasmida)

Theoretical and empirical studies have shown differential management of transposable elements in organisms with different reproductive strategies. To investigate this issue, we analysed the R2 retroelement structure and variability in parthenogenetic and bisexual populations of Bacillus rossius stick insects, as well as insertions inheritance in the offspring of parthenogenetic isolates and of crosses. The B. rossius genome hosts a functional (R2Br(fun) ) and a degenerate (R2Br(deg) ) element, their presence correlating with neither reproductive strategies nor population distribution. The median-joining network method indicated that R2Br(fun) duplicates through a multiple source model, while R2Br(deg) is apparently still duplicating via a master gene model. Offspring analyses showed that unisexual and bisexual offspring have a similar number of R2Br-occupied sites. Multiple or recent shifts from gonochoric to parthenogenetic reproduction may explain the observed data. Moreover, insertion frequency spectra show that higher-frequency insertions in unisexual offspring significantly outnumber those in bisexual offspring. This suggests that unisexual offspring eliminate insertions with lower efficiency. A comparison with simulated insertion frequencies shows that inherited insertions in unisexual and bisexual offspring are significantly different from the expectation. On the whole, different mechanisms of R2 elimination in unisexual vs bisexual offspring and a complex interplay between recombination effectiveness, natural selection and time can explain the observed data.

Preferential occupancy of R2 retroelements on the B chromosomes of the grasshopper Eyprepocnemis plorans

R2 non-LTR retrotransposons exclusively insert into the 28S rRNA genes of their host, and are expressed by co-transcription with the rDNA unit. The grasshopper Eyprepocnemis plorans contains transcribed rDNA clusters on most of its A chromosomes, as well as non-transcribed rDNA clusters on the parasitic B chromosomes found in many populations. Here the structure of the E. plorans R2 element, its abundance relative to the number of rDNA units and its retrotransposition activity were determined. Animals screened from five populations contained on average over 12,000 rDNA units on their A chromosomes, but surprisingly only about 100 R2 elements. Monitoring the patterns of R2 insertions in individuals from these populations revealed only low levels of retrotransposition. The low rates of R2 insertion observed in E. plorans differ from the high levels of R2 insertion previously observed in insect species that have many fewer rDNA units. It is proposed that high levels of R2 are strongly selected against in E. plorans, because the rDNA transcription machinery in this species is unable to differentiate between R2-inserted and uninserted units. The B chromosomes of E. plorans contain an additional 7,000 to 15,000 rDNA units, but in contrast to the A chromosomes, from 150 to over 1,500 R2 elements. The higher concentration of R2 in the inactive B chromosomes rDNA clusters suggests these chromosomes can act as a sink for R2 insertions thus further reducing the level of insertions on the A chromosomes. These studies suggest an interesting evolutionary relationship between the parasitic B chromosomes and R2 elements.

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.

28S junctions and chimeric elements of the rDNA targeting non-LTR retrotransposon R2 in crustacean living fossils (Branchiopoda, Notostraca)

The 28S rRNA genes of several metazoans are interrupted by site-specific targeting non-LTR retrotransposons, such as R2. R2 elements have been deeply analyzed but aspects of their retrotransposition mechanism and the origin of the wide diversity observed are still debated. We characterized six new R2 lineages in four tadpole shrimp species (Notostraca), samples deriving from a parthenogenetic population of Triops cancriformis (R2Tc_it) and from bisexual Lepidurus populations of L. lubbocki (R2Ll), L. couesii (R2LcA, R2LcB, R2LcC) and L. arcticus (R2La). All elements fit the canonical R2 structure but R2Ll which turned out to be a chimera with an additional ORF originating from another R2. Consistently with data on LINEs, R2Ll could be the result of recombination due to reverse transcriptase template jump. The analysis of 28S/R2 5' end junctions further suggests aberrant homologous recombination, as observed in RNA viruses.

Retrotransposition of R2 elements in somatic nuclei during the early development of Drosophila

Background

R2 retrotransposable elements exclusively insert in the 28S rRNA genes of their host. Their RNA transcripts are produced by self-processing from a 28S R2 cotranscript. Because full-length R2 transcripts are found in most tissues of R2-active animals, we tested whether new R2 insertions occurred in somatic tissues even though such events would be an evolutionary dead end.

Findings

PCR assays were used to identify somatic R2 insertions in isolated adult tissues and larval imaginal discs of Drosophila simulans. R2 somatic mosaics were detected encompassing cells from individual tissues as well as tissues from multiple body segments. The somatic insertions had 5' junction sequences characteristic of germline insertions suggesting they represented authentic retrotransposition events.

Conclusions

Body segments are specified early in Drosophila development, thus the detection of the same somatic insertion in cells from multiple tissues suggested that the R2 retrotransposition events had occurred before the blastoderm stage of Drosophila development. R2 activity at this stage, when embryonic nuclei are rapidly dividing in a common cytoplasm, suggests that some retrotransposition events appearing as germline events may correspond to germline mosaicism.

R2 dynamics in Triops cancriformis (Bosc, 1801) (Crustacea, Branchiopoda, Notostraca): turnover rate and 28S concerted evolution

The R2 retrotransposon is here characterized in bisexual populations of the European crustacean Triops cancriformis. The isolated element matches well with the general aspects of the R2 family and it is highly differentiated from that of the congeneric North American Triops longicaudatus. The analysis of 5' truncations indicates that R2 dynamics in T. cancriformis populations show a high turnover rate as observed in Drosophila simulans. For the first time in the literature, though, individuals harboring truncation variants, but lacking the complete element, are found. Present results suggest that transposition-mediated deletion mechanisms, possibly involving genomic turnover processes acting on rDNAs, can dramatically decrease the copy number or even delete R2 from the ribosomal locus. The presence of R2 does not seem to impact on the nucleotide variation of inserted 28S rDNA with respect to the uninserted genes. On the other hand, a low level of polymorphism characterizes rDNA units because new 28S variants continuously spread across the ribosomal array. Again, the interplay between transposition-mediated deletion and molecular drive may explain this pattern.

A subtelomeric non-LTR retrotransposon Hebe in the bdelloid rotifer Adineta vaga is subject to inactivation by deletions but not 5' truncations

BACKGROUND

Rotifers of the class Bdelloidea are microscopic freshwater invertebrates best known for: their capacity for anhydrobiosis; the lack of males and meiosis; and for the ability to capture genes from other non-metazoan species. Although genetic exchange between these animals might take place by non-canonical means, the overall lack of meiosis and syngamy should greatly impair the ability of transposable elements (TEs) to spread in bdelloid populations. Previous studies demonstrated that bdelloid chromosome ends, in contrast to gene-rich regions, harbour various kinds of TEs, including specialized telomere-associated retroelements, as well as DNA TEs and retrovirus-like retrotransposons which are prone to horizontal transmission. Vertically-transmitted retrotransposons have not previously been reported in bdelloids and their identification and studies of the patterns of their distribution and evolution could help in the understanding of the high degree of TE compartmentalization within bdelloid genomes.

RESULTS

We identified and characterized a non-long terminal repeat (LTR) retrotransposon residing primarily in subtelomeric regions of the genome in the bdelloid rotifer Adineta vaga. Contrary to the currently prevailing views on the mode of proliferation of non-LTR retrotransposons, which results in frequent formation of 5'-truncated ('dead-on-arrival') copies due to the premature disengagement of the element-encoded reverse transcriptase from its template, this non-LTR element, Hebe, is represented only by non-5'-truncated copies. Most of these copies, however, were subject to internal deletions associated with microhomologies, a hallmark of non-homologous end-joining events.

CONCLUSIONS

The non-LTR retrotransposon Hebe from the bdelloid rotifer A. vaga was found to undergo frequent microhomology-associated deletions, rather than 5'-terminal truncations characteristic of this class of retrotransposons, and to exhibit preference for telomeric localization. These findings represent the first example of a vertically transmitted putatively deleterious TE in bdelloids, and may indicate the involvement of microhomology-mediated non-homologous end-joining in desiccation-induced double-strand break repair at the genome periphery.

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.

The effect of transposon Pokey insertions on sequence variation in the 28S rRNA gene of Daphnia pulex

The goal of this study was to determine the impact of breeding system and the presence of the transposon Pokey on intraindividual variation in 28S rRNA genes. We PCR-amplified, cloned, and sequenced 1000 nucleotides downstream of the Pokey insertion site in genes with and without insertions from 10 obligately and 10 cyclically parthenogenetic isolates of Daphnia pulex. Variation among genes with Pokey insertions was higher than variation among genes without insertions in both cyclic and obligate isolates. Although the differences were not quite significant (p = 0.06 in both cases), the results suggest that Pokey insertions are likely to inhibit the homogenization of their host genes to some extent. We also observed that the complement of 28S rRNA alleles differed between genes with and without inserts in some isolates, suggesting that a particular inserted gene can persist for substantial periods of time and even spread within the rDNA array, despite the fact that insertions are deleterious. This apparently contradictory pattern can be explained if homogenization of rRNA genes occurs primarily by gene conversion, but copies with Pokey inserts can occasionally increase in frequency within arrays owing to unequal crossing over events that do not originate in the inserted genes themselves.

The pattern of R2 retrotransposon activity in natural populations of Drosophila simulans reflects the dynamic nature of the rDNA locus

The pattern and frequency of insertions that enable transposable elements to remain active in a population are poorly understood. The retrotransposable element R2 exclusively inserts into the 28S rRNA genes where it establishes long-term, stable relationships with its animal hosts. Previous studies with laboratory stocks of Drosophila simulans have suggested that control over R2 retrotransposition resides within the rDNA loci. In this report, we sampled 180 rDNA loci of animals collected from two natural populations of D. simulans. The two populations were found to have similar patterns of R2 activity. About half of the rDNA loci supported no or very low levels of R2 transcripts with no evidence of R2 retrotransposition. The remaining half of the rDNA loci had levels of R2 transcripts that varied in a continuous manner over almost a 100-fold range and did support new retrotransposition events. Structural analysis of the rDNA loci in 18 lines that spanned the range of R2 transcript levels in these populations revealed that R2 number and rDNA locus size varied 2-fold; however, R2 activity was not readily correlated with either of these parameters. Instead R2 activity was best correlated with the distribution of elements within the rDNA locus. Loci with no activity had larger contiguous blocks of rDNA units free of R2-insertions. These data suggest a model in which frequent recombination within the rDNA locus continually redistributes R2-inserted units resulting in changing levels of R2 activity within individual loci and persistent R2 activity within the population.

High Ag-NOR-site variation associated to a secondary contact in brown trout from the Iberian Peninsula

The analysis of nucleolar organizer regions (NORs) using silver (Ag-) staining and in situ hybridization (ISH) in brown trout (Salmo trutta) from various river basins in the Iberian Peninsula revealed high variation in the number and location of NORs. A total of 17 different Ag-NOR sites were revealed in 10 different chromosome pairs. Three different Ag-NOR patterns clustered by river basins and strongly associated to the internal transcribed spacer 1 (ITS1) variation were detected. The main variability in NOR-sites was found in a secondary contact between two divergent lineages of brown trouts at Duero basin. Our results confirmed the abrupt break in the spatial distribution of genetic variation of brown trout populations previously reported at Duero basin. We hypothesize that NOR-site variation might be a consequence of hybridization between divergent lineages of brown trouts and that NORs could play a major role in the maintenance of a hybrid zone in Duero basin via post-zygotic isolation mechanisms.

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.

Epigenetic regulation of retrotransposons within the nucleolus of Drosophila

R2 retrotransposable elements exclusively insert into a conserved region of the tandemly organized 28S rRNA genes. Despite inactivating a subset of these genes, R2 elements have persisted in the ribosomal DNA (rDNA) loci of insects for hundreds of millions of years. Controlling R2 proliferation was addressed in this study using lines of Drosophila simulans previously shown to have either active or inactive R2 retrotransposition. Lines with active retrotransposition were shown to have high R2 transcript levels, which nuclear run-on transcription experiments revealed were due to increased transcription of R2-inserted genes. Crosses between R2 active and inactive lines indicated that an important component of this transcriptional control is linked to or near the rDNA locus, with the R2 transcription level of the inactive parent being dominant. Pulsed-field gel analysis suggested that the R2 active and inactive states were determined by R2 distribution within the locus. Molecular and cytological analyses further suggested that the entire rDNA locus from the active line can be silenced in favor of the locus from the inactive line. This silencing of entire rDNA loci represents an example of the large-scale epigenetic control of transposable elements and shares features with the nucleolar dominance frequently seen in interspecies hybrids.

Sequence variation within the rRNA gene loci of 12 Drosophila species

Concerted evolution maintains at near identity the hundreds of tandemly arrayed ribosomal RNA (rRNA) genes and their spacers present in any eukaryote. Few comprehensive attempts have been made to directly measure the identity between the rDNA units. We used the original sequencing reads (trace archives) available through the whole-genome shotgun sequencing projects of 12 Drosophila species to locate the sequence variants within the 7.8-8.2 kb transcribed portions of the rDNA units. Three to 18 variants were identified in >3% of the total rDNA units from 11 species. Species where the rDNA units are present on multiple chromosomes exhibited only minor increases in sequence variation. Variants were 10-20 times more abundant in the noncoding compared with the coding regions of the rDNA unit. Within the coding regions, variants were three to eight times more abundant in the expansion compared with the conserved core regions. The distribution of variants was largely consistent with models of concerted evolution in which there is uniform recombination across the transcribed portion of the unit with the frequency of standing variants dependent upon the selection pressure to preserve that sequence. However, the 28S gene was found to contain fewer variants than the 18S gene despite evolving 2.5-fold faster. We postulate that the fewer variants in the 28S gene is due to localized gene conversion or DNA repair triggered by the activity of retrotransposable elements that are specialized for insertion into the 28S genes of these species.

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.

Characteristics of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) rRNA genes of Apis mellifera (Insecta: Hymenoptera): structure, organization, and retrotransposable elements

As an accompanying manuscript to the release of the honey bee genome, we report the entire sequence of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) ribosomal RNA (rRNA)-encoding gene sequences (rDNA) and related internally and externally transcribed spacer regions of Apis mellifera (Insecta: Hymenoptera: Apocrita). Additionally, we predict secondary structures for the mature rRNA molecules based on comparative sequence analyses with other arthropod taxa and reference to recently published crystal structures of the ribosome. In general, the structures of honey bee rRNAs are in agreement with previously predicted rRNA models from other arthropods in core regions of the rRNA, with little additional expansion in non-conserved regions. Our multiple sequence alignments are made available on several public databases and provide a preliminary establishment of a global structural model of all rRNAs from the insects. Additionally, we provide conserved stretches of sequences flanking the rDNA cistrons that comprise the externally transcribed spacer regions (ETS) and part of the intergenic spacer region (IGS), including several repetitive motifs. Finally, we report the occurrence of retrotransposition in the nuclear large subunit rDNA, as R2 elements are present in the usual insertion points found in other arthropods. Interestingly, functional R1 elements usually present in the genomes of insects were not detected in the honey bee rRNA genes. The reverse transcriptase products of the R2 elements are deduced from their putative open reading frames and structurally aligned with those from another hymenopteran insect, the jewel wasp Nasonia (Pteromalidae). Stretches of conserved amino acids shared between Apis and Nasonia are illustrated and serve as potential sites for primer design, as target amplicons within these R2 elements may serve as novel phylogenetic markers for Hymenoptera. Given the impending completion of the sequencing of the Nasonia genome, we expect our report eventually to shed light on the evolution of the hymenopteran genome within higher insects, particularly regarding the relative maintenance of conserved rDNA genes, related variable spacer regions and retrotransposable elements.

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.

Organization of chromosome ends in the rice blast fungus, Magnaporthe oryzae

Eukaryotic pathogens of humans often evade the immune system by switching the expression of surface proteins encoded by subtelomeric gene families. To determine if plant pathogenic fungi use a similar mechanism to avoid host defenses, we sequenced the 14 chromosome ends of the rice blast pathogen, Magnaporthe oryzae. One telomere is directly joined to ribosomal RNA-encoding genes, at the end of the ∼2 Mb rDNA array. Two are attached to chromosome-unique sequences, and the remainder adjoin a distinct subtelomere region, consisting of a telomere-linked RecQ-helicase (TLH) gene flanked by several blocks of tandem repeats. Unlike other microbes, M.oryzae exhibits very little gene amplification in the subtelomere regions—out of 261 predicted genes found within 100 kb of the telomeres, only four were present at more than one chromosome end. Therefore, it seems unlikely that M.oryzae uses switching mechanisms to evade host defenses. Instead, the M.oryzae telomeres have undergone frequent terminal truncation, and there is evidence of extensive ectopic recombination among transposons in these regions. We propose that the M.oryzae chromosome termini play more subtle roles in host adaptation by promoting the loss of terminally-positioned genes that tend to trigger host defenses.

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.

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.

Long-term inheritance of the 28S rDNA-specific retrotransposon R2

R2 is a non-long-terminal-repeat (LTR) retrotransposon that inserts specifically into 28S rDNA. R2 has been identified in many species of arthropods and three species of chordates. R2 may be even more widely distributed in animals, and its origin may be traceable to early animal evolution. In this study, we identified R2 elements in medaka fish, White Cloud Mountain minnow, Reeves' turtle, hagfish, sea lilies, and some arthropod species, using degenerate polymerase chain reaction methods. We also identified two R2 elements from the public genomic sequence database of the bloodfluke Schistosoma mansoni. One of the two bloodfluke R2 elements has two zinc-finger motifs at the N-terminus; this differs from other known R2 elements, which have one or three zinc-finger motifs. Phylogenetic analysis revealed that the whole phylogeny of R2 can be divided into 11 parts (subclades), in which the local R2 phylogeny and the corresponding host phylogeny are consistent. Divergence-versus-age analysis revealed that there is no reliable evidence for the horizontal transfer of R2 but supports the proposition that R2 has been vertically transferred since before the divergence of the deuterostomes and protostomes. The seeming inconsistency between the R2 phylogeny and the phylogeny of their hosts is due to the existence of paralogous lineages. The number of N-terminal zinc-finger motifs is consistent with the deep phylogeny of R2 and indicates that the common ancestor of R2 had three zinc-finger motifs at the N-terminus. This study revealed the long-term vertical inheritance and the ancient origin of sequence specificity of R2, both of which seem applicable to some other non-LTR retrotransposons.

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.

Integration of the 5' end of the retrotransposon, R2Bm, can be complemented by homologous recombination

R2Bm is a non-long-terminal-repeat (non-LTR) retrotransposon that was identified at a specific target site in the 28S rRNA genes of the silkworm, Bombyx mori. Although in vitro analysis has revealed that the 3' end of R2Bm is integrated into the target site by means of target-primed reverse transcription (TPRT), the mechanism of the 5' end integration is not well understood. We established a novel in vivo system to assay the insertion mechanism of R2Bm using a cultured cell line, C65, and a baculovirus, AcNPV, as host and vector, respectively. The 3' end of R2Bm integrated at the target site in the rRNA genes of C65 cells when an AcNPV containing both the full-length 3' UTR and the entire open reading frame (ORF) of R2Bm was introduced while the 5' end integration was incorrect. The 5' end of R2Bm was integrated, however, when the 28S gene sequence upstream of the R2Bm target site was added to the R2Bm sequence. Thus, in our assay, homologous sequences were likely essential for the successful integration of the entire R2Bm into the host cell genome. We also demonstrated that the failure to integrate caused by a frame-shifted ORF was rescued by co-infection with a helper virus that contained only the R2Bm ORF. This indicates that R2 retrotransposition can be complemented in trans. These findings suggest that the host's mechanism for DNA repair may be necessary for the integration of the 5' end of R2Bm and that R2Bm protein may only have the ability to integrate the 3' end of the element by TPRT.

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.

R2 retrotransposition on assembled nucleosomes depends on the translational position of the target site

R2 retrotransposons insert into the 28S rRNA genes of insects. Integration occurs by specific cleavage of the target site and utilization of the released DNA end to prime reverse transcription of the RNA transcript. Specificity of the protein to the target site is dependent upon nucleotide sequence recognition extending from 35 bp upstream to 15 bp downstream of the cleavage site. In this report, we show that sequence recognition and cleavage by the R2 protein can occur while the target site is assembled into nucleosomes. Reconstitution of DNA fragments containing the 28S gene sequence into a set of nucleosomes with different translational frames revealed that the R2 site adopted the same rotational orientation with respect to the histone octamer. Binding and cleavage by the R2 protein were most efficient when the upstream binding site for the R2 protein was near a nucleosome end. Interaction of the R2 protein with the nucleosome disrupted the histone:DNA contacts in the 50 bp region directly bound by R2, but did not modify the remainder of the nucleosome structure.

Distribution of rRNA introns in the three-dimensional structure of the ribosome

More than 1200 introns have been documented at over 150 unique sites in the small and large subunit ribosomal RNA genes (as of February 2002). Nearly all of these introns are assigned to one of four main types: group I, group II, archaeal and spliceosomal. This sequence information has been organized into a relational database that is accessible through the Comparative RNA Web Site (http://www.rna.icmb.utexas.edu/) While the rRNA introns are distributed across the entire tree of life, the majority of introns occur within a few phylogenetic groups. We analyzed the distributions of rRNA introns within the three-dimensional structures of the 30S and 50S ribosomes. Most sites in rRNA genes that contain introns contain only one type of intron. While the intron insertion sites occur at many different coordinates, the majority are clustered near conserved residues that form tRNA binding sites and the subunit interface. Contrary to our expectations, many of these positions are not accessible to solvent in the mature ribosome. The correlation between the frequency of intron insertions and proximity of the insertion site to functionally important residues suggests an association between intron evolution and rRNA function.

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.

SINEs and LINEs: the art of biting the hand that feeds you

SINEs and LINEs are short and long interspersed retrotransposable elements, respectively, that invade new genomic sites using RNA intermediates. SINEs and LINEs are found in almost all eukaryotes (although not in Saccharomyces cerevisiae) and together account for at least 34% of the human genome. The noncoding SINEs depend on reverse transcriptase and endonuclease functions encoded by partner LINEs. With the completion of many genome sequences, including our own, the database of SINEs and LINEs has taken a great leap forward. The new data pose new questions that can only be answered by detailed studies of the mechanism of retroposition. Current work ranges from the biochemistry of reverse transcription and integration invitro, target site selection in vivo, nucleocytoplasmic transport of the RNA and ribonucleoprotein intermediates, and mechanisms of genomic turnover. Two particularly exciting new ideas are that SINEs may help cells survive physiological stress, and that the evolution of SINEs and LINEs has been shaped by the forces of RNA interference. Taken together, these studies promise to explain the birth and death of SINEs and LINEs, and the contribution of these repetitive sequence families to the evolution of genomes.