Programmed genome rearrangements in ciliates

Remarkable chromosomes and karyotypes: A top 10 list

Chromosomes and karyotypes are particularly rich in oddities and extremes. Described below are 10 remarkable chromosomes and karyotypes sprinkled throughout the tree of life. These include variants in chromosome number, structure, and dynamics both natural and engineered. This versatility highlights the robustness and tolerance of the mitotic and meiotic machinery to dramatic changes in chromosome and karyotype architecture. These examples also illustrate that the robustness comes at a cost, enabling the evolution of chromosomes that subvert mitosis and meiosis.

End resection and telomere healing of DNA double-strand breaks during nematode programmed DNA elimination

Most DNA double-strand breaks (DSBs) are harmful to genome integrity. However, some forms of DSBs are essential to biological processes, such as meiotic recombination and V(D)J recombination. DSBs are also required for programmed DNA elimination (PDE) in ciliates and nematodes. In nematodes, the DSBs are healed with telomere addition. While telomere addition sites have been well-characterized, little is known regarding the DSBs that fragment nematode chromosomes. Here, we used embryos from the nematode Ascaris to study the timing of PDE breaks and examine the DSBs and their end processing. Using END-seq, we characterize the DSB ends and demonstrate that DNA breaks are introduced before mitosis, followed by extensive end resection. The resection profile is unique for each break site, and the resection generates 3' overhangs before the addition of telomeres. Interestingly, telomere healing occurs much more frequently on retained DSB ends than on eliminated ends. This biased repair of the DSB ends in Ascaris may be due to the sequestration of the eliminated DNA into micronuclei, preventing their ends from telomere healing. Additional DNA breaks occur within the eliminated DNA in both Ascaris and Parascaris, ensuring chromosomal breakage and providing a fail-safe mechanism for nematode PDE.

Extensive DNA methylome rearrangement during early lamprey embryogenesis

DNA methylation (5mC) is a repressive gene regulatory mark widespread in vertebrate genomes, yet the developmental dynamics in which 5mC patterns are established vary across species. While mammals undergo two rounds of global 5mC erasure, teleosts, for example, exhibit localized maternal-to-paternal 5mC remodeling. Here, we studied 5mC dynamics during the embryonic development of sea lamprey, a jawless vertebrate which occupies a critical phylogenetic position as the sister group of the jawed vertebrates. We employed 5mC quantification in lamprey embryos and tissues, and discovered large-scale maternal-to-paternal epigenome remodeling that affects ~30% of the embryonic genome and is predominantly associated with partially methylated domains. We further demonstrate that sequences eliminated during programmed genome rearrangement (PGR), are hypermethylated in sperm prior to the onset of PGR. Our study thus unveils important insights into the evolutionary origins of vertebrate 5mC reprogramming, and how this process might participate in diverse developmental strategies.

MicroRNAs in Tetrahymena thermophila: An epigenetic regulatory mechanism in the response to cadmium stress

Among the epigenetic mechanisms based on non-coding RNA are microRNAs (miRNAs) that are involved in the post-transcriptional regulation of mRNAs. In many organisms, the expression of genes involved in the cellular response to biotic or abiotic stress depends on the regulation, generally inhibitory, performed by miRNAs. For the first time in the eukaryotic microorganism (ciliate-model) Tetrahymena thermophila, miRNAs involved in the post-transcriptional regulation of transcripts linked to the response to cadmium have been isolated and analyzed. Forty de novo miRNAs (we named tte-miRNAs) have been isolated from control and Cd-treated populations (1 or 24 h exposures). An exhaustive comparative analysis of the features of these mature tte-miRNAs and their precursor sequences (pre-tte-miRNAs) confirms that they are true miRNAs. In addition to the three types of miRNA isoforms previously described in other organisms, two new types are also described among the tte-miRNAs studied. A certain percentage of the pre-tte-miRNA sequences are in introns from genes with many introns, and have been defined as 5', 3'-tailed mirtrons. A qRT-PCR analysis of selected tte-miRNAs together with some of their targets has validated them. Cd is one of the most toxic metals for the cell, which must defend itself against its toxicity by various mechanisms, such as expulsion by membrane pumps, chelation by metallothioneins, among others. Like other toxic metals, Cd also causes a well-known series of cellular effects such as intense proteotoxicity. Many of the targets that are regulated by the tte-miRNAs are transcripts encoding proteins that fit into these defense mechanisms and toxic metal effects.

