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

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

MITE infestation accommodated by genome editing in the germline genome of the ciliate Blepharisma.

During their development following sexual conjugation, ciliates excise numerous internal eliminated sequences (IESs) from a copy of the germline genome to produce the functional somatic genome. Most IESs are thought to have originated from transposons, but the presumed homology is often obscured by sequence decay. To obtain more representative perspectives on the nature of IESs and ciliate genome editing, we assembled 40,000 IESs of Blepharisma stoltei, a species belonging to a lineage (Heterotrichea) that diverged early from those of the intensively studied model ciliate species. About a quarter of IESs were short (<115 bp), largely nonrepetitive, and with a pronounced ~10 bp periodicity in length; the remainder were longer (up to 7 kbp) and nonperiodic and contained abundant interspersed repeats. Contrary to the expectation from current models, the assembled Blepharisma germline genome encodes few transposases. Instead, its most abundant repeat (8,000 copies) is a Miniature Inverted-repeat Transposable Element (MITE), apparently a deletion derivative of a germline-limited Pogo-family transposon. We hypothesize that MITEs are an important source of IESs whose proliferation is eventually self-limiting and that rather than defending the germline genomes against mobile elements, transposase domestication actually facilitates the accumulation of junk DNA.

Chromatin remodeling is required for sRNA-guided DNA elimination in Paramecium.

Small RNAs mediate the silencing of transposable elements and other genomic loci, increasing nucleosome density and preventing undesirable gene expression. The unicellular ciliate Paramecium is a model to study dynamic genome organization in eukaryotic cells, given its unique feature of nuclear dimorphism. Here, the formation of the somatic macronucleus during sexual reproduction requires eliminating thousands of transposon remnants (IESs) and transposable elements scattered throughout the germline micronuclear genome. The elimination process is guided by Piwi-associated small RNAs and leads to precise cleavage at IES boundaries. Here we show that IES recognition and precise excision are facilitated by recruiting ISWI1, a Paramecium homolog of the chromatin remodeler ISWI. ISWI1 knockdown substantially inhibits DNA elimination, quantitatively similar to development-specific sRNA gene knockdowns but with much greater aberrant IES excision at alternative boundaries. We also identify key development-specific sRNA biogenesis and transport proteins, Ptiwi01 and Ptiwi09, as ISWI1 cofactors in our co-immunoprecipitation studies. Nucleosome profiling indicates that increased nucleosome density correlates with the requirement for ISWI1 and other proteins necessary for IES excision. We propose that chromatin remodeling together with small RNAs is essential for efficient and precise DNA elimination in Paramecium.

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.

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

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

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

SUMMARY: Ciliates are single-celled eukaryotes that eliminate specific, interspersed DNA sequences (internally eliminated sequences, IESs) from their genomes during development. These are challenging to annotate and assemble because IES-containing sequences are typically much less abundant in the cell than those without, and IES sequences themselves often contain repetitive and low-complexity sequences. Long-read sequencing technologies from Pacific Biosciences and Oxford Nanopore have the potential to reconstruct longer IESs than has been possible with short reads but require a different assembly strategy. Here we present BleTIES, a software toolkit for detecting, assembling, and analyzing IESs using mapped long reads. AVAILABILITY AND IMPLEMENTATION: BleTIES is implemented in Python 3. Source code is available at https://github.com/Swart-lab/bleties (MIT license) and also distributed via Bioconda. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The macronuclear genome of the Antarctic psychrophilic marine ciliate Euplotes focardii reveals new insights on molecular cold adaptation.

