Functional diversification of Dicer-like proteins and small RNAs required for genome sculpting

Functions and mechanisms of eukaryotic RNA-guided programmed DNA elimination

Many eukaryotic organisms, from ciliates to mammals, employ programmed DNA elimination during their postmeiotic reproduction. The process removes specific regions from the somatic DNA and has broad functions, including the irreversible silencing of genes, sex determination, and genome protection from transposable elements or integrating viruses. Multiple mechanisms have evolved that explain the sequence selectivity of the process. In some cases, the eliminated sequences lack centromeres and are flanked by conserved sequence motifs that are specifically recognized and cleaved by designated nucleases. Upon cleavage, all DNA fragments that lack centromeres are lost during the following mitosis. Alternatively, specific sequences can be destined for elimination by complementary small RNAs (sRNAs) as in some ciliates. These sRNAs enable a PIWI-mediated recruitment of chromatin remodelers, followed up by the precise positioning of a cleavage complex formed from a transposase like PiggyBac or Tc1. Here, we review the known molecular interplay of the cellular machinery that is involved in precise sRNA-guided DNA excision, and additionally, we highlight prominent knowledge gaps. We focus on the modes through which sRNAs enable the precise localization of the cleavage complex, and how the nuclease activity is controlled to prevent off-target cleavage. A mechanistic understanding of this process could enable the development of novel eukaryotic genome editing tools.

The linker region of a development-specific DNA polymerase X ensures efficient repair of programmed DNA double-strand breaks in Parameciumtetraurelia

During the sexual cycle, programmed genome rearrangement in Paramecium tetraurelia involves the non-homologous end joining (NHEJ) DNA repair pathway to eliminate specific germinal internal eliminated sequences (IESs) from the newly developing somatic nucleus. Besides the core NHEJ factors Ku70/80 and Xrcc4/Lig4, additional enzymes are required to process the 4-base 5'-protruding ends generated following DNA cleavage at IES boundaries, prior to their ligation. Here, we report that PolXa,b,c,d, four P. tetraurelia distant orthologs of the human Polλ DNA polymerase, are involved in the repair of IES excision junctions. During rearrangements, PolX-depleted cells accumulate genome-wide errors, such as unrepaired double-strand breaks, one-nucleotide deletions, and IES retention. Although all PolX paralogs can process DNA ends, two of them (PolXa&b) are induced during rearrangements and have evolved a specific linker sequence downstream of their BRCT domain, which provides them with tight nuclear anchoring properties. We show that PolXa accumulates in nuclear foci together with other NHEJ actors and the Dicer-like enzyme Dcl5, which is involved in the biogenesis of IES-specific small RNAs. We propose that these 'DNA repair foci' correspond to the sites where IES concatemers, a by-product of IES excision, are ligated together to produce the precursors of these small RNAs.

ISWI1 complex proteins facilitate developmental genome editing in Paramecium

One of the most extensive forms of natural genome editing occurs in ciliates, a group of microbial eukaryotes. Ciliate germline and somatic genomes are contained in distinct nuclei within the same cell. During the massive reorganization process of somatic genome development, ciliates eliminate tens of thousands of DNA sequences from a germline genome copy. Recently, we showed that the chromatin remodeler ISWI1 is required for somatic genome development in the ciliate Paramecium tetraurelia Here, we describe two high similarity paralogous proteins, ICOPa and ICOPb, essential for their genome editing. ICOPa and ICOPb are highly divergent from known proteins; the only domain detected showed distant homology with the WSD (WHIM2 + WHIM3) motif. We show that both ICOPa and ICOPb interact with the chromatin remodeler ISWI1. Upon ICOP knockdown, changes in alternative DNA excision boundaries and nucleosome densities are similar to those observed for ISWI1 knockdown. We thus propose that a complex comprising ISWI1 and either or both ICOPa and ICOPb are needed for Paramecium's precise genome editing.

The PIWI-interacting protein Gtsf1 controls the selective degradation of small RNAs in Paramecium

Ciliates undergo developmentally programmed genome elimination, in which small RNAs direct the removal of transposable elements (TEs) during the development of the somatic nucleus. Twenty-five nucleotide scanRNAs (scnRNAs) are produced from the entire germline genome and transported to the maternal somatic nucleus, where selection of scnRNAs corresponding to germline-specific sequences is thought to take place. Selected scnRNAs then guide the elimination of TEs in the developing somatic nucleus. How germline-specific scnRNAs are selected remains to be determined. Here, we provide important mechanistic insights into the scnRNA selection pathway by identifying a Paramecium homolog of Gtsf1 as essential for the selective degradation of scnRNAs corresponding to retained somatic sequences. Consistently, we also show that Gtsf1 is localized in the maternal somatic nucleus where it associates with the scnRNA-binding protein Ptiwi09. Furthermore, we demonstrate that the scnRNA selection process is critical for genome elimination. We propose that Gtsf1 is required for the coordinated degradation of Ptiwi09-scnRNA complexes that pair with target RNA via the ubiquitin pathway, similarly to the mechanism suggested for microRNA target-directed degradation in metazoans.

Heterochromatin-dependent transcription links the PRC2 complex to small RNA-mediated DNA elimination

Facultative heterochromatin is marked by the repressive histone modification H3K27me3 in eukaryotes. Deposited by the PRC2 complex, H3K27me3 is essential for regulating gene expression during development, and chromatin bearing this mark is generally considered transcriptionally inert. The PRC2 complex has also been linked to programmed DNA elimination during development in ciliates such as Paramecium. Due to a lack of mechanistic insight, a direct involvement has been questioned as most eliminated DNA segments in Paramecium are shorter than the size of a nucleosome. Here, we identify two sets of histone methylation readers essential for PRC2-mediated DNA elimination in Paramecium: Firefly1/2 and Mayfly1-4. The chromodomain proteins Firefly1/2 act in tight association with TFIIS4, a transcription elongation factor required for noncoding RNA transcription. These noncoding transcripts act as scaffolds for sequence-specific targeting by PIWI-bound sRNAs, resulting in local nucleosome depletion and DNA elimination. Our findings elucidate the molecular mechanism underlying the role of PRC2 in PIWI-mediated DNA elimination and suggest that its role in IES elimination may be to activate rather than repress transcription.

Genes and proteins expressed at different life cycle stages in the model protist Euplotes vannus revealed by both transcriptomic and proteomic approaches

Sexual reproduction first appeared in unicellular protists and has continued to be an essential biological process in almost all eukaryotes. Ciliated protists, which contain both germline and somatic genomes within a single cell, have evolved a special form of sexual reproduction called conjugation that involves mitosis, meiosis, fertilization, nuclear differentiation, genome rearrangement, and the development of unique cellular structures. The molecular basis and mechanisms of conjugation vary dramatically among ciliates, and many details of the process and its regulation are still largely unknown. In order to better comprehend these processes and mechanisms from an evolutionary perspective, this study provides the first comprehensive overview of the transcriptome and proteome profiles during the entire life cycle of the newly-established marine model ciliate Euplotes vannus. Transcriptome analyses from 14 life cycle stages (three vegetative stages and 11 sexual stages) revealed over 26,000 genes that are specifically expressed at different stages, many of which are related to DNA replication, transcription, translation, mitosis, meiosis, nuclear differentiation, and/or genome rearrangement. Quantitative proteomic analyses identified 338 proteins with homologs associated with conjugation and/or somatic nuclear development in other ciliates, including dicer-like proteins, Hsp90 proteins, RNA polymerase II and transcription elongation factors, ribosomal-associated proteins, and ubiquitin-related proteins. Four of these homologs belong to the PIWI family, each with different expression patterns identified and confirmed by RT-qPCR, which may function in small RNA-mediated genome rearrangement. Proteins involved in the nonhomologous end-joining pathway are induced early during meiosis and accumulate in the developing new somatic nucleus, where more than 80% of the germline sequences are eliminated from the somatic genome. A number of new candidate genes and proteins likely to play roles in conjugation and its related genome rearrangements have also been revealed. The gene expression profiles reported here will be valuable resources for further studies of the origin and evolution of sexual reproduction in this new model species.

