Developmentally regulated chromosome fragmentation linked to imprecise elimination of repeated sequences in paramecia

The genome and comparative transcriptome of the euryhaline model ciliate Paramecium duboscqui reveal adaptations to environmental salinity

BACKGROUND

As a potential model organism for studies of environmental and cell biology, Paramecium duboscqui is a special euryhaline species of Paramecium that can be found in fresh, brackish, or marine water in natural salinity ranges between 0‰ and 33‰. However, the genome information as well as molecular mechanisms that account for its remarkable halotolerant traits remain extremely unknown. To characterize its genome feature, we combined PacBio and Illumina sequencing to assemble the first high-quality and near-complete macronuclear genome of P. duboscqui. Meanwhile, comparative transcriptomic profiles under different salinities gave underlying insight into the molecular mechanism of its adaptations to environmental salinity.

RESULTS

The results showed that the MAC genome of P. duboscqui comprises 160 contigs, with 113 of them possessing telomere (~ 28.82 Mb haploid genome size). Through comparative genomic analyses with the other ciliate, we found that gene families encoding transmembrane transporter proteins have been expanded in P. duboscqui, showing enormous potential in salinity adaptation. Like other Paramecium, P. duboscqui utilizes TGA as its only termination codon and has reassigned TAA and TAG to encode glutamine. P. duboscqui showed different growth rates under different salinities, with an optimum growth rate in 5‰ salinity. A comparison of the transcriptomic profiles among P. duboscqui grown under different concentrations showed that genes involved in protein folding, oxygen respiration, and glutathione-dependent detoxification were upregulated in the high-salt group, whereas genes encoding DNA-binding proteins and transcription factors were upregulated in the low-salt group, suggesting distinct mechanisms for responding to low and high salinity. Weighted gene coexpression network analysis (WGCNA) linked the hub genes expressed at 30‰ salinity to cysteine-type peptidase activity, lipid transfer, sodium hydrogen exchange, and cell division, with the hub genes expressed at 0‰ salinity involved in transmembrane transport and protein localization.

CONCLUSIONS

This study characterizes a new euryhaline model Paramecium, provides novel insights into Paramecium evolution, and describes the molecular mechanisms that have allow P. duboscqui to adapt to different osmotic environments.

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.

Uncoupling programmed DNA cleavage and repair scrambles the Paramecium somatic genome

In the ciliate Paramecium, precise excision of numerous internal eliminated sequences (IESs) from the somatic genome is essential at each sexual cycle. DNA double-strands breaks (DSBs) introduced by the PiggyMac endonuclease are repaired in a highly concerted manner by the non-homologous end joining (NHEJ) pathway, illustrated by complete inhibition of DNA cleavage when Ku70/80 proteins are missing. We show that expression of a DNA-binding-deficient Ku70 mutant (Ku70-6E) permits DNA cleavage but leads to the accumulation of unrepaired DSBs. We uncoupled DNA cleavage and repair by co-expressing wild-type and mutant Ku70. High-throughput sequencing of the developing macronucleus genome in these conditions identifies the presence of extremities healed by de novo telomere addition and numerous translocations between IES-flanking sequences. Coupling the two steps of IES excision ensures that both extremities are held together throughout the process, suggesting that DSB repair proteins are essential for assembly of a synaptic precleavage complex.

Programmed chromosome fragmentation in ciliated protozoa: multiple means to chromosome ends

SUMMARYCiliated protozoa undergo large-scale developmental rearrangement of their somatic genomes when forming a new transcriptionally active macronucleus during conjugation. This process includes the fragmentation of chromosomes derived from the germline, coupled with the efficient healing of the broken ends by de novo telomere addition. Here, we review what is known of developmental chromosome fragmentation in ciliates that have been well-studied at the molecular level (Tetrahymena, Paramecium, Euplotes, Stylonychia, and Oxytricha). These organisms differ substantially in the fidelity and precision of their fragmentation systems, as well as in the presence or absence of well-defined sequence elements that direct excision, suggesting that chromosome fragmentation systems have evolved multiple times and/or have been significantly altered during ciliate evolution. We propose a two-stage model for the evolution of the current ciliate systems, with both stages involving repetitive or transposable elements in the genome. The ancestral form of chromosome fragmentation is proposed to have been derived from the ciliate small RNA/chromatin modification process that removes transposons and other repetitive elements from the macronuclear genome during development. The evolution of this ancestral system is suggested to have potentiated its replacement in some ciliate lineages by subsequent fragmentation systems derived from mobile genetic elements.

Inter-generational nuclear crosstalk links the control of gene expression to programmed genome rearrangement during the Paramecium sexual cycle

