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

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.

Genomic insights into the cellular specialization of predation in raptorial protists

BACKGROUND

Predation is a fundamental mechanism for organisms to acquire energy, and various species have evolved diverse tools to enhance their hunting abilities. Among protozoan predators, raptorial Haptorian ciliates are particularly fascinating as they possess offensive extrusomes known as toxicysts, which are rapidly discharged upon prey contact. However, our understanding of the genetic processes and specific toxins involved in toxicyst formation and discharge is still limited.

RESULTS

In this study, we investigated the predation strategies and subcellular structures of seven Haptoria ciliate species and obtained their genome sequences using single-cell sequencing technology. Comparative genomic analysis revealed distinct gene duplications related to membrane transport proteins and hydrolytic enzymes in Haptoria, which play a crucial role in the production and discharge of toxicysts. Transcriptomic analysis further confirmed the abundant expression of genes related to membrane transporters and cellular toxins in Haptoria compared to Trichostomatia. Notably, polyketide synthases (PKS) and L-amino acid oxidases (LAAO) were identified as potentially toxin genes that underwent extensive duplication events in Haptoria.

CONCLUSIONS

Our results shed light on the evolutionary and genomic adaptations of Haptorian ciliates for their predation strategies in evolution and provide insights into their toxic mechanisms.

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.

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

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

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.

De novo transcriptome assembly and gene annotation for the toxic dinoflagellate Dinophysis

Species within the dinoflagellate genus Dinophysis can produce okadiac acid and dinophysistoxins leading to diarrhetic shellfish poisoning. Since the first report of D. ovum from the Gulf of Mexico in 2008, reports of other Dinophysis species across US have increased. Members of the D. cf. acuminata complex (D. acuminata, D. acuta, D. ovum, D. sacculus) are difficult to differentiate due to their morphological similarities. Dinophysis feeds on and steals the chloroplasts from the ciliate, Mesodinium rubrum, which in turn has fed on and captured the chloroplasts of its prey, the cryptophyte Teleaulax amphioxeia. The objective of this study was to generate de novo transcriptomes for new isolates of these mixotrophic organisms. The transcriptomes obtained will serve as a reference for future experiments to assess the effect of different abiotic and biotic conditions and will also provide a useful resource for screening potential marker genes to differentiate among the closely related species within the D. cf. acuminata-complex. The complete comprehensive detailed workflow and links to obtain the transcriptome data are provided.

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

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

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.

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.

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.

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.