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

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.

Study on Analysis of Several Molecular Identification Methods for Ciliates of Colpodea (Protista, Ciliophora)

The application of molecular techniques to accurately identify protozoan species can correct previous misidentifications based on traditional morphological identification. Colpodea ciliates have many toxicological and cytological applications, but their subtle morphological differences and small body size hinder species delineation. Herein, we used Cox I and β-tubulin genes, alongside fluorescence in situ hybridization (FISH), to evaluate each method in delineating Colpodea species. For this analysis, Colpoda harbinensis n. sp., C. reniformis, two populations of C. inflata, Colpoda compare grandis, and five populations of Paracolpoda steinii, from the soil in northeastern China, were used. We determined that (1) the Cox I gene was more suitable than the β-tubulin gene as a molecular marker for defining intra- and interspecific level relationships of Colpoda. (2) FISH probes designed for Colpoda sp., C. inflata, Colpoda compare grandis, and Paracolpoda steinii, provided rapid interspecific differentiation of Colpodea species. (3) Colpoda harbinensis n. sp. was established and mainly characterized by its size in vivo (approximately $80\times 60\mu \text{m}$ ), a reniform body in outline, one macronucleus, its spherical shape, a sometimes nonexistent micronucleus, 11–15 somatic kineties, and five or six postoral kineties. In conclusion, combining oligonucleotide probes, DNA barcoding, and morphology for the first time, we have greatly improved the delineation of Colpodea and confirmed that Cox I gene was a promising DNA barcoding marker for species of Colpodea, and FISH could provide useful morphological information as complementing traditional techniques such as silver carbonate.

Selection and paucity of phylogenetic signal challenge the utility of alpha-tubulin in reconstruction of evolutionary history of free-living litostomateans (Protista, Ciliophora)

The class Litostomatea represents a highly diverse but monophyletic group, uniting both free-living and endosymbiotic ciliates. Ribosomal RNA genes and ITS-region sequences helped to recognize and define the main litostomatean lineages, but did not provide enough phylogenetic signal to unambiguously resolve their interrelationships. In this study, we attempted to improve the resolution among main free-living predatory lineages by adding the gene coding for alpha-tubulin. However, our phylogenetic analyses challenged the performance of alpha-tubulin in reconstruction of evolutionary history of free-living litostomateans. We identified several mutually interconnected problems associated with the ciliate alpha-tubulin gene: the paucity of phylogenetic signal, molecular homoplasies and non-neutral evolution. Positive selection may generate molecular homoplasies (parallel evolution), while negative selection may cause a small number of changes and hence little phylogenetic informativness. Both problems were encountered in nucleotide and amino acid alpha-tubulin alignments, indicating an action of various selective pressures. Taking into account the involvement of alpha-tubulin in many essential biological processes, this protein could be so strongly affected by purifying selection that it even might have become an inappropriate molecular marker for reconstruction of phylogenetic relationships. Therefore, a great caution should be paid when tubulin genes are included in phylogenetic and/or phylogenomic analyses.

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.

Novel insights into the genetic diversity of Balantidium and Balantidium-like cyst-forming ciliates

Balantidiasis is considered a neglected zoonotic disease with pigs serving as reservoir hosts. However, Balantidium coli has been recorded in many other mammalian species, including primates. Here, we evaluated the genetic diversity of B. coli in non-human primates using two gene markers (SSrDNA and ITS1-5.8SDNA-ITS2). We analyzed 49 isolates of ciliates from fecal samples originating from 11 species of captive and wild primates, domestic pigs and wild boar. The phylogenetic trees were computed using Bayesian inference and Maximum likelihood. Balantidium entozoon from edible frog and Buxtonella sulcata from cattle were included in the analyses as the closest relatives of B. coli, as well as reference sequences of vestibuliferids. The SSrDNA tree showed the same phylogenetic diversification of B. coli at genus level as the tree constructed based on the ITS region. Based on the polymorphism of SSrDNA sequences, the type species of the genus, namely B. entozoon, appeared to be phylogenetically distinct from B. coli. Thus, we propose a new genus Neobalantidium for the homeothermic clade. Moreover, several isolates from both captive and wild primates (excluding great apes) clustered with B. sulcata with high support, suggesting the existence of a new species within this genus. The cysts of Buxtonella and Neobalantidium are morphologically indistinguishable and the presence of Buxtonella-like ciliates in primates opens the question about possible occurrence of these pathogens in humans.

