RNA-mediated programming of developmental genome rearrangements in Paramecium tetraurelia

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.

Programmed chromosome fragmentation in ciliated protozoa: multiple means to chromosome ends

SUMMARYCiliated protozoa undergo large-scale developmental rearrangement of their somatic genomes when forming a new transcriptionally active macronucleus during conjugation. This process includes the fragmentation of chromosomes derived from the germline, coupled with the efficient healing of the broken ends by de novo telomere addition. Here, we review what is known of developmental chromosome fragmentation in ciliates that have been well-studied at the molecular level (Tetrahymena, Paramecium, Euplotes, Stylonychia, and Oxytricha). These organisms differ substantially in the fidelity and precision of their fragmentation systems, as well as in the presence or absence of well-defined sequence elements that direct excision, suggesting that chromosome fragmentation systems have evolved multiple times and/or have been significantly altered during ciliate evolution. We propose a two-stage model for the evolution of the current ciliate systems, with both stages involving repetitive or transposable elements in the genome. The ancestral form of chromosome fragmentation is proposed to have been derived from the ciliate small RNA/chromatin modification process that removes transposons and other repetitive elements from the macronuclear genome during development. The evolution of this ancestral system is suggested to have potentiated its replacement in some ciliate lineages by subsequent fragmentation systems derived from mobile genetic elements.

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.

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.

Membrane Traffic and Ca2+ -Signals in Ciliates

A Paramecium cell has as many types of membrane interactions as mammalian cells, as established with monoclonal antibodies by R. Allen and A. Fok. Since then, we have identified key-players, such as SNARE-proteins, Ca²⁺ -regulating proteins, including Ca²⁺ -channels, Ca²⁺ -pumps, Ca²⁺ -binding proteins of different affinity etc. at the molecular level, probed their function and localized them at the light and electron microscopy level. SNARE-proteins, in conjunction with a synaptotagmin-like Ca²⁺ -sensor protein, mediate membrane fusion. This interaction is additionally regulated by monomeric GTPases whose spectrum in Tetrahymena and Paramecium has been established by A. Turkewitz. As known from mammalian cells, GTPases are activated on membranes in conjunction with lumenal acidification by an H⁺ -ATPase. For these complex molecules we found in Paramecium an unsurpassed number of 17 a-subunit paralogs which connect the polymeric head and basis part, V1 and V0. (This multitude may reflect different local functional requirements.) Together with plasmalemmal Ca²⁺ -influx-channels, locally enriched intracellular InsP₃ -type (InsP₃ R, mainly in osmoregulatory system) and ryanodine receptor-like Ca²⁺ -release channels (ryanodine receptor-like proteins, RyR-LP), this complexity mediates Ca²⁺ signals for most flexible local membrane-to-membrane interactions. As we found, the latter channel types miss a substantial portion of the N-terminal part. Caffeine and 4-chloro-meta-cresol (the agent used to probe mutations of RyRs in man during surgery in malignant insomnia patients) initiate trichocyst exocytosis by activating Ca²⁺ -release channels type CRC-IV in the peripheral part of alveolar sacs. This is superimposed by Ca²⁺ -influx, i.e. a mechanism called "store-operated Ca²⁺ -entry" (SOCE). For the majority of key players, we have mapped paralogs throughout the Paramecium cell, with features in common or at variance in the different organelles participating in vesicle trafficking. Local values of free Ca²⁺ -concentration, [Ca²⁺ ]_i , and their change, e.g. upon exocytosis stimulation, have been registered by flurochromes and chelator effects. In parallel we have registered release of Ca²⁺ from alveolar sacs by quenched-flow analysis combined with cryofixation and x-ray microanalysis.

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.

Evolutionary Plasticity of Mating-Type Determination Mechanisms in Paramecium aurelia Sibling Species

The Paramecium aurelia complex, a group of morphologically similar but sexually incompatible sibling species, is a unique example of the evolutionary plasticity of mating-type systems. Each species has two mating types, O (Odd) and E (Even). Although O and E types are homologous in all species, three different modes of determination and inheritance have been described: genetic determination by Mendelian alleles, stochastic developmental determination, and maternally inherited developmental determination. Previous work in three species of the latter kind has revealed the key roles of the E-specific transmembrane protein mtA and its highly specific transcription factor mtB: type O clones are produced by maternally inherited genome rearrangements that inactivate either mtA or mtB during development. Here we show, through transcriptome analyses in five additional species representing the three determination systems, that mtA expression specifies type E in all cases. We further show that the Mendelian system depends on functional and nonfunctional mtA alleles, and identify novel developmental rearrangements in mtA and mtB which now explain all cases of maternally inherited mating-type determination. Epistasis between these genes likely evolved from less specific interactions between paralogs in the P. aurelia common ancestor, after a whole-genome duplication, but the mtB gene was subsequently lost in three P. aurelia species which appear to have returned to an ancestral regulation mechanism. These results suggest a model accounting for evolutionary transitions between determination systems, and highlight the diversity of molecular solutions explored among sibling species to maintain an essential mating-type polymorphism in cell populations.

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.

Functional diversification of Paramecium Ku80 paralogs safeguards genome integrity during precise programmed DNA elimination

Gene duplication and diversification drive the emergence of novel functions during evolution. Because of whole genome duplications, ciliates from the Paramecium aurelia group constitute a remarkable system to study the evolutionary fate of duplicated genes. Paramecium species harbor two types of nuclei: a germline micronucleus (MIC) and a somatic macronucleus (MAC) that forms from the MIC at each sexual cycle. During MAC development, ~45,000 germline Internal Eliminated Sequences (IES) are excised precisely from the genome through a 'cut-and-close' mechanism. Here, we have studied the P. tetraurelia paralogs of KU80, which encode a key DNA double-strand break repair factor involved in non-homologous end joining. The three KU80 genes have different transcription patterns, KU80a and KU80b being constitutively expressed, while KU80c is specifically induced during MAC development. Immunofluorescence microscopy and high-throughput DNA sequencing revealed that Ku80c stably anchors the PiggyMac (Pgm) endonuclease in the developing MAC and is essential for IES excision genome-wide, providing a molecular explanation for the previously reported Ku-dependent licensing of DNA cleavage at IES ends. Expressing Ku80a under KU80c transcription signals failed to complement a depletion of endogenous Ku80c, indicating that the two paralogous proteins have distinct properties. Domain-swap experiments identified the α/β domain of Ku80c as the major determinant for its specialized function, while its C-terminal part is required for excision of only a small subset of IESs located in IES-dense regions. We conclude that Ku80c has acquired the ability to license Pgm-dependent DNA cleavage, securing precise DNA elimination during programmed rearrangements. The present study thus provides novel evidence for functional diversification of genes issued from a whole-genome duplication.

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.

