Diversity estimates of microeukaryotes below the chemocline of the anoxic Mariager Fjord, Denmark

Microbial communities of extreme environments have often been assumed to have low species richness. We analysed 18S rRNA gene signatures in a sample collected below the chemocline of the anoxic Mariager Fjord in Denmark, and from these data we computed novel parametric and standard nonparametric estimates of protistan phylotype richness. Our results indicate unexpectedly high richness in this environment: at the 99.5% phylotype definition, our most conservative estimate was 568 phylotypes (+/-114, standard error). Phylogenetic analyses revealed that the sequences collected cover the majority of described lineages in the eukaryotic domain. Out of 384 sequences analysed, 307 were identified as protistan targets, none of which was identical to known sequences. However, based on what is known about species that are phylogenetically related to the Mariager sequences, most of the latter seem to belong to strictly or facultative anaerobe organisms. We also found signatures that together with other environmental 18S rRNA gene sequences represent environmental clades of possibly high taxonomic levels (class to kingdom level). One of these clades, consisting exclusively of sequences from anoxic sampling sites, branches at the base of the eukaryotic evolutionary tree among the earliest eukaryotic lineages. Assuming eukaryotic evolution under oxygen-depleted conditions, these sequences may represent immediate descendants of early eukaryotic ancestors.     

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.     

Molecular Identification of Nanoplanktonic Protists Based on Small Subunit Ribosomal RNA Gene Sequences for Ecological Studies1

Nanoplanktonic protists are comprised of a diverse assemblage of species which are responsible for a variety of trophic processes in marine and freshwater ecosystems. Current methods for identifying small protists by electron microscopy do not readily permit both identification and enumeration of nanoplanktonic protists in field samples. Thus, one major goal in the application of molecular approaches in protistan ecology has been the detection and quantification of individual species in natural water samples. Sequences of small subunit ribosomal RNA (SSU rRNA) genes have proven to be useful towards achieving this goal. Comparison of sequences from clone libraries of protistan SSU rRNA genes amplified from natural assemblages of protists by the polymerase chain reaction (PCR) can be used to examine protistan diversity. Furthermore, oligonucleotide probes complementary to short sequence regions unique to species of small protists can be designed by comparative analysis of rRNA gene sequences. These probes may be used to either detect the RNA of particular species of protists in total nucleic acid extracts immobilized on membranes, or the presence of target species in water samples via in situ hybridization of whole cells. Oligonucleotide probes may also serve as primers for the selective amplification of target sequences from total population DNA by PCR. Thus, molecular sequence information is becoming increasingly useful for identifying and enumerating protists, and for studying their spatial and temporal distribution in nature. Knowledge of protistan species composition, abundance and variability in an environment can ultimately be used to relate community structure to various aspects of community function and biogeochemical activity.     

Microeukaryote Community Patterns along an O2/H2S Gradient in a Supersulfidic Anoxic Fjord (Framvaren, Norway)

To resolve the fine-scale architecture of anoxic protistan communities, we conducted a cultivation-independent 18S rRNA survey in the superanoxic Framvaren Fjord in Norway. We generated three clone libraries along the steep O₂/H₂S gradient, using the multiple-primer approach. Of 1,100 clones analyzed, 753 proved to be high-quality protistan target sequences. These sequences were grouped into 92 phylotypes, which displayed high protistan diversity in the fjord (17 major eukaryotic phyla). Only a few were closely related to known taxa. Several sequences were dissimilar to all previously described sequences and occupied a basal position in the inferred phylogenies, suggesting that the sequences recovered were derived from novel, deeply divergent eukaryotes. We detected sequence clades with evolutionary importance (for example, clades in the euglenozoa) and clades that seem to be specifically adapted to anoxic environments, challenging the hypothesis that the global dispersal of protists is uniform. Moreover, with the detection of clones affiliated with jakobid flagellates, we present evidence that primitive descendants of early eukaryotes are present in this anoxic environment. To estimate sample coverage and phylotype richness, we used parametric and nonparametric statistical methods. The results show that although our data set is one of the largest published inventories, our sample missed a substantial proportion of the protistan diversity. Nevertheless, statistical and phylogenetic analyses of the three libraries revealed the fine-scale architecture of anoxic protistan communities, which may exhibit adaptation to different environmental conditions along the O₂/H₂S gradient.     

Microeukaryote community patterns along an O2/H2S gradient in a supersulfidic anoxic fjord (Framvaren, Norway)