Multinucleate cells are found in many eukaryotes, but how multiple nuclei coordinate their functions is still poorly understood. In the cytoplasm of the ciliate Paramecium tetraurelia, two micronuclei (MIC) serving sexual reproduction coexist with a somatic macronucleus (MAC) dedicated to gene expression. During sexual processes, the MAC is progressively destroyed while still ensuring transcription, and new MACs develop from copies of the zygotic MIC. Several gene clusters are successively induced and switched off before vegetative growth resumes. Concomitantly, programmed genome rearrangement (PGR) removes transposons and their relics from the new MACs. Development of the new MACs is controlled by the old MAC, since the latter expresses genes involved in PGR, including the PGM gene encoding the essential PiggyMac endonuclease that cleaves the ends of eliminated sequences. Using RNA deep sequencing and transcriptome analysis, we show that impairing PGR upregulates key known PGR genes, together with ∼600 other genes possibly also involved in PGR. Among these genes, 42% are no longer induced when no new MACs are formed, including 180 genes that are co-expressed with PGM under all tested conditions. We propose that bi-directional crosstalk between the two coexisting generations of MACs links gene expression to the progression of MAC development.

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.

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.

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.

Absolute quantification of chromosome copy numbers in the polyploid macronucleus of Tetrahymena thermophila at the single-cell level

Amitosis is widespread among eukaryotes, but the underlying mechanisms are poorly understood. The polyploid macronucleus (MAC) of unicellular ciliates divides by amitosis, making ciliates a potentially valuable model system to study this process. However, a method to accurately quantify the copy number of MAC chromosomes has not yet been established. Here we used droplet digital PCR (ddPCR) to quantify the absolute copy number of the MAC chromosomes in Tetrahymena thermophila. We first confirmed that ddPCR is a sensitive and reproducible method to determine accurate chromosome copy numbers at the single-cell level. We then used ddPCR to determine the copy number of different MAC chromosomes by analyzing individual T. thermophila cells in the G1 and the amitotic (AM) phases. The average copy number of MAC chromosomes was 90.9 at G1 phase, approximately half the number at AM phase (189.8). The copy number of each MAC chromosome varied among individual cells in G1 phase and correlated with cell size, suggesting that amitosis accompanied by unequal cytokinesis causes copy number variability. Furthermore, the fact that MAC chromosome copy number is less variable among AM-phase cells suggests that the copy number is standardized by regulating DNA replication. We also demonstrated that copy numbers differ among different MAC chromosomes and that interchromosomal variations in copy number are consistent across individual cells. Our findings demonstrate that ddPCR can be used to model amitosis in T. thermophila and possibly in other ciliates.

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.

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.

Cancer cells employ an evolutionarily conserved polyploidization program to resist therapy

Unusually large cancer cells with abnormal nuclei have been documented in the cancer literature since 1858. For more than 100 years, they have been generally disregarded as irreversibly senescent or dying cells, too morphologically misshapen and chromatin too disorganized to be functional. Cell enlargement, accompanied by whole genome doubling or more, is observed across organisms, often associated with mitigation strategies against environmental change, severe stress, or the lack of nutrients. Our comparison of the mechanisms for polyploidization in other organisms and non-transformed tissues suggest that cancer cells draw from a conserved program for their survival, utilizing whole genome doubling and pausing proliferation to survive stress. These polyaneuploid cancer cells (PACCs) are the source of therapeutic resistance, responsible for cancer recurrence and, ultimately, cancer lethality.

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.

Programmed genome rearrangements in ciliates

Ciliates are a highly divergent group of unicellular eukaryotes with separate somatic and germline genomes found in distinct dimorphic nuclei. This characteristic feature is tightly linked to extremely laborious developmentally regulated genome rearrangements in the development of a new somatic genome/nuclei following sex. The transformation from germline to soma genome involves massive DNA elimination mediated by non-coding RNAs, chromosome fragmentation, as well as DNA amplification. In this review, we discuss the similarities and differences in the genome reorganization processes of the model ciliates Paramecium and Tetrahymena (class Oligohymenophorea), and the distantly related Euplotes, Stylonychia, and Oxytricha (class Spirotrichea).

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.

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.

A meiosis-specific Spt5 homolog involved in non-coding transcription

