Epigenetics of ciliates

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.

The Conjusome-A Transient Organelle Linking Genome Rearrangements in the Parental and Developing Macronuclei

The conjusome plays an important role in the conjugation events that occur in Tetrahymena thermophila. The conjusome appears in the anterior of conjugant pairs during the early stages of new macronuclei (anlagen) development. It lacks a membrane, and is composed of a network of fibrous, electron dense material, containing background cytoplasm and ribosomes. Several proteins localize to this organelle, including Pdd1p, a chromodomain protein that participates in the formation of chromatin-containing structures in developing macronuclear anlagen, and is associated with the elimination of specific germ-line sequences from developing macronuclei. Conjugants lacking the PDD1 allele in the parental macronucleus do not show Pdd1p antibody staining in conjusomes. Investigations were performed using mutant cell lines, uniparental cytogamy and drug treatment, and show that the conjusome appears to be dependent on parental macronuclei condensation, and is a transitory organelle that traffics nuclear determinants from the parental macronucleus to the developing anlagen. These data, taken together with Pdd1p knockout experiments, suggest the conjusome is involved in the epigenetic phenomena that occur during conjugation and sexual reorganization. This is likely a conserved organelle. Conjusome-like structures were also observed in another Ciliate, Stylonichia. In general, conjusomes have features that resemble germ line P-granules.

Origins of genome-editing excisases as illuminated by the somatic genome of the ciliate Blepharisma

Massive DNA excision occurs regularly in ciliates, ubiquitous microbial eukaryotes with somatic and germline nuclei in the same cell. Tens of thousands of internally eliminated sequences (IESs) scattered throughout the ciliate germline genome are deleted during the development of the streamlined somatic genome. The genus Blepharisma represents one of the two high-level ciliate clades (subphylum Postciliodesmatophora) and, unusually, has dual pathways of somatic nuclear and genome development. This makes it ideal for investigating the functioning and evolution of these processes. Here we report the somatic genome assembly of Blepharisma stoltei strain ATCC 30299 (41 Mbp), arranged as numerous telomere-capped minichromosomal isoforms. This genome encodes eight PiggyBac transposase homologs no longer harbored by transposons. All appear subject to purifying selection, but just one, the putative IES excisase, has a complete catalytic triad. We hypothesize that PiggyBac homologs were ancestral excisases that enabled the evolution of extensive natural genome editing.

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.

Transcriptome analysis of the binucleate ciliate Tetrahymena thermophila with asynchronous nuclear cell cycles

Tetrahymena thermophila harbors two functionally and physically distinct nuclei within a shared cytoplasm. During vegetative growth, the "cell cycles" of the diploid micronucleus and polyploid macronucleus are offset. Micronuclear S phase initiates just before cytokinesis and is completed in daughter cells before onset of macronuclear DNA replication. Mitotic micronuclear division occurs mid-cell cycle, while macronuclear amitosis is coupled to cell division. Here we report the first RNA-seq cell cycle analysis of a binucleated ciliated protozoan. RNA was isolated across 1.5 vegetative cell cycles, starting with a macronuclear G1 population synchronized by centrifugal elutriation. Using MetaCycle, 3244 of the 26,000+ predicted genes were shown to be cell cycle regulated. Proteins present in both nuclei exhibit a single mRNA peak that always precedes their macronuclear function. Nucleus-limited genes, including nucleoporins and importins, are expressed before their respective nucleus-specific role. Cyclin D and A/B gene family members exhibit different expression patterns that suggest nucleus-restricted roles. Periodically expressed genes cluster into seven cyclic patterns. Four clusters have known PANTHER gene ontology terms associated with G1/S and G2/M phase. We propose that these clusters encode known and novel factors that coordinate micro- and macronuclear-specific events such as mitosis, amitosis, DNA replication, and cell division.

The micronuclear histone H3 clipping in the unicellular eukaryote Tetrahymena thermophila

Clipping of the histone H3 N-terminal tail has been implicated in multiple fundamental biological processes for a growing list of eukaryotes. H3 clipping, serving as an irreversible process to permanently remove some post-translational modifications (PTMs), may lead to noticeable changes in chromatin dynamics or gene expression. The eukaryotic model organism Tetrahymena thermophila is among the first few eukaryotes that exhibits H3 clipping activity, wherein the first six amino acids of H3 are cleaved off during vegetative growth. Clipping only occurs in the transcriptionally silent micronucleus of the binucleated T. thermophila, thus offering a unique opportunity to reveal the role of H3 clipping in epigenetic regulation. However, the physiological functions of the truncated H3 and its protease(s) for clipping remain elusive. Here, we review the major findings of H3 clipping in T. thermophila and highlight its association with histone modifications and cell cycle regulation. We also summarize the functions and mechanisms of H3 clipping in other eukaryotes, focusing on the high diversity in terms of protease families and cleavage sites. Finally, we predict several protease candidates in T. thermophila and provide insights for future studies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-022-00151-0.

