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

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.

Integrative analyses on the ciliates Colpoda illuminate the life history evolution of soil microorganisms

Microorganisms play a central role in sustaining soil ecosystems and agriculture, and these functions are usually associated with their complex life history. Yet, the regulation and evolution of life history have remained enigmatic and poorly understood, especially in protozoa, the third most abundant group of organisms in the soil. Here, we explore the life history of a cosmopolitan species-Colpoda steinii. Our analysis has yielded a high-quality macronuclear genome for C. steinii, with size of 155 Mbp and 37,123 protein-coding genes, as well as mean intron length of ~93 bp, longer than most other studied ciliates. Notably, we identify two possible whole-genome duplication events in C. steinii, which may account for its genome being about twice the size of C. inflata's, another co-existing species. We further resolve the gene expression profiles in diverse life stages of C. steinii, which are also corroborated in C. inflata. During the resting cyst stage, genes associated with cell death and vacuole formation are upregulated, and translation-related genes are downregulated. While the translation-related genes are upregulated during the excystment of resting cysts. Reproductive cysts exhibit a significant reduction in cell adhesion. We also demonstrate that most genes expressed in specific life stages are under strong purifying selection. This study offers a deeper understanding of the life history evolution that underpins the extraordinary success and ecological functions of microorganisms in soil ecosystems.IMPORTANCEColpoda species, as a prominent group among the most widely distributed and abundant soil microorganisms, play a crucial role in sustaining soil ecosystems and promoting plant growth. This investigation reveals their exceptional macronuclear genomic features, including significantly large genome size, long introns, and numerous gene duplications. The gene expression profiles and the specific biological functions associated with the transitions between various life stages are also elucidated. The vast majority of genes linked to life stage transitions are subject to strong purifying selection, as inferred from multiple natural strains newly isolated and deeply sequenced. This substantiates the enduring and conservative nature of Colpoda's life history, which has persisted throughout the extensive evolutionary history of these highly successful protozoa in soil. These findings shed light on the evolutionary dynamics of microbial eukaryotes in the ever-fluctuating soil environments. This integrative research represents a significant advancement in understanding the life histories of these understudied single-celled eukaryotes.

Morphology and molecular phylogeny of three Parasonderia species including a new species (Ciliophora, Plagiopylea)

Ciliates of the class Plagiopylea play a vital role in various anaerobic environments as consumers of prokaryotes. Yet, the diversity and phylogeny of this group of ciliates, especially marine representatives, remain poorly known. In this study, three Parasonderia species, viz., Parasonderia elongata spec. nov., and the already known P. cyclostoma and P. vestita, discovered in anaerobic sediments from various intertidal zones in China, were investigated based on their living morphology, infraciliature, and small subunit ribosomal rRNA gene sequences. Parasonderia elongata can be recognized by its larger body size, elongated body shape, oval oral opening, number of oral kineties, and significantly shortened leftmost postbuccal polykineties on the cell surface. Improved diagnosis and redescription of P. cyclostoma is provided for the first time, including data on infraciliature and molecular sequence. Phylogenetic analyses revealed that the three species cluster together and with the sequence of a Chinese population of P. vestita already present in the GenBank database, forming a robust clade.

Exploring the landscape of symbiotic diversity and distribution in unicellular ciliated protists

BACKGROUND

The eukaryotic-bacterial symbiotic system plays an important role in various physiological, developmental, and evolutionary processes. However, our current understanding is largely limited to multicellular eukaryotes without adequate consideration of diverse unicellular protists, including ciliates.

RESULTS

To investigate the bacterial profiles associated with unicellular organisms, we collected 246 ciliate samples spanning the entire Ciliophora phylum and conducted single-cell based metagenome sequencing. This effort has yielded the most extensive collection of bacteria linked to unicellular protists to date. From this dataset, we identified 883 bacterial species capable of cohabiting with ciliates, unveiling the genomes of 116 novel bacterial cohabitants along with 7 novel archaeal cohabitants. Highlighting the intimate relationship between ciliates and their cohabitants, our study unveiled that over 90% of ciliates coexist with bacteria, with individual hosts fostering symbiotic relationships with multiple bacteria concurrently, resulting in the observation of seven distinct symbiotic patterns among bacteria. Our exploration of symbiotic mechanisms revealed the impact of host digestion on the intracellular diversity of cohabitants. Additionally, we identified the presence of eukaryotic-like proteins in bacteria as a potential contributing factor to their resistance against host digestion, thereby expanding their potential host range.

CONCLUSIONS