Spt5 is a conserved and essential transcriptional regulator that binds directly to RNA polymerase and is involved in transcription elongation, polymerase pausing and various co-transcriptional processes. To investigate the role of Spt5 in non-coding transcription, we used the unicellular model Paramecium tetraurelia. In this ciliate, development is controlled by epigenetic mechanisms that use different classes of non-coding RNAs to target DNA elimination. We identified two SPT5 genes. One (STP5v) is involved in vegetative growth, while the other (SPT5m) is essential for sexual reproduction. We focused our study on SPT5m, expressed at meiosis and associated with germline nuclei during sexual processes. Upon Spt5m depletion, we observed absence of scnRNAs, piRNA-like 25 nt small RNAs produced at meiosis. The scnRNAs are a temporal copy of the germline genome and play a key role in programming DNA elimination. Moreover, Spt5m depletion abolishes elimination of all germline-limited sequences, including sequences whose excision was previously shown to be scnRNA-independent. This suggests that in addition to scnRNA production, Spt5 is involved in setting some as yet uncharacterized epigenetic information at meiosis. Our study establishes that Spt5m is crucial for developmental genome rearrangements and necessary for scnRNA production.

Multimerization properties of PiggyMac, a domesticated piggyBac transposase involved in programmed genome rearrangements

During sexual processes, the ciliate Paramecium eliminates 25–30% of germline DNA from its somatic genome. DNA elimination includes excision of ∼45 000 short, single-copy internal eliminated sequences (IESs) and depends upon PiggyMac (Pgm), a domesticated piggyBac transposase that is essential for DNA cleavage at IES ends. Pgm carries a core transposase region with a putative catalytic domain containing three conserved aspartic acids, and a downstream cysteine-rich (CR) domain. A C-terminal extension of unknown function is predicted to adopt a coiled-coil (CC) structure. To address the role of the three domains, we designed an in vivo complementation assay by expressing wild-type or mutant Pgm-GFP fusions in cells depleted for their endogenous Pgm. The DDD triad and the CR domain are essential for Pgm activity and mutations in either domain have a dominant-negative effect in wild-type cells. A mutant lacking the CC domain is partially active in the presence of limiting Pgm amounts, but inactive when Pgm is completely absent, suggesting that presence of the mutant protein increases the overall number of active complexes. We conclude that IES excision involves multiple Pgm subunits, of which at least a fraction must contain the CC domain.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

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

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

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.

Phylogenomic analysis reveals genome-wide purifying selection on TBE transposons in the ciliate Oxytricha

Background

Transposable elements are a major player contributing to genetic variation and shaping genome evolution. Multiple independent transposon domestication events have occurred in ciliates, recruiting transposases to key roles in cellular processes. In the ciliate Oxytricha trifallax, the telomere-bearing elements (TBE), a Tc1/mariner transposon, occupy a significant portion of the germline genome and are involved in programmed genome rearrangements that produce a transcriptionally active somatic nucleus from a copy of the germline nucleus during development.

Results

Here we provide a thorough characterization of the distribution and sequences of TBE transposons in the Oxytricha germline genome. We annotate more than 10,000 complete and 24,000 partial TBE sequences. TBEs cluster into four major families and display a preference for either insertion into DNA segments that are retained in the somatic genome or their maintenance at such sites. The three TBE-encoded genes in all four families display dN/dS ratios much lower than 1, suggesting genome-wide purifying selection. We also identify TBE homologs in other ciliate species for phylogenomic analysis.

Conclusions

This paper provides genome-wide characterization of a major class of ciliate transposons. Phylogenomic analysis reveals selective constraints on transposon-encoded genes, shedding light on the evolution and domesticated functions of these transposons.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-016-0057-9) contains supplementary material, which is available to authorized users.

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.

Genomic characterization of variable surface antigens reveals a telomere position effect as a prerequisite for RNA interference-mediated silencing in Paramecium tetraurelia

Antigenic or phenotypic variation is a widespread phenomenon of expression of variable surface protein coats on eukaryotic microbes. To clarify the mechanism behind mutually exclusive gene expression, we characterized the genetic properties of the surface antigen multigene family in the ciliate Paramecium tetraurelia and the epigenetic factors controlling expression and silencing. Genome analysis indicated that the multigene family consists of intrachromosomal and subtelomeric genes; both classes apparently derive from different gene duplication events: whole-genome and intrachromosomal duplication. Expression analysis provides evidence for telomere position effects, because only subtelomeric genes follow mutually exclusive transcription. Microarray analysis of cultures deficient in Rdr3, an RNA-dependent RNA polymerase, in comparison to serotype-pure wild-type cultures, shows cotranscription of a subset of subtelomeric genes, indicating that the telomere position effect is due to a selective occurrence of Rdr3-mediated silencing in subtelomeric regions. We present a model of surface antigen evolution by intrachromosomal gene duplication involving the maintenance of positive selection of structurally relevant regions. Further analysis of chromosome heterogeneity shows that alternative telomere addition regions clearly affect transcription of closely related genes. Consequently, chromosome fragmentation appears to be of crucial importance for surface antigen expression and evolution. Our data suggest that RNAi-mediated control of this genetic network by trans-acting RNAs allows rapid epigenetic adaptation by phenotypic variation in combination with long-term genetic adaptation by Darwinian evolution of antigen genes.