Application of RNA interference and protein localization to investigate housekeeping and developmentally regulated genes in the emerging model protozoan Paramecium caudatum

Unicellular eukaryotes represent tremendous evolutionary diversity. However, the molecular mechanisms underlying this diversity remain largely unexplored, partly due to a limitation of genetic tools to only a few model species. Paramecium caudatum is a well-known unicellular eukaryote with an unexpectedly large germline genome, of which only two percent is retained in the somatic genome following sexual processes, revealing extensive DNA elimination. However, further progress in understanding the molecular mechanisms governing this process is hampered by a lack of suitable genetic tools. Here, we report the successful application of gene knockdown and protein localization methods to interrogate the function of both housekeeping and developmentally regulated genes in P. caudatum. Using these methods, we achieved the expected phenotypes upon RNAi by feeding, and determined the localization of these proteins by microinjection of fusion constructs containing fluorescent protein or antibody tags. Lastly, we used these methods to reveal that P. caudatum PiggyMac, a domesticated piggyBac transposase, is essential for sexual development, and is likely to be an active transposase directly involved in DNA cleavage. The application of these methods lays the groundwork for future studies of gene function in P. caudatum and can be used to answer important biological questions in the future.

Single cell transcriptomics reveals UAR codon reassignment in Palmarella salina (Metopida, Armophorea) and confirms Armophorida belongs to APM clade

Anaerobes have emerged in several major lineages of ciliates, but the number of independent transitions to anaerobiosis among ciliates is unknown. The APM clade (Armophorea, Muranotrichea, Parablepharismea) represents the largest clade of obligate anaerobes among ciliates and contains free-living marine and freshwater representatives as well as gut endobionts of animals. The evolution of APM group has only recently started getting attention, and our knowledge on its phylogeny and genetics is still limited to a fraction of taxa. While ciliates portray a wide array of alternatives to the standard genetic code across numerous classes, the APM ciliates were considered to be the largest group using exclusively standard nuclear genetic code. In this study, we present a pan-ciliate phylogenomic analysis with emphasis on the APM clade, bringing the first phylogenomic analysis of the family Tropidoatractidae (Armophorea) and confirming the position of Armophorida within Armophorea. We include five newly sequenced single cell transcriptomes from marine, freshwater, and endobiotic APM ciliates - Palmarella salina, Anteclevelandella constricta, Nyctotherus sp., Caenomorpha medusula, and Thigmothrix strigosa. We report the first discovery of an alternative nuclear genetic code among APM ciliates, used by Palmarella salina (Tropidoatractidae, Armophorea), but not by its close relative, Tropidoatractus sp., and provide a comparative analysis of stop codon identity and frequency indicating the precedency to the UAG codon loss/reassignment over the UAA codon reassignment in the specific ancestor of Palmarella. Comparative genomic and proteomic studies of this group may help explain the constraints that underlie UAR stop-to-sense reassignment, the most frequent type of alternative nuclear genetic code, not only in ciliates, but eukaryotes in general.

From germline genome to highly fragmented somatic genome: genome-wide DNA rearrangement during the sexual process in ciliated protists

Genomes are incredibly dynamic within diverse eukaryotes and programmed genome rearrangements (PGR) play important roles in generating genomic diversity. However, genomes and chromosomes in metazoans are usually large in size which prevents our understanding of the origin and evolution of PGR. To expand our knowledge of genomic diversity and the evolutionary origin of complex genome rearrangements, we focus on ciliated protists (ciliates). Ciliates are single-celled eukaryotes with highly fragmented somatic chromosomes and massively scrambled germline genomes. PGR in ciliates occurs extensively by removing massive amounts of repetitive and selfish DNA elements found in the silent germline genome during development of the somatic genome. We report the partial germline genomes of two spirotrich ciliate species, namely Strombidium cf. sulcatum and Halteria grandinella, along with the most compact and highly fragmented somatic genome for S. cf. sulcatum. We provide the first insights into the genome rearrangements of these two species and compare these features with those of other ciliates. Our analyses reveal: (1) DNA sequence loss through evolution and during PGR in S. cf. sulcatum has combined to produce the most compact and efficient nanochromosomes observed to date; (2) the compact, transcriptome-like somatic genome in both species results from extensive removal of a relatively large number of shorter germline-specific DNA sequences; (3) long chromosome breakage site motifs are duplicated and retained in the somatic genome, revealing a complex model of chromosome fragmentation in spirotrichs; (4) gene scrambling and alternative processing are found throughout the core spirotrichs, offering unique opportunities to increase genetic diversity and regulation in this group.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00213-x.