The Transmembrane Protein Semi1 Positions Gamete Nuclei for Reciprocal Fertilization in Tetrahymena

During sexual reproduction in the ciliate, Tetrahymena thermophila, cells of complementary mating type pair ("conjugate") undergo simultaneous meiosis and fertilize each other. In both mating partners only one of the four meiotic products is "selected" to escape autophagy, and this nucleus divides mitotically to produce two pronuclei. The migrating pronucleus of one cell translocates to the mating partner and fuses with its stationary pronucleus and vice versa. Selection of the designated gametic nucleus was thought to depend on its position within the cell because it always attaches to the junction with the partner cell. Here we show that a transmembrane protein, Semi1, is crucial for attachment. Loss of Semi1 causes failure to attach and consequent infertility. However, a nucleus is selected and gives rise to pronuclei regardless of Semi1 expression, indicating that attachment of a nucleus to the junction is not a precondition for selection but follows the selection process.

Programmed Chromosome Deletion in the Ciliate Oxytricha trifallax

The ciliate Oxytricha trifallax contains two nuclei: a germline micronucleus and a somatic macronucleus. These two nuclei diverge significantly in genomic structure. The micronucleus contains approximately 100 chromosomes of megabase scale, while the macronucleus contains 16,000 gene-sized, high ploidy "nanochromosomes." During its sexual cycle, a copy of the zygotic germline micronucleus develops into a somatic macronucleus via DNA excision and rearrangement. The rearrangement process is guided by multiple RNA-based pathways that program the epigenetic inheritance of sequences in the parental macronucleus of the subsequent generation. Here, we show that the introduction of synthetic DNA molecules homologous to a complete native nanochromosome during the rearrangement process results in either loss or heavy copy number reduction of the targeted nanochromosome in the macronucleus of the subsequent generation. This phenomenon was tested on a variety of nanochromosomes with different micronuclear structures, with deletions resulting in all cases. Deletion of the targeted nanochromosome results in the loss of expression of the targeted genes, including gene knockout phenotypes that were phenocopied using alternative knockdown approaches. Further investigation of the chromosome deletion showed that, although the full length nanochromosome was lost, remnants of the targeted chromosome remain. We were also able to detect the presence of telomeres on these remnants. The chromosome deletions and remnants are epigenetically inherited when backcrossed to wild type strains, suggesting that an undiscovered mechanism programs DNA elimination and cytoplasmically transfers to both daughter cells during conjugation. Programmed deletion of targeted chromosomes provides a novel approach to investigate genome rearrangement and expands the available strategies for gene knockout in Oxytricha trifallax.

Loss of a Fragile Chromosome Region leads to the Screwy Phenotype in Paramecium tetraurelia

A conspicuous cell-shape phenotype known as “screwy” was reported to result from mutations at two or three uncharacterized loci in the ciliate Paramecium tetraurelia. Here, we describe a new screwy mutation, Spinning Top, which appeared spontaneously in the cross of an unrelated mutant with reference strain 51. The macronuclear (MAC) genome of the Spinning Top mutant is shown to lack a ~28.5-kb segment containing 18 genes at the end of one chromosome, which appears to result from a collinear deletion in the micronuclear (MIC) genome. We tested several candidate genes from the deleted locus by dsRNA-induced silencing in wild-type cells, and identified a single gene responsible for the phenotype. This gene, named Spade, encodes a 566-aa glutamine-rich protein with a C₂HC zinc finger. Its silencing leads to a fast phenotype switch during vegetative growth, but cells recover a wild-type phenotype only 5–6 divisions after silencing is stopped. We analyzed 5 independently-obtained mutant alleles of the Sc1 locus, and concluded that all of them also lack the Spade gene and a number of neighboring genes in the MAC and MIC genomes. Mapping of the MAC deletion breakpoints revealed two different positions among the 5 alleles, both of which differ from the Spinning Top breakpoint. These results suggest that this MIC chromosome region is intrinsically unstable in strain 51.

Evolutionary origins and impacts of genome architecture in ciliates

Genome architecture is well diversified among eukaryotes in terms of size and content, with many being radically shaped by ancient and ongoing genome conflicts with transposable elements (e.g., the large transposon‐rich genomes common among plants). In ciliates, a group of microbial eukaryotes with distinct somatic and germ‐line genomes present in a single cell, the consequences of these genome conflicts are most apparent in their developmentally programmed genome rearrangements. This complicated developmental phenomenon has largely overshadowed and outpaced our understanding of how germ‐line and somatic genome architectures have influenced the evolutionary dynamism and potential in these taxa. In our review, we highlight three central concepts: how the evolution of atypical ciliate germ‐line genome architectures is linked to ancient genome conflicts; how the complex, epigenetically guided transformation of germline to soma during development can generate widespread genetic variation; and how these features, coupled with their unusual life cycle, have increased the rate of molecular evolution linked to genome architecture in these taxa.

Six domesticated PiggyBac transposases together carry out programmed DNA elimination in Paramecium

The domestication of transposable elements has repeatedly occurred during evolution and domesticated transposases have often been implicated in programmed genome rearrangements, as remarkably illustrated in ciliates. In Paramecium, PiggyMac (Pgm), a domesticated PiggyBac transposase, carries out developmentally programmed DNA elimination, including the precise excision of tens of thousands of gene-interrupting germline Internal Eliminated Sequences (IESs). Here, we report the discovery of five groups of distant Pgm-like proteins (PgmLs), all able to interact with Pgm and essential for its nuclear localization and IES excision genome-wide. Unlike Pgm, PgmLs lack a conserved catalytic site, suggesting that they rather have an architectural function within a multi-component excision complex embedding Pgm. PgmL depletion can increase erroneous targeting of residual Pgm-mediated DNA cleavage, indicating that PgmLs contribute to accurately position the complex on IES ends. DNA rearrangements in Paramecium constitute a rare example of a biological process jointly managed by six distinct domesticated transposases.

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.

The fitness cost of mis-splicing is the main determinant of alternative splicing patterns

BACKGROUND

Most eukaryotic genes are subject to alternative splicing (AS), which may contribute to the production of protein variants or to the regulation of gene expression via nonsense-mediated messenger RNA (mRNA) decay (NMD). However, a fraction of splice variants might correspond to spurious transcripts and the question of the relative proportion of splicing errors to functional splice variants remains highly debated.

