Towards a Molecular Taxonomy for Protists: Benefits, Risks, and Applications in Plankton Ecology

Nuclear and Mitochondrial SSU rRNA Genes Reveal Hidden Diversity of Haptophrya Endosymbionts in Freshwater Planarians and Challenge Their Traditional Classification in Astomatia

Like many other aquatic animals, freshwater planarians have also become partners of symbiotic ciliates from the class Oligohymenophorea. In the present study, we explored the hidden diversity and addressed the questionable systematic position of mouthless obligatory gut endosymbionts of freshwater planarians, using the nuclear and mitochondrial SSU rRNA genes. Although all isolated ciliates morphologically corresponded to a single species, molecular analyses suggested the existence of three genetically distinct entities: Haptophrya planariarum, Haptophrya dugesiarum nov. spec., and Haptophrya schmidtearum nov. spec. The two former species share the same planarian host, which indicates a speciation model involving one duplication event without host switching. Such a diversification pattern was recognized also in astome ciliates inhabiting megascolecid and glossoscolecid earthworms. The present multi-gene phylogenies along with the secondary structure of the mitochondrial 16S rRNA molecule, however, challenge the traditional classification of Haptophrya within the subclass Astomatia. Haptophrya very likely evolved from an orphan scuticociliate lineage by the loss of oral apparatus and by the transformation of the thigmotactic field into an adhesive sucker. Since astomy evolved multiple times independently within the Oligohymenophorea, the loss of cell mouth cannot be used as a sole argument for the assignment of Haptophrya to the Astomatia anymore.

Paramecium bursaria—A Complex of Five Cryptic Species: Mitochondrial DNA COI Haplotype Variation and Biogeographic Distribution

Ciliates are a diverse protistan group and many consist of cryptic species complexes whose members may be restricted to particular biogeographic locations. Mitochondrial genes, characterized by a high resolution for closely related species, were applied to identify new species and to distinguish closely related morphospecies. In the current study, we analyzed 132 sequences of COI mtDNA fragments obtained from P. bursaria species collected worldwide. The results allowed, for the first time, to generate a network of COI haplotypes and demonstrate the relationships between P. bursaria strains, as well as to confirm the existence of five reproductively isolated haplogroups. The P. bursaria haplogroups identified in the present study correspond to previously reported syngens (R1, R2, R3, R4, and R5), thus we decided to propose the following binominal names for each of them: P. primabursaria, P. bibursaria, P. tribursaria, P. tetrabursaria, and P. pentabursaria, respectively. The phylogeographic distribution of P. bursaria species showed that P. primabursaria and P. bibursaria were strictly Eurasian, except for two South Australian P. bibursaria strains. P. tribursaria was found mainly in Eastern Asia, in two stands in Europe and in North America. In turn, P. tetrabursaria was restricted to the USA territory, whereas P. pentabursaria was found in two European localities.

Phylogenetic analysis of eight species of Anomopoda based on transcriptomic and mitochondrial DNA sequences

Anomopoda is the widespread planktonic microcrustacean, which plays a crucial role in aquatic ecosystem. There are few studies about the evolutionary relationships among various Anomopoda basing on molecular data. In the present study, phylogenetic analysis of eight Anomopoda was carried out. Firstly, the culture system was developed to breed cladocerans. By using this system, eight species (Daphnia magna, D. pulex, D. sinensis, Ceriodaphnia reticulata, Moina micrura, Scapholeberis kingi, Simocephalus vetulus and Eurycercus lamellatus) were purified and cultured stably in the laboratory. Then, transcriptomic sequences and partial mitochondrial DNA sequences were both used to reconstruct the phylogenetic tree among 8 species. Transcriptomic sequences were sequenced on Illumina Hiseq 2500 platform. After assembly and annotation, transcriptomic sequences were spliced together and aligned for phylogenetic analysis. Basing on the orthologous genes derived from transcriptomic sequences, the phylogenetic analysis showed that 4 genera of Daphniidae were clustered into one group, and among the 4 genera, Ceriodaphnia was closer to Daphnia than Simocephalus, while Scapholeberis was farthest from other species. In addition, Eurycercidae was closer to Daphniidae than Moinidae. The phylogenetic trees based on both 12S rRNA and 16S rRNA sequences were similar with that based on transcriptomic sequences. Meanwhile, the phylogenetic tree based on 16S rRNA sequences was more suitable than that based on 12S rRNA sequences. These results suggested that the phylogenetic analysis basing on the transcriptomic sequences was available in cladocerans, which will help us to effectively understand the phylogenetic relationships among various cladocerans.

