Large-scale phylogenomic analysis reveals the phylogenetic position of the problematic taxon Protocruzia and unravels the deep phylogenetic affinities of the ciliate lineages

Dynamics and timing of diversification events of ciliated eukaryotes from a large phylogenomic perspective

Ciliophora, an exceptionally diverse lineage of unicellular eukaryotes, exhibits a remarkable range of species richness across classes in the ciliate Tree of Life. In this study, we have acquired transcriptome and genome data from 40 representative species in seven ciliate classes. Utilizing 247 genes and 105 taxa, we devised a comprehensive phylogenomic tree for Ciliophora, encompassing over 60 % of orders and constituting the most extensive dataset of ciliate species to date. We established a robust phylogenetic framework that encompasses ambiguous taxa and the major classes within the phylum. Our findings support the monophyly of each of two subphyla (Postciliodesmatophora and Intramacronucleata), along with three subclades (Protocruzia, CONTHREEP, and SAPML) nested within Intramacronucleata, and elucidate evolutionary positions among the major classes within the phylum. Drawing on the robust ciliate Tree of Life and three constraints, we estimated the radiation of Ciliophora around 1175 Ma during the middle of the Proterozoic Eon, and most of the ciliate classes diverged from their sister lineage during the latter half of this period. Additionally, based on the time-calibrated tree and species richness pattern, we investigated net diversification rates of Ciliophora and its classes. The global net diversification rate for Ciliophora was estimated at 0.004979 species/Ma. Heterogeneity in net diversification rates was evident at the class level, with faster rates observed in Oligohymenophorea and Spirotrichea than other classes within the subclades CONTHREEP and SAPML, respectively. Notably, our analysis suggests that variations in net diversification rates, rather than clade ages, appear to contribute to the differences in species richness in Ciliophora at the class level.

Giants' cooperation: a draft genome of the giant ciliate Muniziella cunhai suggests its ecological role in the capybara's digestive metabolism

Several hundred ciliate species live in animals' guts as a part of their microbiome. Among them, Muniziella cunhai (Trichostomatia, Pycnotrichidae), the largest described ciliate, is found exclusively associated with Hydrochoerus hydrochaeris (capybara), the largest known rodent reaching up to 90 kg. Here, we present the sequence, structural and functional annotation of this giant microeukaryote macronuclear genome and discuss its phylogenetic placement. The 85 Mb genome is highly AT rich (GC content 25.71 %) and encodes a total of 11 397 protein-coding genes, of which 2793 could have their functions predicted with automated functional assignments. Functional annotation showed that M. cunhai can digest recalcitrant structural carbohydrates, non-structural carbohydrates, and microbial cell walls, suggesting a role in diet metabolization and in microbial population control in the capybara's intestine. Moreover, the phylogenetic placement of M. cunhai provides insights on the origins of gigantism in the subclass Trichostomatia.

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

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

Rediscovery of Remarkably Rare Anaerobic Tentaculiferous Ciliate Genera Legendrea and Dactylochlamys (Ciliophora: Litostomatea)

Free-living anaerobic ciliates are of considerable interest from an ecological and an evolutionary standpoint. Extraordinary tentacle-bearing predatory lineages have evolved independently several times within the phylum Ciliophora, including two rarely encountered anaerobic litostomatean genera, Legendrea and Dactylochlamys. In this study, we significantly extend the morphological and phylogenetic characterization of these two poorly known groups of predatory ciliates. We provide the first phylogenetic analysis of the monotypic genus Dactylochlamys and the three valid species of Legendrea based on the 18S rRNA gene and ITS-28S rRNA gene sequences. Prior to this study, neither group had been studied using silver impregnation methods. We provide the first protargol-stained material and also a unique video material including documentation, for the first time, of the hunting and feeding behavior of a Legendrea species. We briefly discuss the identity of methanogenic archaeal and bacterial endosymbionts of both genera based on 16S rRNA gene sequences, and the importance of citizen science for ciliatology from a historical and contemporary perspective.

Comprehensive phylogenomic analyses reveal that order Armophorida is most closely related to class Armophorea (Protista, Ciliophora)

Ciliate species within the class Armophorea are widely distributed in various anaerobic environments, hence they are of great interest to researchers studying evolution and adaptation of eukaryotes to extreme habitats. However, phylogenetic relationships within the class remain largely elusive, most especially assignment of the order Armophorida and classification within the family Metopidae. In this study, we newly sequenced transcriptomes and the SSU rDNA of five armophorean species, Sulfonecta cf. uniserialis (order Armophorida), Nyctotheroides sp. (order Clevelandellida), and Metopus major, M. paraes, and Brachonella contorta (order Metopida). Comprehensive phylogenomic analyses revealed that Armophorea was most closely related to Muranotrichea and Parablepharismea. Our results indicate that the order Armophorida either belongs to Armophorea or represents a new class within APM (Armophorea-Parablepharismea-Muranotrichea). Analyses combining ecological niches and molecular trees showed that APM species might descend from an anaerobic free-living ciliate species. Existing molecular phylogenomic/phylogenetic and morphological evidence indicate that the family Metopidae is non-monophyletic and should be further classified with inclusion of additional lines of evidences. Our results provide new insights into the long-debated relationships within Armophorea.

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.

