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

Copy number and sequence variation in rDNA of Daphnia pulex from natural populations: insights from whole-genome sequencing

Ribosomal DNA (rDNA) has a vital role in ribosome biogenesis as it contains the genes that encode ribosomal RNA (rRNA) separated by intergenic spacers (IGSs). The rRNA genes occur in hundreds to tens of thousands of copies per haploid genome in eukaryotes and are generally highly conserved with low variation within species. Due to the repetitive nature and large size of rDNA arrays, detecting intraindividual variation can be difficult. In this study, we use whole-genome sequences of 169 Daphnia pulex individuals from 10 natural populations to measure the copy number and sequence variation in rDNA. This revealed that variation in rDNA copy number between individuals spans an order of magnitude. We further observed a substantial level of sequence variation within individual genomes. As expected, single-nucleotide polymorphisms occurred in regions of lower functional constraint such as the IGS and expansion segments of the rRNA genes. The presence of strong linkage disequilibrium among variants facilitated identification of haplotypes within each population. Although there was evidence of recombination among haplotypes from different populations, it is insufficient to eliminate linkage disequilibrium within populations. Estimating copy number and haplotype diversity within individuals revealed that the level of intraindividual sequence variation is not strongly correlated with copy number. The observed patterns of variation highlight a complex evolutionary history of rDNA in D. pulex. Future research should explore the functional implications of rDNA copy number and sequence variation on organismal phenotypes.

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.

Conserved and divergent DNA recognition specificities and functions of R2 retrotransposon N-terminal domains

R2 non-long terminal repeat (non-LTR) retrotransposons are among the most extensively distributed mobile genetic elements in multicellular eukaryotes and show promise for applications in transgene supplementation of the human genome. They insert new gene copies into a conserved site in 28S ribosomal DNA with exquisite specificity. R2 clades are defined by the number of zinc fingers (ZFs) at the N terminus of the retrotransposon-encoded protein, postulated to additively confer DNA site specificity. Here, we illuminate general principles of DNA recognition by R2 N-terminal domains across and between clades, with extensive, specific recognition requiring only one or two compact domains. DNA-binding and protection assays demonstrate broadly shared as well as clade-specific DNA interactions. Gene insertion assays in cells identify the N-terminal domains sufficient for target-site insertion and reveal roles in second-strand cleavage or synthesis for clade-specific ZFs. Our results have implications for understanding evolutionary diversification of non-LTR retrotransposon insertion mechanisms and the design of retrotransposon-based gene therapies.

Substantial rDNA copy number reductions alter timing of development and produce variable tissue-specific phenotypes in C. elegans

The genes that encode ribosomal RNAs are present in several hundred copies in most eukaryotes. These vast arrays of repetitive ribosomal DNA (rDNA) have been implicated not just in ribosome biogenesis, but also aging, cancer, genome stability, and global gene expression. rDNA copy number is highly variable among and within species; this variability is thought to associate with traits relevant to human health and disease. Here we investigate the phenotypic consequences of multicellular life at the lower bounds of rDNA copy number. We use the model Caenorhabditis elegans, which has previously been found to complete embryogenesis using only maternally provided ribosomes. We find that individuals with rDNA copy number reduced to ∼5% of wild type are capable of further development with variable penetrance. Such individuals are sterile and exhibit severe morphological defects, particularly in post-embryonically dividing tissues such as germline and vulva. Developmental completion and fertility are supported by an rDNA copy number ∼10% of wild type, with substantially delayed development. Worms with rDNA copy number reduced to ∼33% of wild type display a subtle developmental timing defect that was absent in worms with higher copy numbers. Our results support the hypothesis that rDNA requirements vary across tissues and indicate that the minimum rDNA copy number for fertile adulthood is substantially less than the lowest naturally observed total copy number. The phenotype of individuals with severely reduced rDNA copy number is highly variable in penetrance and presentation, highlighting the need for continued investigation into the biological consequences of rDNA copy number variation.

Ribosomal DNA replication time coordinates completion of genome replication and anaphase in yeast

Timely completion of genome replication is a prerequisite for mitosis, genome integrity, and cell survival. A challenge to this timely completion comes from the need to replicate the hundreds of untranscribed copies of rDNA that organisms maintain in addition to the copies required for ribosome biogenesis. Replication of these rDNA arrays is relegated to late S phase despite their large size, repetitive nature, and essentiality. Here, we show that, in Saccharomyces cerevisiae, reducing the number of rDNA repeats leads to early rDNA replication, which results in delaying replication elsewhere in the genome. Moreover, cells with early-replicating rDNA arrays and delayed genome-wide replication aberrantly release the mitotic phosphatase Cdc14 from the nucleolus and enter anaphase prematurely. We propose that rDNA copy number determines the replication time of the rDNA locus and that the release of Cdc14 upon completion of rDNA replication is a signal for cell cycle progression.

