The extent of protist diversity: insights from molecular ecology of freshwater eukaryotes

Classical studies on protist diversity of freshwater environments worldwide have led to the idea that most species of microbial eukaryotes are known. One exemplary case would be constituted by the ciliates, which have been claimed to encompass a few thousands of ubiquitous species, most of them already described. Recently, molecular methods have revealed an unsuspected protist diversity, especially in oceanic as well as some extreme environments, suggesting the occurrence of a hidden diversity of eukaryotic lineages. In order to test if this holds also for freshwater environments, we have carried out a molecular survey of small subunit ribosomal RNA genes in water and sediment samples of two ponds, one oxic and another suboxic, from the same geographic area. Our results show that protist diversity is very high. The majority of phylotypes affiliated within a few well established eukaryotic kingdoms or phyla, including alveolates, cryptophytes, heterokonts, Cercozoa, Centroheliozoa and haptophytes, although a few sequences did not display a clear taxonomic affiliation. The diversity of sequences within groups was very large, particularly that of ciliates, and a number of them were very divergent from known species, which could define new intra-phylum groups. This suggests that, contrary to current ideas, the diversity of freshwater protists is far from being completely described.     

Protist Diversity and Seasonal Dynamics in Skagerrak Plankton Communities as Revealed by Metabarcoding and Microscopy

Protist community composition and seasonal dynamics are of major importance for the production of higher trophic levels, such as zooplankton and fish. Our aim was to reveal how the protist community in the Skagerrak changes through the seasons by combining high‐throughput sequencing and microscopy of plankton collected monthly over two years. The V4 region of the 18S rRNA gene was amplified by eukaryote universal primers from the total RNA/cDNA. We found a strong seasonal variation in protist composition and proportional abundances, and a difference between two depths within the euphotic zone. Highest protist richness was found in late summer‐early autumn, and lowest in winter. Temperature was the abiotic factor explaining most of the variation in diversity. Dinoflagellates was the most abundant and diverse group followed by ciliates and diatoms. We found about 70 new taxa recorded for the first time in the Skagerrak. The seasonal pattern in relative read abundance of major phytoplankton groups was well in accordance with microscopical biovolumes. This is the first metabarcoding study of the protist plankton community of all taxonomic groups and through seasons in the Skagerrak, which may serve as a baseline for future surveys to reveal effects of climate and environmental changes.     

Protist community composition during spring in an Arctic flaw lead polynya

The overwintering deployment of an icebreaker during the Canadian Flaw Lead study provided an opportunity to evaluate how protist communities (phytoplankton and other single-celled eukaryotes) respond to changing spring irradiance conditions in flaw lead polynyas, where open water persists between the central pack ice and land fast ice. We combined microscopic analysis of the protist communities (all cell sizes) with clone libraries of 18S rRNA genes and 18S rRNA (from RNA converted to cDNA) of size-fractionated seawater (0.2–3.0 μm) from 10 to 12 m depth in the surface mixed layer. The rRNA gene analysis provided information on the presence of organisms, while the rRNA analysis provided information on the most active members of the community. There was little overlap between the two types of clone libraries, and there were large community shifts over time. Heterotrophic dinoflagellates and ciliates were the most common sequences recovered. The relative proportion of photosynthetic protist sequences increased in March and April, and there was greater representation of Bacillariophyta, Prasinophyta, Haptophyta, and Cryptophyta in the rRNA compared to rRNA gene libraries. Microscopy indicated that large-celled diatoms dominated the community in May, when chlorophyll concentrations were greatest. However, the RNA sequencing showed that heterotrophic and putative parasitic protists were proportionately more active, and the concomitant decrease in nutrients suggested that the spring phytoplankton bloom had begun to decline by this time. These observations provide evidence of substantial changes in protist community structure and function during the spring transition.     

Multiple genes of apparent algal origin suggest ciliates may once have been photosynthetic