RESULTS

We propose a test to quantify the fraction of AS events corresponding to errors. This test is based on the fact that the fitness cost of splicing errors increases with the number of introns in a gene and with expression level. We analyzed the transcriptome of the intron-rich eukaryote Paramecium tetraurelia. We show that in both normal and in NMD-deficient cells, AS rates strongly decrease with increasing expression level and with increasing number of introns. This relationship is observed for AS events that are detectable by NMD as well as for those that are not, which invalidates the hypothesis of a link with the regulation of gene expression. Our results show that in genes with a median expression level, 92-98% of observed splice variants correspond to errors. We observed the same patterns in human transcriptomes and we further show that AS rates correlate with the fitness cost of splicing errors.

CONCLUSIONS

These observations indicate that genes under weaker selective pressure accumulate more maladaptive substitutions and are more prone to splicing errors. Thus, to a large extent, patterns of gene expression variants simply reflect the balance between selection, mutation, and drift.

A germline-limited piggyBac transposase gene is required for precise excision in Tetrahymena genome rearrangement

Developmentally programmed genome rearrangement accompanies differentiation of the silent germline micronucleus into the transcriptionally active somatic macronucleus in the ciliated protozoan Tetrahymena thermophila. Internal eliminated sequences (IES) are excised, followed by rejoining of MAC-destined sequences, while fragmentation occurs at conserved chromosome breakage sequences, generating macronuclear chromosomes. Some macronuclear chromosomes, referred to as non-maintained chromosomes (NMC), are lost soon after differentiation. Large NMC contain genes implicated in development-specific roles. One such gene encodes the domesticated piggyBac transposase TPB6, required for heterochromatin-dependent precise excision of IES residing within exons of functionally important genes. These conserved exonic IES determine alternative transcription products in the developing macronucleus; some even contain free-standing genes. Examples of precise loss of some exonic IES in the micronucleus and retention of others in the macronucleus of related species suggest an evolutionary analogy to introns. Our results reveal that germline-limited sequences can encode genes with specific expression patterns and development-related functions, which may be a recurring theme in eukaryotic organisms experiencing programmed genome rearrangement during germline to soma differentiation.

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.

Multimerization properties of PiggyMac, a domesticated piggyBac transposase involved in programmed genome rearrangements

During sexual processes, the ciliate Paramecium eliminates 25–30% of germline DNA from its somatic genome. DNA elimination includes excision of ∼45 000 short, single-copy internal eliminated sequences (IESs) and depends upon PiggyMac (Pgm), a domesticated piggyBac transposase that is essential for DNA cleavage at IES ends. Pgm carries a core transposase region with a putative catalytic domain containing three conserved aspartic acids, and a downstream cysteine-rich (CR) domain. A C-terminal extension of unknown function is predicted to adopt a coiled-coil (CC) structure. To address the role of the three domains, we designed an in vivo complementation assay by expressing wild-type or mutant Pgm-GFP fusions in cells depleted for their endogenous Pgm. The DDD triad and the CR domain are essential for Pgm activity and mutations in either domain have a dominant-negative effect in wild-type cells. A mutant lacking the CC domain is partially active in the presence of limiting Pgm amounts, but inactive when Pgm is completely absent, suggesting that presence of the mutant protein increases the overall number of active complexes. We conclude that IES excision involves multiple Pgm subunits, of which at least a fraction must contain the CC domain.

Evolutionary Cell Biology of Proteins from Protists to Humans and Plants

During evolution, the cell as a fine-tuned machine had to undergo permanent adjustments to match changes in its environment, while "closed for repair work" was not possible. Evolution from protists (protozoa and unicellular algae) to multicellular organisms may have occurred in basically two lineages, Unikonta and Bikonta, culminating in mammals and angiosperms (flowering plants), respectively. Unicellular models for unikont evolution are myxamoebae (Dictyostelium) and increasingly also choanoflagellates, whereas for bikonts, ciliates are preferred models. Information accumulating from combined molecular database search and experimental verification allows new insights into evolutionary diversification and maintenance of genes/proteins from protozoa on, eventually with orthologs in bacteria. However, proteins have rarely been followed up systematically for maintenance or change of function or intracellular localization, acquirement of new domains, partial deletion (e.g. of subunits), and refunctionalization, etc. These aspects are discussed in this review, envisaging "evolutionary cell biology." Protozoan heritage is found for most important cellular structures and functions up to humans and flowering plants. Examples discussed include refunctionalization of voltage-dependent Ca²⁺ channels in cilia and replacement by other types during evolution. Altogether components serving Ca²⁺ signaling are very flexible throughout evolution, calmodulin being a most conservative example, in contrast to calcineurin whose catalytic subunit is lost in plants, whereas both subunits are maintained up to mammals for complex functions (immune defense and learning). Domain structure of R-type SNAREs differs in mono- and bikonta, as do Ca²⁺ -dependent protein kinases. Unprecedented selective expansion of the subunit a which connects multimeric base piece and head parts (V0, V1) of H⁺ -ATPase/pump may well reflect the intriguing vesicle trafficking system in ciliates, specifically in Paramecium. One of the most flexible proteins is centrin when its intracellular localization and function throughout evolution is traced. There are many more examples documenting evolutionary flexibility of translation products depending on requirements and potential for implantation within the actual cellular context at different levels of evolution. From estimates of gene and protein numbers per organism, it appears that much of the basic inventory of protozoan precursors could be transmitted to highest eukaryotic levels, with some losses and also with important additional "inventions."

Unusual features of non-dividing somatic macronuclei in the ciliate class Karyorelictea

Genome structure and nuclear organization have been intensely studied in model ciliates such as Tetrahymena and Paramecium, yet few studies have focused on nuclear features of other ciliate clades including the class Karyorelictea. In most ciliates, both the somatic macronuclei and germline micronuclei divide during cell division and macronuclear development only occurs after conjugation. However, the macronuclei of Karyorelictea are non-dividing (i.e. division minus (Div-)) and develop anew from micronuclei during each asexual division. As macronuclei age within Karyorelictea, they undergo changes in morphology and DNA content until they are eventually degraded and replaced by newly developed macronuclei. No less than two macronuclei and one micronucleus are present in karyorelictid species, which suggests that a mature macronucleus 1) might be needed to sustain the cell while a new macronucleus is developing and 2) likely plays a role in guiding the development of the new macronucleus. Here we use a phylogenetic framework to compile information on the morphology and development of nuclei in Karyorelictea, largely relying on the work of Dr. Igor Raikov (1932-1998). We synthesize data to speculate on the functional implications of key features of Karyorelictea including the presence of at least two macronuclei in each cell and the inability for macronuclei to divide.