A practical reference for studying meiosis in the model ciliate Tetrahymena thermophila

Meiosis is a critical cell division program that produces haploid gametes for sexual reproduction. Abnormalities in meiosis are often causes of infertility and birth defects (e.g., Down syndrome). Most organisms use a highly specialized zipper-like protein complex, the synaptonemal complex (SC), to guide and stabilize pairing of homologous chromosomes in meiosis. Although the SC is critical for meiosis in many eukaryotes, there are organisms that perform meiosis without a functional SC. However, such SC-less meiosis is poorly characterized. To understand the features of SC-less meiosis and its adaptive significance, the ciliated protozoan Tetrahymena was selected as a model. Meiosis research in Tetrahymena has revealed intriguing aspects of the regulatory programs utilized in its SC-less meiosis, yet additional efforts are needed for obtaining an in-depth comprehension of mechanisms that are associated with the absence of SC. Here, aiming at promoting a wider application of Tetrahymena for meiosis research, we introduce basic concepts and core techniques for studying meiosis in Tetrahymena and then suggest future directions for expanding the current Tetrahymena meiosis research toolbox. These methodologies could be adopted for dissecting meiosis in poorly characterized ciliates that might reveal novel features. Such data will hopefully provide insights into the function of the SC and the evolution of meiosis from a unique perspective.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-022-00149-8.

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.

Genome architecture used to supplement species delineation in two cryptic marine ciliates

The purpose of this study is to determine which taxonomic methods can elucidate clear and quantifiable differences between two cryptic ciliate species, and to test the utility of genome architecture as a new diagnostic character in the discrimination of otherwise indistinguishable taxa. Two cryptic tintinnid ciliates, Schmidingerella arcuata and Schmidingerella meunieri, are compared via traditional taxonomic characters including lorica morphometrics, ribosomal RNA (rRNA) gene barcodes and ecophysiological traits. In addition, single-cell 'omics analyses (single-cell transcriptomics and genomics) are used to elucidate and compare patterns of micronuclear genome architecture between the congeners. The results include a highly similar lorica that is larger in S. meunieri, a 0%-0.5% difference in rRNA gene barcodes, two different and nine indistinguishable growth responses among 11 prey treatments, and distinct patterns of micronuclear genomic architecture for genes detected in both ciliates. Together, these results indicate that while minor differences exist between S. arcuata and S. meunieri in common indices of taxonomic identification (i.e., lorica morphology, DNA barcode sequences and ecophysiology), differences exist in their genomic architecture, which suggests potential genetic incompatibility. Different patterns of micronuclear architecture in genes shared by both isolates also enable the design of species-specific primers, which are used in this study as unique "architectural barcodes" to demonstrate the co-occurrence of both ciliates in samples collected from a NW Atlantic estuary. These results support the utility of genomic architecture as a tool in species delineation, especially in ciliates that are cryptic or otherwise difficult to differentiate using traditional methods of identification.

Ciliate Research. From Myth to Trendsetting Science

This special issue of the Journal of Eukaryotic Microbiology (JEM) summarizes achievements obtained by generations of researchers with ciliates in widely different disciplines. In fact, ciliates range among the first cells seen under the microscope centuries ago. Their beauty made them an object of scientia amabilis and their manifold reactions made them attractive for college experiments and finally challenged causal analyses at the cellular level. Some of this work was honored by a Nobel Prize. Some observations yielded a baseline for additional novel discoveries, occasionally facilitated by specific properties of some ciliates. This also offers some advantage in the exploration of closely related parasites (malaria). Articles contributed here by colleagues from all over the world encompass a broad spectrum of ciliate life, from genetics to evolution, from molecular cell biology to ecology, from intercellular signaling to epigenetics etc. This introductory chapter, largely based on my personal perception, aims at integrating work presented in this special issue of JEM into a broader historical context up to current research.

Macronuclear development in ciliates, with a focus on nuclear architecture