A Spontaneous Inversion of the X Chromosome Heterochromatin Provides a Tool for Studying the Structure and Activity of the Nucleolus in Drosophila melanogaster

The pericentromeric heterochromatin is largely composed of repetitive sequences, making it difficult to analyze with standard molecular biological methods. At the same time, it carries many functional elements with poorly understood mechanisms of action. The search for new experimental models for the analysis of heterochromatin is an urgent task. In this work, we used the Rif1 mutation, which suppresses the underreplication of all types of repeated sequences, to analyze heterochromatin regions in polytene chromosomes of Drosophila melanogaster. In the Rif1 background, we discovered and described in detail a new inversion, In(1)19EHet, which arose on a chromosome already carrying the In(1)sc⁸ inversion and transferred a large part of X chromosome heterochromatin, including the nucleolar organizer to a new euchromatic environment. Using nanopore sequencing and FISH, we have identified the eu- and heterochromatin breakpoints of In(1)19EHet. The combination of the new inversion and the Rif1 mutation provides a promising tool for studies of X chromosome heterochromatin structure, nucleolar organization, and the nucleolar dominance phenomenon. In particular, we found that, with the complete polytenization of rDNA repeats, the nucleolus consists of a cloud-like structure corresponding to the classical nucleolus of polytene chromosomes, as well as an unusual intrachromosomal structure containing alternating transcriptionally active and inactive regions.

First discovered, long out of sight, finally visible: ribosomal DNA

With the advent of long-read sequencing, previously unresolvable genomic elements are being revisited in an effort to generate fully complete reference genomes. One such element is ribosomal DNA (rDNA), the highly conserved genomic region that encodes rRNAs. Genomic structure and content of the rDNA are variable in both prokarya and eukarya, posing interesting questions about the biology of rDNA. Here, we consider the types of variation observed in rDNA - including locus structure and number, copy number, and sequence variation - and their known phenotypic consequences. With recent advances in long-read sequencing technology, incorporating the full rDNA sequence into reference genomes is within reach. This knowledge will have important implications for understanding rDNA biology within the context of cell physiology and whole-organism phenotypes.

Impaired function of rDNA transcription initiation machinery leads to derepression of ribosomal genes with insertions of R2 retrotransposon

Eukaryotic genomes harbor hundreds of rRNA genes, many of which are transcriptionally silent. However, little is known about selective regulation of individual rDNA units. In Drosophila melanogaster, some rDNA repeats contain insertions of the R2 retrotransposon, which is capable to be transcribed only as part of pre-rRNA molecules. rDNA units with R2 insertions are usually inactivated, although R2 expression may be beneficial in cells with decreased rDNA copy number. Here we found that R2-inserted rDNA units are enriched with HP1a and H3K9me3 repressive mark, whereas disruption of the heterochromatin components slightly affects their silencing in ovarian germ cells. Surprisingly, we observed a dramatic upregulation of R2-inserted rRNA genes in ovaries lacking Udd (Under-developed) or other subunits (TAF1b and TAF1c-like) of the SL1-like complex, which is homologues to mammalian Selective factor 1 (SL1) involved in rDNA transcription initiation. Derepression of rRNA genes with R2 insertions was accompanied by a reduction of H3K9me3 and HP1a enrichment. We suggest that the impairment of the SL1-like complex affects a mechanism of selective activation of intact rDNA units which competes with heterochromatin formation. We also propose that R2 derepression may serve as an adaptive response to compromised rRNA synthesis.

Repression of interrupted and intact rDNA by the SUMO pathway in Drosophila melanogaster