The macronuclear (MAC) genomes of ciliates belonging to the genus Euplotes species are comprised of numerous small DNA molecules, nanochromosomes, each typically encoding a single gene. These genomes are responsible for all gene expression during vegetative cell growth. Here, we report the analysis of the MAC genome from the Antarctic psychrophile Euplotes focardii. Nanochromosomes containing bacterial sequences were not found, suggesting that phenomena of horizontal gene transfer did not occur recently, even though this ciliate species has a substantial associated bacterial consortium. As in other euplotid species, E. focardii MAC genes are characterized by a high frequency of translational frameshifting. Furthermore, in order to characterize differences that may be consequent to cold adaptation and defense to oxidative stress, the main constraints of the Antarctic marine microorganisms, we compared E. focardii MAC genome with those available from mesophilic Euplotes species. We focussed mainly on the comparison of tubulin, antioxidant enzymes and heat shock protein (HSP) 70 families, molecules which possess peculiar characteristic correlated with cold adaptation in E. focardii. We found that α-tubulin genes and those encoding SODs and CATs antioxidant enzymes are more numerous than in the mesophilic Euplotes species. Furthermore, the phylogenetic trees showed that these molecules are divergent in the Antarctic species. In contrast, there are fewer hsp70 genes in E. focardii compared to mesophilic Euplotes and these genes do not respond to thermal stress but only to oxidative stress. Our results suggest that molecular adaptation to cold and oxidative stress in the Antarctic environment may not only be due to particular amino acid substitutions but also due to duplication and divergence of paralogous genes.

Two Sets of Piwi Proteins Are Involved in Distinct sRNA Pathways Leading to Elimination of Germline-Specific DNA.

Piwi proteins and piRNAs protect eukaryotic germlines against the spread of transposons. During development in the ciliate Paramecium, two Piwi-dependent sRNA classes are involved in the elimination of transposons and transposon-derived DNA: scan RNAs (scnRNAs), associated with Ptiwi01 and Ptiwi09, and iesRNAs, whose binding partners we now identify as Ptiwi10 and Ptiwi11. scnRNAs derive from the maternal genome and initiate DNA elimination during development, whereas iesRNAs continue DNA targeting until the removal process is complete. Here, we show that scnRNAs and iesRNAs are processed by distinct Dicer-like proteins and bind Piwi proteins in a mutually exclusive manner, suggesting separate biogenesis pathways. We also demonstrate that the PTIWI10 gene is transcribed from the developing nucleus and that its transcription depends on prior DNA excision, suggesting a mechanism of gene expression control triggered by the removal of short DNA segments interrupting the gene.

The Macronuclear Genome of Stentor coeruleus Reveals Tiny Introns in a Giant Cell.

The giant, single-celled organism Stentor coeruleus has a long history as a model system for studying pattern formation and regeneration in single cells. Stentor [1, 2] is a heterotrichous ciliate distantly related to familiar ciliate models, such as Tetrahymena or Paramecium. The primary distinguishing feature of Stentor is its incredible size: a single cell is 1 mm long. Early developmental biologists, including T.H. Morgan [3], were attracted to the system because of its regenerative abilities-if large portions of a cell are surgically removed, the remnant reorganizes into a normal-looking but smaller cell with correct proportionality [2, 3]. These biologists were also drawn to Stentor because it exhibits a rich repertoire of behaviors, including light avoidance, mechanosensitive contraction, food selection, and even the ability to habituate to touch, a simple form of learning usually seen in higher organisms [4]. While early microsurgical approaches demonstrated a startling array of regenerative and morphogenetic processes in this single-celled organism, Stentor was never developed as a molecular model system. We report the sequencing of the Stentor coeruleus macronuclear genome and reveal key features of the genome. First, we find that Stentor uses the standard genetic code, suggesting that ciliate-specific genetic codes arose after Stentor branched from other ciliates. We also discover that ploidy correlates with Stentor's cell size. Finally, in the Stentor genome, we discover the smallest spliceosomal introns reported for any species. The sequenced genome opens the door to molecular analysis of single-cell regeneration in Stentor.

The eukaryotic way to defend and edit genomes by sRNA-targeted DNA deletion.

While there is currently burgeoning interest in the application of the CRISPR/Cas (clustered regularly interspaced short palindromic repeats/CRISPR-associated genes) to genome editing, it is perhaps not widely appreciated that this is the second discovery of a small RNA (sRNA)-targeted DNA-deletion system. The first sRNA-targeted DNA-deletion system to be discovered, which we call IES/Ias (internal eliminated sequence/IES-associated genes) to contrast with CRISPR/Cas, is found in ciliates, and, like CRISPR/Cas, is thought to serve as a form of immune defense against invasive DNAs. The manner in which the ciliate IES/Ias system functions is distinct from that of the CRISPR/Cas system in archaea and bacteria, and arose independently through a synthesis of RNA interference-derived and DNA-specific molecular components. Despite the major differences between CRISPR/Cas and IES/Ias, both systems face similar conceptual challenges in targeting invasive DNAs. In this review, we focus on the discovery, effects, function, and evolutionary consequences of the IES/Ias system.

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.

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.

