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

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.

A SUMO E3 ligase promotes long non-coding RNA transcription to regulate small RNA-directed DNA elimination

Small RNAs target their complementary chromatin regions for gene silencing through nascent long non-coding RNAs (lncRNAs). In the ciliated protozoan Tetrahymena, the interaction between Piwi-associated small RNAs (scnRNAs) and the nascent lncRNA transcripts from the somatic genome has been proposed to induce target-directed small RNA degradation (TDSD), and scnRNAs not targeted for TDSD later target the germline-limited sequences for programmed DNA elimination. In this study, we show that the SUMO E3 ligase Ema2 is required for the accumulation of lncRNAs from the somatic genome and thus for TDSD and completing DNA elimination to make viable sexual progeny. Ema2 interacts with the SUMO E2 conjugating enzyme Ubc9 and enhances SUMOylation of the transcription regulator Spt6. We further show that Ema2 promotes the association of Spt6 and RNA polymerase II with chromatin. These results suggest that Ema2-directed SUMOylation actively promotes lncRNA transcription, which is a prerequisite for communication between the genome and small RNAs.

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.

MITE infestation accommodated by genome editing in the germline genome of the ciliate Blepharisma

During their development following sexual conjugation, ciliates excise numerous internal eliminated sequences (IESs) from a copy of the germline genome to produce the functional somatic genome. Most IESs are thought to have originated from transposons, but the presumed homology is often obscured by sequence decay. To obtain more representative perspectives on the nature of IESs and ciliate genome editing, we assembled 40,000 IESs of Blepharisma stoltei, a species belonging to a lineage (Heterotrichea) that diverged early from those of the intensively studied model ciliate species. About a quarter of IESs were short (<115 bp), largely nonrepetitive, and with a pronounced ~10 bp periodicity in length; the remainder were longer (up to 7 kbp) and nonperiodic and contained abundant interspersed repeats. Contrary to the expectation from current models, the assembled Blepharisma germline genome encodes few transposases. Instead, its most abundant repeat (8,000 copies) is a Miniature Inverted-repeat Transposable Element (MITE), apparently a deletion derivative of a germline-limited Pogo-family transposon. We hypothesize that MITEs are an important source of IESs whose proliferation is eventually self-limiting and that rather than defending the germline genomes against mobile elements, transposase domestication actually facilitates the accumulation of junk DNA.

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.

Phylogenetic profiling resolves early emergence of PRC2 and illuminates its functional core

This study strengthens the support for PRC2 emergence before the diversification of eukaryotes, detects a common presence of E(z) and ESC, indicating a conserved core, identifies diverse VEFS-Box Su(z)12 candidate proteins, and proposes a substrate specificity shift during E(z) evolution.

Polycomb repressive complex 2 (PRC2) is involved in maintaining transcriptionally silent chromatin states through methylating lysine 27 of histone H3 by the catalytic subunit enhancer of zeste [E(z)]. Here, we report the diversity of PRC2 core subunit proteins in different eukaryotic supergroups with emphasis on the early-diverged lineages and explore the molecular evolution of PRC2 subunits by phylogenetics. For the first time, we identify the putative ortholog of E(z) in Discoba, a lineage hypothetically proximal to the eukaryotic root, strongly supporting emergence of PRC2 before the diversification of eukaryotes. Analyzing 283 species, we robustly detect a common presence of E(z) and ESC, indicating a conserved functional core. Full-length Su(z)12 orthologs were identified in some lineages and species only, indicating, nonexclusively, high divergence of VEFS-Box–containing Su(z)12-like proteins, functional convergence of sequence-unrelated proteins, or Su(z)12 dispensability. Our results trace E(z) evolution within the SET-domain protein family, proposing a substrate specificity shift during E(z) evolution based on SET-domain and H3 histone interaction prediction.

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.

GC content, but not nucleosome positioning, directly contributes to intron splicing efficiency in Paramecium

Eukaryotic genes are interrupted by introns that must be accurately spliced from mRNA precursors. With an average length of 25 nt, the more than 90,000 introns of Paramecium tetraurelia stand among the shortest introns reported in eukaryotes. The mechanisms specifying the correct recognition of these tiny introns remain poorly understood. Splicing can occur cotranscriptionally, and it has been proposed that chromatin structure might influence splice site recognition. To investigate the roles of nucleosome positioning in intron recognition, we determined the nucleosome occupancy along the P. tetraurelia genome. We show that P. tetraurelia displays a regular nucleosome array with a nucleosome repeat length of ∼151 bp, among the smallest periodicities reported. Our analysis has revealed that introns are frequently associated with inter-nucleosomal DNA, pointing to an evolutionary constraint favoring introns at the AT-rich nucleosome edge sequences. Using accurate splicing efficiency data from cells depleted for nonsense-mediated decay effectors, we show that introns located at the edge of nucleosomes display higher splicing efficiency than those at the center. However, multiple regression analysis indicates that the low GC content of introns, rather than nucleosome positioning, is associated with high splicing efficiency. Our data reveal a complex link between GC content, nucleosome positioning, and intron evolution in Paramecium.

