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

Phylogeny and species delimitation of ciliates in the genus Spirostomum (Class, Heterotrichea) using single-cell transcriptomes

Ciliates are single-celled microbial eukaryotes that diverged from other eukaryotic lineages over a billion years ago. The extensive evolutionary timespan of ciliate has led to enormous genetic and phenotypic changes, contributing significantly to their high level of diversity. Recent analyses based on molecular data have revealed numerous cases of cryptic species complexes in different ciliate lineages, demonstrating the need for a robust approach to delimit species boundaries and elucidate phylogenetic relationships. Heterotrich ciliate species of the genus Spirostomum are abundant in freshwater and brackish environments and are commonly used as biological indicators for assessing water quality. However, some Spirostomum species are difficult to identify due to a lack of distinguishable morphological characteristics, and the existence of cryptic species in this genus remains largely unexplored. Previous phylogenetic studies have focused on only a few loci, namely the ribosomal RNA genes, alpha-tubulin, and mitochondrial CO1. In this study, we obtained single-cell transcriptome of 25 Spirostomum species populations (representing six morphospecies) sampled from South Korea and the USA, and used concatenation- and coalescent-based methods for species tree inference and delimitation. Phylogenomic analysis of 37 Spirostomum populations and 265 protein-coding genes provided a robustious insight into the evolutionary relationships among Spirostomum species and confirmed that species with moniliform and compact macronucleus each form a distinct monophyletic lineage. Furthermore, the multispecies coalescent (MSC) model suggests that there are at least nine cryptic species in the Spirostomum genus, three in S. minus, two in S. ambiguum, S. subtilis, and S. teres each. Overall, our fine sampling of closely related Spirostomum populations and wide scRNA-seq allowed us to demonstrate the hidden crypticity of species within the genus Spirostomum, and to resolve and provide much stronger support than hitherto to the phylogeny of this important ciliate genus.

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.

First Records of 'Flagship' Soil Ciliates in North America

'Flagship' ciliates were investigated from soil samples collected in Florida, USA. This was undertaken to determine if species thought to be restricted to a given world region could be uncovered from similar habitats in a novel location, e.g. another continent. Two species of Condylostomides were discovered, and recorded from the North American continent for the first time. Condylostomides etoschensis was known only from Africa, but was found to be thriving in a Florida study site. An 18S rDNA sequence for this species was determined for the first time. Also discovered from the same study site was the ciliate Condylostomides coeruleus, previously known only from Central and South America. These two 'flagship' ciliates were found in the same habitat, from a continent well outside of their previously recorded biogeographies. Molecular sequencing and microscopy investigations were conducted to form the baseline for future work within this genus. Soil ciliates can obtain large population numbers and form cysts and are therefore likely able to disperse globally. These new records provide additional evidence that large distances, even between continents, do not hinder microbes from thriving globally. The absence of these conspicuously-colored gold and blue ciliates from previous studies is likely due to undersampling, rather than to any physical barriers.

Morphology and Phylogenetic Position of Gruberia lanceolata (Gruber 1884) (Ciliophora, Heterotrichea) from Rio de Janeiro, Brazil

Ciliates of the genus Gruberia are poorly studied. Consequently, most species lack detailed morphological descriptions, and all gene sequences in GenBank are not classified at the species level. In this study, a detailed morphological description of a population of G. lanceolata from Brazil is presented, based on live and protargol-stained organisms. We also present the 18S rRNA gene sequence and the phylogenetic position of this species. The primary characteristics of G. lanceolata from the Maricá Lagoon are as follows: an elongate fusiform body 280-870 × 40-160 μm in size; rosy cortical granules; a peristome occupying approximately 1/3-1/2 of body length; an adoral zone comprising 115-330 membranelles; a paroral membrane in 35-50 fragments; and a moniliform macronucleus with 11-16 nodules. Based on our observations and data from pertinent literature, we suggest G. beninensis to be a junior synonym of G. lanceolata.

New contributions to Gruberia lanceolata (Gruber, 1884) Kahl, 1932 based on analyses of multiple populations and genes (Ciliophora, Heterotrichea, Gruberiidae)

