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

Complex evolutionary patterns within the tubulin gene family of ciliates, unicellular eukaryotes with diverse microtubular structures

BACKGROUND

Tubulins are major components of the eukaryotic cytoskeletons that are crucial in many cellular processes. Ciliated protists comprise one of the oldest eukaryotic lineages possessing cilia over their cell surface and assembling many diverse microtubular structures. As such, ciliates are excellent model organisms to clarify the origin and evolution of tubulins in the early stages of eukaryote evolution. Nonetheless, the evolutionary history of the tubulin subfamilies within and among ciliate classes is unclear.

RESULTS

We analyzed the evolutionary pattern of ciliate tubulin gene family based on genomes/transcriptomes of 60 species covering 10 ciliate classes. Results showed: (1) Six tubulin subfamilies (α_Tub, β_Tub, γ_Tub, δ_Tub, ε_Tub, and ζ_Tub) originated from the last eukaryotic common ancestor (LECA) were observed within ciliates. Among them, α_Tub, β_Tub, and γ_Tub were present in all ciliate species, while δ_Tub, ε_Tub, and ζ_Tub might be independently lost in some species. (2) The evolutionary history of the tubulin subfamilies varied. Evolutionary history of ciliate γ_Tub, δ_Tub, ε_Tub, and ζ_Tub showed a certain degree of consistency with the phylogeny of species after the divergence of ciliate classes, while the evolutionary history of ciliate α_Tub and β_Tub varied among different classes. (3) Ciliate α- and β-tubulin isoforms could be classified into an "ancestral group" present in LECA and a "divergent group" containing only ciliate sequences. Alveolata-specific expansion events probably occurred within the "ancestral group" of α_Tub and β_Tub. The "divergent group" might be important for ciliate morphological differentiation and wide environmental adaptability. (4) Expansion events of the tubulin gene family appeared to be consistent with whole genome duplication (WGD) events in some degree. More Paramecium-specific tubulin expansions were detected than Tetrahymena-specific ones. Compared to other Paramecium species, the Paramecium aurelia complex underwent a more recent WGD which might have experienced more tubulin expansion events.

CONCLUSIONS

Evolutionary history among different tubulin gene subfamilies seemed to vary within ciliated protists. And the complex evolutionary patterns of tubulins among different ciliate classes might drive functional diversification. Our investigation provided meaningful information for understanding the evolution of tubulin gene family in the early stages of eukaryote evolution.

Exploring the biogeography, morphology, and phylogeny of the condylostomatid ciliates (Alveolata, Ciliophora, Heterotrichea), with establishment of four new Condylostoma species and a revision including redescriptions of five species found in China

Species of the ciliate class Heterotrichea Stein, 1859 are a cosmopolitan group of unicellular eukaryotic microorganisms, many of which have been widely used as models in various fields of research such as regenerative biology, functional ecology, environmental toxicology, and symbiotic behavior. However, species identification in the heterotrich family Condylostomatidae, especially the most species-rich and type genus Condylostoma Bory de Saint-Vincent, 1824, remains challenging due to incomplete original descriptions, few reliable distinguishing characters, and overlapping features between different species. This study presents an updated revision of Condylostoma and its related genus Condylostomides da Silva Neto, 1994 based on descriptions of five species, including nine populations collected from China, using both morphological and molecular methods. The main findings are as follows: (1) 43 nominal species and about 130 populations are reviewed, resulting in the recognition of 30 valid species of Condylostoma and eight valid species of Condylostomides; (2) keys, synonyms, biogeographic distributions and amended/improved diagnoses of all valid species are provided; (3) based on the available data, four new Condylostoma species (C. marinum sp. nov., C. petzi sp. nov., C. villeneuvei sp. nov., and C. microstomum sp. nov.), one new combination (Condylostomides minimus (Dragesco, 1954) comb. nov. & nom. corr.), and two corrected names (Condylostoma ancestrale Villeneuve-Brachon, 1940 nom. corr. and Condylostomides nigrus (Dragesco, 1960) nom. corr.) are suggested; (4) cryptic species are detected and proposed for the first time to form the Condylostoma curvum species complex; (5) three highly confusing Condylostoma species, C. kris, C. spatiosum, and C. minutum, are redefined for the first time based on modern taxonomic methods; (6) a 'flagship' species, Condylostomides coeruleus, is recorded for the first time from the continent of Asia, substantially expanding its biogeography; (7) ciliature adjacent to the distal end of the paroral membrane within the family Condylostomatidae is uniformly defined as frontal membranelles and is classified into three patterns according to the arrangement of kinetosomes, which serve as important key features.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00223-3.