To resolve the fine-scale architecture of anoxic protistan communities, we conducted a cultivation-independent 18S rRNA survey in the superanoxic Framvaren Fjord in Norway. We generated three clone libraries along the steep O(2)/H(2)S gradient, using the multiple-primer approach. Of 1,100 clones analyzed, 753 proved to be high-quality protistan target sequences. These sequences were grouped into 92 phylotypes, which displayed high protistan diversity in the fjord (17 major eukaryotic phyla). Only a few were closely related to known taxa. Several sequences were dissimilar to all previously described sequences and occupied a basal position in the inferred phylogenies, suggesting that the sequences recovered were derived from novel, deeply divergent eukaryotes. We detected sequence clades with evolutionary importance (for example, clades in the euglenozoa) and clades that seem to be specifically adapted to anoxic environments, challenging the hypothesis that the global dispersal of protists is uniform. Moreover, with the detection of clones affiliated with jakobid flagellates, we present evidence that primitive descendants of early eukaryotes are present in this anoxic environment. To estimate sample coverage and phylotype richness, we used parametric and nonparametric statistical methods. The results show that although our data set is one of the largest published inventories, our sample missed a substantial proportion of the protistan diversity. Nevertheless, statistical and phylogenetic analyses of the three libraries revealed the fine-scale architecture of anoxic protistan communities, which may exhibit adaptation to different environmental conditions along the O(2)/H(2)S gradient.     

Broad range DNA probes for detecting and amplifying eubacterial nucleic acids

In this report we describe and characterize two oligomer probes that are broadly homologous to conserved eubacterial 16S ribosomal RNA (rRNA) sequences not present in human 18 rRNA or human mitochondrial 12S rRNA. One or both of the probes can detect all of 23 phylogenetically diverse eubacterial nucleic acids against which they were tested by dot blot hybridization. A sensitivity of about 1 bacterium per 10 eukaryotic cells was achieved. By using these oligomer sequences or their complements as primers in the polymerase chain reaction (PCR), the equivalent of 1 pg of E. coli DNA was detected in the presence of a large excess of eukaryotic DNA. Information useful for partial phylogenetic classification of detected organisms may be obtained by direct sequence analysis of the amplified DNA and comparison with known sequences or catalogs. Such broadly homologous probes offer advantages over more narrowly specific probes for detecting organisms whose identity is unknown. They could thus be employed for recognizing infection by organisms that cannot be cultured as may occur, for example, in tissue culture or in plant or animal diseases of unknown cause, provided the probes fail to hybridize with host nucleic acids.     

Culturing Heterotrophic Protists from the Baltic Sea: Mostly the “Usual Suspects” but a Few Novelties as Well

The study of cultured strains has a long tradition in protistological research and has greatly contributed to establishing the morphology, taxonomy, and ecology of many protist species. However, cultivation‐independent techniques, based on 18S rRNA gene sequences, have demonstrated that natural protistan assemblages mainly consist of hitherto uncultured protist lineages. This mismatch impedes the linkage of environmental diversity data with the biological features of cultured strains. Thus, novel taxa need to be obtained in culture to close this knowledge gap. In this study, traditional cultivation techniques were applied to samples from coastal surface waters and from deep oxygen‐depleted waters of the Baltic Sea. Based on 18S rRNA gene sequencing, 126 monoclonal cultures of heterotrophic protists were identified. The majority of the isolated strains were affiliated with already cultured and described taxa, mainly chrysophytes and bodonids. This was likely due to “culturing bias” but also to the eutrophic nature of the Baltic Sea. Nonetheless, ~ 12% of the isolates in our culture collection showed highly divergent 18S rRNA gene sequences compared to those of known organisms and thus may represent novel taxa, either at the species level or at the genus level. Moreover, we also obtained evidence that some of the isolated taxa are ecologically relevant, under certain conditions, in the Baltic Sea.     

A multiple PCR-primer approach to access the microeukaryotic diversity in environmental samples