As the first large-scale analysis of prokaryotic associations with ciliate protists, this study provides a valuable resource for future research on eukaryotic-bacterial symbioses. Video Abstract.

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.

Morphogenesis of an anaerobic ciliate Heterometopus palaeformis (Kahl, 1927) Foissner, 2016 (Ciliophora, Armophorea) with notes on its morphological and molecular characterization

The morphology, morphogenesis, and molecular phylogeny of Heterometopus palaeformis (Kahl, 1927) Foissner, 2016 were studied using microscopical observations on live and protargol-stained specimens as well SSU rRNA gene sequencing. The morphogenetic data for the genus are presented for the first time. Compared to other metopids, the morphogenesis of H. palaeformis is distinct since its (1) perizonal stripe rows 4 and 5 are involved in the formation of the opisthe's adoral polykinetids; (2) perizonal stripe rows 3-5 and two adjacent preoral dome kineties contribute to most of the opisthe's paroral membrane while perizonal stripe rows 1 and 2 contribute very little; (3) four kinety rows are formed to the left of the opisthe's adoral zone of polykinetids. The Chinese population resembles the original and neotype populations well in terms of general morphology - characterized by a life size of 55-120 × 10-20 μm, an elongate ellipsoidal body with a hardly spiralized flat preoral dome, about 18 somatic kineties and 20 adoral polykinetids. The SSU rDNA sequence of the present population exhibits a disparity of 1.33%-2.22% divergence from sequences of other populations. Nevertheless, phylogenetic analysis reveals that populations of H. palaeformis form a separate, stable cluster within the paraphyletic Metopidae clade.

Divergent marine anaerobic ciliates harbor closely related Methanocorpusculum endosymbionts

Ciliates are a diverse group of protists known for their ability to establish various partnerships and thrive in a wide variety of oxygen-depleted environments. Most anaerobic ciliates harbor methanogens, one of the few known archaea living intracellularly. These methanogens increase the metabolic efficiency of host fermentation via syntrophic use of host end-product in methanogenesis. Despite the ubiquity of these symbioses in anoxic habitats, patterns of symbiont specificity and fidelity are not well known. We surveyed two unrelated, commonly found groups of anaerobic ciliates, the Plagiopylea and Metopida, isolated from anoxic marine sediments. We sequenced host 18S rRNA and symbiont 16S rRNA marker genes as well as the symbiont internal transcribed spacer region from our cultured ciliates to identify hosts and their associated methanogenic symbionts. We found that marine ciliates from both of these co-occurring, divergent groups harbor closely related yet distinct intracellular archaea within the Methanocorpusculum genus. The symbionts appear to be stable at the host species level, but at higher taxonomic levels, there is evidence that symbiont replacements have occurred. Gaining insight into this unique association will deepen our understanding of the complex transmission modes of marine microbial symbionts, and the mutualistic microbial interactions occurring across domains of life.

The energy metabolism of Balantidium polyvacuolum inhabiting the hindgut of Xenocypris davidi

Anaerobic parasitic ciliates are a specialized group of ciliates that are adapted to anoxic and oxygen-depleted habitats. Among them, Balantidium polyvacuolum, which inhabits the hindgut of Xenocyprinae fishes, has received very limited scientific attention, so the molecular mechanism of its adaptation to the digestive tract microenvironment is still unclear. In this study, transmission electron microscopy (TEM) and single-cell transcriptome analysis were used to uncover the metabolism of B. polyvacuolum. Starch granules, endosymbiotic bacteria, and multiple specialized mitochondrion-related organelles (MROs) of various shapes were observed. The MROs may have completely lost the electron transport chain (ETC) complexes I, III, IV, and V and only retained succinate dehydrogenase subunit A (SDHA) of complex II. The tricarboxylic acid (TCA) cycle was also incomplete. It can be inferred that the hypoxic intestinal environment has led to the specialization of the mitochondria in B. polyvacuolum. Moreover, carbohydrate-active enzymes (CAZymes), including carbohydrate esterases, enzymes with a carbohydrate-binding module, glycoside hydrolases, and glycosyltransferases, were identified, which may constitute evidence that B. polyvacuolum is able to digest carbohydrates and starch. These findings can improve our knowledge of the energy metabolism and adaptive mechanisms of B. polyvacuolum.

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.

Extreme environments offer an unprecedented opportunity to understand microbial eukaryotic ecology, evolution, and genome biology