The architecture of a scrambled genome reveals massive levels of genomic rearrangement during development.

Programmed DNA rearrangements in the single-celled eukaryote Oxytricha trifallax completely rewire its germline into a somatic nucleus during development. This elaborate, RNA-mediated pathway eliminates noncoding DNA sequences that interrupt gene loci and reorganizes the remaining fragments by inversions and permutations to produce functional genes. Here, we report the Oxytricha germline genome and compare it to the somatic genome to present a global view of its massive scale of genome rearrangements. The remarkably encrypted genome architecture contains >3,500 scrambled genes, as well as >800 predicted germline-limited genes expressed, and some posttranslationally modified, during genome rearrangements. Gene segments for different somatic loci often interweave with each other. Single gene segments can contribute to multiple, distinct somatic loci. Terminal precursor segments from neighboring somatic loci map extremely close to each other, often overlapping. This genome assembly provides a draft of a scrambled genome and a powerful model for studies of genome rearrangement.

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.

The draft assembly of the radically organized Stylonychia lemnae macronuclear genome.

Stylonychia lemnae is a classical model single-celled eukaryote, and a quintessential ciliate typified by dimorphic nuclei: A small, germline micronucleus and a massive, vegetative macronucleus. The genome within Stylonychia's macronucleus has a very unusual architecture, comprised variably and highly amplified "nanochromosomes," each usually encoding a single gene with a minimal amount of surrounding noncoding DNA. As only a tiny fraction of the Stylonychia genes has been sequenced, and to promote research using this organism, we sequenced its macronuclear genome. We report the analysis of the 50.2-Mb draft S. lemnae macronuclear genome assembly, containing in excess of 16,000 complete nanochromosomes, assembled as less than 20,000 contigs. We found considerable conservation of fundamental genomic properties between S. lemnae and its close relative, Oxytricha trifallax, including nanochromosomal gene synteny, alternative fragmentation, and copy number. Protein domain searches in Stylonychia revealed two new telomere-binding protein homologs and the presence of linker histones. Among the diverse histone variants of S. lemnae and O. trifallax, we found divergent, coexpressed variants corresponding to four of the five core nucleosomal proteins (H1.2, H2A.6, H2B.4, and H3.7) suggesting that these ciliates may possess specialized nucleosomes involved in genome processing during nuclear differentiation. The assembly of the S. lemnae macronuclear genome demonstrates that largely complete, well-assembled highly fragmented genomes of similar size and complexity may be produced from one library and lane of Illumina HiSeq 2000 shotgun sequencing. The provision of the S. lemnae macronuclear genome sets the stage for future detailed experimental studies of chromatin-mediated, RNA-guided developmental genome rearrangements.

Functional diversification of Dicer-like proteins and small RNAs required for genome sculpting.

In eukaryotes, small RNAs (sRNAs) have key roles in development, gene expression regulation, and genome integrity maintenance. In ciliates, such as Paramecium, sRNAs form the heart of an epigenetic system that has evolved from core eukaryotic gene silencing components to selectively target DNA for deletion. In Paramecium, somatic genome development from the germline genome accurately eliminates the bulk of typically gene-interrupting, noncoding DNA. We have discovered an sRNA class (internal eliminated sequence [IES] sRNAs [iesRNAs]), arising later during Paramecium development, which originates from and precisely delineates germline DNA (IESs) and complements the initial sRNAs ("scan" RNAs [scnRNAs]) in targeting DNA for elimination. We show that whole-genome duplications have facilitated successive differentiations of Paramecium Dicer-like proteins, leading to cooperation between Dcl2 and Dcl3 to produce scnRNAs and to the production of iesRNAs by Dcl5. These innovations highlight the ability of sRNA systems to acquire capabilities, including those in genome development and integrity.

The Oxytricha trifallax macronuclear genome: a complex eukaryotic genome with 16,000 tiny 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.