Revealing gene expression heterogeneity in a clonal population of Tetrahymena thermophila through single-cell RNA sequencing

We performed single-cell RNA sequencing (scRNA-seq) on a population of 5,000 Tetrahymena thermophila, using the 10x Genomics 3' gene expression analysis, to investigate gene expression variability within this clonal population. Initially, we estimated the 3'-untranslated regions (3' UTRs), which were absent in existing annotation files but are crucial for the 10x Genomics 3' gene expression analysis, using the peaks2utr method. This allowed us to create a modified annotation file, which was then utilized in our scRNA-seq analysis. Our analysis revealed significant gene expression variability within the population, even after removing the effect of cell phase-related features. This variability predominantly appeared in six distinct clusters. Through gene ontology and KEGG pathway enrichment analyses, we identified that these were primarily associated with ribosomal proteins, proteins specific to mitochondria, proteins involved in peroxisome-specific carbon metabolism, cytoskeletal proteins, motor proteins, and immobilized antigens.

Single-cell transcriptomic analysis reveals genome evolution in predatory litostomatean ciliates

Many ciliated protists prey on other large microbial organisms, including other protists and microscopic metazoans. The ciliate class Litostomatea unites both predatory and endosymbiotic species. The evolution of predation ability in ciliates remains poorly understood, in part, due to a lack of genomic data. To fill this gap, we acquired the transcriptome profiles of six predatory litostomateans using single-cell sequencing technology and investigated their transcriptomic features. Our results show that: (1) in contrast to non-predatory ciliates, the predatory litostomateans have expanded gene families associated with transmembrane activity and reactive oxidative stress response pathways, potentially as a result of cellular behaviors such as fast contraction and extension; (2) the expansion of the calcium-activated BK potassium channel gene family, which hypothetically regulates cell contractility, is an ancient evolutionary event for the class Litostomatea, suggesting a rewired metabolism associated with the hunting behavior of predatory ciliates; and (3) three whole genome duplication (WGD) events have been detected in litostomateans, with genes associated with biosynthetic processes, transmembrane activity, and calcium-activated potassium channel activity being retained during the WGD events. In addition, we explored the evolutionary relationships among 17 ciliate species, including eight litostomateans, and provided a rich foundational dataset for future in-depth phylogenomic studies of Litostomatea. Our comprehensive analyses suggest that the rewired cellular metabolism via expanded gene families and WGD events might be the potential genetic basis for the predation ability of raptorial ciliates.

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.

Single-cell transcriptomics supports presence of cryptic species and reveals low levels of population genetic diversity in two testate amoebae morphospecies with large population sizes

The enormous population sizes and wide biogeographical distribution of many microbial eukaryotes set the expectation of high levels of intraspecific genetic variation. However, studies investigating protist populations remain scarce, mostly due to limited 'omics data. Instead, most genetics studies of microeukaryotes have thus far relied on single loci, which can be misleading and do not easily allow for detection of recombination, a hallmark of sexual reproduction. Here, we analyze >40 genes from 72 single-cell transcriptomes from two morphospecies-Hyalosphenia papilio and Hyalosphenia elegans-of testate amoebae (Arcellinida, Amoebozoa) to assess genetic diversity in samples collected over four years from New England bogs. We confirm the existence of cryptic species based on our multilocus dataset, which provides evidence of recombination within and high levels of divergence between the cryptic species. At the same time, total levels of genetic diversity within cryptic species are low, suggesting that these abundant organisms have small effective population sizes, perhaps due to extinction and repopulation events coupled with efficient modes of dispersal. This study is one of the first to investigate population genetics in uncultivable heterotrophic protists using transcriptomics data and contributes towards understanding cryptic species of nonmodel microeukaryotes.

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.

Single-cell transcriptome reveals cell division-regulated hub genes in the unicellular eukaryote Paramecium

The transition from growth to division during the cell cycle encompasses numerous conserved processes such as large-scale DNA replication and protein synthesis. In ciliate cells, asexual cell division is accompanied by additional cellular changes including amitotic nuclear division, extensive ciliogenesis, and trichocyst replication. However, the molecular mechanisms underlying these processes remain elusive. In this study, we present single-cell gene expression profiles of Paramecium cf. multimicronucleatum cells undergoing cell division. Our results reveal that the most up-regulated genes in dividing cells compared to growing cells are associated with 1) cell cycle signaling pathways including transcription, DNA replication, chromosome segregation and protein degradation; 2) microtubule proteins and tubulin glycylases which are essential for ciliogenesis, nuclei separation and structural differentiation signaling; and 3) trichocyst matrix proteins involved in trichocyst synthesis and reproduction. Furthermore, weighted gene co-expression network analysis identified hub genes that may play crucial roles during cell division. Our findings provide insights into cell cycle regulators, microtubules and trichocyst matrix proteins that may exert influence on this process in ciliates.