Anaerobic Ciliates as a Model Group for Studying Symbioses in Oxygen-depleted Environments

Anaerobiosis has independently evolved in multiple lineages of ciliates, allowing them to colonize a variety of anoxic and oxygen-depleted habitats. Anaerobic ciliates commonly form symbiotic relationships with various prokaryotes, including methanogenic archaea and members of several bacterial groups. The hypothesized functions of these ecto- and endosymbionts include the symbiont utilizing the ciliate's fermentative end-products to increase host's anaerobic metabolic efficiency, or the symbiont directly providing the host with energy by denitrification or photosynthesis. The host, in turn, may protect the symbiont from competition, the environment, and predation. Despite rapid advances in sampling, molecular, and microscopy methods, as well as the associated broadening of the known diversity of anaerobic ciliates, many aspects of these ciliate symbioses, including host-specificity and co-evolution, remain largely unexplored. Nevertheless, with the number of comparative genomic and transcriptomic analyses targeting anaerobic ciliates and their symbionts on the rise, insights into the nature of these symbioses and the evolution of the ciliate transition to obligate anaerobiosis continue to deepen. This review summarizes the current body of knowledge regarding the complex nature of symbioses in anaerobic ciliates, the diversity of these symbionts, their role in the evolution of ciliate anaerobiosis and their significance in ecosystem-level processes.

Distribution of Membrane Trafficking System Components Across Ciliate Diversity Highlights Heterogenous Organelle-Associated Machinery

The ciliate phylum is a group of protists noted for their unusual membrane trafficking system and apparent environmental ubiquity; as highly successful microbial predators, they are found in all manner of environments and the ability for specific species to adapt to extremely challenging conditions makes them valued as bioindicators. Ciliates have also been used for many years as cell biological models due to their large cell size and ease of culturing, and for many fundamental cell structures, particularly membrane-bound organelles, ciliates were some of the earliest organisms in which these were observed via microscopy. In this study, we carried out a comparative genomic survey of selected membrane trafficking proteins in a pan-ciliate transcriptome and genome dataset. We observed considerable loss of membrane trafficking system (MTS) proteins that would indicate a loss of machinery that is generally conserved across eukaryotic diversity, even after controlling for potentially incomplete genome representation. In particular, the DSL1 complex was missing in all surveyed ciliates. This protein complex has been shown as involved in peroxisome biogenesis in some model systems, and a paucity of DSL1 components has been indicative of degenerate peroxisome. However, Tetrahymena thermophila (formerly Tetrahymena pyroformis) was one of the original models for visualising peroxisomes. Conversely, the AP3 complex essential for mucocyst maturation in T. thermophila, is poorly conserved despite the presence of secretory lysosome-related organelles across ciliate diversity. We discuss potential resolutions for these apparent paradoxes in the context of the heterogenous distribution of MTS machinery across the diversity of ciliates.

A New Protocruzia Species (Ciliophora: Protocruziida) Isolated From the Mariana Trench Area

A new species of Protocruzia, isolated from the deep-sea Pacific Ocean (>3,000-m depth) in the vicinity of the Mariana Trench, is described based on morphological and molecular data. The systematic status of the ciliate genus Protocruzia has long been highly ambiguous, and Protocruzia species have been assigned to an independent class until recently. In the present study, we described Protocruzia marianaensis sp. n. as a small (25–32 × 14–17 μm in vivo) drop-shaped ciliate, with longitudinal furrows along the ciliary rows on the right side, six adoral membranelles, eight somatic kineties, and one macronucleus comprising 7–11 nuclear globules. Phylogenetic analyses inferred from small subunit rRNA gene sequences revealed that seven Protocruzia species in the phylogenetic tree formed a fully supported clade representing an independent class. Protocruzia marianaensis sp. n. was established to be most closely related to Protocruzia adhaerens, with a sequence similarity of 96.64%, and was found to be able to survive at both atmospheric pressure and hydrostatic pressure of 320 bar, thereby indicating effective barotolerance.

Ciliate SSU-rDNA reference alignments and trees for phylogenetic placements of metabarcoding data