A Primer for Microbiome Time-Series Analysis

Time-series can provide critical insights into the structure and function of microbial communities. The analysis of temporal data warrants statistical considerations, distinct from comparative microbiome studies, to address ecological questions. This primer identifies unique challenges and approaches for analyzing microbiome time-series. In doing so, we focus on (1) identifying compositionally similar samples, (2) inferring putative interactions among populations, and (3) detecting periodic signals. We connect theory, code and data via a series of hands-on modules with a motivating biological question centered on marine microbial ecology. The topics of the modules include characterizing shifts in community structure and activity, identifying expression levels with a diel periodic signal, and identifying putative interactions within a complex community. Modules are presented as self-contained, open-access, interactive tutorials in R and Matlab. Throughout, we highlight statistical considerations for dealing with autocorrelated and compositional data, with an eye to improving the robustness of inferences from microbiome time-series. In doing so, we hope that this primer helps to broaden the use of time-series analytic methods within the microbial ecology research community.

Species delimitation for the molecular taxonomy and ecology of the widely distributed microbial eukaryote genus Euplotes (Alveolata, Ciliophora)

Recent advances in high-throughput sequencing and metabarcoding technologies are revolutionizing our understanding of the diversity and ecology of microbial eukaryotes (protists). The interpretation of protist diversity and the elucidation of their ecosystem function are, however, impeded by problems with species delimitation, especially as it applies to molecular taxonomy. Here, using the ciliate Euplotes as an example, we describe approaches for species delimitation based on integrative taxonomy by using evolutionary and ecological perspectives and selecting the most appropriate metabarcoding gene markers as proxies for species units. Our analyses show that: Euplotes (sensu lato) comprises six distinct clades, mainly as result of ecological speciation; the validity of the genera Euplotes (sensu stricto), Euplotoides, Euplotopsis and Moneuplotes are not supported; the vannus-type group, which includes species without distinct morphological differences, seems to be undergoing incipient speciation and contains cryptic species; the hypervariable V4 region of the small subunit rDNA and D1-D2 region of the large subunit rDNA are the promising candidates for general species delimitation in Euplotes.

Evaluation of the molecular variability and characteristics of Paramecium polycaryum and Paramecium nephridiatum, within subgenus Cypriostomum (Ciliophora, Protista)

Although some Paramecium species are suitable research objects in many areas of life sciences, the biodiversity structure of other species is almost unknown. In the current survey, we present a molecular analysis of 60 Cypriostomum strains, which for the first time allows for the study of intra- and interspecific relationships within that subgenus, as well as the assessment of the biogeography patterns of its morphospecies. Analysis of COI mtDNA variation revealed three main clades (separated from each other by approximately 130 nucleotide substitutions), each one with internal sub-clusters (differing by 30 to 70 substitutions - a similar range found between P. aurelia cryptic species and P. bursaria syngens). The first clade is represented exclusively by P. polycaryum; the second one includes only four strains identified as P. calkinsi. The third cluster seems to be paraphyletic, as it includes P. nephridiatum, P. woodruffi, and Eucandidatus P. hungarianum. Some strains, previously identified as P. calkinsi, had COI sequences identical or very similar to P. nephridiatum ones. Morphological reinvestigation of several such strains revealed common morphological features with P. nephridiatum. The paper contains new information concerning speciation within particular species, i.e. existence of cryptic species within P. polycaryum (three) and in P. nephridiatum (six).

Are We Overestimating Protistan Diversity in Nature?