Not all Is SET for Methylation: Evolution of Eukaryotic Protein Methyltransferases

Dynamic posttranslational modifications to canonical histones that constitute the nucleosome (H2A, H2B, H3, and H4) control all aspects of enzymatic transactions with DNA. Histone methylation has been studied heavily for the past 20 years, and our mechanistic understanding of the control and function of individual methylation events on specific histone arginine and lysine residues has been greatly improved over the past decade, driven by excellent new tools and methods. Here, we will summarize what is known about the distribution and some of the functions of protein methyltransferases from all major eukaryotic supergroups. The main conclusion is that protein, and specifically histone, methylation is an ancient process. Many taxa in all supergroups have lost some subfamilies of both protein arginine methyltransferases (PRMT) and the heavily studied SET domain lysine methyltransferases (KMT). Over time, novel subfamilies, especially of SET domain proteins, arose. We use the interactions between H3K27 and H3K36 methylation as one example for the complex circuitry of histone modifications that make up the "histone code," and we discuss one recent example (Paramecium Ezl1) for how extant enzymes that may resemble more ancient SET domain KMTs are able to modify two lysine residues that have divergent functions in plants, fungi, and animals. Complexity of SET domain KMT function in the well-studied plant and animal lineages arose not only by gene duplication but also acquisition of novel DNA- and histone-binding domains in certain subfamilies.

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.

Structural insights into the interactions of Polycomb Repressive Complex 2 with chromatin

Polycomb repressive complexes are a family of chromatin modifier enzymes which are critical for regulating gene expression and maintaining cell-type identity. The reversible chemical modifications of histone H3 and H2A by the Polycomb proteins are central to its ability to function as a gene silencer. PRC2 is both a reader and writer of the tri-methylation of histone H3 lysine 27 (H3K27me3) which serves as a marker for transcription repression, and heterochromatin boundaries. Over the last few years, several studies have provided key insights into the mechanisms regulating the recruitment and activation of PRC2 at Polycomb target genes. In this review, we highlight the recent structural studies which have elucidated the roles played by Polycomb cofactor proteins in mediating crosstalk between histone post-translational modifications and the recruitment of PRC2 and the stimulation of PRC2 methyltransferase activity.

Transposon debris in ciliate genomes

Ciliate genomes are replete with "junk" DNA insertions that require DNA splicing for their removal. What is the evolutionary source of these insertions? This Primer explores a new study in PLOS Biology that supports a classic model of transposon origins.

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.

Role of Polycomb in the control of transposable elements

It is generally considered that Polycomb Repressive Complex (PRC)2 deposits the histone mark H3K27me3 on silent protein-coding genes, while transposable elements are repressed by DNA and/or H3K9 methylation. Yet, there is increasing evidence that PRC2 also targets and even silences transposable elements in representatives of several distantly related eukaryotic lineages. In plants and animals, H3K27me3 is present on transposable elements in mutants and specific cell types devoid of DNA methylation. In this Opinion, we summarize the experimental evidence for this phenomenon across the eukaryotic kingdom, and discuss its functional and evolutionary significance. We hypothesize that an ancestral role of Polycomb group (PcG) proteins was to silence transposable elements.

The unusual structure of the PiggyMac cysteine-rich domain reveals zinc finger diversity in PiggyBac-related transposases

BACKGROUND

Transposons are mobile genetic elements that colonize genomes and drive their plasticity in all organisms. DNA transposon-encoded transposases bind to the ends of their cognate transposons and catalyze their movement. In some cases, exaptation of transposon genes has allowed novel cellular functions to emerge. The PiggyMac (Pgm) endonuclease of the ciliate Paramecium tetraurelia is a domesticated transposase from the PiggyBac family. It carries a core catalytic domain typical of PiggyBac-related transposases and a short cysteine-rich domain (CRD), flanked by N- and C-terminal extensions. During sexual processes Pgm catalyzes programmed genome rearrangements (PGR) that eliminate ~ 30% of germline DNA from the somatic genome at each generation. How Pgm recognizes its DNA cleavage sites in chromatin is unclear and the structure-function relationships of its different domains have remained elusive.

