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

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.

The Evolutionary Transition of the RNA World to Obcells to Cellular-Based Life

The obcell hypothesis is a proposed route for the RNA world to develop into a primitive cellular one. It posits that this transition began with the emergence of the proto-ribosome which enabled RNA to colonise the external surface of lipids by the synthesis of amphipathic peptidyl-RNAs. The obcell hypothesis also posits that the emergence of a predation-based ecosystem provided a selection mechanism for continued sophistication amongst early life forms. Here, I argue for this hypothesis owing to its significant explanatory power; it offers a rationale why a ribosome which initially was capable only of producing short non-coded peptides was advantageous and it forgoes issues related to maintaining a replicating RNA inside a lipid enclosure. I develop this model by proposing that the evolutionary selection for improved membrane anchors resulted in the emergence of primitive membrane pores which enabled obcells to gradually evolve into a cellular morphology. Moreover, I introduce a model of obcell production which advances that tRNAs developed from primers of the RNA world.

Decryption of the survival "black box": gene family expansion promotes the encystment in ciliated protists

BACKGROUND

Encystment is an important survival strategy extensively employed by microbial organisms to survive unfavorable conditions. Single-celled ciliated protists (ciliates) are popular model eukaryotes for studying encystment, whereby these cells degenerate their ciliary structures and develop cyst walls, then reverse the process under more favorable conditions. However, to date, the evolutionary basis and mechanism for encystment in ciliates is largely unknown. With the rapid development of high-throughput sequencing technologies, genome sequencing and comparative genomics of ciliates have become effective methods to provide insights into above questions.

RESULTS

Here, we profiled the MAC genome of Pseudourostyla cristata, a model hypotrich ciliate for encystment studies. Like other hypotrich MAC genomes, the P. cristata MAC genome is extremely fragmented with a single gene on most chromosomes, and encodes introns that are generally small and lack a conserved branch point for pre-mRNA splicing. Gene family expansion analyses indicate that multiple gene families involved in the encystment are expanded during the evolution of P. cristata. Furthermore, genomic comparisons with other five representative hypotrichs indicate that gene families of phosphorelay sensor kinase, which play a role in the two-component signal transduction system that is related to encystment, show significant expansion among all six hypotrichs. Additionally, cyst wall-related chitin synthase genes have experienced structural changes that increase them from single-exon to multi-exon genes during evolution. These genomic features potentially promote the encystment in hypotrichs and enhance their ability to survive in adverse environments during evolution.

CONCLUSIONS

We systematically investigated the genomic structure of hypotrichs and key evolutionary phenomenon, gene family expansion, for encystment promotion in ciliates. In summary, our results provided insights into the evolutionary mechanism of encystment in ciliates.

Global and local functions of the Fused kinase ortholog CdaH in intracellular patterning in Tetrahymena

Ciliates assemble numerous microtubular structures into complex cortical patterns. During ciliate division, the pattern is duplicated by intracellular segmentation that produces a tandem of daughter cells. In Tetrahymena thermophila, the induction and positioning of the division boundary involves two mutually antagonistic factors: posterior CdaA (cyclin E) and anterior CdaI (Hippo kinase). Here, we characterized the related cdaH-1 allele, which confers a pleiotropic patterning phenotype including an absence of the division boundary and an anterior-posterior mispositioning of the new oral apparatus. CdaH is a Fused or Stk36 kinase ortholog that localizes to multiple sites that correlate with the effects of its loss, including the division boundary and the new oral apparatus. CdaH acts downstream of CdaA to induce the division boundary and drives asymmetric cytokinesis at the tip of the posterior daughter. CdaH both maintains the anterior-posterior position of the new oral apparatus and interacts with CdaI to pattern ciliary rows within the oral apparatus. Thus, CdaH acts at multiple scales, from induction and positioning of structures on the cell-wide polarity axis to local organelle-level patterning.

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.

MarFERReT, an open-source, version-controlled reference library of marine microbial eukaryote functional genes

Metatranscriptomics generates large volumes of sequence data about transcribed genes in natural environments. Taxonomic annotation of these datasets depends on availability of curated reference sequences. For marine microbial eukaryotes, current reference libraries are limited by gaps in sequenced organism diversity and barriers to updating libraries with new sequence data, resulting in taxonomic annotation of about half of eukaryotic environmental transcripts. Here, we introduce Marine Functional EukaRyotic Reference Taxa (MarFERReT), a marine microbial eukaryotic sequence library designed for use with taxonomic annotation of eukaryotic metatranscriptomes. We gathered 902 publicly accessible marine eukaryote genomes and transcriptomes and assessed their sequence quality and cross-contamination issues, selecting 800 validated entries for inclusion in MarFERReT. Version 1.1 of MarFERReT contains reference sequences from 800 marine eukaryotic genomes and transcriptomes, covering 453 species- and strain-level taxa, totaling nearly 28 million protein sequences with associated NCBI and PR2 Taxonomy identifiers and Pfam functional annotations. The MarFERReT project repository hosts containerized build scripts, documentation on installation and use case examples, and information on new versions of MarFERReT.

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.

Three closely-related subclasses Phacodiniidia Small & Lynn, 1985, Protohypotrichia Shi et al., 1999, and Euplotia Jankowski, 1979 (Protista, Ciliophora): A new contribution to their phylogeny with reconsiderations on the evolutionary hypotheses