Documenting the immense diversity of single-celled, eukaryotic organisms (protists) has been a formidable challenge for ecologists. These species were originally defined by morphological criteria, but shortcomings of the morphospecies concept, and a bewildering array of sizes and cellular attributes, has made constructing a taxonomy that is useful for ecologists nearly impossible. Consequently, physiological and genetic information has been integrated to address these shortcomings, and to develop the framework of a unifying taxonomy. DNA sequence information, in particular, has revolutionized studies of protistan diversity. However, the exponential increase in sequence-based protistan species richness published from field surveys in recent years raises the question of whether we have moved beyond characterizing species-level diversity and begun to reveal intraspecies diversity. The answer to that question appears to be 'yes', at least for some protistan lineages. The need to document such microdiversity may be justified, but it is important for protistologists to recognize and acknowledge that possibility, and its consequences.

Comparative genomics and the nature of placozoan species

Placozoans are a phylum of nonbilaterian marine animals currently represented by a single described species, Trichoplax adhaerens, Schulze 1883. Placozoans arguably show the simplest animal morphology, which is identical among isolates collected worldwide, despite an apparently sizeable genetic diversity within the phylum. Here, we use a comparative genomics approach for a deeper appreciation of the structure and causes of the deeply diverging lineages in the Placozoa. We generated a high-quality draft genome of the genetic lineage H13 isolated from Hong Kong and compared it to the distantly related T. adhaerens. We uncovered substantial structural differences between the two genomes that point to a deep genomic separation and provide support that adaptation by gene duplication is likely a crucial mechanism in placozoan speciation. We further provide genetic evidence for reproductively isolated species and suggest a genus-level difference of H13 to T. adhaerens, justifying the designation of H13 as a new species, Hoilungia hongkongensis nov. gen., nov. spec., now the second described placozoan species and the first in a new genus. Our multilevel comparative genomics approach is, therefore, likely to prove valuable for species distinctions in other cryptic microscopic animal groups that lack diagnostic morphological characters, such as some nematodes, copepods, rotifers, or mites.

Diversity and potential activity patterns of planktonic eukaryotic microbes in a mesoeutrophic coastal area (eastern English Channel)

The diversity of planktonic eukaryotic microbes was studied at a coastal station of the eastern English Channel (EEC) from March 2011 to July 2015 (77 samples) using high throughput sequencing (454-pyrosequencing and Illumina) of the V2-V3 hypervariable region of the 18S SSU rDNA gene. Similar estimations of OTU relative abundance and taxonomic distribution for the dominant higher taxonomic groups (contributing >1% of the total number of OTUs) were observed with the two methods (Kolmogorov-Smirnov p-value = 0.22). Eight super-groups were identified throughout all samples: Alveolata, Stramenopiles, Opisthokonta, Hacrobia, Archeaplastida, Apusozoa, Rhizaria, and Amoebozoa (ordered by decreasing OTU richness). To gain further insight into microbial activity in the EEC, ribosomal RNA was extracted for samples from 2013–2015 (30 samples). Analysis of 18S rDNA and rRNA sequences led to the detection of 696 and 700 OTUs, respectively. Cluster analysis based on OTUs’ abundance indicated three major seasonal groups that were associated to spring, winter/autumn, and summer conditions. The clusters inferred from rRNA data showed a clearer seasonal representation of the community succession than the one based on rDNA. The rRNA/rDNA ratio was used as a proxy for relative cell activity. When all OTUs were considered, the average rRNA:rDNA ratio showed a linear trend around the 1:1 line, suggesting a linear relation between OTU abundance (rDNA) and activity (rRNA). However, this ratio was highly variable over time when considering individual OTUs. Interestingly, the OTU affiliated with P. globosa displayed rRNA:rDNA ratio that allowed to delimit high vs low abundance and high vs low activity periods. It unveiled quite well the Phaeocystis bloom dynamic regarding cell proliferation and activity, and could even be used as early indicator of an upcoming bloom.

Feeding by the harmful phototrophic dinoflagellate Takayama tasmanica (Family Kareniaceae)

