A broad molecular phylogeny of ciliates: identification of major evolutionary trends and radiations within the phylum

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.

Calmodulin in Paramecium: Focus on Genomic Data

Calcium (Ca²⁺) is a universal second messenger that plays a key role in cellular signaling. However, Ca²⁺ signals are transduced with the help of Ca²⁺-binding proteins, which serve as sensors, transducers, and elicitors. Among the collection of these Ca²⁺-binding proteins, calmodulin (CaM) emerged as the prototypical model in eukaryotic cells. This is a small protein that binds four Ca²⁺ ions and whose functions are multiple, controlling many essential aspects of cell physiology. CaM is universally distributed in eukaryotes, from multicellular organisms, such as human and land plants, to unicellular microorganisms, such as yeasts and ciliates. Here, we review most of the information gathered on CaM in Paramecium, a group of ciliates. We condense the information here by mentioning that mature Paramecium CaM is a 148 amino acid-long protein codified by a single gene, as in other eukaryotic microorganisms. In these ciliates, the protein is notoriously localized and regulates cilia function and can stimulate the activity of some enzymes. When Paramecium CaM is mutated, cells show flawed locomotion and/or exocytosis. We further widen this and additional information in the text, focusing on genomic data.

Integrative studies on the taxonomy and molecular phylogeny of four new Pleuronema species (Protozoa, Ciliophora, Scuticociliatia)

The ciliate genus Pleuronema comprises approximately 30 nominal species and has been reported in freshwater, brackish water, and marine habitats. Nevertheless, recent studies have indicated that there might be a large undiscovered species diversity. In the present work, four new Pleuronema species, namely P. foissneri sp. nov., P. parasmalli sp. nov., P. parasalmastra sp. nov., and P. paraorientale sp. nov., collected from Shenzhen, southern China, was investigated using taxonomic methods. The diagnosis, description, comparisons with morphologically related species and detailed morphometric data are supplied for each. The small subunit ribosomal RNA (SSU rRNA) gene of the four new species is sequenced and their molecular phylogeny is analyzed. The SSU rRNA gene tree shows that Pleuronema is polyphyletic comprising several separate clades. All four new species cluster consistently with P. orientale KF206429, P. puytoraci KF840520 and P. setigerum FJ848874 within the core Pleuronematidae + Peniculistomatidae clade. Phylogenies of Pleuronematidae-related taxa are also discussed.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-022-00130-5.

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.

Diluted Seawater and Ammonia-N Tolerance of Two Mangrove Crab Species. New Insights to Understand the Vulnerability of Pristine Islands Ecosystems Organisms

Mangrove ecosystems are the primary receptors of anthropogenic pollution in tropical areas. Assessing the vulnerability of these ecosystems can be expressed, among other indicators, by studying the health of ‘ecosystem engineers’. In this study, mangrove forests facing opposing anthropogenic pressures were studied (i) in the uninhabited island of Europa (Mozambique Channel), considered as a pristine ecosystem, and, (ii) on the island of Mayotte, facing regular domestic wastewater discharges. Using an ecophysiological approach, the effects of diluted seawater (DSW) and increased ammonia-N were studied for two fiddler crab species: Gelasimus tetragonon (GT) on the island of Europa and Paraleptuca chlorophthalmus (PC) on the island of Mayotte. Osmoregulation curves and osmoregulatory capacity were determined along with O2 consumption rates after a 96 h exposure period. Histological analyses were also carried out on two important metabolic organs: the hepatopancreas and the posterior gills. Results indicate that both crab species are good hyper-hypo-osmoregulators but only PC can maintain its osmoregulatory capacity when exposed to ammonia-N. Oxygen consumption is increased in GT after 96 h of exposure to ammonia-N but this does not occur in PC. Finally, a thickening of the gill osmoregulatory epithelium was observed after 96 h in PC when exposed to ammonium but not in GT. Therefore, the two species do not have the same tolerance to DSW and increased ammonia-N. PC shows physiological acclimation capacities in order to better manage nitrogenous enrichments. GT did not show the same physiological plasticity when exposed to ammonia-N and could be more at risk by this kind of stress. These results along with those from other studies regarding the effects of domestic effluents on mangrove crabs are discussed. Therefore, the greater vulnerability of organisms occupying pristine ecosystems could induce major changes in mangrove functioning if crabs, that are engineer species of the ecosystem, are about to reduce their bioturbation activity or, even, disappear from the mangrove forests.

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.

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

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

Discrepancies Between Molecular and Morphological Databases of Soil Ciliates Studied for Temperate Grasslands of Central Europe