Gruberia Kahl, 1932 is a species-poor genus comprising only seven named species. Most of these species have not been reinvestigated since the original reports. In the present work, we investigated the taxonomy and phylogeny of Gruberia lanceolata (Gruber, 1884) Kahl, 1932 based on analyses of morphology and multiple gene sequences from four South Korean populations. This species is mainly characterized by a well-developed peristome region, segmented paroral membrane, and moniliform macronucleus. Some morphological features were not stable among the four populations investigated, such as body shape and size, cell color, and the ratio of oral length to body length. However, our molecular analyses of four different genetic markers - three nuclear DNA markers (18S rDNA, ITS1-5.8S-ITS2 region, D1D2 of 28S rDNA) and one mitochondrial (mt) marker (CO1 gene) - indicated that all Korean populations examined were the same species. Based on our present findings and historic works, we propose that G. calkinsi, G. aculeata, and G. beninensis are junior synonyms of G. lanceolata.

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.

Morphology and phylogenetic position of an unusual Stentor polymorphus (Ciliophora: Heterotrichea) without symbiotic algae

We investigated the live morphology, infraciliature, and small subunit rRNA gene sequences of an unusual population of Stentor polymorphus without symbiotic algae that was isolated from the southeastern region of Brazil. The morphological and molecular data confirmed the identity of this strain as S. polymorphus. The Brazilian S. polymorphus organism is 850-2,000 μm in length in vivo and has colorless cortical granules, a moniliform macronucleus with 6-12 nodules, somatic ciliature composed of 50-60 kineties, a single contractile vacuole located to the left of the cytostome, and a conspicuous oral pouch, and it does not build a lorica. Based on the phylogenetic analyses, the Brazilian S. polymorphus was located within a cluster consisting of four other S. polymorphus sequences, with high support values using both the Bayesian inference and maximum likelihood algorithms. Our study presents the first report of a S. polymorphus population without its symbionts under natural conditions. On the basis of our findings, we propose that the presence or absence of symbiotic algae should not be used as a taxonomic character for the identification of Stentor species.

Placing environmental next-generation sequencing amplicons from microbial eukaryotes into a phylogenetic context

Nucleotide positions in the hypervariable V4 and V9 regions of the small subunit (SSU)-rDNA locus are normally difficult to align and are usually removed before standard phylogenetic analyses. Yet, with next-generation sequencing data, amplicons of these regions are all that are available to answer ecological and evolutionary questions that rely on phylogenetic inferences. With ciliates, we asked how inclusion of the V4 or V9 regions, regardless of alignment quality, affects tree topologies using distinct phylogenetic methods (including PairDist that is introduced here). Results show that the best approach is to place V4 amplicons into an alignment of full-length Sanger SSU-rDNA sequences and to infer the phylogenetic tree with RAxML. A sliding window algorithm as implemented in RAxML shows, though, that not all nucleotide positions in the V4 region are better than V9 at inferring the ciliate tree. With this approach and an ancestral-state reconstruction, we use V4 amplicons from European nearshore sampling sites to infer that rather than being primarily terrestrial and freshwater, colpodean ciliates may have repeatedly transitioned from terrestrial/freshwater to marine environments.

Morphology and 18S rDNA gene sequence of Blepharisma sinuosum Sawaya, 1940 (Ciliophora: Heterotrichea) from Brazil

The morphology and morphometric data of seven populations of Blepharisma sinuosum from southeastern Brazil were investigated. The description is based on live observations, protargol impregnation, and scanning electron microscopy. Blepharisma sinuosum measures 75-255μm in length and 25-93μm in width and has a spindle-shaped body, pink color, a single contractile vacuole located at the posterior end, 50 adoral membranelles, a conspicuous paroral, 17-35 somatic kineties, a moniliform macronucleus with 2-7 connected nodules, and 3-20 micronuclei. Morphological comparisons with similar species were performed and suggest that B. americanum is the junior synonym of B. sinuosum. The 18S rDNA gene sequence of B. sinuosum was obtained and compared with that of other Blepharisma species. The length and GC content of the obtained sequence is 1652bp and 47.03%, respectively, and has a very high structural similarity (99.9%) with the B. undulans sequence. The validity of the classification of Blepharisma species in morphonuclear subgenera is also discussed.

Morphological and molecular phylogeny of dileptid and tracheliid ciliates: resolution at the base of the class Litostomatea (Ciliophora, Rhynchostomatia)