Ribosomal RNAs (rRNAs) are essential components of the ribosome and are among the most abundant macromolecules in the cell. To ensure high rRNA level, eukaryotic genomes contain dozens to hundreds of rDNA genes, however, only a fraction of the rRNA genes seems to be active, while others are transcriptionally silent. We found that individual rDNA genes have high level of cell-to-cell heterogeneity in their expression in Drosophila melanogaster. Insertion of heterologous sequences into rDNA leads to repression associated with reduced expression in individual cells and decreased number of cells expressing rDNA with insertions. We found that SUMO (Small Ubiquitin-like Modifier) and SUMO ligase Ubc9 are required for efficient repression of interrupted rDNA units and variable expression of intact rDNA. Disruption of the SUMO pathway abolishes discrimination of interrupted and intact rDNAs and removes cell-to-cell heterogeneity leading to uniformly high expression of individual rDNA in single cells. Our results suggest that the SUMO pathway is responsible for both repression of interrupted units and control of intact rDNA expression.

R1 retrotransposons in the nucleolar organizers of Drosophila melanogaster are transcribed by RNA polymerase I upon heat shock

The ribosomal RNA genes (rDNA) of Drosophila melanogaster reside within centromere-proximal nucleolar organizers on both the X and Y chromosomes. Each locus contains between 200-300 tandem repeat rDNA units that encode 18S, 5.8S, 2S, and 28S ribosomal RNAs (rRNAs) necessary for ribosome biogenesis. In arthropods like Drosophila, about 60% of the rDNA genes have R1 and/or R2 retrotransposons inserted at specific sites within their 28S regions; these units likely fail to produce functional 28S rRNA. We showed earlier that R2 expression increases upon nucleolar stress caused by the loss of the ribosome assembly factor, Nucleolar Phosphoprotein of 140 kDa (Nopp140). Here we show that R1 expression is selectively induced by heat shock. Actinomycin D, but not α-amanitin, blocked R1 expression in S2 cells upon heat shock, indicating that R1 elements are transcribed by Pol I. A series of RT-PCRs established read-through transcription by Pol I from the 28S gene region into R1. Sequencing the RT-PCR products confirmed the 28S-R1 RNA junction and the expression of R1 elements within nucleolar rDNA rather than R1 elements known to reside in centromeric heterochromatin. Using a genome-wide precision run-on sequencing (PRO-seq) data set available at NCBI-GEO, we show that Pol I activity on R1 elements is negligible under normal non-heat shock conditions but increases upon heat shock. We propose that prior to heat shock Pol I pauses within the 5' end of R1 where we find a consensus "pause button", and that heat shock releases Pol I for read-through transcription farther into R1.

Transgenerational dynamics of rDNA copy number in Drosophila male germline stem cells

rDNA loci, composed of hundreds of tandemly duplicated arrays of rRNA genes, are known to be among the most unstable genetic elements due to their repetitive nature. rDNA instability underlies aging (replicative senescence) in yeast cells, however, its contribution to the aging of multicellular organisms is poorly understood. In this study, we investigate the dynamics of rDNA loci during aging in the Drosophila male germline stem cell (GSC) lineage, and show that rDNA copy number decreases during aging. Our study further reveals that this age-dependent decrease in rDNA copy number is heritable from generation to generation, yet GSCs in young animals that inherited reduced rDNA copy number are capable of recovering normal rDNA copy number. Based on these findings, we propose that rDNA loci are dynamic genetic elements, where rDNA copy number changes dynamically yet is maintained through a recovery mechanism in the germline.

Site-specific non-LTR retrotransposons

Although most of non-long terminal repeat (non-LTR) retrotransposons are incorporated in the host genome almost randomly, some non-LTR retrotransposons are incorporated into specific sequences within a target site. On the basis of structural and phylogenetic features, non-LTR retrotransposons are classified into two large groups, restriction enzyme-like endonuclease (RLE)-encoding elements and apurinic/apyrimidinic endonuclease (APE)-encoding elements. All clades of RLE-encoding non-LTR retrotransposons include site-specific elements. However, only two of more than 20 APE-encoding clades, Tx1 and R1, contain site-specific non-LTR elements. Site-specific non-LTR retrotransposons usually target within multi-copy RNA genes, such as rRNA gene (rDNA) clusters, or repetitive genomic sequences, such as telomeric repeats; this behavior may be a symbiotic strategy to reduce the damage to the host genome. Site- and sequence-specificity are variable even among closely related non-LTR elements and appeared to have changed during evolution. In the APE-encoding elements, the primary determinant of the sequence- specific integration is APE itself, which nicks one strand of the target DNA during the initiation of target primed reverse transcription (TPRT). However, other factors, such as interaction between mRNA and the target DNA, and access to the target region in the nuclei also affect the sequence-specificity. In contrast, in the RLE-encoding elements, DNA-binding motifs appear to affect their sequence-specificity, rather than the RLE domain itself. Highly specific integration properties of these site-specific non-LTR elements make them ideal alternative tools for sequence-specific gene delivery, particularly for therapeutic purposes in human diseases.

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.