Ciliates are defined by the presence of dimorphic nuclei as they have both a somatic macronucleus and germline micronucleus within each individual cell. The size and structure of both germline micronuclei and somatic macronuclei varies tremendously among ciliates. Except just after conjugation (i.e. the nuclear exchange in sexual cycle), the germline micronucleus is transcriptionally-inactive and contains canonical chromosomes that will be inherited between generations. In contrast, the transcriptionally-active macronucleus contains chromosomes that vary in size in different classes of ciliates, with some lineages having extensively-fragmented gene-sized somatic chromosomes while others contain longer multigene chromosomes. Here, we describe the variation in somatic macronuclear architecture in lineages sampled across the ciliate tree of life, specifically focusing on lineages with extensively fragmented chromosomes (e.g. the classes Phyllopharyngea and Spirotrichea). Further, we synthesize information from the literature on the development of ciliate macronuclei, focusing on changes in nuclear architecture throughout life cycles. These data highlight the tremendous diversity among ciliate nuclear cycles, extend our understanding of patterns of genome evolution, and provide insight into different germline and somatic nuclear features (e.g. nuclear structure and development) among eukaryotes.

Polycomb Repressive Complex 2 in Eukaryotes-An Evolutionary Perspective

Polycomb repressive complex 2 (PRC2) represents a group of evolutionarily conserved multi-subunit complexes that repress gene transcription by introducing trimethylation of lysine 27 on histone 3 (H3K27me3). PRC2 activity is of key importance for cell identity specification and developmental phase transitions in animals and plants. The composition, biochemistry, and developmental function of PRC2 in animal and flowering plant model species are relatively well described. Recent evidence demonstrates the presence of PRC2 complexes in various eukaryotic supergroups, suggesting conservation of the complex and its function. Here, we provide an overview of the current understanding of PRC2-mediated repression in different representatives of eukaryotic supergroups with a focus on the green lineage. By comparison of PRC2 in different eukaryotes, we highlight the possible common and diverged features suggesting evolutionary implications and outline emerging questions and directions for future research of polycomb repression and its evolution.

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.

Phylogenomics of the Epigenetic Toolkit Reveals Punctate Retention of Genes across Eukaryotes

Epigenetic processes in eukaryotes play important roles through regulation of gene expression, chromatin structure, and genome rearrangements. The roles of chromatin modification (e.g., DNA methylation and histone modification) and non-protein-coding RNAs have been well studied in animals and plants. With the exception of a few model organisms (e.g., Saccharomyces and Plasmodium), much less is known about epigenetic toolkits across the remainder of the eukaryotic tree of life. Even with limited data, previous work suggested the existence of an ancient epigenetic toolkit in the last eukaryotic common ancestor. We use PhyloToL, our taxon-rich phylogenomic pipeline, to detect homologs of epigenetic genes and evaluate their macroevolutionary patterns among eukaryotes. In addition to data from GenBank, we increase taxon sampling from understudied clades of SAR (Stramenopila, Alveolata, and Rhizaria) and Amoebozoa by adding new single-cell transcriptomes from ciliates, foraminifera, and testate amoebae. We focus on 118 gene families, 94 involved in chromatin modification and 24 involved in non-protein-coding RNA processes based on the epigenetics literature. Our results indicate 1) the presence of a large number of epigenetic gene families in the last eukaryotic common ancestor; 2) differential conservation among major eukaryotic clades, with a notable paucity of genes within Excavata; and 3) punctate distribution of epigenetic gene families between species consistent with rapid evolution leading to gene loss. Together these data demonstrate the power of taxon-rich phylogenomic studies for illuminating evolutionary patterns at scales of >1 billion years of evolution and suggest that macroevolutionary phenomena, such as genome conflict, have shaped the evolution of the eukaryotic epigenetic toolkit.

BleTIES: Annotation of natural genome editing in ciliates using long read sequencing

SUMMARY

Ciliates are single-celled eukaryotes that eliminate specific, interspersed DNA sequences (internally eliminated sequences, IESs) from their genomes during development. These are challenging to annotate and assemble because IES-containing sequences are typically much less abundant in the cell than those without, and IES sequences themselves often contain repetitive and low-complexity sequences. Long read sequencing technologies from Pacific Biosciences and Oxford Nanopore have the potential to reconstruct longer IESs than has been possible with short reads, but require a different assembly strategy. Here we present BleTIES, a software toolkit for detecting, assembling, and analyzing IESs using mapped long reads.

AVAILABILITY AND IMPLEMENTATION

BleTIES is implemented in Python 3. Source code is available at https://github.com/Swart-lab/bleties (MIT license), and also distributed via Bioconda.

SUPPLEMENTARY INFORMATION

Benchmarking of BleTIES with published sequence data.

The Evolution of the Hallmarks of Aging