GTSF1 is required for transposon silencing in the unicellular eukaryote Paramecium tetraurelia

The PIWI-interacting RNA (piRNA) pathway is crucial for transposon repression and the maintenance of genomic integrity. Gametocyte-specific factor 1 (GTSF1), a PIWI-associated protein indispensable for transposon repression, has been recently shown to potentiate the catalytic activity of PIWI in many metazoans. Whether the requirement of GTSF1 extends to PIWI proteins beyond metazoans is unknown. In this study, we identified a homolog of GTSF1 in the unicellular eukaryote Paramecium tetraurelia (PtGtsf1) and found that its role as a PIWI-cofactor is conserved. PtGtsf1 interacts with PIWI (Ptiwi09) and Polycomb Repressive Complex 2 and is essential for PIWI-dependent DNA elimination of transposons during sexual development. PtGtsf1 is crucial for the degradation of PIWI-bound small RNAs that recognize the organism's own genomic sequences. Without PtGtsf1, self-matching small RNAs are not degraded and results in an accumulation of H3K9me3 and H3K27me3, which may disturb transposon recognition. Our results demonstrate that the PIWI-GTSF1 interaction also exists in unicellular eukaryotes with a role in transposon silencing.

Nuclear dualism without extensive DNA elimination in the ciliate Loxodes magnus

Most eukaryotes have one nucleus and nuclear genome per cell. Ciliates have instead evolved distinct nuclei that coexist in each cell: a silent germline vs. transcriptionally active somatic nuclei. In the best-studied model species, both nuclei can divide asexually, but only germline nuclei undergo meiosis and karyogamy during sex. Thereafter, thousands of DNA segments, called internally eliminated sequences (IESs), are excised from copies of the germline genomes to produce the streamlined somatic genome. In Loxodes, however, somatic nuclei cannot divide but instead develop from germline copies even during asexual cell division, which would incur a huge overhead cost if genome editing was required. Here, we purified and sequenced both genomes in Loxodes magnus to see whether their nondividing somatic nuclei are associated with differences in genome architecture. Unlike in other ciliates studied to date, we did not find canonical germline-limited IESs, implying Loxodes does not extensively edit its genomes. Instead, both genomes appear large and equivalent, replete with retrotransposons and repetitive sequences, unlike the compact, gene-rich somatic genomes of other ciliates. Two other hallmarks of nuclear development in ciliates-domesticated DDE-family transposases and editing-associated small RNAs-were also not found. Thus, among the ciliates, Loxodes genomes most resemble those of conventional eukaryotes. Nonetheless, base modifications, histone marks, and nucleosome positioning of vegetative Loxodes nuclei are consistent with functional differentiation between actively transcribed somatic vs. inactive germline nuclei. Given their phylogenetic position, it is likely that editing was present in the ancestral ciliate but secondarily lost in the Loxodes lineage.

Two paralogous PHD finger proteins participate in natural genome editing in Paramecium tetraurelia

The unicellular eukaryote Paramecium tetraurelia contains functionally distinct nuclei: germline micronuclei (MICs) and a somatic macronucleus (MAC). During sex, the MIC genome is reorganized into a new MAC genome and the old MAC is lost. Almost 45,000 unique internal eliminated sequences (IESs) distributed throughout the genome require precise excision to guarantee a functional new MAC genome. Here, we characterize a pair of paralogous PHD finger proteins involved in DNA elimination. DevPF1, the early-expressed paralog, is present in only some of the gametic and post-zygotic nuclei during meiosis. Both DevPF1 and DevPF2 localize in the new developing MACs, where IES excision occurs. Upon DevPF2 knockdown (KD), long IESs are preferentially retained and late-expressed small RNAs decrease; no length preference for retained IESs was observed in DevPF1-KD and development-specific small RNAs were abolished. The expression of at least two genes from the new MAC with roles in genome reorganization seems to be influenced by DevPF1- and DevPF2-KD. Thus, both PHD fingers are crucial for new MAC genome development, with distinct functions, potentially via regulation of non-coding and coding transcription in the MICs and new MACs.

Small-RNA-guided histone modifications and somatic genome elimination in ciliates

Transposable elements and other repeats are repressed by small-RNA-guided histone modifications in fungi, plants and animals. The specificity of silencing is achieved through base-pairing of small RNAs corresponding to the these genomic loci to nascent noncoding RNAs, which allows the recruitment of histone methyltransferases that methylate histone H3 on lysine 9. Self-reinforcing feedback loops enhance small RNA production and ensure robust and heritable repression. In the unicellular ciliate Paramecium tetraurelia, small-RNA-guided histone modifications lead to the elimination of transposable elements and their remnants, a definitive form of repression. In this organism, germline and somatic functions are separated within two types of nuclei with different genomes. At each sexual cycle, development of the somatic genome is accompanied by the reproducible removal of approximately a third of the germline genome. Instead of recruiting a H3K9 methyltransferase, small RNAs corresponding to eliminated sequences tether Polycomb Repressive Complex 2, which in ciliates has the unique property of catalyzing both lysine 9 and lysine 27 trimethylation of histone H3. These histone modifications that are crucial for the elimination of transposable elements are thought to guide the endonuclease complex, which triggers double-strand breaks at these specific genomic loci. The comparison between ciliates and other eukaryotes underscores the importance of investigating small-RNAs-directed chromatin silencing in a diverse range of organisms. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > RNAi: Mechanisms of Action.

Structural maintenance of chromosomes (SMC) proteins are required for DNA elimination in Paramecium

Chromosome (SMC) proteins are a large family of ATPases that play important roles in the organization and dynamics of chromatin. They are central regulators of chromosome dynamics and the core component of condensin. DNA elimination during zygotic somatic genome development is a characteristic feature of ciliated protozoa such as Paramecium This process occurs after meiosis, mitosis, karyogamy, and another mitosis, which result in the formation of a new germline and somatic nuclei. The series of nuclear divisions implies an important role of SMC proteins in Paramecium sexual development. The relationship between DNA elimination and SMC has not yet been described. Here, we applied RNA interference, genome sequencing, mRNA sequencing, immunofluorescence, and mass spectrometry to investigate the roles of SMC components in DNA elimination. Our results show that SMC4-2 is required for genome rearrangement, whereas SMC4-1 is not. Functional diversification of SMC4 in Paramecium led to a formation of two paralogues where SMC4-2 acquired a novel, development-specific function and differs from SMC4-1. Moreover, our study suggests a competitive relationship between these two proteins.

Inter-generational nuclear crosstalk links the control of gene expression to programmed genome rearrangement during the Paramecium sexual cycle