Research in extreme environments has substantially expanded our understanding of the ecology and evolution of life on Earth, but a major group of organisms has been largely overlooked: microbial eukaryotes (i.e., protists). In this Perspective, we summarize data from over 80 studies of protists in extreme environments and identify focal lineages that are of significant interest for further study, including clades within Echinamoebida, Heterolobosea, Radiolaria, Haptophyta, Oomycota, and Cryptophyta. We argue that extreme environments are prime sampling targets to fill gaps in the eukaryotic tree of life and to increase our understanding of the ecology, metabolism, genome architecture, and evolution of eukaryotic life.

Anaerocyclidiidae fam. nov. (Oligohymenophorea, Scuticociliatia): A newly recognized major lineage of anaerobic ciliates hosting prokaryotic symbionts

The research on anaerobic ciliates, to date, has mainly been focused on representatives of obligately anaerobic classes such as Armophorea or Plagiopylea. In this study, we focus on the anaerobic representatives of the subclass Scuticociliatia, members of the class Oligohymenophorea, which is mainly composed of aerobic ciliates. Until now, only a single anaerobic species, Cyclidium porcatum (here transferred to the genus Anaerocyclidium gen. nov.), has been described both molecularly and morphologically. Our broad sampling of anoxic sediments together with cultivation and single cell sequencing approaches have shown that scuticociliates are common and diversified in anoxic environments. Our results show that anaerobic scuticociliates represent a distinctive evolutionary lineage not closely related to the family Cyclidiidae (order Pleuronematida), as previously suggested. However, the phylogenetic position of the newly recognized lineage within the subclass Scuticociliatia remains unresolved. Based on molecular and morphological data, we establish the family Anaerocyclidiidae fam. nov. to accommodate members of this clade. We further provide detailed morphological descriptions and 18S rRNA gene sequences for six new Anaerocyclidium species and significantly broaden the described diversity of anaerobic scuticociliates.

Culture-Independent Single-Cell PacBio Sequencing Reveals Epibiotic Variovorax and Nucleus Associated Mycoplasma in the Microbiome of the Marine Benthic Protist Geleia sp. YT (Ciliophora, Karyorelictea)

Microbes in marine sediments constitute up to five-sixths of the planet's total biomass, but their diversity is little explored, especially for those forming associations with unicellular protists. Heterotrophic ciliates are among the most dominant and diversified marine benthic protists and comprise hotspot niches of bacterial colonization. To date, studies using culture-independent single-cell approaches to explore microbiomes of marine benthic ciliates in nature are almost absent, even for the most ubiquitous species. Here, we characterize the major bacterial groups associated with a representative marine benthic ciliate, Geleia sp. YT, collected directly from the coastal zone of Yantai, China. PacBio sequencing of the nearly full-length 16Sr RNA genes was performed on single cells of Geleia. Fluorescence in situ hybridization (FISH) analysis with genus-specific probes was further applied to locate the dominant bacterial groups. We identified a Variovorax-like bacterium as the major epibiotic symbiont residing in the kineties of the ciliate host. We provide evidence of a nucleus-associated bacterium related to the human pathogen Mycoplasma, which appeared prevalently in the local populations of Geleia sp. YT for 4 months. The most abundant bacterial taxa associated with Geleia sp. YT likely represent its core microbiome, hinting at the important roles of the ciliate-bacteria consortium in the marine benthos. Overall, this work has contributed to the knowledge of the diversity of life in the enigmatic marine benthic ciliate and its symbioses.

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.

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

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

Mitochondrial metabolism of the facultative parasite Chilodonella uncinata (Alveolata, Ciliophora)

BACKGROUND

Chilodonella uncinata is an aerobic ciliate capable of switching between being free-living and parasitic on fish fins and gills, causing tissue damage and host mortality. It is widely used as a model organism for genetic studies, but its mitochondrial metabolism has never been studied. Therefore, we aimed to describe the morphological features and metabolic characteristics of its mitochondria.

METHODS

Fluorescence staining and transmission electron microscopy (TEM) were used to observe the morphology of mitochondria. Single-cell transcriptome data of C. uncinata were annotated by the Clusters of Orthologous Genes (COG) database. Meanwhile, the metabolic pathways were constructed based on the transcriptomes. The phylogenetic analysis was also made based on the sequenced cytochrome c oxidase subunit 1 (COX1) gene.

RESULTS

Mitochondria were stained red using Mito-tracker Red staining and were stained slightly blue by DAPI dye. The cristae and double membrane structures of the mitochondria were observed by TEM. Besides, many lipid droplets were evenly distributed around the macronucleus. A total of 2594 unigenes were assigned to 23 functional classifications of COG. Mitochondrial metabolic pathways were depicted. The mitochondria contained enzymes for the complete tricarboxylic acid (TCA) cycle, fatty acid metabolism, amino acid metabolism, and cytochrome-based electron transport chain (ETC), but only partial enzymes involved in the iron-sulfur clusters (ISCs).