piRNAs and their diverse roles: a transposable element-driven tactic for gene regulation?

P-element-induced wimpy testis (PIWI)-interacting RNAs (piRNAs) are small, noncoding RNAs known for silencing transposable elements (TEs) in the germline of animals. Most genomes host TEs, which are notorious for mobilizing themselves and endangering survival of the host if not controlled. By silencing TEs in the germline, piRNAs prevent harmful mutations from being passed on to the next generation. How piRNAs are generated and how they silence TEs were the focus of researchers ever since their discovery. Now a spate of recent papers are beginning to tell us that piRNAs can play roles beyond TE silencing and are involved in diverse cellular processes from mRNA regulation to development or genome rearrangement. In this review, we discuss some of these recently reported roles. Data on these new roles are often rudimentary, and the involvement of piRNAs in these processes is yet to be definitely established. What is interesting is that the reports are on animals widely separated on the phylogenetic tree of life and that piRNAs were also found outside the gonadal tissues. Some of these piRNAs map to TE sequences, prompting us to hypothesize that genomes may have co-opted the TE-derived piRNA system for their own regulation.-Sarkar, A., Volff, J.-N., Vaury, C. piRNAs and their diverse roles: a transposable element-driven tactic for gene regulation?

Classical non-homologous end-joining pathway utilizes nascent RNA for error-free double-strand break repair of transcribed genes

DNA double-strand breaks (DSBs) leading to loss of nucleotides in the transcribed region can be lethal. Classical non-homologous end-joining (C-NHEJ) is the dominant pathway for DSB repair (DSBR) in adult mammalian cells. Here we report that during such DSBR, mammalian C-NHEJ proteins form a multiprotein complex with RNA polymerase II and preferentially associate with the transcribed genes after DSB induction. Depletion of C-NHEJ factors significantly abrogates DSBR in transcribed but not in non-transcribed genes. We hypothesized that nascent RNA can serve as a template for restoring the missing sequences, thus allowing error-free DSBR. We indeed found pre-mRNA in the C-NHEJ complex. Finally, when a DSB-containing plasmid with several nucleotides deleted within the E. coli lacZ gene was allowed time to repair in lacZ-expressing mammalian cells, a functional lacZ plasmid could be recovered from control but not C-NHEJ factor-depleted cells, providing important mechanistic insights into C-NHEJ-mediated error-free DSBR of the transcribed genome.

Most adult mammalian cells prefer to repair double-strand DNA breaks though the classical nonhomologous end-joining pathway. Here the authors present evidence that a nascent RNA transcript can serve as a template to facilitate error-free repair.

Chromosome fusions triggered by noncoding RNA

Chromosomal fusions are common in normal and cancer cells and can produce aberrant gene products that promote transformation. The mechanisms driving these fusions are poorly understood, but recurrent fusions are widespread. This suggests an underlying mechanism, and some authors have proposed a possible role for RNA in this process. The unicellular eukaryote Oxytricha trifallax displays an exorbitant capacity for natural genome editing, when it rewrites its germline genome to form a somatic epigenome. This developmental process provides a powerful model system to directly test the influence of small noncoding RNAs on chromosome fusion events during somatic differentiation. Here we show that small RNAs are capable of inducing chromosome fusions in 4 distinct cases (out of 4 tested), including one fusion of 3 chromosomes. We further show that these RNA-mediated chromosome fusions are heritable over multiple sexual generations and that transmission of the acquired fusion is associated with endogenous production of novel piRNA molecules that target the fused junction. We also demonstrate the capacity of a long noncoding RNA (lncRNA) to induce chromosome fusion of 2 distal germline loci. These results underscore the ability of short-lived, aberrant RNAs to act as drivers of chromosome fusion events that can be stably transmitted to future generations.

Two sets of RNAi components are required for heterochromatin formation in trans triggered by truncated transgenes

Across kingdoms, RNA interference (RNAi) has been shown to control gene expression at the transcriptional- or the post-transcriptional level. Here, we describe a mechanism which involves both aspects: truncated transgenes, which fail to produce intact mRNA, induce siRNA accumulation and silencing of homologous loci in trans in the ciliate Paramecium. We show that silencing is achieved by co-transcriptional silencing, associated with repressive histone marks at the endogenous gene. This is accompanied by secondary siRNA accumulation, strictly limited to the open reading frame of the remote locus. Our data shows that in this mechanism, heterochromatic marks depend on a variety of RNAi components. These include RDR3 and PTIWI14 as well as a second set of components, which are also involved in post-transcriptional silencing: RDR2, PTIWI13, DCR1 and CID2. Our data indicates differential processing of nascent un-spliced and long, spliced transcripts thus suggesting a hitherto-unrecognized functional interaction between post-transcriptional and co-transcriptional RNAi. Both sets of RNAi components are required for efficient trans-acting RNAi at the chromatin level and our data indicates similar mechanisms contributing to genome wide regulation of gene expression by epigenetic mechanisms.

Programmed Rearrangement in Ciliates: Paramecium

Programmed genome rearrangements in the ciliate Paramecium provide a nice illustration of the impact of transposons on genome evolution and plasticity. During the sexual cycle, development of the somatic macronucleus involves elimination of ∼30% of the germline genome, including repeated DNA (e.g., transposons) and ∼45,000 single-copy internal eliminated sequences (IES). IES excision is a precise cut-and-close process, in which double-stranded DNA cleavage at IES ends depends on PiggyMac, a domesticated piggyBac transposase. Genome-wide analysis has revealed that at least a fraction of IESs originate from Tc/mariner transposons unrelated to piggyBac. Moreover, genomic sequences with no transposon origin, such as gene promoters, can be excised reproducibly as IESs, indicating that genome rearrangements contribute to the control of gene expression. How the system has evolved to allow elimination of DNA sequences with no recognizable conserved motif has been the subject of extensive research during the past two decades. Increasing evidence has accumulated for the participation of noncoding RNAs in epigenetic control of elimination for a subset of IESs, and in trans-generational inheritance of alternative rearrangement patterns. This chapter summarizes our current knowledge of the structure of the germline and somatic genomes for the model species Paramecium tetraurelia, and describes the DNA cleavage and repair factors that constitute the IES excision machinery. We present an overview of the role of specialized RNA interference machineries and their associated noncoding RNAs in the control of DNA elimination. Finally, we discuss how RNA-dependent modification and/or remodeling of chromatin may guide PiggyMac to its cognate cleavage sites.