Copy number variations of 11 macronuclear chromosomes and their gene expression in Oxytricha trifallax.

Ciliated protozoa are peculiar for their nuclear dimorphism, wherein two types of nuclei divide nuclear functions: a germline micronucleus (MIC) is transcriptionally inert during vegetative growth, but serves as the genetic blueprint for the somatic macronucleus (MAC), which is responsible for all transcripts supporting cell growth and reproduction. While all the advantages/disadvantages associated with nuclear dimorphism are not clear, an essential advantage seems to be the ability to produce a highly polyploid MAC, which then allows for the maintenance of extremely large single cells - many ciliate cells are larger than small metazoa. In some ciliate classes, chromosomes in the MAC are extensively fragmented to create extremely short chromosomes that often carry single genes, and these chromosomes may be present in different copy numbers, resulting in different ploidies. While using gene copy number to regulate gene expression is limited in most eukaryotic systems, the extensive fragmentation in some ciliate classes provides this opportunity to every MAC gene. However, it is still unclear if this mechanism is in fact used extensively in these ciliates. To address this, we have quantified copy numbers of 11 MAC chromosomes and their gene expression in Oxytricha trifallax (CI: Spirotrichea). We compared copy numbers between two subpopulations of O. trifallax, and copy numbers of 7 orthologous genes between O. trifallax and the closely related Stylonychia lemnae. We show that copy numbers of MAC chromosomes are variable, dynamic, and positively correlated to gene expression. These features might be conserved in all spirotrichs, and might exist in other classes of ciliates with heavily fragmented MAC chromosomes.

Characterization and taxonomic validity of the ciliate Oxytricha trifallax (Class Spirotrichea) based on multiple gene sequences: limitations in identifying genera solely by morphology.

Oxytricha trifallax - an established model organism for studying genome rearrangements, chromosome structure, scrambled genes, RNA-mediated epigenetic inheritance, and other phenomena - has been the subject of a nomenclature controversy for several years. Originally isolated as a sibling species of O. fallax, O. trifallax was reclassified in 1999 as Sterkiella histriomuscorum, a previously identified species, based on morphological similarity. The proper identification of O. trifallax is crucial to resolve in order to prevent confusion in both the comparative genomics and the general scientific communities. We analyzed nine conserved nuclear gene sequences between the two given species and several related ciliates. Phylogenetic analyses suggest that O. trifallax and a bona fide S. histriomuscorum have accumulated significant evolutionary divergence from each other relative to other ciliates such that they should be unequivocally classified as separate species. We also describe the original isolation of O. trifallax, including its comparison to O. fallax, and we provide criteria to identify future isolates of O. trifallax.

The Oxytricha trifallax mitochondrial genome.

The Oxytricha trifallax mitochondrial genome contains the largest sequenced ciliate mitochondrial chromosome (~70 kb) plus a ~5-kb linear plasmid bearing mitochondrial telomeres. We identify two new ciliate split genes (rps3 and nad2) as well as four new mitochondrial genes (ribosomal small subunit protein genes: rps- 2, 7, 8, 10), previously undetected in ciliates due to their extreme divergence. The increased size of the Oxytricha mitochondrial genome relative to other ciliates is primarily a consequence of terminal expansions, rather than the retention of ancestral mitochondrial genes. Successive segmental duplications, visible in one of the two Oxytricha mitochondrial subterminal regions, appear to have contributed to the genome expansion. Consistent with pseudogene formation and decay, the subtermini possess shorter, more loosely packed open reading frames than the remainder of the genome. The mitochondrial plasmid shares a 251-bp region with 82% identity to the mitochondrial chromosome, suggesting that it most likely integrated into the chromosome at least once. This region on the chromosome is also close to the end of the most terminal member of a series of duplications, hinting at a possible association between the plasmid and the duplications. The presence of mitochondrial telomeres on the mitochondrial plasmid suggests that such plasmids may be a vehicle for lateral transfer of telomeric sequences between mitochondrial genomes. We conjecture that the extreme divergence observed in ciliate mitochondrial genomes may be due, in part, to repeated invasions by relatively error-prone DNA polymerase-bearing mobile elements.