CONCLUSIONS

Our results showed that C. uncinata possess typical mitochondria. Stored lipid droplets inside mitochondria may be the energy storage of C. uncinata that helps its transmission from a free-living to a parasitic lifestyle. These findings also have improved our knowledge of the mitochondrial metabolism of C. uncinata and increased the volume of molecular data for future studies of this facultative parasite.

Morphology and phylogenetic position of three anaerobic ciliates from the classes Odontostomatea and Muranotrichea (Ciliophora)

The diversity of the classes Odontostomatea and Muranotrichea, which contain solely obligate anaerobes, is poorly understood. We studied two populations of Mylestoma sp., one of Saprodinium dentatum (Odontostomatea), two of Muranothrix felix sp. nov., and one of Muranothrix sp. (Muranotrichea) employing live observation, protargol impregnation, scanning electron microscopy, and 18S rRNA gene sequencing. Conspecificity of Mylestoma sp., described here, with a previously described species of this genus cannot be excluded since no species have been studied with modern methods. Phylogenetically, the genus Mylestoma is closely related to the odontostomatid Discomorphella pedroeneasi, although the phylogenetic position of class Odontostomatea itself remains unresolved. The newly described muranotrichean species, Muranothrix felix sp. nov., is morphologically similar to M. gubernata but can be distinguished by its fewer macronuclear nodules and fewer adoral membranelles; moreover, it is clearly distinguished from M. gubernata by its 18S rRNA gene sequence. Another population, designated here as Muranothrix sp., most likely represents a separate species.

Energetics and evolution of anaerobic microbial eukaryotes

Mitochondria and aerobic respiration have been suggested to be required for the evolution of eukaryotic cell complexity. Aerobic respiration is several times more energetically efficient than fermentation. Moreover, aerobic respiration occurs at internalized mitochondrial membranes that are not constrained by a sublinear scaling with cell volume. However, diverse and complex anaerobic eukaryotes (for example, free-living and parasitic unicellular, and even small multicellular, eukaryotes) that exclusively rely on fermentation for energy generation have evolved repeatedly from aerobic ancestors. How do fermenting eukaryotes maintain their cell volumes and complexity while relying on such a low energy-yielding process? Here I propose that reduced rates of ATP generation in fermenting versus respiring eukaryotes are compensated for by longer cell cycles that satisfy lifetime energy demands. A literature survey and growth efficiency calculations show that fermenting eukaryotes divide approximately four to six times slower than aerobically respiring counterparts with similar cell volumes. Although ecological advantages such as competition avoidance offset lower growth rates and yields in the short term, fermenting eukaryotes inevitably have fewer physiological and ecological possibilities, which ultimately constrain their long-term evolutionary trajectories.

Eukaryotic evolution: Spatial proteomics sheds light on mitochondrial reduction

Multi-organelle spatial proteomics has revolutionized animal cell biology, but its use in protists has so far been limited. A new study delivers the first such proteome of a free-living protist, uncovering a previously overlooked function of highly reduced mitochondria.

Reduced mitochondria provide an essential function for the cytosolic methionine cycle

The loss of mitochondria in oxymonad protists has been associated with the redirection of the essential Fe-S cluster assembly to the cytosol. Yet as our knowledge of diverse free-living protists broadens, the list of functions of their mitochondrial-related organelles (MROs) expands. We revealed another such function in the closest oxymonad relative, Paratrimastix pyriformis, after we solved the proteome of its MRO with high accuracy, using localization of organelle proteins by isotope tagging (LOPIT). The newly assigned enzymes connect to the glycine cleavage system (GCS) and produce folate derivatives with one-carbon units and formate. These are likely to be used by the cytosolic methionine cycle involved in S-adenosyl methionine recycling. The data provide consistency with the presence of the GCS in MROs of free-living species and its absence in most endobionts, which typically lose the methionine cycle and, in the case of oxymonads, the mitochondria.

The Genome of the Mitochondrion-Related Organelle in Cepedea longa, a Large Endosymbiotic Opalinid Inhabiting the Recta of Frogs

