Epigenetic regulation of retrotransposons within the nucleolus of Drosophila

Causes and Consequences of Varying Transposable Element Activity: An Evolutionary Perspective

Transposable elements (TEs) are genomic parasites found in nearly all eukaryotes, including humans. This evolutionary success of TEs is due to their replicative activity, involving insertion into new genomic locations. TE activity varies at multiple levels, from between taxa to within individuals. The rapidly accumulating evidence of the influence of TE activity on human health, as well as the rapid growth of new tools to study it, motivated an evaluation of what we know about TE activity thus far. Here, we discuss why TE activity varies, and the consequences of this variation, from an evolutionary perspective. By studying TE activity in nonhuman organisms in the context of evolutionary theories, we can shed light on the factors that affect TE activity. While the consequences of TE activity are usually deleterious, some have lasting evolutionary impacts by conferring benefits on the host or affecting other evolutionary processes.

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

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

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

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

Structural RNA components supervise the sequential DNA cleavage in R2 retrotransposon

Retroelements are the widespread jumping elements considered as major drivers for genome evolution, which can also be repurposed as gene-editing tools. Here, we determine the cryo-EM structures of eukaryotic R2 retrotransposon with ribosomal DNA target and regulatory RNAs. Combined with biochemical and sequencing analysis, we reveal two essential DNA regions, Drr and Dcr, required for recognition and cleavage. The association of 3' regulatory RNA with R2 protein accelerates the first-strand cleavage, blocks the second-strand cleavage, and initiates the reverse transcription starting from the 3'-tail. Removing 3' regulatory RNA by reverse transcription allows the association of 5' regulatory RNA and initiates the second-strand cleavage. Taken together, our work explains the DNA recognition and RNA supervised sequential retrotransposition mechanisms by R2 machinery, providing insights into the retrotransposon and application reprogramming.

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.

Evolution and dosage compensation of nucleolar organizing regions (NORs) mediated by mobile elements in turtles with female (ZZ/ZW) but not with male (XX/XY) heterogamety

Understanding the evolution and regulation of nucleolar organizing regions (NORs) is important to elucidate genome structure and function. This is because ribosomal gene (rDNA) copy number and activity mediate protein biosynthesis, stress response, ageing, disease, dosage compensation and genome stability. Here, we found contrasting dosage compensation of sex-linked NORs in turtles with male and female heterogamety. Most taxa examined exhibit homomorphic rRNA gene clusters in a single autosome pair (determined by 28S rDNA fluorescence in situ hybridization), whereas NORs are sex-linked in Apalone spinifera, Pelodiscus sinensis and Staurotypus triporcatus. Full-dosage compensation upregulates the male X-NOR (determined via silver staining-AgNOR) in Staurotypus (who lacks Y-NOR) compared with female X-AgNORs. In softshell Apalone and Pelodiscus, who share homologous ZZ/ZW micro-chromosomes, their enlarged W-NOR is partially active (due to 28S rDNA invasion by R2 retroelements), whereas their smaller Z-NOR is silent in females but active in both male-Zs (presumably because the W-NOR meets cellular demands and excessive NOR activity is costly). We hypothesize that R2 disruption favoured W enlargement to add intact 28S-units, perhaps facilitated by reduced recombination during sex chromosome evolution. The molecular basis of the potentially adaptive female Z-silencing is likely intricate and perhaps epigenetic, as non-ribosomal Z genes are active in Apalone females. Yet, Emydura maquarii exhibit identical heteromorphism in their autosomal NOR (R2 invaded 28S-units and the small-autosome NOR is silent), suggesting that the softshell turtle pattern can evolve independent of sex chromosome evolution. Our study illuminates the complex sex chromosome evolution and dosage compensation of non-model systems that challenges classic paradigms.

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.

The function of twister ribozyme variants in non-LTR retrotransposition in Schistosoma mansoni