By measuring the change in soil protist communities, the effect of human land use on grasslands can be monitored to promote sustainable ecosystem functioning. Protists form the active link in the rhizosphere between the plant roots and higher trophic organisms; however, only few morphological species and their ecological values have yet been described in this context. To investigate the communicability between morphological and molecular databases used in the molecular barcoding of protists and in the biomonitoring of grassland soil, the present high-throughput sequencing (HTS) study (N=150) covered the area of central Europe (mesoscale) known to be well studied for ciliated protists. HTS delivered 2,404 unique reads identifying taxa in all major ciliophoran classes but exact reference matches were few. The study identified clear discrepancies between databases for well-studied taxa, where molecular databases contained multiple gene variants for single morphospecies of dominant taxa. Gene variants presented own biogeography - the eukaryotic microdiversity along gradients (e.g., land-use intensity, soil water). It is possible that many of the so called novel phylogenetic lineages and hidden diversity pointed out in environmental surveys could be evidence for the severe lack of molecular data for already known and morphologically described species, present in morphological databases.

Genomics: Stentor's Trumpet Sounds Anew

The amazing regenerative abilities of the giant ciliate Stentor coeruleus made it a favorite subject for classical embryologists. Now, its genome has been sequenced, enabling renewed experimental study and revealing unexpected surprises in mRNA splicing and the genetic code.

Structure of the germline genome of Tetrahymena thermophila and relationship to the massively rearranged somatic genome

The germline genome of the binucleated ciliate Tetrahymena thermophila undergoes programmed chromosome breakage and massive DNA elimination to generate the somatic genome. Here, we present a complete sequence assembly of the germline genome and analyze multiple features of its structure and its relationship to the somatic genome, shedding light on the mechanisms of genome rearrangement as well as the evolutionary history of this remarkable germline/soma differentiation. Our results strengthen the notion that a complex, dynamic, and ongoing interplay between mobile DNA elements and the host genome have shaped Tetrahymena chromosome structure, locally and globally. Non-standard outcomes of rearrangement events, including the generation of short-lived somatic chromosomes and excision of DNA interrupting protein-coding regions, may represent novel forms of developmental gene regulation. We also compare Tetrahymena's germline/soma differentiation to that of other characterized ciliates, illustrating the wide diversity of adaptations that have occurred within this phylum.

Methods for taxonomic studies of rumen ciliates (alveolata: ciliophora): a brief review

This review presents the principal methods used in taxonomic studies of rumen ciliates: live observation, Lugol staining, fixation and staining with methyl-green formalin saline (MFS) solution, protargol staining, silver carbonate impregnation, scanning electron microscopy and molecular techniques. Mastering these techniques is essential for successful research on the taxonomy of rumen ciliates. No single technique reveals all of the characteristics required for a complete description of a rumen ciliate; therefore, it is necessary to combine the use of these techniques as appropriate to the rumen ciliate group under study. Tables are provided to summarize: 1) morphological methods more appropriate for revealing morphological structures of interest, 2) morphological methods indicated for each group of rumen ciliates, and 3) main primers used for PCR amplification of the 18S rDNA of rumen ciliates.

Fate of engineered cerium oxide nanoparticles in an aquatic environment and their toxicity toward 14 ciliated protist species

The potential environmental impacts of engineered cerium oxide nanoparticles (CeO2 NPs) on aquatic organisms have remained largely unknown. Therefore, the laboratory study featured herein was performed to determine the fate of CeO2 NPs in an aquatic environment and their toxicity towards 14 different ciliated protist species at a specified population level. An investigation of 48 h aggregation kinetics in the Dryl's solution showed the CeO2 NPs to be relatively stable. The pH values in three test medium were too far away from PZC, which explained the stability of CeO2 NPs. CeO2 NPs generally elicited more toxicity with increasing NP concentration, following certain dose-response relationships. Nano-CeO2 resulted in greater toxicity in a particle state than when added as bulk material. LC50 values showed a negative correlation with the surface-to-volume ratio for these protists, suggesting that surface adsorption of CeO2 NPs might contribute to the observed toxicity. Additionally, acute toxic responses of 14 ciliated protist species to CeO2 NPs were not significantly phylogenetically conserved. The results of these observations provide a better insight into the potential risks of CeO2 NPs in an aquatic environment.

The All-Data-Based Evolutionary Hypothesis of Ciliated Protists with a Revised Classification of the Phylum Ciliophora (Eukaryota, Alveolata)

The phylum Ciliophora plays important roles in a wide range of biological studies. However, the evolutionary relationships of many groups remain unclear due to a lack of sufficient molecular data. In this study, molecular dataset was expanded with representatives from 55 orders and all major lineages. The main findings are: (1) 14 classes were recovered including one new class, Protocruziea n. cl.; (2) in addition to the two main branches, Postciliodesmatophora and Intramacronucleata, a third branch, the Mesodiniea, is identified as being basal to the other two subphyla; (3) the newly defined order Discocephalida is revealed to be a sister clade to the euplotids, strongly suggesting the separation of discocephalids from the hypotrichs; (4) the separation of mobilids from the peritrichs is not supported; (5) Loxocephalida is basal to the main scuticociliate assemblage, whereas the thigmotrichs are placed within the order Pleuronematida; (6) the monophyly of classes Phyllopharyngea, Karyorelictea, Armophorea, Prostomatea, Plagiopylea, Colpodea and Heterotrichea are confirmed; (7) ambiguous genera Askenasia, CyclotrichiumParaspathidium and Plagiocampa show close affiliation to the well known plagiopyleans; (8) validity of the subclass Rhynchostomatia is supported, and (9) the systematic positions of Halteriida and Linconophoria remain unresolved and are thus regarded as incertae sedis within Spirotrichea.