Although ciliates are one of the most dominant microbial eukaryotic groups in many environments, there is a lack of updated global ciliate alignments and reference trees that can be used for phylogenetic placement methods to analyze environmental metabarcoding data. Here we fill this gap by providing reference alignments and trees for those ciliates taxa with available SSU-rDNA sequences derived from identified species. Each alignment contains 478 ciliate and six outgroup taxa, and they were made using different masking strategies for alignment positions (unmasked, masked and masked except the hypervariable V4 region). We constrained the monophyly of the major ciliate groups based on the recently updated classification of protists and based on phylogenomic data. Taxa of uncertain phylogenetic position were kept unconstrained, except for Mesodinium species that we constrained to form a clade with the Litostomatea. These ciliate reference alignments and trees can be used to perform taxonomic assignments of metabarcoding data, discover novel ciliate clades, estimate species richness, and overlay measured ecological parameters onto the phylogenetic placements.

Large-scale phylogenomic analysis provides new insights into the phylogeny of the class Oligohymenophorea (Protista, Ciliophora) with establishment of a new subclass Urocentria nov. subcl

The class Oligohymenophorea is one of the most diverse assemblage of ciliated protists, which are particularly important in fundamental biological studies including understanding the evolutionary relationships among the lineages. Phylogenetic relationships within the class remain largely elusive, especially within the subclass Peniculia, which contains the long-standing problematic taxa Urocentrum and Paranassula. In the present study, we sequenced the genomes and/or transcriptomes of six non-culturable oligohymenophoreans using single-cell sequencing techniques. Phylogenomic analysis was performed based on expanded taxon sampling of 85 taxa, including 157 nuclear genes encoding 36,953 amino acids. The results indicate that: (1) urocentrids form an independent branch that is sister to the clade formed by Scuticociliatia and Hymenostomatia, which, together with the morphological data, supports the establishment of a new subclass, Urocentria n. subcl., within Oligohymenophorea; (2) phylogenomic analysis and ortholog comparison reveal a close relationship between Paranassula and peniculines, providing corroborative evidence for removing Paranassula from Nassulida and elevating it as an order, Paranassulida, within the subclass Peniculia; (3) based on the phylogenomic analyses and morphological data, we hypothesize that Peritrichia is the earliest diverging clade within Oligohymenophorea while Scuticociliatia and Hymenostomatia share the most common ancestor, followed successively by Urocentria and Peniculia. In addition, stop codon analyses indicate that oligohymenophoreans widely use UGA as the stop codon, while UAR are reassigned to glutamate (peritrichs) or glutamine (others), supporting the evolutionary hypothesis.

Incorporating mitogenome sequencing into integrative taxonomy: The multidisciplinary redescription of the ciliate Thuricola similis (Peritrichia, Vaginicolidae) provides new insights into the evolutionary relationships among Oligohymenophorea subclasses

The evolutionary relationships among Oligohymenophorea subclasses are under debate as the phylogenomic analysis using a large dataset of nuclear coding genes is significantly different to the 18S rDNA phylogeny, and it is unfortunately not stable within and across different published studies. In addition to nuclear genes, the faster-evolving mitochondrial genes have also shown the ability to solve phylogenetic problems in many ciliated taxa. However, due to the paucity of mitochondrial data, the corresponding work is scarce, let alone the phylogenomic analysis based on mitochondrial gene dataset. In this work, we presented the characterization on Thuricola similis Bock, 1963, a loricate peritrich (Oligohymenophorea), incorporating mitogenome sequencing into integrative taxonomy. As the first mitogenome for the subclass Peritrichia, it is linear, 38,802 bp long, and contains two rRNAs, 12 tRNAs, and 43 open reading frames (ORFs). As a peculiarity, it includes a central repeated region composed of tandemly repeated A-T rich units working as a bi-transcriptional start. Moreover, taking this opportunity, the phylogenomic analyses based on a set of mitochondrial genes were also performed, revealing that T. similis, as a representative of Peritrichia subclass, branches basally to other three Oligohymenophorea subclasses, namely Hymenostomatia, Peniculia, and Scuticociliatia. Evolutionary relationships among those Oligohymenophorea subclasses were discussed, also in the light of recent phylogenomic reconstructions based on a set of nuclear genes. Besides, as a little-known species, T. similis was also redescribed and neotypified based on data from two populations collected from wastewater treatment plants (WWTPs) in Brazil and Italy, by means of integrative methods (i.e., living observation, silver staining methods, scanning and transmission electron microscopy, and 18S rDNA phylogeny). After emended diagnosis, it is characterized by: (1) the sewage habitat; (2) the lorica with a single valve and small undulations; (3) the 7-22 µm-long inner stalk; and (4) the presence of only a single postciliary microtubule on the left side of the aciliferous row in the haplokinety. Among Vaginicolidae family, our 18S rRNA gene-based phylogenetic analysis revealed that Thuricola and Cothurnia are monophyletic genera, and Vaginicola could be a polyphyletic genus.

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.