Genomic insights into the phylogeny and biomass-degrading enzymes of rumen ciliates

Understanding the biodiversity and genetics of gut microbiomes has important implications for host physiology and industrial enzymes, whereas most studies have been focused on bacteria and archaea, and to a lesser extent on fungi and viruses. One group, still underexplored and elusive, is ciliated protozoa, despite its importance in shaping microbiota populations. Integrating single-cell sequencing and an assembly-and-identification pipeline, we acquired 52 high-quality ciliate genomes of 22 rumen morphospecies from 11 abundant morphogenera. With these genomes, we resolved the taxonomic and phylogenetic framework that revised the 22 morphospecies into 19 species spanning 13 genera and reassigned the genus Dasytricha from Isotrichidae to a new family Dasytrichidae. Comparative genomic analyses revealed that extensive horizontal gene transfers and gene family expansion provided rumen ciliate species with a broad array of carbohydrate-active enzymes (CAZymes) to degrade all major kinds of plant and microbial carbohydrates. In particular, the genomes of Diplodiniinae and Ophryoscolecinae species encode as many CAZymes as gut fungi, and ~80% of their degradative CAZymes act on plant cell-wall. The activities of horizontally transferred cellulase and xylanase of ciliates were experimentally verified and were 2-9 folds higher than those of the inferred corresponding bacterial donors. Additionally, the new ciliate dataset greatly facilitated rumen metagenomic analyses by allowing ~12% of the metagenomic sequencing reads to be classified as ciliate sequences.

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.

Transcriptomic Differences between Free-Living and Parasitic Chilodonella uncinata (Alveolata, Ciliophora)

Chilodonella uncinata is a facultatively parasitic ciliate, which can opportunistically parasitize on fish gills and fins, and sometimes even cause host mortality. Previous molecular studies of C. uncinata mainly focused on genetic diversity and molecular evolution. There are currently no transcriptome reports studying differences between free-living and parasitic C. uncinata. We addressed this by sequencing transcriptomes of these two C. uncinata lifestyle types using Smart-seq2 and Illumina HiSeq technologies. In total, 1040 differentially expressed genes (DEGs) were identified. Compared with the free-living type, 494 genes of the parasitic type were downregulated and 546 genes were upregulated. These DEGs were identified through BLAST with NCBI-nr, Swiss-Port, and Pfam databases and then annotated by GO enrichment and KEGG pathway analysis. The results showed that parasitism-related genes such as heat shock proteins (HSPs), actin I, and leishmanolysin were significantly upregulated in parasitic C. uncinata. The ciliary-related dynein heavy chain also had a higher expression in parasitic C. uncinata. Furthermore, there were significant differences in the amino acid metabolism, fatty acid metabolism, lipid metabolism, and TCA cycle. This study increases the volume of molecular data available for C. uncinata and contributes to our understanding of the mechanisms underlying the transition from a free-living to a parasitic lifestyle.