The twister ribozyme is widely distributed over numerous organisms and is especially abundant in Schistosoma mansoni, but has no confirmed biological function. Of the 17 non-LTR retrotransposons known in S. mansoni, none have thus far been associated with ribozymes. Here we report the identification of novel twister variant (T-variant) ribozymes and their function in S. mansoni non-LTR retrotransposition. We show that T-variant ribozymes are located at the 5′ end of Perere-3 non-LTR retrotransposons in the S. mansoni genome. T-variant ribozymes were demonstrated to be catalytically active in vitro. In reporter constructs, T-variants were shown to cleave in vivo, and cleavage of T-variants was sufficient for the translation of downstream reporter genes. Our analysis shows that the T-variants and Perere-3 are transcribed together. Target site duplications (TSDs); markers of target-primed reverse transcription (TPRT) and footmarks of retrotransposition, are located adjacent to the T-variant cleavage site and suggest that T-variant cleavage has taken place inS. mansoni. Sequence heterogeneity in the TSDs indicates that Perere-3 retrotransposition is not site-specific. The TSD sequences contribute to the 5′ end of the terminal ribozyme helix (P1 stem). Based on these results we conclude that T-variants have a functional role in Perere-3 retrotransposition.

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

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

Novel ribozymes: discovery, catalytic mechanisms, and the quest to understand biological function

Small endonucleolytic ribozymes promote the self-cleavage of their own phosphodiester backbone at a specific linkage. The structures of and the reactions catalysed by members of individual families have been studied in great detail in the past decades. In recent years, bioinformatics studies have uncovered a considerable number of new examples of known catalytic RNA motifs. Importantly, entirely novel ribozyme classes were also discovered, for most of which both structural and biochemical information became rapidly available. However, for the majority of the new ribozymes, which are found in the genomes of a variety of species, a biological function remains elusive. Here, we concentrate on the different approaches to find catalytic RNA motifs in sequence databases. We summarize the emerging principles of RNA catalysis as observed for small endonucleolytic ribozymes. Finally, we address the biological functions of those ribozymes, where relevant information is available and common themes on their cellular activities are emerging. We conclude by speculating on the possibility that the identification and characterization of proteins that we hypothesize to be endogenously associated with catalytic RNA might help in answering the ever-present question of the biological function of the growing number of genomically encoded, small endonucleolytic ribozymes.

The processivity factor Pol32 mediates nuclear localization of DNA polymerase delta and prevents chromosomal fragile site formation in Drosophila development

The Pol32 protein is one of the universal subunits of DNA polymerase δ (Pol δ), which is responsible for genome replication in eukaryotic cells. Although the role of Pol32 in DNA repair has been well-characterized, its exact function in genome replication remains obscure as studies in single cell systems have not established an essential role for Pol32 in the process. Here we characterize Pol32 in the context of Drosophila melanogaster development. In the rapidly dividing embryonic cells, loss of Pol32 halts genome replication as it specifically disrupts Pol δ localization to the nucleus. This function of Pol32 in facilitating the nuclear import of Pol δ would be similar to that of accessory subunits of DNA polymerases from mammalian Herpes viruses. In post-embryonic cells, loss of Pol32 reveals mitotic fragile sites in the Drosophila genome, a defect more consistent with Pol32’s role as a polymerase processivity factor. Interestingly, these fragile sites do not favor repetitive sequences in heterochromatin, with the rDNA locus being a striking exception. Our study uncovers a possibly universal function for DNA polymerase ancillary factors and establishes a powerful system for the study of chromosomal fragile sites in a non-mammalian organism.

Cancer etiological studies suggest that the majority of pathological mutations occurred under near normal DNA replication conditions, emphasizing the importance of understanding replication regulation under non-lethal conditions. To gain such a better understanding, we investigated the function of Pol32, a conserved ancillary subunit of the essential DNA polymerase Delta complex, through the development of the fruit fly Drosophila. We uncovered a previously unappreciated function of Pol32 in regulating the nuclear import of the polymerase complex, and this function is developmentally regulated. By utilizing mutations in pol32 and other replication factors, we have started to define basic features of Chromosome Fragile Sites (CFS) in Drosophila somatic cells. CFS is a major source of genome instability associated with replication stresses, and has been an important topic of cancer biology. We discovered that CFS formation does not favor genomic regions with repetitive sequences except the highly transcribed locus encoding ribosomal RNA. Our work lays the groundwork for future studies using Drosophila as an alternative system to uncover the most fundamental features of CFS.

