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

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

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

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.

Paramecium Genetics, Genomics, and Evolution

The ciliate genus Paramecium served as one of the first model systems in microbial eukaryotic genetics, contributing much to the early understanding of phenomena as diverse as genome rearrangement, cryptic speciation, cytoplasmic inheritance, and endosymbiosis, as well as more recently to the evolution of mating types, introns, and roles of small RNAs in DNA processing. Substantial progress has recently been made in the area of comparative and population genomics. Paramecium species combine some of the lowest known mutation rates with some of the largest known effective populations, along with likely very high recombination rates, thereby harboring a population-genetic environment that promotes an exceptionally efficient capacity for selection. As a consequence, the genomes are extraordinarily streamlined, with very small intergenic regions combined with small numbers of tiny introns. The subject of the bulk of Paramecium research, the ancient Paramecium aurelia species complex, is descended from two whole-genome duplication events that retain high degrees of synteny, thereby providing an exceptional platform for studying the fates of duplicate genes. Despite having a common ancestor dating to several hundred million years ago, the known descendant species are morphologically indistinguishable, raising significant questions about the common view that gene duplications lead to the origins of evolutionary novelties.

A Review for the Special Issue on Paramecium as a Modern Model Organism

This review provides background and perspective for the articles contributing to the Special Issue of MDPI Micro-organisms on Paramecium as a Modern Model Organism. The six articles cover a variety of topics, each taking advantage of an important aspect of Paramecium biology: peripheral surface proteins that are developmentally regulated, endosymbiont algae and bacteria, ion channel regulation by calmodulin, regulation of cell mating reactivity and senescence, and the introns that dwell in the large genome. Each article highlights a significant aspect of Paramecium and its versatility.

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

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

A Population-Genetic Lens into the Process of Gene Loss Following Whole-Genome Duplication

Whole-genome duplications (WGDs) have occurred in many eukaryotic lineages. However, the underlying evolutionary forces and molecular mechanisms responsible for the long-term retention of gene duplicates created by WGDs are not well understood. We employ a population-genomic approach to understand the selective forces acting on paralogs and investigate ongoing duplicate-gene loss in multiple species of Paramecium that share an ancient WGD. We show that mutations that abolish protein function are more likely to be segregating in retained WGD paralogs than in single-copy genes, most likely because of ongoing nonfunctionalization post-WGD. This relaxation of purifying selection occurs in only one WGD paralog, accompanied by the gradual fixation of nonsynonymous mutations and reduction in levels of expression, and occurs over a long period of evolutionary time, “marking” one locus for future loss. Concordantly, the fitness effects of new nonsynonymous mutations and frameshift-causing indels are significantly more deleterious in the highly expressed copy compared with their paralogs with lower expression. Our results provide a novel mechanistic model of gene duplicate loss following WGDs, wherein selection acts on the sum of functional activity of both duplicate genes, allowing the two to wander in expression and functional space, until one duplicate locus eventually degenerates enough in functional efficiency or expression that its contribution to total activity is too insignificant to be retained by purifying selection. Retention of duplicates by such mechanisms predicts long times to duplicate-gene loss, which should not be falsely attributed to retention due to gain/change in function.

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

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

Broad domains of histone marks in the highly compact Paramecium macronuclear genome

The unicellular ciliate Paramecium contains a large vegetative macronucleus with several unusual characteristics, including an extremely high coding density and high polyploidy. As macronculear chromatin is devoid of heterochromatin, our study characterizes the functional epigenomic organization necessary for gene regulation and proper Pol II activity. Histone marks (H3K4me3, H3K9ac, H3K27me3) reveal no narrow peaks but broad domains along gene bodies, whereas intergenic regions are devoid of nucleosomes. Our data implicate H3K4me3 levels inside ORFs to be the main factor associated with gene expression, and H3K27me3 appears in association with H3K4me3 in plastic genes. Silent and lowly expressed genes show low nucleosome occupancy, suggesting that gene inactivation does not involve increased nucleosome occupancy and chromatin condensation. Because of a high occupancy of Pol II along highly expressed ORFs, transcriptional elongation appears to be quite different from that of other species. This is supported by missing heptameric repeats in the C-terminal domain of Pol II and a divergent elongation system. Our data imply that unoccupied DNA is the default state, whereas gene activation requires nucleosome recruitment together with broad domains of H3K4me3. In summary, gene activation and silencing in Paramecium run counter to the current understanding of chromatin biology.

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

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

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

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

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

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

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

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

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

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

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

What's Genetic Variation Got to Do with It? Starvation-Induced Self-Fertilization Enhances Survival in Paramecium

The pervasiveness of sex despite its well-known costs is a long-standing puzzle in evolutionary biology. Current explanations for the success of sex in nature largely rely on the adaptive significance of the new or rare genotypes that sex may generate. Less explored is the possibility that sex-underlying molecular mechanisms can enhance fitness and convey benefits to the individuals that bear the immediate costs of sex. Here, we show that the molecular environment associated with self-fertilization can increase stress resistance in the ciliate Paramecium tetraurelia. This advantage is independent of new genetic variation, coupled with a reduced nutritional input, and offers fresh insights into the mechanistic origin of sex. In addition to providing evidence that the molecular underpinnings of sexual reproduction and the stress response are linked in P. tetraurelia, these findings supply an integrative explanation for the persistence of self-fertilization in this ciliate.

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.