Multinucleate cells are found in many eukaryotes, but how multiple nuclei coordinate their functions is still poorly understood. In the cytoplasm of the ciliate Paramecium tetraurelia, two micronuclei (MIC) serving sexual reproduction coexist with a somatic macronucleus (MAC) dedicated to gene expression. During sexual processes, the MAC is progressively destroyed while still ensuring transcription, and new MACs develop from copies of the zygotic MIC. Several gene clusters are successively induced and switched off before vegetative growth resumes. Concomitantly, programmed genome rearrangement (PGR) removes transposons and their relics from the new MACs. Development of the new MACs is controlled by the old MAC, since the latter expresses genes involved in PGR, including the PGM gene encoding the essential PiggyMac endonuclease that cleaves the ends of eliminated sequences. Using RNA deep sequencing and transcriptome analysis, we show that impairing PGR upregulates key known PGR genes, together with ∼600 other genes possibly also involved in PGR. Among these genes, 42% are no longer induced when no new MACs are formed, including 180 genes that are co-expressed with PGM under all tested conditions. We propose that bi-directional crosstalk between the two coexisting generations of MACs links gene expression to the progression of MAC development.

Developmental mRNA clearance by PIWI-bound endo-siRNAs in Paramecium

The clearance of untranslated mRNAs by Argonaute proteins is essential for embryonic development in metazoans. However, it is currently unknown whether similar processes exist in unicellular eukaryotes. The ciliate Paramecium tetraurelia harbors a vast array of PIWI-clade Argonautes involved in various small RNA (sRNA) pathways, many of which have not yet been investigated. Here, we investigate the function of a PIWI protein, Ptiwi08, whose expression is limited to a narrow time window during development, concomitant with the start of zygotic transcription. We show that Ptiwi08 acts in an endogenous small interfering RNA (endo-siRNA) pathway involved in the clearance of untranslated mRNAs. These endo-siRNAs are found in clusters that are strictly antisense to their target mRNAs and are a subset of siRNA-producing clusters (SRCs). Furthermore, the endo-siRNAs are 2'-O-methylated by Hen1 and require Dcr1 for their biogenesis. Our findings suggest that sRNA-mediated developmental mRNA clearance extends beyond metazoans and may be a more widespread mechanism than previously anticipated.

A potential role for SARS-CoV-2 small viral RNAs in targeting host microRNAs and modulating gene expression

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes coronavirus disease (COVID-19) in humans, with symptoms ranging from mild to severe, including fatality. The molecular mechanisms surrounding the effects of viral infection on the host RNA machinery remain poorly characterized. We used a comparative transcriptomics approach to investigate the effects of SARS-CoV-2 infection on the host mRNA and sRNA expression machinery in a human lung epithelial cell line (Calu-3) and an African green monkey kidney cell line (Vero-E6). Upon infection, we observed global changes in host gene expression and differential expression of dozens of host miRNAs, many with known links to viral infection and immune response. Additionally, we discovered an expanded landscape of more than a hundred SARS-CoV-2-derived small viral RNAs (svRNAs) predicted to interact with differentially expressed host mRNAs and miRNAs. svRNAs are derived from distinct regions of the viral genome and sequence signatures suggest they are produced by a non-canonical biogenesis pathway. 52 of the 67 svRNAs identified in Calu-3 cells are predicted to interact with differentially expressed miRNAs, with many svRNAs having multiple targets. Accordingly, we speculate that these svRNAs may play a role in SARS-CoV-2 propagation by modulating post-transcriptional gene regulation, and that methods for antagonizing them may have therapeutic value.

Chromatin remodeling is required for sRNA-guided DNA elimination in Paramecium

Small RNAs mediate the silencing of transposable elements and other genomic loci, increasing nucleosome density and preventing undesirable gene expression. The unicellular ciliate Paramecium is a model to study dynamic genome organization in eukaryotic cells, given its unique feature of nuclear dimorphism. Here, the formation of the somatic macronucleus during sexual reproduction requires eliminating thousands of transposon remnants (IESs) and transposable elements scattered throughout the germline micronuclear genome. The elimination process is guided by Piwi-associated small RNAs and leads to precise cleavage at IES boundaries. Here we show that IES recognition and precise excision are facilitated by recruiting ISWI1, a Paramecium homolog of the chromatin remodeler ISWI. ISWI1 knockdown substantially inhibits DNA elimination, quantitatively similar to development-specific sRNA gene knockdowns but with much greater aberrant IES excision at alternative boundaries. We also identify key development-specific sRNA biogenesis and transport proteins, Ptiwi01 and Ptiwi09, as ISWI1 cofactors in our co-immunoprecipitation studies. Nucleosome profiling indicates that increased nucleosome density correlates with the requirement for ISWI1 and other proteins necessary for IES excision. We propose that chromatin remodeling together with small RNAs is essential for efficient and precise DNA elimination in Paramecium.

A small RNA-guided PRC2 complex eliminates DNA as an extreme form of transposon silencing

In animal germlines, transposons are silenced at the transcriptional or post-transcriptional level to prevent deleterious expression. Ciliates employ a more direct approach by physically eliminating transposons from their soma, utilizing piRNAs to recognize transposons and imprecisely excise them. Ancient, mutated transposons often do not require piRNAs and are precisely eliminated. Here, we characterize the Polycomb Repressive Complex 2 (PRC2) in Paramecium and demonstrate its involvement in the removal of transposons and transposon-derived DNA. Our results reveal a striking difference between the elimination of new and ancient transposons at the chromatin level and show that the complex may be guided by Piwi-bound small RNAs (sRNAs). We propose that imprecise elimination in ciliates originates from an ancient transposon silencing mechanism, much like in plants and metazoans, through sRNAs, repressive methylation marks, and heterochromatin formation. However, it is taken a step further by eliminating DNA as an extreme form of transposon silencing.

Identification of novel, functional, long noncoding RNAs involved in programmed, large-scale genome rearrangements

Noncoding RNAs (ncRNAs) make up to ∼98% percent of the transcriptome of a given organism. In recent years, one relatively new class of ncRNAs, long noncoding RNAs (lncRNAs), were shown to be more than mere by-products of gene expression and regulation. The unicellular eukaryote Paramecium tetraurelia is a member of the ciliate phylum, an extremely heterogeneous group of organisms found in most bodies of water across the globe. A hallmark of ciliate genetics is nuclear dimorphism and programmed elimination of transposons and transposon-derived DNA elements, the latter of which is essential for the maintenance of the somatic genome. Paramecium and ciliates in general harbor a plethora of different ncRNA species, some of which drive the process of large-scale genome rearrangements, including DNA elimination, during sexual development. Here, we identify and validate the first known functional lncRNAs in ciliates to date. Using deep-sequencing and subsequent bioinformatic processing and experimental validation, we show that Paramecium expresses at least 15 lncRNAs. These candidates were predicted by a highly conservative pipeline, and informatic analyses hint at differential expression during development. Depletion of two lncRNAs, lnc1 and lnc15, resulted in clear phenotypes, decreased survival, morphological impairment, and a global effect on DNA elimination.

Early developmental, meiosis-specific proteins - Spo11, Msh4-1, and Msh5 - Affect subsequent genome reorganization in Paramecium tetraurelia