One Cell at a Time: Advances in Single-Cell Methods and Instrumentation for Discovery in Aquatic Microbiology

Studying microbes from a single-cell perspective has become a major theme and interest within the field of aquatic microbiology. One emerging trend is the unfailing observation of heterogeneity in activity levels within microbial populations. Wherever researchers have looked, intra-population variability in biochemical composition, growth rates, and responses to varying environmental conditions has been evident and probably reflect coexisting genetically distinct strains of the same species. Such observations of heterogeneity require a shift away from bulk analytical approaches and development of new methods or adaptation of existing techniques, many of which were first pioneered in other, unrelated fields, e.g., material, physical, and biomedical sciences. Many co-opted approaches were initially optimized using model organisms. In a field with so few cultivable models, method development has been challenging but has also contributed tremendous insights, breakthroughs, and stimulated curiosity. In this perspective, we present a subset of methods that have been effectively applied to study aquatic microbes at the single-cell level. Opportunities and challenges for innovation are also discussed. We suggest future directions for aquatic microbiological research that will benefit from open access to sophisticated instruments and highly interdisciplinary collaborations.

Molecular evolutionary analysis of the SM and SNARE vesicle fusion machinery in ciliates shows concurrent expansions in late secretory machinery

Protists in the phylum Ciliophora possess a complex membrane-trafficking system, including osmoregulatory Contractile Vacuoles and specialized secretory organelles. Molecular cell biological investigations in Tetrahymena thermophila have identified components of the protein machinery associated with the secretory organelles, mucocysts. The Qa-SNARE Syn7lp plays a role in mucocyst biogenesis as do subunits of the CORVET tethering complex (specifically Vps8). Indeed, Tetrahymena thermophila possesses expanded gene complements of several CORVET components, including Vps33 which is also a Sec1/Munc18 (SM) protein that binds Qa-SNAREs. Moreover, the Qa-SNAREs in Paramecium tetraurelia have been localized to various endomembrane organelles. Here, we use comparative genomics and phylogenetics to determine the evolutionary history of the SM and Qa-SNARE proteins across the Ciliophora. We identify that the last ciliate common ancestor possessed the four SM proteins and six Qa-SNAREs common to eukaryotes, including the uncommonly retained Syntaxin 17. We furthermore identify independent expansion of these protein families in several ciliate classes, including concurrent expansions of the SM protein-Qa SNARE partners Sec1:SynPM in the oligohymenophorean ciliates lineage, consistent with novel Contractile Vacuole specific innovations. Overall, these data are consistent with SM proteins and Qa-SNAREs being a common set of components for endomembrane modulation in the ciliates.

Paramecium epigenetics in development and proliferation

The term epigenetics is used for any layer of genetic information aside from the DNA base-sequence information. Mammalian epigenetic research increased our understanding of chromatin dynamics in terms of cytosine methylation and histone modification during differentiation, aging, and disease. Instead, ciliate epigenetics focused more on small RNA-mediated effects. On the one hand, these do concern the transport of RNA from parental to daughter nuclei, representing a regulated transfer of epigenetic information across generations. On the other hand, studies of Paramecium, Tetrahymena, Oxytricha, and Stylonychia revealed an almost unique function of transgenerational RNA. Rather than solely controlling chromatin dynamics, they control sexual progeny's DNA content quantitatively and qualitatively. Thus epigenetics seems to control genetics, at least genetics of the vegetative macronucleus. This combination offers ciliates, in particular, an epigenetically controlled genetic variability. This review summarizes the epigenetic mechanisms that contribute to macronuclear heterogeneity and relates these to nuclear dimorphism. This system's adaptive and evolutionary possibilities raise the critical question of whether such a system is limited to unicellular organisms or binuclear cells. We discuss here the relevance of ciliate genetics and epigenetics to multicellular organisms.

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.