The evolutionary theory of aging has set the foundations for a comprehensive understanding of aging. The biology of aging has listed and described the "hallmarks of aging," i.e., cellular and molecular mechanisms involved in human aging. The present paper is the first to infer the order of appearance of the hallmarks of bilaterian and thereby human aging throughout evolution from their presence in progressively narrower clades. Its first result is that all organisms, even non-senescent, have to deal with at least one mechanism of aging - the progressive accumulation of misfolded or unstable proteins. Due to their cumulation, these mechanisms are called "layers of aging." A difference should be made between the first four layers of unicellular aging, present in some unicellular organisms and in all multicellular opisthokonts, that stem and strike "from the inside" of individual cells and span from increasingly abnormal protein folding to deregulated nutrient sensing, and the last four layers of metacellular aging, progressively appearing in metazoans, that strike the cells of a multicellular organism "from the outside," i.e., because of other cells, and span from transcriptional alterations to the disruption of intercellular communication. The evolution of metazoans and eumetazoans probably solved the problem of aging along with the problem of unicellular aging. However, metacellular aging originates in the mechanisms by which the effects of unicellular aging are kept under control - e.g., the exhaustion of stem cells that contribute to replace damaged somatic cells. In bilaterians, additional functions have taken a toll on generally useless potentially limited lifespan to increase the fitness of organisms at the price of a progressively less efficient containment of the damage of unicellular aging. In the end, this picture suggests that geroscience should be more efficient in targeting conditions of metacellular aging rather than unicellular aging itself.

Tetrahymena meiosis: Simple yet ingenious

The presence of meiosis, which is a conserved component of sexual reproduction, across organisms from all eukaryotic kingdoms, strongly argues that sex is a primordial feature of eukaryotes. However, extant meiotic structures and processes can vary considerably between organisms. The ciliated protist Tetrahymena thermophila, which diverged from animals, plants, and fungi early in evolution, provides one example of a rather unconventional meiosis. Tetrahymena has a simpler meiosis compared with most other organisms: It lacks both a synaptonemal complex (SC) and specialized meiotic machinery for chromosome cohesion and has a reduced capacity to regulate meiotic recombination. Despite this, it also features several unique mechanisms, including elongation of the nucleus to twice the cell length to promote homologous pairing and prevent recombination between sister chromatids. Comparison of the meiotic programs of Tetrahymena and higher multicellular organisms may reveal how extant meiosis evolved from proto-meiosis.

New contribution to epigenetic studies: Isolation of micronuclei with high purity and DNA integrity in the model ciliated protist, Tetrahymena thermophila

The ciliated protist Tetrahymena thermophila is a well-known model organism with typical nuclear dimorphism containing a somatic macronucleus (MAC) and a germline micronucleus (MIC). The presence in the same cell compartment of two nuclei with distinctly different structural and functional properties provides an ideal model system to explore mechanisms of genome maintenance. Although methods for the isolation of MIC have been available for many years, cross-contamination and DNA degradation remain unresolved. Here, we describe a reliable and quick method to isolate MIC with high purity and DNA integrity in T. thermophila. Different factors are examined to optimize the MIC purification. The MAC contamination ratio in purified MIC is about 0.19% and DNA integrity of purified MIC is maintained. We also establish a more accurate method to detect the contamination rate of nuclei including microscopic observation and PCR detection. This study will facilitate further epigenetic research in Tetrahymena.

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.

A Polycomb repressive complex is required for RNAi-mediated heterochromatin formation and dynamic distribution of nuclear bodies

Polycomb group (PcG) proteins are widely utilized for transcriptional repression in eukaryotes. Here, we characterize, in the protist Tetrahymena thermophila, the EZL1 (E(z)-like 1) complex, with components conserved in metazoan Polycomb Repressive Complexes 1 and 2 (PRC1 and PRC2). The EZL1 complex is required for histone H3 K27 and K9 methylation, heterochromatin formation, transposable element control, and programmed genome rearrangement. The EZL1 complex interacts with EMA1, a helicase required for RNA interference (RNAi). This interaction is implicated in co-transcriptional recruitment of the EZL1 complex. Binding of H3K27 and H3K9 methylation by PDD1-another PcG protein interacting with the EZL1 complex-reinforces its chromatin association. The EZL1 complex is an integral part of Polycomb bodies, which exhibit dynamic distribution in Tetrahymena development: Their dispersion is driven by chromatin association, while their coalescence by PDD1, likely via phase separation. Our results provide a molecular mechanism connecting RNAi and Polycomb repression, which coordinately regulate nuclear bodies and reorganize the genome.

Programmed DNA Elimination in Vertebrates