Developmental DNA elimination in Paramecium tetraurelia occurs through a trans-nuclear comparison of the genomes of two distinct types of nuclei: the germline micronucleus (MIC) and the somatic macronucleus (MAC). During sexual reproduction, which starts with meiosis of the germline nuclei, MIC-limited sequences including Internal Eliminated Sequences (IESs) and transposons are eliminated from the developing MAC in a process guided by noncoding RNAs (scnRNAs and iesRNAs). However, our current understanding of this mechanism is still very limited. Therefore, studying both genetic and epigenetic aspects of these processes is a crucial step to understand this phenomenon in more detail. Here, we describe the involvement of homologs of classical meiotic proteins, Spo11, Msh4-1, and Msh5 in this phenomenon. Based on our analyses, we propose that proper functioning of Spo11, Msh4-1, and Msh5 during Paramecium sexual reproduction are necessary for genome reorganization and viable progeny. Also, we show that double-strand breaks (DSBs) in DNA induced during meiosis by Spo11 are crucial for proper IESs excision. In summary, our investigations show that early sexual reproduction processes may significantly influence later somatic genome integrity.

Paramecium epigenetics in development and proliferation

The term epigenetics is used for any layer of genetic information aside from the DNA base-sequence information. Mammalian epigenetic research increased our understanding of chromatin dynamics in terms of cytosine methylation and histone modification during differentiation, aging, and disease. Instead, ciliate epigenetics focused more on small RNA-mediated effects. On the one hand, these do concern the transport of RNA from parental to daughter nuclei, representing a regulated transfer of epigenetic information across generations. On the other hand, studies of Paramecium, Tetrahymena, Oxytricha, and Stylonychia revealed an almost unique function of transgenerational RNA. Rather than solely controlling chromatin dynamics, they control sexual progeny's DNA content quantitatively and qualitatively. Thus epigenetics seems to control genetics, at least genetics of the vegetative macronucleus. This combination offers ciliates, in particular, an epigenetically controlled genetic variability. This review summarizes the epigenetic mechanisms that contribute to macronuclear heterogeneity and relates these to nuclear dimorphism. This system's adaptive and evolutionary possibilities raise the critical question of whether such a system is limited to unicellular organisms or binuclear cells. We discuss here the relevance of ciliate genetics and epigenetics to multicellular organisms.

The transient Spt4-Spt5 complex as an upstream regulator of non-coding RNAs during development

The Spt4-Spt5 complex is conserved and essential RNA polymerase elongation factor. To investigate the role of the Spt4-Spt5 complex in non-coding transcription during development, we used the unicellular model Paramecium tetraurelia. In this organism harboring both germline and somatic nuclei, massive transcription of the entire germline genome takes place during meiosis. This phenomenon starts a series of events mediated by different classes of non-coding RNAs that control developmentally programmed DNA elimination. We focused our study on Spt4, a small zinc-finger protein encoded in P. tetraurelia by two genes expressed constitutively and two genes expressed during meiosis. SPT4 genes are not essential in vegetative growth, but they are indispensable for sexual reproduction, even though genes from both expression families show functional redundancy. Silencing of the SPT4 genes resulted in the absence of double-stranded ncRNAs and reduced levels of scnRNAs - 25 nt-long sRNAs produced from these double-stranded precursors in the germline nucleus. Moreover, we observed that the presence of a germline-specific Spt4-Spt5m complex is necessary for transfer of the scnRNA-binding PIWI protein between the germline and somatic nucleus. Our study establishes that Spt4, together with Spt5m, is essential for expression of the germline genome and necessary for developmental genome rearrangements.

Developmental timing of programmed DNA elimination in Paramecium tetraurelia recapitulates germline transposon evolutionary dynamics

With its nuclear dualism, the ciliate Paramecium constitutes a unique model to study how host genomes cope with transposable elements (TEs). P. tetraurelia harbors two germline micronuclei (MICs) and a polyploid somatic macronucleus (MAC) that develops from one MIC at each sexual cycle. Throughout evolution, the MIC genome has been continuously colonized by TEs and related sequences that are removed from the somatic genome during MAC development. Whereas TE elimination is generally imprecise, excision of approximately 45,000 TE-derived internal eliminated sequences (IESs) is precise, allowing for functional gene assembly. Programmed DNA elimination is concomitant with genome amplification. It is guided by noncoding RNAs and repressive chromatin marks. A subset of IESs is excised independently of this epigenetic control, raising the question of how IESs are targeted for elimination. To gain insight into the determinants of IES excision, we established the developmental timing of DNA elimination genome-wide by combining fluorescence-assisted nuclear sorting with high-throughput sequencing. Essentially all IESs are excised within only one endoreplication round (32C to 64C), whereas TEs are eliminated at a later stage. We show that DNA elimination proceeds independently of replication. We defined four IES classes according to excision timing. The earliest excised IESs tend to be independent of epigenetic factors, display strong sequence signals at their ends, and originate from the most ancient integration events. We conclude that old IESs have been optimized during evolution for early and accurate excision by acquiring stronger sequence determinants and escaping epigenetic control.

One Cell, Two Gears: Extensive Somatic Genome Plasticity Accompanies High Germline Genome Stability in Paramecium

Mutation accumulation (MA) experiments are conventionally employed to study spontaneous germline mutations. However, MA experiments can also shed light on somatic genome plasticity in a habitual and genetic drift-maximizing environment. Here, we revisit an MA experiment that uncovered extraordinary germline genome stability in Paramecium tetraurelia, a single-celled eukaryote with nuclear dimorphism. Our re-examination of isogenic P. tetraurelia MA lines propagated in nutrient-rich medium for >40 sexual cycles reveals that their polyploid somatic genome accrued hundreds of intervening DNA segments (IESs), which are normally eliminated during germline-soma differentiation. These IESs frequently occupy a fraction of the somatic DNA copies of a given locus, producing IES excision/retention polymorphisms, and preferentially fall into a class of epigenetically controlled sequences. Relative to control lines, retained IESs are flanked by stronger cis-acting signals and interrupt an excess of highly expressed coding exons. These findings suggest that P. tetraurelia’s elevated germline DNA replication fidelity is associated with pervasive somatic genome plasticity. They show that MA regimes are powerful tools for investigating the role that developmental plasticity, somatic mutations, and epimutations have in ecology and evolution.

Two Piwis with Ago-like functions silence somatic genes at the chromatin level

Most sRNA biogenesis mechanisms involve either RNAse III cleavage or ping-pong amplification by different Piwi proteins harbouring slicer activity. Here, we follow the question why the mechanism of transgene-induced silencing in the ciliate Paramecium needs both Dicer activity and two Ptiwi proteins. This pathway involves primary siRNAs produced from non-translatable transgenes and secondary siRNAs from targeted endogenous loci. Our data does not indicate any signatures from ping-pong amplification but Dicer cleavage of long dsRNA. Ptiwi13 and 14 prefer different sub-cellular localizations and different preferences for primary and secondary siRNAs but do not load them mutually exclusive. Both Piwis enrich for antisense RNAs and show a general preference for uridine-rich sRNAs along the entire sRNA length. In addition, Ptiwi14-loaded siRNAs show a 5´-U signature. Our data indicates both Ptiwis and 2´-O-methylation contributing to strand selection of Dicer cleaved siRNAs. This unexpected function of the two distinct vegetative Piwis extends the increasing knowledge of the diversity of Piwi functions in diverse silencing pathways. We describe an unusual mode of action of Piwi proteins extending not only the great variety of Piwi-associated RNAi pathways but moreover raising the question whether this could have been the primordial one.

Emergent RNA-RNA interactions can promote stability in a facultative phototrophic endosymbiosis