The trophic mode of a phototrophic dinoflagellate is a critical factor in the dynamics of its harmful algal bloom. Recent discoveries of the mixotrophic capabilities of phototrophic dinoflagellates have changed the traditional view of bloom dynamics and prediction models. Here, mixotrophy in the harmful phototrophic dinoflagellate Takayama tasmanica was examined. Moreover, growth and ingestion rates of T. tasmanica on each of Alexandrium minutum CCMP1888 and Alexandrium tamarense CCMP1493, suitable prey, were determined as a function of prey concentration. This study reported for the first time that T. tasmanica is a mixotrophic species. Among the phytoplankton species offered as prey, T. tasmanica fed on all prey species whose equivalent spherical diameter (ESD) was greater than 30 μm, but also A. minutum whose ESD was 19 μm. In contrast, T. tasmanica did not feed on the phototrophic dinoflagellates Heterocapsa triquetra, Gymnodinium aureolum, Scrippsiella acuminata (previously S. trochoidea), Cochlodinium polykrikoides, Alexandrium affine, Alexandrium insuetum, and Alexandrium pacificum that its sister species Takayama helix is able to feed on. With increasing mean prey concentration, ingestion rates of T. tasmanica on A. minutum increased, but became saturated at the prey concentrations of >2130 cells mL^-1 (1070 ng C mL^-1). The maximum ingestion rate (MIR) of T. tasmanica on A. minutum was 0.5 ng C predator^-1 d^-1 (1.0 cells predator^-1 d^-1) which is only 64% of the body carbon of a T. tasmanica cell. Growth rates of T. tasmanica on A. minutum were not affected by prey concentrations. Thus, the low maximum ingestion rate is likely to be responsible for the small increases of its growth rate through mixotrophy. In addition, neither growth nor ingestion rates of T. tasmanica feeding on Alexandrium tamarense were affected by prey concentrations. The maximum ingestion rate of T. tasmanica on A. minutum was considerably lower than that of T. helix on the same prey species. Therefore, the mixotrophic ability of T. tasmanica is weaker than that of T. helix, and also T. tasmanica may have an ecological niche different from that of T. helix in marine ecosystems.

Unravel the hidden protistan diversity: application of blocking primers to suppress PCR amplification of metazoan DNA

Planktonic protists, including both autotroph and heterotroph, have been recognized as a major contributor to primary production and consumers of bacteria, archaea and picophytoplankton. However, the understanding of protistan diversity is typically impeded by the large amount of metazoans when employing universal primers to environmental samples. In this study, we developed 3 blocking primers to inhibit the amplification of metazoan DNA in PCR. First, we optimized the design and concentration of blocking primers by using 2 metazoans and 3 protists in quantitative PCR. No significant difference (q > 0.05) was found in protistan community structure at phylum, family or OTU levels and genetic diversity between samples amplified with and without adding blocking primers, indicating that the blocking primer does not alter the composition of protistan community. An application to samples with a high abundance of metazoans showed that the blocking primer can reduce 42.1-72.4% of metazoan sequences, resulting the retrieval of a higher protistan richness and diversity. The use of the blocking primer helped to identify protistan community composition from a large size fraction, which is usually dominated by metazoans. Protistan community from the small (0.2-10 μm) and large (> 10 μm) size fractions exhibited a low similarity of 36.6% and shared 14.7% OTUs. About 63.8% of the OTUs were unique to the large fraction. Species from groups, such as Lepidodinium, Warnowia, Kareniaceae, Torodiniales, Phaeocystis, Chrysochromulina and Chrysophyceae, were enriched in the large fraction, indicating that they could be largely underestimated in studies that exclude the large-sized cells. Blocking primers are a promising tool to increase the efficiency for the characterization of protistan diversity in aquatic planktonic ecosystems.

Kin Discrimination in Protists: From Many Cells to Single Cells and Backwards