Mitochondrion-related organelles (MROs) are loosely defined as degenerated mitochondria in anaerobic and microaerophilic lineages. Opalinids are commonly regarded as commensals in the guts of cold-blooded amphibians. It may represent an intermediate adaptation stage between the conventional aerobic mitochondria and derived anaerobic MROs. In the present study, we sequenced and analyzed the MRO genome of Cepedea longa. It has a linear MRO genome with large inverted repeat gene regions at both ends. Compared to Blastocystis and Proteromonas lacertae, the MRO genome of C. longa has a higher G + C content and repeat sequences near the central region. Although three Opalinata species have different morphological characteristics, phylogenetic analyses based on eight concatenated nad genes indicate that they are close relatives. The phylogenetic analysis showed that C. longa clustered with P. lacertae with strong support. The 18S rRNA gene-based phylogeny resolved the Opalinea clade as a sister clade to Karotomorpha, which then further grouped with Proteromonas. The paraphyly of Proteromonadea needs to be verified due to the lack of MRO genomes for key species, such as Karotomorpha, Opalina and Protoopalina. Besides, our dataset and analyses offered slight support for the paraphyly of Bigyra.

Group-specific functional patterns of mitochondrion-related organelles shed light on their multiple transitions from mitochondria in ciliated protists

Adaptations of ciliates to hypoxic environments have arisen independently several times. Studies on mitochondrion-related organelle (MRO) metabolisms from distinct anaerobic ciliate groups provide evidence for understanding the transitions from mitochondria to MROs within eukaryotes. To deepen our knowledge about the evolutionary patterns of ciliate anaerobiosis, mass-culture and single-cell transcriptomes of two anaerobic species, Metopus laminarius (class Armophorea) and Plagiopyla cf. narasimhamurtii (class Plagiopylea), were sequenced and their MRO metabolic maps were compared. In addition, we carried out comparisons using publicly available predicted MRO proteomes from other ciliate classes (i.e., Armophorea, Litostomatea, Muranotrichea, Oligohymenophorea, Parablepharismea and Plagiopylea). We found that single-cell transcriptomes were similarly comparable to their mass-culture counterparts in predicting MRO metabolic pathways of ciliates. The patterns of the components of the MRO metabolic pathways might be divergent among anaerobic ciliates, even among closely related species. Notably, our findings indicate the existence of group-specific functional relics of electron transport chains (ETCs). Detailed group-specific ETC functional patterns are as follows: full oxidative phosphorylation in Oligohymenophorea and Muranotrichea; only electron-transfer machinery in Armophorea; either of these functional types in Parablepharismea; and ETC functional absence in Litostomatea and Plagiopylea. These findings suggest that adaptation of ciliates to anaerobic conditions is group-specific and has occurred multiple times. Our results also show the potential and the limitations of detecting ciliate MRO proteins using single-cell transcriptomes and improve the understanding of the multiple transitions from mitochondria to MROs within ciliates.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-022-00147-w.

Phylogeny of the anaerobic ciliate genus Sonderia (Protista: Ciliophora: Plagiopylea), including the description of three novel species and a brief revision of the genus

Anaerobic protists in general, and ciliates in particular, are important components of anoxic or hypoxic environments, however, their diversity remains underestimated. Sonderia is a poorly studied genus that is distributed worldwide and is commonly found in anaerobic environments. In the present study, the taxonomy and phylogeny of three new species, namely Sonderia aposinuata sp. nov., Sonderia paramacrochilus sp. nov. and Sonderia steini sp. nov., collected from China, were investigated based on microscopic observations and SSU rRNA gene sequencing methods. Sonderia aposinuata sp. nov. is diagnosed mainly by having a relatively large body size, a crescent-shaped oral opening, numerous slender extrusomes, one suture on the ventral side and two on the dorsal side, and a buccal cavity that occupies the anterior third of the cell. Sonderia paramacrochilus sp. nov. closely resembles S. macrochilus but differs mainly by its oral opening being located closer to the anterior cell margin and its spindle-shaped extrusomes. Sonderia steini sp. nov. is a freshwater species that can be recognized by its shallow buccal cavity, sparsely distributed rod-shaped extrusomes, and having 68-79 monokinetidal somatic kineties that form sutures on both sides of the body. Phylogenetic analyses based on small subunit ribosomal RNA (SSU rRNA) gene sequence data support the monophyly of the family Sonderiidae, however, Sonderia is paraphyletic. The genus Sonderia is briefly revised and a key to the identification of species belonging to this genus is supplied.

Evidence for an Independent Hydrogenosome-to-Mitosome Transition in the CL3 Lineage of Fornicates

