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

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.

Noncanonical usage of stop codons in ciliates expands proteins with structurally flexible Q-rich motifs

Serine(S)/threonine(T)-glutamine(Q) cluster domains (SCDs), polyglutamine (polyQ) tracts and polyglutamine/asparagine (polyQ/N) tracts are Q-rich motifs found in many proteins. SCDs often are intrinsically disordered regions that mediate protein phosphorylation and protein-protein interactions. PolyQ and polyQ/N tracts are structurally flexible sequences that trigger protein aggregation. We report that due to their high percentages of STQ or STQN amino acid content, four SCDs and three prion-causing Q/N-rich motifs of yeast proteins possess autonomous protein expression-enhancing activities. Since these Q-rich motifs can endow proteins with structural and functional plasticity, we suggest that they represent useful toolkits for evolutionary novelty. Comparative Gene Ontology (GO) analyses of the near-complete proteomes of 26 representative model eukaryotes reveal that Q-rich motifs prevail in proteins involved in specialized biological processes, including Saccharomyces cerevisiae RNA-mediated transposition and pseudohyphal growth, Candida albicans filamentous growth, ciliate peptidyl-glutamic acid modification and microtubule-based movement, Tetrahymena thermophila xylan catabolism and meiosis, Dictyostelium discoideum development and sexual cycles, Plasmodium falciparum infection, and the nervous systems of Drosophila melanogaster, Mus musculus and Homo sapiens. We also show that Q-rich-motif proteins are expanded massively in 10 ciliates with reassigned TAAQ and TAGQ codons. Notably, the usage frequency of CAGQ is much lower in ciliates with reassigned TAAQ and TAGQ codons than in organisms with expanded and unstable Q runs (e.g. D. melanogaster and H. sapiens), indicating that the use of noncanonical stop codons in ciliates may have coevolved with codon usage biases to avoid triplet repeat disorders mediated by CAG/GTC replication slippage.

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.

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 analysis of anterior-posterior mRNA localization in Stentor coeruleus reveals a role for the microtubule cytoskeleton

Cells have complex and beautiful structures that are important for their function, but understanding the molecular mechanisms that produce these structures is a challenging problem due to the gap in size scales between molecular interactions and cellular structures. The giant ciliate Stentor coeruleus is a unicellular model organism whose large size, reproducible structure, and ability to heal wounds and regenerate has historically allowed the formation of structure in a single cell to be addressed using methods of experimental embryology. Such studies have shown that specific cellular structures, such as the oral apparatus, always form in specific regions of the cell, which raises the question: what is the source of positional information within this organism? By analogy with embryonic development, in which localized mRNA is often used to mark position, we asked whether position along the anterior-posterior axis of Stentor might be marked by specific regionalized mRNAs. By physically bisecting cells and conducting half-cell RNA sequencing, we were able to identify sets of messages enriched in either the anterior or posterior half. We repeated this analysis in cells in which a set of longitudinal microtubule bundles running down the whole length of the cell, known as KM-fibers, were disrupted by RNAi of b-tubulin. We found that many messages either lost their regionalized distribution or switched to an opposite distribution, such that anterior-enriched messages in control became posterior-enriched in the RNAi cells, or vice versa. This study indicates that mRNA can be regionalized within a single giant cell and that microtubules may play a role, possibly by serving as tracks for the movement of the messages.

The nuclear transport factor CSE1 drives macronuclear volume increase and macronuclear node coalescence in Stentor coeruleus

Stentor coeruleus provides a unique opportunity to study how cells regulate nuclear shape because its macronucleus undergoes a rapid, dramatic, and developmentally regulated shape change. We found that the volume of the macronucleus increases during coalescence, suggesting an inflation-based mechanism. When the nuclear transport factor, CSE1, is knocked down by RNAi, the shape and volume changes of the macronucleus are attenuated, and nuclear morphology is altered. CSE1 protein undergoes a dynamic relocalization correlated with nuclear shape changes, being mainly cytoplasmic prior to nuclear coalescence, and accumulating inside the macronucleus during coalescence. At the end of regeneration, CSE1 protein levels are reduced as the macronucleus returns to its pre-coalescence volume. We propose a model in which nuclear transport via CSE1 is required to increase the volume of the macronucleus, thereby decreasing the surface-to-volume ratio and driving coalescence of the nodes into a single mass.

Contrasting outcomes of genome reduction in mikrocytids and microsporidians

BACKGROUND

Intracellular symbionts often undergo genome reduction, losing both coding and non-coding DNA in a process that ultimately produces small, gene-dense genomes with few genes. Among eukaryotes, an extreme example is found in microsporidians, which are anaerobic, obligate intracellular parasites related to fungi that have the smallest nuclear genomes known (except for the relic nucleomorphs of some secondary plastids). Mikrocytids are superficially similar to microsporidians: they are also small, reduced, obligate parasites; however, as they belong to a very different branch of the tree of eukaryotes, the rhizarians, such similarities must have evolved in parallel. Since little genomic data are available from mikrocytids, we assembled a draft genome of the type species, Mikrocytos mackini, and compared the genomic architecture and content of microsporidians and mikrocytids to identify common characteristics of reduction and possible convergent evolution.

RESULTS

At the coarsest level, the genome of M. mackini does not exhibit signs of extreme genome reduction; at 49.7 Mbp with 14,372 genes, the assembly is much larger and gene-rich than those of microsporidians. However, much of the genomic sequence and most (8075) of the protein-coding genes code for transposons, and may not contribute much of functional relevance to the parasite. Indeed, the energy and carbon metabolism of M. mackini share several similarities with those of microsporidians. Overall, the predicted proteome involved in cellular functions is quite reduced and gene sequences are extremely divergent. Microsporidians and mikrocytids also share highly reduced spliceosomes that have retained a strikingly similar subset of proteins despite having reduced independently. In contrast, the spliceosomal introns in mikrocytids are very different from those of microsporidians in that they are numerous, conserved in sequence, and constrained to an exceptionally narrow size range (all 16 or 17 nucleotides long) at the shortest extreme of known intron lengths.

CONCLUSIONS

Nuclear genome reduction has taken place many times and has proceeded along different routes in different lineages. Mikrocytids show a mix of similarities and differences with other extreme cases, including uncoupling the actual size of a genome with its functional reduction.

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.