During four decades (1960-1990s), the conceptualization and experimental design of studies in kin recognition relied on work with multicellular eukaryotes, particularly Unikonta (including invertebrates and vertebrates) and some Bikonta (including plants). This pioneering research had an animal behavior approach. During the 2000s, work on taxa-, clone- and kin-discrimination and recognition in protists produced genetic and molecular evidence that unicellular organisms (e.g. Saccharomyces, Dictyostelium, Polysphondylium, Tetrahymena, Entamoeba and Plasmodium) could distinguish between same (self or clone) and different (diverse clones), as well as among conspecifics of close or distant genetic relatedness. Here, we discuss some of the research on the genetics of kin discrimination/recognition and highlight the scientific progress made by switching emphasis from investigating multicellular to unicellular systems (and backwards). We document how studies with protists are helping us to understand the microscopic, cellular origins and evolution of the mechanisms of kin discrimination/recognition and their significance for the advent of multicellularity. We emphasize that because protists are among the most ancient organisms on Earth, belong to multiple taxonomic groups and occupy all environments, they can be central to reexamining traditional hypotheses in the field of kin recognition, reformulating concepts, and generating new knowledge.

Nomenclature for the Nameless: A Proposal for an Integrative Molecular Taxonomy of Cryptic Diversity Exemplified by Planktonic Foraminifera

Investigations of biodiversity, biogeography, and ecological processes rely on the identification of "species" as biologically significant, natural units of evolution. In this context, morphotaxonomy only provides an adequate level of resolution if reproductive isolation matches morphological divergence. In many groups of organisms, morphologically defined species often disguise considerable genetic diversity, which may be indicative of the existence of cryptic species. The diversity hidden by morphological species can be disentangled through genetic surveys, which also provide access to data on the ecological distribution of genetically circumscribed units. These units can be identified by unique DNA sequence motifs and allow studies of evolutionary and ecological processes at different levels of divergence. However, the nomenclature of genetically circumscribed units within morphological species is not regulated and lacks stability. This represents a major obstacle to efforts to synthesize and communicate data on genetic diversity for multiple stakeholders. We have been confronted with such an obstacle in our work on planktonic foraminifera, where the stakeholder community is particularly diverse, involving geochemists, paleoceanographers, paleontologists, and biologists, and the lack of stable nomenclature beyond the level of formal morphospecies prevents effective transfer of knowledge. To circumvent this problem, we have designed a stable, reproducible, and flexible nomenclature system for genetically circumscribed units, analogous to the principles of a formal nomenclature system. Our system is based on the definition of unique DNA sequence motifs collocated within an individual, their typification (in analogy with holotypes), utilization of their hierarchical phylogenetic structure to define levels of divergence below that of the morphospecies, and a set of nomenclature rules assuring stability. The resulting molecular operational taxonomic units remain outside the domain of current nomenclature codes, but are linked to formal morphospecies as regulated by the codes. Subsequently, we show how this system can be applied to classify genetically defined units using the SSU rDNA marker in planktonic foraminifera and we highlight its potential use for other groups of organisms where similarly high levels of connectivity between molecular and formal taxonomies can be achieved.

Amplicon-Based Pyrosequencing Reveals High Diversity of Protistan Parasites in Ships' Ballast Water: Implications for Biogeography and Infectious Diseases

Ships' ballast water (BW) commonly moves macroorganisms and microorganisms across the world's oceans and along coasts; however, the majority of these microbial transfers have gone undetected. We applied high-throughput sequencing methods to identify microbial eukaryotes, specifically emphasizing the protistan parasites, in ships' BW collected from vessels calling to the Chesapeake Bay (Virginia and Maryland, USA) from European and Eastern Canadian ports. We utilized tagged-amplicon 454 pyrosequencing with two general primer sets, amplifying either the V4 or V9 domain of the small subunit (SSU) of the ribosomal RNA (rRNA) gene complex, from total DNA extracted from water samples collected from the ballast tanks of bulk cargo vessels. We detected a diverse group of protistan taxa, with some known to contain important parasites in marine systems, including Apicomplexa (unidentified apicomplexans, unidentified gregarines, Cryptosporidium spp.), Dinophyta (Blastodinium spp., Euduboscquella sp., unidentified syndinids, Karlodinium spp., Syndinium spp.), Perkinsea (Parvilucifera sp.), Opisthokonta (Ichthyosporea sp., Pseudoperkinsidae, unidentified ichthyosporeans), and Stramenopiles (Labyrinthulomycetes). Further characterization of groups with parasitic taxa, consisting of phylogenetic analyses for four taxa (Cryptosporidium spp., Parvilucifera spp., Labyrinthulomycetes, and Ichthyosporea), revealed that sequences were obtained from both known and novel lineages. This study demonstrates that high-throughput sequencing is a viable and sensitive method for detecting parasitic protists when present and transported in the ballast water of ships. These data also underscore the potential importance of human-aided dispersal in the biogeography of these microbes and emerging diseases in the world's oceans.