Over the last few decades, an increasing number of vertebrate taxa have been identified that undergo programmed genome rearrangement, or programmed DNA loss, during development. In these organisms, the genome of germ cells is often reproducibly different from the genome of all other cells within the body. Although we clearly have not identified all vertebrate taxa that undergo programmed genome loss, the list of species known to undergo loss now represents ∼10% of vertebrate species, including several basally diverging lineages. Recent studies have shed new light on the targets and mechanisms of DNA loss and their association with canonical modes of DNA silencing. Ultimately, expansion of these studies into a larger collection of taxa will aid in reconstructing patterns of shared/independent ancestry of programmed DNA loss in the vertebrate lineage, as well as more recent evolutionary events that have shaped the structure and content of eliminated DNA.

Functional chimeric genes in ciliates: An instructive case from Euplotes raikovi

In ciliates, with every sexual event the transcriptionally active genes of the sub-chromosomic somatic genome that resides in the cell macronucleus are lost. They are de novo assembled starting from 'Macronuclear Destined Sequences' that arise from the fragmentation of transcriptionally silent DNA sequences of the germline chromosomic genome enclosed in the cell micronucleus. The RNA-mediated epigenetic mechanism that drives the assembly of these sequences is subject to errors which result in the formation of chimeric genes. Studying a gene family that in Euplotes raikovi controls the synthesis of protein signal pheromones responsible for a self/not-self recognition mechanism, we identified the chimeric structure of an 851-bp macronuclear gene previously known to specify soluble and membrane-bound pheromone molecules through an intron-splicing mechanism. This chimeric gene, designated mac-er-1*, conserved the native pheromone-gene structure throughout its coding and 3' regions. Instead, its 5' region is completely unrelated to the pheromone gene structure at the level of a 360-bp sequence, which derives from the assembly with a MDS destined to compound a 2417-bp gene encoding a 696-amino acid protein with unknown function. This mac-er-1* gene characterization provides further evidence that ciliates rely on functional chimeric genes that originate in non-programmed phenomena of somatic MDS recombination to increase the species genetic variability independently of gene reshuffling phenomena of the germline genome.

Comparative Studies on the Polymorphism and Copy Number Variation of mtSSU rDNA in Ciliates (Protista, Ciliophora): Implications for Phylogenetic, Environmental, and Ecological Research

While nuclear small subunit ribosomal DNA (nSSU rDNA) is the most commonly-used gene marker in studying phylogeny, ecology, abundance, and biodiversity of microbial eukaryotes, mitochondrial small subunit ribosomal DNA (mtSSU rDNA) provides an alternative. Recently, both copy number variation and sequence variation of nSSU rDNA have been demonstrated for diverse organisms, which can contribute to misinterpretation of microbiome data. Given this, we explore patterns for mtSSU rDNA among 13 selected ciliates (representing five classes), a major component of microbial eukaryotes, estimating copy number and sequence variation and comparing to that of nSSU rDNA. Our study reveals: (1) mtSSU rDNA copy number variation is substantially lower than that for nSSU rDNA; (2) mtSSU rDNA copy number ranges from 1.0 × 104 to 8.1 × 105; (3) a most common sequence of mtSSU rDNA is also found in each cell; (4) the sequence variation of mtSSU rDNA are mainly indels in poly A/T regions, and only half of species have sequence variation, which is fewer than that for nSSU rDNA; and (5) the polymorphisms between haplotypes of mtSSU rDNA would not influence the phylogenetic topology. Together, these data provide more insights into mtSSU rDNA as a powerful marker especially for microbial ecology studies.

The control of gene expression and cell identity by H3K9 trimethylation

Histone 3 lysine 9 trimethylation (H3K9me3) is a conserved histone modification that is best known for its role in constitutive heterochromatin formation and the repression of repetitive DNA elements. More recently, it has become evident that H3K9me3 is also deposited at certain loci in a tissue-specific manner and plays important roles in regulating cell identity. Notably, H3K9me3 can repress genes encoding silencing factors, pointing to a fundamental principle of repressive chromatin auto-regulation. Interestingly, recent studies have shown that H3K9me3 deposition requires protein SUMOylation in different contexts, suggesting that the SUMO pathway functions as an important module in gene silencing and heterochromatin formation. In this Review, we discuss the role of H3K9me3 in gene regulation in various systems and the molecular mechanisms that guide the silencing machinery to target loci.

N(6)-Methyladenine in eukaryotes

DNA modifications are a major form of epigenetic regulation that eukaryotic cells utilize in concert with histone modifications. While much work has been done elucidating the role of 5-methylcytosine over the past several decades, only recently has it been recognized that N(6)-methyladenine (N6-mA) is present in quantifiable and biologically active levels in the DNA of eukaryotic cells. Unlike prokaryotes which utilize N6-mA to recognize "self" from "foreign" DNA, eukaryotes have been found to use N6-mA in varying ways, from regulating transposable elements to gene regulation in response to hypoxia and stress. In this review, we examine the current state of the N6-mA in research field, and the current understanding of the biochemical mechanisms which deposit and remove N6-mA from the eukaryotic genome.