Plantae (as defined by Cavalier-Smith, 1981) plastids evolved via primary endosymbiosis whereby a heterotrophic protist enslaved a photosynthetic cyanobacterium. This "primary" plastid spread into other eukaryotes via secondary endosymbiosis. An important but contentious theory in algal evolution is the chromalveolate hypothesis that posits chromists (cryptophytes, haptophytes, and stramenopiles) and alveolates (ciliates, apicomplexans, and dinoflagellates) share a common ancestor that contained a red-algal-derived "secondary" plastid. Under this view, the existence of several later-diverging plastid-lacking chromalveolates such as ciliates and oomycetes would be explained by plastid loss in these lineages. To test the idea of a photosynthetic ancestry for ciliates, we used the 27,446 predicted proteins from the macronuclear genome of Tetrahymena thermophila to query prokaryotic and eukaryotic genomes. We identified 16 proteins of possible algal origin in the ciliates Tetrahymena and Paramecium tetraurelia. Fourteen of these are present in other chromalveolates. Here we compare and contrast the likely scenarios for algal-gene origin in ciliates either via multiple rounds of horizontal gene transfer (HGT) from algal prey or symbionts, or through endosymbiotic gene transfer (EGT) during a putative photosynthetic phase in their evolution.     

Chrysophytes and other protists in High Arctic lakes: molecular gene surveys,pigment signatures and microscopy

Microscopic analysis of the phytoplankton and other protist communities in High Arctic lakes has shown that they often contain taxa in the Chrysophyceae. Such studies have been increasingly supported by pigment analysis using high-performance liquid chromatography (HPLC) to identify the major algal groups. However, the use of 18S rRNA gene surveys in other systems indicates that many protists, especially small heterotrophs, are underreported or missed by microscopy and HPLC. Here, we investigated the late summer protist community structure of three contrasting lakes in High Arctic polar desert catchments (Char Lake at 74°42′ N, Lake A at 83°00′ N and Ward Hunt Lake at 83°05′ N) with a combination of microscopy, pigment analysis and small subunit 18S ribosomal RNA gene surveys. All three methods showed that chrysophytes were well represented, accounting for 50–70% of total protist community biomass and 25–50% of total 18S rRNA gene sequences. HPLC analysis supported these observations by showing the ubiquitous presence of chrysophyte pigments. The clone libraries revealed a greater contribution of heterotrophs to the protist communities than suggested by microscopy. The flagellate Telonema and ciliates were common in all three lakes, and one fungal sequence was recovered from Char Lake. The approaches yielded complementary information about the protist community structure in the three lakes and underscored the importance of chrysophytes, suggesting that they are well adapted to cope with the low nutrient supply and strong seasonality that characterize the High Arctic environment.     

Comparative analysis between protist communities from the deep‐sea pelagic ecosystem and specific deep hydrothermal habitats

Protist communities associated with deep seawater and bivalves from six hydrothermal sites in the Pacific Ocean were characterized by microscopy and molecular rRNA gene surveys (18S rRNA) and compared with planktonic communities from Pacific deep‐pelagic seawater (from 500 to 3000 m in depth). Genetic libraries from larger size fractions (> 3 µm) of deep‐pelagic water were mainly dominated by Dinophyceae, whereas small size fractions (< 3 µm) mainly revealed radiolarians and Syndiniales. In contrast, more specific opportunistic detritivores and grazers, mostly belonging to Stramenopiles and Cercozoa, were detected from water surrounding vent chimneys. Protist communities were different in the pallial cavity of the giant hydrothermal bivalves Bathymodiolus thermophilus and Calyptogena magnifica, dominated by Ciliophora (primarily belonging to Phyllopharyngea, Oligohymenophorea and Oligotrichea) and Cercozoa. Interestingly, protist communities retrieved from the pallial cavity liquid of hydrothermal bivalves were remarkably homogeneous along the Southern East Pacific Rise, in contrast to bivalves collected on the Mid‐Atlantic Ridge hydrothermal vents and cold seeps from the Gulf of Mexico. Hence, complex protist communities seem to occur inside hydrothermal bivalves, and these metazoa may constitute a stable micro‐niche for micro‐eukaryotes, including grazers, detritivores, symbionts and potential parasites. From these communities, new lineages within the ciliates may emerge.     

Interrelationships among major protistan groups based on a parsimony network of 5S rRNA sequences

To test the validity of the maximum parsimony approach to discern protistan interrelationships, we have derived an optimal network of 16S-like rRNA sequences using our parsimony algorithm and compared it with those reported using the distance matrix method. We have also derived an optimal network topology of 50 5S rRNA sequences through an interactive search using our algorithm. In both these networks, the kinetoplastids and euglenoids form a linkage group with Dictyostelium emerging from its neighbourhood. The cryptophytes, dinoflagellates and chromophytes and green algae emerge as independent lines suggesting that plastids arose more than once during protistan evolution. The large 5S rRNA tree further indicates independent origins of mesozoa and metazoa; kinetoplastids and ciliates; and diphyletic origin of fungi. Comparatively close positions of charales and land plants, chytrids and Zygomycetes, Physarum and amoeba, and red algae and green algae are also seen in this network.     