Expanded phylogenetic analyses of the class Heterotrichea (Ciliophora, Postciliodesmatophora) using five molecular markers and morphological data

Most studies of the molecular evolution of Heterotrichea have been based solely on the 18S-rDNA gene, which were inconsistent with morphological classification. Because of the limitations of single locus phylogenies and the recurring problem of lack of resolution of deeper nodes found in previous studies, we present hypotheses of the evolution of internal groups of the class Heterotrichea based on multi-loci analyses (18S-rDNA, 28S-rDNA, ITS1-5.8S-ITS2 region, COI and alpha-tubulin) and morphological data. Phylogenetic trees from protein coding gene data are presented for Heterotrichea for the first time. Phylogenetic analyses included Bayesian inference, maximum likelihood, maximum parsimony methods, and optimal trees were statistically compared to alternative topologies from the literature. Additionally, the Bayesian concordance approach (BCA algorithm) was used to assess the concordance factor between topologies obtained from isolated analyses. Because different loci may evolve at different rates, resulting in different gene topologies, we also estimated a species tree for Heterotrichea using the STAR coalescence-based method. The results show that: (1) single gene trees are inconsistent regarding the position of some heterotrichean families; (2) the concatenation of all data in a total-evidence tree improved the resolution of deep nodes among the heterotrichean families and genera; (3) the coalescent-based species tree is consistent with phylogenies based on the 18S-rDNA gene and shows Spirostomidae as the stem group of Heterotrichea; (4) however, the total-evidence tree suggests that the large Heterotrichea cluster is divided into nine lineages in which Peritromidae diverges at the base of the Heterotrichea tree.

Cooperativity and evolution of Tetrahymena two-domain arginine kinase

Tetrahymena pyriformis contains two arginine kinases, a 40-kDa enzyme (AK1) with a myristoylation signal sequence at the N-terminus and a two-domain 80-kDa enzyme (AK2). The former is localized mainly in cilia and the latter is in the cytoplasm. AK1 was successfully synthesized using an insect cell-free protein synthesis system and subjected to peptide mass fingerprinting (PMF) analysis. The masses corresponding to unmodified N-terminal tryptic peptide or N-terminal myristoylated peptide were not observed, suggesting that N-terminal peptides were not ionized in this analysis. We performed PMF analyses for two other phosphagen kinases (PKs) with myristoylation signals, an AK from Nematostella vectensis and a PK from Ectocarpus siliculosus. In both cases, the myristoylated, N-terminal peptides were clearly identified. The differences between the experimental and theoretical masses were within 0.0165-0.0583 Da, supporting the accuracy of the identification. Domains 1 and 2 of Tetrahymena two-domain AK2 were expressed separately in Escherichia coli and the extent of cooperativity was estimated on the basis of their kinetic constants. The results suggested that each of the domains functions independently, namely no cooperativity is displayed between the two domains. This is in sharp contrast to the two-domain AK from Anthopleura.

Transposon domestication versus mutualism in ciliate genome rearrangements

Ciliated protists rearrange their genomes dramatically during nuclear development via chromosome fragmentation and DNA deletion to produce a trimmer and highly reorganized somatic genome. The deleted portion of the genome includes potentially active transposons or transposon-like sequences that reside in the germline. Three independent studies recently showed that transposase proteins of the DDE/DDD superfamily are indispensible for DNA processing in three distantly related ciliates. In the spirotrich Oxytricha trifallax, high copy-number germline-limited transposons mediate their own excision from the somatic genome but also contribute to programmed genome rearrangement through a remarkable transposon mutualism with the host. By contrast, the genomes of two oligohymenophorean ciliates, Tetrahymena thermophila and Paramecium tetraurelia, encode homologous PiggyBac-like transposases as single-copy genes in both their germline and somatic genomes. These domesticated transposases are essential for deletion of thousands of different internal sequences in these species. This review contrasts the events underlying somatic genome reduction in three different ciliates and considers their evolutionary origins and the relationships among their distinct mechanisms for genome remodeling.

Insights into the phylogeny of systematically controversial haptorian ciliates (Ciliophora, Litostomatea) based on multigene analyses