The hidden genomic diversity of ciliated protists revealed by single-cell genome sequencing

Background

Ciliated protists are a widely distributed, morphologically diverse, and genetically heterogeneous group of unicellular organisms, usually known for containing two types of nuclei: a transcribed polyploid macronucleus involved in gene expression and a silent diploid micronucleus responsible for transmission of genetic material during sexual reproduction and generation of the macronucleus. Although studies in a few species of culturable ciliated protists have revealed the highly dynamic nature of replicative and recombination events relating the micronucleus to the macronucleus, the broader understanding of the genomic diversity of ciliated protists, as well as their phylogenetic relationships and metabolic potential, has been hampered by the inability to culture numerous other species under laboratory conditions, as well as the presence of symbiotic bacteria and microalgae which provide a challenge for current sequencing technologies. Here, we optimized single-cell sequencing methods and associated data analyses, to effectively remove contamination by commensal bacteria, and generated high-quality genomes for a number of Euplotia species.

Results

We obtained eight high-quality Euplotia genomes by using single-cell genome sequencing techniques. The genomes have high genomic completeness, with sizes between 68 and 125 M and gene numbers between 14K and 25K. Through comparative genomic analysis, we found that there are a large number of gene expansion events in Euplotia genomes, and these expansions are closely related to the phenotypic evolution and specific environmental adaptations of individual species. We further found four distinct subgroups in the genus Euplotes, which exhibited considerable genetic distance and relative lack of conserved genomic syntenies. Comparative genomic analyses of Uronychia and its relatives revealed significant gene expansion associated with the ciliary movement machinery, which may be related to the unique and strong swimming ability.

Conclusions

We employed single-cell genomics to obtain eight ciliate genomes, characterized the underestimated genomic diversity of Euplotia, and determined the divergence time of representative species in this subclass for the first time. We also further investigated the extensive duplication events associated with speciation and environmental adaptation. This study provides a unique and valuable resource for understanding the evolutionary history and genetic diversity of ciliates.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01202-1.

Evolutionary insights and brief review of Loxodes Ehrenberg, 1830 (Ciliophora, Karyorelictea, Loxodidae) with description of a new species from Mexico

Karyorelictids are a group of ciliates inhabiting marine and freshwater biotopes and possessing a non-dividing macronucleus. We describe a new freshwater species based on morphology and the 18S rRNA gene sequence data. Loxodes tziscaensis n. sp. can be easily distinguished from other Loxodes species by the arrangement of the nuclear apparatus and features of the buccal and somatic ciliature. The current proposed 18S rRNA phylogeny of Loxodes, including seven Loxodes species, shows two morphologically well-supported groups. Group A (L. rostrum, type species; L. vorax and L. tziscaensis n. sp.) includes species with a single nuclear group (two macronuclei and one micronucleus), in contrast to species of group B, which possess more than one nuclear group (L. striatus, L. magnus, L. kahli, L. penardi, and L. rex). We propose that the last common ancestor of Loxodes was a marine Remanella-like species possessing a single nuclear group. The division and differentiation of the micronucleus into a new macronucleus and the retention of the old macronuclei, independently of cell division, may have been two crucial processes during the evolution and diversification of Loxodes species with one nuclear group into species with multiple nuclear groups.

Chromosome organization and gene expansion in the highly fragmented genome of the ciliate Strombidium stylifer