Fornicata, a lineage of a broader and ancient anaerobic eukaryotic clade Metamonada, contains diverse taxa that are ideally suited for evolutionary studies addressing various fundamental biological questions, such as the evolutionary trajectory of mitochondrion-related organelles (MROs), the transition between free-living and endobiotic lifestyles, and the derivation of alternative genetic codes. To this end, we conducted detailed microscopic and transcriptome analyses in a poorly documented strain of an anaerobic free-living marine flagellate, PCS, in the so-called CL3 fornicate lineage. Fortuitously, we discovered that the original culture contained two morphologically similar and closely related CL3 representatives, which doubles the taxon representation within this lineage. We obtained a monoeukaryotic culture of one of them and formally describe it as a new member of the family Caviomonadidae, Euthynema mutabile gen. et sp. nov. In contrast to previously studied caviomonads, the endobiotic Caviomonas mobilis and Iotanema spirale, E. mutabile possesses an ultrastructurally discernible MRO. We sequenced and assembled the transcriptome of E. mutabile, and by sequence subtraction, obtained transcriptome data from the other CL3 clade representative present in the original PCS culture, denoted PCS-ghost. Transcriptome analyses showed that the reassignment of only one of the UAR stop codons to encode Gln previously reported from I. spirale does not extend to its free-living relatives and is likely due to a unique amino acid substitution in I. spirale’s eRF1 protein domain responsible for termination codon recognition. The backbone fornicate phylogeny was robustly resolved in a phylogenomic analysis, with the CL3 clade amongst the earliest branching lineages. Metabolic and MRO functional reconstructions of CL3 clade members revealed that all three, including I. spirale, encode homologs of key components of the mitochondrial protein import apparatus and the ISC pathway, indicating the presence of a MRO in all of them. In silico evidence indicates that the organelles of E. mutabile and PCS-ghost host ATP and H₂ production, unlike the cryptic MRO of I. spirale. These data suggest that the CL3 clade has experienced a hydrogenosome-to-mitosome transition independent from that previously documented for the lineage leading to Giardia.

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.

Morphology and phylogeny of two anaerobic freshwater ciliates: Brachonella comma sp. nov. and the widely-distributed but little-known caenomorphid, Ludio parvulus Penard, 1922

Hypoxic, sulfidic freshwater sediments typically support a diffuse consortium of distinctive ciliated protists, including caenomorphids, metopids, and odontostomatids among others. A recent resurgence of interest in these important members of sapropelic food webs has resulted in the description of many new species and an effort, still in its infancy, to characterize them from a morphologic, molecular, and metabolic standpoint and to determine their phylogenetic relationships. Their seemingly invariable association with prokaryotic endosymbionts and, less commonly, ectosymbionts, has become a focus for many researchers. In this report, based on morphologic and molecular data we describe a Brachonella species (Ciliophora, Metopida) new to science and analyze its phylogeny. We also provide a morphologic and molecular characterization of the smallest representative of the Caenomorphidae Poche, 1913, Ludio parvulus Penard, 1922. The phylogenetic analysis confirms the inclusion of this species in the Caenomorphidae.

Orthology Prediction and Phylogenetic Analysis Methods in Plants

In this chapter, we outline a pipeline for ortholog prediction and phylogenetic analysis in plants. This computational pipeline uses algorithms from different software to enable bioinformatic-beginner biologists to predict orthologs that can be shared with many distinct plant nonmodel and model species and identify gene loss events. Prediction of orthologs allows (1) investigation of the evolutionary relationships of plant genomes, (2) discovery of their origin, function, and (3) the impact of their adaptability to the environment.We developed a pipeline to fit, not only eukaryote but also prokaryote organisms, with small or large genomes. All results acquired from the orthologs predication will enable phylogenetic tree construction, using gene and species (phylogenomic) phylogeny approaches.

A putatively extinct higher taxon of Spirotrichea (Ciliophora) from the Lower Cretaceous of Brazil

Fossil microeukaryotes are key elements for understanding ancient ecosystems at microscopic level and improving the knowledge on the diversification of microbial life as a whole. We describe Palaeohypothrix bahiensis gen. et sp. nov., an exceptionally well-preserved Lower Cretaceous (Berriasian-Barremian; 145-125 Mya) amber-entrapped microeukaryote, identified as a spirotrich ciliate. The preservation of structures interpreted as the nuclear apparatus and remains of the ciliature revealed a novel ground plan, not found in modern Spirotrichea, thus representing a putatively extinct higher taxon lineage, viz. the Palaeohypotricha nov. tax. Based on cladistic analysis, the new taxon is hypothesized as phylogenetically related to the Protohypotrichia.

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.