Assessing whether alpha-tubulin sequences are suitable for phylogenetic reconstruction of Ciliophora with insights into its evolution in euplotids

The current understanding of ciliate phylogeny is mainly based on analyses of a single gene, the small subunit ribosomal RNA (SSU-rDNA). However, phylogenetic trees based on single gene sequence are not reliable estimators of species trees, and SSU-rDNA genealogies are not useful for resolution of some branches within Ciliophora. Since congruence between multiple loci is the best tool to determine evolutionary history, we assessed the usefulness of alpha-tubulin gene, a protein-coding gene that is frequently sequenced, for ciliate phylogeny. Here, we generate alpha-tubulin gene sequences of 12 genera and 30 species within the order Euplotida, one of the most frequently encountered ciliate clades with numerous apparently cosmopolitan species, as well as four genera within its putative sister order Discocephalida. Analyses of the resulting data reveal that: 1) the alpha-tubulin gene is suitable phylogenetic marker for euplotids at the family level, since both nucleotide and amino acid phylogenies recover all monophyletic euplotid families as defined by both morphological criteria and SSU-rDNA trees; however, alpha-tubulin gene is not a good marker for defining species, order and subclass; 2) for seven out of nine euplotid species for which paralogs are detected, gene duplication appears recent as paralogs are monophyletic; 3) the order Euplotida is non-monophyletic, and the family Uronychiidae with sequences from four genera, is non-monophyletic; and 4) there is more genetic diversity within the family Euplotidae than is evident from dargyrome (geometrical pattern of dorsal “silverline system” in ciliates) patterns, habit and SSU-rDNA phylogeny, which indicates the urgent need for taxonomic revision in this area.

Evolutionary dynamics of a conserved sequence motif in the ribosomal genes of the ciliate Paramecium

Background

In protozoa, the identification of preserved motifs by comparative genomics is often impeded by difficulties to generate reliable alignments for non-coding sequences. Moreover, the evolutionary dynamics of regulatory elements in 3' untranslated regions (both in protozoa and metazoa) remains a virtually unexplored issue.

Results

By screening Paramecium tetraurelia's 3' untranslated regions for 8-mers that were previously found to be preserved in mammalian 3' UTRs, we detect and characterize a motif that is distinctly conserved in the ribosomal genes of this ciliate. The motif appears to be conserved across Paramecium aurelia species but is absent from the ribosomal genes of four additional non-Paramecium species surveyed, including another ciliate, Tetrahymena thermophila. Motif-free ribosomal genes retain fewer paralogs in the genome and appear to be lost more rapidly relative to motif-containing genes. Features associated with the discovered preserved motif are consistent with this 8-mer playing a role in post-transcriptional regulation.

Conclusions

Our observations 1) shed light on the evolution of a putative regulatory motif across large phylogenetic distances; 2) are expected to facilitate the understanding of the modulation of ribosomal genes expression in Paramecium; and 3) reveal a largely unexplored--and presumably not restricted to Paramecium--association between the presence/absence of a DNA motif and the evolutionary fate of its host genes.

Alpha-tubulin and small subunit rRNA phylogenies of peritrichs are congruent and do not support the clustering of mobilids and sessilids (Ciliophora, Oligohymenophorea)