Deletion of Drosophila Nopp140 induces subcellular ribosomopathies

The nucleolar and Cajal body phosphoprotein of 140 kDa (Nopp140) is considered a ribosome assembly factor, but its precise functions remain unknown. To approach this problem, we deleted the Nopp140 gene in Drosophila using FLP-FRT recombination. Genomic PCR, reverse transcriptase-PCR (RT-PCR), and immunofluorescence microscopy confirmed the loss of Nopp140, its messenger RNA (mRNA), and protein products from all tissues examined. Nopp140-/- larvae arrested in the second instar stage and most died within 8 days. While nucleoli appeared intact in Nopp140-/- cells, the C/D small nucleolar ribonucleoprotein (snoRNP) methyltransferase, fibrillarin, redistributed to the nucleoplasm in variable amounts depending on the cell type; RT-PCRs showed that 2'-O-methylation of ribosomal RNA (rRNA) in Nopp140-/- cells was reduced at select sites within both the 18S and 28S rRNAs. Ultrastructural analysis showed that Nopp140-/- cells were deficient in cytoplasmic ribosomes, but instead contained abnormal electron-dense cytoplasmic granules. Immunoblot analysis showed a loss of RpL34, and metabolic labeling showed a significant drop in protein translation, supporting the loss of functional ribosomes. Northern blots showed that pre-RNA cleavage pathways were generally unaffected by the loss of Nopp140, but that R2 retrotransposons that naturally reside within the 28S region of normally silent heterochromatic Drosophila ribosomal DNA (rDNA) genes were selectively expressed in Nopp140-/- larvae. Unlike copia elements and the related R1 retrotransposon, R2 expression appeared to be preferentially dependent on the loss of Nopp140 and not on environmental stresses. We believe the phenotypes described here define novel intracellular ribosomopathies resulting from the loss of Nopp140.

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.

In and out of the rRNA genes: characterization of Pokey elements in the sequenced Daphnia genome

Background

Only a few transposable elements are known to exhibit site-specific insertion patterns, including the well-studied R-element retrotransposons that insert into specific sites within the multigene rDNA. The only known rDNA-specific DNA transposon, Pokey (superfamily: piggyBac) is found in the freshwater microcrustacean, Daphnia pulex. Here, we present a genome-wide analysis of Pokey based on the recently completed whole genome sequencing project for D. pulex.

Results

Phylogenetic analysis of Pokey elements recovered from the genome sequence revealed the presence of four lineages corresponding to two divergent autonomous families and two related lineages of non-autonomous miniature inverted repeat transposable elements (MITEs). The MITEs are also found at the same 28S rRNA gene insertion site as the Pokey elements, and appear to have arisen as deletion derivatives of autonomous elements. Several copies of the full-length Pokey elements may be capable of producing an active transposase. Surprisingly, both families of Pokey possess a series of 200 bp repeats upstream of the transposase that is derived from the rDNA intergenic spacer (IGS). The IGS sequences within the Pokey elements appear to be evolving in concert with the rDNA units. Finally, analysis of the insertion sites of Pokey elements outside of rDNA showed a target preference for sites similar to the specific sequence that is targeted within rDNA.

Conclusions

Based on the target site preference of Pokey elements and the concerted evolution of a segment of the element with the rDNA unit, we propose an evolutionary path by which the ancestors of Pokey elements have invaded the rDNA niche. We discuss how specificity for the rDNA unit may have evolved and how this specificity has played a role in the long-term survival of these elements in the subgenus Daphnia.

Endonuclease domain of the Drosophila melanogaster R2 non-LTR retrotransposon and related retroelements: a new model for transposition

The molecular mechanisms of the transposition of non-long terminal repeat (non-LTR) retrotransposons are not well understood; the key questions of how the 3′-ends of cDNA copies integrate and how site-specific integration occurs remain unresolved. Integration depends on properties of the endonuclease (EN) domain of retrotransposons. Using the EN domain of the Drosophila R2 retrotransposon as a model for other, closely related non-LTR retrotransposons, we investigated the EN domain and found that it resembles archaeal Holliday-junction resolvases. We suggest that these non-LTR retrotransposons are co-transcribed with the host transcript. Combined with the proposed resolvase activity of the EN domain, this model yields a novel mechanism for site-specific retrotransposition within this class of retrotransposons, with resolution proceeding via a Holliday junction intermediate.

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

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

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