iGDP: An integrated genome decontamination pipeline for wild ciliated microeukaryotes

Ciliates are a large group of ubiquitous and highly diverse single-celled eukaryotes that play an essential role in the functioning of microbial food webs. However, their genomic diversity is far from clear due to the need to develop cultivation methods for most species, so most research is based on wild organisms that almost invariably contain contaminants. Here we establish an integrated Genome Decontamination Pipeline (iGDP) that combines homology search, telomere reads-assisted and clustering approaches to filter contaminated ciliate genome assemblies from wild specimens. We benchmarked the performance of iGDP using genomic data from a contaminated ciliate culture and the results showed that iGDP could recall 91.9% of the target sequences with 96.9% precision. We also used a synthetic dataset to offer guidelines for the application of iGDP in the removal of various groups of contaminants. Compared with several popular metagenome binning tools, iGDP could show better performance. To further validate the effectiveness of iGDP on real-world data, we applied it to decontaminate genome assemblies of three wild ciliate specimens and obtained their genomes with high quality comparable to that of previously well-studied model ciliate genomes. It is anticipated that the newly generated genomes and the established iGDP method will be valuable community resources for detailed studies on ciliate biodiversity, phylogeny, ecology and evolution. The pipeline (https://github.com/GWang2022/iGDP) can be implemented automatically to reduce manual filtering and classification and may be further developed to apply to other microeukaryotes.

Giant proteins in a giant cell: Molecular basis of ultrafast Ca2+-dependent cell contraction

The giant single-celled eukaryote, Spirostomum, exhibits one of the fastest movements in the biological world. This ultrafast contraction is dependent on Ca²⁺ rather than ATP and therefore differs to the actin-myosin system in muscle. We obtained the high-quality genome of Spirostomum minus from which we identified the key molecular components of its contractile apparatus, including two major Ca²⁺ binding proteins (Spasmin 1 and 2) and two giant proteins (GSBP1 and GSBP2), which act as the backbone and allow for the binding of hundreds of spasmins. The evidence suggests that the GSBP-spasmin protein complex is the functional unit of the mesh-like contractile fibrillar system, which, coupled with various other subcellular structures, provides the mechanism for repetitive ultrafast cell contraction and extension. These findings improve our understanding of the Ca²⁺-dependent ultrafast movement and provide a blueprint for future biomimicry, design, and construction of this kind of micromachine.

Developmentally Programmed Switches in DNA Replication: Gene Amplification and Genome-Wide Endoreplication in Tetrahymena

Locus-specific gene amplification and genome-wide endoreplication generate the elevated copy number of ribosomal DNA (rDNA, 9000 C) and non-rDNA (90 C) chromosomes in the developing macronucleus of Tetrahymena thermophila. Subsequently, all macronuclear chromosomes replicate once per cell cycle during vegetative growth. Here, we describe an unanticipated, programmed switch in the regulation of replication initiation in the rDNA minichromosome. Early in development, the 21 kb rDNA minichromosome is preferentially amplified from 2 C to ~800 C from well-defined origins, concurrent with genome-wide endoreplication (2 C to 8-16 C) in starved mating Tetrahymena (endoreplication (ER) Phase 1). Upon refeeding, rDNA and non-rDNA chromosomes achieve their final copy number through resumption of just the endoreplication program (ER Phase 2). Unconventional rDNA replication intermediates are generated primarily during ER phase 2, consistent with delocalized replication initiation and possible formation of persistent RNA-DNA hybrids. Origin usage and replication fork elongation are affected in non-rDNA chromosomes as well. Despite the developmentally programmed 10-fold reduction in the ubiquitous eukaryotic initiator, the Origin Recognition Complex (ORC), active initiation sites are more closely spaced in ER phases 1 and 2 compared to vegetative growing cells. We propose that initiation site selection is relaxed in endoreplicating macronuclear chromosomes and may be less dependent on ORC.

Genome-Wide Identification of G Protein-Coupled Receptors in Ciliated Eukaryotes

G protein-coupled receptors (GPCRs) are the largest family of transmembrane receptors and play important roles in many physiological processes. As a representative group of protozoa, ciliates represent the highest stage of eukaryotic cell differentiation and evolution in terms of their reproductive mode, two-state karyotype, and extremely diverse cytogenesis patterns. GPCRs have been poorly reported in ciliates. In this study, we identified 492 GPCRs in 24 ciliates. Using the existing classification system for animals, GPCRs in ciliates can be assigned to four families, including families A, B, E, and F. Most (377 members) belong to family A. The number of GPCRs is extremely different in different ciliates; the Heterotrichea ciliates usually have more GPCRs than other ciliates. Parasitic or symbiotic ciliates usually have only a few GPCRs. Gene/genome duplication events seem to play important roles in the expansion of the GPCR superfamily in ciliates. GPCRs in ciliates displayed seven typical domain organizations. GPCRs in an ortholog group are common and conserved in all ciliates. The gene expression analysis of the members in this conserved ortholog group in the model ciliate, Tetrahymena thermophila, suggested that these GPCRs play important roles in the life cycle of ciliates. In summary, this study provides the first comprehensive genome-wide identification of GPCRs in ciliates, improving our understanding of the evolution and function of GPCR in ciliates.

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.

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.

Cooperative hydrodynamics accompany multicellular-like colonial organization in the unicellular ciliate Stentor

Evolution of multicellularity from early unicellular ancestors is arguably one of the most important transitions since the origin of life1,2. Multicellularity is often associated with higher nutrient uptake3, better defense against predation, cell specialization and better division of labor4. While many single-celled organisms exhibit both solitary and colonial existence3,5,6, the organizing principles governing the transition and the benefits endowed are less clear. Using the suspension-feeding unicellular protist Stentor coeruleus, we show that hydrodynamic coupling between proximal neighbors results in faster feeding flows that depend on the separation between individuals. Moreover, we find that the accrued benefits in feeding current enhancement are typically asymmetric- individuals with slower solitary currents gain more from partnering than those with faster currents. We find that colony-formation is ephemeral in Stentor and individuals in colonies are highly dynamic unlike other colony-forming organisms like Volvox carteri 3. Our results demonstrate benefits endowed by the colonial organization in a simple unicellular organism and can potentially provide fundamental insights into the selective forces favoring early evolution of multicellular organization.

Studying Habituation in Stentor coeruleus

Learning is usually associated with a complex nervous system, but there is increasing evidence that life at all levels, down to single cells, can display intelligent behaviors. In both natural and artificial systems, learning is the adaptive updating of system parameters based on new information, and intelligence is a measure of the computational process that facilitates learning. Stentor coeruleus is a unicellular pond-dwelling organism that exhibits habituation, a form of learning in which a behavioral response decreases following a repeated stimulus. Stentor contracts in response to mechanical stimulation, which is an apparent escape response from aquatic predators. However, repeated low-force perturbations induce habituation, demonstrated by a progressive reduction in contraction probability. Here, we introduce a method for quantifying Stentor habituation using a microcontroller board-linked apparatus that can deliver mechanical pulses at a specified force and frequency, including methods for building the apparatus and setting up the experiment in a way that minimizes external perturbations. In contrast to the previously described approaches for mechanically stimulating Stentor, this device allows the force of stimulation to be varied under computer control during the course of a single experiment, thus greatly increasing the variety of input sequences that can be applied. Understanding habituation at the level of a single cell will help characterize learning paradigms that are independent of complex circuitry.

Biodiversity of freshwater ciliates (Protista, Ciliophora) in the Lake Weishan Wetland, China: the state of the art

Ciliates are core components of the structure of and function of aquatic microbial food webs. They play an essential role in the energy flow and material circulation within aquatic ecosystems. However, studies on the taxonomy and biodiversity of freshwater ciliates, especially those in wetlands in China are limited. To address this issue, a project to investigate the freshwater ciliates of the Lake Weishan Wetland, Shandong Province, commenced in 2019. Here, we summarize our findings to date on the diversity of ciliates. A total of 187 ciliate species have been found, 94 of which are identified to species-level, 87 to genus-level, and six to family-level. These species show a high morphological diversity and represent five classes, i.e., Heterotrichea, Litostomatea, Prostomatea, Oligohymenophorea, and Spirotrichea. The largest number of species documented are oligohymenophoreans. A comprehensive database of these ciliates, including morphological data, gene sequences, microscope slide specimens and a DNA bank, has been established. In the present study, we provide an annotated checklist of retrieved ciliates as well as information on the sequences of published species. Most of these species are recorded in China for the first time and more than 20% are tentatively identified as new to science. Additionally, an investigation of environmental DNA revealed that the ciliate species diversity in Lake Weishan Wetland is higher than previously supposed.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-022-00154-x.

Insights into the Mechanism of Pre-mRNA Splicing of Tiny Introns from the Genome of a Giant Ciliate Stentor coeruleus

Stentor coeruleus is a ciliate known for its regenerative ability. Recent genome sequencing reveals that its spliceosomal introns are exceptionally small. We wondered whether the multimegadalton spliceosome has any unique characteristics for removal of the tiny introns. First, we analyzed intron features and identified spliceosomal RNA/protein components. We found that all snRNAs are present, whereas many proteins are conserved but slightly reduced in size. Some regulators, such as Serine/Arginine-rich proteins, are noticeably undetected. Interestingly, while most parts of spliceosomal proteins, including Prp8′s positively charged catalytic cavity, are conserved, regions of branching factors projecting to the active site are not. We conjecture that steric-clash avoidance between spliceosomal proteins and a sharply looped lariat might occur, and splicing regulation may differ from other species.

Hydrodynamic dissection of Stentor coeruleus in a microfluidic cross junction

Stentor coeruleus, a single-cell ciliated protozoan, is a model organism for wound healing and regeneration studies. Despite Stentor's large size (up to 2 mm in extended state), microdissection of Stentor remains challenging. In this work, we describe a hydrodynamic cell splitter, consisting of a microfluidic cross junction, capable of splitting Stentor cells in a non-contact manner at a high throughput of ∼500 cells per minute under continuous operation. Introduction of asymmetry in the flow field at the cross junction leads to asymmetric splitting of the cells to generate cell fragments as small as ∼8.5 times the original cell size. Characterization of cell fragment viability shows reduced 5-day survival as fragment size decreases and as the extent of hydrodynamic stress imposed on the fragments increases. Our results suggest that cell fragment size and composition, as well as mechanical stress, play important roles in the long-term repair of Stentor cells and warrant further investigations. Nevertheless, the hydrodynamic splitter can be useful for studying phenomena immediately after cell splitting, such as the closure of wounds in the plasma membrane which occurs on the order of 100-1000 seconds in Stentor.

Modular, cascade-like transcriptional program of regeneration in Stentor

The giant ciliate Stentor coeruleus is a classical model system for studying regeneration and morphogenesis in a single cell. The anterior of the cell is marked by an array of cilia, known as the oral apparatus, which can be induced to shed and regenerate in a series of reproducible morphological steps, previously shown to require transcription. If a cell is cut in half, each half regenerates an intact cell. We used RNA sequencing (RNAseq) to assay the dynamic changes in Stentor’s transcriptome during regeneration, after both oral apparatus shedding and bisection, allowing us to identify distinct temporal waves of gene expression including kinases, RNA -binding proteins, centriole biogenesis factors, and orthologs of human ciliopathy genes. By comparing transcriptional profiles of different regeneration events, we identified distinct modules of gene expression corresponding to oral apparatus regeneration, posterior holdfast regeneration, and recovery after wounding. By measuring gene expression after blocking translation, we show that the sequential waves of gene expression involve a cascade mechanism in which later waves of expression are triggered by translation products of early-expressed genes. Among the early-expressed genes, we identified an E2F transcription factor and the RNA-binding protein Pumilio as potential regulators of regeneration based on the expression pattern of their predicted target genes. RNAi-mediated knockdown experiments indicate that Pumilio is required for regenerating oral structures of the correct size. E2F is involved in the completion of regeneration but is dispensable for earlier steps. This work allows us to classify regeneration genes into groups based on their potential role for regeneration in distinct cell regeneration paradigms, and provides insight into how a single cell can coordinate complex morphogenetic pathways to regenerate missing structures.

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.

Determining protein polarization proteome-wide using physical dissection of individual Stentor coeruleus cells

Cellular components are non-randomly arranged with respect to the shape and polarity of the whole cell.^1-4 Patterning within cells can extend down to the level of individual proteins and mRNA.^5,6 But how much of the proteome is actually localized with respect to cell polarity axes? Proteomics combined with cellular fractionation^7-11 has shown that most proteins localize to one or more organelles but does not tell us how many proteins have a polarized localization with respect to the large-scale polarity axes of the intact cell. Genome-wide localization studies in yeast^12-15 found that only a few percent of proteins have a localized position relative to the cell polarity axis defined by sites of polarized cell growth. Here, we describe an approach for analyzing protein distribution within a cell with a visibly obvious global patterning-the giant ciliate Stentor coeruleus.^16,17 Ciliates, including Stentor, have highly polarized cell shapes with visible surface patterning.^1,18 A Stentor cell is roughly 2 mm long, allowing a "proteomic dissection" in which microsurgery is used to separate cellular fragments along the anterior-posterior axis, followed by comparative proteomic analysis. In our analysis, 25% of the proteome, including signaling proteins, centrin/SFI proteins, and GAS2 orthologs, shows a polarized location along the cell's anterior-posterior axis. We conclude that a large proportion of all proteins are polarized with respect to global cell polarity axes and that proteomic dissection provides a simple and effective approach for spatial proteomics.

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.

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.

The Codon Usage Bias Analysis of Free-Living Ciliates' Macronuclear Genomes and Clustered Regularly Interspaced Short Palindromic Repeats/Cas9 Vector Construction of Stylonychia lemnae

Ciliates represent higher unicellular animals, and several species are also important model organisms for molecular biology research. Analyses of codon usage bias (CUB) of the macronuclear (MAC) genome in ciliates can not only promote a better understanding of the genetic mode and evolution history of these organisms but also help optimize codons to improve the gene editing efficiency of model ciliates. In this study, macronuclear genome sequences of nine free-living ciliates were analyzed with CodonW software to calculate the following indices: the guanine-cytosine content (GC); the frequency of the nucleotides U, C, A, and G at the third position of codons (U3s, C3s, A3s, G3s); the effective number of codons (ENC); the correlation between GC at the first and second positions (GC12); the frequency of the nucleotides G + C at the third position of synonymous codons (GC3s); the relative synonymous codon usage (RSCU). Parity rule 2 plot analysis, neutrality plot analysis, and correlation analysis were performed to explore the factors that influence codon preference. The results showed that the GC contents in nine ciliates' MAC genomes were lower than 50% and appeared AT-rich. The base compositions of GC12 and GC3s are markedly distinct and the codon usage pattern and evolution of ciliates are affected by genetic mutation and natural selection. According to the synonymous codon analysis, the codons of most ciliates ended with A or U and eight codons were the general optimal codons of nine ciliates. A clustered regularly interspaced short palindromic repeats/Cas9 (CRISPR/Cas9) expression vector of Stylonychia lemnae was constructed by optimizing the macronuclear genome codon and was successfully used to knock out the Adss gene. This is the first such extensive investigation of the MAC genome CUB of ciliates and the initial successful application of the CRISPR/Cas9 technique in free-living ciliates.

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.

Regeneration in Stentor coeruleus

We often think about regeneration in terms of replacing missing structures, such as organs or tissues, with new structures generated via cell proliferation and differentiation. But at a smaller scale, single cells, themselves, are capable of regenerating when part of the cell has been removed. A classic model organism that facilitates the study of cellular regeneration in the giant ciliate Stentor coeruleus. These cells, which can grow to more than a millimeter in size, have the ability to survive after extensive wounding of their surface, and are able to regenerate missing structures. Even a small piece of a cell can regenerate a whole cell with normal geometry, in a matter of hours. Such regeneration requires cells to be able to trigger organelle biogenesis in response to loss of structures. But subcellular regeneration also relies on intracellular mechanisms to create and maintain global patterning within the cell. These mechanisms are not understood, but at a conceptual level they involve processes that resemble those seen in animal development and regeneration. Here we discuss single-celled regeneration in Stentor from the viewpoint of standard regeneration paradigms in animals. For example, there is evidence that regeneration of the oral apparatus in Stentor follows a sender-receiver model similar to crustacean eyestalk regeneration. By drawing these analogies, we find that many of the concepts already known from the study of animal-scale regeneration and development can be applied to the study of regeneration at the cellular level, such as the concepts of determination, induction, mosaic vs. regulative development, and epimorphosis vs. morphallaxis. We propose that the similarities may go beyond analogy, and that some aspects of animal development and regeneration may have evolved by exploiting pre-existing subcellular developmental strategies from unicellular ancestors.

The Cilioprotist Cytoskeleton , a Model for Understanding How Cell Architecture and Pattern Are Specified: Recent Discoveries from Ciliates and Comparable Model Systems

The cytoskeletons of eukaryotic, cilioprotist microorganisms are complex, highly patterned, and diverse, reflecting the varied and elaborate swimming, feeding, reproductive, and sensory behaviors of the multitude of cilioprotist species that inhabit the aquatic environment. In the past 10-20 years, many new discoveries and technologies have helped to advance our understanding of how cytoskeletal organelles are assembled in many different eukaryotic model systems, in relation to the construction and modification of overall cellular architecture and function. Microtubule organizing centers, particularly basal bodies and centrioles, have continued to reveal their central roles in architectural engineering of the eukaryotic cell, including in the cilioprotists. This review calls attention to (1) published resources that illuminate what is known of the cilioprotist cytoskeleton; (2) recent studies on cilioprotists and other model organisms that raise specific questions regarding whether basal body- and centriole-associated nucleic acids, both DNA and RNA, should continue to be considered when seeking to employ cilioprotists as model systems for cytoskeletal research; and (3) new, mainly imaging, technologies that have already proven useful for, but also promise to enhance, future cytoskeletal research on cilioprotists.

Evolution of plant telomerase RNAs: farther to the past, deeper to the roots

The enormous sequence heterogeneity of telomerase RNA (TR) subunits has thus far complicated their characterization in a wider phylogenetic range. Our recent finding that land plant TRs are, similarly to known ciliate TRs, transcribed by RNA polymerase III and under the control of the type-3 promoter, allowed us to design a novel strategy to characterize TRs in early diverging Viridiplantae taxa, as well as in ciliates and other Diaphoretickes lineages. Starting with the characterization of the upstream sequence element of the type 3 promoter that is conserved in a number of small nuclear RNAs, and the expected minimum TR template region as search features, we identified candidate TRs in selected Diaphoretickes genomes. Homologous TRs were then used to build covariance models to identify TRs in more distant species. Transcripts of the identified TRs were confirmed by transcriptomic data, RT-PCR and Northern hybridization. A templating role for one of our candidates was validated in Physcomitrium patens. Analysis of secondary structure demonstrated a deep conservation of motifs (pseudoknot and template boundary element) observed in all published TRs. These results elucidate the evolution of the earliest eukaryotic TRs, linking the common origin of TRs across Diaphoretickes, and underlying evolutionary transitions in telomere repeats.

Ciliate behavior: blueprints for dynamic cell biology and microscale robotics

Place a drop of pond water under the microscope, and you will likely find an ocean of extraordinary and diverse single-celled organisms called ciliates. This remarkable group of single-celled organisms wield microtubules, active systems, electrical signaling, and chemical sensors to build intricate geometrical structures and perform complex behaviors that can appear indistinguishable from those of macroscopic animals. Advances in computer vision and machine learning are making it possible to completely digitize and track the dynamics of complex ciliates and mine these data for the hidden structure, patterns, and motifs that are responsible for their behaviors. By deconstructing the diversity of ciliate behaviors in the natural world, themes for organizing and controlling matter at the microscale are beginning to take hold, suggesting new modular approaches for the design of autonomous molecular machines that emulate nature’s finest examples.

Microfluidic guillotine reveals multiple timescales and mechanical modes of wound response in Stentor coeruleus

Background

Wound healing is one of the defining features of life and is seen not only in tissues but also within individual cells. Understanding wound response at the single-cell level is critical for determining fundamental cellular functions needed for cell repair and survival. This understanding could also enable the engineering of single-cell wound repair strategies in emerging synthetic cell research. One approach is to examine and adapt self-repair mechanisms from a living system that already demonstrates robust capacity to heal from large wounds. Towards this end, Stentor coeruleus, a single-celled free-living ciliate protozoan, is a unique model because of its robust wound healing capacity. This capacity allows one to perturb the wounding conditions and measure their effect on the repair process without immediately causing cell death, thereby providing a robust platform for probing the self-repair mechanism.

Results

Here we used a microfluidic guillotine and a fluorescence-based assay to probe the timescales of wound repair and of mechanical modes of wound response in Stentor. We found that Stentor requires ~ 100–1000 s to close bisection wounds, depending on the severity of the wound. This corresponds to a healing rate of ~ 8–80 μm²/s, faster than most other single cells reported in the literature. Further, we characterized three distinct mechanical modes of wound response in Stentor: contraction, cytoplasm retrieval, and twisting/pulling. Using chemical perturbations, active cilia were found to be important for only the twisting/pulling mode. Contraction of myonemes, a major contractile fiber in Stentor, was surprisingly not important for the contraction mode and was of low importance for the others.

Conclusions

While events local to the wound site have been the focus of many single-cell wound repair studies, our results suggest that large-scale mechanical behaviors may be of greater importance to single-cell wound repair than previously thought. The work here advances our understanding of the wound response in Stentor and will lay the foundation for further investigations into the underlying components and molecular mechanisms involved.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-00970-0.

Genome Size Covaries More Positively with Propagule Size than Adult Size: New Insights into an Old Problem

The body size and (or) complexity of organisms is not uniformly related to the amount of genetic material (DNA) contained in each of their cell nuclei ('genome size'). This surprising mismatch between the physical structure of organisms and their underlying genetic information appears to relate to variable accumulation of repetitive DNA sequences, but why this variation has evolved is little understood. Here, I show that genome size correlates more positively with egg size than adult size in crustaceans. I explain this and comparable patterns observed in other kinds of animals and plants as resulting from genome size relating strongly to cell size in most organisms, which should also apply to single-celled eggs and other reproductive propagules with relatively few cells that are pivotal first steps in their lives. However, since body size results from growth in cell size or number or both, it relates to genome size in diverse ways. Relationships between genome size and body size should be especially weak in large organisms whose size relates more to cell multiplication than to cell enlargement, as is generally observed. The ubiquitous single-cell 'bottleneck' of life cycles may affect both genome size and composition, and via both informational (genotypic) and non-informational (nucleotypic) effects, many other properties of multicellular organisms (e.g., rates of growth and metabolism) that have both theoretical and practical significance.

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.

Characterization of Pre-mRNA Splicing and Spliceosomal Machinery in Porphyridium purpureum and Evolutionary Implications for Red Algae

Pre-mRNA splicing is a highly conserved eukaryotic process, but our understanding of it is limited by a historical focus on well-studied organisms such as humans and yeast. There is considerable diversity in mechanisms and components of pre-mRNA splicing, especially in lineages that have evolved under the pressures of genome reduction. The ancestor of red algae is thought to have undergone genome reduction prior to the lineage's radiation, resulting in overall gene and intron loss in extant groups. Previous studies on the extremophilic red alga Cyanidioschyzon merolae revealed an intron-sparse genome with a highly reduced spliceosome. To determine whether these features applied to other red algae, we investigated multiple aspects of pre-mRNA splicing in the mesophilic red alga Porphyridium purpureum. Through strand-specific RNA-Seq, we observed high levels of intron retention across a large number of its introns, and nearly half of the transcripts for these genes are not spliced at all. We also discovered a relationship between variability of 5' splice site sequences and levels of splicing. To further investigate the connections between intron retention and splicing machinery, we bioinformatically assembled the P. purpureum spliceosome, and biochemically verified the presence of snRNAs. While most other core spliceosomal components are present, our results suggest highly divergent or missing U1 snRNP proteins, despite the presence of an uncharacteristically long U1 snRNA. These unusual aspects highlight the diverse nature of pre-mRNA splicing that can be seen in lesser-studied eukaryotes, raising the importance of investigating fundamental eukaryotic processes outside of model organisms.

Characterization and Comparative Analyses of Mitochondrial Genomes in Single-Celled Eukaryotes to Shed Light on the Diversity and Evolution of Linear Molecular Architecture

Determination and comparisons of complete mitochondrial genomes (mitogenomes) are important to understand the origin and evolution of mitochondria. Mitogenomes of unicellular protists are particularly informative in this regard because they are gene-rich and display high structural diversity. Ciliates are a highly diverse assemblage of protists and their mitogenomes (linear structure with high A+T content in general) were amongst the first from protists to be characterized and have provided important insights into mitogenome evolution. Here, we report novel mitogenome sequences from three representatives (Strombidium sp., Strombidium cf. sulcatum, and Halteria grandinella) in two dominant ciliate lineages. Comparative and phylogenetic analyses of newly sequenced and previously published ciliate mitogenomes were performed and revealed a number of important insights. We found that the mitogenomes of these three species are linear molecules capped with telomeric repeats that differ greatly among known species. The genomes studied here are highly syntenic, but larger in size and more gene-rich than those of other groups. They also all share an AT-rich tandem repeat region which may serve as the replication origin and modulate initiation of bidirectional transcription. More generally we identified a split version of ccmf, a cytochrome c maturation-related gene that might be a derived character uniting taxa in the subclasses Hypotrichia and Euplotia. Finally, our mitogenome comparisons and phylogenetic analyses support to reclassify Halteria grandinella from the subclass Oligotrichia to the subclass Hypotrichia. These results add to the growing literature on the unique features of ciliate mitogenomes, shedding light on the diversity and evolution of their linear molecular architecture.

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.

RNA Interference by Cyanobacterial Feeding Demonstrates the SCSG1 Gene Is Essential for Ciliogenesis during Oral Apparatus Regeneration in Stentor

In the giant ciliate Stentor coeruleus, oral apparatus (OA) regeneration is an experimentally tractable regeneration paradigm that occurs via a series of morphological steps. OA regeneration is thought to be driven by a complex regulatory system that orchestrates the temporal expression of conserved and specific genes. We previously identified a S. coeruleus-specific gene (named SCSG1) that was significantly upregulated during the ciliogenesis stages of OA regeneration, with an expression peak at the stage of the first OA cilia appearance. We established a novel RNA interference (RNAi) method through cyanobacteria Synechocystis sp. PCC6803 feeding in S. coeruleus. The expression of SCSG1 gene was significantly knocked down by using this method and induced abnormal ciliogenesis of OA regeneration in S. coeruleus, suggesting that SCSG1 is essential for OA regeneration in S. coeruleus. This novel RNAi method by cyanobacterial feeding has potential utility for studying other ciliates.

SIGAR: Inferring Features of Genome Architecture and DNA Rearrangements by Split-Read Mapping

Ciliates are microbial eukaryotes with distinct somatic and germline genomes. Postzygotic development involves extensive remodeling of the germline genome to form somatic chromosomes. Ciliates therefore offer a valuable model for studying the architecture and evolution of programed genome rearrangements. Current studies usually focus on a few model species, where rearrangement features are annotated by aligning reference germline and somatic genomes. Although many high-quality somatic genomes have been assembled, a high-quality germline genome assembly is difficult to obtain due to its smaller DNA content and abundance of repetitive sequences. To overcome these hurdles, we propose a new pipeline, SIGAR (Split-read Inference of Genome Architecture and Rearrangements) to infer germline genome architecture and rearrangement features without a germline genome assembly, requiring only short DNA sequencing reads. As a proof of principle, 93% of rearrangement junctions identified by SIGAR in the ciliate Oxytricha trifallax were validated by the existing germline assembly. We then applied SIGAR to six diverse ciliate species without germline genome assemblies, including Ichthyophthirius multifilii, a fish pathogen. Despite the high level of somatic DNA contamination in each sample, SIGAR successfully inferred rearrangement junctions, short eliminated sequences, and potential scrambled genes in each species. This pipeline enables pilot surveys or exploration of DNA rearrangements in species with limited DNA material access, thereby providing new insights into the evolution of chromosome rearrangements.

Ciliates as model organisms for the ecotoxicological risk assessment of heavy metals: A meta-analysis

Ciliates are key components of aquatic ecosystems, significantly contributing to the decomposition of organic matter and energy transfer to higher trophic levels. They are considered good biological indicators of chemical pollution and relatively sensitive to heavy metal contamination. In this study, we performed a meta-analysis of the available toxicity data of heavy metals and ciliates to assess: (1) the sensitivity of freshwater ciliates to different heavy metals, (2) the relative sensitivity of ciliates in comparison to the standard test species used in ecotoxicological risk assessment, and (3) the difference in sensitivity across ciliate taxa. Our study shows that the tolerance of ciliates to heavy metals varies notably, which is partly influenced by differences in methodological conditions across studies. Ciliates are, in general, sensitive to Mercury > Cadmium > Copper > Zinc > Lead > Chromium. Also, this study shows that most ciliates are more tolerant to heavy metal pollution than the standard test species used in ecotoxicological risk assessments, i.e., Raphidocelis subcapitata, Daphnia magna, and Onchornyncus mykiss. Threshold concentrations derived from toxicity data for these species is expected to confer sufficient protection for the vast majority of ciliate species. Our data analysis also shows that the most commonly tested ciliate species, Paramecium caudatum and Tetrahymena thermophila, are not necessarily the most sensitive ones to heavy metal pollution. Finally, this study stresses the importance of developing standard toxicity test protocols for ciliates, which could lead to a better comprehension of the toxicological impact of heavy metals and other contaminants to ciliate species.

Proteomic identification and expression of oral apparatus constituents in cell regeneration of giant ciliate Stentor coeruleus (strain WHEL)

The giant ciliate Stentor coeruleus (S. coeruleus) is a suitable model organism for studying morphogenesis and regeneration at the single-cell level. It contains a prominent structure on the anterior end of the cell, known as the oral apparatus (OA). OA can be induced to shed by urea treatment and then new OA regenerates via a series of defined morphological events and the cell resumes normal feeding activity. We identified OA constituents in S. coeruleus by mass spectrometry. A total of 882 OA-associated proteins were identified; the homologs of 181 of these are known OA constituents in other organisms. The expression pattern of OA-associated genes during regeneration was investigated using single-cell transcriptome sequencing. The expression of most OA-associated genes was high during regeneration, indicating their stable expression after OA shedding. We also identified OA-associated differentially expressed genes that may be involved in regulating OA reconstruction. In summary, this study gives preliminary insight into the molecular basis of OA in S. coeruleus.

Single-Cell Transcriptomics Reveal a Correlation between Genome Architecture and Gene Family Evolution in Ciliates

Ciliates, a eukaryotic clade that is over 1 billion years old, are defined by division of genome function between transcriptionally inactive germline micronuclei and functional somatic macronuclei. To date, most analyses of gene family evolution have been limited to cultivable model lineages (e.g., Tetrahymena, Paramecium, Oxytricha, and Stylonychia). Here, we focus on the uncultivable Karyorelictea and its understudied sister class Heterotrichea, which represent two extremes in genome architecture. Somatic macronuclei within the Karyorelictea are described as nearly diploid, while the Heterotrichea have hyperpolyploid somatic genomes. Previous analyses indicate that genome architecture impacts ciliate gene family evolution as the most diverse and largest gene families are found in lineages with extensively processed somatic genomes (i.e., possessing thousands of gene-sized chromosomes). To further assess ciliate gene family evolution, we analyzed 43 single-cell transcriptomes from 33 ciliate species representing 10 classes. Focusing on conserved eukaryotic genes, we use estimates of transcript diversity as a proxy for the number of paralogs in gene families among four focal clades: Karyorelictea, Heterotrichea, extensive fragmenters (with gene-size somatic chromosomes), and non-extensive fragmenters (with more traditional somatic chromosomes), the latter two within the subphylum Intramacronucleata. Our results show that (i) the Karyorelictea have the lowest average transcript diversity, while Heterotrichea are highest among the four groups; (ii) proteins in Karyorelictea are under the highest functional constraints, and the patterns of selection in ciliates may reflect genome architecture; and (iii) stop codon reassignments vary among members of the Heterotrichea and Spirotrichea but are conserved in other classes.IMPORTANCE To further our understanding of genome evolution in eukaryotes, we assess the relationship between patterns of molecular evolution within gene families and variable genome structures found among ciliates. We combine single-cell transcriptomics with bioinformatic tools, focusing on understudied and uncultivable lineages selected from across the ciliate tree of life. Our analyses show that genome architecture correlates with patterns of protein evolution as lineages with more canonical somatic genomes, such as the class Karyorelictea, have more conserved patterns of molecular evolution compared to other classes. This study showcases the power of single-cell transcriptomics for investigating genome architecture and evolution in uncultivable microbial lineages and provides transcriptomic resources for further research on genome evolution.

Comparative Genomics Analysis of Ciliates Provides Insights on the Evolutionary History Within "Nassophorea-Synhymenia-Phyllopharyngea" Assemblage

Ciliated protists (ciliates) are widely used for investigating evolution, mostly due to their successful radiation after their early evolutionary branching. In this study, we employed high-throughput sequencing technology to reveal the phylogenetic position of Synhymenia, as well as two classes Nassophorea and Phyllopharyngea, which have been a long-standing puzzle in the field of ciliate systematics and evolution. We obtained genomic and transcriptomic data from single cells of one synhymenian (Chilodontopsis depressa) and six other species of phyllopharyngeans (Chilodochona sp., Dysteria derouxi, Hartmannula sinica, Trithigmostoma cucullulus, Trochilia petrani, and Trochilia sp.). Phylogenomic analysis based on 157 orthologous genes comprising 173,835 amino acid residues revealed the affiliation of C. depressa within the class Phyllopharyngea, and the monophyly of Nassophorea, which strongly support the assignment of Synhymenia as a subclass within the class Phyllopharyngea. Comparative genomic analyses further revealed that C. depressa shares more orthologous genes with the class Nassophorea than with Phyllopharyngea, and the stop codon usage in C. depressa resembles that of Phyllopharyngea. Functional enrichment analysis demonstrated that biological pathways in C. depressa are more similar to Phyllopharyngea than Nassophorea. These results suggest that genomic and transcriptomic data can be used to provide insights into the evolutionary relationships within the "Nassophorea-Synhymenia-Phyllopharyngea" assemblage.

A Complex Hierarchy of Avoidance Behaviors in a Single-Cell Eukaryote

Complex behavior is associated with animals with nervous systems, but decision-making and learning also occur in non-neural organisms [1], including singly nucleated cells [2-5] and multi-nucleate synctia [6-8]. Ciliates are single-cell eukaryotes, widely dispersed in aquatic habitats [9], with an extensive behavioral repertoire [10-13]. In 1906, Herbert Spencer Jennings [14, 15] described in the sessile ciliate Stentor roeseli a hierarchy of responses to repeated stimulation, which are among the most complex behaviors reported for a singly nucleated cell [16, 17]. These results attracted widespread interest [18, 19] and exert continuing fascination [7, 20-22] but were discredited during the behaviorist orthodoxy by claims of non-reproducibility [23]. These claims were based on experiments with the motile ciliate Stentor coeruleus. We acquired and maintained the correct organism in laboratory culture and used micromanipulation and video microscopy to confirm Jennings' observations. Despite significant individual variation, not addressed by Jennings, S. roeseli exhibits avoidance behaviors in a characteristic hierarchy of bending, ciliary alteration, contractions, and detachment, which is distinct from habituation or conditioning. Remarkably, the choice of contraction versus detachment is consistent with a fair coin toss. Such behavioral complexity may have had an evolutionary advantage in protist ecosystems, and the ciliate cortex may have provided mechanisms for implementing such behavior prior to the emergence of multicellularity. Our work resurrects Jennings' pioneering insights and adds to the list of exceptional features, including regeneration [24], genome rearrangement [25], codon reassignment [26], and cortical inheritance [27], for which the ciliate clade is renowned.

Selection for tandem stop codons in ciliate species with reassigned stop codons

The failure of mRNA translation machinery to recognize a stop codon as a termination signal and subsequent translation of the 3' untranslated region (UTR) is referred to as stop codon readthrough, the frequency of which is related to the length, composition, and structure of mRNA sequences downstream of end-of-gene stop codons. Secondary in-frame stop codons within a few positions downstream of the primary stop codons, so-called tandem stop codons (TSCs), serve as backup termination signals, which limit the effects of readthrough: polypeptide product degradation, mislocalization, and aggregation. In this study, ciliate species with UAA and UAG stop codons reassigned to code for glutamine are found to possess statistical excesses of TSCs at the beginning of their 3' UTRs. The overrepresentation of TSCs in these species is greater than that observed in standard code organisms. Though the overall numbers of TSCs are lower in most species with alternative stop codons because they use fewer than three unique stop codons, the relatively great overrepresentation of TSCs in alternative-code ciliate species suggests that there exist stronger selective pressures to maintain TSCs in these organisms compared to standard code organisms.

Population Genetics of Paramecium Mitochondrial Genomes: Recombination, Mutation Spectrum, and Efficacy of Selection

The evolution of mitochondrial genomes and their population-genetic environment among unicellular eukaryotes are understudied. Ciliate mitochondrial genomes exhibit a unique combination of characteristics, including a linear organization and the presence of multiple genes with no known function or detectable homologs in other eukaryotes. Here we study the variation of ciliate mitochondrial genomes both within and across 13 highly diverged Paramecium species, including multiple species from the P. aurelia species complex, with four outgroup species: P. caudatum, P. multimicronucleatum, and two strains that may represent novel related species. We observe extraordinary conservation of gene order and protein-coding content in Paramecium mitochondria across species. In contrast, significant differences are observed in tRNA content and copy number, which is highly conserved in species belonging to the P. aurelia complex but variable among and even within the other Paramecium species. There is an increase in GC content from ∼20% to ∼40% on the branch leading to the P. aurelia complex. Patterns of polymorphism in population-genomic data and mutation-accumulation experiments suggest that the increase in GC content is primarily due to changes in the mutation spectra in the P. aurelia species. Finally, we find no evidence of recombination in Paramecium mitochondria and find that the mitochondrial genome appears to experience either similar or stronger efficacy of purifying selection than the nucleus.

UAA and UAG may Encode Amino Acid in Cathepsin B Gene of Euplotes octocarinatus

The ciliate Euplotes deviates from the universal genetic code by translating UGA as cysteine and using UAA and UAG as the termination codon. Here, we cloned and sequenced the Cathepsin B gene of Euplotes octocarinatus (Eo-CTSB) which containing several in-frame stop codons throughout the coding sequence. We provide evidences, based on 3'-RACE method and Western blot, that the Eo-CTSB gene is actively expressed. Comparison of the derived amino acid sequence with the homologs in other eukaryotes revealed that UAA and UAG may code for glutamine in Eo-CTSB. These findings imply an evolutionary complexity of stop codon reassignment in eukaryotes.

Tetrahymena Comparative Genomics Database (TCGD): a community resource for Tetrahymena

Ciliates are a large and diverse group of unicellular organisms characterized by having the following two distinct type of nuclei within a single cell: micronucleus (MIC) and macronucleus (MAC). Although the genomes of several ciliates in different groups have been sequenced, comparative genomics data for multiple species within a ciliate genus are not yet available. Here we collected the genome information and comparative genomics analysis results for 10 species in the Tetrahymena genus, including the previously sequenced model organism Tetrahymena thermophila and 9 newly sequenced species, and constructed a genus-level comparative analysis platform, the Tetrahymena Comparative Genomics Database (TCGD). Genome sequences, transcriptomic data, gene models, functional annotation, ortholog groups and synteny maps were built into this database and a user-friendly interface was developed for searching, visualizing and analyzing these data. In summary, the TCGD (http://ciliate.ihb.ac.cn) will be an important and useful resource for the ciliate research community.

Genome analyses of the new model protist Euplotes vannus focusing on genome rearrangement and resistance to environmental stressors

As a model organism for studies of cell and environmental biology, the free-living and cosmopolitan ciliate Euplotes vannus shows intriguing features like dual genome architecture (i.e., separate germline and somatic nuclei in each cell/organism), "gene-sized" chromosomes, stop codon reassignment, programmed ribosomal frameshifting (PRF) and strong resistance to environmental stressors. However, the molecular mechanisms that account for these remarkable traits remain largely unknown. Here we report a combined analysis of de novo assembled high-quality macronuclear (MAC; i.e., somatic) and partial micronuclear (MIC; i.e., germline) genome sequences for E. vannus, and transcriptome profiling data under varying conditions. The results demonstrate that: (a) the MAC genome contains more than 25,000 complete "gene-sized" nanochromosomes (~85 Mb haploid genome size) with the N50 ~2.7 kb; (b) although there is a high frequency of frameshifting at stop codons UAA and UAG, we did not observe impaired transcript abundance as a result of PRF in this species as has been reported for other euplotids; (c) the sequence motif 5'-TA-3' is conserved at nearly all internally-eliminated sequence (IES) boundaries in the MIC genome, and chromosome breakage sites (CBSs) are duplicated and retained in the MAC genome; (d) by profiling the weighted correlation network of genes in the MAC under different environmental stressors, including nutrient scarcity, extreme temperature, salinity and the presence of ammonia, we identified gene clusters that respond to these external physical or chemical stimulations, and (e) we observed a dramatic increase in HSP70 gene transcription under salinity and chemical stresses but surprisingly, not under temperature changes; we link this temperature-resistance to the evolved loss of temperature stress-sensitive elements in regulatory regions. Together with the genome resources generated in this study, which are available online at Euplotes vannus Genome Database (http://evan.ciliate.org), these data provide molecular evidence for understanding the unique biology of highly adaptable microorganisms.

Spliceosome Profiling Visualizes Operations of a Dynamic RNP at Nucleotide Resolution

Tools to understand how the spliceosome functions in vivo have lagged behind advances in the structural biology of the spliceosome. Here, methods are described to globally profile spliceosome-bound pre-mRNA, intermediates, and spliced mRNA at nucleotide resolution. These tools are applied to three yeast species that span 600 million years of evolution. The sensitivity of the approach enables the detection of canonical and non-canonical events, including interrupted, recursive, and nested splicing. This application of statistical modeling uncovers independent roles for the size and position of the intron and the number of introns per transcript in substrate progression through the two catalytic stages. These include species-specific inputs suggestive of spliceosome-transcriptome coevolution. Further investigations reveal the ATP-dependent discard of numerous endogenous substrates after spliceosome assembly in vivo and connect this discard to intron retention, a form of splicing regulation. Spliceosome profiling is a quantitative, generalizable global technology used to investigate an RNP central to eukaryotic gene expression.