Specificity of qPCR assays.
For all TNAs, the assessment of primer and probe specificities was performed previously with respective controls (
6,
40,
50); however, in order to confirm their specificities with other isolates, additional
in silico and experimental amplification tests were performed in our study. The results revealed that
Planktothrix-specific primers and TaqMan probes were specific for
P. rubescens and
P. agardhii but not for other species of cyanobacteria. The PC-IGS and 16S rRNA gene fragments were amplified from all
P. agardhii and
P. rubescens strains tested, whereas the
mcyA gene fragment was amplified only from MC-producing
Planktothrix strains (see Table S1 in the supplemental material). All nontarget DNAs could be amplified by targeting the cyanobacterium-specific 16S rRNA-ITS gene region but gave no products when those DNAs were amplified with
Planktothrix-specific primer-probe sets. The presence of as much as 0.0125 to 1.25 ng of
Microcystis,
Lyngbya,
Synechococcus, and
Anabaena DNAs in the qPCR mixtures did not show a significant effect on
Cq values of the samples which had only
P. rubescens DNA (
Table 1). For
mcyA, the majority of nonspecific templates did not show any amplification, but in some cases, late amplification signals (
Cq > 31) were observed for nontarget species, which could be abundant in the lake sediments (e.g.,
Synechococcus,
Microcystis, and
Anabaena). Nevertheless, the sequencing results for the
mcyA amplicons from the qPCR runs with sediment samples from 1-, 4-, and 24-cm core depths (16 sequences for each depth) confirmed the specificity of TNA-MAP (
6), demonstrating the specific amplification of the target
mcyA gene of
P. rubescens and
P. agardhii only from sediment samples.
The application of SYBR green with optimal primer sets can be as sensitive as qPCR using TaqMan probes (
35), but it is much cheaper to operate. In our study, SYBR green assay with the CYA371F-CYA783R primer set did not show sufficient specificity over the wide range of annealing temperatures examined (52 to 60°C). The alternative primer set, CSIF-373R, appeared to be highly specific for cyanobacteria, and a single amplification product with a characteristic melting temperature was observed when both genomic DNAs from a variety of algal cultures and preserved DNAs from sediment samples were analyzed. Importantly, the characteristic melting temperatures of different cyanobacterial standard strains, including
P. rubescens strains TCC 29 and TCC 38 (86 ± 0.1°C),
P. agardhii strains SAG 5.81 and CCAP 1459/15 (86.1 ± 0.1°C),
Anabaena strain TCC 79 (85.8 ± 0.2°C),
Microcystis strains HUB 5.3 and HUB 524 (85.9 ± 0.1°C), and
Synechococcus strains TCC 32 (86.2 ± 0.3°C) and TCC 173 and 175 (87.3 ± 0.6°C), corresponded to those for DNAs amplified from sediment samples (86.1 ± 0.2°C). This result confirms that specific amplification occurred with preserved DNAs from sediment samples. The sizes of the PCR products varied among the different cyanobacterial genera (275 bp to 350 bp) (
25), which resulted in small differences (<1.5°C) in the characteristic melting temperatures. When a DNA mix containing known amounts of specific and nonspecific templates was analyzed by qPCR, only specific amplification was observed, based on the presence of a single melting peak. Post-PCR analyses by inspection of an agarose gel containing each sediment sample showed no unspecific products, which confirmed the results of melting curve analyses. Primer dimers were not formed during qPCR experiments, as shown by the negative-control samples and postamplification melting curve analyses for every run. Moreover, since the application of SYBR green with relatively large (∼300-bp) qPCR products might be questionable for environmental samples, the resulting rRNA-ITS gene amplification products of qPCR runs were sequenced to obtain the ultimate proof of the specificity of this assay for sediment samples. The sequencing results (44 sequences) obtained for depths of 1, 4, 10, and 24 cm revealed a wide range of cyanobacterial diversity, and no unspecific sequences were observed. The majority of clones were affiliated with
Synechococcus spp.,
Planktothrix,
Anabaena, and
Microcystis spp. Thus, qPCR with SYBR green chemistry was found to be specific even with sediment samples, was reliable and robust enough to prevent false-positive results, and therefore could be applied successfully to enumerate total cyanobacteria in sediment samples.
Sensitivity and amplification efficiency of qPCR assays performed with culture and preserved DNA samples.
Standard curves based on serial dilutions of genomic DNA were constructed for all four qPCR assays, and they were not extrapolated beyond the ranges of dilution, which were defined as 23.2 × 10
5 to 23 and 27.4 × 10
5 to 27 copies per template for each of the target regions detected by TNAs and SYBR green reactions, respectively. The amplification efficiencies (
E) for all qPCR assays were in the range of 92 to 103% (
Table 2). Since quantification errors typically occur toward both ends of a calibration curve, cell quantification for sediment samples was achieved by direction toward the central region of the standard curves, which is found to be the most resistant against background effects (
28). A significant relationship between the
Planktothrix cell numbers estimated by microscopic counting and those estimated by TNA-16S rRNA was found (
Fig. 1C). The cell number per milliliter estimated by qPCR with DNA extracts was 1.11 × 10
7 ± 9.81 × 10
5, compared to 2.08 × 10
7 ± 0.56 × 10
6 for enumeration by direct microscopic counting. The regression equation was as follows:
y = 4.38
x − 1,184 (
R2 = 0.97;
P < 10
−4;
n = 18), where
y and
x are the log cell numbers determined by microscopic and TNA-16S rRNA methods, respectively.
Controls for evaluating the performances of the multiple assays, amplifying each target individually and comparing the results with those for the multiplex assays, were performed for both cultures and sediment samples. Resulting Cq values for all assays were highly comparable, independent of the sample type and target gene (TCC 29 had Cq values of 16.66 ± 0.45 [singleplex] and 16.39 ± 0.18 [multiplex], and B5-1 had values of 27.14 ± 0.40 [singleplex] and 26.52 ± 0.09 [multiplex], as determined by TNA-PlPC). The different qPCR mixes tested did not reveal significant differences in performance, demonstrating similar Cq values for both culture and sediment samples as determined by either TNA-PlPC or TNA-16S rRNA (data not shown). However, Rotor Gene multiplex PCR master mix, optimized specifically for Rotor Gene cyclers, had a shorter analysis time and was therefore used for all further qPCRs.
When two alternative TNAs for enumeration of total
Planktothrix populations via the PC-IGS region of
Planktothrix, i.e., TNA-PC and TNA-PIPC, were compared, the latter demonstrated a higher sensitivity of detection, showing lower
Cq values (12.80 ± 0.35 for TCC 29 and 10.86 ± 0.24 for TCC 83-1 with TNA-PlPC versus 15.38 ± 0.37 for TCC 29 and 13.52 ± 0.07 for TCC 83-1 with TNA-PC). Thus, all qPCR assays with sediment samples were performed by TNA-PlPC. Total
Planktothrix enumeration by either 16S rRNA or PC-IGS gene-targeting primer-probe sets revealed comparable results for both sediment and culture samples, showing the high reproducibility and reliability of these genetic markers. Generally,
CQ values obtained by TNA-16S rRNA were the lowest among the other TNAs for all preserved DNA samples (
Table 2). This was likely due to the multiple copy number of the 16S rRNA gene (four copies of the 16S rRNA gene are found in
P. agardhii CYA126/8). It was feasible to detect small numbers of
Planktothrix, down to 3 copies per template, by TNA-16S rRNA assay, while such detection was not reliable by TNA-PlPC. Thus, qPCR targeting the 16S rRNA gene of
Planktothrix spp. has certain advantages over the application of PC-IGS-based qPCR, since it has increased sensitivity and enables the detection of rare organisms in a sample (
34,
35,
50).
Assessment of inhibition due to coextracted materials from sediment samples.
Due to the potentially high humic acid content in sediment samples, an absence of inhibition should be demonstrated when performing qPCR with such samples. The approach applied in this case was based on the assumption that inhibitors are diluted out when a log-linear relationship is achieved between
Cq and the dilution factor (
33). No inhibition effect was noticed when preserved DNAs from sediment samples from different depths along the sediment core were mixed with
P. rubescens strain TCC 29 at any ratio, as assessed by either TNA-16S rRNA or TNA-PlPC. Further ANOVA showed no significant effect of sediment DNA addition on
Cq values at any dilution rate (
P = 0.935 and
P = 0.584 for TNA-16S rRNA and TNA-PlPC, respectively).
Cq values of all mixes were 15.62 ± 0.14 and 11.06 ± 0.33 by TNA-PlPC and TNA-16S rRNA, respectively, while
Cq values of strain TCC 29 in the absence of sediment samples were 15.46 ± 0.27 and 10.86 ± 0.51, respectively. When both culture and preserved DNAs were diluted 10-, 100-, and 1,000-fold and mixed at a ratio of 1:1,
Cq values of those mixes were not changed compared with those of TCC 29 alone at the respective dilution rate. Additionally, a
CQ value of 11.19 ± 0.73 was achieved for all nondiluted culture and sediment samples mixed at a ratio of 1:1, as determined by TNA-16S rRNA, which was close to the
CQ value for TCC 29 alone (10.86 ± 0.51). Importantly, the amplification efficiencies of all qPCR assays performed with preserved DNAs from sediment samples (0.87 to 1.16) (
Table 2;
Fig. 1B) were as high as those with the respective standards (0.91 to 1.05) (
Table 2;
Fig. 1A). Similar qPCR amplification efficiencies ensured that no PCR-inhibiting contaminants were present in the sediment samples. Based on the obtained results, we concluded that qPCRs were not inhibited, which is a prerequisite for reliable enumeration of both toxic and total
Planktothrix spp. and total cyanobacteria in preserved DNA.
Quantification of total/toxic Planktothrix spp. and total cyanobacteria in preserved DNAs from lake sediments.
Sediment sampling sites were selected based on the fact that Lake Geneva, Lake Annecy, and Lake Bourget are deep large lakes (with anoxic bottom waters over at least the last few decades) from the same region but show different histories in terms of their trophic status and, consequently, different diversity and quantities of cyanobacteria (
16,
24). Thus, if DNA is preserved in the laminated sediments of these lakes, DNA-based analyses of sediment samples should reveal these differences regarding planktonic communities, and more particularly, qPCR analyses are likely to reflect the abundance and dominance of a wide range of cyanobacterial species in water columns.
Lake sediments contained 0.24 (L1) to 5.95 (B5-10) μg DNA per gram of dry sediment. Lake Annecy sediments had smaller amounts of genomic DNA extracted than those of Lake Bourget and Lake Geneva, probably due to low productivity. The DNAs from lake sediments resulted in amplification of 1,500-bp PCR products from the 16S rRNA-ITS gene region; thus, as reported before by Coolen et al. (
9) for other ecosystems, fragmentation of this region of DNA during burial in the sediments did not occur, or at least was rare enough to favor the preservation of a sufficient quantity and quality of DNAs from several cyanobacterial groups. Moreover, the amplification products of the 16S rRNA-ITS gene region were sequenced, and around 100 sequences per core depth were obtained (O. Savichtcheva et al., unpublished data). To date, sequencing analyses with preserved DNAs from sediments have not been performed often enough, indicating the necessity and novelty of our approach. In this study, we focused mainly on demonstration of the quality and integrity of a preserved DNA fragment, rather than discussing the diversity of various cyanobacterial groups, by sequencing the 16S rRNAITS gene region. The sequencing results obtained for depths of 1, 4, 10, and 24 cm revealed a wide range of cyanobacterial diversity, and no chimeric sequences were observed. The majority of clones were affiliated with
Synechococcus spp. but also included
Cyanobium,
Anabaena,
Nostoc,
Microcystis, and
Planktothrix spp. Thus, these results demonstrate sufficient quality and fragment integrity of the preserved DNA.
We determined the amounts of DNA derived from
Planktothrix and cyanobacteria in sediment cores by qPCR. Overall, qPCR assays performed with preserved DNAs from different lake sediments revealed clear differences in the quantitative compositions of total and toxic
Planktothrix and total cyanobacteria (
Fig. 2). Thus, qPCR enumeration showed high cyanobacterial counts in Lake Bourget sediments and very low ones in Lake Annecy sediments, which is in good correspondence with available data on lake water monitoring (
16,
24,
43). Lake Bourget and Lake Geneva were indeed affected by eutrophication in the second half of the 20th century, and these systems are now recovering toward an oligotrophic status, while Lake Annecy has not been impacted by such strong eutrophication, and its oligotrophy has been preserved over the past decade (
42). Lake Geneva underwent strong eutrophication in the early 1960s, but following the implementation of phosphorus reduction measures in the 1970s, the phosphorus concentration decreased from ∼90 μg liter
−1 in the mid-1970s to ∼30 μg liter
−1 in 2005 in Lake Geneva (
53) and from 120 to 26 μg liter
−1 between 1980 and 2001 in Lake Bourget (
24).
Among all the lakes, Lake Bourget had the highest counts for total (as determined by both 16S rRNA and PC-IGS gene-based qPCR) and toxic (sample B5-1)
Planktothrix populations as well as for total cyanobacteria (sample B5-1). Total
Planktothrix quantifications based on the 16S rRNA and PC-IGS genes showed comparable results and demonstrated similar trends for all sediment cores (
Fig. 2A). Generally,
mcyA-based counts exceeded those of 16S rRNA or PC-IGS gene-based assays (
Fig. 2A). The fact that sequencing results for the
mcyA amplicons confirmed the amplification of only
Planktothrix spp. rejects the idea that the primers targeting this conserved condensation domain (21) might be able to coamplify other MC-producing taxa. At present, the highest
mcyA proportions (>100%) in sediment samples cannot be explained. From these preliminary results, it was concluded that the proportion of the
mcyA gene has remained high in Lake Bourget since the 1950s and in Lake Geneva since the 1980s. In the future, alternative qPCR assays to enumerate toxic
Planktothrix spp. should be applied to validate those results.
For Lake Bourget, the proliferation of MC-producing
P. rubescens has been reported previously (
7,
24); the dominance of this taxon in Lake Bourget has been known at least since 1996 and is associated mainly with the process of restoring this ecosystem. In Lake Annecy sediments, no
Planktothrix spp. were detected by qPCR, and
Planktothrix has never been reported for planktonic assemblages in this lake. The presence of
P. rubescens has been detected in the pelagic compartment of Lake Geneva, but generally in small amounts (compared to those in Lake Bourget).
Planktothrix spp. have been detected in Lake Geneva, especially during the years 1963 to 1985 (
16). Importantly, we observed the same trend by performing qPCR assays with lake sediments, which confirmed the lower
Planktothrix counts in Lake Geneva than in Lake Bourget.
Within the Lake Geneva sediment core, total and toxic
Planktothrix spp. were detected at 1-, 4-, and 10-cm depths (for dating, see Table S2 in the supplemental material) but not at the 24-cm core depth, while other cyanobacteria were detected down to the 24-cm layer. Interestingly, total cyanobacterial qPCR counts increased gradually from 1 to 10 cm, consistent with the period of strong eutrophication observed in the late 1980s (
53), and decreased again at the 24-cm core depth (consistent with a period of pre-eutrophication according to data from the long-term monitoring of deep perialpine lakes).
Within the Lake Bourget sediment core, particularly high
Planktothrix and cyanobacterial counts were observed at a 1-cm core depth (
Fig. 2), and they decreased with the depth but remained at high levels at 4-, 10-, and 24-cm core depths. This suggests that the presence of
P. rubescens is ancient in this lake (the 24-cm depth in sediment corresponds to the period from 1956 to 1961 [see Table S2 in the supplemental material]). Some older data on planktonic counts (performed occasionally) confirm the presence of this taxon in Lake Bourget in 1913 (
32).
In the Lake Annecy sediment core, cyanobacteria were detected by qPCR at 1- and 4-cm depths but not at 10- and 24-cm core depths (
Fig. 2). The values were generally low, with the highest quantification recorded for the 4-cm layer (corresponding to the period from 1992 to 1996).
The data obtained in our study are consistent with the past dynamics of cyanobacteria and likely reflect the temporal variations in Planktothrix and cyanobacterial species in these deep lakes. Thus, combining paleolimnology and molecular tools, the amounts of cyanobacterial and Planktothrix compositional turnover in sediments from different lakes could be analyzed, and distinct regions of the sediment core, corresponding to different trophic periods in the lake (oligotrophic to eutrophic), could be identified, inferring the recent history of the lake ecosystems.
Limitations of the use of qPCR for paleolimnology, with future perspectives.
One of the major limitations in the use of the qPCR approach for paleolimnology is the efficiency of DNA preservation, which could vary from 1 year to another, especially depending on the length or stability of anoxia. Even though it is difficult to prove that no shifts in degradation processes occurred from 1 year to the other during the past century in the three studied lakes, it could be expected that rather good preservation conditions have occurred during the last century in these three cold, anoxic aquatic sediments (
19). Paleolimnological studies performed recently on Lake Bourget and Lake Annecy (
19,
39,
42) provided converging information about the efficient storage of organic material in these lakes during the last century. Additionally, we showed that high-quality DNA could be extracted from these sediments, and since a short DNA region is targeted in performing qPCR, the qPCR method should be considered even with partial degradation of preserved DNA. Consequently, we assume that full analyses of these sediment cores by using the qPCR approach and sequencing in order to target cyanobacteria (and potentially other phytoplanktonic taxa, such as diatoms) are of special interest to reconstruct the comparative history of these lakes (
19). Moreover, to avoid the obstacles related to the different DNA preservation levels within sediment layers, it would be useful to apply the relative counts (e.g., the ratios of
Planktothrix counts to total cyanobacterial counts) instead of their absolute numbers.
On the other hand, not all cellular material is transported equally well to sediments, and there are species- or cell-specific variations in the level of (post)depositional degradation of intracellular DNA (
10). We assume that obtained DNAs represent preserved DNAs from
Planktothrix and other cyanobacteria residing in the water column rather than DNAs derived from physiologically active microorganisms residing in the sediments.
P. rubescens is characterized by its location in the metalimnic layer of the water column, and internal waves appear to have a major impact on its spatial distribution and proliferation, influencing the growth of this species by a direct impact on light availability (
12). This species is particularly well adapted for growing at low light intensities and also at low temperatures (
56). In contrast, other cyanobacterial colonies (especially
Microcystis spp.) have a benthic life cycle, where they form blooms during the summer, sink and reach sediments in autumn, remain in the sediments as vegetative cells during winter, and provide an inoculum for the water column in spring (
30). For
Planktothrix, such a benthic life cycle has not been reported. All of these factors could introduce a certain level of discrepancy and error to the qPCR-based approach performed with preserved DNA, especially where comparisons between different phytoplanktonic groups (e.g., diatoms, cyanobacteria, and chlorophyceae) are concerned. We should thus be cautious in using DNA as a quantitative indicator of past phytoplankton community structure. The combined application of DNA and traditional biomarkers, such as lipids or pigments, would also provide converging information about the importance of phytoplankton groups stored in sediments.
Further studies on the application of qPCR to sediment samples, as well as development of the data set by sequencing of the 16S rRNA-ITS gene region, should be conducted. Only a few characteristics about cell viability of cyanobacteria in sediments have been studied yet; thus, the identification of the active microbial population within sediment samples can be addressed by targeting rRNA content directly rather than targeting DNA. Moreover, the preservation (during sedimentation processes and after deposition on sediments), origin, and size of DNA in sediments should not be neglected. Even though it is expected that similar preservation conditions might occur in various anoxic aquatic sediments, further work examining whether the DNAs from different species are equally well preserved should be performed.