Programmed Genome Rearrangements in Tetrahymena

Ciliates are champions in programmed genome rearrangements. They carry out extensive restructuring during differentiation to drastically alter the complexity, relative copy number, and arrangement of sequences in the somatic genome. This chapter focuses on the model ciliate Tetrahymena, perhaps the simplest and best-understood ciliate studied. It summarizes past studies on various genome rearrangement processes and describes in detail the remarkable progress made in the past decade on the understanding of DNA deletion and other processes. The process occurs at thousands of specific sites to remove defined DNA segments that comprise roughly one-third of the genome including all transposons. Interestingly, this DNA rearranging process is a special form of RNA interference. It involves the production of double-stranded RNA and small RNA that guides the formation of heterochromatin. A domesticated piggyBac transposase is believed to cut off the marked chromatin, and the retained sequences are joined together through nonhomologous end-joining processes. Many of the proteins and DNA players involved have been analyzed and are described. This link provides possible explanations for the evolution, mechanism, and functional roles of the process. The article also discusses the interactions between parental and progeny somatic nuclei that affect the selection of sequences for deletion, and how the specific deletion boundaries are determined after heterochromatin marking.

DNA deletion as a mechanism for developmentally programmed centromere loss

A hallmark of active centromeres is the presence of the histone H3 variant CenH3 in the centromeric chromatin, which ensures faithful genome distribution at each cell division. A functional centromere can be inactivated, but the molecular mechanisms underlying the process of centromere inactivation remain largely unknown. Here, we describe the loss of CenH3 protein as part of a developmental program leading to the formation of the somatic nucleus in the eukaryote Paramecium. We identify two proteins whose depletion prevents developmental loss of CenH3: the domesticated transposase Pgm involved in the formation of DNA double strand cleavages and the Polycomb-like lysine methyltransferase Ezl1 necessary for trimethylation of histone H3 on lysine 9 and lysine 27. Taken together, our data support a model in which developmentally programmed centromere loss is caused by the elimination of DNA sequences associated with CenH3.

TFIIS-Dependent Non-coding Transcription Regulates Developmental Genome Rearrangements

Because of their nuclear dimorphism, ciliates provide a unique opportunity to study the role of non-coding RNAs (ncRNAs) in the communication between germline and somatic lineages. In these unicellular eukaryotes, a new somatic nucleus develops at each sexual cycle from a copy of the zygotic (germline) nucleus, while the old somatic nucleus degenerates. In the ciliate Paramecium tetraurelia, the genome is massively rearranged during this process through the reproducible elimination of repeated sequences and the precise excision of over 45,000 short, single-copy Internal Eliminated Sequences (IESs). Different types of ncRNAs resulting from genome-wide transcription were shown to be involved in the epigenetic regulation of genome rearrangements. To understand how ncRNAs are produced from the entire genome, we have focused on a homolog of the TFIIS elongation factor, which regulates RNA polymerase II transcriptional pausing. Six TFIIS-paralogs, representing four distinct families, can be found in P. tetraurelia genome. Using RNA interference, we showed that TFIIS4, which encodes a development-specific TFIIS protein, is essential for the formation of a functional somatic genome. Molecular analyses and high-throughput DNA sequencing upon TFIIS4 RNAi demonstrated that TFIIS4 is involved in all kinds of genome rearrangements, including excision of ~48% of IESs. Localization of a GFP-TFIIS4 fusion revealed that TFIIS4 appears specifically in the new somatic nucleus at an early developmental stage, before IES excision. RT-PCR experiments showed that TFIIS4 is necessary for the synthesis of IES-containing non-coding transcripts. We propose that these IES+ transcripts originate from the developing somatic nucleus and serve as pairing substrates for germline-specific short RNAs that target elimination of their homologous sequences. Our study, therefore, connects the onset of zygotic non coding transcription to the control of genome plasticity in Paramecium, and establishes for the first time a specific role of TFIIS in non-coding transcription in eukaryotes.

Paramecium tetraurelia provides an excellent model for studying the mechanisms involved in the production of non-coding transcripts and their mode of action. Different types of non-coding RNAs (ncRNAs) were shown to be implicated in the programmed DNA elimination process that occurs in this organism. At each sexual cycle, during development of the somatic nucleus from the germline nucleus, the genome is massively rearranged through the reproducible elimination of germline-specific sequences including thousands of short, single copy, non-coding Internal Eliminated Sequences (IES). Here, we demonstrate, using RNA interference, that the TFIIS4 gene encoding a development-specific homolog of RNA polymerase II elongation factor TFIIS, is indispensable for ncRNA synthesis in the new somatic nucleus. TFIIS4 depletion impairs the assembly of a functional somatic genome and affects excision of a large fraction of IESs, which leads to strong lethality in the sexual progeny. We propose that TFIIS4-dependent ncRNAs provide an important component of the molecular machinery that is responsible for developmental genome remodeling in Paramecium.

Pdsg1 and Pdsg2, novel proteins involved in developmental genome remodelling in Paramecium

The epigenetic influence of maternal cells on the development of their progeny has long been studied in various eukaryotes. Multicellular organisms usually provide their zygotes not only with nutrients but also with functional elements required for proper development, such as coding and non-coding RNAs. These maternally deposited RNAs exhibit a variety of functions, from regulating gene expression to assuring genome integrity. In ciliates, such as Paramecium these RNAs participate in the programming of large-scale genome reorganization during development, distinguishing germline-limited DNA, which is excised, from somatic-destined DNA. Only a handful of proteins playing roles in this process have been identified so far, including typical RNAi-derived factors such as Dicer-like and Piwi proteins. Here we report and characterize two novel proteins, Pdsg1 and Pdsg2 (Paramecium protein involved in Development of the Somatic Genome 1 and 2), involved in Paramecium genome reorganization. We show that these proteins are necessary for the excision of germline-limited DNA during development and the survival of sexual progeny. Knockdown of PDSG1 and PDSG2 genes affects the populations of small RNAs known to be involved in the programming of DNA elimination (scanRNAs and iesRNAs) and chromatin modification patterns during development. Our results suggest an association between RNA-mediated trans-generational epigenetic signal and chromatin modifications in the process of Paramecium genome reorganization.

Genomic characterization of variable surface antigens reveals a telomere position effect as a prerequisite for RNA interference-mediated silencing in Paramecium tetraurelia