The ciliate subclass Haptoria is a diverse taxon that includes most of the free-living predators in the class Litostomatea. Phylogenetic study of this group was initially conducted using a single molecular marker small-subunit ribosomal RNA (SSU rRNA genes). Multi-gene analysis has been limited because very few other sequences were available. We performed phylogenetic analyses of Haptoria incorporating new SSU rRNA gene sequences from several debated members of the taxon, in particular, the first molecular data from Cyclotrichium. We also provided nine large-subunit ribosomal RNA (LSU rRNA) gene sequences and 10 alpha-tubulin sequences from diverse haptorians, and two possible relatives of controversial haptorians (Plagiopylea, Prostomatea). Phylogenies inferred from the different molecules showed the following: (i) Cyclotrichium and Paraspathidium were clearly separated from the haptorids and even from class Litostomatea, rejecting their high-level taxonomic assignments based on morphology. Both genera branch instead with the classes Plagiopylea, Prostomatea and Oligohymenophora. This raises the possibility that the well-known but phylogenetically problematic cyclotrichiids Mesodinium and Myrionecta may also have affinities here, rather than with litostomes; (ii) the transfer of Trachelotractus to Litostomatea is supported, especially by the analyses of SSU rRNA and LSU rRNA genes, however, Trachelotractus and Chaenea (more uncertainly) generally form the two deepest lineages within litostomes; and (iii) phylogenies of the new molecular markers are consistent with SSU rRNA gene information in recovering order Pleurostomatida as monophyletic. However, Pleurostomatida branches cladistically within order Haptorida, as does subclass Trichostomatia (on the basis of SSU rRNA phylogenies). Our results suggest that the class-level taxonomy of ciliates is still not resolved, and also that a systematic revision of litostomes is required, beginning at high taxonomic levels (taxa currently ranked as subclasses and orders).

Congruence and indifference between two molecular markers for understanding oral evolution in the Marynidae sensu lato (Ciliophora, Colpodea)

Our understanding of the evolution of oral structures within the Colpodida is confounded by the low number of morphological characters that can be used in constructing hypotheses, and by the low taxon and character sampling in molecular phylogenetic analyses designed to assess these hypotheses. Here we increase character sampling by sequencing the mitochondrial SSU-rDNA locus for three isolates of the Marynidae sensu lato. We show that the inferred mitochondrial and nuclear SSU-rDNA trees, as well as concatenated and constrained analyses, are congruent in not recovering a monophyletic Marynidae. However, due to low node support, the trees are indifferent to whether the morphological characters used to unite the Marynidae are the result of retention of ancestral states or convergence. In light of this indifference and an increased amount of nuclear and mitochondrial SSU-rDNA data, alternative hypotheses of oral evolution in the Colpodida are presented.

DNA elimination in ciliates: transposon domestication and genome surveillance

Ciliated protozoa extensively remodel their somatic genomes during nuclear development, fragmenting their chromosomes and removing large numbers of internal eliminated sequences (IESs). The sequences eliminated are unique and repetitive DNAs, including transposons. Recent studies have identified transposase proteins that appear to have been domesticated and are used by these cells to eliminate DNA not wanted in the somatic macronucleus. This DNA elimination process is guided by meiotically produced small RNAs, generated in the germline nucleus, that recognize homologous sequences leading to their removal. These scan RNAs are found in complexes with PIWI proteins. Before they search the developing genome for IESs to eliminate, they scan the parental somatic nucleus and are removed from the pool if they match homologous sequences in that previously reorganized genome. In Tetrahymena, the scan RNAs target heterochromatin modifications to mark IESs for elimination. This DNA elimination pathway in ciliates shares extensive similarity with piRNA-mediated silencing of metazoans and highlights the remarkable ability of homologous RNAs to shape developing genomes.

Small subunit rRNA phylogenies show that the class nassophorea is not monophyletic (Phylum Ciliophora)

The hypostome ciliates have been generally classified into two classes, Phyllopharyngea and Nassophorea. The status of Nassophorea and its relationship with Phyllopharyngea is one of the most controversial issues in ciliate systematics. Here we focus on the phylogenetic interrelationships of Nassophorea and Phyllopharyngea based on small subunit ribosomal RNA gene sequences. The three nassophorean subgroups, synhymeniids, microthoracids, and nassulids, each emerged as monophyletic, with synhymeniids as a sister group of Phyllopharyngea, and microthoracids as a sister of the synhymeniids+Phyllopharyngea clade in all phylogenies. The exact placement of the nassulids, however, remains uncertain. Following a detailed analysis of phenotypic characters, we hypothesize that: (1) the Phyllopharyngea could have evolved from synhymeniids, with the further development of their subkinetal microtubules as one of the major events; and (2) the development of monokinetid structures, as well as the reduction and specialization of the cyrtos and cortex, might have occurred during the diversifications of the microthoracids, synhymeniids, and Phyllopharyngea from a common ancestor. Expanding the class Phyllopharyngea to include the synhymeniids as a subclass, and designating a new subclass Subkinetalia n. subcl. for the group comprising cyrtophorians, chonotrichians, rhynchodians, and suctorians, are proposed.

Further consideration of the phylogeny of some "traditional" heterotrichs (Protista, Ciliophora) of uncertain affinities, based on new sequences of the small subunit rRNA gene

