The SUMO pathway is developmentally regulated and required for programmed DNA elimination in Paramecium tetraurelia

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.

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 Rearrangement in Ciliates: Paramecium

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

Signalling in ciliates: long- and short-range signals and molecular determinants for cellular dynamics

In ciliates, unicellular representatives of the bikont branch of evolution, inter- and intracellular signalling pathways have been analysed mainly in Paramecium tetraurelia, Paramecium multimicronucleatum and Tetrahymena thermophila and in part also in Euplotes raikovi. Electrophysiology of ciliary activity in Paramecium spp. is a most successful example. Established signalling mechanisms include plasmalemmal ion channels, recently established intracellular Ca²⁺ -release channels, as well as signalling by cyclic nucleotides and Ca²⁺ . Ca²⁺ -binding proteins (calmodulin, centrin) and Ca²⁺ -activated enzymes (kinases, phosphatases) are involved. Many organelles are endowed with specific molecules cooperating in signalling for intracellular transport and targeted delivery. Among them are recently specified soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs), monomeric GTPases, H⁺ -ATPase/pump, actin, etc. Little specification is available for some key signal transducers including mechanosensitive Ca²⁺ -channels, exocyst complexes and Ca²⁺ -sensor proteins for vesicle-vesicle/membrane interactions. The existence of heterotrimeric G-proteins and of G-protein-coupled receptors is still under considerable debate. Serine/threonine kinases dominate by far over tyrosine kinases (some predicted by phosphoproteomic analyses). Besides short-range signalling, long-range signalling also exists, e.g. as firmly installed microtubular transport rails within epigenetically determined patterns, thus facilitating targeted vesicle delivery. By envisaging widely different phenomena of signalling and subcellular dynamics, it will be shown (i) that important pathways of signalling and cellular dynamics are established already in ciliates, (ii) that some mechanisms diverge from higher eukaryotes and (iii) that considerable uncertainties still exist about some essential aspects of signalling.

Depletion of UBC9 Causes Nuclear Defects during the Vegetative and Sexual Life Cycles in Tetrahymena thermophila

Ubc9p is the sole E2-conjugating enzyme for SUMOylation, and its proper function is required for regulating key nuclear events such as transcription, DNA repair, and mitosis. In Tetrahymena thermophila, the genome is separated into a diploid germ line micronucleus (MIC) that divides by mitosis and a polyploid somatic macronucleus (MAC) that divides amitotically. This unusual nuclear organization provides novel opportunities for the study of SUMOylation and Ubc9p function. We identified the UBC9 gene and demonstrated that its complete deletion from both MIC and MAC genomes is lethal. Rescue of the lethal phenotype with a GFP-UBC9 fusion gene driven by a metallothionein promoter generated a cell line with CdCl2-dependent expression of green fluorescent protein (GFP)-Ubc9p. Depletion of Ubc9p in vegetative cells resulted in the loss of MICs, but MACs continued to divide. In contrast, expression of catalytically inactive Ubc9p resulted in the accumulation of multiple MICs. Critical roles for Ubc9p were also identified during the sexual life cycle of Tetrahymena. Cell lines that were depleted for Ubc9p did not form mating pairs and therefore could not complete any of the subsequent stages of conjugation, including meiosis and macronuclear development. Mating between cells expressing catalytically inactive Ubc9p resulted in arrest during macronuclear development, consistent with our observation that Ubc9p accumulates in the developing macronucleus.

SUMOylation is developmentally regulated and required for cell pairing during conjugation in Tetrahymena thermophila