The macronuclear genomic landscape within Tetrahymena thermophila

The extent of intraspecific genomic variation is key to understanding species evolutionary history, including recent adaptive shifts. Intraspecific genomic variation remains poorly explored in eukaryotic micro-organisms, especially in the nuclear dimorphic ciliates, despite their fundamental role as laboratory model systems and their ecological importance in many ecosystems. We sequenced the macronuclear genome of 22 laboratory strains of the oligohymenophoran Tetrahymena thermophila, a model species in both cellular biology and evolutionary ecology. We explored polymorphisms at the junctions of programmed eliminated sequences, and reveal their utility to barcode very closely related cells. As for other species of the genus Tetrahymena, we confirm micronuclear centromeres as gene diversification centres in T. thermophila, but also reveal a two-speed evolution in these regions. In the rest of the genome, we highlight recent diversification of genes coding for extracellular proteins and cell adhesion. We discuss all these findings in relation to this ciliate's ecology and cellular characteristics.

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.

Programmed DNA elimination in Mesorhabditis nematodes

Some species undergo programmed DNA elimination (PDE), whereby portions of the genome are systematically destroyed in somatic cells. PDE has emerged independently in several phyla, but its function is unknown. Although the mechanisms are partially solved in ciliates, PDE remains mysterious in metazoans because the study species were not yet amenable to functional approaches. We fortuitously discovered massive PDE in the free-living nematode genus Mesorhabditis, from the same family as C. elegans. As such, these species offer many experimental advantages to start elucidating the PDE mechanisms in an animal. Here, we used cytology to describe the dynamics of chromosome fragmentation and destruction in early embryos. Elimination occurs once in development, at the third embryonic cell division in the somatic blastomeres. Chromosomes are first fragmented during S phase. Next, some of the fragments fail to align on the mitotic spindle and remain outside the re-assembled nuclei after mitosis. These fragments are gradually lost after a few cell cycles. The retained fragments form new mini chromosomes, which are properly segregated in the subsequent cell divisions. With genomic approaches, we found that Mesorhabditis mainly eliminate repeated regions and also about a hundred genes. Importantly, none of the eliminated protein-coding genes are shared between closely related Mesorhabditis species. Our results strongly suggest PDE has not been selected for regulating genes with important biological functions in Mesorhabditis but rather mainly to irreversibly remove repeated sequences in the soma. We propose that PDE may target genes, provided their elimination in the soma is invisible to selection.

From oral structure to molecular evidence: new insights into the evolutionary phylogeny of the ciliate order Sessilida (Protista, Ciliophora), with the establishment of two new families and new contributions to the poorly studied family Vaginicolidae

Ciliated protists represent one of the most primitive and diverse lineages of eukaryotes, with nuclear dimorphism, a distinctive sexual process (conjugation), and extensive genome rearrangements. Among divergent ciliate lineages, the peritrich order Sessilida includes members with a colonial lifestyle, which may hint to an independent evolutionary attempt for multicellularity, although they are still single-celled organisms. To date, the evolution and phylogeny of this group are still far from clear, in part due to the paucity of molecular and/or morphological data for many taxa. In this study, we extend taxon sampling of a loricate group of sessilids by obtaining 69 new rDNA (SSU rDNA, ITS1-5.8S rDNA-ITS2, and LSU rDNA) sequences from 20 well-characterized representative species and analyze the phylogenetic relationships within Sessilida. The main findings are: (i) the genera Rhabdostyla and Campanella each represents a unique taxon at family level, supporting the establishment of two new families, i.e., Rhabdostylidae n. fam. and Campanellidae n. fam., respectively, the former being sister to a morphologically heterogeneous clade comprising Astylozoidae and several incertae sedis species and the latter occupying the basal position within the Sessilida clade; (ii) the structure of infundibular polykinety 3 is likely to be a phylogenetically informative character for resolving evolutionary relationships among sessilids; (iii) differences between sparsely and the densely arranged silverline systems could be a suprageneric taxonomic character; (iv) the monophyly of Vaginicolidae is confirmed, which is consistent with its specialized morphology, i.e., the possession of a typical peritrich lorica which might be an apomorphy for this group; (v) within Vaginicolidae, the monotypic Cothurniopsis sensu Stokes, 1893 is a synonym of Cothurnia Ehrenberg, 1831, and a new combination is created, i.e., Cothurnia valvata nov. comb.; (vi) Vaginicola sensu lato comprises at least two distinctly divergent clades, one affiliated with Thuricola and the other with a systematically puzzling clade represented by Vaginicola tincta.

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

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

