Epigenetics of ciliates

A specialized TFIIB is required for transcription of transposon-targeting noncoding RNAs

Transposable elements (TEs) pose threats to genome stability. Therefore, small RNA-mediated heterochromatinization suppresses the transcription and hence the mobility of TEs. Paradoxically, transcription of noncoding RNA (ncRNA) from TEs is needed for the production of TE-targeting small RNAs and/or recruiting the silencing machinery to TEs. Hence, specialized RNA polymerase II (Pol II) regulators are required for such unconventional transcription in different organisms, including the developmental stage-specific Mediator complex (Med)-associated proteins in the ncRNA transcription from TE-related sequences in Tetrahymena. Yet it remains unclear how the Pol II transcriptional machinery is assembled at TE-related sequences for the ncRNA transcription. Here, we report that Pol II is regulated by Emit3, a stage-specific TFIIB-like protein specialized in TE transcription. Emit3 interacts with the TFIIH complex and localizes to TE-dense regions, especially at sites enriched with a G-rich sequence motif. Deletion of Emit3 globally abolishes Pol II-chromatin association in the meiotic nucleus, disrupts the chromatin binding of Med, and impairs the TE-biased localization of TFIIH. Conversely, Emit3's preferential localization to TE-rich loci relies in part on Med-associated proteins. These findings suggest that Emit3, TFIIH, and Med-associated proteins work together to initiate Pol II ncRNA transcription from TE-dense regions, possibly in a sequence-dependent manner.

Genome content reorganization in the non-model ciliate Chilodonella uncinata: insights into nuclear architecture, DNA content, and chromosome fragmentation during macronuclear development

Ciliates are a model lineage for studies of genome architecture given their unusual genome structures. All ciliates have both somatic macronuclei (MAC) and germline micronuclei (MIC), both of which develop from a zygotic nucleus following sex (i.e., conjugation). Nuclear developmental stages are not well documented among non-model ciliates, including Chilodonella uncinata (class Phyllopharyngea), the focus of our work. Here, we characterize nuclear architecture and genome dynamics in C. uncinata by combining 4',6-diamidino-2-phenylindole (DAPI) staining and fluorescence in situ hybridization (FISH) techniques with confocal microscopy. We developed a telomere probe for staining, which alongside DAPI allows for the identification of fragmented somatic chromosomes among the total DNA in the nuclei. We quantify both total DNA and telomere-bound signals from more than 250 nuclei sampled from 116 individual cells, and analyze changes in DNA content and nuclear architecture across Chilodonella's nuclear life cycle. Specifically, we find that MAC developmental stages in the ciliate C. uncinata are different from those reported from other ciliate species. These data provide insights into nuclear dynamics during development and enrich our understanding of genome evolution in non-model ciliates.

IMPORTANCE

Ciliates are a clade of diverse single-celled eukaryotic microorganisms that contain at least one somatic macronucleus (MAC) and germline micronucleus (MIC) within each cell/organism. Ciliates rely on complex genome rearrangements to generate somatic genomes from a zygotic nucleus. However, the development of somatic nuclei has only been documented for a few model ciliate genera, including Paramecium, Tetrahymena, and Oxytricha. Here, we study the MAC developmental process in the non-model ciliate, C. uncinata. We analyze both total DNA and the generation of gene-sized somatic chromosomes using a laser scanning confocal microscope to describe C. uncinata's nuclear life cycle. We show that DNA content changes dramatically during their life cycle and in a manner that differs from previous studies on model ciliates. Our study expands knowledge of genome dynamics in ciliates and among eukaryotes more broadly.

Functional involvement of RNAs and intrinsically disordered proteins in the assembly of heterochromatin

Heterochromatin is a highly condensed chromatin structure observed in the nuclei of eukaryotic cells. It plays a pivotal role in repressing undesired gene expression and establishing functional chromosomal domains, including centromeres and telomeres. Heterochromatin is characterized by specific histone modifications and the formation of higher-order chromatin structures mediated by proteins, such as HP1 and Polycomb repressive complexes (PRCs), which recognize the specific histone modifications. Recent studies have identified the involvement of non-coding RNAs (ncRNAs) and intrinsically disordered proteins (IDPs) in heterochromatin, leading to the proposal of a new model in which liquid-liquid phase separation (LLPS) contributes to heterochromatin formation and function. This emerging model not only broadens our understanding of heterochromatin's molecular mechanisms but also provides insights into its dynamic regulation depending on cellular context. Such advancements pave the way for exploring heterochromatin's role in genome organization and stability, as well as its implications in development and disease.

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

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

Genome content in the non-model ciliate Chilodonella uncinata: insights into nuclear architecture, gene-sized chromosomes among the total DNA in their somatic macronuclei during their development

Ciliates are a model lineage for studies of genome architecture given their unusual genome structures. All ciliates have both somatic macronuclei (MAC) and germline micronuclei (MIC), both of which develop from a zygotic nucleus following sex (i.e., conjugation). Nuclear developmental stages are not as well explored among non-model ciliate genera, including Chilodonella uncinata (Class- Phyllopharyngea), the focus of our work. Here, we characterize nuclear architecture and genome dynamics in C. uncinata by combining DAPI (4',6-diamidino-2-phenylindole) staining and fluorescence in situ hybridization (FISH) techniques with confocal microscopy. We developed a telomere probe for staining alongside DAPI, which allows for the identification of fragmented somatic chromosomes among the total DNA in the nuclei. We quantify both total DNA and telomere-bound signals to explore changes in DNA content and chromosome maturation across Chilodonella's nuclear life cycle. Specifically, we find that MAC developmental stages in the ciliate C. uncinata are different than the data reported from other ciliate species. These data provide insights into nuclear dynamics during nuclear development and enrich our understanding of genome evolution in non-model ciliates.

Nuclear dualism without extensive DNA elimination in the ciliate Loxodes magnus

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

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.

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

Supplementary data are available at Bioinformatics online.

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.

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.

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.

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.

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

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.

Transgenerational Inheritance: Parental Guidance Suggested

RNAs have been attractive candidates to transmit epigenetic information over multiple generations. In Tetrahymena, a new study demonstrates that the selective degradation of small RNAs that occurs by interaction with the parental genome can communicate altered patterns of heterochromatin formation and DNA elimination in offspring.

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.