Phylogenomics solves a long-standing evolutionary puzzle in the ciliate world: The subclass Peritrichia is monophyletic

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.

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.

On four epibiotic peritrichous ciliates (Protozoa, Ciliophora) found in Lake Weishan Wetland: morphological and molecular data support the establishment of a new genus, Parapiosoma gen. nov., and two new species

During a study on the diversity of ciliated protists in Lake Weishan Wetland, the largest wetland in northern China, four epibiotic sessilid peritrichs were isolated from aquatic host animals. Two of them, i.e., Epistylis cambari Kellicott, 1885 and Epistylis lwoffi Fauré-Fremiet, 1943, were known species whereas the other two, i.e., Parapiosoma typicum gen. nov., sp. nov. and Orborhabdostyla gracilis sp. nov., are new to science. The new genus Parapiosoma gen. nov. is characterized by its branched non-contractile stalk, everted peristomial lip, obconical macronucleus and transverse silverlines. Two species are assigned to the new genus, namely Parapiosoma typicum sp. nov. and Parapiosoma gasterostei (Fauré-Fremiet, 1905) comb. nov. Morphologically, P. typicum sp. nov. is recognized by its goblet-shaped zooids, single-layered peristomial lip, dichotomously branched stalk, and infundibular polykinety 3 (P3) containing three equal-length rows. Orborhabdostyla gracilis sp. nov. is characterized by its slender zooid, curved macronucleus, and three equal-length rows in infundibular P3. Improved diagnoses and redescriptions of E. cambari and E. lwoffi are provided including, for the first time, data on the ciliature of E. cambari. Phylogenetic analyses based on SSU rDNA, ITS1-5.8S rDNA -ITS2, and LSU rDNA sequence data strongly support the assertion that the family Epistylididae comprises morphospecies with different evolutionary lineages and indicate that Parapiosoma gen. nov. may represent a new taxon at family level.

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.

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.

Large-scale phylogenetic analysis provides insights into the diversification and evolution of sessilid peritrich ciliates (Protista: Ciliophora)

The Peritrichia is a speciose and morphologically distinctive assemblage of ciliated protists that was first observed by Antonie van Leeuwenhoek over 340 years ago. In the last two decades, the phylogenetic relationships of this group have been increasingly debated as morphological and molecular analyses have generated contrasting conclusions, mainly owing to limited sampling. In the present study, we performed expanded phylogenetic analyses of 152 sessilid peritrichs collected from 14 different provinces of China and 141 SSU rDNA peritrich sequences from GenBank. The results of the analyses revealed new divergent relationships between and within major clades that challenge the morphological classification of this group including, (1) the recovery of four major phylogenetically divergent clades in the monophyletic order Sessilida, (2) aboral structures such as the stalk and spasmoneme were evolutionary labile, (3) the stalk or/and spasmoneme was lost in each divergent clade indicating that parallel evolution occurred in sessilid peritrichs and (4) the life cycle and habit drive the diversity of aboral structures as well as diversification and evolution in peritrichs.

New considerations on the phylogeny of Sessilida (Protista: Ciliophora: Peritrichia) based on multiple-gene information, with emphasis on colonial taxa

The subclass Peritrichia, containing two orders Sessilida and Mobilida, is a major group of ciliates with worldwide distribution and high species diversity. Several studies have investigated the phylogeny of peritrichs; however, the evolutionary relationships and classification of some families and genera within the Sessilida remain unclear. In the present study, we isolated and identified 22 peritrich populations representing four families and six genera and obtained 64 rDNA sequences to perform phylogenetic analyses and assess their systematic relationships. Ancestral character reconstruction was also carried out to infer evolutionary routes within the Sessilida. The results indicate: (1) family Vaginicolidae is monophyletic and acquisition of the typical peritrich lorica represents a single evolutionary divergence; (2) core epistylidids evolved from a Zoothamnium-like ancestor and experienced spasmoneme loss during evolution; (3) Campanella clusters with species in the basal clade and shows stable morphological differences with other epistylidids, supporting its assignment to a separate family; (4) the structure of the peristomial lip may be a genus-level character rather than a diagnostic character for discriminating Epistylididae and Operculariidae, thus a redefinition of Operculariidae should be carried out when more species have been investigated; (5) some characters, such as lifestyle (solitary or colonial), spasmoneme and living habit (sessile or free-swimming), evolved repeatedly among sessilids indicating that species with non-contractile stalks or that are free-swimming have multiple evolutionary routes and might derive from any sessilid lineage without a lorica. The close phylogenetic relationships of some morphologically distinct sessilids imply that the diagnoses of some genera and families should be improved.

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.