Genome plasticity in Paramecium bursaria revealed by population genomics

Background

Ciliates are an ancient and diverse eukaryotic group found in various environments. A unique feature of ciliates is their nuclear dimorphism, by which two types of nuclei, the diploid germline micronucleus (MIC) and polyploidy somatic macronucleus (MAC), are present in the same cytoplasm and serve different functions. During each sexual cycle, ciliates develop a new macronucleus in which newly fused genomes are extensively rearranged to generate functional minichromosomes. Interestingly, each ciliate species seems to have its way of processing genomes, providing a diversity of resources for studying genome plasticity and its regulation. Here, we sequenced and analyzed the macronuclear genome of different strains of Paramecium bursaria, a highly divergent species of the genus Paramecium which can stably establish endosymbioses with green algae.

Results

We assembled a high-quality macronuclear genome of P. bursaria and further refined genome annotation by comparing population genomic data. We identified several species-specific expansions in protein families and gene lineages that are potentially associated with endosymbiosis. Moreover, we observed an intensive chromosome breakage pattern that occurred during or shortly after sexual reproduction and contributed to highly variable gene dosage throughout the genome. However, patterns of copy number variation were highly correlated among genetically divergent strains, suggesting that copy number is adjusted by some regulatory mechanisms or natural selection. Further analysis showed that genes with low copy number variation among populations tended to function in basic cellular pathways, whereas highly variable genes were enriched in environmental response pathways.

Conclusions

We report programmed DNA rearrangements in the P. bursaria macronuclear genome that allow cells to adjust gene copy number globally according to individual gene functions. Our results suggest that large-scale gene copy number variation may represent an ancient mechanism for cells to adapt to different environments.

Supplementary information

The online version contains supplementary material available at 10.1186/s12915-020-00912-2.

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

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

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

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

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

ParameciumDB 2019: integrating genomic data across the genus for functional and evolutionary biology