Assessing the utility of Hsp90 gene for inferring evolutionary relationships within the ciliate subclass Hypotricha (Protista, Ciliophora)

BACKGROUND

Although phylogenomic analyses are increasingly used to reveal evolutionary relationships among ciliates, relatively few nuclear protein-coding gene markers have been tested for their suitability as candidates for inferring phylogenies within this group. In this study, we investigate the utility of the heat-shock protein 90 gene (Hsp90) as a marker for inferring phylogenetic relationships among hypotrich ciliates.

RESULTS

A total of 87 novel Hsp90 gene sequences of 10 hypotrich species were generated. Of these, 85 were distinct sequences. Phylogenetic analyses based on these data showed that: (1) the Hsp90 gene amino acid trees are comparable to the small subunit rDNA tree for recovering phylogenetic relationships at the rank of class, but lack sufficient phylogenetic signal for inferring evolutionary relationships at the genus level; (2) Hsp90 gene paralogs are recent and therefore unlikely to pose a significant problem for recovering hypotrich clades; (3) definitions of some hypotrich orders and families need to be revised as their monophylies are not supported by various gene markers; (4) The order Sporadotrichida is paraphyletic, but the monophyly of the "core" Urostylida is supported; (5) both the subfamily Oxytrichinae and the genus Urosoma seem to be non-monophyletic, but monophyly of Urosoma is not rejected by AU tests.

CONCLUSIONS

Our results for the first time demonstrate that the Hsp90 gene is comparable to SSU rDNA for recovering phylogenetic relationships at the rank of class, and its paralogs are unlikely to pose a significant problem for recovering hypotrich clades. This study shows the value of careful gene marker selection for phylogenomic analyses of ciliates.

Genomics of New Ciliate Lineages Provides Insight into the Evolution of Obligate Anaerobiosis

Oxygen plays a crucial role in energetic metabolism of most eukaryotes. Yet adaptations to low-oxygen concentrations leading to anaerobiosis have independently arisen in many eukaryotic lineages, resulting in a broad spectrum of reduced and modified mitochondrion-related organelles (MROs). In this study, we present the discovery of two new class-level lineages of free-living marine anaerobic ciliates, Muranotrichea, cl. nov. and Parablepharismea, cl. nov., that, together with the class Armophorea, form a major clade of obligate anaerobes (APM ciliates) within the Spirotrichea, Armophorea, and Litostomatea (SAL) group. To deepen our understanding of the evolution of anaerobiosis in ciliates, we predicted the mitochondrial metabolism of cultured representatives from all three classes in the APM clade by using transcriptomic and metagenomic data and performed phylogenomic analyses to assess their evolutionary relationships. The predicted mitochondrial metabolism of representatives from the APM ciliates reveals functional adaptations of metabolic pathways that were present in their last common ancestor and likely led to the successful colonization and diversification of the group in various anoxic environments. Furthermore, we discuss the possible relationship of Parablepharismea to the uncultured deep-sea class Cariacotrichea on the basis of single-gene analyses. Like most anaerobic ciliates, all studied species of the APM clade host symbionts, which we propose to be a significant accelerating factor in the transitions to an obligately anaerobic lifestyle. Our results provide an insight into the evolutionary mechanisms of early transitions to anaerobiosis and shed light on fine-scale adaptations in MROs over a relatively short evolutionary time frame.

Why the -omic future of Apicomplexa should include gregarines

Gregarines, a polyphyletic group of apicomplexan parasites infecting mostly non-vertebrates hosts, remains poorly known at taxonomic, phylogenetic and genomic levels. However, it represents an essential group for understanding evolutionary history and adaptive capacities of apicomplexan parasites to the remarkable diversity of their hosts. Because they have a mostly extracellular lifestyle, gregarines have developed other cellular developmental forms and host-parasite interactions, compared with their much better studied apicomplexan cousins, intracellular parasites of vertebrates (Hemosporidia, Coccidia, Cryptosporidia). This review highlights the promises offered by the molecular exploration of gregarines, that have been until now left on the side of the road of the comparative -omic exploration of apicomplexan parasites. Elucidating molecular bases for both their ultrastructural, functional and behavioural similarities and differences, compared with those of the typical apicomplexan models, is expected to provide entirely novel clues on the adaptive capacities developed by Apicomplexa over evolution. A challenge remains to identify which gregarines should be explored in priority, as recent metadata from open and host-associated environments have confirmed how underestimated is our current view on true gregarine biodiversity. It is now time to turn to gregarines to widen the currently highly skewed view we have of adaptive mechanisms developed by Apicomplexa.

Documentation of a new hypotrich species in the family Amphisiellidae, Lamtostyla gui n. sp. (Protista, Ciliophora) using a multidisciplinary approach