The systematic relationships and taxonomic positions of the traditional heterotrich genera Condylostentor, Climacostomum, Fabrea, Folliculina, Peritromus, and Condylostoma, as well as the licnophorid genus Licnophora, were re-examined using new data from sequences of the gene coding for small subunit ribosomal RNA. Trees constructed using distance-matrix, Bayesian inference, and maximum-parsimony methods all showed the following relationships: (1) the "traditional" heterotrichs consist of several paraphyletic groups, including the current classes Heterotrichea, Armophorea and part of the Spirotrichea; (2) the class Heterotrichea was confirmed as a monophyletic assemblage based on our analyses of 31 taxa, and the genus Peritromus was demonstrated to be a peripheral group; (3) the genus Licnophora occupied an isolated branch on one side of the deepest divergence in the subphylum Intramacronucleata and was closely affiliated with spirotrichs, armophoreans, and clevelandellids; (4) Condylostentor, a recently defined genus with several truly unique morphological features, is more closely related to Condylostoma than to Stentor; (5) Folliculina, Eufolliculina, and Maristentor always clustered together with high bootstrap support; and (6) Climacostomum occupied a paraphyletic position distant from Fabrea, showing a close relationship with Condylostomatidae and Chattonidiidae despite of modest support.

Refined annotation and assembly of the Tetrahymena thermophila genome sequence through EST analysis, comparative genomic hybridization, and targeted gap closure

BACKGROUND

Tetrahymena thermophila, a widely studied model for cellular and molecular biology, is a binucleated single-celled organism with a germline micronucleus (MIC) and somatic macronucleus (MAC). The recent draft MAC genome assembly revealed low sequence repetitiveness, a result of the epigenetic removal of invasive DNA elements found only in the MIC genome. Such low repetitiveness makes complete closure of the MAC genome a feasible goal, which to achieve would require standard closure methods as well as removal of minor MIC contamination of the MAC genome assembly. Highly accurate preliminary annotation of Tetrahymena's coding potential was hindered by the lack of both comparative genomic sequence information from close relatives and significant amounts of cDNA evidence, thus limiting the value of the genomic information and also leaving unanswered certain questions, such as the frequency of alternative splicing.

RESULTS

We addressed the problem of MIC contamination using comparative genomic hybridization with purified MIC and MAC DNA probes against a whole genome oligonucleotide microarray, allowing the identification of 763 genome scaffolds likely to contain MIC-limited DNA sequences. We also employed standard genome closure methods to essentially finish over 60% of the MAC genome. For the improvement of annotation, we have sequenced and analyzed over 60,000 verified EST reads from a variety of cellular growth and development conditions. Using this EST evidence, a combination of automated and manual reannotation efforts led to updates that affect 16% of the current protein-coding gene models. By comparing EST abundance, many genes showing apparent differential expression between these conditions were identified. Rare instances of alternative splicing and uses of the non-standard amino acid selenocysteine were also identified.

CONCLUSION

We report here significant progress in genome closure and reannotation of Tetrahymena thermophila. Our experience to date suggests that complete closure of the MAC genome is attainable. Using the new EST evidence, automated and manual curation has resulted in substantial improvements to the over 24,000 gene models, which will be valuable to researchers studying this model organism as well as for comparative genomics purposes.

Molecular phylogeny of the Heterotrichea (Ciliophora, Postciliodesmatophora) based on small subunit rRNA gene sequences

A comprehensive molecular analysis of the phylogenetic relationships within the Heterotrichea including all described families is still lacking. For this reason, the complete nuclear small subunit (SSU) rDNA was sequenced from further representatives of the Blepharismidae and the Stentoridae. In addition, the SSU rDNA of a new, undescribed species of the genus Condylostomides (Condylostomatidae) was sequenced. The detailed phylogenetic analyses revealed a consistent branching pattern: while the terminal branches are generally well resolved, the basal relationships remain unsolved. Moreover, the data allow some conclusions about the macronuclear evolution within the genera Blepharisma, Stentor, and Spirostomum suggesting that a single, compact macronucleus represents the ancestral state.

Molecular phylogeny of Stentor (Ciliophora: Heterotrichea) based on small subunit ribosomal RNA sequences

To determine the phylogenetic position of Stentor within the Class Heterotrichea, the complete small subunit rRNA genes of three Stentor species, namely Stentor polymorphus, Stentor coeruleus, and Stentor roeseli, were sequenced and used to construct phylogenetic trees using the maximum parsimony, neighbor joining, and Bayesian analysis. With all phylogenetic methods, the genus Stentor was monophyletic, with S. roeseli branching basally.

Molecular systematics and ultrastructural characterization of a forgotten species: Chattonidium setense (Ciliophora, Heterotrichea)

In the present paper we redescribe the ciliate Chattonidium setense Villeneuve 1937 combining morphological observations (live, stained, scanning, and transmission electron microscope) with behavioral notes and molecular data. Ultrastructural analysis revealed remarkable similarities between Chattonidium and representative members of the class Heterotrichea in cortical structure and cytoplasmic organization. The most similar genus for these aspects appears to be Condylostoma. To verify this relatedness, 18S rRNA genes from Chattonidium and from one Condylostoma species were sequenced. Phylogenetic analysis indicates Chattonidium belongs to the class Heterotrichea defined according to the modern taxonomy, and confirms its relatedness with Condylostoma already hypothesized by Villeneuve-Brachon (1940). The presence of the aboral cavity complex, a unique feature never described in other ciliates, and its peculiar organization revealed by ultrastructural analysis fully justify, in our opinion, the maintenance of Chattonidium in the separate family Chattonidiidae, established by Villeneuve-Brachon in 1940.