Comparative analysis of single-cell genome sequencing techniques toward the characterization of germline and somatic genomes in ciliated protists

Ciliated protists contain both germline micronucleus (MIC) and somatic macronucleus (MAC) in a single cytoplasm. Programmed genome rearrangements occur in ciliates during sexual processes, and the extent of rearrangements varies dramatically among species, which lead to significant differences in genomic architectures. However, genomic sequences remain largely unknown for most ciliates due to the difficulty in culturing and in separating the germline from the somatic genome in a single cell. Single-cell whole genome amplification (WGA) has emerged as a powerful technology to characterize the genomic heterogeneity at the single-cell level. In this study, we compared two single-cell WGA, multiple displacement amplification (MDA) and multiple annealing and looping-based amplification cycles (MALBAC) in characterizing the germline and somatic genomes in ciliates with different genomic architectures. Our results showed that: 1) MALBAC exhibits strong amplification bias towards MAC genome while MDA shows bias towards MIC genome of ciliates with extensively fragmented MAC genome; 2) both MDA and MALBAC could amplify MAC genome more efficiently in ciliates with moderately fragmented MAC genome. Moreover, we found that more sample replicates could help to obtain more genomic data. Our work provides a reference for selecting the appropriate method to characterize germline and somatic genomes of ciliates.

Programmed DNA elimination in the parasitic nematode Ascaris

In most organisms, the whole genome is maintained throughout the life span. However, exceptions occur in some species where the genome is reduced during development through a process known as programmed DNA elimination (PDE). In the human and pig parasite Ascaris, PDE occurs during the 4 to 16 cell stages of embryogenesis, when germline chromosomes are fragmented and specific DNA sequences are reproducibly lost in all somatic cells. PDE was identified in Ascaris over 120 years ago, but little was known about its molecular details until recently. Genome sequencing revealed that approximately 1,000 germline-expressed genes are eliminated in Ascaris, suggesting PDE is a gene silencing mechanism. All germline chromosome ends are removed and remodeled during PDE. In addition, PDE increases the number of chromosomes in the somatic genome by splitting many germline chromosomes. Comparative genomics indicates that these germline chromosomes arose from fusion events. PDE separates these chromosomes at the fusion sites. These observations indicate that PDE plays a role in chromosome karyotype and evolution. Furthermore, comparative analysis of PDE in other parasitic and free-living nematodes illustrates conserved features of PDE, suggesting it has important biological significance. We summarize what is known about PDE in Ascaris and its relatives. We also discuss other potential functions, mechanisms, and the evolution of PDE in these parasites of humans and animals of veterinary importance.

Small RNA-mediated genome rearrangement pathways in ciliates

Most eukaryotes employ a combination of transcriptional and post-transcriptional silencing mechanisms to suppress transposons, yet ciliates employ a more extreme approach. They separate germline and somatic functions into distinct nuclei, enabling the elimination of transposons from the active somatic genome through diverse small RNA-mediated genome rearrangement pathways during sexual processes.

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.

The nematode Oscheius tipulae as a genetic model for programmed DNA elimination

Programmed DNA elimination (PDE) is a notable exception to the paradigm of genome integrity. In metazoa, PDE often occurs coincident with germline to somatic cell differentiation. During PDE, portions of genomic DNA are lost, resulting in reduced somatic genomes. Prior studies have described the sequences lost, as well as chromosome behavior, during metazoan PDE. However, a system for studying the mechanisms and consequences of PDE in metazoa is lacking. Here, we present a functional and genetic model for PDE in the free-living Rhabditidae nematode Oscheius tipulae, a family that also includes Caenorhabditis elegans. O. tipulae was recently suggested to eliminate DNA. Using staged embryos and DNA FISH, we showed that O. tipulae PDE occurs during embryogenesis at the 8-16 cell stages. We identified a conserved motif, named Sequence For Elimination (SFE), for all 12 break sites on the six chromosomes at the junctions of retained and eliminated DNA. SFE mutants exhibited a "fail-to-eliminate" phenotype only at the modified sites. END-seq revealed that breaks can occur at multiple positions within the SFE, with extensive end resection followed by telomere addition to both retained and eliminated ends. We identified many functional SFEs at the chromosome ends through END-seq in the wild-type embryos, genome sequencing of SFE mutants, and comparative genomics of 23 wild isolates. We suggest that these alternative SFEs provide flexibility in the sequences eliminated and a fail-safe mechanism for PDE. These studies establish O. tipulae as a new, attractive model for studying the mechanisms and consequences of PDE in a metazoan.