An integrated approach considering both morphologic and molecular data is now required to improve biodiversity estimations and provide more robust systematics interpretations in hypotrichs, a highly differentiated group of ciliates. In present study, we document a new hypotrich species, Lamtostyla gui n. sp., collected from Chongming wetland, Shanghai, China, based on investigations using living observation, protargol staining, scanning and transmission electron microscopy, and gene sequencing. The new species is mainly recognized by having a short amphisiellid median cirral row composed of four cirri, three frontoventral cirri, three dorsal kinetids, four to eight macronuclear nodules, and small colorless cortical granules distributed as rosettes around dorsal bristles. Transmission electron microscope observation finds the associated microtubules of cirri and pharyngeal discs of L. gui are distinct from those in other hypotrichs. Morphogenesis of this species indicates that parental adoral membranelles retained intact or partial renewed is a potential feature to separate Lamtostyla granulifera-group and Lamtostyla lamottei-group. Phylogenetic analysis based on small subunit ribosomal RNA (rRNA) gene shows that this molecular marker is not useful to resolve phylogenetic relationships of the genus Lamtostyla, as well as many other hypotrichous taxa. We additionally characterize the internal transcribed spacers (ITS) region and the almost complete large subunit rRNA, which will be essential for future studies aimed at solving phylogenetic problems of Lamtostyla, or even the family Amphisiellidae. As a final remark, the critical screening of GenBank using ITS genes of our organism allows us to recognize a large amount of hypotrichous sequences have been misclassified as fungi. This observation suggests that hypotrichs could be frequently found in fungi-rich environment and overlooked by fungal specialists.

Comparative Studies on the Polymorphism and Copy Number Variation of mtSSU rDNA in Ciliates (Protista, Ciliophora): Implications for Phylogenetic, Environmental, and Ecological Research

While nuclear small subunit ribosomal DNA (nSSU rDNA) is the most commonly-used gene marker in studying phylogeny, ecology, abundance, and biodiversity of microbial eukaryotes, mitochondrial small subunit ribosomal DNA (mtSSU rDNA) provides an alternative. Recently, both copy number variation and sequence variation of nSSU rDNA have been demonstrated for diverse organisms, which can contribute to misinterpretation of microbiome data. Given this, we explore patterns for mtSSU rDNA among 13 selected ciliates (representing five classes), a major component of microbial eukaryotes, estimating copy number and sequence variation and comparing to that of nSSU rDNA. Our study reveals: (1) mtSSU rDNA copy number variation is substantially lower than that for nSSU rDNA; (2) mtSSU rDNA copy number ranges from 1.0 × 104 to 8.1 × 105; (3) a most common sequence of mtSSU rDNA is also found in each cell; (4) the sequence variation of mtSSU rDNA are mainly indels in poly A/T regions, and only half of species have sequence variation, which is fewer than that for nSSU rDNA; and (5) the polymorphisms between haplotypes of mtSSU rDNA would not influence the phylogenetic topology. Together, these data provide more insights into mtSSU rDNA as a powerful marker especially for microbial ecology studies.

A Phylogenomic Approach to Clarifying the Relationship of Mesodinium within the Ciliophora: A Case Study in the Complexity of Mixed-Species Transcriptome Analyses

Recent high-throughput sequencing endeavors have yielded multigene/protein phylogenies that confidently resolve several inter- and intra-class relationships within the phylum Ciliophora. We leverage the massive sequencing efforts from the Marine Microbial Eukaryote Transcriptome Sequencing Project, other SRA submissions, and available genome data with our own sequencing efforts to determine the phylogenetic position of Mesodinium and to generate the most taxonomically rich phylogenomic ciliate tree to date. Regardless of the data mining strategy, the multiprotein data set, or the molecular models of evolution employed, we consistently recovered the same well-supported relationships among ciliate classes, confirming many of the higher-level relationships previously identified. Mesodinium always formed a monophyletic group with members of the Litostomatea, with mixotrophic species of Mesodinium-M. rubrum, M. major, and M. chamaeleon-being more closely related to each other than to the heterotrophic member, M. pulex. The well-supported position of Mesodinium as sister to other litostomes contrasts with previous molecular analyses including those from phylogenomic studies that exploited the same transcriptomic databases. These topological discrepancies illustrate the need for caution when mining mixed-species transcriptomes and indicate that identifying ciliate sequences among prey contamination-particularly for Mesodinium species where expression from stolen prey nuclei appears to dominate-requires thorough and iterative vetting with phylogenies that incorporate sequences from a large outgroup of prey.

The New Tree of Eukaryotes