Non-coding RNA Transcription in Tetrahymena Meiotic Nuclei Requires Dedicated Mediator Complex-Associated Proteins

To preserve genome integrity, eukaryotic cells use small RNA-directed mechanisms to repress transposable elements (TEs). Paradoxically, in order to silence TEs, precursors of the small RNAs must be transcribed from TEs. However, it is still poorly understood how these precursors are transcribed from TEs under silenced conditions. In the otherwise transcriptionally silent germline micronucleus (MIC) of Tetrahymena, a burst of non-coding RNA (ncRNA) transcription occurs during meiosis. The transcripts are processed into small RNAs that serve to identify TE-related sequences for elimination. The Mediator complex (Med) has an evolutionarily conserved role for transcription by bridging gene-specific transcription factors and RNA polymerase II. Here, we report that three Med-associated factors, Emit1, Emit2, and Rib1, are required for the biogenesis of small ncRNAs. Med localizes to the MIC only during meiosis, and both Med localization and MIC ncRNA transcription require Emit1 and Emit2. In the MIC, Med occupies TE-rich pericentromeric and telomeric regions in a Rib1-dependent manner. Rib1 is dispensable for ncRNA transcription but is required for the accumulation of double-stranded ncRNAs. Nuclear and sub-nuclear localization of the three Med-associated proteins is interdependent. Hence, Emit1 and Emit2 act coordinately to import Med into the MIC, and Rib1 recruits Med to specific chromosomal locations to quantitatively or qualitatively promote the biogenesis of functional ncRNA. Our results underscore that the transcription machinery can be regulated by a set of specialized Med-associated proteins to temporally transcribe TE-related sequences from a silent genome for small RNA biogenesis and genome defense.

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.

RNAi-dependent Polycomb repression controls transposable elements in Tetrahymena

In this study, Zhao et al. show that in the protozoan Tetrahymena thermophila, germline-specific internally eliminated sequences (IESs) become transcriptionally activated in mutants deficient in the RNAi-dependent Polycomb repression pathway. Their findings suggest that the interplay between RNAi and Polycomb repression is a widely conserved phenomenon whose ancestral role is epigenetic silencing of TEs.

RNAi and Polycomb repression play evolutionarily conserved and often coordinated roles in transcriptional silencing. Here, we show that, in the protozoan Tetrahymena thermophila, germline-specific internally eliminated sequences (IESs)—many related to transposable elements (TEs)—become transcriptionally activated in mutants deficient in the RNAi-dependent Polycomb repression pathway. Germline TE mobilization also dramatically increases in these mutants. The transition from noncoding RNA (ncRNA) to mRNA production accompanies transcriptional activation of TE-related sequences and vice versa for transcriptional silencing. The balance between ncRNA and mRNA production is potentially affected by cotranscriptional processing as well as RNAi and Polycomb repression. We posit that interplay between RNAi and Polycomb repression is a widely conserved phenomenon, whose ancestral role is epigenetic silencing of TEs.

Nuclear localization signal targeting to macronucleus and micronucleus in binucleated ciliate Tetrahymena thermophila

Ciliated protozoa possess two morphologically and functionally distinct nuclei: a macronucleus (MAC) and a micronucleus (MIC). The MAC is transcriptionally active and functions in all cellular events. The MIC is transcriptionally inactive during cell growth, but functions in meiotic events to produce progeny nuclei. Thus, these two nuclei must be distinguished by the nuclear proteins required for their distinct functions during cellular events such as cell proliferation and meiosis. To understand the mechanism of the nuclear transport specific to either MAC or MIC, we identified specific nuclear localization signals (NLSs) in two MAC- and MIC-specific nuclear proteins, macronuclear histone H1 and micronuclear linker histone-like protein (Mlh1), respectively. By expressing GFP-fused fragments of these proteins in Tetrahymena thermophila cells, two distinct regions in macronuclear histone H1 protein were assigned as independent MAC-specific NLSs and two distinct regions in Mlh1 protein were assigned as independent MIC-specific NLSs. These NLSs contain several essential lysine residues responsible for the MAC- and MIC-specific nuclear transport, but neither contains any consensus sequence with known monopartite or bipartite NLSs in other model organisms. Our findings contribute to understanding how specific nuclear targeting is achieved to perform distinct nuclear functions in binucleated ciliates.

Whats, hows and whys of programmed DNA elimination in Tetrahymena

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

The protective function of noncoding DNA in genome defense of eukaryotic male germ cells