The huge variety of species and worldwide distribution of ciliated protists in class Spirotrichea continue to make it one of the most complicated and confused groups in Ciliophora, despite significant research interest in the unique molecular genetics of these organisms. In this study, the morphological and molecular information were integrated, and it is inferred from a new perspective for the evolutionary relationship among Phacodiniidia, Protohypotrichia, Hypotrichia and Euplotia. Our results indicate that Kiitricha and Caryotricha, two members in Protohypotrichia, may represent two parallel branches of evolution; Euplotidae and Aspidiscidae represent the most recently diverged taxa within Euplotida, followed by Certesiidae, Gastrocirrhidae, and Uronychidae. Further, representative morphological characters (e.g. fronto-ventral-transverse cirral anlagen, undulating membranes, marginal cirri and caudal cirri) were stochastically mapped on phylogenies to speculate evolutionary path and morphological characters of the evolutionary transition node groups were assumed.

SDRAP for annotating scrambled or rearranged genomes

Genomes sometimes undergo large-scale rearrangements. Programmed genome rearrangements in ciliates offer an extreme example, making them a compelling model system to study DNA rearrangements. Currently, available methods for genome annotation are not adequate for highly scrambled genomes. We present a theoretical framework and software implementation for the systematic extraction and analysis of DNA rearrangement annotations from pairs of genome assemblies corresponding to precursor and product versions. The software makes no assumptions about the structure of the rearrangements, and permits the user to select parameters to suit the data. Compared to previous approaches, this work achieves more complete precursor-product mappings, allows for full transparency and reproducibility, and can be adapted to genomic data from different sources.

Identification of a non-canonical ciliate nuclear genetic code where UAA and UAG code for different amino acids

The genetic code is one of the most highly conserved features across life. Only a few lineages have deviated from the "universal" genetic code. Amongst the few variants of the genetic code reported to date, the codons UAA and UAG virtually always have the same translation, suggesting that their evolution is coupled. Here, we report the genome and transcriptome sequencing of a novel uncultured ciliate, belonging to the Oligohymenophorea class, where the translation of the UAA and UAG stop codons have changed to specify different amino acids. Genomic and transcriptomic analyses revealed that UAA has been reassigned to encode lysine, while UAG has been reassigned to encode glutamic acid. We identified multiple suppressor tRNA genes with anticodons complementary to the reassigned codons. We show that the retained UGA stop codon is enriched in the 3'UTR immediately downstream of the coding region of genes, suggesting that there is functional drive to maintain tandem stop codons. Using a phylogenomics approach, we reconstructed the ciliate phylogeny and mapped genetic code changes, highlighting the remarkable number of independent genetic code changes within the Ciliophora group of protists. According to our knowledge, this is the first report of a genetic code variant where UAA and UAG encode different amino acids.

Genome-wide identification of ATP-binding cassette transporter B subfamily, focusing on its structure, evolution and rearrangement in ciliates

ATP-binding cassette subfamily B (ABCB) has been implicated in various essential functions such as multidrug resistance, auxin transport and heavy metal tolerance in animals and plants. However, the functions, the genomic distribution and the evolutionary history have not been characterized systematically in lower eukaryotes. As a lineage of highly specialized unicellular eukaryotes, ciliates have extremely diverse genomic features including nuclear dimorphism. To further understand the genomic structure and evolutionary history of this gene family, we investigated the ABCB gene subfamily in 11 ciliates. The results demonstrate that there is evidence of substantial gene duplication, which has occurred by different mechanisms in different species. These gene duplicates show consistent purifying selection, suggesting functional constraint, in all but one species, where positive selection may be acting to generate novel function. We also compare the gene structures in the micronuclear and macronuclear genomes and find no gene scrambling during genome rearrangement, despite the abundance of such scrambling in two of our focal species. These results lay the foundation for future analyses of the function of these genes and the mechanisms responsible for their evolution across diverse eukaryotic lineages.

Comparative Transcriptomics of Chilodonella hexasticha and C. uncinata Provide New Insights into Adaptations to a Parasitic Lifestyle and Mdivi-1 as a Potential Agent for Chilodonellosis Control

Chilodonella hexasticha is a harmful parasitic ciliate that can cause severe damage to fish and high mortalities worldwide. Its congeneric species, C. uncinata, is a facultative parasite that not only can be free-living but also can parasitize on fish gills and fins. In this study, single-cell transcriptomes of these two species were assembled and characterized. Numerous enzymes related to energy metabolism and parasitic adaption were identified through annotation in the Non-Redundant (NR), Clusters of Orthologous Genes (COG), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. The expression of isocitrate dehydrogenase (IDH), cytochrome c oxidase subunit 1 (Cox1) and ATP synthase F1, delta subunit (ATP5D) was up-regulated in C. hexasticha compared with C. uncinata. The oxidative phosphorylation process was also enriched in C. hexasticha. The main mitochondrial metabolic pathways in C. hexasticha were depicted and enzymes related to energy metabolism pathways were compared between these two species. More importantly, mitochondrial division inhibitor 1 (mdivi-1) proved to be very effective in killing both C. hexasticha and C. uncinata, which could be a novel drug for Chilodonellosis control. This study can help us better understand the energy metabolisms of C. hexasticha and C. uncinata and provide new insight into novel targets for chilodonellosis control. Meanwhile, the transcriptome data can also facilitate genomic studies of these two species in the future.

Comparative genome analysis of three euplotid protists provides insights into the evolution of nanochromosomes in unicellular eukaryotic organisms