Individual Genetic Heterogeneity

Genetic variation has been widely covered in literature, however, not from the perspective of an individual in any species. Here, a synthesis of genetic concepts and variations relevant for individual genetic constitution is provided. All the different levels of genetic information and variation are covered, ranging from whether an organism is unmixed or hybrid, has variations in genome, chromosomes, and more locally in DNA regions, to epigenetic variants or alterations in selfish genetic elements. Genetic constitution and heterogeneity of microbiota are highly relevant for health and wellbeing of an individual. Mutation rates vary widely for variation types, e.g., due to the sequence context. Genetic information guides numerous aspects in organisms. Types of inheritance, whether Mendelian or non-Mendelian, zygosity, sexual reproduction, and sex determination are covered. Functions of DNA and functional effects of variations are introduced, along with mechanism that reduce and modulate functional effects, including TARAR countermeasures and intraindividual genetic conflict. TARAR countermeasures for tolerance, avoidance, repair, attenuation, and resistance are essential for life, integrity of genetic information, and gene expression. The genetic composition, effects of variations, and their expression are considered also in diseases and personalized medicine. The text synthesizes knowledge and insight on individual genetic heterogeneity and organizes and systematizes the central concepts.

Identification of novel, functional, long noncoding RNAs involved in programmed, large-scale genome rearrangements

Noncoding RNAs (ncRNAs) make up to ∼98% percent of the transcriptome of a given organism. In recent years, one relatively new class of ncRNAs, long noncoding RNAs (lncRNAs), were shown to be more than mere by-products of gene expression and regulation. The unicellular eukaryote Paramecium tetraurelia is a member of the ciliate phylum, an extremely heterogeneous group of organisms found in most bodies of water across the globe. A hallmark of ciliate genetics is nuclear dimorphism and programmed elimination of transposons and transposon-derived DNA elements, the latter of which is essential for the maintenance of the somatic genome. Paramecium and ciliates in general harbor a plethora of different ncRNA species, some of which drive the process of large-scale genome rearrangements, including DNA elimination, during sexual development. Here, we identify and validate the first known functional lncRNAs in ciliates to date. Using deep-sequencing and subsequent bioinformatic processing and experimental validation, we show that Paramecium expresses at least 15 lncRNAs. These candidates were predicted by a highly conservative pipeline, and informatic analyses hint at differential expression during development. Depletion of two lncRNAs, lnc1 and lnc15, resulted in clear phenotypes, decreased survival, morphological impairment, and a global effect on DNA elimination.

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.

Ciliate Research. From Myth to Trendsetting Science

This special issue of the Journal of Eukaryotic Microbiology (JEM) summarizes achievements obtained by generations of researchers with ciliates in widely different disciplines. In fact, ciliates range among the first cells seen under the microscope centuries ago. Their beauty made them an object of scientia amabilis and their manifold reactions made them attractive for college experiments and finally challenged causal analyses at the cellular level. Some of this work was honored by a Nobel Prize. Some observations yielded a baseline for additional novel discoveries, occasionally facilitated by specific properties of some ciliates. This also offers some advantage in the exploration of closely related parasites (malaria). Articles contributed here by colleagues from all over the world encompass a broad spectrum of ciliate life, from genetics to evolution, from molecular cell biology to ecology, from intercellular signaling to epigenetics etc. This introductory chapter, largely based on my personal perception, aims at integrating work presented in this special issue of JEM into a broader historical context up to current research.

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

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

Macronuclear development in ciliates, with a focus on nuclear architecture