Long interspersed nuclear elements (LINEs) show tissue-specific, mosaic genome and methylation-unrestricted, widespread expression of noncoding RNAs in somatic tissues of the rat

We report strong somatic and germ line expression of LINE RNAs in eight different tissues of rat by using a novel ~2.8 kb genomic PstI-LINE DNA (P1-LINE) isolated from the rat brain. P1-LINE is present in a 93 kb LINE-SINE-cluster in sub-telomeric region of chromosome 12 (12p12) and as multiple truncated copies interspersed in all rat chromosomes. P1-LINEs occur as inverted repeats at multiple genomic loci in tissue-specific and mosaic patterns. P1-LINE RNAs are strongly expressed in brain, liver, lungs, heart, kidney, testes, spleen and thymus into large to small heterogeneous RNAs (~5.0 to 0.2 kb) in tissue-specific and dynamic patterns in individual rats. P1-LINE DNA is strongly methylated at CpG-dinucleotides in most genomic copies in all the tissues and weakly hypomethylated in few copies in some tissues. Small (700-75 nt) P1-LINE RNAs expressed in all tissues may be possible precursors for small regulatory RNAs (PIWI-interacting/piRNAs) bioinformatically derived from P1-LINE. The strong and dynamic expression of LINE RNAs from multiple chromosomal loci and the putative piRNAs in somatic tissues of rat under normal physiological conditions may define functional chromosomal domains marked by LINE RNAs as long noncoding RNAs (lncRNAs) unrestricted by DNA methylation. The tissue-specific, dynamic RNA expression and mosaic genomic distribution of LINEs representing a steady-state genomic flux of retrotransposon RNAs suggest for biological role of LINE RNAs as long ncRNAs and small piRNAs in mammalian tissues independent of their cellular fate for translation, reverse-transcription and retrotransposition. This may provide evolutionary advantages to LINEs and mammalian genomes.

Chromatin structure of ribosomal RNA genes in dipterans and its relationship to the location of nucleolar organizers

Nucleoli, nuclear organelles in which ribosomal RNA is synthesized and processed, emerge from nucleolar organizers (NORs) located in distinct chromosomal regions. In polytene nuclei of dipterans, nucleoli of some species can be observed under light microscopy exhibiting distinctive morphology: Drosophila and chironomid species display well-formed nucleoli in contrast to the fragmented and dispersed nucleoli seen in sciarid flies. The available data show no apparent relationship between nucleolar morphology and location of NORs in Diptera. The regulation of rRNA transcription involves controlling both the transcription rate per gene as well as the proportion of rRNA genes adopting a proper chromatin structure for transcription, since active and inactive rRNA gene copies coexist in NORs. Transcription units organized in nucleosomes and those lacking canonical nucleosomes can be analyzed by the method termed psoralen gel retarding assay (PGRA), allowing inferences on the ratio of active to inactive rRNA gene copies. In this work, possible connections between chromosomal location of NORs and proportion of active rRNA genes were studied in Drosophila melanogaster, and in chironomid and sciarid species. The data suggested a link between location of NORs and proportion of active rRNA genes since the copy number showing nucleosomal organization predominates when NORs are located in the pericentric heterochromatin. The results presented in this work are in agreement with previous data on the chromatin structure of rRNA genes from distantly related eukaryotes, as assessed by the PGRA.

Nucleolar dominance of the Y chromosome in Drosophila melanogaster

The rDNA genes are transcribed by RNA polymerase I to make structural RNAs for ribosomes. Hundreds of rDNA genes are typically arranged in an array that spans megabase pairs of DNA. These arrays are the major sites of transcription in growing cells, accounting for as much as 50% of RNA synthesis. The repetitive rDNA arrays are thought to use heterochromatic gene silencing as a mechanism for metabolic regulation, since repeated sequences nucleate heterochromatin formation in eukaryotes. Drosophila melanogaster carries an rDNA array on the X chromosome and on the Y chromosome, and genetic analysis has suggested that both are transcribed. However, using a chromatin-marking assay, we find that the entire X chromosome rDNA array is normally silenced in D. melanogaster males, while the Y chromosome rDNA array is dominant and expressed. This resembles "nucleolar dominance," a phenomenon that occurs in interspecific hybrids where an rDNA array from one parental species is silenced, and that from the other parent is preferentially transcribed. Interspecies nucleolar dominance is thought to result from incompatibilities between species-specific transcription factors and the rDNA promoters in the hybrid, but our results show that nucleolar dominance is a normal feature of rDNA regulation. Nucleolar dominance within D. melanogaster is only partially dependent on known components of heterochromatic gene silencing, implying that a distinctive chromatin regulatory system may act at rDNA genes. Finally, we isolate variant Y chromosomes that allow X chromosome array expression and suggest that the large-scale organization of rDNA arrays contribute to nucleolar dominance. This is the first example of allelic inactivation in D. melanogaster.