IMPORTANCE

Alternating surface protein structures have been described for almost all eukaryotic microbes, and a broad variety of functions have been described, such as virulence factors, adhesion molecules, and molecular camouflage. Mechanisms controlling gene expression of variable surface proteins therefore represent a powerful tool for rapid phenotypic variation across kingdoms in pathogenic as well as free-living eukaryotic microbes. However, the epigenetic mechanisms controlling synchronous expression and silencing of individual genes are hardly understood. Using the ciliate Paramecium tetraurelia as a (epi)genetic model, we showed that a subtelomeric gene position effect is associated with the selective occurrence of RNAi-mediated silencing of silent surface protein genes, suggesting small interfering RNA (siRNA)-mediated epigenetic cross talks between silent and active surface antigen genes. Our integrated genomic and molecular approach discloses the correlation between gene position effects and siRNA-mediated trans-silencing, thus providing two new parameters for regulation of mutually exclusive gene expression and the genomic organization of variant gene families.

Paramecium tetraurelia chromatin assembly factor-1-like protein PtCAF-1 is involved in RNA-mediated control of DNA elimination

Genome-wide DNA remodelling in the ciliate Paramecium is ensured by RNA-mediated trans-nuclear crosstalk between the germline and the somatic genomes during sexual development. The rearrangements include elimination of transposable elements, minisatellites and tens of thousands non-coding elements called internally eliminated sequences (IESs). The trans-nuclear genome comparison process employs a distinct class of germline small RNAs (scnRNAs) that are compared against the parental somatic genome to select the germline-specific subset of scnRNAs that subsequently target DNA elimination in the progeny genome. Only a handful of proteins involved in this process have been identified so far and the mechanism of DNA targeting is unknown. Here we describe chromatin assembly factor-1-like protein (PtCAF-1), which we show is required for the survival of sexual progeny and localizes first in the parental and later in the newly developing macronucleus. Gene silencing shows that PtCAF-1 is required for the elimination of transposable elements and a subset of IESs. PTCAF-1 depletion also impairs the selection of germline-specific scnRNAs during development. We identify specific histone modifications appearing during Paramecium development which are strongly reduced in PTCAF-1 depleted cells. Our results demonstrate the importance of PtCAF-1 for the epigenetic trans-nuclear cross-talk mechanism.

Genome-wide analysis of genetic and epigenetic control of programmed DNA deletion

During the development of the somatic genome from the Paramecium germline genome the bulk of the copies of ∼45 000 unique, internal eliminated sequences (IESs) are deleted. IES targeting is facilitated by two small RNA (sRNA) classes: scnRNAs, which relay epigenetic information from the parental nucleus to the developing nucleus, and iesRNAs, which are produced and used in the developing nucleus. Why only certain IESs require sRNAs for their removal has been enigmatic. By analyzing the silencing effects of three genes: PGM (responsible for DNA excision), DCL2/3 (scnRNA production) and DCL5 (iesRNA production), we identify key properties required for IES elimination. Based on these results, we propose that, depending on the exact combination of their lengths and end bases, some IESs are less efficiently recognized or excised and have a greater requirement for targeting by scnRNAs and iesRNAs. We suggest that the variation in IES retention following silencing of DCL2/3 is not primarily due to scnRNA density, which is comparatively uniform relative to IES retention, but rather the genetic properties of IESs. Taken together, our analyses demonstrate that in Paramecium the underlying genetic properties of developmentally deleted DNA sequences are essential in determining the sensitivity of these sequences to epigenetic control.

Epigenetics of ciliates

Research using ciliates revealed early examples of epigenetic phenomena and continues to provide novel findings. These protozoans maintain separate germline and somatic nuclei that carry transcriptionally silent and active genomes, respectively. Examining the differences in chromatin within distinct nuclei of Tetrahymena identified histone variants and established that transcriptional regulators act by modifying histones. Formation of somatic nuclei requires both transcriptional activation of silent chromatin and large-scale DNA elimination. This somatic genome remodeling is directed by homologous RNAs, acting with an RNA interference (RNAi)-related machinery. Furthermore, the content of the parental somatic genome provides a homologous template to guide this genome restructuring. The mechanisms regulating ciliate DNA rearrangements reveal the surprising power of homologous RNAs to remodel the genome and transmit information transgenerationally.

Non-coding RNAs mediate the rearrangements of genomic DNA in ciliates