The covalent attachment of small ubiquitin-like modifier (SUMO) to target proteins regulates numerous nuclear events in eukaryotes, including transcription, mitosis and meiosis, and DNA repair. Despite extensive interest in nuclear pathways within the field of ciliate molecular biology, there have been no investigations of the SUMO pathway in Tetrahymena. The developmental program of sexual reproduction of this organism includes cell pairing, micronuclear meiosis, and the formation of a new somatic macronucleus. We identified the Tetrahymena thermophila SMT3 (SUMO) and UBA2 (SUMO-activating enzyme) genes and demonstrated that the corresponding green fluorescent protein (GFP) tagged gene products are found predominantly in the somatic macronucleus during vegetative growth. Use of an anti-Smt3p antibody to perform immunoblot assays with whole-cell lysates during conjugation revealed a large increase in SUMOylation that peaked during formation of the new macronucleus. Immunofluorescence using the same antibody showed that the increase was localized primarily within the new macronucleus. To initiate functional analysis of the SUMO pathway, we created germ line knockout cell lines for both the SMT3 and UBA2 genes and found both are essential for cell viability. Conditional Smt3p and Uba2p cell lines were constructed by incorporation of the cadmium-inducible metallothionein promoter. Withdrawal of cadmium resulted in reduced cell growth and increased sensitivity to DNA-damaging agents. Interestingly, Smt3p and Uba2p conditional cell lines were unable to pair during sexual reproduction in the absence of cadmium, consistent with a function early in conjugation. Our studies are consistent with multiple roles for SUMOylation in Tetrahymena, including a dynamic regulation associated with the sexual life cycle.

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

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

Highly precise and developmentally programmed genome assembly in Paramecium requires ligase IV-dependent end joining

During the sexual cycle of the ciliate Paramecium, assembly of the somatic genome includes the precise excision of tens of thousands of short, non-coding germline sequences (Internal Eliminated Sequences or IESs), each one flanked by two TA dinucleotides. It has been reported previously that these genome rearrangements are initiated by the introduction of developmentally programmed DNA double-strand breaks (DSBs), which depend on the domesticated transposase PiggyMac. These DSBs all exhibit a characteristic geometry, with 4-base 5′ overhangs centered on the conserved TA, and may readily align and undergo ligation with minimal processing. However, the molecular steps and actors involved in the final and precise assembly of somatic genes have remained unknown. We demonstrate here that Ligase IV and Xrcc4p, core components of the non-homologous end-joining pathway (NHEJ), are required both for the repair of IES excision sites and for the circularization of excised IESs. The transcription of LIG4 and XRCC4 is induced early during the sexual cycle and a Lig4p-GFP fusion protein accumulates in the developing somatic nucleus by the time IES excision takes place. RNAi–mediated silencing of either gene results in the persistence of free broken DNA ends, apparently protected against extensive resection. At the nucleotide level, controlled removal of the 5′-terminal nucleotide occurs normally in LIG4-silenced cells, while nucleotide addition to the 3′ ends of the breaks is blocked, together with the final joining step, indicative of a coupling between NHEJ polymerase and ligase activities. Taken together, our data indicate that IES excision is a “cut-and-close” mechanism, which involves the introduction of initiating double-strand cleavages at both ends of each IES, followed by DSB repair via highly precise end joining. This work broadens our current view on how the cellular NHEJ pathway has cooperated with domesticated transposases for the emergence of new mechanisms involved in genome dynamics.

Double-strand breaks (DSBs) are among the most deleterious lesions that may occur on DNA. Some physiological processes, however, involve the introduction of DSBs and their subsequent repair. In the ciliate Paramecium, programmed DSBs initiate the extensive genome rearrangements that take place at each sexual cycle, during the development of the somatic nucleus. In particular, short intervening germline sequences (one every 1–2 kb along the genome) are spliced out from coding and non-coding regions. In this study, we present evidence that this process is a two-step mechanism and involves DNA cleavage at both ends of each excised sequence, followed by DSB repair. We demonstrate that cellular end-joining proteins, Ligase IV and its partner, Xrcc4p, are essential for the closure of broken excision sites, which has to be precise at the nucleotide level to allow the assembly of functional genes. This precision stands in sharp contrast to the notion that end joining is an error-prone DSB repair pathway. Therefore, Paramecium provides an excellent model for analysis of an intrinsically precise end joining pathway that has been recruited for genome-wide DSB repair.

Functional study of genes essential for autogamy and nuclear reorganization in Paramecium