The taxonomy and molecular phylogeny of two epibiotic colonial peritrich ciliates (Ciliophora, Peritrichia)

Two colonial sessilid peritrichs, Epistylis qingdaoensis n. sp. and Carchesium cyclopidarum Nenninger, 1948, were isolated from a marine crustacean and a freshwater mayfly, respectively. Morphological characters for each species were revealed by in vivo observations and silver staining methods. Epistylis qingdaoensis n. sp. is characterized by the apperance of the colony which is up to 250 μm high and usually contains fewer than eight zooids, the single-layered peristomial lip, the conspicuously conical peristomial disc, and the structure of infundibular polykinety 3 which comprises three isometric ciliary rows. Carchesium cyclopidarum is recognized by the colony height of about 600 μm, the small zooid size of 35-50 × 20-30 μm, and the low number of silverlines. Phylogenetic analyses based on small subunit ribosomal DNA sequences were performed to reveal their evolutionary relationships. Surprisingly, neither species clustered with its congeners. Epistylis qingdaoensis n. sp. nested within a clade of Zoothamnium species that was sister to a clade comprising core Epistylis species and several Zoothamnium species. Carchesium cyclopidarum nested within a clade comprising solitary peritrichs and was distantly related to the type species of Carchesium, C. polypinum. These findings imply for the first time that Carchesium is non-monophyletic.

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.

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

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

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 species of Trichodina lishuiensis n. sp. (Ciliophora: Trichodinidae) in urinary bladder of Odorrana schmackeri (Amphibia: Ranidae) from Zhejiang, China

Endoparasitic trichodinids are rather rare ciliates. In this study we describe a new species named Trichodina lishuiensis from the bladder of Odorrana schmackeri collected in Zhejiang, China, with the prevalence of 20% (9/45). We identified T. lishuiensis as a new species by morphological comparison and molecular analysis. The ciliates were observed using the dry-silver and protargol staining methods, as well as SEM (scanning electron microscopy). Trichodina lishuiensis is a small species (cell diameter 31.8-43.9 μm), with incompact denticles connection, medium-wide blades and thick rays. We also sequenced a 1712 bp-long fragment of the small subunit ribosomal RNA gene (SSU rRNA). Phylogenetic analyses showed that the new species clustered with Trichodina unionis. The route of transmission of Trichodina species in the urinary bladder remains a mystery. We hypothesize that the transmisison takes place during the amplexus, with eggs and sperm discharged from the cloaca, and that trichodinids 'accompany' the amphibian through its whole life cycle, but further studies are needed to test this hypothesis.

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.

Integrative Studies on a New Ciliate Campanella sinica n. sp. (Protista, Ciliophora, Peritrichia) Based on the Morphological and Molecular Data, With Notes on the Phylogeny and Systematics of the Family Epistylididae

During an investigation on freshwater peritrichs, a new colonial sessilid ciliate, Campanella sinica n. sp., was isolated from aquatic plants in an artificial freshwater pond in Qingdao, China. Specimen observations of this species were performed both in vivo and using silver staining. C. sinica n. sp. is characterized by the appearance of the mature colony, which is up to 2 cm high and contains more than 1,000 zooids, the asymmetric horn-shaped zooids, strongly everted and multi-layered peristomial lip, the slightly convex peristomial disc, and the well-developed haplokinety and polykinety, which make more than four circuits of the peristome before descending into the infundibulum. The small subunit ribosomal DNA (SSU rDNA), 5.8s rDNA and its flank internal transcribed spacers (ITS1-5.8s rDNA-ITS2), and large subunit ribosomal DNA (LSU rDNA) are sequenced and used for phylogenetic analyses which reveal that the family Epistylididae Kahl, 1933 is non-monophyletic whereas the genus Campanella is monophyletic and nests within the basal clade of the sessilids. The integrative results support the assertion that the genus Campanella represents a separate lineage from other epistylidids, suggesting a further revision of the family Epistylididae is needed. We revise Campanella including the transfer into this genus of a taxon formerly assigned to Epistylis, which we raise to species rank, i.e., Campanella ovata (Nenninger, 1948) n. grad. & n. comb. (original combination Epistylis purneri f. ovata Nenninger, 1948). In addition, we provide a key to the identification of the species of Campanella.