Antigenic or phenotypic variation is a widespread phenomenon of expression of variable surface protein coats on eukaryotic microbes. To clarify the mechanism behind mutually exclusive gene expression, we characterized the genetic properties of the surface antigen multigene family in the ciliate Paramecium tetraurelia and the epigenetic factors controlling expression and silencing. Genome analysis indicated that the multigene family consists of intrachromosomal and subtelomeric genes; both classes apparently derive from different gene duplication events: whole-genome and intrachromosomal duplication. Expression analysis provides evidence for telomere position effects, because only subtelomeric genes follow mutually exclusive transcription. Microarray analysis of cultures deficient in Rdr3, an RNA-dependent RNA polymerase, in comparison to serotype-pure wild-type cultures, shows cotranscription of a subset of subtelomeric genes, indicating that the telomere position effect is due to a selective occurrence of Rdr3-mediated silencing in subtelomeric regions. We present a model of surface antigen evolution by intrachromosomal gene duplication involving the maintenance of positive selection of structurally relevant regions. Further analysis of chromosome heterogeneity shows that alternative telomere addition regions clearly affect transcription of closely related genes. Consequently, chromosome fragmentation appears to be of crucial importance for surface antigen expression and evolution. Our data suggest that RNAi-mediated control of this genetic network by trans-acting RNAs allows rapid epigenetic adaptation by phenotypic variation in combination with long-term genetic adaptation by Darwinian evolution of antigen genes.

IMPORTANCE

Alternating surface protein structures have been described for almost all eukaryotic microbes, and a broad variety of functions have been described, such as virulence factors, adhesion molecules, and molecular camouflage. Mechanisms controlling gene expression of variable surface proteins therefore represent a powerful tool for rapid phenotypic variation across kingdoms in pathogenic as well as free-living eukaryotic microbes. However, the epigenetic mechanisms controlling synchronous expression and silencing of individual genes are hardly understood. Using the ciliate Paramecium tetraurelia as a (epi)genetic model, we showed that a subtelomeric gene position effect is associated with the selective occurrence of RNAi-mediated silencing of silent surface protein genes, suggesting small interfering RNA (siRNA)-mediated epigenetic cross talks between silent and active surface antigen genes. Our integrated genomic and molecular approach discloses the correlation between gene position effects and siRNA-mediated trans-silencing, thus providing two new parameters for regulation of mutually exclusive gene expression and the genomic organization of variant gene families.

Local effect of enhancer of zeste-like reveals cooperation of epigenetic and cis-acting determinants for zygotic genome rearrangements

In the ciliate Paramecium tetraurelia, differentiation of the somatic nucleus from the zygotic nucleus is characterized by massive and reproducible deletion of transposable elements and of 45,000 short, dispersed, single-copy sequences. A specific class of small RNAs produced by the germline during meiosis, the scnRNAs, are involved in the epigenetic regulation of DNA deletion but the underlying mechanisms are poorly understood. Here, we show that trimethylation of histone H3 (H3K27me3 and H3K9me3) displays a dynamic nuclear localization that is altered when the endonuclease required for DNA elimination is depleted. We identified the putative histone methyltransferase Ezl1 necessary for H3K27me3 and H3K9me3 establishment and show that it is required for correct genome rearrangements. Genome-wide analyses show that scnRNA-mediated H3 trimethylation is necessary for the elimination of long, repeated germline DNA, while single copy sequences display differential sensitivity to depletion of proteins involved in the scnRNA pathway, Ezl1- a putative histone methyltransferase and Dcl5- a protein required for iesRNA biogenesis. Our study reveals cis-acting determinants, such as DNA length, also contribute to the definition of germline sequences to delete. We further show that precise excision of single copy DNA elements, as short as 26 bp, requires Ezl1, suggesting that development specific H3K27me3 and H3K9me3 ensure specific demarcation of very short germline sequences from the adjacent somatic sequences.

The unicellular eukaryote Paramecium tetraurelia provides an extraordinary model for studying the mechanisms involved in zygotic genome rearrangements. At each sexual cycle, differentiation of the somatic nucleus from the zygotic nucleus is characterized by extensive remodeling of the entire somatic genome, which includes the precise excision of 45,000 short noncoding germline DNA segments to reconstitute functional open reading frames. Exploiting the unique properties of the Paramecium genome, we show that the enhancer of zeste like protein Ezl1 is necessary for histone H3 trimethylation on lysines 27 and 9 and is required for the precise excision of 31,000 of these single copy, dispersed germline DNA segments that can be as short as 26 bp in length. This implies that histone marks usually associated with heterochromatin may contribute to the precise demarcation of segments that are even shorter than the length of DNA wrapped around a single nucleosome. A quantitative analysis of high throughput sequencing datasets further shows that the underlying genetic properties of the germline DNA segments might act in concert with epigenetic signals to define germline specific sequences.

Ku-mediated coupling of DNA cleavage and repair during programmed genome rearrangements in the ciliate Paramecium tetraurelia

During somatic differentiation, physiological DNA double-strand breaks (DSB) can drive programmed genome rearrangements (PGR), during which DSB repair pathways are mobilized to safeguard genome integrity. Because of their unique nuclear dimorphism, ciliates are powerful unicellular eukaryotic models to study the mechanisms involved in PGR. At each sexual cycle, the germline nucleus is transmitted to the progeny, but the somatic nucleus, essential for gene expression, is destroyed and a new somatic nucleus differentiates from a copy of the germline nucleus. In Paramecium tetraurelia, the development of the somatic nucleus involves massive PGR, including the precise elimination of at least 45,000 germline sequences (Internal Eliminated Sequences, IES). IES excision proceeds through a cut-and-close mechanism: a domesticated transposase, PiggyMac, is essential for DNA cleavage, and DSB repair at excision sites involves the Ligase IV, a specific component of the non-homologous end-joining (NHEJ) pathway. At the genome-wide level, a huge number of programmed DSBs must be repaired during this process to allow the assembly of functional somatic chromosomes. To understand how DNA cleavage and DSB repair are coordinated during PGR, we have focused on Ku, the earliest actor of NHEJ-mediated repair. Two Ku70 and three Ku80 paralogs are encoded in the genome of P. tetraurelia: Ku70a and Ku80c are produced during sexual processes and localize specifically in the developing new somatic nucleus. Using RNA interference, we show that the development-specific Ku70/Ku80c heterodimer is essential for the recovery of a functional somatic nucleus. Strikingly, at the molecular level, PiggyMac-dependent DNA cleavage is abolished at IES boundaries in cells depleted for Ku80c, resulting in IES retention in the somatic genome. PiggyMac and Ku70a/Ku80c co-purify as a complex when overproduced in a heterologous system. We conclude that Ku has been integrated in the Paramecium DNA cleavage factory, enabling tight coupling between DSB introduction and repair during PGR.