Molecular Identification of Soil Eukaryotes and Focused Approaches Targeting Protist and Faunal Groups Using High-Throughput Metabarcoding

While until recently the application of high-throughput sequencing approaches has mostly been restricted to bacteria and fungi, these methods have now also become available to less often studied (eukaryotic) groups, such as fauna and protists. Such approaches allow routine diversity screening for large numbers of samples via DNA metabarcoding. Given the enormous taxonomic diversity within the eukaryote tree of life, metabarcoding approaches targeting a single specific DNA region do not allow to discriminate members of all eukaryote clades at high taxonomic resolution. Here, we report on protocols that enable studying the diversity of soil eukaryotes and, at high taxonomic resolution, of individual faunal and protist groups therein using a tiered approach: first, the use of a general eukaryotic primer set targeting a wide range of eukaryotes provides a rough impression on the entire diversity of protists and faunal groups. Second, more focused approaches enable deciphering subsets of soil eukaryotes in higher taxonomic detail. We provide primers and protocols for two examples: soil microarthropods and cercozoan protists.

Delimiting Species Boundaries within a Paraphyletic Species Complex: Insights from Morphological, Genetic, and Molecular Data on Paramecium sonneborni (Paramecium aurelia species complex, Ciliophora, Protozoa)

The demarcation of boundaries between protist species is often problematic because of the absence of a uniform species definition, the abundance of cryptic diversity, and the occurrence of convergent morphology. The ciliates belonging to the Paramecium aurelia complex, consisting of 15 species, are a good model for such systematic and evolutionary studies. One member of the complex is P. sonneborni, previously known only from one stand in Texas (USA), but recently found in two new sampling sites in Cyprus (creeks running to Salt Lake and Oroklini Lake near Larnaca). The studied Paramecium sonneborni strains (from the USA and Cyprus) reveal low viability in the F1 and F2 generations of interstrain hybrids and may be an example of ongoing allopatric speciation. Despite its molecular distinctiveness, we postulate that P. sonneborni should remain in the P. aurelia complex, making it a paraphyletic taxon. Morphological studies have revealed that some features of the nuclear apparatus of P. sonneborni correspond to the P. aurelia spp. complex, while others are similar to P. jenningsi and P. schewiakoffi. The observed discordance indicates rapid splitting of the P. aurelia-P. jenningsi-P. schewiakoffi group, in which genetic, morphological, and molecular boundaries between species are not congruent.

Estimating Protistan Diversity Using High-Throughput Sequencing

Sequencing hypervariable regions from the 18S rRNA gene is commonly employed to characterize protistan biodiversity, yet there are concerns that short reads do not provide the same taxonomic resolution as full-length sequences. A total of 7,432 full-length sequences were used to perform an in silico analysis of how sequences of various lengths and target regions impact downstream ecological interpretations. Sequences that were longer than 400 nucleotides and included the V4 hypervariable region generated results similar to those derived from full-length 18S rRNA gene sequences. Present high-throughput sequencing capabilities are approaching protistan diversity estimation comparable to whole gene sequences.

Evidence of Taxa-, Clone-, and Kin-discrimination in Protists: Ecological and Evolutionary Implications