Chromosomes are well-organized carriers of genetic information in eukaryotes and are usually quite long, carrying hundreds and thousands of genes. Intriguingly, a clade of single-celled ciliates, Spirotrichea, feature nanochromosomes-also called "gene-sized chromosomes". These chromosomes predominantly carry only one gene, flanked by short telomere sequences. However, the organization and copy number variation of the chromosomes in these highly fragmented genomes remain unexplored in many groups of Spirotrichea, including the marine Strombidium. Using deep genome sequencing, we assembled the macronuclear genome of Strombidium stylifer into more than 18,000 nanochromosomes (~2.4 Kb long on average). Our results show that S. stylifer occupies an intermediate position during the evolutionary history of Strombidium lineage and experienced significant expansions in several gene families related to guanyl ribonucleotide binding. Based on the nucleotide distribution bias analysis and conserved motifs search in non-genic regions, we found that the subtelomeric regions have a conserved adenine-thymine (AT)-rich sequence motif. We also found that the copy number of nanochromosomes lacks precise regulation. This work sheds light on the unique features of chromosome structure in eukaryotes with highly fragmented genomes and reveals that a rather specialized evolutionary strategy at the genomic level has resulted in great diversity within the ciliated lineages.

DAPI staining and DNA content estimation of nuclei in uncultivable microbial eukaryotes (Arcellinida and Ciliates)