Stable endosymbiosis between eukaryotic microbes has driven the evolution of further cellular complexity. Yet the mechanisms that can act to stabilize an emergent eukaryote–eukaryote endosymbiosis are unclear. Using the model facultative endosymbiotic system, Paramecium bursaria, we demonstrate that endosymbiont–host RNA–RNA interactions can drive a cost to host growth upon endosymbiont digestion. These RNA–RNA interactions are facilitated by the host RNA-interference system. For endosymbiont messenger RNA sharing a high level of sequence identity with host transcripts, this process can result in host gene knockdown. We propose that these endosymbiont–host RNA–RNA interactions—“RNA-interference collisions”—represent an emergent mechanism to sanction the host for breakdown of the endosymbiosis, promoting the stability of the facultative endosymbiotic interaction.

Eukaryote–eukaryote endosymbiosis was responsible for the spread of chloroplast (plastid) organelles. Stability is required for the metabolic and genetic integration that drives the establishment of new organelles, yet the mechanisms that act to stabilize emergent endosymbioses—between two fundamentally selfish biological organisms—are unclear. Theory suggests that enforcement mechanisms, which punish misbehavior, may act to stabilize such interactions by resolving conflict. However, how such mechanisms can emerge in a facultative endosymbiosis has yet to be explored. Here, we propose that endosymbiont–host RNA–RNA interactions, arising from digestion of the endosymbiont population, can result in a cost to host growth for breakdown of the endosymbiosis. Using the model facultative endosymbiosis between Paramecium bursaria and Chlorella spp., we demonstrate that this mechanism is dependent on the host RNA-interference (RNAi) system. We reveal through small RNA (sRNA) sequencing that endosymbiont-derived messenger RNA (mRNA) released upon endosymbiont digestion can be processed by the host RNAi system into 23-nt sRNA. We predict multiple regions of shared sequence identity between endosymbiont and host mRNA, and demonstrate through delivery of synthetic endosymbiont sRNA that exposure to these regions can knock down expression of complementary host genes, resulting in a cost to host growth. This process of host gene knockdown in response to endosymbiont-derived RNA processing by host RNAi factors, which we term “RNAi collisions,” represents a mechanism that can promote stability in a facultative eukaryote–eukaryote endosymbiosis. Specifically, by imposing a cost for breakdown of the endosymbiosis, endosymbiont–host RNA–RNA interactions may drive maintenance of the symbiosis across fluctuating ecological conditions.

Massive colonization of protein-coding exons by selfish genetic elements in Paramecium germline genomes

Ciliates are unicellular eukaryotes with both a germline genome and a somatic genome in the same cytoplasm. The somatic macronucleus (MAC), responsible for gene expression, is not sexually transmitted but develops from a copy of the germline micronucleus (MIC) at each sexual generation. In the MIC genome of Paramecium tetraurelia, genes are interrupted by tens of thousands of unique intervening sequences called internal eliminated sequences (IESs), which have to be precisely excised during the development of the new MAC to restore functional genes. To understand the evolutionary origin of this peculiar genomic architecture, we sequenced the MIC genomes of 9 Paramecium species (from approximately 100 Mb in Paramecium aurelia species to >1.5 Gb in Paramecium caudatum). We detected several waves of IES gains, both in ancestral and in more recent lineages. While the vast majority of IESs are single copy in present-day genomes, we identified several families of mobile IESs, including nonautonomous elements acquired via horizontal transfer, which generated tens to thousands of new copies. These observations provide the first direct evidence that transposable elements can account for the massive proliferation of IESs in Paramecium. The comparison of IESs of different evolutionary ages indicates that, over time, IESs shorten and diverge rapidly in sequence while they acquire features that allow them to be more efficiently excised. We nevertheless identified rare cases of IESs that are under strong purifying selection across the aurelia clade. The cases examined contain or overlap cellular genes that are inactivated by excision during development, suggesting conserved regulatory mechanisms. Similar to the evolution of introns in eukaryotes, the evolution of Paramecium IESs highlights the major role played by selfish genetic elements in shaping the complexity of genome architecture and gene expression.

A comparative genomics study of nine Paramecium species reveals successful invasion of genes by transposable elements in their germline genomes, showing that the internal eliminated sequences (IESs) followed an evolutionary trajectory remarkably similar to that of spliceosomal introns.

Characterization of the RNA-interference pathway as a tool for reverse genetic analysis in the nascent phototrophic endosymbiosis, Paramecium bursaria

Endosymbiosis was fundamental for the evolution of eukaryotic complexity. Endosymbiotic interactions can be dissected through forward- and reverse-genetic experiments, such as RNA-interference (RNAi). However, distinguishing small (s)RNA pathways in a eukaryote-eukaryote endosymbiotic interaction is challenging. Here, we investigate the repertoire of RNAi pathway protein-encoding genes in the model nascent endosymbiotic system, Paramecium bursaria-Chlorella spp. Using comparative genomics and transcriptomics supported by phylogenetics, we identify essential proteome components of the small interfering (si)RNA, scan (scn)RNA and internal eliminated sequence (ies)RNA pathways. Our analyses reveal that copies of these components have been retained throughout successive whole genome duplication (WGD) events in the Paramecium clade. We validate feeding-induced siRNA-based RNAi in P. bursaria via knock-down of the splicing factor, u2af1, which we show to be crucial to host growth. Finally, using simultaneous knock-down 'paradox' controls to rescue the effect of u2af1 knock-down, we demonstrate that feeding-induced RNAi in P. bursaria is dependent upon a core pathway of host-encoded Dcr1, Piwi and Pds1 components. Our experiments confirm the presence of a functional, host-derived RNAi pathway in P. bursaria that generates 23-nt siRNA, validating the use of the P. bursaria-Chlorella spp. system to investigate the genetic basis of a nascent endosymbiosis.

Cross-Generational Effects and Non-random Developmental Response to Temperature Variation in Paramecium

Unicellular organisms such as ciliates are largely neglected in research on adaptive developmental plasticity, although their nuclear dualism offers ideal circumstances to study development outside an embryonic context. Here, we gain first insights into the ability of the ciliate Paramecium to develop potentially adaptive phenotypic changes in response to early-life adversity. We show that, upon exposure to unconventional culture temperatures, germ line-to-soma differentiation gives rise to coordinated molecular changes that may help attune the number of functional gene copies to the new external conditions. The non-random somatic heterogeneity that developmental plasticity generates is largely epigenetically controlled, shaped by the parental experience, and may prompt a stress response. These findings establish Paramecium as a new model system to study the molecular basis and evolutionary significance of developmental plasticity. In echoing previous indications in mammals, they call for an incorporation of intergenerational effects in adaptation studies.

Programmed genome rearrangements in ciliates

Ciliates are a highly divergent group of unicellular eukaryotes with separate somatic and germline genomes found in distinct dimorphic nuclei. This characteristic feature is tightly linked to extremely laborious developmentally regulated genome rearrangements in the development of a new somatic genome/nuclei following sex. The transformation from germline to soma genome involves massive DNA elimination mediated by non-coding RNAs, chromosome fragmentation, as well as DNA amplification. In this review, we discuss the similarities and differences in the genome reorganization processes of the model ciliates Paramecium and Tetrahymena (class Oligohymenophorea), and the distantly related Euplotes, Stylonychia, and Oxytricha (class Spirotrichea).

The Role of Small Noncoding RNA in DNA Double-Strand Break Repair

