Abstract
Ciliates represent central nodes in freshwater planktonic food webs, and many species show pronounced seasonality, with short-lived maxima of a few dominant taxa while many being rare or ephemeral. These observations are primarily based on morphospecies counting methods, which, however, have limitations concerning the amount and volume of samples that can be processed. For high sampling frequencies at large scales, high throughput sequencing (HTS) of freshwater ciliates seems to be a promising tool. However, several studies reported large discrepancy between species abundance determinations by molecular compared to morphological means. Therefore, we compared ciliate DNA metabarcodes (V9 regions of the 18S rRNA gene) with morphospecies counts for a 3-year study (Lake Zurich, Switzerland; biweekly sampling, n = 74). In addition, we isolated, cultivated and sequenced the 18S rRNA gene of twelve selected ciliate species that served as seeds for HTS analyses. This workflow allowed for a detailed comparison of V9 data with microscopic analyses by quantitative protargol staining (QPS). The dynamics of V9 read abundances over the seasonal cycle corroborated well with morphospecies population patterns. Annual successions of rare and ephemeral species were more adequately characterized by V9 reads than by QPS. However, numbers of species specific sequence reads only partly reflected rank orders seen by counts. In contrast, biomass-based assemblage compositions showed higher similarity to V9 read numbers, probably indicating a relation between cell sizes and numbers / sizes of macronuclei (or 18S rRNA operons). Full-length 18S rRNA sequences of ciliates assigned to certain morphospecies are urgently needed for barcoding approaches as planktonic taxa are still poorly represented in public databases and the interpretation of HTS data depends on profound reference sequences. Through linking operational taxonomic units (OTUs) with known morphospecies, we can use the deep knowledge about the autecology of these species.
Introduction
In temperate lakes, seasonal successions of phyto- and metazooplankton often show predictive dynamics, which were described in a conceptual framework, the Plankton Ecology Group (PEG) model (). During the last 30 years, ecological research has widened the understanding of mechanisms responsible for seasonality with a key finding that dynamics of heterotrophic protists should be regarded separately from metazoan zooplankton patterns (). Microzooplankton may account for a major proportion of total grazing on phytoplankton, consuming more than 50% of primary production during distinct periods in the year (; ). The algivorous protist fraction is often dominated by ciliates (Ciliophora), representing, e.g., the first and most effective grazers of phytoplankton spring blooms (; ; ; ). However, feeding modes of ciliate species are quite diverse (bacterivores, algivorous, omnivorous, predacious), thus, they form a major fraction of protistan plankton throughout the year (; ). Ciliates may even be the dominant bacterivores (; ) or consumers of autotrophic picoplankton () in certain lakes. Seasonal successions of ciliate assemblages and therein involved species have been well described for various temperate lakes (; ; ; ; ; see also the extensive literature summary in ). These studies have one common thread – their outcome is based on the identification and quantification of morphologically well-defined species (morphospecies, ) via various microscopy techniques. Morphospecies counting allows for the quantification of abundances and biomasses per water volume, which is a prerequisite to study energy flows between trophic entities in lakes.
Since 20 years, the traditional morphology based approaches were either supplemented or even replaced by molecular methods in ecological surveys (; ; ; ). Nowadays, high throughput sequencing (HTS) technologies have revolutionized the study of protistan communities in freshwater ecosystems (; ; ; ). Among others, environmental HTS offers the potential to discover cryptic and rare species which are most likely not recognized by microscopy (; ). Nevertheless, molecular methods for analyzing pro- and eukaryotic microbes have to cope with multiple sources of bias. Above all, PCR is a major source of errors due to primer bias and preferential amplification (; ; ; ). One major issue emerging from PCR and primer bias is the critical translation of amplicon abundance into organismic abundance (; ; ), but for solutions see also . Further sources of error include sequencing (), amplicon clustering (; ), and incompleteness and errors in reference databases (). Such errors can have extensive consequences for ciliate diversity estimations, often leading to an overestimation of operational taxonomic units (OTUs) richness in an environment. Beyond, many OTUs still remain unassigned and numerous ciliate morphospecies are not yet represented in reference databases (e.g., NCBI). Moreover, the high copy number of 18S rRNA genes in ciliates () and the fact that copy numbers vary greatly among species, could make quantitative results even more unreliable (). Thus, it is questionable if HTS approaches alone can be used to answer ecological questions about ciliate species such as succession patterns and their relevance for ecosystem functioning. Several surveys used a combination of molecular and morphological methods to study the diversity and dynamics of planktonic ciliate assemblages, both from marine (, ; , ) as well as from freshwater systems (; ; ). The simultaneous identification of morphospecies besides environmental sequencing helps to differentiate meaningful OTUs from potential sequencing errors. Finally, isolation and subsequent sequencing of taxonomic marker genes from identified ciliate species will enable us to link sequencing data (OTUs) to organisms and their known autecology.
In this study, we conducted a 3-year sampling campaign in Lake Zurich (biweekly sampling, n = 74) to establish a substantial data set for the evaluation and improvement of ciliate specific HTS data. Both methods, the classical morphospecies counting via quantitative silver impregnation and 18S rRNA amplicon sequencing, were used to investigate the variability in the ciliate assemblage over time. To test if the two different approaches corroborated with each other, we focused on seasonal patterns of twelve specific ciliate morphospecies. These species were isolated directly from lake water and most of them were successfully cultivated. Subsequently, their V9 regions of the 18S rRNA gene were sequenced and used as seeds for HTS analyses. This workflow resulted in new sequence information for several well-known ciliate morphospecies living in freshwater systems. We intended to address the following aims: (a) Are seasonal patterns of abundant as well as rare and ephemeral ciliate species determined via HTS correlated with successions observed by morphospecies counting? (b) Does HTS allow for a more precise description of the seasonality of rare and ephemeral ciliate species? (c) Is there a relation between numbers of sequence reads and cellular abundance or biomass values for the ciliate species? (d) Do relative numbers of sequence reads mirror the quantitative assemblage structure determined via morphospecies counting?
Materials and Methods
Investigation Site, Sampling and Environmental Parameters
Lake Zurich (47°19.3′N, 8°33.9′E; Switzerland) is a pre-alpine, oligo-mesotrophic, and monomictic lake which is the object of intensive limnological research since five decades (; ). The major morphometric descriptors are: altitude = 406 m a.s.l., maximum depth = 136 m, mean depth = 51 m, area = 65.06 km2, volume = 3.3 km3, and theoretical water renewal time = 1.2 years. The lake serves as drinking water reservoir for more than 1 million people. Plankton samples were taken biweekly from March 2014 to 2017 (n = 74 samples). Lake water was collected from a depth of 5 m with a 5 L Ruttner water sampler. The 5 L water volume of each sample was split in 300 mL for the morphospecies counting method and 4 L for HTS (see below), i.e., the analyzed water came from the same starting volume for both methods (see also the workflow scheme in Supplementary Figure S1). Additionally, net hauls (10 μm mesh size) were collected during the whole sampling period from the upper 20 m to identify, isolate and cultivate living specimens. Accompanying abiotic parameters, i.e., water temperature and oxygen concentration were measured in situ for the sampling depth (5 m) with a multiparameter probe (6600 V2, Yellow Springs Instruments, United States). Data for the description of the trophic status, i.e., total phosphorus (TP) and nitrate (NO3-N) concentrations, and total abundances of phytoplankton were obtained from the Water Supply Company Zurich. For bacterial counts, 20 mL of each 5 L sample were preserved with formaldehyde (2% final concentration), stained with SYBR Green I (Sigma Aldrich) and evaluated via flow cytometry (Cytoflex S, Beckman Coulter; ).
Identification, Abundance and Biomass of Ciliates
Determination of living ciliates from net hauls was done subsequently after sampling in the lab using a Zeiss Axio ImagerM1 microscope (magnification: 100× to 840×) and a Zeiss Discovery.V8 stereo microscope (magnification: 10× to 80×). Single or multiple cells of identified ciliates (Table 1 and Supplementary Figure S1) were isolated using a drawn glass micropipette and washed in 3–5 drops of sterile filtered lake water or medium (Volvic mineral water). Cultures were sustained by sterilized wheat grains supporting bacterial growth for bacterivorous ciliates or by a Cryptomonas strain SAG 26.80 (Culture Collection of Algae at University of Göttingen, Germany) serving as food source for algivorous / omnivorous species. Ciliate cultures were kept at 18°C and a 14 h light (20 μmol m-2 s-1)/8 h dark cycle.
Table 1
| Taxonomy | Cultivation successful | Food | Newly sequenced taxa? | Accession number GenBank | Sequence length in bp | Closest known relative in Genbank (similarity; accession number) |
|---|---|---|---|---|---|---|
| Intramacronucleata | ||||||
| Litostomatea | ||||||
| Pelagodileptus trachelioides | no | – | no | LS999902 | 757 | Pelagodileptus trachelioides (100%; AB558117) |
| Oligohymenophorea | ||||||
| Cinetochilum margaritaceum | yes | wheat grainsc | yes | LR025746 | 2832 | Urocentrum turbo (90%; AF255357) |
| Histiobalantium bodamicum | yes | Cryptomonas sp. | yes | LS999901 | 1724 | Histiobalantium comosa (99%; KU665372) |
| Stokesia vernalis | no | – | no | LS999907 | 300 | Stokesia vernalis (100%; HM030738) |
| Prostomatea | ||||||
| Balanion planctonicuma | yes | Cryptomonas sp. | yes | LS999896 | 826 | Balanion masanensis (92%; AM412525) |
| Coleps sp. mixotroph | yes | Cryptomonas sp. | no | LS999899 | 1269 | Coleps spetai (100%; AM292312) |
| Spirotrichea | ||||||
| Codonella craterab | no | – | no | LS999898 | 1362 | Tintinnopsis lacustris (100%; JQ408161) |
| Halteria bifurcataa | yes | wheat grainsc | yes | LS999900 | 1618 | Halteria grandinella (98%; AF194410) |
| Pelagostrombidium mirabile | yes | Cryptomonas sp. | yes | LS999903 | 1676 | Strombidium paracalkinsi (96%; KJ737432) |
| Rimostrombidium lacustris | yes | Cryptomonas sp. | no | LS999904 | 1582 | Rimostrombidium lacustris (99%; DQ986131) |
| Uroleptus willii | yes | Cryptomonas sp. | no | LS999908 | 1619 | Uroleptus willii (100%; EU399543) |
| Postciliodesmatophora | ||||||
| Heterotrichea | ||||||
| Stentor roeselii | yes | Cryptomonas sp. | no | LS999906 | 1289 | Stentor roeselii (100%; KP970248) |
List of the twelve ciliate species which have been isolated and partly cultivated during the investigation period (March 2014–2017).
(a) Incertae sedis, and , respectively. (b) Synonym of Tintinnopsis lacustris, . (c) Mixed bacterial suspension growing on wheat grains. Specific sequences of the V9 regions (18S rRNA gene) were determined for all listed ciliates. Five of these species were not represented in databases up to now. We intended to get as complete sequence information (up to 2832 bp) about the rRNA genes as possible.
For counting and biomass calculations, 300 mL of lake water were preserved with freshly prepared 15 mL Bouin’s solution [containing 10.7 mL saturated picric acid, 3.6 mL formaldehyde (37% stock solution) and 0.7 mL glacial acetic acid; ]. Subsamples (100–150 mL) were filtered through 0.8 μm cellulose nitrate filters equipped with counting grids (Sartorius, Germany). Then, the quantitative protargol (silver proteinate) staining (QPS) method was applied following the protocol of with slight modifications after . After the staining procedure, silver impregnated ciliates on filters were embedded in Canada balsam (Merck, Germany) providing permanent slides which keep their quality for years. Preparations were analyzed at 200× to 1,600× magnification with a Zeiss Axio ImagerM1 microscope. The identification of living and stained ciliate species was done according to , , , , ) and . We followed the classification systems by and for taxonomic affiliation of detected species. For each sample we inspected an equivalent filter area until at least 400 cells were determined and counted. Finally, the total filter area was checked for rare species which would have been overlooked by our routine counting procedure. Biomass (i.e., fresh weight) for each identified ciliate species was calculated by multiplying ciliate abundances with species-specific conversion factors published by , , ). These published factors are based on average cell size of a species, a geometric approximation of the cell form to calculate cell volume and a specific density of 1 pg μm-3 (for details see chapter 2 in ).
Sequencing of Isolated Ciliate Species
Partial 18S rRNA genes including the V9 region were sequenced either from cultivated ciliates or directly from living individuals gathered from net hauls (Table 1 and Supplementary Figure S1). Cells were washed as described above and starved for ∼24 h in sterile filtrated lake water or medium to ensure that all potential food particles were digested. Subsequently, ciliates were transferred to 1–2 μL PCR grade water which was then directly used as template for a nested PCR using the primers from ; see Supplementary Figure S2). Alternatively, the primers Euk360F () and Univ1492RE () were used for our standard PCR protocol (Supplementary Figure S2).
Conditions for the nested PCR were as follows: First denaturation at 96°C for 1 min followed by 30 cycles, each consisting of 1 min at 96°C, 2 min at 55°C, 3 min at 72°C, followed by a final elongation of 10 min at 72°C using the GoTaq® Green Master Mix (Promega). This procedure was successfully applied for the species Cinetochilum margaritaceum (Table 1). The partial 18S rDNA of the remaining eleven species was amplified using the standard PCR protocol. The conditions were as follows: First denaturation at 94°C for 3 min followed by 35 cycles, each consisting of 1 min at 94°C, 1 min at 52°C, and 2 min at 72°C, followed by a final elongation of 5 min at 72°C using the GoTaq® Green Master Mix (Promega). PCR products were purified with Agencourt AMPure XP PCR Purification Kit (Beckman Coulter) and Sanger sequenced with ABI BigDye chemistry on an ABI 3130x Genetic Analyzer (Applied Biosystems). In some cases, additional sequencing primers E528F () and SR10f () were used in order to sequence the entire V9 region.
DNA Extraction, Amplification and Sequencing
Raw water samples (2 × 2 L) were pre-filtered through a 150 μm net to remove larger zooplankton and subsequently, filtered onto a 0.65 μm membrane filter (Durapore, Millipore) using a peristaltic pump with a low flowrate of 50 mL min-1. The filters (duplicates) were directly transferred to a cryovial containing 1.5 mL RNAlater (QIAGEN), placed in the refrigerator (∼5°C) overnight and finally stored at -80°C until further processing. For DNA extraction, each filter was transferred into a lysing matrix tube (Lysing Matrix E, MP Biomedicals) and 600 μL RLT buffer and 6 μL β-Mercaptoethanol were added. The remaining liquid in each cryovial was centrifuged and discarded. To the residual pellet, 200 μL RLT buffer and 2 μL β-Mercaptoethanol were added and mixed. This mix was added to the lysing matrix tube. To break and remove cells from the filter, each matrix tube was subjected to bead-beating for 45 s by 30 Hz. Total environmental DNA was then extracted using AllPrep DNA/RNA Mini Kit (QIAGEN).
The hypervariable V9 region (about 150 bp long) of the 18S rDNA was amplified from the extracted DNA following the protocol of . The forward primer was 1391F (5′-GTACACACCGCCCGTC-3′; referring to position 1629–1644 in the Saccharomyces cerevisiae reference with accession number U53879 in NCBI’s GenBank; ) and the reverse primer EukB (5′-TGATCCTTCTGCAGGTTCACCTAC-3′; referring to position 1774–1797 in S. cerevisiae; ). The PCR protocol for V9 amplification employed an initial denaturation step at 98°C for 30 s, followed by 30 cycles of 10 s at 98°C, 20 s at 61°C, 25 s at 72°C and a final 5-min extension at 72°C. The reactions were run in volumes of 50 μL using 0.5 μL Phusion polymerase (Biolabs), 10 μL 5× Phusion GC buffer (Biolabs), 1 μL 10 mM dNTPs, 0.5 μL template DNA, 32.5 μL pure water, and 0.5 μL forward, and 0.5 μL reverse primers. Successful amplification was checked with agarose gel electrophoresis using 1.0 g agarose (Carl Roth GmbH), 100 mL TAE buffer (1×) and 5 μL ethidium bromide. From each DNA extract, triplicate PCR reactions were conducted to minimize PCR bias. Prior to purification (Qiagen’s MinElute Kit), PCR replicates obtained from the same DNA extract were combined.
From the resulting PCR products, sequencing libraries were constructed using the NEB Next® UltraTM DNA Library Prep Kit for Illumina (NEB). The quality of the libraries was assessed with an Agilent Bioanalyzer 2100 system. Eighteen samples of the time-series dataset (March 2014 – October 2014) were sequenced on an Illumina NextSeq platform, generating 250-bp-paired-end reads. All other samples (November 2014 – March 2017) were sequenced on an Illumina MiSeq platform, generating 300-bp paired-end reads. High-throughput sequencing was conducted by SeqIT GmbH & Co. KG (Kaiserslautern, Germany). We used at least 250 bp-NextSeq reads to guarantee the highest possible sequence quality (100% overlap of read 1 and read 2 sequences; please note: sequences were exceeding 150 bp in length due to multiplexing and adapters).
Pre- and Postprocessing of High Throughput Sequencing Data
Paired-end reads were merged using a custom script. The quality of the HTS data was evaluated as following: In an initial step, excessive primer overhangs were clipped with CUTADAPT version 1.18 (). Reads were then quality-filtered using the split.libraries.py command in QIIME version 1.8.0 (). Only such reads were retained which had exactly matching barcodes and primers, contained exclusively unambiguous nucleotides and had a minimum length of 90 basepairs. In the final quality-filtering step, all reads underwent de novo chimera analysis in UCHIME v5.2.236 ().
All high-quality reads were eventually de-replicated into amplicons and clustered in SWARM version 2.2.2 () using d = 1. This value d refers to a local clustering threshold instead of an arbitrary global clustering threshold. Additionally, SWARM does not depend on the input order. OTU grow iteratively by comparing each generation of assigned amplicons to the remaining reads in the dataset. An OTU is closed when no new amplicon with d or fewer differences can be assigned to the OTU. Using a custom script, an OTU contingency table was created based on the output files of SWARM (the script can be directly obtained upon request from Dominik Forster, University of Kaiserslautern).
Assignment of HTS Data to Cultured Ciliate Species
Prior to comparing HTS data to sequence data of cultured ciliate species, the OTU contingency table was normalized to the smallest number of sequences in all samples (110,100 sequences in sample March 31, 2014). Normalization was done in R version 3.5.1 (R Development Core Team, 2008) using the rrarefy command in the package vegan1. From each OTU, we extracted the seed sequence as a representative. The representatives of the HTS dataset were then compared against the sequences obtained from the twelve selected ciliate species using blastn in BLAST version 2.6.0 (). HTS sequences had to share a fragment of at least 90 consecutive nucleotides and a similarity of 97% in order to be assigned to one of these twelve cultured ciliates. Only those OTUs whose representative sequences matched our twelve selected ciliates (Table 1) were extracted from the normalized contingency table for further analyses.
Statistical Analyses: Comparison of Molecular and Morphological Data
To test if succession patterns of morphospecies counts were resembled by HTS data, a Spearman’s rank-order correlation coefficient was calculated for each of the twelve investigated ciliate species. As we performed several correlation tests on the same data set, a Bonferroni correction was applied to protect from an inflation of the alpha level. In addition, cell counts, biomass values and numbers of sequence reads per sampling date were used to compile three dissimilarity matrices (Bray Curtis) – one for each measurement. Correlation of sequence versus cell counts or biomass matrices was tested using a simple Mantel test with 10,000 permutations. All statistical analyses were performed using R version 3.4.4 (R Development Core Team, 2008).
Results
Seasonality of the Ciliate Assemblage With Special Focus on Twelve Selected Species
Lake Zurich is monomictic and water turnover usually occurs in early springtime, i.e., between February and April. Environmental parameters showed typical seasonal patterns, i.e., there were no extreme climatic years during the investigation period (Figure 1). Water temperature in the sampling depth (5 m) reached a minimum of 4.7°C in winter and a maximum of 23.7°C in summer. Maximal oxygen concentrations were measured during the phytoplankton spring bloom phases (April–May). Toward the end of the year, oxygen concentrations slightly declined, however, values in 5 m depth were still above 8.7 mg L-1 (>70% saturation) even during winter. Total phosphorus and nitrate (NO3-N) concentrations decreased owing to uptake by autotrophs toward the end of the year and increased as a result of water turnover in early spring. Maximal abundances of phytoplankton were reached during summer months, whereas numbers of total heterotrophic bacteria increased during spring and remained at a high level (ca. 4 × 109 bacteria L-1) until autumn. Total ciliate abundances followed typical repetitive annual patterns with slightly pronounced peaks during springtime and maximal values (ca. 60 × 103 cells L-1) during summer at highest water temperatures (Figure 1). During the 3 years of investigation, we detected 48 different ciliate morphotypes based on microscopic counts, and 31 of these morphotypes could be clearly assigned to morphospecies level (Supplementary Table S1). Seventeen morphotypes were assigned only to genus level, partly (i) representing yet undescribed species and (ii) morphotypes which could not be linked to a known species owing to methodological restrictions in the evaluation of QPS preparations.
FIGURE 1
With our isolation and cultivation workflow (Supplementary Figure S1), we determined specific sequences (including V9 regions of the 18S rRNA gene) for twelve ciliate species assigned to five classes from Lake Zurich that were used as seeds for HTS analyses (Table 1, accession numbers: LR025746, LS999896, LS999898-LS999904, LS999906-LS999908). Five of these species were not yet represented in public genetic databases. Several highly abundant and very frequent species could be characterized by morphospecies counting as well as HTS (Figure 2). The comparison of both methodological approaches resulted in the following outcomes: (i) Seasonal successions for the most abundant morphospecies based on counts showed remarkable correlations with HTS data (see seven top panels in Figure 2, Table 2, and Supplementary Figure S3). Even for species with striking abundance fluctuations, e.g., Cinetochilum margaritaceum, Histiobalantium bodamicum, Pelagostrombidium mirabile, a mixotrophic Coleps sp., and Rimostrombidium lacustris both methods gave highly comparable temporal dynamics. However, the correlation of reads specific for Balanion planctonicum with cells counts was rather low (see discussion for possible reasons). A moderate correlation between the two methods was observed for the second most abundant species Halteria bifurcata. Halteria counts probably comprised two barely distinguishable species, namely H. bifurcata and H. grandinella, leading to potential errors in morphospecies counts. Nevertheless, there were statistically significant positive correlations (Spearman’s rank-order) between microscopic counts and HTS data for eleven of the twelve selected ciliates, but high coefficients >0.6 were determined for only five species (Table 2). For eight species, significant correlations were found for the whole 3 years’ dataset, but not within each year (Table 2). (ii) Divergence between HTS and morphospecies counting was noticed for ephemeral and rare species, based on morphological investigations (e.g., Uroleptus willii, Stokesia vernalis). Here, HTS seemed to be a more precise method to describe their seasonality than the classical approach (Figure 2, 3 and Supplementary Figure S3). For example, HTS demonstrated the annual appearance for Codonella cratera and Pelagodileptus trachelioides in all 3 years of the study, although these species could be found only in 2 of 3 years by microscopy. The largest discrepancy between both methods was observed for Stentor roeselii. We detected this species microscopically only in 1 year (i.e., in one out of 74 analyzed protargol impregnated slides), whereas amplicon reads demonstrated the appearance in all 3 years of investigation (Figure 2).
FIGURE 2
Table 2
 |