Protist diversity in suboxic and sulfidic waters of the Black Sea

The oxic-anoxic transition zone of the Black Sea comprises a large suboxic zone as well as anoxic and sulfidic waters. While the prokaryotes and biogeochemical cycles that characterize this zone have been frequently studied, little is known about the diversity or ecology of its microbial eukaryotes. Here, we present the first broad qualitative report of the protist species composition in the Black Sea redoxcline using molecular tools. Fingerprint analysis from the whole redoxcline revealed a complex community structure of metabolically active protists with distinct shifts along the redox gradient. Additionally, 18S rRNA gene clone libraries were used to compare protist species composition of suboxic and sulfidic water layers. Among the ciliates, sequences related to Pleuronema and Strombidium were dominant in both water layers whereas sequences affiliated with anaerobic plagiopylids and Cyclidium were detected only in the sulfidic zone. Among the flagellates, mainly stramenopiles (mostly bicosoecids and chrysophytes) occurred throughout the redoxcline. In the sulfidic zone we found stramenopile sequences but also euglenozoans, jakobids and choanoflagellates that were related to clonal sequences from other anoxic marine habitats, thus indicating the existence of globally distributed groups of anoxic flagellates. Higher species diversity in the sulfidic zone and about twice as many novel sequence types of ciliates and stramenopiles compared with the suboxic layer emphasizes the importance of anoxic, sulfidic waters as habitat for high protist diversity although the function of these organisms is yet unknown.     

Early evolution and the origin of eukaryotes

Our understanding of evolutionary relationships in the eukaryotic world has been revolutionized by molecular systematics. Phylogenies based upon comparisons of rRNAs define five major eukaryotic assemblages plus a series of paraphyletic protist lineages. Comparison of conserved genes that were duplicated prior to the divergence of eubacteria, archaebacteria, and eukaryotes, positions the root of the universal tree within the eubacterial line of descent. In this review a novel model is presented which uses the rRNA and protein based phylogenies to describe the evolutionary origins of eukaryotes.     

Relationship between the flagellates and the ciliates.

The flagellates and the ciliates have long been considered to be closely related because of their unicellular nature and the similarity in the structures of the axoneme of the flagella and cilia in both groups. Most protozoologists believe that the ciliates arose from a flagellate. The flagellates that are most similar in structure to the ciliates are the dinoflagellates and two genera of uncertain taxonomic position, Colponema and Katablepharis. Structurally, dinoflagellates have a number of similarities with ciliates. These include the similarity of the cortical alveoli in the ciliates to the thecal vesicles in the dinoflagellates, the possession of tubular cristae, the similarity of the parasomal sac of the ciliates to the pusule of the dinoflagellates, the possession of similar trichocysts and mucocysts, and some similarity in the feeding apparatus. Colponema spp. are probably related to the dinoflagellates and have many of the same similarities with the ciliates. Katablepharis spp. are very similar in structure to the swarmer (embryo) of the suctorian ciliates. Indeed, reduction in the number of cilia to two in the suctorian swarmer and elimination of the macronucleus would result in a cell that is very similar to the Katablepharis cell. The feeding apparatus of Katablepharis spp. and the rest of the ciliates consists of two concentric microtubular arrays associated with vesicles. Information available from nucleotide sequencing of rRNA places the dinoflagellates in an ancestral position to the ciliates. The rRNA of Colponema and Katablepharis spp. has not yet been investigated. The use of stop codons in mRNA is discussed in relation to phylogeny.     

ALGAE CONTAINING CHLOROPHYLLS a + c ARE PARAPHYLETIC: MOLECULAR EVOLUTIONARY ANALYSIS OF THE CHROMOPHYTA

Sequence comparisons of small subunit ribosomal RNA coding regions from 12 chlorophylls a + c-containing algae were used to infer phylogenetic relationships within the Chromophyta. Three chromophyte lines of descent, delineated by the Bacillariophyceae, the Phaeophyceae/Xanthophyceae, and the Chrysophyceae/Eustigmatophyceae/Synurophyceae are members of a complex evolutionary assemblage, which also includes representatives of the Oomycota (“lower” fungi). Maximum parsimony and distance matrix methods demonstrate a common evolutionary history for these lineages but their relative branching order could not be determined. Other algal species with chlorophylls a + c, including dinoflagellates and prymnesiophytes, are not members of this complex assemblage. Dinoflagellates are specifically related to apicomplexans and ciliates, and the prymnesiophyte, Emiliania huxleyi, represents an independent photosynthetic lineage that separated from other eukaryotes during the nearly simultaneous divergence of plants, animals, fungi, and a number of other protist lineages. The small subunit rRNA phylogenies of chromophytes/oomycetes were compared to those derived from comparisons of ultrastructural characters. Only tubular, tripartite mastigonemes (flagellar hairs) characterized all studied taxa of chromophytes/oomycetes as a monophyletic assemblage.     