DNA damage is a common phenomenon promoted through a variety of exogenous and endogenous factors. The DNA damage response (DDR) pathway involves a wide range of proteins, and as was indicated, small noncoding RNAs (sncRNAs). These are double-strand break-induced RNAs (diRNAs) and DNA damage response small RNA (DDRNA). Moreover, RNA binding proteins (RBPs) and RNA modifications have also been identified to modulate diRNA and DDRNA function in the DDR process. Several theories have been formulated regarding the synthesis and function of these sncRNAs during DNA repair; nevertheless, these pathways’ molecular details remain unclear. Here, we review the current knowledge regarding the mechanisms of diRNA and DDRNA biosynthesis and discuss the role of sncRNAs in maintaining genome stability.

The Paramecium histone chaperone Spt16-1 is required for Pgm endonuclease function in programmed genome rearrangements

In Paramecium tetraurelia, a large proportion of the germline genome is reproducibly removed from the somatic genome after sexual events via a process involving small (s)RNA-directed heterochromatin formation and DNA excision and repair. How germline limited DNA sequences are specifically recognized in the context of chromatin remains elusive. Here, we use a reverse genetics approach to identify factors involved in programmed genome rearrangements. We have identified a P. tetraurelia homolog of the highly conserved histone chaperone Spt16 subunit of the FACT complex, Spt16-1, and show its expression is developmentally regulated. A functional GFP-Spt16-1 fusion protein localized exclusively in the nuclei where genome rearrangements take place. Gene silencing of Spt16-1 showed it is required for the elimination of all germline-limited sequences, for the survival of sexual progeny, and for the accumulation of internal eliminated sequence (ies)RNAs, an sRNA population produced when elimination occurs. Normal accumulation of 25 nt scanRNAs and deposition of silent histone marks H3K9me3 and H3K27me3 indicated that Spt16-1 does not regulate the scanRNA-directed heterochromatin pathway involved in the early steps of DNA elimination. We further show that Spt16-1 is required for the correct nuclear localization of the PiggyMac (Pgm) endonuclease, which generates the DNA double-strand breaks required for DNA elimination. Thus, Spt16-1 is essential for Pgm function during programmed genome rearrangements. We propose a model in which Spt16-1 mediates interactions between the excision machinery and chromatin, facilitating endonuclease access to DNA cleavage sites during genome rearrangements.

The genome is generally similar in all the cells of an organism. However, in the ciliate Paramecium tetraurelia, massive and reproducible programmed DNA elimination leads to a highly streamlined somatic genome. In eukaryotes, DNA is packaged into nucleosomes, which ensure genome integrity but act as a barrier to enzymes acting on DNA. How the endonuclease PiggyMac gains access to the genome to initiate DNA elimination remains elusive. Here, we identified four P. tetraurelia genes encoding homologs of the conserved histone chaperone Spt16, which can modulate access to DNA by promoting nucleosome assembly and disassembly. We demonstrated that the most divergent gene, SPT16-1, has a highly specialized expression pattern, similar to that of PiggyMac, and a specific role in programmed DNA elimination. We show that the Spt16-1 protein, like PiggyMac, is exclusively localized in the differentiating somatic nucleus, and is also required for the dramatic elimination of germline-limited sequences. We further show that Spt16-1 directs the correct nuclear localization of the PiggyMac endonuclease. Thus, Spt16-1 is essential for PiggyMac function during programmed DNA elimination. We propose that Spt16-1 mediates the interaction between PiggyMac and chromatin or DNA, facilitating endonuclease access to DNA cleavage sites.

Roles of Noncoding RNAs in Ciliate Genome Architecture

Ciliates are an interesting model system for investigating diverse functions of noncoding RNAs, especially in genome defence pathways. During sexual development, the ciliate somatic genome undergoes massive rearrangement and reduction through removal of transposable elements and other repetitive DNA. This is guided by a multitude of noncoding RNAs of different sizes and functions, the extent of which is only recently becoming clear. The genome rearrangement pathways evolved as a defence against parasitic DNA, but interestingly also use the transposable elements and transposases to execute their own removal. Thus, ciliates are also a good model for the coevolution of host and transposable element, and the mutual dependence between the two. In this review, we summarise the genome rearrangement pathways in three diverse species of ciliate, with focus on recent discoveries and the roles of noncoding RNAs.

•

Ciliate genomes undergo massive rearrangement and reduction during development.

•

Transposon elimination is guided by small RNAs and carried out by transposases.

•

New pathways for noncoding RNA production have recently been discovered in ciliates.

•

Diverse ciliate species have different mechanisms for RNA-guided genome remodeling.

Conserved chromosomal functions of RNA interference

RNA interference (RNAi), a cellular process through which small RNAs target and regulate complementary RNA transcripts, has well-characterized roles in post-transcriptional gene regulation and transposon repression. Recent studies have revealed additional conserved roles for RNAi proteins, such as Argonaute and Dicer, in chromosome function. By guiding chromatin modification, RNAi components promote chromosome segregation during both mitosis and meiosis and regulate chromosomal and genomic dosage response. Small RNAs and the RNAi machinery also participate in the resolution of DNA damage. Interestingly, many of these lesser-studied functions seem to be more strongly conserved across eukaryotes than are well-characterized functions such as the processing of microRNAs. These findings have implications for the evolution of RNAi since the last eukaryotic common ancestor, and they provide a more complete view of the functions of RNAi.

Exogenous RNAi mechanisms contribute to transcriptome adaptation by phased siRNA clusters in Paramecium

Extensive research has characterized distinct exogenous RNAi pathways interfering in gene expression during vegetative growth of the unicellular model ciliate Paramecium. However, role of RNAi in endogenous transcriptome regulation, and environmental adaptation is unknown. Here, we describe the first genome-wide profiling of endogenous sRNAs in context of different transcriptomic states (serotypes). We developed a pipeline to identify, and characterize 2602 siRNA producing clusters (SRCs). Our data show no evidence that SRCs produce miRNAs, and in contrast to other species, no preference for strand specificity of siRNAs. Interestingly, most SRCs overlap coding genes and a separate group show siRNA phasing along the entire open reading frame, suggesting that the mRNA transcript serves as a source for siRNAs. Integrative analysis of siRNA abundance and gene expression levels revealed surprisingly that mRNA and siRNA show negative as well as positive associations. Two RNA-dependent RNA Polymerase mutants, RDR1 and RDR2, show a drastic loss of siRNAs especially in phased SRCs accompanied with increased mRNA levels. Importantly, most SRCs depend on both RDRs, reminiscent to primary siRNAs in the RNAi against exogenous RNA, indicating mechanistic overlaps between exogenous and endogenous RNAi contributing to flexible transcriptome adaptation.

Environmentally induced plasticity of programmed DNA elimination boosts somatic variability in Paramecium tetraurelia

Can ecological changes impact somatic genome development? Efforts to resolve this question could reveal a direct link between environmental changes and somatic variability, potentially illuminating our understanding of how variation can surface from a single genotype under stress. Here, we tackle this question by leveraging the biological properties of ciliates. When Paramecium tetraurelia reproduces sexually, its polyploid somatic genome regenerates from the germline genome through a developmental process that involves the removal of thousands of ORF-interrupting sequences known as internal eliminated sequences (IESs). We show that exposure to nonstandard culture temperatures impacts the efficiency of this process of programmed DNA elimination, prompting the emergence of hundreds of incompletely excised IESs in the newly developed somatic genome. These alternative DNA isoforms display a patterned genomic topography, impact gene expression, and might be inherited transgenerationally. On this basis, we conclude that environmentally induced developmental thermoplasticity contributes to genotypic diversification in Paramecium.