Peripheral and abundant noncoding DNA has been hypothesized to protect the genome and the central protein-coding sequences against DNA damage in somatic genome. In the cytosol, invading exogenous nucleic acids may first be deactivated by small RNAs encoded by noncoding DNA via mechanisms similar to the prokaryotic CRISPR-Cas system. In the nucleus, the radicals generated by radiation in the cytosol, radiation energy and invading exogenous nucleic acids are absorbed, blocked and/or reduced by peripheral heterochromatin, and damaged DNA in heterochromatin is removed and excluded from the nucleus to the cytoplasm through nuclear pore complexes. To further strengthen the hypothesis, this review summarizes the experimental evidence supporting the protective function of noncoding DNA in the genome of male germ cells. Based on these data, this review provides evidence supporting the protective role of noncoding DNA in the genome defense of sperm genome through similar mechanisms to those of the somatic genome.

Replicating and Cycling Stores of Information Perpetuate Life

Life is perpetuated through a single-cell bottleneck between generations in many organisms. Here, I highlight that this cell holds information in two distinct stores: in the linear DNA sequence that is replicated during cell divisions, and in the three-dimensional arrangement of molecules that can change during development but is recreated at the start of each generation. These two interdependent stores of information - one replicating with each cell division and the other cycling with a period of one generation - coevolve while perpetuating an organism. Unlike the genome sequence, the arrangement of molecules, including DNA, RNAs, proteins, sugars, lipids, etc., is not well understood. Because this arrangement and the genome sequence are transmitted together from one generation to the next, analysis of both is necessary to understand evolution and origins of inherited diseases. Recent developments suggest that tools are in place to examine how all the information to build an organism is encoded within a single cell, and how this cell code is reproduced in every generation. See also the video abstract here: https://youtu.be/IdWEL-T6TPU.

Intergenerational Transmission of Gene Regulatory Information in Caenorhabditis elegans

Epigenetic mechanisms can stably maintain gene expression states even after the initiating conditions have changed. Often epigenetic information is transmitted only to daughter cells, but evidence is emerging, in both vertebrate and invertebrate systems, for transgenerational epigenetic inheritance (TEI), the transmission of epigenetic gene regulatory information across generations. Each new description of TEI helps uncover the properties, molecular mechanisms and biological roles for TEI. The nematode Caenorhabditis elegans has been particularly instrumental in the effort to understand TEI, as multiple environmental and genetic triggers can initiate an epigenetic signal that can alter the expression of both transgenes and endogenous loci. Here, we review recent studies of TEI in C. elegans.

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.

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.

Improved methods and resources for paramecium genomics: transcription units, gene annotation and gene expression

Background

The 15 sibling species of the Paramecium aurelia cryptic species complex emerged after a whole genome duplication that occurred tens of millions of years ago. Given extensive knowledge of the genetics and epigenetics of Paramecium acquired over the last century, this species complex offers a uniquely powerful system to investigate the consequences of whole genome duplication in a unicellular eukaryote as well as the genetic and epigenetic mechanisms that drive speciation. High quality Paramecium gene models are important for research using this system. The major aim of the work reported here was to build an improved gene annotation pipeline for the Paramecium lineage.

Results

We generated oriented RNA-Seq transcriptome data across the sexual process of autogamy for the model species Paramecium tetraurelia. We determined, for the first time in a ciliate, candidate P. tetraurelia transcription start sites using an adapted Cap-Seq protocol. We developed TrUC, multi-threaded Perl software that in conjunction with TopHat mapping of RNA-Seq data to a reference genome, predicts transcription units for the annotation pipeline. We used EuGene software to combine annotation evidence. The high quality gene structural annotations obtained for P. tetraurelia were used as evidence to improve published annotations for 3 other Paramecium species. The RNA-Seq data were also used for differential gene expression analysis, providing a gene expression atlas that is more sensitive than the previously established microarray resource.

Conclusions