Drosophila melanogaster linker histone dH1 is required for transposon silencing and to preserve genome integrity

Histone H1 is an intrinsic component of chromatin, whose important contribution to chromatin structure is well-established in vitro. Little is known, however, about its functional roles in vivo. Here, we have addressed this question in Drosophila, a model system offering many advantages since it contains a single dH1 variant. For this purpose, RNAi was used to efficiently deplete dH1 in flies. Expression-profiling shows that dH1 depletion affects expression of a relatively small number of genes in a regional manner. Furthermore, depletion up-regulates inactive genes, preferentially those located in heterochromatin, while active euchromatic genes are down-regulated, suggesting that the contribution of dH1 to transcription regulation is mainly structural, organizing chromatin for proper gene-expression regulation. Up-regulated genes are remarkably enriched in transposons. In particular, R1/R2 retrotransposons, which specifically integrate in the rDNA locus, are strongly up-regulated. Actually, depletion increases expression of transposon-inserted rDNA copies, resulting in synthesis of aberrant rRNAs and enlarged nucleolus. Concomitantly, dH1-depleted cells accumulate extra-chromosomal rDNA, show increased γH2Av content, stop proliferation and activate apoptosis, indicating that depletion causes genome instability and affects proliferation. Finally, the contributions to maintenance of genome integrity and cell proliferation appear conserved in human hH1s, as their expression rescues proliferation of dH1-depleted cells.

Host RNAs, including transposons, are encapsidated by a eukaryotic single-stranded RNA virus

Next-generation sequencing is a valuable tool in our growing understanding of the genetic diversity of viral populations. Using this technology, we have investigated the RNA content of a purified nonenveloped single-stranded RNA virus, flock house virus (FHV). We have also investigated the RNA content of virus-like particles (VLPs) of FHV and the related Nudaurelia capensis omega virus. VLPs predominantly package ribosomal RNA and transcripts of their baculoviral expression vectors. In addition, we find that 5.3% of the packaged RNAs are transposable elements derived from the Sf21 genome. This observation may be important when considering the therapeutic use of VLPs. We find that authentic FHV virions also package a variety of host RNAs, accounting for 1% of the packaged nucleic acid. Significant quantities of host messenger RNAs, ribosomal RNA, noncoding RNAs, and transposable elements are readily detected. The packaging of these host RNAs elicits the possibility of horizontal gene transfer between eukaryotic hosts that share a viral pathogen. We conclude that the genetic content of nonenveloped RNA viruses is variable, not just by genome mutation, but also in the diversity of RNA transcripts that are packaged.

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.

Targeting novel sites: The N-terminal DNA binding domain of non-LTR retrotransposons is an adaptable module that is implicated in changing site specificities

Restriction-like endonuclease (RLE) bearing non-LTR retrotransposons are site-specific elements that integrate into the genome through target primed reverse transcription (TPRT). RLE-bearing elements have been used as a model system for investigating non-LTR retrotransposon integration. R2 elements target a specific site in the 28S rDNA gene. We previously demonstrated that the two major sub-classes of R2 (R2-A and R2-D) target the R2 insertion site in an opposing manner with regard to the pairing of known DNA binding domains and bound sequences-indicating that the A- and D-clades represent independently derived modes of targeting that site. Elements have been discovered that group phylogenetically with R2 but do not target the canonical R2 site. Here we extend our earlier studies to show that a separate R2-A clade element, which targets a site other than the canonical R2 site, does so by using the N-terminal zinc fingers and Myb motifs. We further extend our targeting studies beyond R2 clade elements by investigating the ability of the N-terminal zinc fingers from the nematode NeSL-1 element to target its integration site. Our data are consistent with the use of an N-terminal DNA binding domain as one of the major targeting determinants used by RLE-bearing non-LTR retrotransposons to secure a protein subunit near the insertion site. This N-terminal DNA binding domain can undergo modifications, allowing the element to target novel sites. The binding orientation of the N-terminal domain relative to the insertion site is quite variable.