The Cariaco Basin off the Venezuelan coast in the Caribbean Sea is the world's largest truly marine body of anoxic water. The first rRNA survey of microbial eukaryotes in this environment revealed a number of novel lineages, but sampled only a fraction of the entire diversity. The goal of this study was to significantly improve recovery of protistan rRNA from the Basin. This was achieved by a systematic application of multiple PCR primer sets and substantially larger sequencing efforts. We focused on the most diverse habitat in the basin, anoxic waters approximately 100m below the oxic-anoxic interface, and detected novel lineages that escaped the single PCR primer approach. All clones obtained proved unique. A 99% sequence similarity cut-off value combined these clones into operational taxonomic units (OTUs), over 75% of which proved novel. Some of these OTUs form deep branches within established protistan groups. Others signify discovery of novel protistan lineages that appear unrelated to any known microeukaryote. Surprisingly, even this large-scale multi-primer rRNA approach still missed a substantial part of the samples' rRNA diversity. The overlap between the species lists obtained with different primers is low, with only 4% of OTUs shared by all three libraries, and the number of species detected only once is large (55%). This strongly indicates that, at least in anoxic environments, protistan diversity may be much larger than is commonly thought. A single sample appears to contain thousands of largely novel protistan species. Multiple PCR primer combinations may be needed to capture these species.     

Haptophyte Diversity and Vertical Distribution Explored by 18S and 28S Ribosomal RNA Gene Metabarcoding and Scanning Electron Microscopy

Haptophyta encompasses more than 300 species of mostly marine pico‐ and nanoplanktonic flagellates. Our aims were to investigate the Oslofjorden haptophyte diversity and vertical distribution by metabarcoding, and to improve the approach to study haptophyte community composition, richness and proportional abundance by comparing two rRNA markers and scanning electron microscopy (SEM). Samples were collected in August 2013 at the Outer Oslofjorden, Norway. Total RNA/cDNA was amplified by haptophyte‐specific primers targeting the V4 region of the 18S, and the D1‐D2 region of the 28S rRNA. Taxonomy was assigned using curated haptophyte reference databases and phylogenetic analyses. Both marker genes showed Chrysochromulinaceae and Prymnesiaceae to be the families with highest number of Operational Taxonomic Units (OTUs), as well as proportional abundance. The 18S rRNA data set also contained OTUs assigned to eight supported and defined clades consisting of environmental sequences only, possibly representing novel lineages from family to class. We also recorded new species for the area. Comparing coccolithophores by SEM with metabarcoding shows a good correspondence with the 18S rRNA gene proportional abundances. Our results contribute to link morphological and molecular data and 28S to 18S rRNA gene sequences of haptophytes without cultured representatives, and to improve metabarcoding methodology.     

Multiple ribosomal RNA cleavage pathways in mammalian cells

The sequence content of mouse L cell pre-rRNA was examined by RNA gel transfer and blot hybridization. Nuclear RNAs were separated by agarose gel electrophoresis, transferred to diazo-paper, and hybridized to twelve different restriction fragments that are complementary to various sections of 45S pre-rRNA. An abundant new 34S pre-rRNA and less abundant new 37S, 26S and 17S pre-rRNAs were detected. The presence of these new pre-rRNAs suggests the existence of at least two new pre-rRNA cleavage pathways. 34S and 26S pre-rRNAs were also detected in HeLa cells suggesting that these new cleavage pathways are characteristic of mammalian cells. Further, an abundant new 12S precursor to 5.8S rRNA was also detected and is common to all the proposed cleavage pathways. The previously identified 45S, 41S, 32S and 20S pre-rRNAs were readily detected and their general structure confirmed. The 20S pre-rRNA is characteristic of the known pathway used by HeLa and other cells, and its presence suggests that growing mouse L cells use this pre-rRNA cleavage pathway. The 36S pre-rRNA characteristic of the previously described mouse L cell cleavage pathway was not detected. In all these cleavage pathways pre-rRNA cleavage sites are apparently identical and occur at or near the termini of the mature 18S, 5.8S and 28S rRNA sequences. The pathways differ only in the temporal order of cleavage at these sites.     

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.     

Protistan microbial observatory in the Cariaco Basin,Caribbean. I. Pyrosequencing vs Sanger insights into species richness

Microbial diversity and distribution are topics of intensive research. In two companion papers in this issue, we describe the results of the Cariaco Microbial Observatory (Caribbean Sea, Venezuela). The Basin contains the largest body of marine anoxic water, and presents an opportunity to study protistan communities across biogeochemical gradients. In the first paper, we survey 18S ribosomal RNA (rRNA) gene sequence diversity using both Sanger- and pyrosequencing-based approaches, employing multiple PCR primers, and state-of-the-art statistical analyses to estimate microbial richness missed by the survey. Sampling the Basin at three stations, in two seasons, and at four depths with distinct biogeochemical regimes, we obtained the largest, and arguably the least biased collection of over 6000 nearly full-length protistan rRNA gene sequences from a given oceanographic regime to date, and over 80 000 pyrosequencing tags. These represent all major and many minor protistan taxa, at frequencies globally similar between the two sequence collections. This large data set provided, via the recently developed parametric modeling, the first statistically sound prediction of the total size of protistan richness in a large and varied environment, such as the Cariaco Basin: over 36 000 species, defined as almost full-length 18S rRNA gene sequence clusters sharing over 99% sequence homology. This richness is a small fraction of the grand total of known protists (over 100 000–500 000 species), suggesting a degree of protistan endemism.     

