Functional expression of a sequence-specific endonuclease encoded by the retrotransposon R2Bm

A fraction of the 28S ribosomal genes in certain insect species is interrupted by the insertion elements R1 and R2. These two elements from the silkworm Bombyx mori (R1Bm and R2Bm) are retrotransposons capable of transposing in a highly sequence-specific manner. We report here the functional expression in E. coli of the entire single open reading frame of R2Bm and show that it encodes a double-stranded endo-nuclease (integrase) that can specifically cleave the 28S gene at the R2 insertion site. The resulting cleavage is a 4 bp staggered 5' overhang. Deletion analysis of the 28S gene revealed that the DNA sequence required for specific cleavage is asymmetric with respect to the actual insertion (cleavage) site, with fewer than 10 bp required at one side and at least 24 bp at the other side of the site. A model is proposed based on these and previous data to account for the sequence-specific integration of the R2 retrotransposon.     

Evolution of R1 and R2 in the rDNA units of the genus Drosophila

R1 and R2 are non-long terminal repeat (non-LTR) retrotransposable elements that specifically insert in the 28S ribosomal RNA (rRNA) genes of insects. Using the Drosophila genus, which includes some of the best characterized insect taxa, we have conducted a number of studies on the evolution of these elements. We find that R1 and R2 are subject to the same recombinational forces that give rise to the concerted evolution of the rDNA units. The turnover of R1 and R2 elements can be readily documented in different strains of D. melanogaster using 5′ truncated elements as restriction-length polymorphisms. This turnover leads to uniform populations of elements with nucleotide sequence divergence of different copies averaging only 0.23% for the R2 and 0.47% for the R1 elements. Molecular phylogenetic analysis of elements from 16 different species of Drosophila suggests that these elements have been stable components of the rDNA locus for the 50–70 million year history of the Drosophila genus. Using changes at synonymous positions within the protein-encoding regions as estimates of the baseline substitution rate, it could be shown that R1 and R2 are evolving at rates similar to that of typical protein encoding genes provided corrections are made for the low codon bias of the elements. R1 and R2 are clearly well-adapted for their existence in the rDNA units of their host. This revised version was published online in August 2006 with corrections to the Cover Date.     

Rapid R2 retrotransposition leads to the loss of previously inserted copies via large deletions of the rDNA locus

R2 non-long terminal repeat retrotransposable elements insert specifically into the 28S rRNA genes of a wide range of animals. These elements maintain long-term stable relationships with the host genome. By scoring the variation present at the 5' ends of individual R2 copies, lines of Drosophila simulans have been identified with high rates of R2 retrotransposition. Comparing the R2 elements present in the parents with that of their progeny after 1 or 30 generations in this report revealed that retrotransposition rates were higher through the female germ line compared with the male germ line. In addition, most events in females occur late in germ line development. Surprisingly, the gain of new R2 insertions by retrotranspositions was counterbalanced by deletions of preexisting R2 insertions. These deletions occurred by the loss of large segments of the rDNA units that contained on average an estimated 15 R2 elements. When monitored over single generations, the rate of loss of preexisting elements was higher than the rate of new insertions. However, the chromosomes with the largest deletions appear to be eliminated from the population because the rates of R2 insertions and deletions after 30 generations were approximately equal. These findings suggest that high rates of R2 retrotransposition do not necessarily lead to dramatic increases in the level of R2 insertions in the rDNA locus but can lead to a more rapid turnover of rDNA units.     

Turnover of R1 (Type I) and R2 (Type Ii) Retrotransposable Elements in the Ribosomal DNA of Drosophila Melanogaster

R1 and R2 are distantly related non-long terminal repeat retrotransposable elements each of which inserts into a specific site in the 28S rRNA genes of most insects. We have analyzed aspects of R1 and R2 abundance and sequence variation in 27 geographical isolates of Drosophila melanogaster. The fraction of 28S rRNA genes containing these elements varied greatly between strains, 17-67% for R1 elements and 2-28% for R2 elements. The total percentage of the rDNA repeats inserted ranged from 32 to 77%. The fraction of the rDNA repeats that contained both of these elements suggested that R1 and R2 exhibit neither an inhibition of nor preference for insertion into a 28S gene already containing the other type of element. Based on the conservation of restriction sites in the elements of all strains, and sequence analysis of individual elements from three strains, nucleotide divergence is very low for R1 and R2 elements within or between strains (less than 0.6%). This sequence uniformity is the expected result of the forces of concerted evolution (unequal crossovers and gene conversion) which act on the rRNA genes themselves. Evidence for the role of retrotransposition in the turnover of R1 and R2 was obtained by using naturally occurring 5' length polymorphisms of the elements as markers for independent transposition events. The pattern of these different length 5' truncations of R1 and R2 was found to be diverse and unique to most strains analyzed. Because recombination can only, with time, amplify or eliminate those length variants already present, the diversity found in each strain suggests that retrotransposition has played a critical role in maintaining these elements in the rDNA repeats of D. melanogaster.     

Evolution of target specificity in R1 clade non-LTR retrotransposons

Although most non-long terminal repeat (non-LTR) retrotransposons are inserted throughout the host genome, many non-LTR elements in the R1 clade are inserted into specific sites within the target sequence. Four R1 clade families have distinct target specificity: R1 and RT insert into specific sites of 28S rDNA, and TRAS and SART insert into different sites within the (TTAGG)(n) telomeric repeats. To study the evolutionary history of target specificity of R1-clade retrotransposons, we have screened extensively novel representatives of the clade from various insects by in silico and degenerate polymerase chain reaction (PCR) cloning. We found four novel sequence-specific elements; Waldo (WaldoAg1, 2, and WaldoFs1) inserts into ACAY repeats, Mino (MinoAg1) into AC repeats, R6 into another specific site of the 28S rDNA, and R7 into a specific site of the 18S rDNA. In contrast, several elements (HOPE, WISHBm1, HidaAg1, NotoAg1, KagaAg1, Ha1Fs1) lost target sequence specificity, although some of them have preferred target sequences. Phylogenetic trees based on the RT and EN domains of each element showed that (1) three rDNA-specific elements, RT, R6, and R7, diverged from Waldo; (2) the elements having similar target sequences are phylogenetically related; and (3) the target specificity in the R1 clade was obtained once and thereafter altered and lost several times independently. These data indicate that the target specificity in R1 clade retroelements has changed during evolution and is more divergent than has been speculated so far.     

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.     

Site-specific ribosomal DNA insertion elements in Anopheles gambiae and A. arabiensis: nucleotide sequence of gene-element boundaries.

The nucleotide sequence of the junctions between the 28S ribosomal gene and site-specific insertion elements from two sibling mosquito species, Anopheles gambiae and A. arabiensis, is reported. In both species, elements insert at the same point within the 28S gene, but this site is 634 basepairs (bp) 3' of the R1 (Type I) insertion site in Drosophila melanogaster. The two mosquito elements each have poly A tails and a polyadenylation signal, but the extreme 3' and 5' ends show no other similarity to each other or to any other insertion element. In both mosquito species, identical target site duplications of 17 bp are generated. The sequence TNTCCCTNT found in this duplication is also found in the 14 bp target site duplications that flank R1 elements in D. melanogaster. Another sequence in this duplication, GGGATAACT, is very similar to the sequence GGGAGTAACT found in the 24 base sequence required by the Bombyx mori R2 endonuclease.     

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.     

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

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