Diversification of small RNA amplification mechanisms for targeting transposon-related sequences in ciliates

The silencing of repetitive transposable elements (TEs) is ensured by signal amplification of the initial small RNA trigger, which occurs at distinct steps of TE silencing in different eukaryotes. How such a variety of secondary small RNA biogenesis mechanisms has evolved has not been thoroughly elucidated. Ciliated protozoa perform small RNA-directed programmed DNA elimination of thousands of TE-related internal eliminated sequences (IESs) in the newly developed somatic nucleus. In the ciliate Paramecium, secondary small RNAs are produced after the excision of IESs. In this study, we show that in another ciliate, Tetrahymena, secondary small RNAs accumulate at least a few hours before their derived IESs are excised. We also demonstrate that DNA excision is dispensable for their biogenesis in this ciliate. Therefore, unlike in Paramecium, small RNA amplification occurs before IES excision in Tetrahymena This study reveals the remarkable diversity of secondary small RNA biogenesis mechanisms, even among ciliates with similar DNA elimination processes, and thus raises the possibility that the evolution of TE-targeting small RNA amplification can be traced by investigating the DNA elimination mechanisms of ciliates.

A mating-type mutagenesis screen identifies a zinc-finger protein required for specific DNA excision events in Paramecium

In the ciliate Paramecium tetraurelia, functional genes are reconstituted during development of the somatic macronucleus through the precise excision of ∼45 000 single-copy Internal Eliminated Sequences (IESs), thought to be the degenerate remnants of ancient transposon insertions. Like introns, IESs are marked only by a weak consensus at their ends. How such a diverse set of sequences is faithfully recognized and precisely excised remains unclear: specialized small RNAs have been implicated, but in their absence up to ∼60% of IESs are still correctly excised. To get further insight, we designed a mutagenesis screen based on the hypersensitivity of a specific excision event in the mtA gene, which determines mating types. Unlike most IES-containing genes, the active form of mtA is the unexcised one, allowing the recovery of hypomorphic alleles of essential IES recognition/excision factors. Such is the case of one mutation recovered in the Piwi gene PTIWI09, a key player in small RNA-mediated IES recognition. Another mutation identified a novel protein with a C2H2 zinc finger, mtGa, which is required for excision of a small subset of IESs characterized by enrichment in a 5-bp motif. The unexpected implication of a sequence-specific factor establishes a new paradigm for IES recognition and/or excision.

Dicer-like Enzymes with Sequence Cleavage Preferences

Dicer proteins are known to produce small RNAs (sRNAs) from long double-stranded RNA (dsRNA) templates. These sRNAs are bound by Argonaute proteins, which select the guide strand, often with a 5' end sequence bias. However, Dicer proteins have never been shown to have sequence cleavage preferences. In Paramecium development, two classes of sRNAs that are required for DNA elimination are produced by three Dicer-like enzymes: Dcl2, Dcl3, and Dcl5. Through in vitro cleavage assays, we demonstrate that Dcl2 has a strict size preference for 25 nt and a sequence preference for 5' U and 5' AGA, while Dcl3 has a sequence preference for 5' UNG. Dcl5, however, has cleavage preferences for 5' UAG and 3' CUAC/UN, which leads to the production of RNAs precisely matching short excised DNA elements with corresponding end base preferences. Thus, we characterize three Dicer-like enzymes that are involved in Paramecium development and propose a biological role for their sequence-biased cleavage products.

The Challenges of Genome-Wide Studies in a Unicellular Eukaryote With Two Nuclear Genomes

We present here methods to study a eukaryotic microorganism with two nuclear genomes, both originating from the same zygotic genome. Paramecium, like other ciliates, is characterized by nuclear dimorphism, which is the presence of two types of nuclei with distinct organization and functions in the same cytoplasm. The two diploid germline micronuclei (MIC) undergo meiosis and fertilization to transmit the genetic information across sexual generations. The highly polyploid somatic macronucleus (MAC) contains a reduced version of the genome optimized for gene expression. Reproducible programmed DNA elimination of about 30% of the complexity of the 100Mb MIC genome occurs during development of the MAC along with endoreplication to 800 copies. Large regions that contain transposable elements and other repeats are eliminated, and short single copy remnants of transposable elements, which often interrupt coding sequences, are precisely excised to restore functional open reading frames. Genome-wide studies of this process require access to MIC DNA which has long been impossible. The breakthrough with respect to this technical obstacle came with development of a MIC purification protocol involving a critical step of flow cytometry to sort nuclei representing only 0.5% of total genomic DNA. Here, we provide a step-by-step protocol and important tips for purifying nuclei, and present the methods developed for downstream analysis of NGS data.

Whats, hows and whys of programmed DNA elimination in Tetrahymena

Programmed genome rearrangements in ciliates provide fascinating examples of flexible epigenetic genome regulations and important insights into the interaction between transposable elements (TEs) and host genomes. DNA elimination in Tetrahymena thermophila removes approximately 12 000 internal eliminated sequences (IESs), which correspond to one-third of the genome, when the somatic macronucleus (MAC) differentiates from the germline micronucleus (MIC). More than half of the IESs, many of which show high similarity to TEs, are targeted for elimination in cis by the small RNA-mediated genome comparison of the MIC to the MAC. Other IESs are targeted for elimination in trans by the same small RNAs through repetitive sequences. Furthermore, the small RNA–heterochromatin feedback loop ensures robust DNA elimination. Here, we review an updated picture of the DNA elimination mechanism, discuss the physiological and evolutionary roles of DNA elimination, and outline the key questions that remain unanswered.

Two Sets of Piwi Proteins Are Involved in Distinct sRNA Pathways Leading to Elimination of Germline-Specific DNA

Piwi proteins and piRNAs protect eukaryotic germlines against the spread of transposons. During development in the ciliate Paramecium, two Piwi-dependent sRNA classes are involved in the elimination of transposons and transposon-derived DNA: scan RNAs (scnRNAs), associated with Ptiwi01 and Ptiwi09, and iesRNAs, whose binding partners we now identify as Ptiwi10 and Ptiwi11. scnRNAs derive from the maternal genome and initiate DNA elimination during development, whereas iesRNAs continue DNA targeting until the removal process is complete. Here, we show that scnRNAs and iesRNAs are processed by distinct Dicer-like proteins and bind Piwi proteins in a mutually exclusive manner, suggesting separate biogenesis pathways. We also demonstrate that the PTIWI10 gene is transcribed from the developing nucleus and that its transcription depends on prior DNA excision, suggesting a mechanism of gene expression control triggered by the removal of short DNA segments interrupting the gene.

•

Identification of two Piwi proteins (Ptiwi10/11) associated with iesRNAs

•

Piwi proteins bind Dicer-produced sRNAs and remove passenger strands

•

Ptiwi10 is expressed from the new somatic macronucleus

•

DNA elimination activates gene transcription

From "Cellular" RNA to "Smart" RNA: Multiple Roles of RNA in Genome Stability and Beyond