Microeukaryotic diversity in marine environments,an analysis of surface layer sediments from the East Sea

Molecular techniques, based on clone library of 18S rRNA gene, were employed to ascertain the diversity of microeukaryotic organisms in sediments from the East Sea. A total of 261 clones were recovered from surface sediments. Most of the clone sequences (90%) were affiliated with protists, dominated by Ciliates (18%) and Dinoflagellates (19%) of Alveolates, phototrophic Stramenopiles (11%), and Cercozoa (20%). Many of the clones were related to uncultivated eukaryotes clones retrieved from anoxic environments with several highly divergent 18S rRNA gene sequences. However, no clones were related to cultivated obligate anaerobic protists. Protistan communities between subsurface layers of 1 and 9 cm shared 23% of total phylotypes which comprised 64% of total clones retrieved. Analysis of diversity indices and rarefaction curve showed that the protistan community within the 1 cm layer exhibited higher diversity than the 9 cm layer. Our results imply that diverse protists remain to be uncovered within marine benthic environments.     

Unveiling trophic functions of uncultured protist taxa by incubation experiments in the brackish baltic sea

Background

Our knowledge of the phylogeny and diversity of aquatic protists is rapidly increasing due to molecular surveys and next-generation sequencing approaches. This has led to a considerable discrepancy between the taxa known from cultures and those known from environmental 18S rRNA gene sequences. Hence, it is generally difficult to assign ecological functions to new taxa detected by culture-independent molecular approaches.

Methodology/Principal Findings

A combination of unamended dark incubations and 18S rRNA sequencing was chosen to link molecular diversity data of uncultured protists with heterotrophic, presumably bacterivorous, growth. The incubations, conducted with Baltic Sea brackish water, resulted in a consistent shift from a protistan community dominated by phototrophs to one in which heterotrophs predominated. This was determined on the basis of cell abundance and 18S rRNA sequences derived from fingerprint analysis and clone libraries. The bulk of enriched phylotypes after incubation were related to hitherto uncultured marine taxa within chrysophytes, ochrophytes, choanoflagellates, cercozoans, and picobiliphytes, mostly represented in recently established or here defined environmental clades. Their growth in the dark, together with coinciding results from studies with a similar objective, provides evidence that these uncultured taxa represent heterotrophic or mixotrophic species.

Conclusions/Significance

These findings shed some light into the trophic role of diverse uncultured protists especially within functionally heterogeneous groups (e.g., chrysophytes, ochrophytes) and groups that appear to be puzzling with regard to their nutrition (picobiliphytes). Additionally, our results indicate that the heterotrophic flagellate community in the southwestern Baltic Sea is dominated by species of marine origin. The combination of unamended incubations with molecular diversity analysis is thus confirmed as a promising approach to explore the trophic mode of environmentally relevant protist taxa for which only sequence data are currently available.     

Distinct protistan assemblages characterize the euphotic zone and deep sea (2500 m) of the western North Atlantic (Sargasso Sea and Gulf Stream)

Protistan diversity was characterized at three locations in the western North Atlantic (Sargasso Sea and Gulf Stream) by sequencing 18S rRNA genes in samples from euphotic (< or = 125 m) and bathypelagic depths (2500 m). A total of 923 partial-length protistan sequences were analysed, revealing 324 distinct operational taxonomic units (OTUs) determined by an automated OTU-calling program set to 95% sequence similarity. Most OTUs were comprised of only one or two sequences suggesting a large but rare pool of protistan diversity. Many OTUs from both depth strata were associated with recently described novel alveolate and stramenopile lineages while many OTUs from the bathypelagic were affiliated with Acantharea, Polycystinea and Euglenozoa and were not observed in euphotic zone libraries. Protistan assemblages from the euphotic zone and the deep sea were largely composed of distinct OTUs; only 28 of the 324 protistan OTUs were detected in both shallow and deep sea clone libraries. The diversity of protistan assemblages in the deep sea was distinctly lower than the diversity of euphotic zone assemblages. Protistan assemblages from the Gulf Stream were the most diverse for either depth strata. Overall, protistan assemblages from different stations but comparable depths were more similar than the assemblages from different depths at the same station. These data suggest that particular groups of protistan OTUs formed distinct 'shallow' and 'deep-sea' assemblages across widely spaced oceanic locales.     