One of the most diverse clades of ciliated protozoa, the class Spirotrichea, displays a series of unique characters in terms of eukaryotic macronuclear (MAC) genome, including high fragmentation that produces nanochromosomes. However, the genomic diversity and evolution of nanochromosomes and gene families for spirotrich MAC genomes are poorly understood. In this study, we assemble the MAC genome of a representative euplotid (a new model organism in Spirotrichea) species, Euplotes aediculatus. Our results indicate that: (a) the MAC genome includes 35,465 contigs with a total length of 97.3 Mb and a contig N50 of 3.4 kb, and contains 13,145 complete nanochromosomes and 43,194 predicted genes, with the majority of these nanochromosomes containing tiny introns and harboring only one gene; (b) genomic comparisons between E. aediculatus and other reported spirotrichs indicate that average GC content and genome fragmentation levels exhibit interspecific variation, and chromosome breaking sites (CBSs) might be lost during evolution, resulting in the increase of multi-gene nanochromosome; (c) gene families associated with chitin metabolism and FoxO signaling pathway are expanded in E. aediculatus, suggesting their potential roles in environment adaptation and survival strategies of E. aediculatus; and (d) a programmed ribosomal frameshift (PRF) with a conservative motif 5'-AAATAR-3' tends to occur in longer genes with more exons, and PRF genes play an important role in many cellular regulation processes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00175-0.

Nontriplet feature of genetic code in Euplotes ciliates is a result of neutral evolution

The triplet nature of the genetic code is considered a universal feature of known organisms. However, frequent stop codons at internal mRNA positions in Euplotes ciliates ultimately specify ribosomal frameshifting by one or two nucleotides depending on the context, thus posing a nontriplet feature of the genetic code of these organisms. Here, we sequenced transcriptomes of eight Euplotes species and assessed evolutionary patterns arising at frameshift sites. We show that frameshift sites are currently accumulating more rapidly by genetic drift than they are removed by weak selection. The time needed to reach the mutational equilibrium is several times longer than the age of Euplotes and is expected to occur after a several-fold increase in the frequency of frameshift sites. This suggests that Euplotes are at an early stage of the spread of frameshifting in expression of their genome. In addition, we find the net fitness burden of frameshift sites to be noncritical for the survival of Euplotes. Our results suggest that fundamental genome-wide changes such as a violation of the triplet character of genetic code can be introduced and maintained solely by neutral evolution.

Four ciliate-specific expansion events occurred during actin gene family evolution of eukaryotes

Actin gene family is a divergent and ancient eukaryotic cellular cytoskeletal gene family, and participates in many essential cellular processes. Ciliated protists offer us an excellent opportunity to investigate gene family evolution, since their gene families evolved faster in ciliates than in other eukaryotes. Nonetheless, actin gene family is well studied in few model ciliate species but little is known about its evolutionary patterns in ciliates. Here, we analyzed the evolutionary pattern of eukaryotic actin gene family based on genomes/transcriptomes of 36 species covering ten ciliate classes, as well as those of nine non-ciliate eukaryotic species. Results showed: (1) Except for conventional actins and actin-related proteins (Arps) shared by various eukaryotes, at least four ciliate-specific subfamilies occurred during evolution of actin gene family. Expansions of Act2 and ArpC were supposed to have happen in the ciliate common ancestor, while expansions of ActI and ActII may have occurred in the ancestor of Armophorea, Muranotrichea, and Spirotrichea. (2) The number of actin isoforms varied greatly among ciliate species. Environmental adaptability, whole genome duplication (WGD) or segmental duplication events, distinct spatial and temporal patterns of expression might play driving forces for the increasement of isoform numbers. (3) The 'birth and death' model of evolution could explain the evolution of actin gene family in ciliates. And actin genes have been generally under strong negative selection to maintain protein structures and physiological functions. Collectively, we provided meaningful information for understanding the evolution of eukaryotic actin gene family.

Stop or Not: Genome-Wide Profiling of Reassigned Stop Codons in Ciliates

Bifunctional stop codons that have both translation and termination functions in the same species are important for understanding the evolution and function of genetic codes in living organisms. Considering the high frequency of bifunctional codons but limited number of available genomes in ciliates, we de novo sequenced seven representative ciliate genomes to explore the evolutionary history of stop codons. We further propose a stop codon reassignment quantification method (stopCR) that can identify bifunctional codons and measure their frequencies in various eukaryotic organisms. Using our newly developed method, we found two previously undescribed genetic codes, illustrating the prevalence of bifunctional stop codons in ciliates. Overall, evolutionary genomic analyses suggest that gain or loss of reassigned stop codons in ciliates is shaped by their living environment, the eukaryotic release factor 1, and suppressor tRNAs. This study provides novel clues about the functional diversity and evolutionary history of stop codons in eukaryotic organisms.

Novel insights into molecular mechanisms of vegetative cell cycle and resting cyst formation in Apodileptus cf. visscheri (Alveolata, Ciliophora)

Ciliates usually with big cell sizes, complex morphological structures, and diverse life cycles, are good model organisms for studying cell proliferation regulation of eukaryotes. Up to date, the molecular regulation mechanisms for the vegetative cell cycle and encystment of these ciliates are poorly understood. Here, transcriptomes of Apodileptus cf. visscheri, which has an asexual vegetative cell cycle and is apt to encyst when environmental conditions become unfavorable, were sequenced to enrich our related knowledge. In this study, three replicates were sequenced for each of four cell stages, including initial period of growth, morphogenesis, cell division, and resting cyst. The significant transcription differences, involving cell cycle, biosynthesis, and energy metabolism pathways, were revealed between the resting cyst and vegetative cell cycle. Further investigations showed that the cell cycle pathway was enriched during morphogenesis stage and cell division stage. Compared to the initial period of growth stage, the differentially expressed genes involved in cellular components and molecular function were significantly enriched during cell division stage, while cellular components and biological processes were significantly enriched during morphogenesis stage. These provide novel insights into a comprehensive understanding at the molecular level of the survival and adaptive mechanism of unicellular eukaryotes.

A proposed metric set for evaluation of genome assembly quality