The Role of Sequence Duplication in Transcriptional Regulation and Genome Evolution

Sequence duplication is nowadays recognized as an important mechanism that underlies the evolution of eukaryote genomes, being indeed one of the most powerful strategies for the generation of adaptive diversity by modulating transcriptional activity. The evolutionary novelties simultaneously associated with sequence duplication and differential gene expression can be collectively referred to as duplication-mediated transcriptional regulation. In the last years, evidence has emerged supporting the idea that sequence duplication and functionalization represent important evolutionary strategies acting at the genome level, and both coding and non-coding sequences have been found to be targets of such events. Moreover, it has been proposed that deleterious effects of sequence duplication might be potentially silenced by endogenous cell machinery (i.e., RNA interference, epigenetic repressive marks, etc). Along these lines, our aim is to highlight the role of sequence duplication on transcriptional activity and the importance of both in genome evolution.

The peculiar genetics of the ribosomal DNA blurs the boundaries of transgenerational epigenetic inheritance

Our goal is to draw a line-hypothetical in its totality but experimentally supported at each individual step-connecting the ribosomal DNA and the phenomenon of transgenerational epigenetic inheritance of induced phenotypes. The reasonableness of this hypothesis is offset by its implication, that many (or most) (or all) of the cases of induced-and-inherited phenotypes that are seen to persist for generations are instead unmapped induced polymorphisms in the ribosomal DNA, and thus are the consequence of the peculiar and enduringly fascinating genetics of the highly transcribed repeat DNA structure at that locus.

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

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

Distribution of the DNA transposon family, Pokey in the Daphnia pulex species complex

BACKGROUND

The Pokey family of DNA transposons consists of two putatively autonomous groups, PokeyA and PokeyB, and two groups of Miniature Inverted-repeat Transposable Elements (MITEs), mPok1 and mPok2. This TE family is unusual as it inserts into a specific site in ribosomal (r)DNA, as well as other locations in Daphnia genomes. The goals of this study were to determine the distribution of the Pokey family in lineages of the Daphnia pulex species complex, and to test the hypothesis that unusally high PokeyA number in some isolates of Daphnia pulicaria is the result of recent transposition. To do this, we estimated the haploid number of Pokey, mPok, and rRNA genes in 45 isolates from five Daphnia lineages using quantitative PCR. We also cloned and sequenced partial copies of PokeyA from four isolates of D. pulicaria.

RESULTS

Haploid PokeyA and PokeyB number is generally less than 20 and tends to be higher outside rDNA in four lineages. Conversely, the number of both groups is much higher outside rDNA (~120) in D. arenata, and PokeyB is also somewhat higher inside rDNA. mPok1 was only detected in D. arenata. mPok2 occurs both outside (~30) and inside rDNA (~6) in D. arenata, but was rare (≤2) outside rDNA in the other four lineages. There is no correlation between Pokey and rRNA gene number (mean = 240 across lineages) in any lineage. Variation among cloned partial PokeyA sequences is significantly higher in isolates with high number compared to isolates with an average number.

CONCLUSIONS

The high Pokey number outside rDNA in D. arenata and inside rDNA in some D. pulicaria isolates is consistent with a recent increase in transposition rate. The D. pulicaria increase may have been triggered by insertion of PokeyA into a region of transcriptionally active rDNA. The expansion in D. arenata (thought to be of hybrid origin) may be a consequence of release from epigenetic repression following hybridization. Previous work found D. obtusa to be very different from the D. pulex complex; mean PokeyA is higher in rDNA (~75), rDNA array size is nearly twice as large (415), and the two are positively correlated. The predominance of Pokey in only one location could be explained by purifying selection against ectopic recombination between elements inside and outside rDNA.

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.