For 15 years, the eukaryote Tree of Life (eToL) has been divided into five to eight major groupings, known as 'supergroups'. However, the tree has been profoundly rearranged during this time. The new eToL results from the widespread application of phylogenomics and numerous discoveries of major lineages of eukaryotes, mostly free-living heterotrophic protists. The evidence that supports the tree has transitioned from a synthesis of molecular phylogenetics and biological characters to purely molecular phylogenetics. Most current supergroups lack defining morphological or cell-biological characteristics, making the supergroup label even more arbitrary than before. Going forward, the combination of traditional culturing with maturing culture-free approaches and phylogenomics should accelerate the process of completing and resolving the eToL at its deepest levels.

A multigene timescale and diversification dynamics of Ciliophora evolution

Ciliophora is one of the most diverse lineages of unicellular eukaryotes. Nevertheless, a robust timescale including all main lineages and employing properly identified ciliate fossils as primary calibrations is lacking. Here, we inferred a time-calibrated multigene phylogeny of Ciliophora evolution, and we used this timetree to investigate the rates and patterns of lineage diversification through time. We implemented a two-step analytical approach that favored both gene and taxon sampling, reducing the uncertainty of time estimates and yielding narrower credibility intervals on the ribosomal-derived chronogram. We estimate the origin of Ciliophora at 1143 Ma, which is substantially younger than previously proposed ages, and the huge diversity explosion occurred during the Paleozoic. Among the current groups recognized as classes, Spirotrichea diverged earlier, its origin was dated at ca. 850 Ma, and Protocruziea was the younger class, with crown age estimated at 56 Ma. Macroevolutionary analysis detected a significant rate shift in diversification dynamics in the spirotrichean clade Hypotrichia + Oligotrichia + Choreotrichia, which had accelerated speciation rate ca. 570 Ma, during the Ediacaran-Cambrian transition. For all crown lineages investigated, speciation rates declined through time, whereas extinction rates remained low and relatively constant throughout the evolutionary history of ciliates.

Ultrastructural Studies on a Model Tintinnid - Schmidingerella meunieri (Kofoid and Campbell, 1929) Agatha and Strüder-Kypke, 2012 (Ciliophora). I. Somatic Kinetids with Unique Ultrastructure

Molecular phylogenies of Oligotrichea currently do not contain all genera and families and display topologies which are often incongruent with morphological findings. In ciliates, the somatic kinetids are rather conserved, i.e., their ultrastructures, particularly the fibrillar associates, often characterise the main groups, except for the choreotrichids. Four different kinetid types are found in protargol-stained choreotrichids and used for reconstructing the taxon's evolution (the "Kinetid Transformation Hypothesis"). Proof for this hypothesis requires transmission electron microscopic studies, which are very rare in the choreotrichids and oligotrichids. Such an approach provides insights into the ultrastructural variability of somatic kinetids in spirotrichs and may also detect apomorphies characterising certain choreotrichid families. In the model tintinnid Schmidingerella meunieri, the ultrastructure of the three kinetid types in the somatic ciliature is studied in cryofixed cells. The data support the "Kinetid Transformation Hypothesis" regarding tintinnids with a ventral kinety. This first detailed study on kinetids in tintinnids and choreotrichids in general reveals totally new kinetid types in ciliates: beyond the three common associates, they are characterised by two or three conspicuous microtubular ribbons extending on the kinetids' left sides. These extraordinary ribbons form together with the overlapping postciliary ribbons a unique network in the cortex of the anterior cell portion. The evolutionary constrains which might have fostered the development of such structures are discussed for the Oligotrichea, the choreotrichids, and tintinnids as their first occurrence is currently uncertain. Additionally, the kinetids in tintinnids, aloricate choreotrichids, oligotrichids, hypotrichs, and euplotids are compared.

Techniques and tools for species identification in ciliates: a review

Ciliates are highly divergent unicellular eukaryotic organisms with nuclear dualism and a highly specialized ciliary pattern. They inhabit all biotopes and play crucial roles in regulating microbial food webs as they prey on bacteria, protists and even on microscopic animals. Nevertheless, subtle morphological differences and tiny sizes hinder proper species identification for many ciliates. In the present review, an attempt has been made to elaborate the various approaches used by modern day ciliate taxonomists for species identification. The different approaches involved in taxonomic characterization of ciliates such as classical (using live-cell observations, staining techniques, etc.), molecular (involving various marker genes) and statistical (delimitation of cryptic species) methods have been reviewed. Ecological and behavioural aspects in species identification have also been discussed. In present-day taxonomy, it is important to use a ‘total evidence’ approach in identifying ciliates, relying on both classical and molecular information whenever possible. This integrative approach will help in the mergence of classical methods with modern-day tools for comprehensive species description in future.

Genetic basis for the establishment of endosymbiosis in Paramecium