Mutations of tubulin glycylation sites reveal cross-talk between the C termini of alpha- and beta-tubulin and affect the ciliary matrix in Tetrahymena

Two types of polymeric post-translational modifications of alpha/beta-tubulin, glycylation and glutamylation, occur widely in cilia and flagella. Their respective cellular functions are poorly understood. Mass spectrometry and immunoblotting showed that two closely related species, the ciliates Tetrahymena and Paramecium, have dramatically different compositions of tubulin post-translational modifications in structurally identical axonemes. Whereas the axonemal tubulin of Paramecium is highly glycylated and has a very low glutamylation content, the axonemal tubulin of Tetrahymena is glycylated and extensively glutamylated. In addition, only the alpha-tubulin of Tetrahymena undergoes detyrosination. Mutations of the known glycylation sites in Tetrahymena tubulin affected the level of each polymeric modification type in both the mutated and nonmutated subunits, revealing cross-talk between alpha- and beta-tubulin. Ultrastructural analyses of glycylation site mutants uncovered defects in the doublet B-subfiber of axonemes and revealed an accumulation of dense material in the ciliary matrix, reminiscent of intraflagellar transport particles seen by others in Chlamydomonas. We propose that polyglycylation and/or polyglutamylation stabilize the B-subfiber of outer doublets and regulate the intraflagellar transport.

A Cladistic Approach for the Classification of Oligotrichid Ciliates (Ciliophora: Spirotricha)

Currently, gene sequence genealogies of the Oligotrichea Bütschli, 1889 comprise only few species. Therefore, a cladistic approach, especially to the Oligotrichida, was made, applying Hennig's method and computer programs. Twenty-three characters were selected and discussed, i.e., the morphology of the oral apparatus (five characters), the somatic ciliature (eight characters), special organelles (four characters), and ontogenetic particulars (six characters). Nine of these characters developed convergently twice. Although several new features were included into the analyses, the cladograms match other morphological trees in the monophyly of the Oligotrichea, Halteriia, Oligotrichia, Oligotrichida, and Choreotrichida. The main synapomorphies of the Oligotrichea are the enantiotropic division mode and the de novo-origin of the undulating membranes. Although the sister group relationship of the Halteriia and the Oligotrichia contradicts results obtained by gene sequence analyses, no morphologic, ontogenetic or ultrastructural features were found, which support a branching of Halteria grandinella within the Stichotrichida. The cladistic approaches suggest paraphyly of the family Strombidiidae probably due to the scarce knowledge. A revised classification of the Oligotrichea is suggested, including all sufficiently known families and genera.

Phylogenetic position of Licnophora, Lechriopyla, and Schizocaryum, three unusual ciliates (phylum Ciliophora) endosymbiotic in echinoderms (phylum Echinodermata)

Various echinoderms are colonized by species from several classes of the Phylum Ciliophora, indicating that the echinoderm "habitat" has been invaded independently on numerous occasions throughout evolutionary history. Two "echinoderm" ciliates whose phylogenetic positions have been problematic are Licnophora macfarlandi Stevens, 1901 and Schizocaryum dogieli Poljansky and Golikova, 1957. Licnophora macfarlandi is an endosymbiont of the respiratory trees of holothuroids, and S. dogieli is found in the esophagus of echinoids. A third species, Lechriopyla mystax Lynch, 1930, is a plagiopylid ciliate found in the intestine of echinoids. Host echinoderms were collected near the Friday Harbor Laboratories, San Juan Island, WA. Specimens of S. dogieli and L. mystax were obtained from the esophagus and intestine, respectively, of the sea urchin Strongylocentrotus pallidus. Specimens of L. macfarlandi were collected from the fluid obtained from the respiratory trees of Parastichopus californicus. Using small subunit ribosomal RNA (SSrRNA) sequences of these three ciliates and a global alignment of SSrRNA sequences of other ciliates, we established the following. 1) Licnophora is a spirotrich ciliate, clearly related to the hypotrichs and stichotrichs; this is corroborated by its possession of macronuclear replication bands. 2) Lechriopyla is the sister genus to Plagiopyla and is a member of the Class Plagiopylea, which was predicted based on its cytology. 3) Schizocaryum clusters in the Class Oligohymenophorea and is most closely related to the scuticociliates; there are currently no morphological features known to relate Schizocaryum to the scuticociliates.

A molecular approach to the tangled intrageneric relationships underlying phylogeny in Euplotes (Ciliophora, Spirotrichea)