Transgenerational inheritance of diet-induced genome rearrangements in Drosophila

Ribosomal RNA gene (rDNA) copy number variation modulates heterochromatin formation and influences the expression of a large fraction of the Drosophila genome. This discovery, along with the link between rDNA, aging, and disease, high-lights the importance of understanding how natural rDNA copy number variation arises. Pursuing the relationship between rDNA expression and stability, we have discovered that increased dietary yeast concentration, emulating periods of dietary excess during life, results in somatic rDNA instability and copy number reduction. Modulation of Insulin/TOR signaling produces similar results, indicating a role for known nutrient sensing signaling pathways in this process. Furthermore, adults fed elevated dietary yeast concentrations produce offspring with fewer rDNA copies demonstrating that these effects also occur in the germline, and are transgenerationally heritable. This finding explains one source of natural rDNA copy number variation revealing a clear long-term consequence of diet.

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

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

Copy number of the transposon, Pokey, in rDNA is positively correlated with rDNA copy number in Daphnia obtuse [corrected]

Pokey is a class II DNA transposon that inserts into 28S ribosomal RNA (rRNA) genes and other genomic regions of species in the subgenus, Daphnia. Two divergent lineages, PokeyA and PokeyB have been identified. Recombination between misaligned rRNA genes changes their number and the number of Pokey elements. We used quantitative PCR (qPCR) to estimate rRNA gene and Pokey number in isolates from natural populations of Daphnia obtusa, and in clonally-propagated mutation accumulation lines (MAL) initiated from a single D. obtusa female. The change in direction and magnitude of Pokey and rRNA gene number did not show a consistent pattern across ∼ 87 generations in the MAL; however, Pokey and rRNA gene number changed in concert. PokeyA and 28S gene number were positively correlated in the isolates from both natural populations and the MAL. PokeyB number was much lower than PokeyA in both MAL and natural population isolates, and showed no correlation with 28S gene number. Preliminary analysis did not detect PokeyB outside rDNA in any isolates and detected only 0 to 4 copies of PokeyA outside rDNA indicating that Pokey may be primarily an rDNA element in D. obtusa. The recombination rate in this species is high and the average size of the rDNA locus is about twice as large as that in other Daphnia species such as D. pulicaria and D. pulex, which may have facilitated expansion of PokeyA to much higher numbers in D. obtusa rDNA than these other species.

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

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

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.

Reduced rDNA copy number does not affect "competitive" chromosome pairing in XYY males of Drosophila melanogaster

The ribosomal DNA (rDNA) arrays are causal agents in X-Y chromosome pairing in meiosis I of Drosophila males. Despite broad variation in X-linked and Y-linked rDNA copy number, polymorphisms in regulatory/spacer sequences between rRNA genes, and variance in copy number of interrupting R1 and R2 retrotransposable elements, there is little evidence that different rDNA arrays affect pairing efficacy. I investigated whether induced rDNA copy number polymorphisms affect chromosome pairing in a "competitive" situation in which complex pairing configurations were possible using males with XYY constitution. Using a common normal X chromosome, one of two different full-length Y chromosomes, and a third chromosome from a series of otherwise-isogenic rDNA deletions, I detected no differences in X-Y or Y-Y pairing or chromosome segregation frequencies that could not be attributed to random variation alone. This work was performed in the context of an undergraduate teaching program at Texas A&M University, and I discuss the pedagogical utility of this and other such experiments.

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.

Ribosomal DNA is active in different B chromosome variants of the grasshopper Eyprepocnemis plorans