Like all ciliates, Paramecium tetraurelia is a unicellular eukaryote that harbors two kinds of nuclei within its cytoplasm. At each sexual cycle, a new somatic macronucleus (MAC) develops from the germ line micronucleus (MIC) through a sequence of complex events, which includes meiosis, karyogamy, and assembly of the MAC genome from MIC sequences. The latter process involves developmentally programmed genome rearrangements controlled by noncoding RNAs and a specialized RNA interference machinery. We describe our first attempts to identify genes and biological processes that contribute to the progression of the sexual cycle. Given the high percentage of unknown genes annotated in the P. tetraurelia genome, we applied a global strategy to monitor gene expression profiles during autogamy, a self-fertilization process. We focused this pilot study on the genes carried by the largest somatic chromosome and designed dedicated DNA arrays covering 484 genes from this chromosome (1.2% of all genes annotated in the genome). Transcriptome analysis revealed four major patterns of gene expression, including two successive waves of gene induction. Functional analysis of 15 upregulated genes revealed four that are essential for vegetative growth, one of which is involved in the maintenance of MAC integrity and another in cell division or membrane trafficking. Two additional genes, encoding a MIC-specific protein and a putative RNA helicase localizing to the old and then to the new MAC, are specifically required during sexual processes. Our work provides a proof of principle that genes essential for meiosis and nuclear reorganization can be uncovered following genome-wide transcriptome analysis.

Gene expression in a paleopolyploid: a transcriptome resource for the ciliate Paramecium tetraurelia

Background

The genome of Paramecium tetraurelia, a unicellular model that belongs to the ciliate phylum, has been shaped by at least 3 successive whole genome duplications (WGD). These dramatic events, which have also been documented in plants, animals and fungi, are resolved over evolutionary time by the loss of one duplicate for the majority of genes. Thanks to a low rate of large scale genome rearrangement in Paramecium, an unprecedented large number of gene duplicates of different ages have been identified, making this organism an outstanding model to investigate the evolutionary consequences of polyploidization. The most recent WGD, with 51% of pre-duplication genes still in 2 copies, provides a snapshot of a phase of rapid gene loss that is not accessible in more ancient polyploids such as yeast.

Results

We designed a custom oligonucleotide microarray platform for P. tetraurelia genome-wide expression profiling and used the platform to measure gene expression during 1) the sexual cycle of autogamy, 2) growth of new cilia in response to deciliation and 3) biogenesis of secretory granules after massive exocytosis. Genes that are differentially expressed during these time course experiments have expression patterns consistent with a very low rate of subfunctionalization (partition of ancestral functions between duplicated genes) in particular since the most recent polyploidization event.

Conclusions

A public transcriptome resource is now available for Paramecium tetraurelia. The resource has been integrated into the ParameciumDB model organism database, providing searchable access to the data. The microarray platform, freely available through NimbleGen Systems, provides a robust, cost-effective approach for genome-wide expression profiling in P. tetraurelia. The expression data support previous studies showing that at short evolutionary times after a whole genome duplication, gene dosage balance constraints and not functional change are the major determinants of gene retention.

The conjugation-specific Die5 protein is required for development of the somatic nucleus in both Paramecium and Tetrahymena

Development in ciliated protozoa involves extensive genome reorganization within differentiating macronuclei, which shapes the somatic genome of the next vegetative generation. Major events of macronuclear differentiation include excision of internal eliminated sequences (IESs), chromosome fragmentation, and genome amplification. Proteins required for these events include those with homology throughout eukaryotes as well as proteins apparently unique to ciliates. In this study, we identified the ciliate-specific Defective in IES Excision 5 (DIE5) genes of Paramecium tetraurelia (PtDIE5) and Tetrahymena thermophila (TtDIE5) as orthologs that encode nuclear proteins expressed exclusively during development. Abrogation of PtDie5 protein (PtDie5p) function by RNA interference (RNAi)-mediated silencing or TtDie5p by gene disruption resulted in the failure of developing macronuclei to differentiate into new somatic nuclei. Tetrahymena DeltaDIE5 cells arrested late in development and failed to complete genome amplification, whereas RNAi-treated Paramecium cells highly amplified new macronuclear DNA before the failure in differentiation, findings that highlight clear differences in the biology of these distantly related species. Nevertheless, IES excision and chromosome fragmentation failed to occur in either ciliate, which strongly supports that Die5p is a critical player in these processes. In Tetrahymena, loss of zygotic expression during development was sufficient to block nuclear differentiation. This observation, together with the finding that knockdown of Die5p in Paramecium still allows genome amplification, indicates that this protein acts late in macronuclear development. Even though DNA rearrangements in these two ciliates look to be quite distinct, analysis of DIE5 establishes the action of a conserved mechanism within the genome reorganization pathway.