Integrated cytogenetics and genomics analysis of transposable elements in the Nile tilapia, Oreochromis niloticus.

Integration of cytogenetics and genomics has become essential to a better view of architecture and function of genomes. Although the advances on genomic sequencing have contributed to study genes and genomes, the repetitive DNA fraction of the genome is still enigmatic and poorly understood. Among repeated DNAs, transposable elements (TEs) are major components of eukaryotic chromatin and their investigation has been hindered even after the availability of whole sequenced genomes. The cytogenetic mapping of TEs in chromosomes has proved to be of high value to integrate information from the micro level of nucleotide sequence to a cytological view of chromosomes. Different TEs have been cytogenetically mapped in cichlids; however, neither details about their genomic arrangement nor appropriated copy number are well defined by these approaches. The current study integrates TEs distribution in Nile tilapia Oreochromis niloticus genome based on cytogenetic and genomics/bioinformatics approach. The results showed that some elements are not randomly distributed and that some are genomic dependent on each other. Moreover, we found extensive overlap between genomics and cytogenetics data and that tandem duplication may be the major mechanism responsible for the genomic dynamics of TEs here analyzed. This paper provides insights in the genomic organization of TEs under an integrated view based on cytogenetics and genomics.

Identification of RNA binding motifs in the R2 retrotransposon-encoded reverse transcriptase.

R2 non-LTR retrotransposons insert at a specific site in the 28S rRNA genes of many animal phyla. R2 elements encode a single polypeptide with reverse transcriptase, endonuclease and nucleic acid binding domains. Integration involves separate cleavage of the two DNA strands at the target site and utilization of the released 3' ends to prime DNA synthesis. Critical to this integration is the ability of the protein to specifically bind 3' and 5' regions of the R2 RNA. In this report, alanine mutations in two conserved motifs N-terminal to the reverse transcriptase domain were generated and shown to result in proteins that retained the ability to cleave the first strand of the DNA target, to reverse transcribe RNA from an annealed primer and to displace annealed RNA when using DNA as a template. However, the mutant proteins had greatly reduced ability to bind 3' and 5' RNA in mobility shift assays, use the DNA target to prime reverse transcription and conduct second-strand DNA cleavage. These motifs thus appear to participate in all activities of the R2 protein known to require specific RNA binding. The similarity of these R2 RNA binding motifs to those of telomerase and group II introns is discussed.

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.

Y chromosome mediates ribosomal DNA silencing and modulates the chromatin state in Drosophila.

Although the Drosophila Y chromosome is degenerated, heterochromatic, and contains few genes, increasing evidence suggests that it plays an important role in regulating the expression of numerous autosomal and X-linked genes. Here we use 15 Y chromosomes originating from a single founder 550 generations ago to study the role of the Y chromosome in regulating rRNA gene transcription, position-effect variegation (PEV), and the link among rDNA copy number, global gene expression, and chromatin regulation. Based on patterns of rRNA gene transcription indicated by transcription of the retrotransposon R2 that specifically inserts into the 28S rRNA gene, we show that X-linked rDNA is silenced in males. The silencing of X-linked rDNA expression by the Y chromosome is consistent across populations and independent of genetic background. These Y chromosomes also vary more than threefold in rDNA locus size and cause dramatically different levels of PEV suppression. The degree of suppression is negatively associated with the number and fraction of rDNA units without transposon insertions, but not with total rDNA locus size. Gene expression profiling revealed hundreds of differentially expressed genes among these Y chromosome introgression lines, as well as a divergent global gene expression pattern between the low-PEV and high-PEV flies. Our findings suggest that the Y chromosome is involved in diverse phenomena related to transcriptional regulation including X-linked rDNA silencing and suppression of PEV phenotype. These results further expand our understanding of the role of the Y chromosome in modulating global gene expression, and suggest a link with modifications of the chromatin state.

Heterochromatin formation promotes longevity and represses ribosomal RNA synthesis.

Organismal aging is influenced by a multitude of intrinsic and extrinsic factors, and heterochromatin loss has been proposed to be one of the causes of aging. However, the role of heterochromatin in animal aging has been controversial. Here we show that heterochromatin formation prolongs lifespan and controls ribosomal RNA synthesis in Drosophila. Animals with decreased heterochromatin levels exhibit a dramatic shortening of lifespan, whereas increasing heterochromatin prolongs lifespan. The changes in lifespan are associated with changes in muscle integrity. Furthermore, we show that heterochromatin levels decrease with normal aging and that heterochromatin formation is essential for silencing rRNA transcription. Loss of epigenetic silencing and loss of stability of the rDNA locus have previously been implicated in aging of yeast. Taken together, these results suggest that epigenetic preservation of genome stability, especially at the rDNA locus, and repression of unnecessary rRNA synthesis, might be an evolutionarily conserved mechanism for prolonging lifespan.

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.

The R2 retrotransposon RNA families.

Analysis of the R2 retrotransposons from multiple silkmoth and fruitfly species have revealed three segments that contain conserved RNA secondary structures. These conserved structures play important roles in the propagation of the R2 element, including R2 RNA processing and transposon integration into the host genome as well as a likely role in translation. Two of the structured regions comprise protein binding sites: one is located in the 3' UTR and the other is in the 5' UTR close to the putative start of the R2 open reading frame (ORF). The 3' structure was deduced from chemical mapping and sequence comparison. The 5' structure was determined using a combination of chemical mapping, oligonucleotide binding, NMR and sequence analysis and contains an unusual pseudoknot structure. The third structure occurs at the 5' end of the R2 RNA and is responsible for self-cleavage of the 5' end of the element from a 28S ribosomal RNA co-transcript. A structure for this fragment was proposed based on motif searching and sequence comparison. There is remarkable similarity in sequence and structure to the hepatitis delta virus (HDV) ribozyme. Seed alignments for the 5' structure and the R2 ribozyme, containing representative sequences and consensus structures, have been submitted to the Rfam database.