The cosmopolitan genus Euplotes is remarkable among ciliates for its species richness. To understand the still tangled taxonomy and phylogenetic relationships within the genus, small subunit rRNA sequences of 11 morphologically defined species colonizing different habitats were determined. Euplotes 18S rRNA is unique among ciliates for its anomalous length and high evolutionary rate. Phylogenetic reconstruction pointed to a high divergence between this genus and the other Spirotrichea, together with a high variability within the genus. Some of the relationships within the Euplotes group were also resolved. Data from the literature, based on morphological features, habitat, and symbiotic relationships, have been compared with our results and are critically discussed. In many cases, the molecular phylogenetic analysis disagreed with species relatedness established on morphological and ecological grounds. The occurrence of a radiation phenomenon in the evolution of the genus is postulated.

Nuclear and cortical regulation in doublets of Paramecium: II. When and how do two cortical domains reorganize to one?

Homopolar doublets with twofold rotational symmetry were generated in Paramecium tetraurelia and in P. undecaurelia by electrofusion or by arrested conjugation. These doublets underwent a complex cortical reorganization over time, which led to their reversion to singlets. This reorganization involved a reduction in number of ciliary rows, a progressive inactivation and loss of one oral meridian, and a reduction and eventual disappearance of one cortical surface (semicell) situated between the two oral meridians. The intermediate steps of this reorganization included some processes that resemble those previously described in regulating doublets of other ciliates, and others that are peculiar to members of the "P. aurelia" species-group and some of its close relatives. The former included a disappearance of one cortical landmark (a contractile vacuole meridian) and transient appearance of another (a third cytoproct) within the narrower semicell. The latter included a reorganization of the paratene zone and the associated invariant (non-duplicating) region to occupy the entire narrower semicell and a redistribution of zones of most active basal-body proliferation within the opposite, wider semicell. The final steps of reorganization involved anterior displacement, invagination, and resorption of one of the two oral apparatuses and eventual disappearance of the associated oral meridian. An oral meridian deprived of its oral apparatus, either by spontaneous resorption or microsurgical removal, could persist for some time in "incomplete doublets" before regulating to the singlet condition. The phylogenetically widespread events encountered in the regulation of doublets to singlets suggest that Paramecium shares some of the global regulatory properties that are likely to be ancestral in ciliates. The more specific events are probably associated with the complex cytoskeletal architecture of this organism and with the frequent occurrence of autogamy that was described in the preceding study (Prajer et al. 1999).

Phylogenetic relationships within the class Spirotrichea (Ciliophora) inferred from small subunit rRNA gene sequences

The small subunit rDNAs of five species belonging to the Euplotidae and eight species of the Oxytrichidae were sequenced to obtain a more detailed picture of the phylogenetic relationships within the Spirotrichea (Ciliophora). Various tree reconstruction algorhythms yielded nearly identical topologies. All Euplotidae were separated from the other Spirotrichea by a deep split. Further, a large genetic distance between the marine genus Moneuplotes and the freshwater species of Euplotoides was found. Differences between the methods used occurred only within the Oxytrichidae. Whereas the monophyly of the Stylonychinae was supported in all trees, the monophyly of the Oxytrichinae was not. However, the molecular data support the morphological and ontogenetic evidence that the pattern of 18 frontal-ventral-transversal cirri evolved in the stemline of the Oxytrichidae and was modified several times independently. Our results are also in agreement with taxonomic revisions: the separation of both Sterkiella nova from Oxytricha and Tetmemena pustulata from Stylonychia.

Phylogenetic position of Sorogena stoianovitchae and relationships within the class Colpodea (Ciliophora) based on SSU rDNA sequences

The ciliate Sorogena stoianovitchae, which can form a multicellular fruiting body, has been classified based upon its ultrastructure and morphology: the oral and somatic infraciliature of S. stoianovitchae most closely resemble those of members of the order Cyrtolophosidida in the class Colpodea. We characterized the small subunit ribosomal DNA (SSU rDNA) gene sequence from S. stoianovitchae and compared this sequence with those from representatives of all ciliate classes. These analyses placed S. stoianovitchae as either sister to members of the class Nassophorea or Colpodea. In an in-group analysis, including all SSU rDNA sequences from members of the classes Nassophorea and Colpodea and representatives of appropriate outgroups, S. stoianovitchae was always sister to Platyophrya vorax (class Colpodea, order Cyrtolophosidida). However, our analyses failed to support the monophyly of the class Colpodea. Instead, our data suggest that there are essentially three unresolved clades: (1) the class Nassophorea; (2) Bresslaua vorax, Colpoda inflata, Pseudoplatyophrya nana, and Bursaria truncatella (class Colpodea); and (3) P. vorax and S. stoianovitchae (class Colpodea).

Polyglycylation of tubulin is essential and affects cell motility and division in Tetrahymena thermophila