Ciliates are defined by the presence of dimorphic nuclei as they have both a somatic macronucleus and germline micronucleus within each individual cell. The size and structure of both germline micronuclei and somatic macronuclei varies tremendously among ciliates. Except just after conjugation (i.e. the nuclear exchange in sexual cycle), the germline micronucleus is transcriptionally-inactive and contains canonical chromosomes that will be inherited between generations. In contrast, the transcriptionally-active macronucleus contains chromosomes that vary in size in different classes of ciliates, with some lineages having extensively-fragmented gene-sized somatic chromosomes while others contain longer multigene chromosomes. Here, we describe the variation in somatic macronuclear architecture in lineages sampled across the ciliate tree of life, specifically focusing on lineages with extensively fragmented chromosomes (e.g. the classes Phyllopharyngea and Spirotrichea). Further, we synthesize information from the literature on the development of ciliate macronuclei, focusing on changes in nuclear architecture throughout life cycles. These data highlight the tremendous diversity among ciliate nuclear cycles, extend our understanding of patterns of genome evolution, and provide insight into different germline and somatic nuclear features (e.g. nuclear structure and development) among eukaryotes.

Epigenetic influences of mobile genetic elements on ciliate genome architecture and evolution

Mobile genetic elements (MGEs) are transient genetic material that can move either within a single organism's genome or between individuals or species. While historically considered 'junk' DNA (i.e. deleterious or at best neutral), more recent studies reveal the adaptive advantages MGEs provide in lineages across the tree of life. Ciliates, a group of single-celled microbial eukaryotes characterized by nuclear dimorphism, exemplify how epigenetic influences from MGEs shape genome architecture and patterns of molecular evolution. Ciliate nuclear dimorphism may have evolved as a response to transposon invasion and ciliates have since co-opted transposons to carry out programmed DNA deletion. Another example of the effect of MGEs is in providing mechanisms for lateral gene transfer from bacteria, which introduces genetic diversity and, in several cases, drives ecological specialization in ciliates. As a third example, the integration of viral DNA, likely through transduction, provides new genetic material and can change the way host cells defend themselves against other viral pathogens. We argue that the acquisition of MGEs through non-Mendelian patterns of inheritance, coupled with their effects on ciliate genome architecture and expression and persistence throughout evolutionary history, exemplify how the transmission of mobile elements should be considered a mechanism of transgenerational epigenetic inheritance.

Polycomb Repressive Complex 2 in Eukaryotes-An Evolutionary Perspective

Polycomb repressive complex 2 (PRC2) represents a group of evolutionarily conserved multi-subunit complexes that repress gene transcription by introducing trimethylation of lysine 27 on histone 3 (H3K27me3). PRC2 activity is of key importance for cell identity specification and developmental phase transitions in animals and plants. The composition, biochemistry, and developmental function of PRC2 in animal and flowering plant model species are relatively well described. Recent evidence demonstrates the presence of PRC2 complexes in various eukaryotic supergroups, suggesting conservation of the complex and its function. Here, we provide an overview of the current understanding of PRC2-mediated repression in different representatives of eukaryotic supergroups with a focus on the green lineage. By comparison of PRC2 in different eukaryotes, we highlight the possible common and diverged features suggesting evolutionary implications and outline emerging questions and directions for future research of polycomb repression and its evolution.

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.

Phenotypic plasticity through disposable genetic adaptation in ciliates

Ciliates have an extraordinary genetic system in which each cell harbors two distinct kinds of nucleus, a transcriptionally active somatic nucleus and a quiescent germline nucleus. The latter undergoes classical, heritable genetic adaptation, while adaptation of the somatic nucleus is only short-term and thus disposable. The ecological and evolutionary relevance of this nuclear dimorphism have never been well formalized, which is surprising given the long history of using ciliates such as Tetrahymena and Paramecium as model organisms. We present a novel, alternative explanation for ciliate nuclear dimorphism which, we argue, should be considered an instrument of phenotypic plasticity by somatic selection on the level of the ciliate clone, as if it were a diffuse multicellular organism. This viewpoint helps to put some enigmatic aspects of ciliate biology into perspective and presents the diversity of ciliates as a large natural experiment that we can exploit to study phenotypic plasticity and organismality.

A telomere-to-telomere assembly of Oscheius tipulae and the evolution of rhabditid nematode chromosomes