The single-celled ciliate Paramecium bursaria is an indispensable model for investigating endosymbiosis between protists and green-algal symbionts. To elucidate the mechanism of this type of endosymbiosis, we combined PacBio and Illumina sequencing to assemble a high-quality and near-complete macronuclear genome of P. bursaria. The genomic characteristics and phylogenetic analyses indicate that P. bursaria is the basal clade of the Paramecium genus. Through comparative genomic analyses with its close relatives, we found that P. bursaria encodes more genes related to nitrogen metabolism and mineral absorption, but encodes fewer genes involved in oxygen binding and N-glycan biosynthesis. A comparison of the transcriptomic profiles between P. bursaria with and without endosymbiotic Chlorella showed differential expression of a wide range of metabolic genes. We selected 32 most differentially expressed genes to perform RNA interference experiment in P. bursaria, and found that P. bursaria can regulate the abundance of their symbionts through glutamine supply. This study provides novel insights into Paramecium evolution and will extend our knowledge of the molecular mechanism for the induction of endosymbiosis between P. bursaria and green algae.

Twisted Tales: Insights into Genome Diversity of Ciliates Using Single-Cell 'Omics

The emergence of robust single-cell 'omics techniques enables studies of uncultivable species, allowing for the (re)discovery of diverse genomic features. In this study, we combine single-cell genomics and transcriptomics to explore genome evolution in ciliates (a > 1 Gy old clade). Analysis of the data resulting from these single-cell 'omics approaches show: 1) the description of the ciliates in the class Karyorelictea as "primitive" is inaccurate because their somatic macronuclei contain loci of varying copy number (i.e., they have been processed by genome rearrangements from the zygotic nucleus); 2) gene-sized somatic chromosomes exist in the class Litostomatea, consistent with Balbiani's (1890) observation of giant chromosomes in this lineage; and 3) gene scrambling exists in the underexplored Postciliodesmatophora (the classes Heterotrichea and Karyorelictea, abbreviated here as the Po-clade), one of two major clades of ciliates. Together these data highlight the complex evolutionary patterns underlying germline genome architectures in ciliates and provide a basis for further exploration of principles of genome evolution in diverse microbial lineages.

EukRef-Ciliophora: a manually curated, phylogeny-based database of small subunit rRNA gene sequences of ciliates

High-throughput sequencing (HTS) surveys, among the most common approaches currently used in environmental microbiology, require reliable reference databases to be correctly interpreted. The EukRef Initiative (eukref.org) is a community effort to manually screen available small subunit (SSU) rRNA gene sequences and produce a public, high-quality and informative framework of phylogeny-based taxonomic annotations. In the context of EukRef, we present a database for the monophyletic phylum Ciliophora, one of the most complex, diverse and ubiquitous protist groups. We retrieved more than 11 500 sequences of ciliates present in GenBank (28% from identified isolates and 72% from environmental surveys). Our approach included the inference of phylogenetic trees for every ciliate lineage and produced the largest SSU rRNA tree of the phylum Ciliophora to date. We flagged approximately 750 chimeric or low-quality sequences, improved the classification of 70% of GenBank entries and enriched environmental and literature metadata by 30%. The performance of EukRef-Ciliophora is superior to the current SILVA database in classifying HTS reads from a global marine survey. Comprehensive outputs are publicly available to make the new tool a useful guide for non-specialists and a quick reference for experts.

Remodeling the Specificity of an Endosomal CORVET Tether Underlies Formation of Regulated Secretory Vesicles in the Ciliate Tetrahymena thermophila

In the endocytic pathway of animals, two related complexes, called CORVET (class C core vacuole/endosome transport) and HOPS (homotypic fusion and protein sorting), act as both tethers and fusion factors for early and late endosomes, respectively. Mutations in CORVET or HOPS lead to trafficking defects and contribute to human disease, including immune dysfunction. HOPS and CORVET are conserved throughout eukaryotes, but remarkably, in the ciliate Tetrahymena thermophila, the HOPS-specific subunits are absent, while CORVET-specific subunits have proliferated. VPS8 (vacuolar protein sorting), a CORVET subunit, expanded to 6 paralogs in Tetrahymena. This expansion correlated with loss of HOPS within a ciliate subgroup, including the Oligohymenophorea, which contains Tetrahymena. As uncovered via forward genetics, a single VPS8 paralog in Tetrahymena (VPS8A) is required to synthesize prominent secretory granules called mucocysts. More specifically, Δvps8a cells fail to deliver a subset of cargo proteins to developing mucocysts, instead accumulating that cargo in vesicles also bearing the mucocyst-sorting receptor Sor4p. Surprisingly, although this transport step relies on CORVET, it does not appear to involve early endosomes. Instead, Vps8a associates with the late endosomal/lysosomal marker Rab7, indicating that target specificity switching occurred in CORVET subunits during the evolution of ciliates. Mucocysts belong to a markedly diverse and understudied class of protist secretory organelles called extrusomes. Our results underscore that biogenesis of mucocysts depends on endolysosomal trafficking, revealing parallels with invasive organelles in apicomplexan parasites and suggesting that a wide array of secretory adaptations in protists, like in animals, depend on mechanisms related to lysosome biogenesis.