DNA double-strand breaks (DSBs) are potential threats for chromosome stability, but they are usually repaired by two major pathways, homologous recombination or non-homologous end joining (NHEJ). DSBs can also be essential during physiological processes, such as the programmed removal of germline sequences that takes place in various eukaryotes, including ciliates, during somatic differentiation. We use the ciliate Paramecium tetraurelia as a unicellular model to study how DNA breakage and DSB repair are coordinated during programmed genome rearrangements. In this organism, assembly of the somatic genome involves the elimination of ∼25% of germline DNA, including the precise excision of thousands of short Internal Eliminated Sequences (IES) scattered along germline chromosomes. A domesticated piggyBac transposase, PiggyMac, is required for double-strand DNA cleavage at IES ends and IES excision sites are very precisely repaired by the NHEJ pathway. Here, we report that a specialized Ku heterodimer, specifically expressed during programmed genome rearrangements, is an essential partner of PiggyMac and activates DNA cleavage. We propose that incorporation of DSB repair proteins in a pre-cleavage complex constitutes a safe and efficient way for Paramecium to direct thousands of programmed DSBs to the NHEJ pathway and make sure that somatic chromosomes are assembled correctly.

A forward genetic screen reveals essential and non-essential RNAi factors in Paramecium tetraurelia

In most eukaryotes, small RNA-mediated gene silencing pathways form complex interacting networks. In the ciliate Paramecium tetraurelia, at least two RNA interference (RNAi) mechanisms coexist, involving distinct but overlapping sets of protein factors and producing different types of short interfering RNAs (siRNAs). One is specifically triggered by high-copy transgenes, and the other by feeding cells with double-stranded RNA (dsRNA)-producing bacteria. In this study, we designed a forward genetic screen for mutants deficient in dsRNA-induced silencing, and a powerful method to identify the relevant mutations by whole-genome sequencing. We present a set of 47 mutant alleles for five genes, revealing two previously unknown RNAi factors: a novel Paramecium-specific protein (Pds1) and a Cid1-like nucleotidyl transferase. Analyses of allelic diversity distinguish non-essential and essential genes and suggest that the screen is saturated for non-essential, single-copy genes. We show that non-essential genes are specifically involved in dsRNA-induced RNAi while essential ones are also involved in transgene-induced RNAi. One of the latter, the RNA-dependent RNA polymerase RDR2, is further shown to be required for all known types of siRNAs, as well as for sexual reproduction. These results open the way for the dissection of the genetic complexity, interconnection, mechanisms and natural functions of RNAi pathways in P. tetraurelia.

Genome-defence small RNAs exapted for epigenetic mating-type inheritance

In the ciliate Paramecium, transposable elements and their single-copy remnants are deleted during the development of somatic macronuclei from germline micronuclei, at each sexual generation. Deletions are targeted by scnRNAs, small RNAs produced from the germ line during meiosis that first scan the maternal macronuclear genome to identify missing sequences, and then allow the zygotic macronucleus to reproduce the same deletions. Here we show that this process accounts for the maternal inheritance of mating types in Paramecium tetraurelia, a long-standing problem in epigenetics. Mating type E depends on expression of the transmembrane protein mtA, and the default type O is determined during development by scnRNA-dependent excision of the mtA promoter. In the sibling species Paramecium septaurelia, mating type O is determined by coding-sequence deletions in a different gene, mtB, which is specifically required for mtA expression. These independently evolved mechanisms suggest frequent exaptation of the scnRNA pathway to regulate cellular genes and mediate transgenerational epigenetic inheritance of essential phenotypic polymorphisms.

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

In eukaryotes, small RNAs (sRNAs) have key roles in development, gene expression regulation, and genome integrity maintenance. In ciliates, such as Paramecium, sRNAs form the heart of an epigenetic system that has evolved from core eukaryotic gene silencing components to selectively target DNA for deletion. In Paramecium, somatic genome development from the germline genome accurately eliminates the bulk of typically gene-interrupting, noncoding DNA. We have discovered an sRNA class (internal eliminated sequence [IES] sRNAs [iesRNAs]), arising later during Paramecium development, which originates from and precisely delineates germline DNA (IESs) and complements the initial sRNAs ("scan" RNAs [scnRNAs]) in targeting DNA for elimination. We show that whole-genome duplications have facilitated successive differentiations of Paramecium Dicer-like proteins, leading to cooperation between Dcl2 and Dcl3 to produce scnRNAs and to the production of iesRNAs by Dcl5. These innovations highlight the ability of sRNA systems to acquire capabilities, including those in genome development and integrity.

Epigenetics of ciliates

Research using ciliates revealed early examples of epigenetic phenomena and continues to provide novel findings. These protozoans maintain separate germline and somatic nuclei that carry transcriptionally silent and active genomes, respectively. Examining the differences in chromatin within distinct nuclei of Tetrahymena identified histone variants and established that transcriptional regulators act by modifying histones. Formation of somatic nuclei requires both transcriptional activation of silent chromatin and large-scale DNA elimination. This somatic genome remodeling is directed by homologous RNAs, acting with an RNA interference (RNAi)-related machinery. Furthermore, the content of the parental somatic genome provides a homologous template to guide this genome restructuring. The mechanisms regulating ciliate DNA rearrangements reveal the surprising power of homologous RNAs to remodel the genome and transmit information transgenerationally.

Molecular functions of small regulatory noncoding RNA

Recently, using large-scale genomic sequencing, a great number of small noncoding RNAs (ncRNA) has been discovered. Short ncRNAs can be classified into three major classes--small interfering RNA (siRNA), microRNA (miRNA), and piwi-interacting RNA (piRNA). These short ncRNAs ranging from 20 to 300 nt in size are now recognized as a new paradigm of gene regulation for controlling many biological processes. In this paper, we review the biogenesis and recent research on the functions of small regulatory non-coding RNAs and aim at understanding their important functions in living organisms.

Genomes on the edge: programmed genome instability in ciliates

Ciliates are an ancient and diverse group of microbial eukaryotes that have emerged as powerful models for RNA-mediated epigenetic inheritance. They possess extensive sets of both tiny and long noncoding RNAs that, together with a suite of proteins that includes transposases, orchestrate a broad cascade of genome rearrangements during somatic nuclear development. This Review emphasizes three important themes: the remarkable role of RNA in shaping genome structure, recent discoveries that unify many deeply diverged ciliate genetic systems, and a surprising evolutionary "sign change" in the role of small RNAs between major species groups.