Unicellular eukaryotes, or protists, are among the most ancient organisms on Earth. Protists belong to multiple taxonomic groups; they are widely distributed geographically and in all environments. Their ability to discriminate among con- and heterospecifics has been documented during the past decade. Here we discuss exemplar cases of taxa-, clone-, and possible kin-discrimination in five major lineages: Mycetozoa (Dictyostelium, Polysphondylium), Dikarya (Saccharomyces), Ciliophora (Tetrahymena), Apicomplexa (Plasmodium) and Archamoebae (Entamoeba). We summarize the proposed genetic mechanisms involved in discrimination-mediated aggregation (self versus different), including the csA, FLO and trg (formerly lag) genes, and the Proliferation Activation Factors (PAFs), which facilitate clustering in some protistan taxa. We caution about the experimental challenges intrinsic to studying recognition in protists, and highlight the opportunities for exploring the ecology and evolution of complex forms of cell-cell communication, including social behavior, in a polyphyletic, still superficially understood group of organisms. Because unicellular eukaryotes are the evolutionary precursors of multicellular life, we infer that their mechanisms of taxa-, clone-, and possible kin-discrimination gave origin to the complex diversification and sophistication of traits associated with species and kin recognition in plants, fungi, invertebrates and vertebrates.

Investigating microbial eukaryotic diversity from a global census: insights from a comparison of pyrotag and full-length sequences of 18S rRNA genes

Next-generation DNA sequencing (NGS) approaches are rapidly surpassing Sanger sequencing for characterizing the diversity of natural microbial communities. Despite this rapid transition, few comparisons exist between Sanger sequences and the generally much shorter reads of NGS. Operational taxonomic units (OTUs) derived from full-length (Sanger sequencing) and pyrotag (454 sequencing of the V9 hypervariable region) sequences of 18S rRNA genes from 10 global samples were analyzed in order to compare the resulting protistan community structures and species richness. Pyrotag OTUs called at 98% sequence similarity yielded numbers of OTUs that were similar overall to those for full-length sequences when the latter were called at 97% similarity. Singleton OTUs strongly influenced estimates of species richness but not the higher-level taxonomic composition of the community. The pyrotag and full-length sequence data sets had slightly different taxonomic compositions of rhizarians, stramenopiles, cryptophytes, and haptophytes, but the two data sets had similarly high compositions of alveolates. Pyrotag-based OTUs were often derived from sequences that mapped to multiple full-length OTUs at 100% similarity. Thus, pyrotags sequenced from a single hypervariable region might not be appropriate for establishing protistan species-level OTUs. However, nonmetric multidimensional scaling plots constructed with the two data sets yielded similar clusters, indicating that beta diversity analysis results were similar for the Sanger and NGS sequences. Short pyrotag sequences can provide holistic assessments of protistan communities, although care must be taken in interpreting the results. The longer reads (>500 bp) that are now becoming available through NGS should provide powerful tools for assessing the diversity of microbial eukaryotic assemblages.

Cryptic diversity in a cosmopolitan diatom known as Asterionellopsis glacialis (Fragilariaceae): Implications for ecology, biogeography, and taxonomy

PREMISE OF THE STUDY

Diatoms have long been known as the most species-rich of algal groups, with a wide range of estimates for species number (20-200 k) due to differing species concepts. The fine valve structure in Asterionellopsis glacialis, a diatom believed cosmopolitan and eurytopic, has never been systematically examined using modern microscopy and is an excellent candidate to genetically test morphology-based conspecificity among its geographically distant culture isolates.

METHODS

Isolates from the Atlantic and Pacific Oceans that were morphologically delineated as A. glacialis were genetically characterized using three nuclear DNA regions (two 18S rDNA fragments and ITS region) and one plastidal (rbcL) and one mitochondrial gene (cox1) and related to SEM-based morphometrics.

KEY RESULTS

Five genetically distinct groupings were found, four of which are new to science. ITS2 RNA transcript secondary structure was species specific as were plastidal and mitochondrial genes, while the 18S gene fragments did not diverge sufficiently to segregate new species efficiently. We genetically circumscribed the A. glacialis epitype.

CONCLUSIONS

The morphological diversification of the species examined in this study lags behind their genetic divergence. The currently accepted 2% cutoff level of operational taxonomic unit (OTU) clustering in 18S rDNA environmental sequencing is too high to recognize genetic diversity in Asterionellopsis and very likely in many other species. Our results support the notion that a considerable number of species and diversity remain to be discovered among diatoms and that species number may be more in line with higher estimates. Molecular signatures of the species discovered here will aid in their globally consistent identification and ultimate understanding of their ecology.