Cryptic host-driven speciation of mobilid ciliates epibiotic on freshwater planarians

Mobilids are among the most taxonomically diverse but morphologically uniform groups of epibiotic ciliates. They attach to their hosts by means of an adhesive disc as harmless commensals such as Urceolaria, or as parasites causing significant economic loss such as some Trichodina species. We investigated the diversity, species boundaries, and phylogenetic relationships of mobilids associated with freshwater planarians, using 114 new sequences of two mitochondrial (16S rRNA gene and cytochrome c oxidase gene) and five nuclear (18S rRNA gene, ITS1-5.8S-ITS2 region, D1/D2 domains of 28S rRNA gene) markers. Although the morphological disparity of the isolated trichodinids and urceolariids was low, Bayesian coalescent analyses revealed the existence of five distinct evolutionary lineages/species given the seven molecular markers. The occurrence of mobilids perfectly correlated with their planarian hosts: Trichodina steinii and two Urceolaria mitra-like taxa were associated exclusively with the planarian Dugesia gonocephala, Trichodina polycelis sp. n. with the planarian Polycelis felina, and Trichodina schmidtea sp. n. with the planarian Schmidtea polychroa. Host organisms thus very likely constitute sharply isolated niches that might permit speciation of their epibiotic ciliates, even though no distinct morphological features appear to be recognizable among ciliates originating from different hosts.

Ciliate Diversity From Aquatic Environments in the Brazilian Atlantic Forest as Revealed by High-Throughput DNA Sequencing

Rainforest aquatic ecosystems include complex habitats with scarce information on their unicellular eukaryote diversity and community structure. We have investigated the diversity of ciliates in freshwater and brackish environments along the Brazilian Atlantic Forest, based on the hypervariable V4 region of the 18S-rDNA obtained by high-throughput DNA sequencing. Our analyses detected 409 ciliate taxonomic units (OTUs), mostly attributed to the classes Oligohymenophorea and Spirotrichea. A total of 11 classes, 12 subclasses, 112 genera, and 144 species were reported. We found the following: (a) the ciliate communities are more diverse in freshwater- than in Atlantic Forest-associated brackish environments; (b) the ciliate communities are composed by a small amount of highly abundant OTUs, but a high number of low-abundant or rare OTUs; (c) nearly one-third of the ciliate OTUs share less than 97% sequence identity to reference sequences and (d) phylogenetic inference supports the hypothesis that the V4 region of the Ciliophora 18S-rDNA is a suitable marker for accurate evolutionary inferences at class level. Our results showed that a considerable fraction of the HTS-detected diversity of ciliates from Brazilian Atlantic Forest is not represented in the currently available molecular databases.

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.

Re-discovery and novel contributions to morphology and multigene phylogeny of Myxophyllum steenstrupi (Ciliophora: Pleuronematida), an obligate symbiont of terrestrial pulmonates

Myxophyllum steenstrupi is a symbiotic ciliate living in the body slime and mantle cavity of terrestrial pulmonates (Gastropoda: Pulmonata). In the present study, M. steenstrupi was re-discovered after almost 30 years and characterized using an integrative morpho-molecular approach for the first time. Myxophyllum is distinguished by a broadly ovate, about 140 × 115 μm-sized body, a nuclear apparatus typically composed of seven macronuclear nodules and a single micronucleus, a central contractile vacuole, a shallow oral cavity situated in the posterior body region and dense somatic ciliature with extensive thigmotactic field. According to the present phylogenetic analyses of two mitochondrial and five nuclear markers, M. steenstrupi is classified in the predominantly free-living order Pleuronematida (Oligohymenophorea: Scuticociliatia). This order also encompasses other taxa isolated from molluscs and traditionally classified along with Myxophyllum in the order Thigmotrichida. The proper classifications of Myxophyllum was hampered by the dramatic remodelling of its oral apparatus (reduction of the paroral membrane and adoral organelles, formation of vestibular kineties), a transformation that was likely induced by its firm association with terrestrial gastropods. The present study also documents that various ciliate lineages independently became commensals or parasites of various aquatic and terrestrial molluscs.

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.