Phylogenomic analysis of Balantidium ctenopharyngodoni (Ciliophora, Litostomatea) based on single-cell transcriptome sequencing

In this paper, we present transcriptome data for Balantidium ctenopharyngodoni Chen, 1955 collected from the hindgut of grass carp (Ctenopharyngodon idella). We evaluated sequence quality and de novo assembled a preliminary transcriptome, including 43.3 megabits and 119,141 transcripts. Then we obtained a final transcriptome, including 17.7 megabits and 35,560 transcripts, by removing contaminative and redundant sequences. Phylogenomic analysis based on a supermatrix with 132 genes comprising 53,873 amino acid residues and phylogenetic analysis based on SSU rDNA of 27 species were carried out herein to reveal the evolutionary relationships among six ciliate groups: Colpodea, Oligohymenophorea, Litostomatea, Spirotrichea, Heterotrichea and Protocruziida. The topologies of both phylogenomic and phylogenetic trees are discussed in this paper. In addition, our results suggest that single-cell sequencing is a sound method of obtaining sufficient omics data for phylogenomic analysis, which is a good choice for uncultivable ciliates. The transcriptome data for Balantidium ctenopharyngodoni are the first omics data within the subclass Trichostomatia, and provide a good basis for ciliate phylogenomic analysis, as well as related omics analysis.

Molecules illuminate morphology: phylogenomics confirms convergent evolution among 'oligotrichous' ciliates

'Oligotrichous' ciliates have been traditionally placed in a presumed monophyletic taxon called the Oligotrichia. However, gene sequences of the small subunit rRNA gene, and several other genes, suggest that the taxon is not monophyletic: although statistical support for this is not strong, the oligotrich Halteria grandinella is associated with the hypotrich ciliates and not with other oligotrich genera, such as Strombidium and Strombidinopsis. This has convinced some taxonomists to emphasize that morphological features strongly support the monophyly of the oligotrichs. To further test this hypothesis of monophyly, we have undertaken a phylogenomic analysis using the transcriptome of H. grandinella cells amplified by a single-cell technique. One hundred and twenty-six of 159 single-gene trees placed H. grandinella as sister to hypotrich species, and phylogenomic analyses based on a subset of 124 genes robustly rejected the monophyly of the Oligotrichia and placed the genus Halteria as sister to the hypotrich genera Stylonychia and Oxytricha. We use these phylogenomic analyses to assess the convergent nature of morphological features of oligotrichous ciliates. A particularly 'strong' morphological feature supporting monophyly of the oligotrichs is enantiotropic cell division, which our results suggest is nevertheless a convergent feature, arising through the need for dividing ciliates to undertake rotokinesis to complete cell division.

Nuclear genetic codes with a different meaning of the UAG and the UAA codon

BACKGROUND

Departures from the standard genetic code in eukaryotic nuclear genomes are known for only a handful of lineages and only a few genetic code variants seem to exist outside the ciliates, the most creative group in this regard. Most frequent code modifications entail reassignment of the UAG and UAA codons, with evidence for at least 13 independent cases of a coordinated change in the meaning of both codons. However, no change affecting each of the two codons separately has been documented, suggesting the existence of underlying evolutionary or mechanistic constraints.

RESULTS

Here, we present the discovery of two new variants of the nuclear genetic code, in which UAG is translated as an amino acid while UAA is kept as a termination codon (along with UGA). The first variant occurs in an organism noticed in a (meta)transcriptome from the heteropteran Lygus hesperus and demonstrated to be a novel insect-dwelling member of Rhizaria (specifically Sainouroidea). This first documented case of a rhizarian with a non-canonical genetic code employs UAG to encode leucine and represents an unprecedented change among nuclear codon reassignments. The second code variant was found in the recently described anaerobic flagellate Iotanema spirale (Metamonada: Fornicata). Analyses of transcriptomic data revealed that I. spirale uses UAG to encode glutamine, similarly to the most common variant of a non-canonical code known from several unrelated eukaryotic groups, including hexamitin diplomonads (also a lineage of fornicates). However, in these organisms, UAA also encodes glutamine, whereas it is the primary termination codon in I. spirale. Along with phylogenetic evidence for distant relationship of I. spirale and hexamitins, this indicates two independent genetic code changes in fornicates.

CONCLUSIONS

Our study documents, for the first time, that evolutionary changes of the meaning of UAG and UAA codons in nuclear genomes can be decoupled and that the interpretation of the two codons by the cytoplasmic translation apparatus is mechanistically separable. The latter conclusion has interesting implications for possibilities of genetic code engineering in eukaryotes. We also present a newly developed generally applicable phylogeny-informed method for inferring the meaning of reassigned codons.