Coding for proteins has been considered the main function of RNA since the "central dogma" of biology was proposed. The discovery of noncoding transcripts shed light on additional roles of RNA, ranging from the support of polypeptide synthesis, to the assembly of subnuclear structures, to gene expression modulation. Cellular RNA has therefore been recognized as a central player in often unanticipated biological processes, including genomic stability. This ever-expanding list of functions inspired us to think of RNA as a "smart" phone, which has replaced the older obsolete "cellular" phone. In this review, we summarize the last two decades of advances in research on the interface between RNA biology and genome stability. We start with an account of the emergence of noncoding RNA, and then we discuss the involvement of RNA in DNA damage signaling and repair, telomere maintenance, and genomic rearrangements. We continue with the depiction of single-molecule RNA detection techniques, and we conclude by illustrating the possibilities of RNA modulation in hopes of creating or improving new therapies. The widespread biological functions of RNA have made this molecule a reoccurring theme in basic and translational research, warranting it the transcendence from classically studied "cellular" RNA to "smart" RNA.

RNA-mediated transgenerational inheritance in ciliates and plants

In the age of next-generation sequencing (NGS) and with the availability of whole sequenced genomes and epigenomes, some attention has shifted from purely sequence-based studies to those of heritable epigenetic modifications. Transgenerational inheritance can be defined as heritable changes to the state of DNA that may be passed on to subsequent generations without alterations to the underlying DNA sequence. Although this phenomenon has been extensively studied in many systems, studies of transgenerational inheritance in mammals and other higher-level eukaryotes may be complicated by the fact that many epigenetic marks are reprogrammed during sexual reproduction. This, by definition, may obscure our interpretation of what is in fact truly transgenerational. Therefore, in this mini review, we discuss what is currently known in the field about transgenerational epigenetic inheritance in ciliates and plants, with a particular emphasis on RNA-mediated processes and changes in chromatin states.

Environmental Temperature Controls Accumulation of Transacting siRNAs Involved in Heterochromatin Formation

Genes or alleles can interact by small RNAs in a homology dependent manner meaning that short interfering (siRNAs) can act in trans at the chromatin level producing stable and heritable silencing phenotypes. Because of the puzzling data on endogenous paramutations, their impact contributing to adaptive evolution in a Lamarckian manner remains unknown. An increasing number of studies characterizes the underlying siRNA accumulation pathways using transgene experiments. Also in the ciliate Paramecium tetraurelia, we induce trans silencing on the chromatin level by injection of truncated transgenes. Here, we characterize the efficiency of this mechanism at different temperatures showing that silencing of the endogenous genes is temperature dependent. Analyzing different transgene constructs at different copy numbers, we dissected whether silencing efficiency is due to varying precursor RNAs or siRNA accumulation. Our data shows that silencing efficiency correlates with more efficient accumulation of primary siRNAs at higher temperatures rather than higher expression of precursor RNAs. Due to higher primary levels, secondary siRNAs also show temperature dependency and interestingly increase their relative proportion to primary siRNAs. Our data shows that efficient trans silencing on the chromatin level in P. tetraurelia depends on environmental parameters, thus being an important epigenetic factor limiting regulatory effects of siRNAs.

Negative Regulators of an RNAi-Heterochromatin Positive Feedback Loop Safeguard Somatic Genome Integrity in Tetrahymena

RNAi-mediated positive feedback loops are pivotal for the maintenance of heterochromatin, but how they are downregulated at heterochromatin-euchromatin borders is not well understood. In the ciliated protozoan Tetrahymena, heterochromatin is formed exclusively on the sequences that are removed from the somatic genome by programmed DNA elimination, and an RNAi-mediated feedback loop is important for assembling heterochromatin on the eliminated sequences. In this study, we show that the heterochromatin protein 1 (HP1)-like protein Coi6p, its interaction partners Coi7p and Lia5p, and the histone demethylase Jmj1p are crucial for confining the production of small RNAs and the formation of heterochromatin to the eliminated sequences. The loss of Coi6p, Coi7p, or Jmj1p causes ectopic DNA elimination. The results provide direct evidence for the existence of a dedicated mechanism that counteracts a positive feedback loop between RNAi and heterochromatin at heterochromatin-euchromatin borders to maintain the integrity of the somatic genome.

•

The HP1-like protein Coi6p confines small RNA and heterochromatin formation

•

Two Coi6p-binding proteins and the histone demethylase Jmj1p likely act with Coi6p

•

Coi6p and Jmj1p are important for preventing ectopic DNA elimination

•

Suppression of RNAi-heterochromatin feedback loop maintains somatic genome integrity

A meiosis-specific Spt5 homolog involved in non-coding transcription

Spt5 is a conserved and essential transcriptional regulator that binds directly to RNA polymerase and is involved in transcription elongation, polymerase pausing and various co-transcriptional processes. To investigate the role of Spt5 in non-coding transcription, we used the unicellular model Paramecium tetraurelia. In this ciliate, development is controlled by epigenetic mechanisms that use different classes of non-coding RNAs to target DNA elimination. We identified two SPT5 genes. One (STP5v) is involved in vegetative growth, while the other (SPT5m) is essential for sexual reproduction. We focused our study on SPT5m, expressed at meiosis and associated with germline nuclei during sexual processes. Upon Spt5m depletion, we observed absence of scnRNAs, piRNA-like 25 nt small RNAs produced at meiosis. The scnRNAs are a temporal copy of the germline genome and play a key role in programming DNA elimination. Moreover, Spt5m depletion abolishes elimination of all germline-limited sequences, including sequences whose excision was previously shown to be scnRNA-independent. This suggests that in addition to scnRNA production, Spt5 is involved in setting some as yet uncharacterized epigenetic information at meiosis. Our study establishes that Spt5m is crucial for developmental genome rearrangements and necessary for scnRNA production.

Thousands of RNA-cached copies of whole chromosomes are present in the ciliate Oxytricha during development

The ciliate Oxytricha trifallax maintains two genomes: a germline genome that is active only during sexual conjugation and a transcriptionally active, somatic genome that derives from the germline via extensive sequence reduction and rearrangement. Previously, we found that long noncoding (lnc) RNA "templates"-telomere-containing, RNA-cached copies of mature chromosomes-provide the information to program the rearrangement process. Here we used a modified RNA-seq approach to conduct the first genome-wide search for endogenous, telomere-to-telomere RNA transcripts. We find that during development, Oxytricha produces long noncoding RNA copies for over 10,000 of its 16,000 somatic chromosomes, consistent with a model in which Oxytricha transmits an RNA-cached copy of its somatic genome to the sexual progeny. Both the primary sequence and expression profile of a somatic chromosome influence the temporal distribution and abundance of individual template RNAs. This suggests that Oxytricha may undergo multiple rounds of DNA rearrangement during development. These observations implicate a complex set of thousands of long RNA molecules in the wiring and maintenance of a highly elaborate somatic genome architecture.

Genetics and Epigenetics of Mating Type Determination in Paramecium and Tetrahymena

While sex is an ancient and highly conserved eukaryotic invention, self-incompatibility systems such as mating types or sexes appear to be derived limitations that show considerable evolutionary plasticity. Within a single class of ciliates, Paramecium and Tetrahymena species have long been known to present a wide variety of mating type numbers and modes of inheritance, but only recently have the genes involved been identified. Although similar transmembrane proteins mediate self/nonself recognition in both ciliates, the mechanisms of mating type determination differ widely, ranging from Mendelian systems to developmental nuclear differentiation, either stochastic or maternally inherited. The non-Mendelian systems rely on programmed editing of the germline genome that occurs during differentiation of the somatic nucleus, and they have co-opted different DNA recombination mechanisms-some previously unknown. Here we review the recent molecular advances and some remaining unsolved questions and discuss the possible implications of these diverse mechanisms for inbreeding/outbreeding balance regulation.