B chromosomes are considered to be genetically inert elements. However, some of them are able to show nucleolus organizer region (NOR) activity, as detected by both cytological and molecular means. The grasshopper Eyprepocnemis plorans shows a B chromosome polymorphism characterized by the existence of many B variants. One of them, B24, shows NOR activity in about half of B-carrying males in the Torrox population. Molecular data have suggested the recent origin for B chromosomes in this species, and on this basis it would be expected that NOR activity was widespread among the different B variants. Here we test this hypothesis in four different B chromosome variants (B1, B2, B5, and B24) from 11 natural populations of the grasshopper E. plorans covering the south and east of the Iberian Peninsula plus the Balearic Islands. We used two different approaches: (1) the cytological observation of nucleoli attached to the distal region of the B chromosome (where the rDNA is located), and (2) the molecular detection of the rDNA transcripts carrying an adenine insertion characteristic of B chromosome ITS2 sequences. The results showed NOR expression not only for B24 but also for the B1 and B2 variants. However, the level of B-NOR expression in these latter variants, measured by the proportion of cells showing nucleoli attached to the B chromosomes, was much lower than that previously reported for B24. This suggests the possibility that structural or genetic background conditions are enhancing the expressivity of the rDNA in the B24 variant.

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

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

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

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

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

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.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Copy number variation of ribosomal DNA and Pokey transposons in natural populations of Daphnia

BACKGROUND

Despite their ubiquity and high diversity in eukaryotic genomes, DNA transposons are rarely encountered in ribosomal DNA (rDNA). In contrast, R-elements, a diverse group of non-LTR retrotransposons, specifically target rDNA. Pokey is a DNA transposon that targets a specific rDNA site, but also occurs in many other genomic locations, unlike R-elements. However, unlike most DNA transposons, Pokey has been a stable component of Daphnia genomes for over 100 million years. Here we use qPCR to estimate the number of 18S and 28S ribosomal RNA genes and Pokey elements in rDNA (rPokey), as well as other genomic locations (gPokey) in two species of Daphnia. Our goals are to estimate the correlation between (1) the number of 18S and 28S rRNA genes, (2) the number of 28S genes and rPokey, and (3) the number of rPokey and gPokey. In addition, we ask whether Pokey number and distribution in both genomic compartments are affected by differences in life history between D. pulex and D. pulicaria.

RESULTS

We found differences in 18S and 28S gene number within isolates that are too large to be explained by experimental variation. In general, Pokey number within isolates is modest (< 20), and most are gPokey. There is no correlation between the number of rRNA genes and rPokey, or between rPokey and gPokey. However, we identified three isolates with unusually high numbers of both rPokey and gPokey, which we infer is a consequence of recent transposition. We also detected other rDNA insertions (rInserts) that could be degraded Pokey elements, R- elements or the divergent PokeyB lineage recently detected in the Daphnia genome sequence. Unlike rPokey, rInserts are positively correlated with rRNA genes, suggesting that they are amplified by the same mechanisms that amplify rDNA units even though rPokey is not. Overall, Pokey frequency and distribution are similar in D. pulex and D. pulicaria suggesting that differences in life history have no impact on Pokey.

CONCLUSIONS

The possibility that many rDNA units do not contain a copy of both 18S and 28S genes suggests that rDNA is much more complicated than once thought, and warrants further study. In addition, the lack of correlation between rPokey, gPokey and rDNA unit numbers suggests that Pokey transposition rate is generally very low, and that recombination, in combination with natural selection, eliminates rPokey much faster than gPokey. Our results suggest that further research to determine the mechanisms by which Pokey has escaped complete inactivation by its host (the usual fate of DNA transposons), would provide important insights into transposon biology.

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.

Genetics: polymorphisms, epigenetics, and something in between

At its broadest sense, to say that a phenotype is epigenetic suggests that it occurs without changes in DNA sequence, yet is heritable through cell division and occasionally from one organismal generation to the next. Since gene regulatory changes are oftentimes in response to environmental stimuli and may be retained in descendent cells, there is a growing expectation that one's experiences may have consequence for subsequent generations and thus impact evolution by decoupling a selectable phenotype from its underlying heritable genotype. But the risk of this overbroad use of “epigenetic” is a conflation of genuine cases of heritable non-sequence genetic information with trivial modes of gene regulation. A look at the term “epigenetic” and some problems with its increasing prevalence argues for a more reserved and precise set of defining characteristics. Additionally, questions arising about how we define the “sequence independence” aspect of epigenetic inheritance suggest a form of genome evolution resulting from induced polymorphisms at repeated loci (e.g., the rDNA or heterochromatin).

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.