ParameciumDB (https://paramecium.i2bc.paris-saclay.fr) is a community model organism database for the genome and genetics of the ciliate Paramecium. ParameciumDB development relies on the GMOD (www.gmod.org) toolkit. The ParameciumDB web site has been publicly available since 2006 when the P. tetraurelia somatic genome sequence was released, revealing that a series of whole genome duplications punctuated the evolutionary history of the species. The genome is linked to available genetic data and stocks. ParameciumDB has undergone major changes in its content and website since the last update published in 2011. Genomes from multiple Paramecium species, especially from the P. aurelia complex, are now included in ParameciumDB. A new modern web interface accompanies this transition to a database for the whole Paramecium genus. Gene pages have been enriched with orthology relationships, among the Paramecium species and with a panel of model organisms across the eukaryotic tree. This update also presents expert curation of Paramecium mitochondrial genomes.

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

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

Feeding exogenous dsRNA interferes with endogenous sRNA accumulation in Paramecium

Supply of exogenous dsRNA (exo-dsRNA), either by injection or by feeding, is a fast and powerful alternative to classical knockout studies. The biotechnical potential of feeding techniques is evident from the numerous studies focusing on oral administration of dsRNA to control pests and viral infection in crops/animal farming. We aimed to dissect the direct and indirect effects of exo-dsRNA feeding on the endogenous short interfering RNA (endo-siRNA) populations of the free-living ciliate Paramecium. We introduced dsRNA fragments against Dicer1 (DCR1), involved in RNA interference (RNAi) against exo- and few endo-siRNAs, and an RNAi unrelated gene, ND169. Any feeding, even the control dsRNA, diminishes genome wide the accumulation of endo-siRNAs and mRNAs. This cannot be explained by direct off-target effects and suggests mechanistic overlaps of the exo- and endo-RNAi mechanisms. Nevertheless, we observe a stronger down-regulation of mRNAs in DCR1 feeding compared with ND169 knockdown. This is likely due to the direct involvement of DCR1 in endo-siRNA accumulation. We further observed a cis-regulatory effect on mRNAs that overlap with phased endo-siRNAs. This interference of exo-dsRNA with endo-siRNAs warrants further investigations into secondary effects in target species/consumers, risk assessment of dsRNA feeding applications, and environmental pollution with dsRNA.

Comparative Genomics Analysis of Ciliates Provides Insights on the Evolutionary History Within "Nassophorea-Synhymenia-Phyllopharyngea" Assemblage

Ciliated protists (ciliates) are widely used for investigating evolution, mostly due to their successful radiation after their early evolutionary branching. In this study, we employed high-throughput sequencing technology to reveal the phylogenetic position of Synhymenia, as well as two classes Nassophorea and Phyllopharyngea, which have been a long-standing puzzle in the field of ciliate systematics and evolution. We obtained genomic and transcriptomic data from single cells of one synhymenian (Chilodontopsis depressa) and six other species of phyllopharyngeans (Chilodochona sp., Dysteria derouxi, Hartmannula sinica, Trithigmostoma cucullulus, Trochilia petrani, and Trochilia sp.). Phylogenomic analysis based on 157 orthologous genes comprising 173,835 amino acid residues revealed the affiliation of C. depressa within the class Phyllopharyngea, and the monophyly of Nassophorea, which strongly support the assignment of Synhymenia as a subclass within the class Phyllopharyngea. Comparative genomic analyses further revealed that C. depressa shares more orthologous genes with the class Nassophorea than with Phyllopharyngea, and the stop codon usage in C. depressa resembles that of Phyllopharyngea. Functional enrichment analysis demonstrated that biological pathways in C. depressa are more similar to Phyllopharyngea than Nassophorea. These results suggest that genomic and transcriptomic data can be used to provide insights into the evolutionary relationships within the "Nassophorea-Synhymenia-Phyllopharyngea" assemblage.

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

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

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

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

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

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

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

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

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

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

Paramecium Biology

Imagine that in 1678 you are Christiaan Huygens or Antonie van Leeuwenhoek seeing paramecia swim gracefully across the field of view of your new microscope. These unicellular, free-living, and swimming cells might have remained a curiosity if not for the ability of H.S. Jennings (Behavior of the lower organisms. Indiana University Press, Bloomington, 1906) and T.M. Sonneborn (Proc Natl Acad Sci USA 23:378-385, 1937) to recognize them for their behavior and genetics, both Mendelian and non-Mendelian. Following many years of painstaking work by Sonneborn and other researchers, Paramecium now serves as a modern model organism that has made specific contributions to cell and molecular biology and development. We will review the continuing usefulness and contributions of Paramecium species in this chapter.Even without a microscope, Paramecium species is visible to the naked eye because of their size (50-300 μ long). Paramecia are holotrichous ciliates, that is, unicellular organisms in the phylum Ciliophora that are covered with cilia. It was the beating of these cilia that propelled them across the slides of the first microscopes and continue to fascinate us today. Over time, Paramecium became a favorite model organism for a large variety of studies. Denis Lyn has called Paramecium the "white rat" of the Ciliophora for their manipulability and amenity to research. We will touch upon the use of Paramecium species to examine swimming behavior, ciliary structure and function, ion channel function, basal body duplication and patterning, non-Mendelian cortical inheritance, programmed DNA rearrangements, regulated secretion and exocytosis, and cell trafficking. In particular, we will focus on the use of P. tetraurelia and P. caudatum.

EuGene: An Automated Integrative Gene Finder for Eukaryotes and Prokaryotes

EuGene is an integrative gene finder applicable to both prokaryotic and eukaryotic genomes. EuGene annotated its first genome in 1999. Starting from genomic DNA sequences representing a complete genome, EuGene is able to predict the major transcript units in the genome from a variety of sources of information: statistical information, similarities with known transcripts and proteins, but also any GFF3 structured information supporting the presence or absence of specific types of elements. EuGene has been used to find genes in the plants Arabidopsis thaliana, Medicago truncatula, and Theobroma cacao; tomato, sunflower, and Rosa genomes; and in the nematode Meloidogyne incognita genome, among many others. The large fraction of plant in this list probably influenced EuGene development, especially in its capacities to withstand a genome with a large number of repeated regions and transposable elements.Depending on the sources of information used for prediction, EuGene can be considered as purely ab initio, purely similarity based, or hybrid. With the general availability of NGS-transcribed sequence data in genome projects, EuGene adopts a default hybrid behavior that strongly relies on similarity information. Initially targeted at eukaryotic genomes, EuGene has also been extended to offer integrative gene prediction for bacteria, allowing for richer and robust predictions than either purely statistical or homology-based prokaryotic gene finders.This text has been written as a practical guide that will give you the capacity to train and execute EuGene on your favorite eukaryotic genome. As the prokaryotic case is simpler and has already been described, only the main differences with the eukaryotic version were reported.

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

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

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

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

The Rab7 subfamily across Paramecium aurelia species; evidence of high conservation in sequence and function

We examined sequence conservation and signatures of selection in Rab7 proteins across 11 Paramecium aurelia species, and determined the localization patterns of two P. tetraurelia Rab7 paralogs when expressed as GFP fusions in live cells. We found that, while there is a variable number of Rab7 paralogs per genome, Rab7 genes are highly conserved in sequence and appear to be under strong purifying selection across aurelias. Additionally, and surprisingly based on earlier studies, we found that two P. tetraurelia Rab7 proteins have virtually identical localization patterns. Consistent with this, when we examined the gene family of a highly conserved Rab binding partner across aurelias (Rab-Interacting Lysosomal Protein, or RILP), we found that residues in key binding sites in RILPs were absolutely conserved in 13 of 21 proteins, representing genes from 9 of the 11 species examined. Of note, RILP gene number appears to be even more constrained than Rab7 gene number per genome. Abbreviation: WGD: Whole genome duplication.