We have developed a gene annotation pipeline tailored for the compact genomes and tiny introns of Paramecium species. A novel component of this pipeline, TrUC, predicts transcription units using Cap-Seq and oriented RNA-Seq data. TrUC could prove useful beyond Paramecium, especially in the case of high gene density. Accurate predictions of 3′ and 5′ UTR will be particularly valuable for studies of gene expression (e.g. nucleosome positioning, identification of cis regulatory motifs). The P. tetraurelia improved transcriptome resource, gene annotations for P. tetraurelia, P. biaurelia, P. sexaurelia and P. caudatum, and Paramecium-trained EuGene configuration are available through ParameciumDB (http://paramecium.i2bc.paris-saclay.fr). TrUC software is freely distributed under a GNU GPL v3 licence (https://github.com/oarnaiz/TrUC).

Electronic supplementary material

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

Insights into transgenerational epigenetics from studies of ciliates

Epigenetics, a term with many meanings, can be broadly defined as the study of dynamic states of the genome. Ciliates, a clade of unicellular eukaryotes, can teach us about the intersection of epigenetics and evolution due to the advantages of working with cultivable ciliate lineages, plus their tendency to express extreme phenotypes such as heritable doublet morphology. Moreover, ciliates provide a powerful model for studying epigenetics given the presence of dimorphic nuclei - a somatic macronucleus and germline micronucleus - within each cell. Here, we exemplify the power of studying ciliates to learn about epigenetic phenomena. We highlight "classical" examples from morphology and physiology including cortical inheritance, mating type determination, and serotype expression. In addition, we detail molecular studies of epigenetic phenomena, including: DNA elimination; alternative processing and unscrambling; and copy number determination. Based on the implications of these studies, we discuss epigenetics as a possible functional mechanism for rapid speciation in ciliates.

Tetrahymena as a Unicellular Model Eukaryote: Genetic and Genomic Tools

Tetrahymena thermophila is a ciliate model organism whose study has led to important discoveries and insights into both conserved and divergent biological processes. In this review, we describe the tools for the use of Tetrahymena as a model eukaryote, including an overview of its life cycle, orientation to its evolutionary roots, and methodological approaches to forward and reverse genetics. Recent genomic tools have expanded Tetrahymena's utility as a genetic model system. With the unique advantages that Tetrahymena provide, we argue that it will continue to be a model organism of choice.

Transgenerational function of Tetrahymena Piwi protein Twi8p at distinctive noncoding RNA loci

Transgenerational transmission of genome-regulatory epigenetic information can determine phenotypes in the progeny of sexual reproduction. Sequence specificity of transgenerational regulation derives from small RNAs assembled into Piwi-protein complexes. Known targets of transgenerational regulation are primarily transposons and transposon-derived sequences. Here, we extend the scope of Piwi-mediated transgenerational regulation to include unique noncoding RNA loci. Ciliates such as Tetrahymena have a phenotypically silent germline micronucleus and an expressed somatic macronucleus, which is differentiated anew from a germline genome copy in sexual reproduction. We show that the nuclear-localized Tetrahymena Piwi protein Twi8p shuttles from parental to zygotic macronuclei. Genetic elimination of Twi8p has no phenotype for cells in asexual growth. On the other hand, cells lacking Twi8p arrest in sexual reproduction with zygotic nuclei that retain the germline genome structure, without the DNA elimination and fragmentation required to generate a functional macronucleus. Twi8p-bound small RNAs originate from long-noncoding RNAs with a terminal hairpin, which become detectable in the absence of Twi8p. Curiously, the loci that generate Twi8p-bound small RNAs are essential for asexual cell growth, even though Twi8 RNPs are essential only in sexual reproduction. Our findings suggest the model that Twi8 RNPs act on silent germline chromosomes to permit their conversion to expressed macronuclear chromosomes. Overall this work reveals that a Piwi protein carrying small RNAs from long-noncoding RNA loci has transgenerational function in establishing zygotic nucleus competence for gene expression.

Structure of the germline genome of Tetrahymena thermophila and relationship to the massively rearranged somatic genome

The germline genome of the binucleated ciliate Tetrahymena thermophila undergoes programmed chromosome breakage and massive DNA elimination to generate the somatic genome. Here, we present a complete sequence assembly of the germline genome and analyze multiple features of its structure and its relationship to the somatic genome, shedding light on the mechanisms of genome rearrangement as well as the evolutionary history of this remarkable germline/soma differentiation. Our results strengthen the notion that a complex, dynamic, and ongoing interplay between mobile DNA elements and the host genome have shaped Tetrahymena chromosome structure, locally and globally. Non-standard outcomes of rearrangement events, including the generation of short-lived somatic chromosomes and excision of DNA interrupting protein-coding regions, may represent novel forms of developmental gene regulation. We also compare Tetrahymena's germline/soma differentiation to that of other characterized ciliates, illustrating the wide diversity of adaptations that have occurred within this phylum.

The molecular hallmarks of epigenetic control

Over the past 20 years, breakthrough discoveries of chromatin-modifying enzymes and associated mechanisms that alter chromatin in response to physiological or pathological signals have transformed our knowledge of epigenetics from a collection of curious biological phenomena to a functionally dissected research field. Here, we provide a personal perspective on the development of epigenetics, from its historical origins to what we define as 'the modern era of epigenetic research'. We primarily highlight key molecular mechanisms of and conceptual advances in epigenetic control that have changed our understanding of normal and perturbed development.