Quality control is essential for genome assemblies; however, a consensus has yet to be reached on what metrics should be adopted for the evaluation of assembly quality. N50 is widely used for contiguity measurement, but its effectiveness is constantly in question. Prevailing metrics for the completeness evaluation focus on gene space, yet challenging areas such as tandem repeats are commonly overlooked. Achieving correctness has become an indispensable dimension for quality control, while prevailing assembly releases lack scores reflecting this aspect. We propose a metric set with a set of statistic indexes for effective, comprehensive evaluation of assemblies and provide a score of a finished assembly for each metric, which can be utilized as a benchmark for achieving high-quality genome assemblies.

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.

Exploration of the Nuclear Proteomes in the Ciliate Oxytricha trifallax

Nuclear dimorphism is a fundamental feature of ciliated protozoa, which have separate somatic and germline genomes in two distinct organelles within a single cell. The transcriptionally active somatic genome, contained within the physically larger macronucleus, is both structurally and functionally different from the silent germline genome housed in the smaller micronucleus. This difference in genome architecture is particularly exaggerated in Oxytricha trifallax, in which the somatic genome comprises tens of thousands of gene-sized nanochromosomes maintained at a high and variable ploidy, while the germline has a diploid set of megabase-scale chromosomes. To examine the compositional differences between the nuclear structures housing the genomes, we performed a proteomic survey of both types of nuclei and of macronuclear histones using quantitative mass spectrometry. We note distinct differences between the somatic and germline nuclei, with many functional proteins being highly enriched in one of the two nuclei. To validate our conclusions and the efficacy of nuclear separation, we used protein localization through a combination of transformations and immunofluorescence. We also note that the macronuclear histones strikingly display only activating marks, consistent with the conclusion that the macronucleus is the hub of transcription. These observations suggest that the compartmentalization of different genome features into separate structures has been accompanied by a similar specialization of nuclear components that maintain and facilitate the functions of the genomes specific to each nucleus.

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.

Genomic insights into the phylogeny and biomass-degrading enzymes of rumen ciliates

Understanding the biodiversity and genetics of gut microbiomes has important implications for host physiology and industrial enzymes, whereas most studies have been focused on bacteria and archaea, and to a lesser extent on fungi and viruses. One group, still underexplored and elusive, is ciliated protozoa, despite its importance in shaping microbiota populations. Integrating single-cell sequencing and an assembly-and-identification pipeline, we acquired 52 high-quality ciliate genomes of 22 rumen morphospecies from 11 abundant morphogenera. With these genomes, we resolved the taxonomic and phylogenetic framework that revised the 22 morphospecies into 19 species spanning 13 genera and reassigned the genus Dasytricha from Isotrichidae to a new family Dasytrichidae. Comparative genomic analyses revealed that extensive horizontal gene transfers and gene family expansion provided rumen ciliate species with a broad array of carbohydrate-active enzymes (CAZymes) to degrade all major kinds of plant and microbial carbohydrates. In particular, the genomes of Diplodiniinae and Ophryoscolecinae species encode as many CAZymes as gut fungi, and ~80% of their degradative CAZymes act on plant cell-wall. The activities of horizontally transferred cellulase and xylanase of ciliates were experimentally verified and were 2-9 folds higher than those of the inferred corresponding bacterial donors. Additionally, the new ciliate dataset greatly facilitated rumen metagenomic analyses by allowing ~12% of the metagenomic sequencing reads to be classified as ciliate sequences.

Comparative genomics reveals insight into the evolutionary origin of massively scrambled genomes

Ciliates are microbial eukaryotes that undergo extensive programmed genome rearrangement, a natural genome editing process that converts long germline chromosomes into smaller gene-rich somatic chromosomes. Three well-studied ciliates include Oxytricha trifallax, Tetrahymena thermophila, and Paramecium tetraurelia, but only the Oxytricha lineage has a massively scrambled genome, whose assembly during development requires hundreds of thousands of precisely programmed DNA joining events, representing the most complex genome dynamics of any known organism. Here we study the emergence of such complex genomes by examining the origin and evolution of discontinuous and scrambled genes in the Oxytricha lineage. This study compares six genomes from three species, the germline and somatic genomes for Euplotes woodruffi, Tetmemena sp., and the model ciliate O. trifallax. We sequenced, assembled, and annotated the germline and somatic genomes of E. woodruffi, which provides an outgroup, and the germline genome of Tetmemena sp. We find that the germline genome of Tetmemena is as massively scrambled and interrupted as Oxytricha's: 13.6% of its gene loci require programmed translocations and/or inversions, with some genes requiring hundreds of precise gene editing events during development. This study revealed that the earlier diverged spirotrich, E. woodruffi, also has a scrambled genome, but only roughly half as many loci (7.3%) are scrambled. Furthermore, its scrambled genes are less complex, together supporting the position of Euplotes as a possible evolutionary intermediate in this lineage, in the process of accumulating complex evolutionary genome rearrangements, all of which require extensive repair to assemble functional coding regions. Comparative analysis also reveals that scrambled loci are often associated with local duplications, supporting a gradual model for the origin of complex, scrambled genomes via many small events of DNA duplication and decay.

Evolutionary bioenergetics of ciliates

Understanding why various organisms evolve alternative ways of living requires information on both the fitness advantages of phenotypic modifications and the costs of constructing and operating cellular features. Although the former has been the subject of a myriad of ecological studies, almost no attention has been given to how organisms allocate resources to alternative structures and functions. We address these matters by capitalizing on an array of observations on diverse ciliate species and from the emerging field of evolutionary bioenergetics. A relatively robust and general estimator for the total cost of a cell per cell cycle (in units of ATP equivalents) is provided, and this is then used to understand how the magnitudes of various investments scale with cell size. Among other things, we examine the costs associated with the large macronuclear genomes of ciliates, as well as ribosomes, various internal membranes, osmoregulation, cilia, and swimming activities. Although a number of uncertainties remain, the general approach taken may serve as blueprint for expanding this line of work to additional traits and phylogenetic lineages.