Symposium on Ciliates in Memory of Denis Lynn

Denis Lynn (1947-2018) was an outstanding protistologist, applying multiple techniques and data sources and thus pioneering an integrative approach in order to investigate ciliate biology. For example, he recognized the importance of the ultrastructure for inferring ciliate phylogeny, based on which he developed his widely accepted classification scheme for the phylum Ciliophora. In this paper, recent findings regarding the evolution and systematics of both peritrichs and the mainly marine planktonic oligotrichean spirotrichs are discussed and compared with the concepts and hypotheses formulated by Denis Lynn. Additionally, the state of knowledge concerning the diversity of ciliates in bromeliad phytotelmata and amitosis in ciliates is reviewed.

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.

Trichostomatid Ciliates (Alveolata, Ciliophora, Trichostomatia) Systematics and Diversity: Past, Present, and Future

The gastrointestinal tracts of most herbivorous mammals are colonized by symbiotic ciliates of the subclass Trichostomatia, which form a well-supported monophyletic group, currently composed by ∼1,000 species, 129 genera, and 21 families, distributed into three orders, Entodiniomorphida, Macropodiniida, and Vestibuliferida. In recent years, trichostomatid ciliates have been playing a part in many relevant functional studies, such as those focusing in host feeding efficiency optimization and those investigating their role in the gastrointestinal methanogenesis, as many trichostomatids are known to establish endosymbiotic associations with methanogenic Archaea. However, the systematics of trichostomatids presents many inconsistencies. Here, we stress the importance of more taxonomic works, to improve classification schemes of this group of organisms, preparing the ground to proper development of such relevant applied works. We will present a historical review of the systematics of the subclass Trichostomatia highlighting taxonomic problems and inconsistencies. Further on, we will discuss possible solutions to these issues and propose future directions to leverage our comprehension about taxonomy and evolution of these symbiotic microeukaryotes.

Morphology and molecular analyses of four epibiotic peritrichs on crustacean and polychaete hosts, including descriptions of two new species (Ciliophora, Peritrichia)

Four epibiotic sessilid peritrichs, i.e., Zoothamnium wilberti n. sp., Baikalonis microdiscus n. sp., Epistylis anastatica (Linnaeus, 1767) Ehrenberg, 1830, and Rhabdostyla commensalisMöbius, 1888, were isolated from one syllid polychaete and three crustacean hosts in Qingdao, China. For each species, specimens were observed both in vivo and following silver staining. Their SSU rDNA was also sequenced for phylogenetic analyses. Zoothamnium wilberti n. sp. is characterized by the appearance of its colony, which is up to 350 μm high, and usually has fewer than 16 zooids, and the dichotomously branched stalk with transverse wrinkles, the conspicuously conical peristomial disc, and infundibular polykinety 3 comprising three isometric ciliary rows. Baikalonis microdiscus n. sp. can be recognized by its barrel-shaped zooid, small peristomial disc, smooth and short stalk, and its unusual infundibular polykinety 3 comprising a long inner row and a short outer row. Two poorly known species, i.e., Epistylis anastatica and Rhabdostyla commensalis, are redescribed and redefined. Phylogenetic analyses reveal that: (i) R. commensalis is closely related to the family Astylozoidae rather than to the morphologically similar Epistylididae; (ii) B. microdiscus n. sp. is sister to the family Scyphidiidae; (iii) E. anastatica groups with vorticellids and ophrydiids, which further supports the polyphyly of the genus Epistylis; and (iv) Z. wilberti n. sp. is nested within the Zoothamniidae, as expected.

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.

Two Urosoma species (Ciliophora, Hypotrichia): A multidisciplinary approach provides new insights into their ultrastructure and systematics

The general morphology and ultrastructure of two soil hypotrichous ciliates, Urosoma emarginata and U. salmastra, were investigated using light microscopy, scanning electron microscopy and transmission electron microscopy. Phylogenetic analyses, based on the newly sequenced small subunit ribosomal (SSU) rRNA genes, were conducted on three U. emarginata populations and one U. salmastra population. Our findings support for the validity of Perilemmaphora Berger, 2008, a rankless taxon comprising spirotrich ciliates having a perilemma. The cortical granules of both species are extrusomes representing a new type of mucocyst in U. emarginata and possibly a new type of pigmentocyst in U. salmastra. Additionally, the lithosomes were revealed as subglobose structures composed of a low electron-dense, homogeneous inner part and an electron-dense outer part. The ultrastructural features of the cortical granules, together with ontogenetic and molecular phylogenetic data, suggest that the genus Urosoma might need to be divided. It is posited that ultrastructural features of hypotrichous ciliates in general may have important taxonomic value warranting further investigation.

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.