Morphology and molecular phylogeny of the anaerobic freshwater ciliate Urostomides spinosus nov. spec. (Ciliophora, Armophorea, Metopida) from China

The morphology and molecular phylogeny of a new metopid ciliate, Urostomides spinosus nov. spec., discovered in a freshwater ditch in Qingdao, China, were investigated using live observation, morphometry and protargol staining as well as molecular phylogenetic methods. Diagnostic features of the new species include a broadly obpyriform body carrying three posterior spines, eight somatic kineties, five preoral dome kineties with specialized row 3, adoral zone composed of about 28 membranelles, making a 270° turn around body axis. Phylogenetic analyses of the SSU rDNA sequence revealed that the genus Urostomides is monophyletic, but its interspecific relationships remained unresolved. Moreover, a closer relationship of the new species with the morphologically similar Urostomides campanula was not supported by the molecular data.

Macroevolutionary analyses of ciliates associated with hosts support high diversification rates

Ciliophora is a phylum that is comprised of extremely diverse microorganisms with regard to their morphology and ecology. They may be found in various environments, as free-living organisms or associated with metazoans. Such associations range from relationships with low metabolic dependence such as epibiosis, to more intimate relationships such as mutualism and parasitism. We know that symbiotic relationships occur along the whole phylogeny of the group, however, little is known about their evolution. Theoretical studies show that there are two routes for the development of parasitism, yet few authors have investigated the evolution of these characteristics using molecular tools. In the present study, we inferred a wide dated molecular phylogeny, based on the 18S rDNA gene, for the entire Ciliophora phylum, mapped life habits throughout the evolutionary time, and evaluated whether symbiotic relationships were linked to the variation in diversification rates and to the mode of evolution of ciliates. Our results showed that the last common ancestor for Ciliophora was likely a free-living organism, and that parasitism is a recent adaptation in ciliates, emerging more than once and independently via two distinct routes: (i) a free-living ciliate evolved into a mutualistic organism and, later, into a parasitic organism, and (ii) a free-living ciliate evolved directly into a parasitic organism. Furthermore, we have found a significant increase in the diversification rate of parasitic and mutualistic ciliates compared with their free-living conspecifics. The evolutionary success in different lineages of symbiont ciliates may be associated with many factors including type and colonization placement on their host, as well as physical and physiological conditions made available by the hosts.

Parallel functional reduction in the mitochondria of apicomplexan parasites

Gregarines are an early-diverging lineage of apicomplexan parasites that hold many clues into the origin and evolution of the group, a remarkable transition from free-living phototrophic algae into obligate parasites of animals.¹ Using single-cell transcriptomics targeting understudied lineages to complement available sequencing data, we characterized the mitochondrial metabolic repertoire across the tree of apicomplexans. In contrast to the large suite of proteins involved in aerobic respiration in well-studied parasites like Toxoplasma or Plasmodium,² we find that gregarine trophozoites have significantly reduced energy metabolism: most lack respiratory complexes III and IV, and some lack the electron transport chains (ETCs) and tricarboxylic acid (TCA) cycle entirely. Phylogenomic analyses show that these reductions took place several times in parallel, resulting in a functional range from fully aerobic organelles to extremely reduced "mitosomes" restricted to Fe-S cluster biosynthesis. The mitochondrial genome has also been lost repeatedly: in species with severe functional reduction simply by gene loss but in one species with a complete ETC by relocating cox1 to the nuclear genome. Severe functional reduction of mitochondria is generally associated with structural reduction, resulting in small, nondescript mitochondrial-related organelles (MROs).³ By contrast, gregarines retain distinctive mitochondria with tubular cristae, even the most functionally reduced cases that also lack genes associated with cristae formation. Overall, the parallel, severe reduction of gregarine mitochondria expands the diversity of organisms that contain MROs and further emphasizes the role of parallel transitions in apicomplexan evolution.

Gregarine single-cell transcriptomics reveals differential mitochondrial remodeling and adaptation in apicomplexans

BACKGROUND

Apicomplexa is a diverse phylum comprising unicellular endobiotic animal parasites and contains some of the most well-studied microbial eukaryotes including the devastating human pathogens Plasmodium falciparum and Cryptosporidium hominis. In contrast, data on the invertebrate-infecting gregarines remains sparse and their evolutionary relationship to other apicomplexans remains obscure. Most apicomplexans retain a highly modified plastid, while their mitochondria remain metabolically conserved. Cryptosporidium spp. inhabit an anaerobic host-gut environment and represent the known exception, having completely lost their plastid while retaining an extremely reduced mitochondrion that has lost its genome. Recent advances in single-cell sequencing have enabled the first broad genome-scale explorations of gregarines, providing evidence of differential plastid retention throughout the group. However, little is known about the retention and metabolic capacity of gregarine mitochondria.