Unusually high evolutionary rate of the elongation factor 1 alpha genes from the Ciliophora and its impact on the phylogeny of eukaryotes

The elongation factor 1 alpha (EF-1 alpha) has become widely employed as a phylogenetic marker for studying eukaryotic evolution. However, a disturbing problem, the artifactual polyphyly of ciliates, is always observed. It has been suggested that the addition of new sequences will help to circumvent this problem. Thus, we have determined 15 new ciliate EF-1 alpha sequences, providing for a more comprehensive taxonomic sampling of this phylum. These sequences have been analyzed together with a representation of eukaryotic sequences using distance-, parsimony-, and likelihood-based phylogenetic methods. Such analyses again failed to recover the monophyly of Ciliophora. A study of the substitution rate showed that ciliate EF-1 alpha genes exhibit a high evolutionary rate, produced in part by an increased number of variable positions. This acceleration could be related to alterations of the accessory functions acquired by this protein, likely to those involving interactions with the cytoskeleton, which is very modified in the Ciliophora. The high evolutionary rate of these sequences leads to an artificial basal emergence of some ciliates in the eukaryotic tree by effecting a long-branch attraction artifact that produces an asymmetric topology for the basal region of the tree. The use of a maximum-likelihood phylogenetic method (which is less sensitive to long-branch attraction) and the addition of sequences to break long branches allow retrieval of more symmetric topologies, which suggests that the asymmetric part of the tree is most likely artifactual. Therefore, the sole reliable part of the tree appears to correspond to the apical symmetric region. These kinds of observations suggest that the general eukaryotic evolution might have consisted of a massive radiation followed by an increase in the evolutionary rates of certain groups that emerge artificially as early branches in the asymmetric base of the tree. Ciliates in the case of the EF-1 alpha genes would offer clear evidence for this hypothesis.     

Genetic code deviations in the ciliates: evidence for multiple and independent events.

In several species of ciliates, the universal stop codons UAA and UAG are translated into glutamine, while in the euplotids, the glutamine codon usage is normal, but UGA appears to be translated as cysteine. Because the emerging position of this monophyletic group in the eukaryotic lineage is relatively late, this deviant genetic code represents a derived state of the universal code. The question is therefore raised as to how these changes arose within the evolutionary pathways of the phylum. Here, we have investigated the presence of stop codons in alpha tubulin and/or phosphoglycerate kinase gene coding sequences from diverse species of ciliates scattered over the phylogenetic tree constructed from 28S rRNA sequences. In our data set, when deviations occur they correspond to in frame UAA and UAG coding for glutamine. By combining these new data with those previously reported, we show that (i) utilization of UAA and UAG codons occurs to different extents between, but also within, the different classes of ciliates and (ii) the resulting phylogenetic pattern of deviations from the universal code cannot be accounted for by a scenario involving a single transition to the unusual code. Thus, contrary to expectations, deviations from the universal genetic code have arisen independently several times within the phylum.     

Endosymbiotic Bacteroidales bacteria of the flagellated protist Pseudotrichonympha grassii in the gut of the termite Coptotermes formosanus

A unique lineage of bacteria belonging to the order Bacteroidales was identified as an intracellular endosymbiont of the protist Pseudotrichonympha grassii (Parabasalia, Hypermastigea) in the gut of the termite Coptotermes formosanus. We identified the 16S rRNA, gyrB, elongation factor Tu, and groEL gene sequences in the endosymbiont and detected a very low level of sequence divergence (<0.9% of the nucleotides) in the endosymbiont population within and among protist cells. The Bacteroidales endosymbiont sequence was affiliated with a cluster comprising only sequences from termite gut bacteria and was not closely related to sequences identified for members of the Bacteroidales attached to the cell surfaces of other gut protists. Transmission electron microscopy showed that there were numerous rod-shaped bacteria in the cytoplasm of the host protist, and we detected the endosymbiont by fluorescence in situ hybridization (FISH) with an oligonucleotide probe specific for the 16S rRNA gene identified. Quantification of the abundance of the Bacteroidales endosymbiont by sequence-specific cleavage of rRNA with RNase H and FISH cell counting revealed, surprisingly, that the endosymbiont accounted for 82% of the total bacterial rRNA and 71% of the total bacterial cells in the gut community. The genetically nearly homogeneous endosymbionts of Pseudotrichonympha were very abundant in the gut symbiotic community of the termite.     