Morphological and molecular identification of epibiontic sessilid Epistylis semiciculus n. sp. (ciliophora, Peritrichia) from Procambarus clarkia (Crustacea, Decapoda) in China

Data on sessilinasis from Procambarus clarkia are limited. This study investigates sessilid diversity in Hubei Province, China in 2016. Procambarus clarkia pereopods were covered by two sessilid morphotypes in April and May, and the gills were adhered by one of these two morphotypes (morphotype I) in January. Identifying the two morphotypes according to their morphological characters through live observations and protargol-stained method is difficult. Both morphotypes show almost identical morphological characteristics: zooids present vase-shaped, the length-to-width ratio is 2:1, the semicircle-shaped peristomial disk is evidently above the peristomial lip, single contractile vacuole is located below the peristomial lip, and the oral infraciliature shows identical arrangement which can be distinguished from other Epistylis. However, the morphometric data of the two morphotypes are significant different: body size of morphotype I is 78.9–103.8 × 32.1–54.6 μm, whereas that of morphotypes II is 136.7–171.5 × 60.9–88.0 μm. To further identify the two morphotypes, molecular regions including small subunit ribosomal DNA (SSU rDNA) sequences, large subunit ribosomal DNA (LSU rDNA) sequences, and ITS1-5.8S-ITS2 sequences were used. Results supported that the two morphotypes are single species of genus Epistylis rather than distinct species based in their distinct size ranges and temporal presence events. We assigned it the name E. semiciculus n. sp. with respect to the semicircular peristomial disk. Overall, these findings emphasized the importance of using molecular data to solve the identification confusion caused by ontogenetic processes. This study is the first morphological and molecular characterization of a sessilid isolated from P. clarkia under aquaculture conditions.

•

The first identification of epibiontic sessilid from Procambarus clarkia.

•

Different attachment sites of sessilids in spring and winter.

•

Two morphotypes of a new Epistylis with significant different morphometric data.

•

Species descriptions include morphological (vivo and stained) and molecular data.

•

Present Epistylis distinguished by body shape, contractile vacuole and infraciliature.

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.

Further analyses on the phylogeny of the subclass Scuticociliatia (Protozoa, Ciliophora) based on both nuclear and mitochondrial data

So far, the phylogenetic studies on ciliated protists have mainly based on single locus, the nuclear ribosomal DNA (rDNA). In order to avoid the limitations of single gene/genome trees and to add more data to systematic analyses, information from mitochondrial DNA sequence has been increasingly used in different lineages of ciliates. The systematic relationships in the subclass Scuticociliatia are extremely confused and largely unresolved based on nuclear genes. In the present study, we have characterized 72 new sequences, including 40 mitochondrial cytochrome oxidase c subunit I (COI) sequences, 29 mitochondrial small subunit ribosomal DNA (mtSSU-rDNA) sequences and three nuclear small subunit ribosomal DNA (nSSU-rDNA) sequences from 47 isolates of 44 morphospecies. Phylogenetic analyses based on single gene as well as concatenated data were performed and revealed: (1) compared to mtSSU-rDNA, COI gene reveals more consistent relationships with those of nSSU-rDNA; (2) the secondary structures of mtSSU-rRNA V4 region are predicted and compared in scuticociliates, which can contribute to discrimination of closely related species; (3) neither nuclear nor mitochondrial data support the monophyly of the order Loxocephalida, which may represent some divergent and intermediate lineages between the subclass Scuticociliatia and Hymenostomatia; (4) the assignments of thigmotrichids to the order Pleuronematida and the confused taxon Sulcigera comosa to the genus Histiobalantium are confirmed by mitochondrial genes; (5) both nuclear and mitochondrial data reveal that the species in the family Peniculistomatidae always group in the genus Pleuronema, suggesting that peniculistomatids are more likely evolved from Pleuronema-like ancestors; (6) mitochondrial genes support the monophyly of the order Philasterida, but the relationships among families of the order Philasterida remain controversial due to the discrepancies between their morphological and molecular data.