Eukaryotic chromosomes have phylogenetic persistence. In many taxa, each chromosome has a single functional centromere with essential roles in spindle attachment and segregation. Fusion and fission can generate chromosomes with no or multiple centromeres, leading to genome instability. Groups with holocentric chromosomes (where centromeric function is distributed along each chromosome) might be expected to show karyotypic instability. This is generally not the case, and in Caenorhabditis elegans, it has been proposed that the role of maintenance of a stable karyotype has been transferred to the meiotic pairing centers, which are found at one end of each chromosome. Here, we explore the phylogenetic stability of nematode chromosomes using a new telomere-to-telomere assembly of the rhabditine nematode Oscheius tipulae generated from nanopore long reads. The 60-Mb O. tipulae genome is resolved into six chromosomal molecules. We find the evidence of specific chromatin diminution at all telomeres. Comparing this chromosomal O. tipulae assembly with chromosomal assemblies of diverse rhabditid nematodes, we identify seven ancestral chromosomal elements (Nigon elements) and present a model for the evolution of nematode chromosomes through rearrangement and fusion of these elements. We identify frequent fusion events involving NigonX, the element associated with the rhabditid X chromosome, and thus sex chromosome-associated gene sets differ markedly between species. Despite the karyotypic stability, gene order within chromosomes defined by Nigon elements is not conserved. Our model for nematode chromosome evolution provides a platform for investigation of the tensions between local genome rearrangement and karyotypic evolution in generating extant genome architectures.

Disruption of a ∼23-24 nucleotide small RNA pathway elevates DNA damage responses in Tetrahymena thermophila

Endogenous RNA interference (RNAi) pathways regulate a wide range of cellular processes in diverse eukaryotes, yet in the ciliated eukaryote, Tetrahymena thermophila, the cellular purpose of RNAi pathways that generate ∼23–24 nucleotide (nt) small (s)RNAs has remained unknown. Here, we investigated the phenotypic and gene expression impacts on vegetatively growing cells when genes involved in ∼23–24 nt sRNA biogenesis are disrupted. We observed slower proliferation and increased expression of genes involved in DNA metabolism and chromosome organization and maintenance in sRNA biogenesis mutants RSP1Δ, RDN2Δ, and RDF2Δ. In addition, RSP1Δ and RDN2Δ cells frequently exhibited enlarged chromatin extrusion bodies, which are nonnuclear, DNA-containing structures that may be akin to mammalian micronuclei. Expression of homologous recombination factor Rad51 was specifically elevated in RSP1Δ and RDN2Δ strains, with Rad51 and double-stranded DNA break marker γ-H2A.X localized to discrete macronuclear foci. In addition, an increase in Rad51 and γ-H2A.X foci was also found in knockouts of TWI8, a macronucleus-localized PIWI protein. Together, our findings suggest that an evolutionarily conserved role for RNAi pathways in maintaining genome integrity may be extended even to the early branching eukaryotic lineage that gave rise to Tetrahymena thermophila.

Boon and Bane of DNA Double-Strand Breaks

DNA double-strand breaks (DSBs), interrupting the genetic information, are elicited by various environmental and endogenous factors. They bear the risk of cell lethality and, if mis-repaired, of deleterious mutation. This negative impact is contrasted by several evolutionary achievements for DSB processing that help maintaining stable inheritance (correct repair, meiotic cross-over) and even drive adaptation (immunoglobulin gene recombination), differentiation (chromatin elimination) and speciation by creating new genetic diversity via DSB mis-repair. Targeted DSBs play a role in genome editing for research, breeding and therapy purposes. Here, I survey possible causes, biological effects and evolutionary consequences of DSBs, mainly for students and outsiders.

An Optimized and Versatile Counter-Flow Centrifugal Elutriation Workflow to Obtain Synchronized Eukaryotic Cells

Cell synchronization is a powerful tool to understand cell cycle events and its regulatory mechanisms. Counter-flow centrifugal elutriation (CCE) is a more generally desirable method to synchronize cells because it does not significantly alter cell behavior and/or cell cycle progression, however, adjusting specific parameters in a cell type/equipment-dependent manner can be challenging. In this paper, we used the unicellular eukaryotic model organism, Tetrahymena thermophila as a testing system for optimizing CCE workflow. Firstly, flow cytometry conditions were identified that reduced nuclei adhesion and improved the assessment of cell cycle stage. We then systematically examined how to achieve the optimal conditions for three critical factors affecting the outcome of CCE, including loading flow rate, collection flow rate and collection volume. Using our optimized workflow, we obtained a large population of highly synchronous G1-phase Tetrahymena as measured by 5-ethynyl-2'-deoxyuridine (EdU) incorporation into nascent DNA strands, bulk DNA content changes by flow cytometry, and cell cycle progression by light microscopy. This detailed protocol can be easily adapted to synchronize other eukaryotic cells.

Constructive Neutral Evolution 20 Years Later