Isoenergetic penta- and hexanucleotide microarray probing and chemical mapping provide a secondary structure model for an RNA element orchestrating R2 retrotransposon protein function.

LNA (locked nucleic acids, i.e. oligonucleotides with a methyl bridge between the 2' oxygen and 4' carbon of ribose) and 2,6-diaminopurine were incorporated into 2'-O-methyl RNA pentamer and hexamer probes to make a microarray that binds unpaired RNA approximately isoenergetically. That is, binding is roughly independent of target sequence if target is unfolded. The isoenergetic binding and short probe length simplify interpretation of binding to a structured RNA to provide insight into target RNA secondary structure. Microarray binding and chemical mapping were used to probe the secondary structure of a 323 nt segment of the 5' coding region of the R2 retrotransposon from Bombyx mori (R2Bm 5' RNA). This R2Bm 5' RNA orchestrates functioning of the R2 protein responsible for cleaving the second strand of DNA during insertion of the R2 sequence into the genome. The experimental results were used as constraints in a free energy minimization algorithm to provide an initial model for the secondary structure of the R2Bm 5' RNA.

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.

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.

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 diversity of retrotransposons and the properties of their reverse transcriptases.

A number of abundant mobile genetic elements called retrotransposons reverse transcribe RNA to generate DNA for insertion into eukaryotic genomes. Four major classes of retrotransposons are described here. First, the long-terminal-repeat (LTR) retrotransposons have similar structures and mechanisms to those of the vertebrate retroviruses. Genes that may enable these retrotransposons to leave a cell have been acquired by these elements in a number of animal and plant lineages. Second, the tyrosine recombinase retrotransposons are similar to the LTR retrotransposons except that they have substituted a recombinase for the integrase and recombine into the host chromosomes. Third, the non-LTR retrotransposons use a cleaved chromosomal target site generated by an encoded endonuclease to prime reverse transcription. Finally, the Penelope-like retrotransposons are not well understood but appear to also use cleaved DNA or the ends of chromosomes as primer for reverse transcription. Described in the second part of this review are the enzymatic properties of the reverse transcriptases (RTs) encoded by retrotransposons. The RTs of the LTR retrotransposons are highly divergent in sequence but have similar enzymatic activities to those of retroviruses. The RTs of the non-LTR retrotransposons have several unique properties reflecting their adaptation to a different mechanism of retrotransposition.

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.

Fruit flies and humans respond differently to retrotransposons.

Retrotransposable element insertions are 20 times more numerous per unit length of DNA in the large human genome compared to the small Drosophila genome. Whereas all Drosophila elements are subject to constant turnover (recent insertion and elimination by selection), this has not generally been the case for human retrotransposons. We suggest that a difference in recombination adopted by these organisms in response to the deleterious effects of interspersed repeated DNA can explain in part this fundamental difference between the evolutionary dynamics of fruit fly and human retrotransposons.

R2 target-primed reverse transcription: ordered cleavage and polymerization steps by protein subunits asymmetrically bound to the target DNA.

R2 elements are non-long terminal repeat retrotransposons that specifically insert into 28S rRNA genes of many animal groups. These elements encode a single protein with reverse transcriptase and endonuclease activities as well as specific DNA and RNA binding properties. In this report, gel shift experiments were conducted to investigate the stoichiometry of the DNA, RNA, and protein components of the integration reaction. The enzymatic functions associated with each of the protein complexes were also determined, and DNase I digests were used to footprint the protein onto the target DNA. Additionally, a short polypeptide containing the N-terminal putative DNA-binding motifs was footprinted on the DNA target site. These combined findings revealed that one protein subunit binds the R2 RNA template and the DNA 10 to 40 bp upstream of the insertion site. This subunit cleaves the first DNA strand and uses that cleavage to prime reverse transcription of the R2 RNA transcript. Another protein subunit(s) uses the N-terminal DNA binding motifs to bind to the 18 bp of target DNA downstream of the insertion site and is responsible for cleavage of the second DNA strand. A complete model for the R2 integration reaction is presented, which with minor modifications is adaptable to other non-LTR retrotransposons.

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.

Role of the Bombyx mori R2 element N-terminal domain in the target-primed reverse transcription (TPRT) reaction.

R2 is a site-specific non-long terminal repeat (non-LTR) retrotransposon encoding a single polypeptide with reverse transcriptase, DNA endonuclease and nucleic acid-binding domains. The current model of R2 retrotransposition involves an ordered series of cleavage and polymerization steps carried out by at least two R2 protein subunits, one bound upstream and one bound downstream of the integration site. The role in the retrotransposition reaction of two conserved DNA-binding motifs, a C2H2 zinc finger (ZF) and a Myb motif, located within the N-terminal domain of the protein are explored in this report. These motifs do not appear to play a role in RT or the ability of the protein to bind the R2 RNA transcript. Methylation and missing nucleoside interference-based DNA footprints using polypeptides to the N-terminal domain suggest the ZF and Myb motifs bind to regions -3 to -1 and +10 to +15 with reference to the insertion site. Mutations in these DNA sites or of the N-terminal protein domain blocked binding and the activity of the downstream subunit. Mutations of the protein domain also affected binding of the upstream subunit but not its function, suggesting the primary path to DNA target recognition by R2 involves both upstream and downstream subunits.

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.

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.