Though representing a major component of eukaryotic biodiversity, many microbial eukaryotes remain poorly studied, including the focus of the present work, testate amoebae of the order Arcellinida (Amoebozoa) and non-model lineages of ciliates (Alveolata). In particular, knowledge of genome structures and changes in genome content over the often-complex life cycles of these lineages remains enigmatic. However, the limited available knowledge suggests that microbial eukaryotes have the potential to challenge our textbook views on eukaryotic genomes and genome evolution. In this study, we developed protocols for DAPI (4',6-diamidino-2-phenylindole) staining of Arcellinida nuclei and adapted protocols for ciliates. In addition, image analysis software was used to estimate the DNA content in the nuclei of Arcellinida and ciliates, and the measurements of target organisms were compared to those  of well-known model organisms.The results demonstrate that the methods we have developed for nuclear staining in these lineages are effective and can be applied to other microbial eukaryotic groups by adjusting certain stages in the protocols.

Biodiversity-based development and evolution: the emerging research systems in model and non-model organisms

Evolutionary developmental biology, or Evo-Devo for short, has become an established field that, broadly speaking, seeks to understand how changes in development drive major transitions and innovation in organismal evolution. It does so via integrating the principles and methods of many subdisciplines of biology. Although we have gained unprecedented knowledge from the studies on model organisms in the past decades, many fundamental and crucially essential processes remain a mystery. Considering the tremendous biodiversity of our planet, the current model organisms seem insufficient for us to understand the evolutionary and physiological processes of life and its adaptation to exterior environments. The currently increasing genomic data and the recently available gene-editing tools make it possible to extend our studies to non-model organisms. In this review, we review the recent work on the regulatory signaling of developmental and regeneration processes, environmental adaptation, and evolutionary mechanisms using both the existing model animals such as zebrafish and Drosophila, and the emerging nonstandard model organisms including amphioxus, ascidian, ciliates, single-celled phytoplankton, and marine nematode. In addition, the challenging questions and new directions in these systems are outlined as well.

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.

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.

The Compact Macronuclear Genome of the Ciliate Halteria grandinella: A Transcriptome-Like Genome with 23,000 Nanochromosomes

How to achieve protein diversity by genome and transcriptome processing is essential for organismal complexity and adaptation. The present work identifies that the macronuclear genome of Halteria grandinella, a cosmopolitan unicellular eukaryote, is composed almost entirely of gene-sized nanochromosomes with extremely short nongenic regions.

How to achieve protein diversity by genome and transcriptome processing is essential for organismal complexity and adaptation. The present work identifies that the macronuclear genome of Halteria grandinella, a cosmopolitan unicellular eukaryote, is composed almost entirely of gene-sized nanochromosomes with extremely short nongenic regions. This challenges our usual understanding of chromosomal structure and suggests the possibility of novel mechvanisms in transcriptional regulation. Comprehensive analysis of multiple data sets reveals that Halteria transcription dynamics are influenced by: (i) nonuniform nanochromosome copy numbers correlated with gene-expression level; (ii) dynamic alterations at both the DNA and RNA levels, including alternative internal eliminated sequence (IES) deletions during macronucleus formation and large-scale alternative splicing in transcript maturation; and (iii) extremely short 5′ and 3′ untranslated regions (UTRs) and universal TATA box-like motifs in the compact 5′ subtelomeric regions of most chromosomes. This study broadens the view of ciliate biology and the evolution of unicellular eukaryotes, and identifies Halteria as one of the most compact known eukaryotic genomes, indicating that complex cell structure does not require complex gene architecture.

New contributions to the phylogeny of the ciliate class Heterotrichea (Protista, Ciliophora): analyses at family-genus level and new evolutionary hypotheses

Heterotrichous ciliates play an important role in aquatic ecosystem energy flow processes and many are model organisms for research in cytology, regenerative biology, and toxicology. In the present study, we combine both morphological and molecular data to infer phylogenetic relationships at family-genus level and propose new evolutionary hypotheses for the class Heterotrichea. The main results include: (1) 96 new ribosomal DNA sequences from 36 populations, representing eight families and 13 genera, including three poorly annotated genera, Folliculinopsis, Ampullofolliculina and Linostomella; (2) the earliest-branching families are Spirostomidae in single-gene trees and Peritromidae in the concatenated tree, but the family Peritromidae probably represents the basal lineage based on its possession of many "primitive" morphological characters; (3) some findings in molecular trees are not supported by morphological evidence, such as the family Blepharismidae is one of the most recent branches and the relationship between Fabreidae and Folliculinidae is very close; (4) the systematic positions of Condylostomatidae, Climacostomidae, and Gruberiidae remain uncertain based either on morphological or molecular data; and (5) the monophyly of each genus included in the present study is supported by the molecular phylogenetic trees, except for Blepharisma in the SSU rDNA tree and Folliculina in the ITS1-5.8S-ITS2 tree.

Combined Genome and Transcriptome Analyses of the Ciliate Schmidingerella arcuata (Spirotrichea) Reveal Patterns of DNA Elimination, Scrambling, and Inversion

Schmidingerella arcuata is an ecologically important tintinnid ciliate that has long served as a model species in plankton trophic ecology. We present a partial micronuclear genome and macronuclear transcriptome resource for S. arcuata, acquired using single-cell techniques, and we report on pilot analyses including functional annotation and genome architecture. Our analysis shows major fragmentation, elimination, and scrambling in the micronuclear genome of S. arcuata. This work introduces a new nonmodel genome resource for the study of ciliate ecology and genomic biology and provides a detailed functional counterpart to ecological research on S. arcuata.

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.

The completed macronuclear genome of a model ciliate Tetrahymena thermophila and its application in genome scrambling and copy number analyses

The ciliate Tetrahymena thermophila has been a powerful model system for molecular and cellular biology. However, some investigations have been limited due to the incomplete closure and sequencing of the macronuclear genome assembly, which for many years has been stalled at 1,158 scaffolds, with large sections of unknown sequences (available in Tetrahymena Genome Database, TGD, http://ciliate.org/ ). Here we completed the first chromosome-level Tetrahymena macronuclear genome assembly, with approximately 300× long Single Molecule, Real-Time reads of the wild-type SB210 cells-the reference strain for the initial macronuclear genome sequencing project. All 181 chromosomes were capped with two telomeres and gaps were entirely closed. The completed genome shows significant improvements over the current assembly (TGD 2014) in both chromosome structure and sequence integrity. The majority of previously identified gene models shown in TGD were retained, with the addition of 36 new genes and 883 genes with modified gene models. The new genome and annotation were incorporated into TGD. This new genome allows for pursuit in some underexplored areas that were far more challenging previously; two of them, genome scrambling and chromosomal copy number, were investigated in this study. We expect that the completed macronuclear genome will facilitate many studies in Tetrahymena biology, as well as multiple lines of research in other eukaryotes.