Exploration of Genetic Variations through Single-cell Whole-genome Sequencing in the Model Ciliate Tetrahymena thermophila

Ciliates are unicellular eukaryotes with separate germline and somatic genomes and diverse life cycles, which make them a unique model to improve our understanding of population genetics through the detection of genetic variations. However, traditional sequencing methods cannot be directly applied to ciliates because the majority are uncultivated. Single-cell whole-genome sequencing (WGS) is a powerful tool for studying genetic variation in microbes, but no studies have been performed in ciliates. We compared the use of single-cell WGS and bulk DNA WGS to detect genetic variation, specifically single nucleotide polymorphisms (SNPs), in the model ciliate Tetrahymena thermophila. Our analyses showed that (i) single-cell WGS has excellent performance regarding mapping rate and genome coverage but lower sequencing uniformity compared with bulk DNA WGS due to amplification bias (which was reproducible); (ii) false-positive SNP sites detected by single-cell WGS tend to occur in genomic regions with particularly high sequencing depth and high rate of C:G to T:A base changes; (iii) SNPs detected in three or more cells should be reliable (an detection efficiency of 83.4-97.4% was obtained for combined data from three cells). This analytical method could be adapted to measure genetic variation in other ciliates and broaden research into ciliate population genetics.

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.

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.

Evolutionary Associations of Endosymbiotic Ciliates Shed Light on the Timing of the Marsupial-Placental Split

Trichostome ciliates are among the most conspicuous protists in the gastrointestinal tract of a large variety of vertebrates. However, little is still known about phylogeny of the trichostome/vertebrate symbiotic systems, evolutionary correlations between trichostome extrinsic traits, and character-dependent diversification of trichostomes. These issues were investigated here, using the relaxed molecular clock technique along with stochastic mapping of character evolution, and binary-state speciation and extinction models. Clock analyses revealed that trichostomes colonized the vertebrate gastrointestinal tract ∼135 Ma, that is, near the paleontological minimum for the split of therian mammals into marsupials and placentals. According to stochastic mapping, the last common ancestor of trichostomes most likely invaded the hindgut of a mammal. Although multiple shifts to fish/amphibian or avian hosts and to the foregut compartments took place during the trichostome phylogeny, only transition to the foregut was recognized as a key innovation responsible for the explosive radiation of ophryoscolecid trichostomes after the Cretaceous/Tertiary boundary, when ungulates began their diversification. Since crown radiations of main trichostome lineages follow those of their mammalian hosts and are in agreement with their historic dispersal routes, the present time-calibrated phylogeny might help to elucidate controversies in the geological and molecular timing of the split between marsupials and placental mammals.

Molecular identification of Trichodina compacta Van As and Basson, 1989 (Ciliophora: Peritrichia) from cultured Oreochromis niloticus in Egypt and its impact on immune responses and tissue pathology

Trichodinids are peritrichous ciliated protozoa that affect both wild and cultured fishes. Several Trichodina species have low host specificity and are morphologically distinct, facilitating their identification based primarily on the presence of adhesive discs and the number of attached denticles. A trichodinid species named Trichodina compacta was first reported by Van As and Basson (1989) (Protozoa: Ciliophora: Peritrichia). However, in trichodinid infestations, morphological characteristics are insufficient for identifying the infesting species. Therefore, molecular and phylogenetic analyses are considered to be promising and useful tools for identifying the infesting species. This study aimed to achieve the molecular identification of a trichodinid infestation in Nile tilapia and to construct the phylogenetic relationships between the identified species and other peritrichous parasites. Moreover, we also aimed to study the pathological and immunological impacts of trichodinids on fry tissue to improve our understanding of the immune responses of teleost fish to trichodinae parasitic infestations and develop a better control method. Here, we used molecular techniques to identify the isolated trichodina species as T. compacta and demonstrated that Trichodina infestation in Nile tilapia is associated with remarkable immunogenic and inflammatory responses (increased il-1β expression and decreased il-8 and tgf-β expression). These findings improve our understanding of the responses of teleost fish to trichodinid parasite infestation and will be helpful for the development of novel control strategies that reverse the inflammatory and immunogenic alterations that occur in infested fish.