Chromosome spreading of associated transposable elements and ribosomal DNA in the fish Erythrinus erythrinus. Implications for genome change and karyoevolution in fish

BACKGROUND

The fish, Erythrinus erythrinus, shows an interpopulation diversity, with four karyomorphs differing by chromosomal number, chromosomal morphology and heteromorphic sex chromosomes. Karyomorph A has a diploid number of 2n = 54 and does not have differentiated sex chromosomes. Karyomorph D has 2n = 52 chromosomes in females and 2n = 51 in males, and it is most likely derived from karyomorph A by the differentiation of a multiple X1X2Y sex chromosome system. In this study, we analyzed karyomorphs A and D by means of cytogenetic approaches to evaluate their evolutionary relationship.

RESULTS

Conspicuous differences in the distribution of the 5S rDNA and Rex3 non-LTR retrotransposon were found between the two karyomorphs, while no changes in the heterochromatin and 18S rDNA patterns were found between them. Rex3 was interstitially dispersed in most chromosomes. It had a compartmentalized distribution in the centromeric regions of only two acrocentric chromosomes in karyomorph A. In comparison, in karyomorph D, Rex3 was found in 22 acrocentric chromosomes in females and 21 in males. All 5S rDNA sites co-localized with Rex3, suggesting that these are associated in the genome. In addition, the origin of the large metacentric Y chromosome in karyomorph D by centric fusion was highlighted by the presence of internal telomeric sites and 5S rDNA/Rex3 sites on this chromosome.

CONCLUSION

We demonstrated that some repetitive DNAs (5S rDNA, Rex3 retroelement and (TTAGGG)n telomeric repeats) were crucial for the evolutionary divergence inside E. erythrinus. These elements were strongly associated with the karyomorphic evolution of this species. Our results indicate that chromosomal rearrangements and genomic modifications were significant events during the course of evolution of this fish. We detected centric fusions that were associated with the differentiation of the multiple sex chromosomes in karyomorph D, as well as a surprising increase of associated 5S rDNA/Rex3 loci, in contrast to karyomorph A. In this sense, E. erythrinus emerges as an excellent model system for better understanding the evolutionary mechanisms underlying the huge genome diversity in fish. This organism can also contribute to understanding vertebrate genome evolution as a whole.

R2 retrotransposons encode a self-cleaving ribozyme for processing from an rRNA cotranscript

The non-long terminal repeat (non-LTR) retrotransposon R2 is inserted into the 28S rRNA genes of many animals. Expression of the element appears to be by cotranscription with the rRNA gene unit. We show here that processing of the rRNA cotranscript at the 5' end of the R2 element in Drosophila simulans is rapid and utilizes an unexpected mechanism. Using RNA synthesized in vitro, the 5' untranslated region of R2 was shown capable of rapid and efficient self-cleavage of the 28S-R2 cotranscript. The 5' end generated in vitro by the R2 ribozyme was at the position identical to that found for in vivo R2 transcripts. The RNA segment corresponding to the R2 ribozyme could be folded into a double pseudoknot structure similar to that of the hepatitis delta virus (HDV) ribozyme. Remarkably, 21 of the nucleotide positions in and around the active site of the HDV ribozyme were identical in R2. R2 elements from other Drosophila species were also shown to encode HDV-like ribozymes capable of self-cleavage. Tracing their sequence evolution in the Drosophila lineage suggests that the extensive similarity of the R2 ribozyme from D. simulans to that of HDV was a result of convergent evolution, not common descent.

The R2 mobile element of Rhynchosciara americana: molecular, cytological and dynamic aspects

Ribosomal RNA genes are encoded by large units clustered (18S, 5S, and 28S) in the nucleolar organizer region in several organisms. Sometimes additional insertions are present in the coding region for the 28S rDNA. These insertions are specific non-long terminal repeat retrotransposons that have very restricted integration targets within the genome. The retrotransposon present in the genome of Rhynchosciara americana, RaR2, was isolated by the screening of a genomic library. Sequence analysis showed the presence of conserved regions, such as a reverse transcriptase domain and a zinc finger motif in the amino terminal region. The insertion site was highly conserved in R. americana and a phylogenetic analysis showed that this element belongs to the R2 clade. The chromosomal localization confirmed that the RaR2 mobile element was inserted into a specific site in the rDNA gene. The expression level of RaR2 in salivary glands during larval development was determined by quantitative RT-PCR, and the increase of relative expression in the 3P of the fourth instar larval could be related to intense gene activity characteristic of this stage. 5'-Truncated elements were identified in different DNA samples. Additionally, in three other Rhynchosciara species, the R2 element was present as a full-length element.

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.