Timing and characteristics of nuclear events during conjugation and genomic exclusion in Paramecium multimicronucleatum

Ciliated protists are ideal material for studying the origin and evolution of sex, because of their nuclear dimorphism (containing both germline micronucleus and somatic macronucleus in the same cytoplasm), special sexual processes (conjugation and autogamy), and high diversity of mating-type systems. However, the study of sexual process is limited to only a few species, due to the difficulties in inducing or observing conjugation. In the present study, we investigate the conjugation process in Paramecium multimicronucleatum: (1) of the three prezygotic divisions, all micronuclei undergo the first two divisions (meiosis I, II), while a variable number of nuclei undergo the third division (mitosis); (2) the synkaryon divides three times after fertilization, giving rise to eight products that differentiate into four macronuclear anlagen and four micronuclei; (3) cells restore the vegetative stage after two successive cell fissions during which the macronuclear anlagen are distributed into daughter cells without division, while micronuclei divide mitotically; (4) the parental macronucleus begins to fragment following the first meiotic division and finally degenerates completely; (5) the entire process takes about 110 h, of which about 85 h are required for macronuclear development. In addition, we describe for the first time the process of genomic exclusion occurring between amicronucleate and micronucleate cells of P. multimicronucleatum, during which the micronucleate cell contributes a pronucleus to the amicronucleate cell, resulting in both exconjugants being homozygotes. These results provide new insights into the diversity of sexual processes and lay an important cytological basis for future in-depth studies of mating systems in ciliates.

Genome architecture used to supplement species delineation in two cryptic marine ciliates

The purpose of this study is to determine which taxonomic methods can elucidate clear and quantifiable differences between two cryptic ciliate species, and to test the utility of genome architecture as a new diagnostic character in the discrimination of otherwise indistinguishable taxa. Two cryptic tintinnid ciliates, Schmidingerella arcuata and Schmidingerella meunieri, are compared via traditional taxonomic characters including lorica morphometrics, ribosomal RNA (rRNA) gene barcodes and ecophysiological traits. In addition, single-cell 'omics analyses (single-cell transcriptomics and genomics) are used to elucidate and compare patterns of micronuclear genome architecture between the congeners. The results include a highly similar lorica that is larger in S. meunieri, a 0%-0.5% difference in rRNA gene barcodes, two different and nine indistinguishable growth responses among 11 prey treatments, and distinct patterns of micronuclear genomic architecture for genes detected in both ciliates. Together, these results indicate that while minor differences exist between S. arcuata and S. meunieri in common indices of taxonomic identification (i.e., lorica morphology, DNA barcode sequences and ecophysiology), differences exist in their genomic architecture, which suggests potential genetic incompatibility. Different patterns of micronuclear architecture in genes shared by both isolates also enable the design of species-specific primers, which are used in this study as unique "architectural barcodes" to demonstrate the co-occurrence of both ciliates in samples collected from a NW Atlantic estuary. These results support the utility of genomic architecture as a tool in species delineation, especially in ciliates that are cryptic or otherwise difficult to differentiate using traditional methods of identification.

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.

Molecular evolutionary analysis of the SM and SNARE vesicle fusion machinery in ciliates shows concurrent expansions in late secretory machinery

Protists in the phylum Ciliophora possess a complex membrane-trafficking system, including osmoregulatory Contractile Vacuoles and specialized secretory organelles. Molecular cell biological investigations in Tetrahymena thermophila have identified components of the protein machinery associated with the secretory organelles, mucocysts. The Qa-SNARE Syn7lp plays a role in mucocyst biogenesis as do subunits of the CORVET tethering complex (specifically Vps8). Indeed, Tetrahymena thermophila possesses expanded gene complements of several CORVET components, including Vps33 which is also a Sec1/Munc18 (SM) protein that binds Qa-SNAREs. Moreover, the Qa-SNAREs in Paramecium tetraurelia have been localized to various endomembrane organelles. Here, we use comparative genomics and phylogenetics to determine the evolutionary history of the SM and Qa-SNARE proteins across the Ciliophora. We identify that the last ciliate common ancestor possessed the four SM proteins and six Qa-SNAREs common to eukaryotes, including the uncommonly retained Syntaxin 17. We furthermore identify independent expansion of these protein families in several ciliate classes, including concurrent expansions of the SM protein-Qa SNARE partners Sec1:SynPM in the oligohymenophorean ciliates lineage, consistent with novel Contractile Vacuole specific innovations. Overall, these data are consistent with SM proteins and Qa-SNAREs being a common set of components for endomembrane modulation in the ciliates.

Paramecium epigenetics in development and proliferation

The term epigenetics is used for any layer of genetic information aside from the DNA base-sequence information. Mammalian epigenetic research increased our understanding of chromatin dynamics in terms of cytosine methylation and histone modification during differentiation, aging, and disease. Instead, ciliate epigenetics focused more on small RNA-mediated effects. On the one hand, these do concern the transport of RNA from parental to daughter nuclei, representing a regulated transfer of epigenetic information across generations. On the other hand, studies of Paramecium, Tetrahymena, Oxytricha, and Stylonychia revealed an almost unique function of transgenerational RNA. Rather than solely controlling chromatin dynamics, they control sexual progeny's DNA content quantitatively and qualitatively. Thus epigenetics seems to control genetics, at least genetics of the vegetative macronucleus. This combination offers ciliates, in particular, an epigenetically controlled genetic variability. This review summarizes the epigenetic mechanisms that contribute to macronuclear heterogeneity and relates these to nuclear dimorphism. This system's adaptive and evolutionary possibilities raise the critical question of whether such a system is limited to unicellular organisms or binuclear cells. We discuss here the relevance of ciliate genetics and epigenetics to multicellular organisms.