Most eukaryotes employ a variety of mechanisms to defend the integrity of their genome by recognizing and silencing parasitic mobile nucleic acids. However, recent studies have shown that genomic DNA undergoes extensive rearrangements, including DNA elimination, fragmentation, and unscrambling, during the sexual reproduction of ciliated protozoa. Non-coding RNAs have been identified to program and regulate genome rearrangement events. In Paramecium and Tetrahymena, scan RNAs (scnRNAs) are produced from micronuclei and transported to vegetative macronuclei, in which scnRNA elicits the elimination of cognate genomic DNA. In contrast, Piwi-interacting RNAs (piRNAs) in Oxytricha enable the retention of genomic DNA that exhibits sequence complementarity in macronuclei. An RNA interference (RNAi)-like mechanism has been found to direct these genomic rearrangements. Furthermore, in Oxytricha, maternal RNA templates can guide the unscrambling process of genomic DNA. The non-coding RNA-directed genome rearrangements may have profound evolutionary implications, for example, eliciting the multigenerational inheritance of acquired adaptive traits.

Transposon domestication versus mutualism in ciliate genome rearrangements

Ciliated protists rearrange their genomes dramatically during nuclear development via chromosome fragmentation and DNA deletion to produce a trimmer and highly reorganized somatic genome. The deleted portion of the genome includes potentially active transposons or transposon-like sequences that reside in the germline. Three independent studies recently showed that transposase proteins of the DDE/DDD superfamily are indispensible for DNA processing in three distantly related ciliates. In the spirotrich Oxytricha trifallax, high copy-number germline-limited transposons mediate their own excision from the somatic genome but also contribute to programmed genome rearrangement through a remarkable transposon mutualism with the host. By contrast, the genomes of two oligohymenophorean ciliates, Tetrahymena thermophila and Paramecium tetraurelia, encode homologous PiggyBac-like transposases as single-copy genes in both their germline and somatic genomes. These domesticated transposases are essential for deletion of thousands of different internal sequences in these species. This review contrasts the events underlying somatic genome reduction in three different ciliates and considers their evolutionary origins and the relationships among their distinct mechanisms for genome remodeling.

A model for the evolution of extremely fragmented macronuclei in ciliates

While all ciliates possess nuclear dimorphism, several ciliates - like those in the classes Phyllopharyngea, Spirotrichea, and Armophorea - have an extreme macronuclear organization. Their extensively fragmented macronuclei contain upwards of 20,000 chromosomes, each with upwards of thousands of copies. These features have evolved independently on multiple occasions throughout ciliate evolutionary history, and currently no models explain these structures in an evolutionary context. In this paper, we propose that competition between two forces - the limitation and avoidance of chromosomal imbalances as a ciliate undergoes successive asexual divisions, and the costs of replicating massive genomes - is sufficient to explain this particular nuclear structure. We present a simulation of ciliate cell evolution under control of these forces, allowing certain features of the population to change over time. Over a wide range of parameters, we observe the repeated emergence of this unusual genomic organization found in nature. Although much remains to be understood about the evolution of macronuclear genome organization, our results show that the proposed model is a plausible explanation for the emergence of these extremely fragmented, highly polyploid genomes.

Genomes on the edge: programmed genome instability in ciliates

Ciliates are an ancient and diverse group of microbial eukaryotes that have emerged as powerful models for RNA-mediated epigenetic inheritance. They possess extensive sets of both tiny and long noncoding RNAs that, together with a suite of proteins that includes transposases, orchestrate a broad cascade of genome rearrangements during somatic nuclear development. This Review emphasizes three important themes: the remarkable role of RNA in shaping genome structure, recent discoveries that unify many deeply diverged ciliate genetic systems, and a surprising evolutionary "sign change" in the role of small RNAs between major species groups.

The Oxytricha trifallax macronuclear genome: a complex eukaryotic genome with 16,000 tiny chromosomes