Evolution has led to a great diversity that ranges from elegant simplicity to ornate complexity. Many complex features are often assumed to be more functional or adaptive than their simpler alternatives. However, in 1999, Arlin Stolzfus published a paper in the Journal of Molecular Evolution that outlined a framework in which complexity can arise through a series of non-adaptive steps. He called this framework Constructive Neutral Evolution (CNE). Despite its two-decade-old roots, many evolutionary biologists still appear to be unaware of this explanatory framework for the origins of complexity. In this perspective piece, we explain the theory of CNE and how it changes the order of events in narratives that describe the evolution of complexity. We also provide an extensive list of cellular features that may have become more complex through CNE. We end by discussing strategies to determine whether complexity arose through neutral or adaptive processes.

Evolutionary Plasticity of Mating-Type Determination Mechanisms in Paramecium aurelia Sibling Species

The Paramecium aurelia complex, a group of morphologically similar but sexually incompatible sibling species, is a unique example of the evolutionary plasticity of mating-type systems. Each species has two mating types, O (Odd) and E (Even). Although O and E types are homologous in all species, three different modes of determination and inheritance have been described: genetic determination by Mendelian alleles, stochastic developmental determination, and maternally inherited developmental determination. Previous work in three species of the latter kind has revealed the key roles of the E-specific transmembrane protein mtA and its highly specific transcription factor mtB: type O clones are produced by maternally inherited genome rearrangements that inactivate either mtA or mtB during development. Here we show, through transcriptome analyses in five additional species representing the three determination systems, that mtA expression specifies type E in all cases. We further show that the Mendelian system depends on functional and nonfunctional mtA alleles, and identify novel developmental rearrangements in mtA and mtB which now explain all cases of maternally inherited mating-type determination. Epistasis between these genes likely evolved from less specific interactions between paralogs in the P. aurelia common ancestor, after a whole-genome duplication, but the mtB gene was subsequently lost in three P. aurelia species which appear to have returned to an ancestral regulation mechanism. These results suggest a model accounting for evolutionary transitions between determination systems, and highlight the diversity of molecular solutions explored among sibling species to maintain an essential mating-type polymorphism in cell populations.

The Role of Small Noncoding RNA in DNA Double-Strand Break Repair

DNA damage is a common phenomenon promoted through a variety of exogenous and endogenous factors. The DNA damage response (DDR) pathway involves a wide range of proteins, and as was indicated, small noncoding RNAs (sncRNAs). These are double-strand break-induced RNAs (diRNAs) and DNA damage response small RNA (DDRNA). Moreover, RNA binding proteins (RBPs) and RNA modifications have also been identified to modulate diRNA and DDRNA function in the DDR process. Several theories have been formulated regarding the synthesis and function of these sncRNAs during DNA repair; nevertheless, these pathways’ molecular details remain unclear. Here, we review the current knowledge regarding the mechanisms of diRNA and DDRNA biosynthesis and discuss the role of sncRNAs in maintaining genome stability.

Comprehensive Chromosome End Remodeling during Programmed DNA Elimination

Germline and somatic genomes are in general the same in a multicellular organism. However, programmed DNA elimination leads to a reduced somatic genome compared to germline cells. Previous work on the parasitic nematode Ascaris demonstrated that programmed DNA elimination encompasses high-fidelity chromosomal breaks and loss of specific genome sequences including a major tandem repeat of 120 bp and ~1,000 germline-expressed genes. However, the precise chromosomal locations of these repeats, breaks regions, and eliminated genes remained unknown. We used PacBio long-read sequencing and chromosome conformation capture (Hi-C) to obtain fully assembled chromosomes of Ascaris germline and somatic genomes, enabling a complete chromosomal view of DNA elimination. We found that all 24 germline chromosomes undergo comprehensive chromosome end remodeling with DNA breaks in their subtelomeric regions and loss of distal sequences including the telomeres at both chromosome ends. All new Ascaris somatic chromosome ends are recapped by de novo telomere healing. We provide an ultrastructural analysis of Ascaris DNA elimination and show that eliminated DNA is incorporated into double membrane-bound structures, similar to micronuclei, during telophase of a DNA elimination mitosis. These micronuclei undergo dynamic changes including loss of active histone marks and localize to the cytoplasm following daughter nuclei formation and cytokinesis where they form autophagosomes. Comparative analysis of nematode chromosomes suggests that chromosome fusions occurred, forming Ascaris sex chromosomes that become independent chromosomes following DNA elimination breaks in somatic cells. These studies provide the first chromosomal view and define novel features and functions of metazoan programmed DNA elimination.