RESULTS

We provide insight into Pgm structure by determining the fold adopted by its CRD, an essential domain required for PGR. Using Nuclear Magnetic Resonance, we show that the Pgm CRD binds two Zn2+ ions and forms an unusual binuclear cross-brace zinc finger, with a circularly permutated treble-clef fold flanked by two flexible arms. The Pgm CRD structure clearly differs from that of several other PiggyBac-related transposases, among which is the well-studied PB transposase from Trichoplusia ni. Instead, the arrangement of cysteines and histidines in the primary sequence of the Pgm CRD resembles that of active transposases from piggyBac-like elements found in other species and of human PiggyBac-derived domesticated transposases. We show that, unlike the PB CRD, the Pgm CRD does not bind DNA. Instead, it interacts weakly with the N-terminus of histone H3, whatever its lysine methylation state.

CONCLUSIONS

The present study points to the structural diversity of the CRD among transposases from the PiggyBac family and their domesticated derivatives, and highlights the diverse interactions this domain may establish with chromatin, from sequence-specific DNA binding to contacts with histone tails. Our data suggest that the Pgm CRD fold, whose unusual arrangement of cysteines and histidines is found in all PiggyBac-related domesticated transposases from Paramecium and Tetrahymena, was already present in the ancestral active transposase that gave rise to ciliate domesticated proteins.

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.

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 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.

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.

PRC2 is high maintenance

As the process that silences gene expression ensues during development, the stage is set for the activity of Polycomb-repressive complex 2 (PRC2) to maintain these repressed gene profiles. PRC2 catalyzes a specific histone posttranslational modification (hPTM) that fosters chromatin compaction. PRC2 also facilitates the inheritance of this hPTM through its self-contained "write and read" activities, key to preserving cellular identity during cell division. As these changes in gene expression occur without changes in DNA sequence and are inherited, the process is epigenetic in scope. Mutants of mammalian PRC2 or of its histone substrate contribute to the cancer process and other diseases, and research into these aberrant pathways is yielding viable candidates for therapeutic targeting. The effectiveness of PRC2 hinges on its being recruited to the proper chromatin sites; however, resolving the determinants to this process in the mammalian case was not straightforward and thus piqued the interest of many in the field. Here, we chronicle the latest advances toward exposing mammalian PRC2 and its high maintenance.

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.

The Polycomb protein Ezl1 mediates H3K9 and H3K27 methylation to repress transposable elements in Paramecium

In animals and plants, the H3K9me3 and H3K27me3 chromatin silencing marks are deposited by different protein machineries. H3K9me3 is catalyzed by the SET-domain SU(VAR)3-9 enzymes, while H3K27me3 is catalyzed by the SET-domain Enhancer-of-zeste enzymes, which are the catalytic subunits of Polycomb Repressive Complex 2 (PRC2). Here, we show that the Enhancer-of-zeste-like protein Ezl1 from the unicellular eukaryote Paramecium tetraurelia, which exhibits significant sequence and structural similarities with human EZH2, catalyzes methylation of histone H3 in vitro and in vivo with an apparent specificity toward K9 and K27. We find that H3K9me3 and H3K27me3 co-occur at multiple families of transposable elements in an Ezl1-dependent manner. We demonstrate that loss of these histone marks results in global transcriptional hyperactivation of transposable elements with modest effects on protein-coding gene expression. Our study suggests that although often considered functionally distinct, H3K9me3 and H3K27me3 may share a common evolutionary history as well as a common ancestral role in silencing transposable elements.

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.

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.

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.

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.

Flow cytometry sorting of nuclei enables the first global characterization of Paramecium germline DNA and transposable elements

Background

DNA elimination is developmentally programmed in a wide variety of eukaryotes, including unicellular ciliates, and leads to the generation of distinct germline and somatic genomes. The ciliate Paramecium tetraurelia harbors two types of nuclei with different functions and genome structures. The transcriptionally inactive micronucleus contains the complete germline genome, while the somatic macronucleus contains a reduced genome streamlined for gene expression. During development of the somatic macronucleus, the germline genome undergoes massive and reproducible DNA elimination events. Availability of both the somatic and germline genomes is essential to examine the genome changes that occur during programmed DNA elimination and ultimately decipher the mechanisms underlying the specific removal of germline-limited sequences.

Results