Analysis of invariant sequences in 266 complete genomes

To date, the complete genome sequences of more than 250 organisms have been determined. This information can now be used to determine whether there exist any invariant sequences that are conserved among all organisms, from bacteria to plants, animals, and humans. The existence of invariant sequences would strongly suggest that these sequences have been inherited unchanged from the last common ancestor of all life, and that they have essential functions. We have developed a new software program to identify invariant sequences conserved among the currently sequenced genomes and applied this analysis to the complete genome sequences of 266 organisms. We have identified 3 invariant DNA sequences longer than or equal to 11 bp and 6 invariant amino acid sequences longer than or equal to 6 aa. The longest invariant DNA sequence, AAGTCGTACAAGGT (15 bp), was found in the 16S/18S rRNA gene. Two 8 aa sequences, GHVDHGKT in IF2 and EF-Tu and DTPGHVDF in EF-G, were the longest invariant amino acid sequences detected. These sequences could be essential elements from the genome of the last common ancestor and may have remained unchanged throughout evolution.     

Barcoding marine nematodes: an improved set of nematode 18S rRNA primers to overcome eukaryotic co-interference

Nematodes form an important component of many benthic marine ecosystems and DNA barcoding approaches could provide an insight into nematode community composition from different environments globally. We have amplified nematode 18S rRNA sequences using standard nematode18S rRNA primers from environmental DNA extracted from intertidal sediment collected from New Jersey coast, USA to test whether the published marine nematode 18S rRNA sequences from GenBank and EMBL databases can effectively assign unknown nematode sequences into genus or species level. Most of the sequenced clones showed some degree of identities with published marine nematode 18S rRNA sequences. However, relatively very few of the sequences could be assigned even to genus level based on sequence assignment rule. In addition, other eukaryotic 18S rRNA sequences were found to be co-amplified with commonly used nematode 18S rRNA primers. We found that the majority of the current nematode 18S rRNA primers will co-amplify other eukaryotes if environmental DNA is the target template. We therefore designed a new set of nematode 18S rRNA primers and evaluated them using environmental DNA in intertidal sediment from the New Jersey coast. In total, 40 clones were screened and subsequently sequenced and all the sequences showed varying degree of identities with published nematode 18S rRNA sequences from GenBank and EMBL databases, and no obvious eukaryotic co-amplicons were detected with new primers. Only 13 out of 40 clones amplified with the new primer set showed 100% identity to published Daptonema and Metachromadora 18S rRNA sequences. The current molecular databases for nematodes are dominated by sequences from NW Europe and need to be more extensively populated with new full length 18S rRNA nematode sequences collected from different biogeographic locations. The new primers developed in this study, in combination with an updated nematode 18S rRNA sequence database, would help us to better investigate and understand the diversity and community composition of free-living marine nematodes based on DNA barcoding approaches during biodiversity or biomonitoring surveys on a global-scale.     

An optimised PCR/T-RFLP fingerprinting approach for the investigation of protistan communities in groundwater environments

Due to the scarcity or complete absence of higher organisms, protists may represent an important higher trophic level (above Prokaryotes) in the food webs of groundwater habitats. Nevertheless, the importance of aquifer protists, especially in contaminated groundwater environments, is poorly understood. Partly, this may be due to a lack of adequate PCR and fingerprinting approaches for protists in aquifers, which can be considered low in protistan or high in non-target rRNA gene copy numbers. Therefore, we have validated the suitability of distinct eukaryote-targeted primer pairs and restriction endonucleases for T-RFLP fingerprinting of protistan communities. By in silico predictions, and by fingerprinting, cloning and sequencing of microeukaryote amplicons from hydrocarbon-contaminated aquifer sediment DNA, we show that the Euk20f/Euk516r primer set in combination with Bsh1236I digestion is best suited for the recovery of diverse protistan 18S rRNA lineages. In contrast to other tested primer sets, a preferred recovery of fungal and archaeal non-target amplicons was not observed. In summary, we present an optimised microeukaryote-targeted PCR/T-RFLP fingerprinting approach which may be of value for the characterisation of protistan communities in groundwater and other habitats.