The Paramecium germline genome provides a niche for intragenic parasitic DNA: evolutionary dynamics of internal eliminated sequences

Insertions of parasitic DNA within coding sequences are usually deleterious and are generally counter-selected during evolution. Thanks to nuclear dimorphism, ciliates provide unique models to study the fate of such insertions. Their germline genome undergoes extensive rearrangements during development of a new somatic macronucleus from the germline micronucleus following sexual events. In Paramecium, these rearrangements include precise excision of unique-copy Internal Eliminated Sequences (IES) from the somatic DNA, requiring the activity of a domesticated piggyBac transposase, PiggyMac. We have sequenced Paramecium tetraurelia germline DNA, establishing a genome-wide catalogue of ∼45,000 IESs, in order to gain insight into their evolutionary origin and excision mechanism. We obtained direct evidence that PiggyMac is required for excision of all IESs. Homology with known P. tetraurelia Tc1/mariner transposons, described here, indicates that at least a fraction of IESs derive from these elements. Most IES insertions occurred before a recent whole-genome duplication that preceded diversification of the P. aurelia species complex, but IES invasion of the Paramecium genome appears to be an ongoing process. Once inserted, IESs decay rapidly by accumulation of deletions and point substitutions. Over 90% of the IESs are shorter than 150 bp and present a remarkable size distribution with a ∼10 bp periodicity, corresponding to the helical repeat of double-stranded DNA and suggesting DNA loop formation during assembly of a transpososome-like excision complex. IESs are equally frequent within and between coding sequences; however, excision is not 100% efficient and there is selective pressure against IES insertions, in particular within highly expressed genes. We discuss the possibility that ancient domestication of a piggyBac transposase favored subsequent propagation of transposons throughout the germline by allowing insertions in coding sequences, a fraction of the genome in which parasitic DNA is not usually tolerated.

Ciliates are unicellular eukaryotes that rearrange their genomes at every sexual generation when a new somatic macronucleus, responsible for gene expression, develops from a copy of the germline micronucleus. In Paramecium, assembly of a functional somatic genome requires precise excision of interstitial DNA segments, the Internal Eliminated Sequences (IES), involving a domesticated piggyBac transposase, PiggyMac. To study IES origin and evolution, we sequenced germline DNA and identified 45,000 IESs. We found that at least some of these unique-copy elements are decayed Tc1/mariner transposons and that IES insertion is likely an ongoing process. After insertion, elements decay rapidly by accumulation of deletions and substitutions. The 93% of IESs shorter than 150 bp display a remarkable size distribution with a periodicity of 10 bp, the helical repeat of double-stranded DNA, consistent with the idea that evolution has only retained IESs that can form a double-stranded DNA loop during assembly of an excision complex. We propose that the ancient domestication of a piggyBac transposase, which provided a precise excision mechanism, enabled transposons to subsequently invade Paramecium coding sequences, a fraction of the genome that does not usually tolerate parasitic DNA.

RNA-guided DNA rearrangements in ciliates: is the best genome defence a good offence?

Genomes, like crazy patchwork quilts, are stitched together over evolutionary time from diverse elements, including some unwelcome invaders. To deal with parasitic mobile elements, most eukaryotes employ a genome self-defensive manoeuvre to recognise and silence such elements by homology-dependent interactions with RNA-protein complexes that alter chromatin. Ciliated protozoa employ more 'offensive' tactics by actually unstitching and reassembling their somatic genomes at every sexual generation to eliminate transposons and their remnants, using as patterns the maternal genomes that were rearranged in the previous cycle. Genetic and genomic studies of the distant relatives Paramecium and Tetrahymena have begun to reveal how such events are carried out with remarkable precision. Whole genome, non-coding transcripts from the maternal genome are compared with transcripts from the zygotic genome that are processed through an RNA interference (RNAi)-related process. Sequences found only in the latter are targeted for elimination by the resulting short 'scanRNAs' in many thousand DNA splicing reactions initiated by a domesticated transposase. The involvement of widely conserved mechanisms and protein factors clearly shows the relatedness of these phenomena to RNAi-mediated heterochromatic gene silencing. Such malleability of the genome on a generational time scale also has profound evolutionary implications, possibly including the epigenetic inheritance of acquired adaptive traits.

RNA-Mediated Epigenetic Programming of Genome Rearrangements

RNA, normally thought of as a conduit in gene expression, has a novel mode of action in ciliated protozoa. Maternal RNA templates provide both an organizing guide for DNA rearrangements and a template that can transport somatic mutations to the next generation. This opportunity for RNA-mediated genome rearrangement and DNA repair is profound in the ciliate Oxytricha, which deletes 95% of its germline genome during development in a process that severely fragments its chromosomes and then sorts and reorders the hundreds of thousands of pieces remaining. Oxytricha's somatic nuclear genome is therefore an epigenome formed through RNA templates and signals arising from the previous generation. Furthermore, this mechanism of RNA-mediated epigenetic inheritance can function across multiple generations, and the discovery of maternal template RNA molecules has revealed new biological roles for RNA and has hinted at the power of RNA molecules to sculpt genomic information in cells.

Exploiting Oxytricha trifallax nanochromosomes to screen for non-coding RNA genes

We took advantage of the unusual genomic organization of the ciliate Oxytricha trifallax to screen for eukaryotic non-coding RNA (ncRNA) genes. Ciliates have two types of nuclei: a germ line micronucleus that is usually transcriptionally inactive, and a somatic macronucleus that contains a reduced, fragmented and rearranged genome that expresses all genes required for growth and asexual reproduction. In some ciliates including Oxytricha, the macronuclear genome is particularly extreme, consisting of thousands of tiny 'nanochromosomes', each of which usually contains only a single gene. Because the organism itself identifies and isolates most of its genes on single-gene nanochromosomes, nanochromosome structure could facilitate the discovery of unusual genes or gene classes, such as ncRNA genes. Using a draft Oxytricha genome assembly and a custom-written protein-coding genefinding program, we identified a subset of nanochromosomes that lack any detectable protein-coding gene, thereby strongly enriching for nanochromosomes that carry ncRNA genes. We found only a small proportion of non-coding nanochromosomes, suggesting that Oxytricha has few independent ncRNA genes besides homologs of already known RNAs. Other than new members of known ncRNA classes including C/D and H/ACA snoRNAs, our screen identified one new family of small RNA genes, named the Arisong RNAs, which share some of the features of small nuclear RNAs.