Peritrich ciliates have been traditionally subdivided into two orders, Sessilida and Mobilida within the subclass Peritrichia. However, all the existing small subunit (SSU) rRNA phylogenetic trees showed that the sessilids and mobilids did not branch together. To shed some light on this disagreement, we tested whether or not the classic Peritrichia is a monophyletic group by assessing the reliability of the SSU rRNA phylogeny in terms of congruency with alpha-tubulin phylogeny. For this purpose, we obtained 10 partial alpha-tubulin sequences from peritrichs and built phylogenetic trees based on alpha-tubulin nucleotide and amino acid data. A phylogenetic tree from the alpha-tubulin and SSU rRNA genes in combination was also constructed and compared with that from the SSU rRNA gene using a similar species sampling. Our results show that the mobilids and sessilids are consistently separated in all trees, which reinforces the idea that the peritrichs do not constitute a monophyletic group. However, in all alpha-tubulin gene trees, the urceolariids and trichodiniids do not group together, suggested mobilids may not be a monophyletic group.

Relationship between two tandemly arranged and light-induced glutathione S-transferase genes from the ciliated protozoa Blepharisma japonicum

Recently we reported a light-induced cDNA encoding glutathione S-transferase (GST) from the ciliated protozoa Blepharisma japonicum, which possessed photosensitive pigments. In this study, a novel cDNA encoding GST was further isolated, and the two GSTs (BjGST1 and BjGST2) showed high sequence identity of 86%. Phylogenetic trees indicated that the BjGSTs were distantly related to known classes of GSTs, and they could form a protozoa-specific class. The recombinant proteins also existed as homo- or heterodimers that exhibited different enzyme activities, appreciating the functional differentiation. Furthermore, the transcription levels of BjGST genes were coordinately regulated in response to light stimulation. In addition, the genomic structure analysis revealed that the two genes were tandemly arranged through an approximately 500-bp spacer region of unusual DNA structure containing cis-acting elements related to oxidative stress response. These results demonstrate that the two BjGSTs are expressed simultaneously and act cooperatively against photooxidative stress.

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.

Reevaluation of the phylogenetic relationship between mobilid and sessilid peritrichs (Ciliophora, Oligohymenophorea) based on small subunit rRNA genes sequences

Based on morphological characters, peritrich ciliates (Class Olygohymenophorea, Subclass Peritrichia) have been subdivided into the Orders Sessilida and Mobilida. Molecular phylogenetic studies on peritrichs have been restricted to members of the Order Sessilida. In order to shed more light into the evolutionary relationships within peritrichs, the complete small subunit rRNA (SSU rRNA) sequences of four mobilid species, Trichodina nobilis, Trichodina heterodentata, Trichodina reticulata, and Trichodinella myakkae were used to construct phylogenetic trees using maximum parsimony, neighbor joining, and Bayesian analyses. Whatever phylogenetic method used, the peritrichs did not constitute a monophyletic group: mobilid and sessilid species did not cluster together. Similarity in morphology but difference in molecular data led us to suggest that the oral structures of peritrichs are the result of evolutionary convergence. In addition, Trichodina reticulata, a Trichodina species with granules in the center of the adhesive disc, branched separately from its congeners, Trichodina nobilis and Trichodina heterodentata, trichodinids without such granules. This indicates that granules in the adhesive disc might be a phylogenetic character of high importance within the Family Trichodinidae.

Phylogenetic Relationships of the Subclass Peritrichia (Oligohymenophorea, Ciliophora) Inferred from Small Subunit rRNA Gene Sequences1

The phylogenetic relationships among peritrichs remain unresolved. In this study, the complete small subunit rRNA (SSrRNA) gene sequences of seven species (Epistylis galea, Campanella umbellaria, Carchesium polypinum, Zoothamnium arbuscula, Vaginicola crystallina, Ophrydium versatile, and Opercularia microdiscum) were determined. Trees were constructed using distance-matrix, maximum-likelihood and maximum-parsimony methods, all of which strongly supported the monophyly of the subclass Peritrichia. Within the peritrichs, 1) E. galea grouped with Opercularia microdiscum and Campanella umbellaria but not the other Epistylis species, which indicates that the genus Epistylis might not be monophyletic; 2) the topological position of Carchesium and Campanella suggested that Carchesium should be placed in the family Zoothamniidae, or be elevated to a higher taxonomic rank, and that Campanella should be independent of the family Epistylididae, and probably be given a new rank; and 3) Opisthonecta grouped strongly with Astylozoon, which suggested that Opisthonecta species were not the ancestors of the stalked peritrichs.