An extremely streamlined macronuclear genome in the free-living protozoan Fabrea salina

Ciliated protists are among the oldest unicellular organisms with a heterotrophic lifestyle and share a common ancestor with Plantae. Unlike any other eukaryotes, there are two distinct nuclei in ciliates with separate germline and somatic cell functions. Here, we assembled a near-complete macronuclear genome of Fabrea salina, which belongs to one of the oldest clades of ciliates. Its extremely minimized genome (18.35 Mb) is the smallest among all free-living heterotrophic eukaryotes and exhibits typical streamlined genomic features, including high gene density, tiny introns, and shrinkage of gene paralogs. Gene families involved in hypersaline stress resistance, DNA replication proteins, and mitochondrial biogenesis are expanded, and the accumulation of phosphatidic acid may play an important role in resistance to high osmotic pressure. We further investigated the morphological and transcriptomic changes in the macronucleus during sexual reproduction and highlighted the potential contribution of macronuclear residuals to this process. We believe that the minimized genome generated in this study provides novel insights into the genome streamlining theory and will be an ideal model to study the evolution of eukaryotic heterotrophs.

Comparative Genome Analysis Reveals Cis-Regulatory Elements on Gene-Sized Chromosomes of Ciliated Protists

Gene-sized chromosomes are a distinct feature of the macronuclear genome in ciliated protists known as spirotrichs. These nanochromosomes are often only several kilobase pairs long and contain a coding region for a single gene. However, the ways in which transcription is regulated on nanochromosomes is still largely unknown. Here, we generated macronuclear genome assemblies for two species of Pseudokeronopsis ciliates to better understand transcription regulation on gene-sized chromosomes. We searched within the short subtelomeric regions for potential cis-regulatory elements and identified distinct AT-rich sequences conserved in both species, at both the 5’ and 3’ end of each gene. We further acquired transcriptomic data for these species, which showed the 5’ cis-regulatory element is associated with active gene expression. Gene family evolution analysis suggests nanochromosomes in spirotrichs may originated approximately 900 million years ago. Together our comparative genomic analyses reveal novel insights into the biological roles of cis-regulatory elements on gene-sized chromosomes.

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.

An approach to determining the nitrofurazone-induced toxic dynamics for ecotoxicity assessment using protozoan periphytons in marine ecosystems

With several observable responses and sensitivity of protozoans to nitrofurazone (NFZ), the toxic effects of NFZ on protozoans can be an early warning signal of NFZ contamination in the aquatic environment. To evaluate the toxic dynamics induced by NFZ, protozoan samples were collected using microscopy glass slides and exposed to the five concentrations of NFZ: 0, 1, 2, 4, and 8 mg ml^-1. Substantial differences in the species composition and toxic-dynamics patterns were observed among all concentrations. Briefly, periphytic euplotids and pleurostomatids were the most prevalent at each concentration level, while dysteriids were less dominant among all treatments. Multivariate analysis revealed significant (P < 0.05) differences in the taxonomic patterns of the test organisms among the five treatments. Furthermore, significant deviation of protozoan communities from the expected taxonomic breadth was observed to occur in a dose-dependent manner. Based on these findings, it is suggested that protozoan periphytons could be used as bioindicators to assess the ecotoxicity of NFZ in the marine environment.

Dependence of genome size and copy number of rRNA gene on cell volume in dinoflagellates

Dinoflagellates are an ecologically important group of protists in aquatic environment and have evolved many unusual and enigmatic genomic features such as immense genome sizes, high repeated genes, and a large portion of hydroxymethyluracil in DNA. Although previous studies have observed positive correlations between the large subunit (LSU) rRNA gene copy number and genome size of a variety of eukaryotic organisms (e.g. higher plants and animals), or between cell volume and LSU rRNA gene copy number, and/or between genome size and cell size, which suggests a possible co-evolution among these three features in different lineages of life, it remains an open question regarding the relationships among these three parameters in dinoflagellates. For the first time, we estimated the copy numbers of the LSU rRNA gene, the genome sizes, and cell volumes within a broad range of dinoflagellates (covering 15 species of 11 genera) using single-cell qPCR-based assay (determining LSU rRNA gene copy number), FlowCAM (cell volume measurement), and ultraviolet spectrophotometry (genome size estimation). The measured copy number of LSU rRNA gene ranged from 398 ± 184 (Prorocentrum minimum) to 152,078 ± 33,555 copies•cell^-1 (Alexandrium pacificum), while the genome size and the cell volume ranged from 5.6 ± 0.2 (Karlodinium veneficum) to 853 ± 19.9 pg•cell^-1 (Pseliodinium pirum), and from 1,070 ± 225 (Kar. veneficum) to 168,474 ± 124,180 μm³ (Ps. pirum), respectively. Together with the three parameters measured in literature, there are significant positive linear correlations between LSU rRNA gene copy numbers and genome sizes, cell volumes and LSU rRNA gene copy numbers, and between genome sizes and cell volumes via comparisons of multi-model regression analyses, suggesting a dependence of genome size and rRNA gene copy number on the cell volumes of dinoflagellates. Validation of the measurement methods was conducted via comparisons between reported data in the literature and that predicted using the linear equations we obtained, and between genome size measured by flow cytometry (FCM) and ultraviolet spectrophotometry (Nanodrop). These results provide insightful understandings of dinoflagellate evolution in terms of the relationships among genomes, gene copy number, and cell volume, and of rRNA gene-based studies in intra-populational and intra-individual genetic diversity, taxonomy, and diversity assessment in the environment of dinoflagellates. The results also provide a dataset useful for reads calibration in environmental metabarcoding studies of dinoflagellates and selection of candidate species for whole genome sequencing.