Environmental drivers of plankton protist communities along latitudinal and vertical gradients in the oldest and deepest freshwater lake

Identifying which abiotic and biotic factors determine microbial community assembly is crucial to understand ecological processes and predict how communities will respond to environmental change. While global surveys aim at addressing this question in the world's oceans, equivalent studies in large freshwater systems are virtually lacking. Being the oldest, deepest and most voluminous freshwater lake on Earth, Lake Baikal offers a unique opportunity to test the effect of horizontal versus vertical gradients in community structure. Here, we characterized the structure of planktonic microbial eukaryotic communities (0.2–30 μm cell size) along a North–South latitudinal gradient (~600 km) from samples collected in coastal and pelagic waters and from surface to the deepest zones (5–1400 m) using an 18S rRNA gene metabarcoding approach. Our results show complex and diverse protist communities dominated by alveolates (ciliates and dinoflagellates), ochrophytes and holomycotan lineages, with cryptophytes, haptophytes, katablepharids and telonemids in moderate abundance and many low-frequency lineages, including several typical marine members, such as diplonemids, syndinians and radiolarians. Depth had a strong significant effect on protist community stratification. By contrast, the effect of the latitudinal gradient was marginal and no significant difference was observed between coastal and surface open water communities. Co-occurrence network analyses showed that epipelagic communities were significantly more interconnected than communities from the dark water column and suggest specific biotic interactions between autotrophic, heterotrophic and parasitic lineages that influence protist community structure. Since climate change is rapidly affecting Siberia and Lake Baikal, our comprehensive protist survey constitutes a useful reference to monitor ongoing community shifts.     

Ellobiopsids of the genus Thalassomyces are alveolates

Ellobiopsids are multinucleate protist parasites of aquatic crustaceans that possess a nutrient absorbing 'root' inside the host and reproductive structures that protrude through the carapace. Ellobiopsids have variously been affiliated with fungi, 'colorless algae', and dinoflagellates, although no morphological character has been identified that definitively allies them with any particular eukaryotic lineage. The arrangement of the trailing and circumferential flagella of the rarely observed bi-flagellated 'zoospore' is reminiscent of dinoflagellate flagellation, but a well-organized 'dinokaryotic nucleus' has never been observed. Using small subunit ribosomal RNA gene sequences from two species of Thalassomyces, phylogenetic analyses robustly place these ellobiopsid species among the alveolates (ciliates, apicomplexans, dinoflagellates and relatives) though without a clear affiliation to any established alveolate lineage. Our trees demonstrate that Thalassomyces fall within a dinoflagellate + apicomplexa + Perkinsidae + "marine alveolate group 1" clade, clustering most closely with dinoflagellates. However, the poor statistical support for branches within this region indicates that additional data will be needed to resolve relationships among these taxa.     

Probing protistan phylogenies with an anti-tubulin antibody

Representatives from most of the protist phyla were probed by immunofluorescence or by immunoblotting with an anti-tubulin antibody of sharp specificity previously raised against Paramecium axonemal tubulin. Excellent intra-phylum homogeneity of results was recorded except for chlorophytes. All ciliates, dinoflagellates and cryptomonads tested were strongly positive while actinopods, Euglenozoa and parabasalids were negative. All representatives of the broad chromophyte assemblage were positive while all rhizopods were negative. This simple binary immunological character was superimposed on a number of published protist phylogenies and seen to fit very well with some of them. Other immunological approaches to protist taxonomy and evolution are briefly reviewed.     

Alpha‐Tubulin and Small Subunit rRNA Phylogenies of Peritrichs Are Congruent and Do Not Support the Clustering of Mobilids and Sessilids (Ciliophora,Oligohymenophorea)

ABSTRACT. 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 α‐tubulin phylogeny. For this purpose, we obtained 10 partial α‐tubulin sequences from peritrichs and built phylogenetic trees based on α‐tubulin nucleotide and amino acid data. A phylogenetic tree from the α‐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 α‐tubulin gene trees, the urceolariids and trichodiniids do not group together, suggested mobilids may not be a monophyletic group.