The CCCTC-binding factor (CTCF) of Drosophila contributes to the regulation of the ribosomal DNA and nucleolar stability

In the repeat array of ribosomal DNA (rDNA), only about half of the genes are actively transcribed while the others are silenced. In arthropods, transposable elements interrupt a subset of genes, often inactivating transcription of those genes. Little is known about the establishment or separation of juxtaposed active and inactive chromatin domains, or preferential inactivation of transposable element interrupted genes, despite identity in promoter sequences. CTCF is a sequence-specific DNA binding protein which is thought to act as a transcriptional repressor, block enhancer-promoter communication, and delimit juxtaposed domains of active and inactive chromatin; one or more of these activities might contribute to the regulation of this repeated gene cluster. In support of this hypothesis, we show that the Drosophila nucleolus contains CTCF, which is bound to transposable element sequences within the rDNA. Reduction in CTCF gene activity results in nucleolar fragmentation and reduced rDNA silencing, as does disruption of poly-ADP-ribosylation thought to be necessary for CTCF nucleolar localization. Our data establish a role for CTCF as a component necessary for proper control of transposable element-laden rDNA transcription and nucleolar stability.

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.

MERE1, a low-copy-number copia-type retroelement in Medicago truncatula active during tissue culture

We have identified an active Medicago truncatula copia-like retroelement called Medicago RetroElement1-1 (MERE1-1) as an insertion in the symbiotic NSP2 gene. MERE1-1 belongs to a low-copy-number family in the sequenced Medicago genome. These copies are highly related, but only three of them have a complete coding region and polymorphism exists between the long terminal repeats of these different copies. This retroelement family is present in all M. truncatula ecotypes tested but also in other legume species like Lotus japonicus. It is active only during tissue culture in both R108 and Jemalong Medicago accessions and inserts preferentially in genes.

Ribosomal DNA contributes to global chromatin regulation

The 35S ribosomal RNA genes (rDNA) are organized as repeated arrays in many organisms. Epigenetic regulation of transcription of the rRNA results in only a subset of copies being transcribed, making rDNA an important model for understanding epigenetic chromatin modification. We have created an allelic series of deletions within the rDNA array of the Drosophila Y chromosome that affect nucleolus size and morphology, but do not limit steady-state rRNA concentrations. These rDNA deletions result in reduced heterochromatin-induced gene silencing elsewhere in the genome, and the extent of the rDNA deletion correlates with the loss of silencing. Consistent with this, chromosomes isolated from strains mutated in genes required for proper heterochromatin formation have very small rDNA arrays, reinforcing the connection between heterochromatin and the rDNA. In wild-type cells, which undergo spontaneous natural rDNA loss, we observed the same correlation between loss of rDNA and loss of heterochromatin-induced silencing, showing that the volatility of rDNA arrays may epigenetically influence gene expression through normal development and differentiation. We propose that the rDNA contributes to a balance between heterochromatin and euchromatin in the nucleus, and alterations in rDNA--induced or natural--affect this balance.

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

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

Background

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

Results

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

Conclusions

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

Secondary structures for 5' regions of R2 retrotransposon RNAs reveal a novel conserved pseudoknot and regions that evolve under different constraints

Sequences from the 5' region of R2 retrotransposons of four species of silk moth are reported. In Bombyx mori, this region of the R2 messenger RNA contains a binding site for R2 protein and mediates interactions critical to R2 element insertion into the host genome. A model of secondary structure for a segment of this RNA is proposed on the basis of binding to oligonucleotide microarrays, chemical mapping, and comparative sequence analysis. Five conserved secondary structures are identified, including a novel pseudoknot. There is an apparent transition from an entirely RNA structure coding function in most of the 5' segment to a protein coding function near the 3' end. This suggests that local regions evolved under separate functional constraints (structural, coding, or both).

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.