We developed a novel experimental approach that uses flow cell imaging and flow cytometry to sort subpopulations of nuclei to high purity. We sorted vegetative micronuclei and macronuclei during development of P. tetraurelia. We validated the method by flow cell imaging and by high throughput DNA sequencing. Our work establishes the proof of principle that developing somatic macronuclei can be sorted from a complex biological sample to high purity based on their size, shape and DNA content. This method enabled us to sequence, for the first time, the germline DNA from pure micronuclei and to identify novel transposable elements. Sequencing the germline DNA confirms that the Pgm domesticated transposase is required for the excision of all ~45,000 Internal Eliminated Sequences. Comparison of the germline DNA and unrearranged DNA obtained from PGM-silenced cells reveals that the latter does not provide a faithful representation of the germline genome.

Conclusions

We developed a flow cytometry-based method to purify P. tetraurelia nuclei to high purity and provided quality control with flow cell imaging and high throughput DNA sequencing. We identified 61 germline transposable elements including the first Paramecium retrotransposons. This approach paves the way to sequence the germline genomes of P. aurelia sibling species for future comparative genomic studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3713-7) contains supplementary material, which is available to authorized users.

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.

Evolution of Prdm Genes in Animals: Insights from Comparative Genomics

Prdm genes encode transcription factors with a subtype of SET domain known as the PRDF1-RIZ (PR) homology domain and a variable number of zinc finger motifs. These genes are involved in a wide variety of functions during animal development. As most Prdm genes have been studied in vertebrates, especially in mice, little is known about the evolution of this gene family. We searched for Prdm genes in the fully sequenced genomes of 93 different species representative of all the main metazoan lineages. A total of 976 Prdm genes were identified in these species. The number of Prdm genes per species ranges from 2 to 19. To better understand how the Prdm gene family has evolved in metazoans, we performed phylogenetic analyses using this large set of identified Prdm genes. These analyses allowed us to define 14 different subfamilies of Prdm genes and to establish, through ancestral state reconstruction, that 11 of them are ancestral to bilaterian animals. Three additional subfamilies were acquired during early vertebrate evolution (Prdm5, Prdm11, and Prdm17). Several gene duplication and gene loss events were identified and mapped onto the metazoan phylogenetic tree. By studying a large number of nonmetazoan genomes, we confirmed that Prdm genes likely constitute a metazoan-specific gene family. Our data also suggest that Prdm genes originated before the diversification of animals through the association of a single ancestral SET domain encoding gene with one or several zinc finger encoding genes.

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.

Cis-acting signals modulate the efficiency of programmed DNA elimination in Paramecium tetraurelia

In Paramecium, the regeneration of a functional somatic genome at each sexual event relies on the elimination of thousands of germline DNA sequences, known as Internal Eliminated Sequences (IESs), from the zygotic nuclear DNA. Here, we provide evidence that IESs’ length and sub-terminal bases jointly modulate IES excision by affecting DNA conformation in P. tetraurelia. Our study reveals an excess of complementary base pairing between IESs’ sub-terminal and contiguous sites, suggesting that IESs may form DNA loops prior to cleavage. The degree of complementary base pairing between IESs’ sub-terminal sites (termed C_in-score) is positively associated with IES length and is shaped by natural selection. Moreover, it escalates abruptly when IES length exceeds 45 nucleotides (nt), indicating that only sufficiently large IESs may form loops. Finally, we find that IESs smaller than 46 nt are favored targets of the cellular surveillance systems, presumably because of their relatively inefficient excision. Our findings extend the repertoire of cis-acting determinants for IES recognition/excision and provide unprecedented insights into the distinct selective pressures that operate on IESs and somatic DNA regions. This information potentially moves current models of IES evolution and of mechanisms of IES recognition/excision forward.

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.

On the path to genetic novelties: insights from programmed DNA elimination and RNA splicing

Understanding how genetic novelties arise is a central goal of evolutionary biology. To this end, programmed DNA elimination and RNA splicing deserve special consideration. While programmed DNA elimination reshapes genomes by eliminating chromatin during organismal development, RNA splicing rearranges genetic messages by removing intronic regions during transcription. Small RNAs help to mediate this class of sequence reorganization, which is not error-free. It is this imperfection that makes programmed DNA elimination and RNA splicing excellent candidates for generating evolutionary novelties. Leveraging a number of these two processes' mechanistic and evolutionary properties, which have been uncovered over the past years, we present recently proposed models and empirical evidence for how splicing can shape the structure of protein-coding genes in eukaryotes. We also chronicle a number of intriguing similarities between the processes of programmed DNA elimination and RNA splicing, and highlight the role that the variation in the population-genetic environment may play in shaping their target sequences.