The 3D Genome: From Structure to Function

The genome is the most functional part of a cell, and genomic contents are organized in a compact three-dimensional (3D) structure. The genome contains millions of nucleotide bases organized in its proper frame. Rapid development in genome sequencing and advanced microscopy techniques have enabled us to understand the 3D spatial organization of the genome. Chromosome capture methods using a ligation approach and the visualization tool of a 3D genome browser have facilitated detailed exploration of the genome. Topologically associated domains (TADs), lamin-associated domains, CCCTC-binding factor domains, cohesin, and chromatin structures are the prominent identified components that encode the 3D structure of the genome. Although TADs are the major contributors to 3D genome organization, they are absent in Arabidopsis. However, a few research groups have reported the presence of TAD-like structures in the plant kingdom.

Chromosome organization and gene expansion in the highly fragmented genome of the ciliate Strombidium stylifer

Chromosomes are well-organized carriers of genetic information in eukaryotes and are usually quite long, carrying hundreds and thousands of genes. Intriguingly, a clade of single-celled ciliates, Spirotrichea, feature nanochromosomes-also called "gene-sized chromosomes". These chromosomes predominantly carry only one gene, flanked by short telomere sequences. However, the organization and copy number variation of the chromosomes in these highly fragmented genomes remain unexplored in many groups of Spirotrichea, including the marine Strombidium. Using deep genome sequencing, we assembled the macronuclear genome of Strombidium stylifer into more than 18,000 nanochromosomes (~2.4 Kb long on average). Our results show that S. stylifer occupies an intermediate position during the evolutionary history of Strombidium lineage and experienced significant expansions in several gene families related to guanyl ribonucleotide binding. Based on the nucleotide distribution bias analysis and conserved motifs search in non-genic regions, we found that the subtelomeric regions have a conserved adenine-thymine (AT)-rich sequence motif. We also found that the copy number of nanochromosomes lacks precise regulation. This work sheds light on the unique features of chromosome structure in eukaryotes with highly fragmented genomes and reveals that a rather specialized evolutionary strategy at the genomic level has resulted in great diversity within the ciliated lineages.

DAPI staining and DNA content estimation of nuclei in uncultivable microbial eukaryotes (Arcellinida and Ciliates)