Spearman rank order correlation analysis of abundances for the twelve selected ciliate species versus the number of reads for their specific V9 sequences.
Values were determined for all 74 samples from March 2014 to 2017, and for each year of investigation separately. Levels of significance: ∗∗∗(p < 0.001), ∗∗(0.001 < p < 0.01), ∗(0.01 < p < 0.05). P-values denote Bonferroni-adjusted significance. Color codes indicate significant correlations and the height of Spearman’s Rho coefficient. ND = not determined.
FIGURE 3
In sum, HTS showed higher percentages of occurrence for all except of one ciliate species (B. planctonicum) than microscopic counts (Figure 3). Differences between the two methods were striking for some ciliates, e.g., specific reads for R. lacustris were recorded for all investigated samples (n = 74), but we found this species only in 50% of inspected QPS preparations. Based on cell counts, we would have defined S. vernalis and C. cratera as truly ephemeral species, however, HTS data showed their appearance in 80 and 77% of samples, respectively (Figure 3).
Abundance and Biomass Values Versus Numbers of Sequence Reads
We first defined rank orders for the seven most abundant species based on numbers, and determined if these classifications were also reflected by biomasses and numbers of species-specific sequence reads (Figure 4). Congruence between number of sequence reads and biomasses was notably higher than between number of sequence reads and morphospecies abundances. However, for some species we found striking under- or overestimations of their quantitative importance in the ciliate assemblage by HTS. For example, counts identified the small B. planctonicum (∼20 × 15 μm) as the numerically dominant species throughout the investigation period, while numbers of sequence reads for this specific taxon were low compared to all other studied ciliates (Figure 4 and Supplementary Table S2). An opposite ratio was detected, e.g., for the larger species R. lacustris (120 × 70 μm), which was characterized by high numbers of sequence reads but low abundances of counted individuals (Figure 2, 4 and Supplementary Table S2). Throughout the study, highest numbers of sequence reads were obtained for H. bodamicum, a species which reached also high total biomass values (Figure 4).
FIGURE 4
For a final overview, we compared ratios of the twelve selected species to each other in relation to abundance, biomass and number of sequence reads for the whole period of investigation (Figure 5). Due to differences in the average cell sizes of these species, the contribution to biomass differed from their proportion of cellular abundances. For example, B. planctonicum contributed with 38% to the cumulative abundance of the twelve selected species, but accounted for only 11% of cumulative biomass and for 0.5% of sequence reads, respectively. An opposite pattern was observed for, e.g., R. lacustris with contributions of 0.9, 7, and 19% to cumulative abundance, biomass and reads values, respectively. In summary, ratios of species specific reads to each other better reflected the assemblage composition based on biomass determination (Figure 5). This trend was supported when comparing Bray Curtis distances between each pair of samples (n = 74) using cell counts, biomass and the number of sequence reads for each species. A Mantel test indicated a stronger correlation between distances based on the number of sequence reads and biomass per species (r = 0.425, p < 0.001) rather than that based on the number of sequence reads and cell counts per species (r = 0.231, p < 0.001).
FIGURE 5
Discussion
Strengths and Limitations of Morphospecies Based Quantifications
Since the concept of the microbial loop and its importance for food web dynamics was first formalized (), research on the role of ciliates as intrinsic components of marine and freshwater pelagic systems has received considerable attention. Today, we can profit from a “long” lasting research history in ciliate morphospecies analyses – with a special focus on plankton. Since the 1990s, new methodological approaches (; ) and excellent taxonomic compilations (; ) allowed for detailed descriptions of ciliate assemblage successions. Thus, seasonal patterns of some key planktonic species are well described, including data on environmental factors (e.g., food sources, temperature effects, top-down grazing pressure) which cause ciliates’ successions (see ; and references therein). Consequently, for any study on planktonic ciliates, a correct morphospecies identification is the prerequisite to make use of this profound knowledge. Here the quantitative protargol staining (QPS) proved as the tool of choice to quantify ciliates in combination with a relatively high taxonomic resolution (; ; ).
Due to our long work experience with this method (), we can define its pros as well as its drawbacks and limitations: (i) As preserved and stained specimens are concentrated on filters and embedded in a layer of agar first and in Canada balsam in the end, preparations have a certain thickness which may negatively affect the optical evaluation during bright-field microscopy. This is especially problematic for the identification of very tiny species, where minute differences, e.g., in their oral or somatic ciliary patterns, need to be recognized. For example, it is difficult to discriminate between small morphospecies within the genera Urotricha and Cyclidium in QPS slides (Supplementary Table S1). Thus, non-quantitative silver impregnation methods are needed additionally for an exact species identification (). (ii) The entire work procedure of QPS is a time-consuming method and commonly 5 h are needed to stain 8 samples, and usually several hours have to be spend for the evaluation of one preparation. As a consequence, we have restrictions in sampling campaigns which aim in high spatial and temporal resolutions of ciliate assemblage dynamics. (iii) All methods for the quantification of planktonic ciliates only allow for the analysis of a limited water volume. For example, we could concentrate a maximum of 150 mL per sample to obtain good preparations for the oligo-mesotrophic Lake Zurich.
The two latter drawbacks (restrictions of evaluable samples and water volumes) seem to hamper the analysis of ephemeral and rare species on the one hand, but also the description of population dynamics of abundant and fast growing representatives on the other hand. In an earlier study on the spring bloom dynamics in Lake Zurich, we found by high sampling frequency (2–4 days sampling intervals) that vernal abundance peaks of numerous species lasted for only a few days (). Such rapid population dynamics were observed in several lakes (e.g., ; ) and were linked to high growth rates of ciliates resulting in doubling times of 1–2 days (; ). Consequently, there is definitely need to investigate larger sample volumes at a high temporal resolution (e.g., daily during the warm season) to elucidate dynamics of planktonic ciliate assemblages. Here, HTS seems to be a promising supplementary approach to morphospecies-based methods.
The Choice of the V9 Region as Marker
Different hypervariable regions have been used for the monitoring of eukaryotic plankton communities in previous studies, and their utility has been debated among researchers. For a review and a comparison, we refer to . Among the three most popular marker regions (V1–V2, V4, and V9), we chose the V9 region for the following reasons: (i) While primer universality is highly similar for all three specific standard primer sets, currently, the V9 region is the best tradeoff between database coverage and taxonomic resolution (sequence variability; ). (ii) Furthermore, this region is notably shorter than the V1–V2 and V4 regions and thus, allows for a much larger sequence depth at lower costs, which is especially relevant in studies with high sample numbers. (iii) Finally, the largest environmental sequencing study for eukaryotic plankton communities available to date () used the V9 region, allowing for a good coverage of most (if not all) major eukaryotic taxonomic lineages. Future studies on environmental sequence metadata from eukaryotic plankton will benefit from surveys, which already used this genetic marker.
Correlations of HTS Data With Morphotype Counting
Translating amplicon abundance of HTS datasets into cell abundance is a critical issue in molecular ecology and environmental sequencing. In earlier studies of protistan mock communities, while all species were detected after HTS sequencing, the relative proportion of sequence types did not correlate with the cell abundances in the cell mixtures (; ). In addition, comparative studies in natural ecosystems showed that sequence abundances and cell abundances are often incongruent (; ; ). However, these studies are not conclusive, because other studies support the quantitative use of HTS data. reported adequate correlations between microscopy counts of planktonic picoeukaryotes and molecular signals after HTS. This is supported by more recent studies of and , showing that assemblage analyses based on ciliate HTS abundance data mirrors patterns obtained from morphological investigations of the same samples. argued that PCR biases (; ; ) and putative sequencing artifacts () do not affect the correlation between relative read and cell abundances. The authors demonstrated that more reads imply a higher proportion of cells. According to these authors, the use of relative read abundance as a proxy of community composition for comparative purposes is justified. The contradicting results of the mentioned studies demonstrate that relative abundances of HTS data may be interpreted quantitatively, but also warn about considerations as direct proxies for cell abundances in natural communities. The circumstances under which quantitative analyses and interpretations of HTS data are possible and reliable are still to be defined.
In our study, both methods showed clear seasonal successions for all selected species and HTS based population dynamics were in most cases remarkably correlated with morphospecies data (Figure 2 and Supplementary Figure S3). Up to now, there were only few attempts to compare temporal HTS datasets with microscopic analyses for freshwater systems. found significant correlations between cell abundance and sequence data only for two out of five investigated protistan species. conducted a survey of the ciliate assemblage from an oligo-mesotrophic lake by microscopy and pyrosequencing of the 18S rRNA gene V4 region. In the latter study, sequencing data did not mirror morphospecies abundances at all. One major reason for this discrepancy was an inadequacy of the reference gene database with extremely low coverage of typical freshwater ciliates. The authors already stated that further barcoding initiatives for freshwater ciliates (not only focused on the V4 region) were urgently needed to fully exploit the potential of HTS in ecological studies of freshwater protists. Also, a deeper sequencing strategy as achievable with pyrosequencing would be beneficial to cover the full extent of protistan plankton diversity in freshwater environments ().
As shown by our study, ciliates’ succession dynamics based on read abundances were very similar to patterns deduced from morphospecies counts. Two steps toward a better congruence of patterns obtained from microscopy and sequencing are (i) a deeper sequence coverage through Illumina sequencing technology and (ii) the availability of ciliate reference sequences obtained from the indigenous assemblage members. By this workflow, HTS helped indeed to shed light on population dynamics of ephemeral and rare species (Figure 2). The advantage of HTS over microscopy in this context was certainly the possibility to screen higher samples volumes (HTS: 4 L, morphospecies analysis: max. 150 mL), which is a prerequisite to address this question (). Since rare species may fall below the limit of microscopic detection during certain periods of the year, amplicon reads will still be found by HTS. Conversely, sequencing does not allow to distinguish between cells appearing as active trophozoites or as resting cysts in the pelagial. Thus, population dynamics of ephemeral and rare ciliate species are still a topic for further research.
Balanion planctonicum represented the only species, which was detected in more samples by microscopy than by HTS (Figure 2, 3). For the interpretation of HTS data, we could make use of only one reference sequence obtained from a clonal B. planctonicum strain (Table 1). However, this morphospecies may represent a species complex harboring cryptic taxa, which would not have been detected by our workflow, and sequencing of more clonal cultures will be needed in future. The following observation supports our thesis: By setting an identity threshold of 99% in comparison to the applied 97%, we received highly similar numbers of specific reads for all species except for B. planctonicum, which could not be detected using the higher threshold level (Supplementary Table S2). In addition, the largest discrepancy between morphological counts and amplicon read abundances was observed during summer months, when total ciliate abundances reached maximal values. As numbers of B. planctonicum specific reads were per se marginal, amplicon signals probably got lost in relation to high numbers of all other ciliate reads owing to normalization routines of HTS analysis.
Relationship Between Numbers of Reads and Abundance or Biomass Values
We examined if the ratios of reads between the twelve chosen species mirrored the assemblage composition deduced from cell counts (Figure 5). Divergence between the two methods could be expected (; ), as probably all ciliate species have an extremely high copy number of the rDNA operon per individual cell. For example, reported a value of ∼316,000 (±7,100 standard deviation) copies per cell for a peritrichous ciliate (Vorticella sp.). The authors also documented, that rDNA copy numbers per cell markedly differed among (from 3,385 to 315,786) and even within 14 investigated isolates. The latter variation may be linked to macronuclear formation or disintegration during different growth and sexual phases of the cell. For our study, species specific variations in rDNA operon copy numbers may be the reason why some, even abundant, species (e.g., B. planctonicum and C. margaritaceum) were characterized by low numbers of reads. In contrast, reads of H. bodamicum reached highest values throughout the study, although being not the most abundant species. Systematic under- or over-representations of species specific reads could also be induced by primer specific amplification patterns. However, the relative high congruence of assemblage compositions based on biomasses and numbers of reads (Figure 5) rather speaks for a relation between species cell size and rDNA copy numbers. Such a relation was previously shown for 18 algal strains representing several eukaryotic classes (). published the only study investigating the linkage between rDNA and rRNA copy numbers and phenotypic traits in ciliates. The relationship between rDNA or rRNA and cell size (volume) of studied ciliates (ribotype scaling) was in agreement with both the metabolic theory of ecology and the growth rate hypothesis, giving a quantitative framework to link cellular rRNA copy numbers with cell size, growth (activity), and biomass stoichiometry. These results may explain the observed congruence between taxonomic marker gene abundances and cellular biomass in natural ciliate assemblages.
The Need for Metabarcoding Planktonic Freshwater Ciliates
In this methodological study, we did not focus on the obvious strengths of HTS for analyzing planktonic ciliates, i.e., (i) to elucidate the beta-diversity during the seasonal cycle, (ii) to get indications of cryptic, closely related species, and (iii) to discover probably yet unknown and undescribed species (). It was our aim to compare HTS data with morphotype counts, in respect to the seasonality of abundant as well as rare species. We were aware that sequence information (18S rDNA) for most of the known morphospecies in Lake Zurich () did not exist at the start of our study. Our selection of the twelve species which were studied in more detail, partly depended on successful isolation and cultivation of specimens. Based on our morphospecies counts, we found that several species within the genus Urotricha () were as abundant as B. planctonicum during the 3 years of investigation (Supplementary Table S1). However, we did not succeed in the isolation and sequencing of Urotricha spp., and unfortunately there is not a single entry in public sequence database concerning the V9 region of the genus. Thus, we could not extract Urotricha specific reads from our HTS data, illustrating the need for further DNA barcoding of morphologically well described planktonic freshwater ciliates. Based on a comprehensive literature survey, assumed that around 180 described morphospecies may be indeed euplanktonic, and that the total number of planktonic species in one lake (alpha diversity) may vary between 50 and 100 (see also ; ). For Lake Zurich, we detected 48 different ciliate morphotypes based on microscopic counts for the epipelagic realm during the 3 years of investigation. This number seems to be rather low in comparison to the expected and partly described diversity in soils or running waters (see Table 3.7 in ). However, ciliates in pelagic realms need special adaptations on planktonic life and unique strategies to encounter sufficient food particles in often nutrient-poor environments. Taking these numbers as approximation of planktonic ciliate diversity, initiatives for DNA barcoding are manageable even when cryptic or still undescribed taxa may be expected. However, it remains challenging, as many ciliate species show very pronounced seasonal changes with distinct population maxima lasting for a few days or weeks only and being rare for most of the year.
Finally, the evaluation of HTS depends on a substantial database of reference sequences. Through linking OTUs with known morphospecies, we can use the vast knowledge about the autecology of these species. This aspect is crucial for science and especially for applied studies that are based on the identification of indicator species.
Statements
Author contributions
TP, TS, BS, and GP conceived and designed the study. GP and EB did the morphospecies counts and the sequencing of isolated/cultivated species. TS, ZQ, and DF conducted and evaluated the HTS analyses. All authors discussed the results and contributed to the manuscript concept. GP and TP wrote the manuscript and created all figures and tables.
Funding
This study was financed by the Swiss National Science Foundation (SNF, Project: 31003A-182489 and 310030E-160603/1), the German Research Foundation (DFG, Project: STO 414/13-1), the Austrian Science Fund (FWF, Project: I2238-B25) and a research grant by the Carl Zeiss Foundation to DF.
Acknowledgments
We thank Thomas Pröschold for developing the nested PCR strategy for the V9 region. We also thank Sebastian Dirren, Deborah Knapp, Eugen Loher, Daniel Marty, and Katya Pushkareva for their help in sampling and laboratory works. We are also thankful to Barbara Kammerlander, Michaela Salcher, Edward Mitchell, and Stefan Andrei for their critical comments on the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.00248/full#supplementary-material
References
1
AltschulS. F.GishW.MillerW.MyersE. W.LipmanD. J. (1990). Basic local alignment search tool.J. Mol. Biol.215403–410. 10.1016/S0022-2836(05)80360-2
2
Amaral-ZettlerL. A.McclimentE. A.DucklowH. W.HuseS. M. (2009). A method for studying protistan diversity using massively parallel sequencing of V9 hypervariable regions of small-subunit ribosomal RNA genes.PLoS One4:e6372. 10.1371/journal.pone.0006372
3
AzamF.FenchelT.FieldJ. G.GrayJ. S.Meyer-ReilL. A.ThingstadT. F. (1983). The ecological role of water-column microbes in the sea.Mar. Ecol. Prog. Ser.10257–263. 10.3354/meps010257
4
BachyC.DolanJ. R.López-GarcíaP.DeschampsP.MoreiraD. (2013). Accuracy of protist diversity assessments: morphology compared with cloning and direct pyrosequencing of 18S rRNA genes and ITS regions using the conspicuous tintinnid ciliates as a case study.ISME J.7244–255. 10.1038/ismej.2012.106
5
BachyC.GómezF.López-GarcíaP.DolanJ. R.MoreiraD. (2012). Molecular phylogeny of tintinnid ciliates (tintinnida, ciliophora).Protist163873–887. 10.1016/j.protis.2012.01.001
6
BachyC.MoreiraD.DolanJ. R.López-GarcíaP. (2014). Seasonal dynamics of free-living tintinnid ciliate communities revealed by environmental sequences from the North-West Mediterranean Sea.FEMS Microbiol. Ecol.87330–342. 10.1111/1574-6941.12224
7
BanerjiA.BagleyM.ElkM.PilgrimE.MartinsonJ.Santo DomingoJ. (2018). Spatial and temporal dynamics of a freshwater eukaryotic plankton community revealed via 18S rRNA gene metabarcoding.Hydrobiologia81871–86. 10.1007/s10750-018-3593-0
8
BeaverJ. R.CrismanT. L. (1989). The role of ciliated protozoa in pelagic freshwater ecosystems.Microb. Ecol.17111–136. 10.1007/BF02011847
9
BerdjebL.ParadaA.NeedhamD. M.FuhrmanJ. A. (2018). Short-term dynamics and interactions of marine protist communities during the spring-summer transition.ISME J.121907–1917. 10.1038/s41396-018-0097-x
10
BockC.SalcherM.JensenM.PandeyR. V.BoenigkJ. (2018). Synchrony of eukaryotic and prokaryotic planktonic communities in three seasonally sampled Austrian lakes.Front. Microbiol.9:1290. 10.3389/fmicb.2018.01290
11
CaporasoJ. G.KuczynskiJ.StombaughJ.BittingerK.BushmanF. D.CostelloE. K.et al (2010). QIIME allows analysis of high-throughput community sequencing data.Nat. Methods7335–336. 10.1038/nmeth.f.303
12
CaronD. A.GastR. J.LimE. L.DennettM. R. (1999). Protistan community structure: molecular approaches for answering ecological questions.Hydrobiologia401215–227. 10.1023/a:1003721923719
13
CarriasJ. F.AmblardC.BourdierG. (1998). Seasonal dynamics and vertical distribution of planktonic ciliates and their relationship to microbial food resources in the oligomesotrophic Lake Pavin.Arch. Hydrobiol.143227–255. 10.1127/archiv-hydrobiol/143/1998/227
14
De VargasC.AudicS.HenryN.DecelleJ.MahéF.LogaresR.et al (2015). Eukaryotic plankton diversity in the sunlit ocean.Science348:6237. 10.1126/science.1261605
15
DolanJ. R.StoeckT. (2011). Repeated sampling reveals differential variability in measures of species richness and community composition in planktonic protists.Environ. Microbiol. Rep.3661–666. 10.1111/j.1758-2229.2011.00250.x
16
DuffR. J.BallH.LavrentyevP. J. (2008). Application of combined morphological-molecular approaches to the identification of planktonic protists from environmental samples.J. Eukaryot. Microbiol.55306–312. 10.1111/j.1550-7408.2008.00328.x
17
EdgarR. C.HaasB. J.ClementeJ. C.QuinceC.KnightR. (2011). UCHIME improves sensitivity and speed of chimera detection.Bioinformatics272194–2200. 10.1093/bioinformatics/btr381
18
EdgcombV.OrsiW.BungeJ.JeonS.ChristenR.LeslinC.et al (2011). Protistan microbial observatory in the Cariaco Basin, Caribbean. I. Pyrosequencing vs Sanger insights into species richness.ISME J.51344–1356. 10.1038/ismej.2011.6
19
EggeE.BittnerL.AndersenT.AudicS.De VargasC.EdvardsenB. (2013). 454 Pyrosequencing to describe microbial eukaryotic community composition, diversity and relative abundance: a test for marine Haptophytes.PLoS One8:e74371. 10.1371/journal.pone.0074371
20
FoissnerW. (2014). An update of ‘basic light and scanning electron microscopic methods for taxonomic studies of ciliated protozoa’.Int. J. Sys. Evol. Microbiol.64271–292. 10.1099/ijs.0.057893-0
21
FoissnerW.BergerH. (1996). A user-friendly guide to the ciliates (Protozoa, Ciliophora) commonly used by hydrobiologists as bioindicators in rivers, lakes, and waste waters, with notes on their ecology.Freshw. Biol.35375–482. 10.1111/j.1365-2427.1996.tb01775.x
22
FoissnerW.BergerH.BlattererH.KohmannF. (1995). Taxonomische und ökologische revision der ciliaten des saprobiensystems. Band IV: gymnostomatea, loxodes, suctoria. informationsberichte des bayer.Landesamtes Wasserwirtschaft1/95:540.
23
FoissnerW.BergerH.KohmannF. (1992). Taxonomische und ökologische revision der ciliaten des saprobiensystems. Band II: peritrichia, heterotrichida, odontostomatida. informationsberichte des bayer.Landesamtes Wasserwirtschaft5/92:502.
24
FoissnerW.BergerH.KohmannF. (1994). Taxonomische und ökologische revision der ciliaten des saprobiensystems. Band III: hymenostomata, prostomatida, nassulida. informationsberichte des bayer.Landesamtes Wasserwirtschaft1/94:548.
25
FoissnerW.BergerH.SchaumburgJ. (1999). Identification and ecology of limnetic plankton ciliates. Informationsberichte des Bayer.Landesamtes Wasserwirtschaft Heft3/99:793.
26
FoissnerW.BlattererH.BergerH.KohmannF. (1991). Taxonomische und ökologische revision der ciliaten des saprobiensystems. Band I: cyrtophorida, oligotrichida, hypotrichia, colpodea. Informationsberichte des bayer.Landesamtes Wasserwirtschaft1/91:478.
27
FoissnerW.PfisterG. (1997). Taxonomic and ecologic revision of urotrichs (Ciliophora, Prostomatida) with three or more caudal cilia, including a user-friendly key.Limnologica27311–347.
28
ForsterD.DunthornM.StoeckT.MahéF. (2016). Comparison of three clustering approaches for detecting novel environmental microbial diversity.PeerJ4:e1692. 10.7717/peerj.1692
29
ForsterD.FilkerS.KochemsR.BreinerH.-W.CordierT.PawlowskiJ.et al (2018). A comparison of different ciliate metabarcode genes as bioindicators for environmental impact assessments of salmon aquaculture.J. Eukaryot. Microbiol.10.1111/jeu.12670 [Epub ahead of print].
30
FuR.GongJ. (2017). Single cell analysis linking ribosomal (r)DNA and rRNA copy numbers to cell size and growth rate provides insights into molecular protistan ecology.J. Eukaryot. Microbiol.64885–896. 10.1111/jeu.12425
31
GaoF.WarrenA.ZhangQ.GongJ.MiaoM.SunP.et al (2016). The all-data-based evolutionary hypothesis of ciliated protists with a revised classification of the phylum Ciliophora (Eukaryota. Alveolata).Sci. Rep.6:24874. 10.1038/srep24874
32
GellerW.BerberovicR.GaedkeU.MüllerH.PauliH. R.TilzerM. M.et al (1991). Relations among the components of autotrophic and heterotrophic plankton during the seasonal cycle 1987 in Lake Constance.Verh. Internat. Verein Limnol.24831–836. 10.1080/03680770.1989.11898860
33
GinerC. R.FornI.RomacS.LogaresR.De VargasC.MassanaR. (2016). Environmental sequencing provides reasonable estimates of the relative abundance of specific picoeukaryotes.Appl. Environ. Microbiol.824757–4766. 10.1128/aem.00560-16
34
GongJ.DongJ.LiuX.MassanaR. (2013). Extremely high copy numbers and polymorphisms of the rDNA operon estimated from single cell analysis of oligotrich and peritrich ciliates.Protist164369–379. 10.1016/j.protis.2012.11.006
35
HuseS. M.HuberJ. A.MorrisonH. G.SoginM. L.WelchD. M. (2007). Accuracy and quality of massively parallel DNA pyrosequencing.Genome Biol.8:R143. 10.1186/gb-2007-8-7-r143
36
HuseS. M.WelchD. M.MorrisonH. G.SoginM. L. (2010). Ironing out the wrinkles in the rare biosphere through improved OTU clustering.Environ. Microbiol.121889–1898. 10.1111/j.1462-2920.2010.02193.x
37
JeromeC. A.MontagnesD. J. S.TaylorF. J. R. (1993). The effect of the quantitative protargol stain and Lugol’s and Bouin’s fixatives on cell size: a more accurate estimate of ciliate species biomass.J. Eukaryot. Microbiol.40254–259. 10.1111/j.1550-7408.1993.tb04913.x
38
KisandV.ZingelP. (2000). Dominance of ciliate grazing on bacteria during spring in a shallow eutrophic lake.Aquat. Microb. Ecol.22135–142. 10.3354/ame022135
39
LaneD.J. (1991). “16S/23S sequencing” in Nucleic Acid Technologies in Bacterial Systematic, edsStackebrandtE.GoodfellowM. (New York, NY: Wiley), 115–175.
40
LiuA.YiZ.LinX.HuX.Al-FarrajS. A.Al-RasheidK. A. S. (2015). Molecular phylogenetic lineage of Plagiopogon and Askenasia (Protozoa, Ciliophora) revealed by their gene sequences.J. Ocean Univ. China14724–730. 10.1007/s11802-015-2559-3
41
LuoW.BockC.LiH. R.PadisákJ.KrienitzL. (2011). Molecular and microscopic diversity of planktonic eukaryotes in the oligotrophic Lake Stechlin (Germany).Hydrobiologia661133–143. 10.1007/s10750-010-0510-6
42
LynnD. H. (2008). The Ciliated Protozoa. Characterization, Classification, and Guide to the Literature, 3rd Edn. Berlin: Springer.
43
MacekM.ŠimekK.PernthalerJ.VyhnálekV.PsennerR. (1996). Growth rates of dominant planktonic ciliates in two freshwater bodies of different trophic degree.J. Plankt. Res.18463–481. 10.1093/plankt/18.4.463
44
MahéF.RognesT.QuinceC.De VargasC.DunthornM. (2015). Swarm v2: highly-scalable and high-resolution amplicon clustering.PeerJ3:e1420. 10.7717/peerj.1420
45
MarinB.PalmA.KlingbergM. A. X.MelkonianM. (2003). Phylogeny and taxonomic revision of plastid-containing Euglenophytes based on SSU rDNA sequence comparisons and synapomorphic signatures in the SSU rRNA secondary structure.Protist15499–145. 10.1078/143446103764928521
46
MartinM. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads.EMBnet.Journal1710–12. 10.14806/ej.17.1.200
47
MathesJ.ArndtH. (1995). Annual cycle of protozooplankton (ciliates, flagellates and sarcodines) in relation to phyto- and metazooplankton in Lake Neumühler See (Mecklenburg, Germany).Arch. Hydrobiol.134337–358.
48
MedingerR.NolteV.PandeyR. V.JostS.OttenwälderB.SchlöttererC.et al (2010). Diversity in a hidden world: potential and limitation of next-generation sequencing for surveys of molecular diversity of eukaryotic microorganisms.Mol. Ecol.1932–40. 10.1111/j.1365-294X.2009.04478.x
49
MedlinL.ElwoodH. J.StickelS.SoginM. L. (1988). The characterization of enzymatically amplified eukaryotic 16S-like rRNA-coding regions.Gene71491–499. 10.1016/0378-1119(88)90066-2
50
MikhailovI. S.ZakharovaY. R.BukinY. S.GalachyantsY. P.PetrovaD. P.SakirkoM. V.et al (2018). Co-occurrence networks among bacteria and microbial eukaryotes of Lake Baikal during a spring phytoplankton bloom.Microb. Ecol.7796–109. 10.1007/s00248-018-1212-2
51
MontagnesD. J. S.LynnD. H. (1987). A quantitative protargol stain (QPS) for ciliates: method description and test of its quantitative nature.Mar. Microb. Food Webs283–93.
52
MüllerH.GellerW. (1993). Maximum growth rates of aquatic ciliated protozoa: the dependence on body size and temperature reconsidered.Arch. Hydrobiol.126315–327.
53
MüllerH.SchöneA.Pinto-CoelhoR. M.SchweizerA.WeisseT. (1991). Seasonal sucession of ciliates in lake Constance.Microb. Ecol.21119–138. 10.1007/BF02539148
54
NakayamaT.MarinB.KranzH. D.SurekB.HussV. A. R.InouyeI.et al (1998). The basal position of scaly green flagellates among the Green Algae (Chlorophyta) is revealed by analyses of nuclear-encoded SSU rRNA sequences.Protist149367–380. 10.1016/S1434-4610(98)70043-4
55
NolteV.PandeyR. V.JostS.MedingerR.OttenwälderB.BoenigkJ.et al (2010). Contrasting seasonal niche separation between rare and abundant taxa conceals the extent of protist diversity.Mol. Ecol.192908–2915. 10.1111/j.1365-294X.2010.04669.x
56
PetrouK.NielsenD. A. (2018). Uptake of dimethylsulphoniopropionate (DMSP) by the diatom Thalassiosira weissflogii: a model to investigate the cellular function of DMSP.Biogeochemistry141265–271. 10.1007/s10533-018-0507-1
57
PfisterG.SonntagB.PoschT. (1999). Comparison of a direct live count and an improved quantitative protargol stain (QPS) in determining abundance and cell volumes of pelagic freshwater protozoa.Aquat. Microb. Ecol.1895–103. 10.3354/ame018095
58
PoschT.EugsterB.PomatiF.PernthalerJ.PitschG.EckertE. M. (2015). Network of interactions between ciliates and phytoplankton during spring.Front. Microbiol.6:1289. 10.3389/fmicb.2015.01289
59
PoschT.KösterO.SalcherM. M.PernthalerJ. (2012). Harmful filamentous cyanobacteria favoured by reduced water turnover with lake warming.Nat. Clim. Change2809–813. 10.1038/nclimate1581
60
SantoferraraL. F.GrattepancheJ.-D.KatzL. A.McmanusG. B. (2014). Pyrosequencing for assessing diversity of eukaryotic microbes: analysis of data on marine planktonic ciliates and comparison with traditional methods.Environ. Microbiol.162752–2763. 10.1111/1462-2920.12380
61
SantoferraraL. F.GrattepancheJ.-D.KatzL. A.McmanusG. B. (2016). Patterns and processes in microbial biogeography: do molecules and morphologies give the same answers?ISME J.101779–1790. 10.1038/ismej.2015.224
62
ShakyaM.QuinceC.CampbellJ. H.YangZ. K.SchadtC. W.PodarM. (2013). Comparative metagenomic and rRNA microbial diversity characterization using archaeal and bacterial synthetic communities.Environ. Microbiol.151882–1899. 10.1111/1462-2920.12086
63
ŠimekK.BobkovaJ.MacekM.NedomaJ.PsennerR. (1995). Ciliate grazing on picoplankton in an eutrophic reservoir during summer phytoplankton maximum: a study at the species and community level.Limnol. Oceanogr.401077–1090. 10.4319/lo.1995.40.6.1077
64
ŠimekK.NedomaJ.ZnachorP.KasalickyV.JezberaJ.HornákK.et al (2014). A finely tuned symphony of factors modulates the microbial food web of a freshwater reservoir in spring.Limnol. Oceanogr.591477–1492. 10.4319/lo.2014.59.5.1477
65
SimonM.López-GarcíaP.DeschampsP.MoreiraD.RestouxG.BertolinoP.et al (2015). Marked seasonality and high spatial variability of protist communities in shallow freshwater systems.ISME J.91941–1953. 10.1038/ismej.2015.6
66
SkibbeO. (1994). An improved quantitative protargol stain for ciliates and other planktonic protists.Arch. Hydrobiol.130339–347.
67
SommarugaR.PsennerR. (1993). Nanociliates of the order prostomatida: their relevance in the microbial food web of a mesotrophic lake.Aquat. Sci.55179–187. 10.1007/BF00877447
68
SommerU.AdrianR.De Senerpont DomisL.ElserJ. J.GaedkeU.IbelingsB.et al (2012). Beyond the plankton ecology group (PEG) model: mechanisms driving plankton succession.Annu. Rev. Ecol. Evol. Syst.43429–448. 10.1146/annurev-ecolsys-110411-160251
69
SommerU.GliwiczZ. M.LampertW.DuncanA. (1986). The PEG-model of seasonal succession of planktonic events in fresh waters.Arch. Hydrobiol.106433–471.
70
SonntagB.PoschT.KlammerS.TeubnerK.PsennerR. (2006). Phagotrophic ciliates and flagellates in an oligotrophic, deep, alpine lake: contrasting variability with seasons and depths.Aquat. Microb. Ecol.43193–207. 10.3354/ame043193
71
SonntagB.Strüder-KypkeM. C.SummererM. (2008). Uroleptus willii nov. sp., a euplanktonic freshwater ciliate (Dorsomarginalia, Spirotrichea, Ciliophora) with algal symbionts: morphological description including phylogenetic data of the small subunit rRNA gene sequence and ecological notes.Denisia23279–288.
72
StoeckT.BassD.NebelM.ChristenR.JonesM. D. M.BreinerH.-W.et al (2010). Multiple marker parallel tag environmental DNA sequencing reveals a highly complex eukaryotic community in marine anoxic water.Mol. Ecol.1921–31. 10.1111/j.1365-294X.2009.04480.x
73
StoeckT.BreinerH.-W.FilkerS.OstermaierV.KammerlanderB.SonntagB. (2014). A morphogenetic survey on ciliate plankton from a mountain lake pinpoints the necessity of lineage-specific barcode markers in microbial ecology.Environ. Microbiol.16430–444. 10.1111/1462-2920.12194
74
StoeckT.HaywardB.TaylorG. T.VarelaR.EpsteinS. S. (2006). A multiple PCR-primer approach to access the microeukaryotic diversity in environmental samples.Protist15731–43. 10.1016/j.protis.2005.10.004
75
StoeckT.KochemsR.ForsterD.LejzerowiczF.PawlowskiJ. (2018). Metabarcoding of benthic ciliate communities shows high potential for environmental monitoring in salmon aquaculture.Ecol. Indic.85153–164. 10.1016/j.ecolind.2017.10.041
76
StraileD. (1998). Biomass allocation and carbon flow in the pelagic food web of Lake Constance.Arch. Hydrobiol. Spec. Issues Advanc. Limnol.53545–563.
77
TanabeA. S.NagaiS.HidaK.YasuikeM.FujiwaraA.NakamuraY.et al (2016). Comparative study of the validity of three regions of the 18S-rRNA gene for massively parallel sequencing-based monitoring of the planktonic eukaryote community.Mol. Ecol. Resour.16402–414. 10.1111/1755-0998.12459
78
TirokK.GaedkeU. (2007). Regulation of planktonic ciliate dynamics and functional composition during spring in Lake Constance.Aquat. Microb. Ecol.4987–100. 10.3354/ame01127
79
TraginM.ZingoneA.VaulotD. (2018). Comparison of coastal phytoplankton composition estimated from the V4 and V9 regions of the 18S rRNA gene with a focus on photosynthetic groups and especially Chlorophyta.Environ. Microbiol.20506–520. 10.1111/1462-2920.13952
80
WearE. K.WilbanksE. G.NelsonC. E.CarlsonC. A. (2018). Primer selection impacts specific population abundances but not community dynamics in a monthly time-series 16S rRNA gene amplicon analysis of coastal marine bacterioplankton.Environ. Microbiol.202709–2726. 10.1111/1462-2920.14091
81
WeberA. A. T.PawlowskiJ. (2013). Can abundance of protists be inferred from sequence data: a case study of Foraminifera.PLoS One8:e56739. 10.1371/journal.pone.0056739
82
WeisseT.MüllerH.Pinto-CoelhoR. M.SchweizerA.SpringmannD.BaldringerG. (1990). Response of the microbial loop to the phytoplankton spring bloom in a large prealpine lake.Limnol. Oceanogr.35781–794. 10.4319/lo.1990.35.4.0781
83
YankovaY.NeuenschwanderS.KösterO.PoschT. (2017). Abrupt stop of deep water turnover with lake warming: drastic consequences for algal primary producers.Sci. Rep.7:13770. 10.1038/s41598-017-13159-9
84
ZhuF.MassanaR.NotF.MarieD.VaulotD. (2005). Mapping of picoeucaryotes in marine ecosystems with quantitative PCR of the 18S rRNA gene.FEMS Microbiol. Ecol.5279–92. 10.1016/j.femsec.2004.10.006
85
ZingelP.AgasildH.NõgesT.KisandV. (2007). Ciliates are the dominant grazers on pico- and nanoplankton in a shallow, naturally highly eutrophic lake.Microb. Ecol.53134–142. 10.1007/s00248-006-9155-4
86
ZingelP.NõgesT. (2010). Seasonal and annual population dynamics of ciliates in a shallow eutrophic lake.Fund. Appl. Limnol.176133–143. 10.1127/1863-9135/2010/0176-0133
Summary
Keywords
freshwater ciliates, morphospecies counts, V9 region of the SSU rRNA gene, Lake Zurich, seasonality, planktonic ciliates
Citation
Pitsch G, Bruni EP, Forster D, Qu Z, Sonntag B, Stoeck T and Posch T (2019) Seasonality of Planktonic Freshwater Ciliates: Are Analyses Based on V9 Regions of the 18S rRNA Gene Correlated With Morphospecies Counts?. Front. Microbiol. 10:248. doi: 10.3389/fmicb.2019.00248
Received
09 November 2018
Accepted
30 January 2019
Published
19 February 2019
Volume
10 - 2019
Edited by
Marcelino T. Suzuki, Sorbonne Universités, France
Reviewed by
David A. Caron, University of Southern California, United States; Veljo Kisand, University of Tartu, Estonia
Updates
Copyright
© 2019 Pitsch, Bruni, Forster, Qu, Sonntag, Stoeck and Posch.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Thomas Posch, posch@limnol.uzh.ch
This article was submitted to Aquatic Microbiology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.