RESULTS

Here, we sequenced transcriptomes from five species of gregarines isolated from cockroaches. We combined these data with those from other apicomplexans, performed detailed phylogenomic analyses, and characterized their mitochondrial metabolism. Our results support the placement of Cryptosporidium as the earliest diverging lineage of apicomplexans, which impacts our interpretation of evolutionary events within the phylum. By mapping in silico predictions of core mitochondrial pathways onto our phylogeny, we identified convergently reduced mitochondria. These data show that the electron transport chain has been independently lost three times across the phylum, twice within gregarines.

CONCLUSIONS

Apicomplexan lineages show variable functional restructuring of mitochondrial metabolism that appears to have been driven by adaptations to parasitism and anaerobiosis. Our findings indicate that apicomplexans are rife with convergent adaptations, with shared features including morphology, energy metabolism, and intracellularity.

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

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

Molecular characterization of the conoid complex in Toxoplasma reveals its conservation in all apicomplexans, including Plasmodium species

The apical complex is the instrument of invasion used by apicomplexan parasites, and the conoid is a conspicuous feature of this apparatus found throughout this phylum. The conoid, however, is believed to be heavily reduced or missing from Plasmodium species and other members of the class Aconoidasida. Relatively few conoid proteins have previously been identified, making it difficult to address how conserved this feature is throughout the phylum, and whether it is genuinely missing from some major groups. Moreover, parasites such as Plasmodium species cycle through 3 invasive forms, and there is the possibility of differential presence of the conoid between these stages. We have applied spatial proteomics and high-resolution microscopy to develop a more complete molecular inventory and understanding of the organisation of conoid-associated proteins in the model apicomplexan Toxoplasma gondii. These data revealed molecular conservation of all conoid substructures throughout Apicomplexa, including Plasmodium, and even in allied Myzozoa such as Chromera and dinoflagellates. We reporter-tagged and observed the expression and location of several conoid complex proteins in the malaria model P. berghei and revealed equivalent structures in all of its zoite forms, as well as evidence of molecular differentiation between blood-stage merozoites and the ookinetes and sporozoites of the mosquito vector. Collectively, we show that the conoid is a conserved apicomplexan element at the heart of the invasion mechanisms of these highly successful and often devastating parasites.

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.

Diversity of electron transport chains in anaerobic protists

Eukaryotic microbes (protists) that occupy low-oxygen environments often have drastically different mitochondrial metabolism compared to their aerobic relatives. A common theme among many anaerobic protists is the serial loss of components of the electron transport chain (ETC). Here, we discuss the diversity of the ETC across the tree of eukaryotes and review hypotheses for how ETCs are modified, and ultimately lost, in protists. We find that while protists have converged to some of the same metabolism as anaerobic animals, there are clear protist-specific strategies to thrive without oxygen.

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.

Anaerobic Fungi: Past, Present, and Future

Anaerobic fungi (AF) play an essential role in feed conversion due to their potent fiber degrading enzymes and invasive growth. Much has been learned about this unusual fungal phylum since the paradigm shifting work of Colin Orpin in the 1970s, when he characterized the first AF. Molecular approaches targeting specific phylogenetic marker genes have facilitated taxonomic classification of AF, which had been previously been complicated by the complex life cycles and associated morphologies. Although we now have a much better understanding of their diversity, it is believed that there are still numerous genera of AF that remain to be described in gut ecosystems. Recent marker-gene based studies have shown that fungal diversity in the herbivore gut is much like the bacterial population, driven by host phylogeny, host genetics and diet. Since AF are major contributors to the degradation of plant material ingested by the host animal, it is understandable that there has been great interest in exploring the enzymatic repertoire of these microorganisms in order to establish a better understanding of how AF, and their enzymes, can be used to improve host health and performance, while simultaneously reducing the ecological footprint of the livestock industry. A detailed understanding of AF and their interaction with other gut microbes as well as the host animal is essential, especially when production of affordable high-quality protein and other animal-based products needs to meet the demands of an increasing human population. Such a mechanistic understanding, leading to more sustainable livestock practices, will be possible with recently developed -omics technologies that have already provided first insights into the different contributions of the fungal and bacterial population in the rumen during plant cell wall hydrolysis.