We analyzed the role of tubulin polyglycylation in Tetrahymena thermophila using in vivo mutagenesis and immunochemical analysis with modification-specific antibodies. Three and five polyglycylation sites were identified at glutamic acids near the COOH termini of alpha- and beta-tubulin, respectively. Mutants lacking all polyglycylation sites on alpha-tubulin have normal phenotype, whereas similar sites on beta-tubulin are essential. A viable mutant with three mutated sites in beta-tubulin showed reduced tubulin glycylation, slow growth and motility, and defects in cytokinesis. Cells in which all five polyglycylation sites on beta-tubulin were mutated were viable if they were cotransformed with an alpha-tubulin gene whose COOH terminus was replaced by the wild-type COOH terminus of beta-tubulin. In this double mutant, beta-tubulin lacked detectable polyglycylation, while the alpha-beta tubulin chimera was hyperglycylated compared with alpha-tubulin in wild-type cells. Thus, the essential function of polyglycylation of the COOH terminus of beta-tubulin can be transferred to alpha-tubulin, indicating it is the total amount of polyglycylation on both alpha- and beta-tubulin that is essential for survival.

Characterization of histone H3/H4 gene region and phylogenetic affinity of Ichthyophthirius multifiliis based on H4 DNA sequence variation

We sequenced the amino-terminal third of the histone H3 and H4 genes and the intergenic region from Ichthyophthirius multifiliis. Fourteen recombinant clones of 646 bp were sequenced and the level of sequence variation detected among these clones was similar to that reported among closely related species of Tetrahymena and to levels of sequence variation detected within other ciliates. The intergenic region is 417 bp and approximately 92% AT rich, making it the longest and most AT-rich ciliate H3/H4 intergenic region yet identified. Similar to Tetrahymena, the intergenic region of Ichthyophthirius contains two CCAAT regions arranged in a complementary orientation. A neighbor-joining tree was constructed based on nucleotide sequence variation among H4 genes to evaluate evolutionary relationships within and among six classes of Ciliophora. The single shortest neighbor-joining tree depicted a sister-group relationship of Ichthyophthirius with taxa of Tetrahymenina, thereby supporting monophyly of Oligohymenophorea.

Characterization of the Euplotes crassus macronuclear rDNA and its potential as a DNA transformation vehicle

We have cloned the macronuclear linear DNA molecule carrying the ribosomal RNA genes from the ciliated protozoan Euplotes crassus. DNA sequence analysis was carried out to locate coding regions and to determine whether sequences that have been mutated to confer antibiotic resistance are conserved in the E. crassus genes. The beginning and end of the primary transcript were mapped. In order to determine whether conserved sequences that might serve as replication origins were present, the 5' and 3' non-coding sequences from E. crassus were compared to the corresponding sequences from the macronuclear linear rDNA molecules from the following euplotid species: Euplotes vannus, Euplotes minuta, Euplotes raikovii and Euplotes rariseta. A DNA transformation construct was made by generating a putative anisomycin resistant mutation along with a mutation generating a restriction site polymorphism. Microinjection of the construct into the developing macronucleus of mated cells resulted in exconjugant cell lines with increased resistance to anisomycin. The injected rDNA with the restriction site polymorphism is detectable in the anisomycin resistant cells and appears to represent a minor fraction of the rDNA.

Protein coding gene trees in ciliates: comparison with rRNA-based phylogenies

We have reexamined the phylogeny of the ciliates using alpha-tubulin and phosphoglycerate kinase gene sequences. For alpha-tubulin, we have compared the amino acid and nucleotide sequences of 20 species representing seven of the nine classes of the phylum (Karyorelictea, Heterotrichea, Hypotrichea, Oligohymenophorea, Colpodea, Nassophorea, and Litostomatea). The phylogenetic tree resembles a bush from which three monophyletic lineages can be distinguished which correspond to the three classes Hypotrichea, Oligohymenophorea, and Litostomatea. For phosphoglycerate kinase, we have compared the amino acid sequences from 7 species representing three classes (Heterotrichea, Hypotrichea, and Oligohymenophorea). The branching pattern is resolved in three deeply separated branches with an early emergence of the heterotrich. Our comparative analysis shows that if alpha-tubulin phylogeny is not informative at the interclass level, the preliminary data from the phosphoglycerate kinase molecule appear more promising. Nevertheless, at low taxonomic level and at the class level, the resolved phylogenetic relationships inferred from both protein and rRNA sequence data are congruent.

Phylogenetic relationships between prostome and colpodean ciliates tested by small subunit rRNA sequences

We analyzed the phylogenetic relationships among the four ciliate classes, Prostomatea, Colpodea, Nassophorea, and Litostomatea, by small subunit ribosomal RNA gene sequence comparison to test several phylogenetic hypotheses that had been proposed for these classes, based on ultrastructural features. The complete coding regions of the prostomes, Prorodon teres and Coleps hirtus, and of the colpodean, Bursaria truncatella, were determined and compared to the database and the literature. The prostomes and colpodeans are strongly supported as monophyletic groups in all analyses. A sister-group relationship between prostomes and colpodeans is strongly supported by parsimony analyses and moderately supported by distance matrix analyses. There was no support for a close relationship between colpodeans and nassophoreans or between colpodeans and litostomes. In all analyses, these four classes branch late in the ciliate tree. This finding further supports the new view of ciliate phylogeny that so-called ancestral taxa, such as the "primitive" prostomes, are in fact derived and that so-called "highly evolved" groups, such as the heterotrichs, branch off very early in the ciliate tree.