A rapid and simple method of detection of Blepharisma japonicum using PCR and immobilisation on FTA paper

Background

The rapid expansion in the availability of genome and DNA sequence information has opened up new possibilities for the development of methods for detecting free-living protozoa in environmental samples. The protozoan Blepharisma japonicum was used to investigate a rapid and simple detection system based on polymerase chain reaction amplification (PCR) from organisms immobilised on FTA paper.

Results

Using primers designed from the α-tubulin genes of Blepharisma, specific and sensitive detection to the equivalent of a single Blepharisma cell could be achieved. Similar detection levels were found using water samples, containing Blepharisma, which were dried onto Whatman FTA paper.

Conclusion

This system has potential as a sensitive convenient detection system for Blepharisma and could be applied to other protozoan organisms.

A new noncanonical nuclear genetic code: translation of UAA into glutamate

Deviant genetic codes reported in ciliates share the same feature: one (UGA) or two (UAR) of the three canonical stop codons are translated into one particular amino acid. In many genera, such as Oxytricha, Paramecium, and Tetrahymena, UAR codons are translated into glutamine. UGA is translated into cysteine in Euplotes or into tryptophan in Colpoda inflata and Blepharisma americanum. Here, we show that three peritrich species (Vorticella microstoma, Opisthonecta henneguyi, and Opisthonecta matiensis) translate UAA into glutamate and that at least UAA in O. matiensis is decoded through a mutant suppressor-like tRNA. This kind of genetic code has never been reported for any living organism. Phylogenetic analysis with alpha-tubulin sequences corroborates that peritrichs, peniculines (Paramecium), and hymenostomates (Tetrahymena) form a monophyletic group (class Oligohymenophorea). The differential translation (glu/gln) of UAR codons, the monophyly of the Oligohymenophorea, and the common evolutionary origin of glutamate and glutamine suggest that deviant genetic codes of present-day oligohymenophoreans could have the same origin.

"Novel cytoskeletal proteins in protists": introductory remarks

Protists provide the opportunity to integrate analyses from a low (molecular) to a high (organism) level of complexity within a broad evolutionary framework. The perpectives they offer in the cytoskeletal field are discussed with respect to emerging concepts of cellular biology.

Evolution of duplicated alpha-tubulin genes in ciliates

Ciliates provide a powerful system to analyze the evolution of duplicated alpha-tubulin genes in the context of single-celled organisms. Genealogical analyses of ciliate alpha-tubulin sequences reveal five apparently recent gene duplications. Comparisons of paralogs in different ciliates implicate differing patterns of substitutions (e.g., ratios of replacement/synonymous nucleotides and radical/conservative amino acids) following duplication. Most substitutions between paralogs in Euplotes crassus, Halteria grandinella and Paramecium tetraurelia are synonymous. In contrast, alpha-tubulin paralogs within Stylonychia lemnae and Chilodonella uncinata are evolving at significantly different rates and have higher ratios of both replacement substitutions to synonymous substitutions and radical amino acid changes to conservative amino acid changes. Moreover, the amino acid substitutions in C. uncinata and S. lemnae paralogs are limited to short stretches that correspond to functionally important regions of the alpha-tubulin protein. The topology of ciliate alpha-tubulin genealogies are inconsistent with taxonomy based on morphology and other molecular markers, which may be due to taxonomic sampling, gene conversion, unequal rates of evolution, or asymmetric patterns of gene duplication and loss.

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.

Phylogenetic relationships of the subclass peritrichia (Oligohymenophorea, Ciliophora) with emphasis on the genus Epistylis, inferred from small subunit rRNA gene sequences