Dileptid and tracheliid ciliates have been traditionally classified within the subclass Haptoria of the class Litostomatea. However, their phylogenetic position among haptorians has been controversial and indicated that they may play a key role in understanding litostomatean evolution. In order to reconstruct the evolutionary history of dileptids and tracheliids, and to unravel their affinity to other haptorians, we have used a cladistic approach based on morphological evidence and a phylogenetic approach based on 18S rRNA gene sequences, including eight new ones. The molecular trees demonstrate that dileptids and tracheliids represent a separate subclass, Rhynchostomatia, that is sister to the subclasses Haptoria and Trichostomatia. The Rhynchostomatia are characterized by a ventrally located oral opening at the base of a proboscis that carries a complex oral ciliature. We have recognized two orders within Rhynchostomatia. The new order Tracheliida is monotypic, while the order Dileptida comprises two families: the new, typically bimacronucleate family Dimacrocaryonidae and the multimacronucleate family Dileptidae. The Haptoria evolved from the last common ancestor of the Litostomatea by polarization of the body, the oral opening locating more or less apically and the oral ciliature simplifying. The Trichostomatia originated from a microaerophylic haptorian by further simplification of the oral ciliature, possibly due to an endosymbiotic lifestyle.

Linkage between light microscopic observations and molecular analysis by single-cell PCR for ciliates

The molecular analysis of ciliates for identification and phylogeny is usually conducted through PCR amplification and DNA sequencing, with DNA extracted from a large number of cells. Therefore, the analysis of ciliates is limited to only those species that can be cultured. We propose a single-cell PCR procedure to overcome the difficulty in the analysis of unculturable species. The procedure has been tested on 6 Stichotrichia strains and uncultured Levicoleps biwae cells, targeting 18S rRNA gene sequences, after carrying out microscopic observations and obtaining photographic records. This procedure enables the systematic analysis of unculturable ciliate strains in natural environments by linking the morphology and genotype of a single cell.

Small subunit ribosomal RNA and mitochondrial cytochrome c oxidase subunit 1 gene sequences of 21 strains of the parasitic scuticociliate Miamiensis avidus (Ciliophora, Scuticociliatia)

The scuticociliate Miamiensis avidus is a histophagous parasite that causes high mortality in cultured marine fishes. Small subunit ribosomal RNA (SSU rRNA) and mitochondrial cytochrome c oxidase subunit 1 (cox1) genes were analyzed for 21 strains of M. avidus isolated from diseased olive flounder (Paralichthys olivaceus), ridged-eye flounder (Pleuronichthys cornutus), and spotted knifejaw (Oplegnathus fasciatus) in Korea and Japan (collected in 2003-2007). Analysis of SSU rRNA gene sequences (1,759 bp) indicates they are very conserved with less than 0.17% (3 nucleotides) differences suggesting that SSU rRNA are useful to identify M. avidus; however, the cox1 gene (900 bp) has higher variations with intraspecific divergences up to 5.67% (51 nucleotides). A distance tree of cox1 gene sequences based on a neighbor-joining analysis can separate 21 strains into five cox1 types (two heterogeneous clusters and three individual branches). The cox1-type matches with serotype of strains but do not reflect geographical origins, host species, or pathogenicity.

Expanding character sampling for ciliate phylogenetic inference using mitochondrial SSU-rDNA as a molecular marker

Molecular systematics of ciliates, particularly at deep nodes, has largely focused on increasing taxon sampling using the nuclear small subunit rDNA (nSSU-rDNA) locus. These previous analyses have generally been congruent with morphologically-based classifications, although there is extensive non-monophyly at many levels. However, caution is needed in interpreting these results as nSSU-rDNA is just a single molecular marker. Here the mitochondrial small subunit rDNA (mtSSU-rDNA) is evaluated for deep ciliate nodes using the Colpodea as an example. Overall, well-supported nodes in the mtSSU-rDNA and concatenated topologies are well supported in the nSSU-rDNA topology; e.g., the non-monophyly of the Cyrtolophosidida. The two moderately- to well-supported incongruences between the loci are the placement of the Sorogenida andColpoda aspera.Our analyses of mtSSU-rDNA support the conclusion, originally derived from nSSU-rDNA, that the morphological characters used in taxonomic circumscriptions of the Colpodea represent a mixture of ancestral and derived states. This demonstration of the efficacy of the mtSSU-rDNA will enable phylogenetic reconstructions of deep nodes in the ciliate tree of life to move from a single-locus to a multi-locus approach.

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.