Though representing a major component of eukaryotic biodiversity, many microbial eukaryotes remain poorly studied, including the focus of the present work, testate amoebae of the order Arcellinida (Amoebozoa) and non-model lineages of ciliates (Alveolata). In particular, knowledge of genome structures and changes in genome content over the often-complex life cycles of these lineages remains enigmatic. However, the limited available knowledge suggests that microbial eukaryotes have the potential to challenge our textbook views on eukaryotic genomes and genome evolution. In this study, we developed protocols for DAPI (4',6-diamidino-2-phenylindole) staining of Arcellinida nuclei and adapted protocols for ciliates. In addition, image analysis software was used to estimate the DNA content in the nuclei of Arcellinida and ciliates, and the measurements of target organisms were compared to those  of well-known model organisms.The results demonstrate that the methods we have developed for nuclear staining in these lineages are effective and can be applied to other microbial eukaryotic groups by adjusting certain stages in the protocols.

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.

Fifty Generations of Amitosis: Tracing Asymmetric Allele Segregation in Polyploid Cells with Single-Cell DNA Sequencing

Amitosis is a widespread form of unbalanced nuclear division whose biomedical and evolutionary significance remain unclear. Traditionally, insights into the genetics of amitosis have been gleaned by assessing the rate of phenotypic assortment. Though powerful, this experimental approach relies on the availability of phenotypic markers. Leveraging Paramecium tetraurelia, a unicellular eukaryote with nuclear dualism and a highly polyploid somatic nucleus, we probe the limits of single-cell whole-genome sequencing to study the consequences of amitosis. To this end, we first evaluate the suitability of single-cell sequencing to study the AT-rich genome of P. tetraurelia, focusing on common sources of genome representation bias. We then asked: can alternative rearrangements of a given locus eventually assort after a number of amitotic divisions? To address this question, we track somatic assortment of developmentally acquired Internal Eliminated Sequences (IESs) up to 50 amitotic divisions post self-fertilization. To further strengthen our observations, we contrast empirical estimates of IES retention levels with in silico predictions obtained through mathematical modeling. In agreement with theoretical expectations, our empirical findings are consistent with a mild increase in variation of IES retention levels across successive amitotic divisions of the macronucleus. The modest levels of somatic assortment in P. tetraurelia suggest that IESs retention levels are largely sculpted at the time of macronuclear development, and remain fairly stable during vegetative growth. In forgoing the requirement for phenotypic assortment, our approach can be applied to a wide variety of amitotic species and could facilitate the identification of environmental and genetic factors affecting amitosis.

Cancer progression as a sequence of atavistic reversions

It has long been recognized that cancer onset and progression represent a type of reversion to an ancestral quasi-unicellular phenotype. This general concept has been refined into the atavistic model of cancer that attempts to provide a quantitative analysis and testable predictions based on genomic data. Over the past decade, support for the multicellular-to-unicellular reversion predicted by the atavism model has come from phylostratigraphy. Here, we propose that cancer onset and progression involve more than a one-off multicellular-to-unicellular reversion, and are better described as a series of reversionary transitions. We make new predictions based on the chronology of the unicellular-eukaryote-to-multicellular-eukaryote transition. We also make new predictions based on three other evolutionary transitions that occurred in our lineage: eukaryogenesis, oxidative phosphorylation and the transition to adaptive immunity. We propose several modifications to current phylostratigraphy to improve age resolution to test these predictions. Also see the video abstract here: https://youtu.be/3unEu5JYJrQ.

Transcribed germline-limited coding sequences in Oxytricha trifallax

The germline-soma divide is a fundamental distinction in developmental biology, and different genes are expressed in germline and somatic cells throughout metazoan life cycles. Ciliates, a group of microbial eukaryotes, exhibit germline-somatic nuclear dimorphism within a single cell with two different genomes. The ciliate Oxytricha trifallax undergoes massive RNA-guided DNA elimination and genome rearrangement to produce a new somatic macronucleus (MAC) from a copy of the germline micronucleus (MIC). This process eliminates noncoding DNA sequences that interrupt genes and also deletes hundreds of germline-limited open reading frames (ORFs) that are transcribed during genome rearrangement. Here, we update the set of transcribed germline-limited ORFs (TGLOs) in O. trifallax. We show that TGLOs tend to be expressed during nuclear development and then are absent from the somatic MAC. We also demonstrate that exposure to synthetic RNA can reprogram TGLO retention in the somatic MAC and that TGLO retention leads to transcription outside the normal developmental program. These data suggest that TGLOs represent a group of developmentally regulated protein-coding sequences whose gene expression is terminated by DNA elimination.

The macronuclear genome of anaerobic ciliate Entodinium caudatum reveals its biological features adapted to the distinct rumen environment

Entodinium caudatum is an anaerobic binucleated ciliate representing the most dominant protozoal species in the rumen. However, its biological features are largely unknown due to the inability to establish an axenic culture. In this study, we primally sequenced its macronucleus (MAC) genome to aid the understanding of its metabolism, physiology, ecology. We isolated the MAC of E. caudatum strain MZG-1 and sequenced the MAC genome using Illumina MiSeq, MinION, and PacBio RSII systems. De novo assembly of the MiSeq sequence reads followed with subsequent scaffolding with MinION and PacBio reads resulted in a draft MAC genome about 117 Mbp. A large number of carbohydrate-active enzymes were likely acquired through horizontal gene transfer. About 8.74% of the E. caudatum predicted proteome was predicted as proteases. The MAC genome of E. caudatum will help better understand its important roles in rumen carbohydrate metabolism, and interaction with other members of the rumen microbiome.

The Compact Macronuclear Genome of the Ciliate Halteria grandinella: A Transcriptome-Like Genome with 23,000 Nanochromosomes

How to achieve protein diversity by genome and transcriptome processing is essential for organismal complexity and adaptation. The present work identifies that the macronuclear genome of Halteria grandinella, a cosmopolitan unicellular eukaryote, is composed almost entirely of gene-sized nanochromosomes with extremely short nongenic regions.

How to achieve protein diversity by genome and transcriptome processing is essential for organismal complexity and adaptation. The present work identifies that the macronuclear genome of Halteria grandinella, a cosmopolitan unicellular eukaryote, is composed almost entirely of gene-sized nanochromosomes with extremely short nongenic regions. This challenges our usual understanding of chromosomal structure and suggests the possibility of novel mechvanisms in transcriptional regulation. Comprehensive analysis of multiple data sets reveals that Halteria transcription dynamics are influenced by: (i) nonuniform nanochromosome copy numbers correlated with gene-expression level; (ii) dynamic alterations at both the DNA and RNA levels, including alternative internal eliminated sequence (IES) deletions during macronucleus formation and large-scale alternative splicing in transcript maturation; and (iii) extremely short 5′ and 3′ untranslated regions (UTRs) and universal TATA box-like motifs in the compact 5′ subtelomeric regions of most chromosomes. This study broadens the view of ciliate biology and the evolution of unicellular eukaryotes, and identifies Halteria as one of the most compact known eukaryotic genomes, indicating that complex cell structure does not require complex gene architecture.

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).

The Function and Evolution of Motile DNA Replication Systems in Ciliates

DNA replication is a ubiquitous and conserved cellular process. However, regulation of DNA replication is only understood in a small fraction of organisms that poorly represent the diversity of genetic systems in nature. Here we used computational and experimental approaches to examine the function and evolution of one such system, the replication band (RB) in spirotrich ciliates, which is a localized, motile hub that traverses the macronucleus while replicating DNA. We show that the RB can take unique forms in different species, from polar bands to a "replication envelope," where replication initiates at the nuclear periphery before advancing inward. Furthermore, we identify genes involved in cellular transport, including calcium transporters and cytoskeletal regulators, that are associated with the RB and may be involved in its function and translocation. These findings highlight the evolution and diversity of DNA replication systems and provide insights into the regulation of nuclear organization and processes.

Morpholino-Mediated Knockdown of Ciliary Genes in Euplotes vannus, a Novel Marine Ciliated Model Organism

Cilia are highly conserved organelles present in almost all types of eukaryotic cells, and defects in cilia structure and/or function are related to many human genetic disorders. Single-celled ciliated protists, which possess diverse types of cilia, are remarkable model organisms for studying cilia structures and functions. Euplotes vannus is a representative ciliate with many intriguing features; for example, it possesses extensively fragmented somatic genomes and a high frequency of + 1 programmed ribosomal frameshifting. However, the molecular mechanisms underlying these remarkable traits remain largely unknown, mainly due to the lack of efficient genetic manipulation tools. Here, we describe the first application of a morpholino-based strategy to knockdown gene expression in E. vannus. Through interfering with the function of two ciliary genes, ZMYND10 and C21ORF59, we showed that these two genes are essential for the ciliary motility and proliferation of E. vannus cells. Strikingly, both ZMYND10- and C21ORF59-knockdown cells developed shorter cilia in the ventral cirri, a special type of ciliary tuft, suggesting a novel role for these genes in the regulation of cilia length. Our data provide a new method to explore gene function in E. vannus, which may help us to understand the functions of evolutionarily conserved cilia-related genes as well as other biological processes in this intriguing model.