The peritrichs have been recognized as a higher taxon of ciliates since 1968. However, the phylogenetic relationships among them are still unsettled, and their placement within the class Oligohymenophorea has only been supported by the analysis of the small subunit rRNA gene sequence of Opisthonecta henneguyi. DNA was isolated directly from field-sampled species for PCR, and was used to resolve relationships within the genus Epistylis and to confirm the stability of the placement of peritrichs. Small subunit rRNA gene sequences of Epistylis plicatilis, Epistylis urceolata, Epistylis chrysemydis, Epistylis hentscheli, Epistylis wenrichi, and Vorticella campanula were sequenced and analyzed using both distance-matrix and maximum-parsimony methods. In phylogenetic trees, the monophyly of both the genus Episrylis and the subclass Peritrichia was strongly supported, while V. campanula clustered with Vorticella microstoma. The topology in which E. plicatilis and E. hentscheli formed a strongly supported sister clade to E. urceolata, E. chrysemydis, and E. wenrichi was consistent with variations in the thickness of the peristomial lip. We concluded that the peristomial area, especially the peristomial lip, might be the important phylogenetic character within the genus Epistylis.

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).

The molecular basis of nuclear genetic code change in ciliates

BACKGROUND

The nuclear genetic code has changed in several lineages of ciliates. These changes, UAR to glutamine and UGA to cysteine, imply that eukaryotic release factor 1 (eRF1), the protein that recognizes stop codons and terminates translation, changes specificity. Here we test whether changes in eRF1 drive genetic code evolution.

RESULTS

Database sequence analysis reveals numerous genetic code alterations in ciliates, including UGA --> tryptophan in Blepharisma americanum and the distantly related Colpoda. We sequenced eRF1 from four ciliates: B. americanum, a heterotrich that independently derived the same eRF1 specificity as Euplotes, and three spirotrichs, Stylonychia lemnae, S. mytilus, and Oxytricha trifallax, that independently derived the same genetic code as Tetrahymena (UAR --> glutamine). Distantly related ciliates with similar codes show characteristic changes in eRF1. We used a sliding window analysis to test associations between changes in specific eRF1 residues and changes in the genetic code. The regions of eRF1 that display convergent substitutions are identical to those identified in a recently reported nonsense suppression mutant screen in yeast.

CONCLUSIONS

Genetic code change by stop codon reassignment is surprisingly frequent in ciliates, with UGA --> tryptophan occurring twice independently. This is the first description of this code, previously found only in bacteria and mitochondria, in a eukaryotic nuclear genome. eRF1 has evolved strikingly convergently in lineages with variant genetic codes. The strong concordance with biochemical data indicates that our methodology may be generally useful for detecting molecular determinants of biochemical changes in evolution.

Phylogenetic relationships of the Subclass Peniculia (Oligohymenophorea, Ciliophora) inferred from small subunit rRNA gene sequences

Peniculine ciliates have been recognized as a distinct higher taxon of ciliates for almost 50 years. However, phylogenetic relationships within the Subclass Peniculia are still unsettled. To contribute to our understanding of their phylogeny and provide evidence for the position of Urocentrum turbo, we sequenced its small subunit (SS) rRNA gene and the SSrRNA genes from Lembadion bullinum, Frontonia sp., Paramecium caudatum, Paramecium multimicronucleatum, Paramecium putrinum, and Paramecium woodruffi. Urocentrum turbo was the only one of these species not to exhibit a shortened Helix E10_1, which we conclude characterizes the "higher" peniculines. Except for U. turbo, the peniculines are strongly supported as a monophyletic clade with Lembadion, Frontonia, and Paramecium species forming separate and strongly supported clades by bootstrap analysis using distance matrix, maximum parsimony, and maximum likelihood methods. Urocentrum turbo is associated with different lineages, depending upon the analysis used. The Paramecium species form at least four clades with the Paramecium aurelia subgroup being the most derived. We conclude that the Subclass Peniculia should be divided into two orders, the Order Urocentrida and Order Peniculida, with the latter order having two suborders, the Suborder Frontoniina and Peniculina. We place U. turbo with the peniculines because of shared morphological and stomatogenetic features.