With more chromosomes than any other sequenced genome, the macronuclear genome of Oxytricha trifallax has a unique and complex architecture, including alternative fragmentation and predominantly single-gene chromosomes.

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor “silent” germline micronuclear genome by a process of “unscrambling” and fragmentation. The tiny macronuclear “nanochromosomes” typically encode single, protein-coding genes (a small portion, 10%, encode 2–8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.

The macronuclear genome of the ciliate Oxytricha trifallax, contained in its somatic nucleus, has a unique genome architecture. Unlike its diploid germline genome, which is transcriptionally inactive during normal cellular growth, the macronuclear genome is fragmented into at least 16,000 tiny (∼3.2 kb mean length) chromosomes, most of which encode single actively transcribed genes and are differentially amplified to a few thousand copies each. The smallest chromosome is just 469 bp, while the largest is 66 kb and encodes a single enormous protein. We found considerable variation in the genome, including frequent alternative fragmentation patterns, generating chromosome isoforms with shared sequence. We also found limited variation in chromosome amplification levels, though insufficient to explain mRNA transcript level variation. Another remarkable feature of Oxytricha's macronuclear genome is its inordinate fondness for telomeres. In conjunction with its possession of tens of millions of chromosome-ending telomeres per macronucleus, we show that Oxytricha has evolved multiple putative telomere-binding proteins. In addition, we identified two new domesticated transposase-like protein classes that we propose may participate in the process of genome rearrangement. The macronuclear genome now provides a crucial resource for ongoing studies of genome rearrangement processes that use Oxytricha as an experimental or comparative model.

The Paramecium germline genome provides a niche for intragenic parasitic DNA: evolutionary dynamics of internal eliminated sequences

Insertions of parasitic DNA within coding sequences are usually deleterious and are generally counter-selected during evolution. Thanks to nuclear dimorphism, ciliates provide unique models to study the fate of such insertions. Their germline genome undergoes extensive rearrangements during development of a new somatic macronucleus from the germline micronucleus following sexual events. In Paramecium, these rearrangements include precise excision of unique-copy Internal Eliminated Sequences (IES) from the somatic DNA, requiring the activity of a domesticated piggyBac transposase, PiggyMac. We have sequenced Paramecium tetraurelia germline DNA, establishing a genome-wide catalogue of ∼45,000 IESs, in order to gain insight into their evolutionary origin and excision mechanism. We obtained direct evidence that PiggyMac is required for excision of all IESs. Homology with known P. tetraurelia Tc1/mariner transposons, described here, indicates that at least a fraction of IESs derive from these elements. Most IES insertions occurred before a recent whole-genome duplication that preceded diversification of the P. aurelia species complex, but IES invasion of the Paramecium genome appears to be an ongoing process. Once inserted, IESs decay rapidly by accumulation of deletions and point substitutions. Over 90% of the IESs are shorter than 150 bp and present a remarkable size distribution with a ∼10 bp periodicity, corresponding to the helical repeat of double-stranded DNA and suggesting DNA loop formation during assembly of a transpososome-like excision complex. IESs are equally frequent within and between coding sequences; however, excision is not 100% efficient and there is selective pressure against IES insertions, in particular within highly expressed genes. We discuss the possibility that ancient domestication of a piggyBac transposase favored subsequent propagation of transposons throughout the germline by allowing insertions in coding sequences, a fraction of the genome in which parasitic DNA is not usually tolerated.

Ciliates are unicellular eukaryotes that rearrange their genomes at every sexual generation when a new somatic macronucleus, responsible for gene expression, develops from a copy of the germline micronucleus. In Paramecium, assembly of a functional somatic genome requires precise excision of interstitial DNA segments, the Internal Eliminated Sequences (IES), involving a domesticated piggyBac transposase, PiggyMac. To study IES origin and evolution, we sequenced germline DNA and identified 45,000 IESs. We found that at least some of these unique-copy elements are decayed Tc1/mariner transposons and that IES insertion is likely an ongoing process. After insertion, elements decay rapidly by accumulation of deletions and substitutions. The 93% of IESs shorter than 150 bp display a remarkable size distribution with a periodicity of 10 bp, the helical repeat of double-stranded DNA, consistent with the idea that evolution has only retained IESs that can form a double-stranded DNA loop during assembly of an excision complex. We propose that the ancient domestication of a piggyBac transposase, which provided a precise excision mechanism, enabled transposons to subsequently invade Paramecium coding sequences, a fraction of the genome that does not usually tolerate parasitic DNA.

Transposon Invasion of the Paramecium Germline Genome Countered by a Domesticated PiggyBac Transposase and the NHEJ Pathway

Sequences related to transposons constitute a large fraction of extant genomes, but insertions within coding sequences have generally not been tolerated during evolution. Thanks to their unique nuclear dimorphism and to their original mechanism of programmed DNA elimination from their somatic nucleus (macronucleus), ciliates are emerging model organisms for the study of the impact of transposable elements on genomes. The germline genome of the ciliate Paramecium, located in its micronucleus, contains thousands of short intervening sequences, the IESs, which interrupt 47% of genes. Recent data provided support to the hypothesis that an evolutionary link exists between Paramecium IESs and Tc1/mariner transposons. During development of the macronucleus, IESs are excised precisely thanks to the coordinated action of PiggyMac, a domesticated piggyBac transposase, and of the NHEJ double-strand break repair pathway. A PiggyMac homolog is also required for developmentally programmed DNA elimination in another ciliate, Tetrahymena. Here, we present an overview of the life cycle of these unicellular eukaryotes and of the developmentally programmed genome rearrangements that take place at each sexual cycle. We discuss how ancient domestication of a piggyBac transposase might have allowed Tc1/mariner elements to spread throughout the germline genome of Paramecium, without strong counterselection against insertion within genes.

An essential role for the DNA breakage-repair protein Ku80 in programmed DNA rearrangements in Tetrahymena thermophila

Programmed DNA rearrangements are important processes present in many organisms. In the ciliated protozoan Tetrahymena thermophila, DNA rearrangements occur during the sexual conjugation process and lead to the deletion of thousands of specific DNA segments and fragmentation of the chromosomes. In this study, we found that the Ku80 homologue, a conserved component of the nonhomologous end-joining process of DNA repair, was essential for these two processes. During conjugation, TKU80 was highly expressed and localized to the new macronucleus, where DNA rearrangements occur. Homokaryon TKU80-knockout mutants are unable to complete conjugation and produce progeny and are arrested at the two-micronuclei/two-macronuclei stage. Analysis of their DNA revealed failure to complete DNA deletion. However, the DNA-cutting step appeared to have occurred, as evidenced by the presence of circularized excised DNA. Moreover, chromosome breakage or de novo telomere addition was affected. The mutant appears to accumulate free DNA ends detectable by terminal deoxynucleotidyl transferase dUTP nick end labeling assays that led to the degradation of most DNA in the developing macronucleus. These findings suggest that Tku80p may serve an end-protective role after DNA cleavage has occurred. Unexpectedly, the large heterochromatin structures that normally associate with DNA rearrangements failed to form without TKU80. Together the results suggest multiple roles for Tku80p and indicate that a Ku-dependent DNA-repair pathway is involved in programmed DNA rearrangements in Tetrahymena.

RNA-guided DNA rearrangements in ciliates: is the best genome defence a good offence?

Genomes, like crazy patchwork quilts, are stitched together over evolutionary time from diverse elements, including some unwelcome invaders. To deal with parasitic mobile elements, most eukaryotes employ a genome self-defensive manoeuvre to recognise and silence such elements by homology-dependent interactions with RNA-protein complexes that alter chromatin. Ciliated protozoa employ more 'offensive' tactics by actually unstitching and reassembling their somatic genomes at every sexual generation to eliminate transposons and their remnants, using as patterns the maternal genomes that were rearranged in the previous cycle. Genetic and genomic studies of the distant relatives Paramecium and Tetrahymena have begun to reveal how such events are carried out with remarkable precision. Whole genome, non-coding transcripts from the maternal genome are compared with transcripts from the zygotic genome that are processed through an RNA interference (RNAi)-related process. Sequences found only in the latter are targeted for elimination by the resulting short 'scanRNAs' in many thousand DNA splicing reactions initiated by a domesticated transposase. The involvement of widely conserved mechanisms and protein factors clearly shows the relatedness of these phenomena to RNAi-mediated heterochromatic gene silencing. Such malleability of the genome on a generational time scale also has profound evolutionary implications, possibly including the epigenetic inheritance of acquired adaptive traits.

RNA-Mediated Epigenetic Programming of Genome Rearrangements

RNA, normally thought of as a conduit in gene expression, has a novel mode of action in ciliated protozoa. Maternal RNA templates provide both an organizing guide for DNA rearrangements and a template that can transport somatic mutations to the next generation. This opportunity for RNA-mediated genome rearrangement and DNA repair is profound in the ciliate Oxytricha, which deletes 95% of its germline genome during development in a process that severely fragments its chromosomes and then sorts and reorders the hundreds of thousands of pieces remaining. Oxytricha's somatic nuclear genome is therefore an epigenome formed through RNA templates and signals arising from the previous generation. Furthermore, this mechanism of RNA-mediated epigenetic inheritance can function across multiple generations, and the discovery of maternal template RNA molecules has revealed new biological roles for RNA and has hinted at the power of RNA molecules to sculpt genomic information in cells.

Keeping the soma free of transposons: programmed DNA elimination in ciliates

Many transposon-related sequences are removed from the somatic macronucleus of ciliates during sexual reproduction. In the ciliate Tetrahymena, an RNAi-related mechanism produces small noncoding RNAs that induce heterochromatin formation, which is followed by DNA elimination. Because RNAi-related mechanisms repress transposon activities in a variety of eukaryotes, the DNA elimination mechanism of ciliates might have evolved from these types of transposon-silencing mechanisms. Nuclear dimorphism allows ciliates to identify any DNA that has invaded the germ-line micronucleus using small RNAs and a whole genome comparison of the micronucleus and the somatic macronucleus.

DNA elimination in ciliates: transposon domestication and genome surveillance

Ciliated protozoa extensively remodel their somatic genomes during nuclear development, fragmenting their chromosomes and removing large numbers of internal eliminated sequences (IESs). The sequences eliminated are unique and repetitive DNAs, including transposons. Recent studies have identified transposase proteins that appear to have been domesticated and are used by these cells to eliminate DNA not wanted in the somatic macronucleus. This DNA elimination process is guided by meiotically produced small RNAs, generated in the germline nucleus, that recognize homologous sequences leading to their removal. These scan RNAs are found in complexes with PIWI proteins. Before they search the developing genome for IESs to eliminate, they scan the parental somatic nucleus and are removed from the pool if they match homologous sequences in that previously reorganized genome. In Tetrahymena, the scan RNAs target heterochromatin modifications to mark IESs for elimination. This DNA elimination pathway in ciliates shares extensive similarity with piRNA-mediated silencing of metazoans and highlights the remarkable ability of homologous RNAs to shape developing genomes.

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.

Functional specialization of Piwi proteins in Paramecium tetraurelia from post-transcriptional gene silencing to genome remodelling

Proteins of the Argonaute family are small RNA carriers that guide regulatory complexes to their targets. The family comprises two major subclades. Members of the Ago subclade, which are present in most eukaryotic phyla, bind different classes of small RNAs and regulate gene expression at both transcriptional and post-transcriptional levels. Piwi subclade members appear to have been lost in plants and fungi and were mostly studied in metazoa, where they bind piRNAs and have essential roles in sexual reproduction. Their presence in ciliates, unicellular organisms harbouring both germline micronuclei and somatic macronuclei, offers an interesting perspective on the evolution of their functions. Here, we report phylogenetic and functional analyses of the 15 Piwi genes from Paramecium tetraurelia. We show that four constitutively expressed proteins are involved in siRNA pathways that mediate gene silencing throughout the life cycle. Two other proteins, specifically expressed during meiosis, are required for accumulation of scnRNAs during sexual reproduction and for programmed genome rearrangements during development of the somatic macronucleus. Our results indicate that Paramecium Piwi proteins have evolved to perform both vegetative and sexual functions through mechanisms ranging from post-transcriptional mRNA cleavage to epigenetic regulation of genome rearrangements.

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.

Extreme heterogeneous composition of the Paramecium caudatum macronuclear genomic DNA between hemoglobin and nucleosome assembly protein-1 genes

The intergenic region between the hemoglobin (hb) and nucleosome assembly protein-1 (nap-1) genes in the Paramecium caudatum macronuclear genome was previously found to be heterogeneously composed. Cloning of this intergenic region from the macronuclear genomic DNA identified four unique DNA fragments of different sizes. Sequencing of the cloned fragments revealed extreme heterogeneity and characteristics of both internal eliminated sequence (IES) and imprecise internal deletion sequences (IIDSs) in the intergenic region. Missing sequences were an AT-rich and direct repeats existed in their boundaries. Southern blotting of the total genomic DNA and polymerase chain reaction (PCR) of the total genomic DNAs indicated that there exist a dozen DNA fragments of different sizes in this intergenic region. It is likely that the heterogeneity found in the P. caudatum macronuclear genome results from the variable removal of an intergenic region.

Subtraction by addition: domesticated transposases in programmed DNA elimination

The ciliate Paramecium tetraurelia must eliminate approximately 60,000 short sequences from its genome to generate uninterrupted coding sequences in its somatic macronucleus. In this issue of Genes & Development, Baudry and colleagues (pp. 2478-2483) identify the protein that excises these noncoding sequences: a domesticated piggyBac transposase that has been adapted to remove what are likely the remnants of transposon insertions. This new study reveals how addition of a transposase to small RNA-directed silencing machinery can guide major genome reorganization.

PiggyMac, a domesticated piggyBac transposase involved in programmed genome rearrangements in the ciliate Paramecium tetraurelia

Programmed genome rearrangements drive functional gene assembly in ciliates during the development of the somatic macronucleus. The elimination of germline sequences is directed by noncoding RNAs and is initiated by DNA double-strand breaks, but the enzymes responsible for DNA cleavage have not been identified. We show here that PiggyMac (Pgm), a domesticated piggyBac transposase, is required for these rearrangements in Paramecium tetraurelia. A GFP-Pgm fusion localizes in developing macronuclei, where rearrangements take place, and RNAi-mediated silencing of PGM abolishes DNA cleavage. This is the first in vivo evidence suggesting an essential endonucleolytic function of a domesticated piggyBac transposase.