Changes in chromatin structure correlate with transcriptional activity of nucleolar rDNA in polytene chromosomes

Ribosomal DNA genes (rDNA) are found in tandem arrays of hundreds of repeated genes, but only a fraction of these genes are actively transcribed. The regulatory mechanism controlling the transition between active and inactive rDNA in higher eukaryotes is vital for cell survival. Here, we show that the nucleolus from Drosophila salivary gland cells contains two levels of chromatin organization reflecting differences in transcriptional activity: Decondensed chromatin is highly occupied with TATA-box-binding protein (TBP), phosphorylated H3S10, and acetylated H3K14, suggesting that rDNA in decondensed nucleolar areas is actively transcribed. Condensed chromatin lacks TBP, phosphorylated H3S10, or acetylated H3K14 and is enriched in the rDNA retrotransposons R1 and R2. The data show that R1 and R2 retrotransposons are not actively transcribed in salivary glands and may lead to the epigenetic silencing of flanking rDNA genes and that the silencing mechanisms of these sequences might be partially independent of heterochromatin formation by methylation of histone H3 at lysine 9 and binding of heterochromatin protein 1.

Chromatin: linking structure and function in the nucleolus

The nucleolus is an informative model structure for studying how chromatin-regulated transcription relates to nuclear organisation. In this review, we describe how chromatin controls nucleolar structure through both the modulation of rDNA activity by convergently-evolved remodelling complexes and by direct effects upon rDNA packaging. This packaging not only regulates transcription but may also be important for suppressing internal recombination between tandem rDNA repeats. The identification of nucleolar histone chaperones and novel chromatin proteins by mass spectrometry suggests that structure-specific chromatin components remain to be characterised and may regulate the nucleolus in novel ways. However, it also suggests that there is considerable overlap between nucleolar and non-nucleolar-chromatin components. We conclude that a fuller understanding of nucleolar chromatin will be essential for understanding how gene organisation is linked with nuclear architecture.

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.

Genomic rearrangements and transcriptional analysis of the spliced leader-associated retrotransposon in RNA interference-deficient Trypanosoma brucei

The Trypanosoma brucei genome is colonized by the site-specific non-LTR retrotransposon SLACS, or spliced leader-associated conserved sequence, which integrates exclusively into the spliced leader (SL) RNA genes. Although there is evidence that the RNA interference (RNAi) machinery regulates SLACS transcript levels, we do not know whether RNAi deficiency affects the genomic stability of SLACS, nor do we understand the mechanism of SLACS transcription. Here, we report that prolonged culturing of RNAi-deficient T. brucei cells, but not wild-type cells, results in genomic rearrangements of SLACS. Furthermore, two populations of SLACS transcripts persist in RNAi-deficient cells: a full-length transcript of approximately 7 kb and a heterogeneous population of small SLACS transcripts ranging in size from 450 to 550 nt. We provide evidence that SLACS transcription initiates at the +1 of the interrupted SL RNA gene and proceeds into the 5' UTR and open reading frame 1 (ORF1). This transcription is carried out by an RNA polymerase with alpha-amanitin sensitivity reminiscent of SL RNA synthesis and is dependent on the SL RNA promoter. Additionally, we show that both sense and antisense small SLACS transcripts originate from ORF1 and that they are associated with proteins in vivo. We speculate that the small SLACS transcripts serve as substrates for the production of siRNAs to regulate SLACS expression.

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.

Expression of protein-coding genes embedded in ribosomal DNA

Ribosomal DNA (rDNA) is a specialised chromosomal location that is dedicated to high-level transcription of ribosomal RNA genes. Interestingly, rDNAs are frequently interrupted by parasitic elements, some of which carry protein genes. These are non-LTR retrotransposons and group II introns that encode reverse transcriptase-like genes, and group I introns and archaeal introns that encode homing endonuclease genes (HEGs). Although rDNA-embedded protein genes are widespread in nuclei, organelles and bacteria, there is surprisingly little information available on how these genes are expressed. Exceptions include a handful of HEGs from group I introns. Recent studies have revealed unusual and essential roles of group I and group I-like ribozymes in the endogenous expression of HEGs. Here we discuss general aspects of rDNA-embedded protein genes and focus on HEG expression from group I introns in the nucleolus.