Skip to main content

Shifts among Eukaryota, Bacteria, and Archaea define the vertical organization of a lake sediment



Lake sediments harbor diverse microbial communities that cycle carbon and nutrients while being constantly colonized and potentially buried by organic matter sinking from the water column. The interaction of activity and burial remained largely unexplored in aquatic sediments. We aimed to relate taxonomic composition to sediment biogeochemical parameters, test whether community turnover with depth resulted from taxonomic replacement or from richness effects, and to provide a basic model for the vertical community structure in sediments.


We analyzed four replicate sediment cores taken from 30-m depth in oligo-mesotrophic Lake Stechlin in northern Germany. Each 30-cm core spanned ca. 170 years of sediment accumulation according to 137Cs dating and was sectioned into layers 1–4 cm thick. We examined a full suite of biogeochemical parameters and used DNA metabarcoding to examine community composition of microbial Archaea, Bacteria, and Eukaryota.


Community β-diversity indicated nearly complete turnover within the uppermost 30 cm. We observed a pronounced shift from Eukaryota- and Bacteria-dominated upper layers (<5 cm) to Bacteria-dominated intermediate layers (5–14 cm) and to deep layers (>14 cm) dominated by enigmatic Archaea that typically occur in deep-sea sediments. Taxonomic replacement was the prevalent mechanism in structuring the community composition and was linked to parameters indicative of microbial activity (e.g., CO2 and CH4 concentration, bacterial protein production). Richness loss played a lesser role but was linked to conservative parameters (e.g., C, N, P) indicative of past conditions.


By including all three domains, we were able to directly link the exponential decay of eukaryotes with the active sediment microbial community. The dominance of Archaea in deeper layers confirms earlier findings from marine systems and establishes freshwater sediments as a potential low-energy environment, similar to deep sea sediments. We propose a general model of sediment structure and function based on microbial characteristics and burial processes. An upper “replacement horizon” is dominated by rapid taxonomic turnover with depth, high microbial activity, and biotic interactions. A lower “depauperate horizon” is characterized by low taxonomic richness, more stable “low-energy” conditions, and a dominance of enigmatic Archaea.


The continuous deposition of organic and inorganic particles to sediments is an important process in all aquatic ecosystems. Approximately one third of the terrestrial organic matter (OM) that enters freshwater is sequestered in sediments [1], although the total amount of OM that reaches the sediments is much greater than the amount that is actually sequestered [2]. This is because microbial activity is responsible for the cycling of carbon, including methane emission [3]. In lake sediments, a proportion of newly settled OM is rapidly recycled and subsequently transformed into secondary compounds, resulting in a distinct uppermost sediment zone of high heterotrophic activity [4, 5]. This is thought to lead to the structuring of microbial communities along environmental gradients that are much steeper than those in marine sediments, with narrower vertical sequences of electron acceptors [6]. The nature of this gradient influences the carbon, nitrogen, and sulfur cycles [79] and potentially affects the microbial community structure [10].

In contrast to the wealth of studies on marine sediments (cf. the 65 studies of [11]), few studies have examined the vertical microbial community structure of freshwater sediments (e.g., [1216]). The community of sediment microbes was thought to be dominated by Bacteria, together with a smaller fraction of methanogenic Archaea (reviews in [6, 17, 18]). This view has been challenged by the recent discovery of an abundance of non-methanogenic Archaea in marine sediments [19, 20]. They are assumed to be adapted to low-energy environments, and at least one lineage seems to be specialized in inter alia amino acid turnover [20]. This discovery has led to a revised perception of microbial communities in marine sediments, where Archaea appear to be as abundant as Bacteria and increase in relative abundance with sediment depth [11]. Data on sediment Archaea in freshwater are scarce, and the causes of the significant variation observed among studies remain largely unknown (e.g., [14, 2123]).

Prokaryotic activity, biomass, and cell numbers decrease with depth in many freshwater and marine sediments (e.g., [24, 25]), although other studies report relatively constant proportions of active cells with depth and find no accumulation of dead cells in deeper sediments [5, 26]. Despite the continuous presence of vegetative cells and resting stages, recent studies of marine systems indicate that the majority of microbial cells in energy-deprived horizons consist of microbial necromass [27, 28] and the proportion of living organisms decreases with the increasing age of the sediment [29]. The vertical, progressive transformation of OM and depletion of electron acceptors may eventually lead to an extremely low-energy environment in deeper sediment layers with very low growth rates similar to sub-seafloor sediments [30].

A basic mechanism thought to explain the vertical distribution of microbes is simply the one-way input of new organisms attached to OM that sinks from the water column. We hypothesize that this process would result in two simplified, competing structural models, wherein the microbial community (1) consists exclusively of sinking colonizers, with the result being a fully nested community structure in which the community gradually changes from a complex and rich community at the surface to an increasingly depauperate community with increasing sediment depth, dominated by progressive cell death, or (2) is structured by niche specialists at various layers that are well adapted to the specific environmental conditions including redox gradients, OM, and nutrient (C, N, P) concentrations. These two models are suitably analogous to the recently developed β-partitioning of the total β-diversity of a community, in which the taxonomic turnover is mathematically separated into richness and nestedness components (see [31] for a theoretical framework and [32] for applications).

The decomposition rate of settled or buried pelagic dead and living organisms is thought to depend primarily on the activity of the indigenous microbial community rather than on chemical processes (e.g., depurination; [33]). The decrease in DNA with depth that has been reported for freshwater sediments (e.g., [34]) is likely to result of nucleic acid degradation of dead organisms, particularly eukaryotes, whose biomass also decreases with depth. As a result, the decomposition of buried organisms should be a function of the active community which is itself buried over time. Vice versa, temporal patterns of sedimentation will also influence the active microbial community, for example by shifting the redox gradient. Historically changing lake conditions are recorded in lake sediments as DNA and as chemical parameters (e.g., [34]). An important question that remains is how decomposition processes within the sediment redox gradient are related to the burial of OM, eukaryotes, and prokaryotes [35, 36].

We examined the biogeochemical properties and microbial community composition (Eukaryota, Bacteria, and Archaea [37]) of sediments in the oligo-mesotrophic hardwater Lake Stechlin in northeast Germany. Our aims were to evaluate (1) whether microbial communities were nested or structured (to test the competing models, above), (2) how sediment parameters reflecting “present” and “past” conditions influence the overall community structure (Table 1), and (3) whether recently reported vertical patterns of marine Archaea [11] can predict those observed in freshwater sediments. We took four replicate 30-cm sediment cores from ca. 30-m water depth. 137Cs dating indicated the cores include sediments deposited over the past ca. 170 years.

Table 1 Definition of “present” and “past” sediment parameters


Sediments in the cores were black in color, with no visible lamination. Water content was 93–97%. No macrozoobenthic organisms were visible, although DNA metabarcoding (see below) detected the presence of nematodes in addition to microbes (Additional file 1). There was an exponential increase of dissolved refractory carbon with sediment depth (fluorescence index (FI) = 1.69–2.01; Fig. 1) across all four cores, indicating the enrichment of fulvic acids [38]. Prokaryotic cell numbers were on average 1.8 ± 0.5×109 ml−1 wet sediment and were highest in the upper sediment layers. Bacterial biomass production as carbon (BPP-C) (range 0–282 μg C ml−1d−1) decreased rapidly with depth, approaching zero below 10 cm. Total DNA concentration (range <0.3–17.6 μg ml−1 sediment) was negatively correlated with FI (r=−0.886) and followed an exponential decay function. DNA half-life was inferred to be t 1/2=22 a (corresponding to 5.4 cm; f(DNA) = 13.9 × e −0.128x,r 2=0.81). RNA content was lower than DNA content in all layers, with DNA:RNA ratios ranging from 2.3 at the surface to 20.8 at 20-cm depth (Fig. 1). The sediment exhibited a typical electron acceptor sequence (Fig. 1) with a mean oxygen penetration depth of 4.6 mm (SD 1.4). Nitrate and nitrite were immediately depleted at the sediment surface, sulfate approached a constant minimum concentration after 5 cm, soluble reactive phosphorous (SRP) and ammonium (NH\(_{4}^{+}\)) increased with sediment depth, N2O gas was not detected, CH4 increased linearly with depth, and CO2 exhibited minima at the surface and at a depth of 10 cm (Fig. 1). More detailed profiles of all measured parameters can be found in the supplemental material (Additional file 2).

Fig. 1

Depth profiles of microbial community clustering and key biological and chemical characteristics of the sediment cores taken from 30-m depth in Lake Stechlin. Microbial communities were clustered by similarity (average clustering) into three groups (a, b, c) corresponding to different depth horizons (upper-left panel). Dates (y-axis) were calculated using 137Cs measurements. Values are means (±1 SE) from four replicate cores. Parameters and units: FI = fluorescence index; cells [ 106 ml−1]; BPP-C = bacterial protein production in carbon [ μg C ml−1d−1]; DNA extract [ng μ l −1]; the shared chao index (R vegan package, [106]); low sim. = proportion of sequences with no close relative [% of sequences]; SRP [ mg l−1]; NH4 [ mg l−1]; SO\(_{4}^{2-}\) [ mg l−1]; DOC [ mg l−1]; CH4 [ μmol l−1]; CO2 [ mmol l−1]; P [ mg g−1 dry weight]; Fe [ mg g−1 dry weight]; Pb [ mg g−1 dry weight]. Additional information for all measured variables per individual core are provided in the Additional files 2 and 12

Total taxon richness across the 60 samples was estimated (Chao) to be 8545 (SE = 173) operational taxonomic units (OTUs). The proportion of sequences with no close relatives in the SILVA reference database (<93% sequence similarity using BLAST) was highest at a depth of 9-10 cm (40±4%). An overview of eukaryotic, bacterial, and archaeal OTUs recovered in each sample layer is provided in Additional file 1.

When pooling all replicates for analysis of the community matrix with a cluster analysis, the microbial communities were grouped into three major clusters corresponding to depths of 0–5, 5–14, and 14–30 cm, with a pronounced separation at 14 cm (Fig. 1). The sediment communities in layers between 14–30-cm depth were more similar to one another (>65% of community structure) than layers in the upper two clusters (<50%). All α-diversity indices (inverse Simpson, evenness, and estimated Chao index) decreased with depth (Fig. 1). This pattern also occurred in the rarefaction analysis of the Hill indices [39] and confirmed a significant separation of the three depth clusters by major taxonomic changes (Fig. 2). Community turnover (distance) increased with depth, following a distance decay curve and approaching a distance of 1 (i.e., no shared taxa) for the comparison of the lowest layer (30-cm depth) with the surface (Table 2). Upon partitioning the β-diversity among sample layers into taxonomic richness and replacement effects [32], taxonomic replacement was consistently high and was significant for multiple sample layers above 12 cm (Table 2). In contrast, the effect of richness increased with depth (R 2=0.96,F=230.8,dF=11, see Additional file 3) and it was significantly elevated in the deepest layer (26–30-cm depth) (Table 2).

Fig. 2

Hill diversities of the three depth horizon clusters. Absolute (upper panel) and average (lower panel) Hill diversity for the three Hill numbers q resembling richness (0), exponential Shannon Index (1), and Inverse Simpson index (2)[39]. The absolute data is based on all sequences obtained for each horizon; the average diversity is based on all sequences from each sampling depth, grouped as horizons, normalized to a sampling coverage of 0.9. The average diversity was significant different between depth horizons for all Hill numbers, ANOVA: F(2,12)=27.9,59.8,67.5 for q=0,1,2; respectively; p<0.001

Table 2 Microbial community turnover and richness for each sampling layer

The OTUs that were the most influential in structuring the microbial community across our 15 sediment layers were identified by calculating species (OTU) contribution to β-diversity [32]. This analysis identified 96 “structuring” OTUs (see Additional file 4) whose identity reflected the interrelationship of domains. The most influential phyla were Euryarchaeota and Thaumarchaeota (Archaea) as well as Chloroflexi, Proteobacteria, and Phycisphaerae (Bacteria) (Additional file 4). Structuring OTUs included both redox-dependent groups (13% of structuring OTUs could be clearly assigned to redox processes by their classification, e.g., Nitrospiraceae, Desulfobacteraceae, and Methylococcales) and redox irrelevant groups (e.g., Eukaryota, Bacteriovoraceae). A number of the 96 structuring OTUs were significantly elevated in one or more zones (Fig. 3). Eukaryotic and bacterial lineages were characteristic for the uppermost cluster (cluster a in Fig. 1), whereas Archaea and Bacteria were elevated in the lowest zone (cluster c) (Fig. 3). For cluster b, only two structuring archaeal OTUs were identified. The residual OTUs from cluster b were significantly elevated either in the upper two clusters (mainly Bacteria) or in the lower two clusters (mainly Archaea). Only one structuring OTU (Methylococcales) was significantly different in its relative abundance in all three clusters (Fig. 3).

Fig. 3

Overview of the sediment structure in Lake Stechlin. The cluster analysis separates three depth horizons: the redox-stratified zone (0–5 cm), which includes a thin layer of oxygen. A few fauna species exist in this zone, i.e., Nematoda, Gastrotricha, and microeukaryotes (e.g., Ciliophora), in addition to large numbers of highly active Bacteria. Below 5 cm, where 50% of the DNA is already decomposed, the system enters the transition zone. This zone is situated below the sulfate-methane transition. Below 14 cm, we find the depauperate horizon, which extents in the deeper sediment, in which Archaea dominate the community. In an extrapolation of the richness component of the community structure, the loss of richness would completely dominate (100%) the microbial community at 1-m depth (approx. 500 a). Following the decay curve of the DNA, 99.99999% of the DNA would be transformed at that depth. On the right side, the ten most structuring OTUs (from Additional file 4) are listed, which were significantly elevated in the corresponding horizon (only results with p<0.01 in the Tukey HSD post hoc test were included). The brackets ab and bc mark those OTUs that were elevated in the upper two or lower two zones, respectively. Only two OTUs were elevated in the transition zone. The gray box marks the single taxon that was significantly different in all three horizons. Taxon names are color coded according to their classification or phototrophy if applicable: phototrophic organism (green), Eukaryota (black), Bacteria (red), and Archaea (blue)

Sequence proportions of Archaea, Bacteria, and Eukaryota (A:B:E) shifted from 10:70:20 at 0 cm to 50:50:0 at 10 cm and 60:40:0 at 30-cm depth (Fig. 4 a). The eukaryotic proportions were correlated with DNA concentration (r=0.869) and decayed exponentially with depth. Multiple linear regression could predict DNA concentration as a function of the occurrence of Eukaryota (75.6% of the variation) and Bacteria (10.0% of the variation; model: R 2=0.856,p<0.001, Additional file 5). A multivariate ordination of all samples confirmed the strong vertical gradient in microbial community structure, reflected in the distance between the surface and deep sediments on axis 1 (Fig. 4 b, Mantel correlation: r=0.735,p<0.001). The three distinct clusters at different depths (above) were recovered using adonis (Fig. 4 b, F=12.3,p=0.0005). The variance in cluster c was reduced compared to the other clusters (betadispersal analysis: Tukey’s honest significant differences between groups, p<0.01), confirming the greater similarities seen in the previous cluster analysis. Microbial community structure was correlated with sediment parameters representative of both “present” (Mantel correlation: r=0.527,p<0.001) and “past” (r=0.459,p<0.001) conditions (see Table 1). These two parameter sets were nearly orthogonal in ordination (Fig. 4 b). Apart from the betadispersal analysis, we came to the same conclusion when we applied weighted phylogenetically based UniFrac distances instead (significant structuring along the depth gradient, significant separation of the three depth clusters, significant correlation with the “present” and “past” principal components with comparable effect sizes; see Additional file 6).

Fig. 4

Microbial proportions and community structure. a Depth profiles of the microbial community (Eukaryota, Bacteria, and Archaea) presented as relative proportions to each other, which was determined by relative pyrosequencing reads per microbial fraction. b NMDS ordination of the vertical sediment microbial community structure. The clusters from Fig. 1 are presented as standard deviation around the group centroid. The color scale of the dots represents sediment depth. Community composition was highly correlated with sediment depths (on axis 1; Mantel correlation: r=0.735,p<0.001). The three distinct clusters at different depths were significantly different (adonis, p=0.0005) with a reduced variance in cluster c (Tukey’s honest significant differences between groups, p<0.01). The principal components of “past” and “present” environmental parameters were correlated with the community matrix (Mantel correlation: “present” r=0.527,p<0.001 and “past” r=0.459,p<0.001) and with the ordination (envfit: “present” R 2=0.658,p<0.001 and “past” R 2=0.547,p<0.001)

In order to examine links between the two parameter types and microbial community structure, we used fuzzy set analysis to test whether past parameters were correlated with the richness component and present parameters with the replacement component of the microbial community. We partitioned the whole dataset into replacement and richness matrices and correlated these with the present and past parameters. The richness community submatrix was strongly correlated with the past parameters (two-dimensional fuzzy set ordination with the first two axes of the PCA, r=0.99), and the replacement community submatrix was correlated with the present parameters (one-dimensional fuzzy set ordination with the first axis of the PCA, r=0.65,p<0.001).


Comprehensive studies of the microbial communities and physical-chemical characteristics of freshwater sediments are scarce, and general concepts are often transferred from marine systems without validation. This is despite the fact that salinity and sulfate concentrations are very different in the two ecosystem types [40, 41] and can have a profound influence on biological communities and biogeochemical processes. We addressed fundamental questions regarding the vertical structure, organization, and inter-relationships among microbial communities and biogeochemical parameters in freshwater sediments and establish a structural model that can potentially be applied to other aquatic sediments.

Vertical organization of Lake Stechlin sediments

The sediment habitat is thought to be autonomous in terms of species richness and community structure [17, 18], despite the constant colonization by microbes that descend from the water column with sinking organic particles (Additional file 4). In Lake Stechlin, we observed a high species (OTU) β-diversity with depth in lake sediment, leading to nearly complete taxonomic turnover of the microbial community within 30-cm depth. Such high turnover may be a common feature of vertical sediment profiles and has been reported for bacterial taxa in coastal marine sediments [42], for marine Archaea and Bacteria [43], and for freshwater Archaea [21]. Previous studies of freshwater sediments that reported moderate species turnover were restricted to low-resolution methods [15, 21, 44]. Our study differed from previous efforts in that most studies have focused on either Bacteria or Archaea and not on all three domains simultaneously, and none of the previous marine or freshwater studies have partitioned β-diversity into richness and replacement components. The former allowed us to examine whole-community patterns and potential interactions, and the latter allowed to distinguish taxonomic changes that result from microbial activity from those related to sediment burial. We found the microbial community to be clearly delineated into three distinct clusters, each spanning multiple sampling layers (details in Additional file 7). Based on the significant contribution of taxonomic replacement (Table 2) and measurable microbial production in the upper two clusters (clusters a, b; Fig. 1), we term this part of the sediment the “replacement horizon” (Fig. 3). We term the lower layers, comprising the deepest cluster (cluster c; Fig. 1) the “depauperate horizon,” based on the importance of taxonomic richness (as opposed to turnover) in structuring the community and based on the constancy of most sediment biogeochemical parameters at these depths (Fig. 3, Additional file 7). Both horizons are discussed in more detail below.

The replacement horizon We define the replacement horizon (Fig. 3) as comprising the two microbial community clusters in which taxonomic turnover was most pronounced (a, b in Fig. 1). This horizon is further subdivided by the sulfate-methane transition into what we term “redox-stratified” and “transition” zones (Information Box, Additional file 7). Bacterial activity was highest in the redox-stratified zone, where most of the settled OM is readily available. Sulfate was depleted below 5-cm depth, and therefore, most redox processes will take place above. The majority of freshwater sediment studies examine this zone in great detail (e.g., [4, 5, 10]), including the identification of redox processes at the millimeter scale [45, 46]. Active decomposition leads to high prokaryotic cell numbers close to the sediment surface [5, 25, 47] and decreasing abundance with depth [15]. We also observed an initial loss of many taxa in the highly active oxycline—indicated in our data by an outlier to the richness community component (Additional file 3)—which may have been intensified by active grazing by ciliates and copepods. Important methane oxidation processes occur in the transition zone (below the sulfate-methane transition)(cluster b, Additional file 7) [48]. At 10 cm, the CO2 minimum could indicate the start of hydrogenotrophic methanogenesis (Fig. 1) as previously described for profundal Lake Stechlin sediments [49]. The archaeal contribution rose constantly in the transition zone, while Bacteria and Eukaryota decreased and the general activity measures declined rapidly.

According to DNA concentrations, more than 85% of settled organisms were decomposed within this replacement horizon, which spanned approximately 60 years. The decay of DNA, in combination with an enrichment of fulvic acids with increasing depth (FI values), confirms the assumption that DNA can serve as a proxy for the buried OM. We conclude that there is a gradient of OM quality in addition to a gradient of electron acceptors in the replacement horizon. This may facilitate the stratification of active microbial taxa with depth because OM can simultaneously serve as an electron donor and acceptor [50, 51]. OM quality can also modulate microbial redox processes [52], with apparent consequences for carbon turnover rates [36, 51].

The depauperate horizon The depauperate horizon (Fig. 3) was characterized by high concentrations of methane (CH4) and CO2, a dominance of Archaea, and low diversity compared to the replacement horizon. Microbial community composition was more nested, with a 9% relative nestedness here compared to only 2% in the replacement horizon. By entering this horizon, the DNA:RNA ratio doubled and Archaea replaced Bacteria as the dominant microorganisms. We believe this reflects an increase in the number of microbes entering a stationary state below this depth, where cell maintenance predominates over cell synthesis due to the low availability of terminal electron acceptors. This is analogous to what has been suggested for cells in low-energy marine environments in the deep sub-seafloor sediment [27, 30, 53]. The variability in community composition was very low across replicates in the depauperate horizon. The nestedness suggests the gradual disappearance of taxa with burying age and a richness component of turnover steadily increasing to more than 20% (Table 2). If the richness component was to further increase in a linear manner, it would be the sole factor structuring the community composition deeper than 1 m (the total sediment depth of Lake Stechlin is 6 m). It is intuitive that the richness component may be a function of the burying time and that it represents the fading signal of preserved organisms. It remains unclear why it does not follow an exponential decay function analogous to that for DNA.

Potential causes of the high taxonomic replacement Many “present” parameters changed rapidly with depth, particularly in the replacement horizon (e.g., DNA, FI, BPP, electron acceptors), and this was a likely driver of the high degree of taxonomic turnover. Several mechanisms could be responsible for these patterns, namely cellular turnover and random appearances. In cellular turnover, taxonomic replacement is potentially caused by cell synthesis, lysis, and recycling of dormant cells, which are assumed to be high in sediments [54, 55], particularly viral lysis [56, 57]. We found indications for cellular recycling caused by the predatory Bacteriovoracaceae (cf. [58]), which was one of the structuring bacterial lineages identified in Additional file 4. Another potential mechanism—one that may be the most important in the depauperate horizon—is differential cell replication. The resources for cell maintenance and growth should depend on cell size and complexity. This means that small cells, such as nano-Archaea (e.g., Candidatus Parvarchaeum), should have a selective advantage because they can continue to grow under conditions in which larger cells must switch to cell maintenance. This could be one explanation for the observed drop in evenness in the depauperate horizon.

In the mechanism of random appearances, the appearance of taxa may be due to the disappearance of others. Because high-throughput sequencing methods produce relative (rather than absolute) data, it may superimpose proportions over quantities. For example, the initial decay of eukaryotes may have opened a niche for previously hidden rare taxa. Further, if there was no growth in the sediment, lineages that are potentially better suited for long-term survival than others would appear, such as spore-forming Bacteria (Firmicutes). However, we (and others: [27, 34]) did not observe an enrichment of this lineage with depth. In addition, the use of replicate cores in combination with our conservative stripping (see the “Methods” section) should have removed most of the random effects. The mostly constant cell numbers with small local maxima and the observed shifts in the evenness support a non-random stratification of sediment communities including cell replication. While the cellular reproduction probably approaches stagnation for most microbes in the depauperate horizon, the slowly shifting redox conditions across seasons and years may be conducive to colonization of the replacement horizon by different niche specialists. The low sedimentation rate of Lake Stechlin (ca. 2 mm per year, as determined by 137Cs dating at the Federal Office for Radiation Protection, Berlin [courtesy of U.-K. Schkade], or 0.4–2.1 g m −2 d −1, as determined by sediment traps [59]) may mitigate this stratification, and the scale of the horizons may differ in systems with higher or lower sedimentation rates.

Burial processes and microbial activities

The microbial sediment community appeared to be highly indigenous, and yet the constant arrival of sinking OM could bury the microbial community. Indeed, buried DNA and organisms preserve historical plankton communities that can be indicative of past conditions of the lake ecosystem [34, 60, 61]. These past environmental conditions are also partly preserved as particulate matter, which is relatively conservative. We found several of these conservative “past” parameters to correlate well with the present-day communities. Although previous studies have found sediment parameters to influence community patterns in marine systems [43], the parameters were not separated into a present-past context, and microbial community turnover was not partitioned into richness and replacement components. Our results revealed that the richness component could be largely explained by the first two principal components of the past parameters. In contrast, the replacement component was not fully explained by the present parameters, indicating that sources of variation other than environmental parameters are important, such as biotic interactions. Strong biotic interactions have been identified in a vertical profile of a meromictic lake with a comparable chain of redox processes as they occur in sediments [62]. Deep sediment layers may offer low-energy niches that favor a large variety of syntrophic microorganisms [63, 64]. The Dehalococcoidales (Chloroflexi [65, 66]) and the Miscellaneous Crenarchaeotic Group (MCG, [67, 68]) are promising candidates for such hybrid forms of energy harvesting and were among the most influential lineages in our data set (Additional file 4). The MCG co-occurred with Dehalococcoidales, similar to what was found in methane hydrate-bearing sediment in Lake Baikal at >1500-m water depth [69]. Another indication that biotic interactions are important in Lake Stechlin sediments is the appearance of the Candidatus Parvarchaeum as a structuring lineage (Fig. 3). This lineage can exhibit cell-to-cell coupling that allows for thermodynamic processes that would otherwise not be possible [70].

Archaea in freshwater sediments

The fact that Archaea can be very numerous in freshwater sediments and even dominate microbial communities is a rather new discovery, and data from comparative studies are lacking. Their recovery rate in relation to bacterial sequences or cell numbers varies between 3–12% ([71], cells), 5–18% ([21], qPCR), and 14–96% ([14], qPCR), depending on the lake, sampled sediment horizons, and methods employed. In most cases, only surface sediment samples have been considered (e.g., [72], 1% of cells), and the few studies involving vertical profiling to date are ambiguous in finding an archaeal depth gradient. Our results and those from cell counts from Lake Biwa (Japan) of [71] suggest an increase in the proportion of Archaea with sediment depth. However, the results obtained by quantitative PCR for Lake Taihu (China; [14]) and Lake Pavin (France; [21]) did not report such a relationship. Archaea have, on average, compact genomes [73] and a lower ribosomal copy number than Bacteria [74], which may lead to underestimates of archaeal abundance. Similar to our results, [21] found three sequential depth clusters in the archaeal community structure within the first 40 cm, defining an intermediate layer between 4 and 12 cm. Next to well-described methanogenic Archaea, we mainly recovered archaeal lineages with no clear functional assignment thus far (similar to [21]), i.e., primarily the MCG (potentially methanogenic, [68]) and the Marine Benthic Group D (MBG-D). Both groups are among the most numerous Archaea in the marine sub-seafloor, and they are thought to metabolize detrital proteins ([20], discussed in more detail below). Interestingly, we also identified a MBG-B as structuring OTU for the transition zone (Fig. 3), a group which was recently described as eukaryotic progenitor from a hydrothermal vent field (Lokiarchaeota, [75]). Several MCG OTUs belonged to the top structuring taxa. MCG was recently named as Bathyarchaeota by [76] for its deep-branching phylogeny and its occurrence in deep subsurface environments—environmental conditions that our cores (30-m water depth and 30-cm length) did not meet.

Our results suggest that the specific niche adaptation of these microbes is not necessarily related or restricted to the deep biosphere but rather to a cellular state of “low activity” [77]. In this context, it is interesting that single MCG OTU sometimes dominated the community in the deep horizons (up to 34% in core D at 26–30 cm), resulting in a reduced overall evenness and a shift of the residual taxa to the rare biosphere, contrasting the potential random effects as discussed above. Another intriguing observation is the considerable overlap of archaeal and partially bacterial lineages between our study and deep-sea environments. Consequently, typical marine lineages (e.g., Archaea in Additional file 4: MGI, MCG [Bathyarchaeota], MHVG, DHVEG-1, DHVEG-6 [Woesearchaeota], DSEG, MBG-A, MBG-B [Lokiarchaeota], MBG-D, MBG-E) are not as “marine” or as “deep-sea” as previously thought. Given the high cost of deep-sea research [30], freshwater sediments might literally pose a row-boat alternative for research questions targeting these “remote” and “extremophile” microorganisms.

Study limitations and perspectives

When we look at systems that sequester carbon, it is important to keep relic DNA in mind. Relic DNA is defined as DNA residuals that remain in the system after cell death. Its presence can inflate richness and misrepresent relative abundances in some types of soils when analyzed with DNA metabarcoding [78]. The very few aquatic studies that investigated relic DNA reported large amounts of extracellular DNA in marine and freshwater sediments, with fragment sizes of up to 10 kb, but with low amplification success [79, 80]. The degradation of extracellular DNA varies widely in different environments and is dependent on the amount of OM and organic clay fractions present [81]. The low RNA content in the deeper layers in our study might indicate the presence of relic DNA; however, cell numbers in these deeper layers (>1.5 × 109 cells ml −1) were similar to those in upper layers (1.9–2.4 × 109 cells ml −1). The high DNA:RNA ratio may therefore result from an abundance of dormant and potentially dying cells [54], rather than extracellular relic DNA. In our lake, the extractable DNA seems to be rapidly decaying (similar to the sequences from eukaryotes as potential progenitors of relic DNA), which points to a short-lived fate of relic DNA in Lake Stechlin. We also found no evidence for fragmented DNA in our sediment samples that would indicate the presence of larger quantities of extracellular DNA (Additional file 8). Relic DNA would have caused an overestimate of the richness, in particular, the deeper layers would appear more species rich than they are. The use of a single SSU primer pair in our study was also a compromise and underestimates the richness of metazoan groups [37, 82], which affects the zooplankton OTUs that we found in the upper layers. This means that the richness component of β-diversity may be more important as a result, leading to a narrower transition zone and greater differences between the horizons. On the other hand, our study was also limited by resolution, since our pyrosequencing efforts could not adequately analyze the rare biosphere, in which we would suspect most signals from relic DNA. Future studies with higher sequencing depths and several group-specific primer pairs will be able to follow the fate of relic DNA in more details as, e.g., of Eukaryota.


Our results indicate that the sediments of Lake Stechlin are a steady-state and chemostat-like environment with a highly stratified indigenous microbial community. Sediments were not a one-way system for the burial of organic matter and were not composed purely of redox-active taxa. We conclude that both processes take place alongside a vertical gradient of electron acceptors and decomposing OM of decreasing quality with depth. The microbial community was structured into distinct groups, and both the microbial community and the sediment parameters could be divided into components relevant for burial and past conditions as well as for recent carbon turnover processes and their context data. Biotic interactions are likely to play an important role, and we were able to identify important sediment taxa for each horizon. We put a spotlight on the largely unexplored freshwater sediments and confirmed earlier findings that were previously described only for marine sediments, such as the importance of marine archaeal lineages and the introduction of a depauperation zone in which the burial process becomes increasingly important.


Sampling site and sampling procedures

Lake Stechlin (latitude 53° 10 N, longitude 13° 02 E) is a dimictic oligo-mesotrophic lake (maximum depth 69.5 m; area 4.23 km2) in northern Germany that has been the subject of more than 55 years of research [83]. Sediment cores were extracted from the southern bay of the lake. Four adjacent sites were sampled to account for spatial heterogeneity in the sediment (sites A–D, Additional file 9). Pore water was collected by four in situ dialysis samplers, so-called peepers [84], which were deployed for 14 days using a frame (1 m2). Shortly before retrieving the peepers, four sediment cores were taken from each site with Perspex tubes (inner diameter 6 or 9 cm; length 60 cm) using gravity corers (UWITEC™, Mondsee, Austria) at 30-m water depth (aphotic depth) on two subsequent days (March 26 and 28, 2012; peepers were retrieved on April 1, 2012). Two sediment cores (6-cm diameter) were stored in the dark at 4 °C until oxygen penetration depth was measured within the next 4 h (see below). The sediments of the 9-cm cores were sliced directly into 1-cm layers for the uppermost 10 cm and then in 4-cm layers for sediment depths of 10–30 cm. One core was used for the analysis of the total sediment, and the other was used for pore water, gas, and microbial analyses (see below). In May 2014, 24 additional cores were taken to determine the age-depth correlation using the cesium 137 technique.

Maximum oxygen penetration depth

Two initial cores were carefully transferred into 20-cm short cores without disturbing the sediment surface. The short cores were kept cool (4 °C) until measurements were taken. Oxygen microprofiles were performed using two Clark-type microelectrodes (OX50 oxygen microsensors, Unisense, Aarhus, Denmark) with a 50- μm glass tip. SensorTracePro 2.0 software (Unisense) was used for data storage. The electrodes were calibrated by two-point calibration. For each core, we measured at least four profiles. The sediment-water interface was defined as the point where the oxygen depletion shifted from linear to non-linear [85].

Pore water analysis

The sampled sediment horizons were centrifuged (13,250gfor 10 min) to retrieve pore water (filtered through rinsed 0.45- μm cellulose acetate membranes, Roth, Germany) for immediate analysis of the dissolved organic carbon (DOC) and FI. DOC was measured as non-purgeable organic carbon with an organic carbon analyzer (multi N/C 3100, Analytic Jena AG, Jena, Germany). FI was measured following the protocol of [38]. Peeper samples were analyzed for concentrations of SRP and ammonium (NH\(_{4}^{+}\)), dissolved iron (Fe 2+/3+), manganese (Mn 2+), chloride (Cl ), nitrate (NO\(_{3}^{-}\)), and sulfate (SO\(_{4}^{2-}\)), following DIN EN ISO 10304-1. SRP and NH\(_{4}^{+}\) were photometrically determined using segmented flow analysis (SFA, Skalar Sanplus, Skalar Analytical B.V., De Breda, Netherlands). Dissolved iron and manganese levels were determined by AAS (PerkinElmer 3300, Rodgau-Juegesheim, Germany), and analyses of the dissolved anions nitrate and sulfate were conducted by ion chromatography (IC, Shimadzu Corporation, Japan).

Total sediment analysis

Sediment water content was analyzed by drying at 85 °C until mass was constant. A subsample was used to determine the organic matter content (4 h at 550 °C) of the sediment. The metal concentrations were determined by ICP-OES (iCAP 6000, Thermo Fisher Scientific, Dreieich, Germany) after aqua regia digestion in a microwave oven (Gigatherm, Grub, Switzerland), and total phosphorus (TP) was determined spectrophotometrically by CARY 1E (Varian Deutschland GmbH, Darmstadt, Germany) after H2SO4/H2O2 digestion (150 C, 16 h). CNHS content was determined using aliquots of dried matter in a vario EL system (Elementar Analysensysteme GmbH, Hanau, Germany).

Gas chromatography

From each depth, 2 ml of sediment was transferred into 10-ml vials filled with 4 ml of distilled water. Samples were fixed with mercury chloride (final conc. 200 mg l −1), sealed, and stored in the dark at 4 °C until analysis. Concentrations of CO2, CH4, and N2O were measured by gas chromatography (Shimadzu GC-14B, Kyoto, Japan) using the headspace technique described in [86].

Bacterial protein production

Bacterial biomass production was determined via 14C leucine incorporation at in situ temperature under anoxic conditions [87] using a modified protocol [88]. Five hundred microliters of sediment was diluted 1:1 with sterile filtered supernatant water and incubated with 14C-leucine (Hartmann Analytics, Braunschweig, Germany; specific activity 306 mCi mmol−1, diluted with cold L-leucine to a final concentration of 50 μmol l−1). Incubations were stopped after 1 h, extracted, and measured in a liquid scintillation analyzer (TriCarb 2810 TR, PerkinElmer Inc., Germany). Disintegrations per minute were converted to pmol leucine ml −1 day −1, assuming a twofold intracellular isotope dilution [89, 90].

Cell counting

Sediment subsamples for cell counting were immediately fixed with ethanol (50% v/vfinal concentration). Prior to analysis, samples were shaken for 1 h at 700 rpm on a thermoshaker and were sonicated three times for 20 s at 5–6 W (Branson Sonifier 150, Danbury, USA). Cells were stained with a SYBR Gold staining solution diluted to 1:1000 (Molecular Probes, Eugene, USA) and were counted with an epifluorescence microscope (Zeiss, Axio Imager. Z1, Jena, Germany).

Nucleic acid extraction and sequencing

To determine the DNA:RNA ratio (as part of the present parameters, Additional file 10), we extracted total nucleic acids using a phenol-chloroform protocol from 200–400 μl sediment, as described by [91]. The DNA:RNA ratio was measured via fluorometry using selectively binding dyes (broad range dsDNA and broad range RNA assay Kit, Life technologies, Darmstadt, Germany) developed for the Qubit 2.0 (Life technologies, Darmstadt, Germany). A second extraction served as template for the sequencing and determination of the total DNA content. A defined sediment subsample (350 μl) from each depth was lyophilized prior to DNA extraction. We used the “Alternative Protocol for Maximum DNA yields” of the UltraClean®; Soil DNA Isolation Kit (MoBio Laboratories Inc., Carlsbad, USA). The quality of the DNA and the presence of putative environmental (small) DNA in 12 representative samples from 12 depths were verified with a microgel electrophoresis system (DNA High Sensitivity Kit, Bioanalyzer, Agilent, USA, see Additional file 8). A total of 5–20 ng of DNA, as measured by NanoPhotometer P300 (Implen, Schatzbogen, Germany), served as the template for PCR amplification (Herculase II system, Life Technologies) using a single universal primer system (926F, 1392R, [92]) targeting the SSU V6-V8 region. The primer pair employed is one important feature of our study in that it detects similar to the primer of [93] all three microbial domains (Archaea, Bacteria, and Eukaryota) in freshwater systems [37, 62, 94]. The variability of the V6–V8 is sufficiently high for all three domains [9597]. PCR products were purified with AMPure Beads (Beckmann Coulter, Brea, USA) and quantified and pooled using a PicoGreen assay (Life Technologies, Carlsbad, USA). High-throughput sequencing was performed in a Roche 454 GS Junior benchtop sequencer (Hoffmann-La Roche, Basel, Switzerland) at the Berlin Center for Genomics in Biodiversity Research.

Data processing

Four hundred fifty-four sequencing data were processed using Mothur (version 1.33.0) following the guidelines of the Mothur SOP (, accessed 02/2014) with the following modifications: (i) for quality trimming, we used a sliding window approach with a relaxed threshold (window size 50, quality cutoff 27), and (ii) the alignment step used SINA (version 1.2.11; [98]) against the SILVA v.111 non-redundant SSU reference database. A total of 396,000 reads (49% of 802,202 raw sequences) were retained. The sequences were clustered into OTUs at 97% sequence similarity (Additional file 10). A representative sequence from each OTU was used for taxonomic classification with the least common ancestor method in SINA, using 0.7 as a setting for minimum similarity as well as for lca-quorum. The classified OTU abundance matrix served as the basis for all subsequent statistical analyses (Additional file 11). The percentage of sequences with low similarity (< 93%) to the next reference sequence was determined by submitting the FASTA files to SILVA NGS [99].


All measured environmental parameters were compiled in a matrix and imported into R (, version 3.2.2; Additional file 12). We replaced two outliers (FI: replicate 4 cm, total phosphorous: replicate 14 cm) with the mean values of the three other sediment cores. Similarly, the 30 cm peeper data from replicate core B were missing and replaced by the mean of the residual replicates. For statistical analysis, the relative proportions of Archaea, Bacteria, and Eukaryota were arc-sin transformed. For the multiple regression analysis on the declining DNA concentrations, we removed one value (replicate 3, 22 cm) to meet the normal distribution criteria of the residuals. Sufficient normal distribution was confirmed by a QQ plot and Shapiro-Wilks test, p=0.183; Cook’s distance was not violated in any case. We categorized the environmental parameters into present (CH4, CO2, DOC, BPP, SRP, NH\(_{4}^{+}\), SO\(_{4}^{2-}\), Cl , Fe 2+/3+, Mn 2+, FI, and RNA:DNA) and past (TC, TN, dry-weight, TP, TS, TH, Al, As, Ca, Cu, Fe, Mg, Mn, Pb, Ti, and Zn) (see Table 1) and assessed the sample variation for each subset by a centered, scaled principal component analysis (PCA, see Additional file 10). The resulting most explanatory PCA axes, which explained 54% (present parameters) and 41% (past parameters) of the sample variation, were used in the community statistics (see below). DNA and cell numbers were not categorized due to their ambiguous nature; Al was used as a substitute for Mg and Ti due to their high degree of correlation (r>0.95); N2O, NO\(_{3}^{-}\), Cd, and Co were excluded due to their very low values, i.e., near or below the detection limit in all the samples.

Community statistics

Random effects were initially excluded from the OTU matrix by removing all OTUs present in only one sample, regardless of the number of reads. The random effects are expected to be very high in sediments due to the burial of random organic matter (e.g., caused by bird droppings, tourist activities, rainfall), and thus, a large number of rare taxa are expected. This reduced the number of OTUs from 29,228 to 9581 but did not influence the sample distances (Mantel test with Hellinger distances: r=0.993,p<0.001). Diversity indices were calculated using the vegan package [100] for R, with a community matrix that was rarefied to the lowest number of reads (642) present in a sample. The rarefied matrix was highly correlated to the initial matrix (Mantel test with Hellinger distances: r=0.923,p<0.001). Nonmetric multidimensional scaling (NMDS) and Mantel tests were calculated based on Hellinger-transformed rarefied OTU matrix with Euclidean distances. The PCA scores from the past and present parameters (see Table 1) were fitted into the NMDS, and their correlation with the underlying distance matrix was tested with a Mantel test. Additionally, we calculated the weighted UniFrac distances with the R package GUniFrac [101] and corresponding phylogenetic distances were based on a maximum likelihood tree calculated with FastTree 2.1 [102]. The weighted UniFrac distances were based on proportional data from the random effects reduced community matrix (9581 OTUs) and were projected as NMDS or MDS. We include the same statistics as in Fig. 4 (Additional file 6). Moreover, we used a fuzzy set ordination [103] to test for the influence of the past and present parameters on the separated richness and replacement community components. For this, we partitioned β-diversity into richness and replacement components using indices from the Jaccard family, following [104] and the functions provided by [32]. In order to identify general vertical patterns, we used a sum table to increase the resolution and thus avoid an artificial increase in turnover versus the richness/nestedness structure due to sampling effects. The sum table was generated by summing up 2000 sequences per depth, if applicable. The final sum table was rarefied to the lowest number of sequences in the depth profile (5987 reads). Hill numbers and rarefaction curves of the Hill numbers were calculated with the iNEXT package [105]. The cluster analysis (UPGMA clustering based on Kulczynski distance, Fig. 1) was also done on the sum table. The depth-dependent nestedness [31], richness component, replacement component, species contribution to beta diversity (SCBD), and LCBD were calculated as described in [32] including significance tests. We note that the nestedness index [31] is dependent on the sample size, and so we refer to it as “relative nestedness”.


  1. 1

    Tranvik LJ, Downing JA, Cotner JB, Loiselle SA, Striegl RG, Ballatore TJ, Dillon P, Finlay K, Fortino K, Knoll LB. Lakes and reservoirs as regulators of carbon cycling and climate. Limnol Oceanography. 2009; 54(6):2298–314. doi:10.4319/lo.2009.54.6_part_2.2298.

    CAS  Article  Google Scholar 

  2. 2

    Dean WE, Gorham E. Magnitude and signifigance of carbon burial in lakes, reservoirs, and peatlands. Geology. 1998; 26(6):535–8. doi:10.1130/0091-7613(1998)026.

    Article  Google Scholar 

  3. 3

    Bastviken D, Tranvik LJ, Downing JA, Crill PM, Enrich-Prast A. Freshwater methane emissions offset the continental carbon sink. Science. 2011; 331(6013):50. doi:10.1126/science.1196808. 3533

    CAS  PubMed  Article  Google Scholar 

  4. 4

    Fischer H, Sachse A, Steinberg CEW, Pusch M. Differential retention and utilization of dissolved organic carbon by bacteria in river sediments. Limnol Oceanography. 2002; 47(6):1702–11. doi:10.4319/lo.2002.47.6.1702. 3184

    Article  Google Scholar 

  5. 5

    Haglund AL, Lantz P, Törnblom E, Tranvik L. Depth distribution of active bacteria and bacterial activity in lake sediment. FEMS Microbiol Ecol. 2003; 46(1):31–8. doi:10.1016/S0168-6496(03)00190-9.

    CAS  PubMed  Article  Google Scholar 

  6. 6

    Nealson KH. Sediment bacteria: who’s there, what are they doing, and what’s new?Annu Rev Earth Planet Sci. 1997; 25(1):403–34. doi:10.1146/

    CAS  PubMed  Article  Google Scholar 

  7. 7

    Megonigal JP, Hines ME, Visscher PT. Anaerobic metabolism: linkages to trace gases and aerobic processes. In: Biogeochemistry. Oxford: Elsevier: 2004. p. 317–424.

    Google Scholar 

  8. 8

    Maerki M, Muller B, Dinkel C, Wehrli B. Mineralization pathways in lake sediments with different oxygen and organic carbon supply. Limnol Oceanography. 2009; 54(2):428.

    CAS  Article  Google Scholar 

  9. 9

    Melton ED, Swanner ED, Behrens S, Schmidt C, Kappler A. The interplay of microbially mediated and abiotic reactions in the biogeochemical Fe cycle. Nat Rev Microbiol. 2014; 12(12):797–808.

    CAS  PubMed  Article  Google Scholar 

  10. 10

    Frindte K, Allgaier M, Grossart HP, Eckert W. Microbial response to experimentally controlled redox transitions at the sediment water interface. PLOS ONE. 2015; 10(11):0143428. doi:10.1371/journal.pone.0143428.

    Article  CAS  Google Scholar 

  11. 11

    Lloyd KG, Schreiber L, Petersen DG, Kjeldsen KU, Lever MA, Steen AD, Stepanauskas R, Richter M, Kleindienst S, Lenk S, Schramm A, Jørgensen BB. Predominant archaea in marine sediments degrade detrital proteins. Nature. 2013; 496(7444):215–8. doi:10.1038/nature12033.

    CAS  PubMed  Article  Google Scholar 

  12. 12

    Miskin I, Rhodes G, Lawlor K, Saunders JR, Pickup RW. Bacteria in post-glacial freshwater sediments. Microbiology. 1998; 144(9):2427–39.

    CAS  PubMed  Article  Google Scholar 

  13. 13

    Li S, Xiao X, Yin X, Wang F. Bacterial community along a historic lake sediment core of Ardley Island, west Antarctica. Extremophiles. 2006; 10(5):461–7. doi:10.1007/s00792-006-0523-2.

    CAS  PubMed  Article  Google Scholar 

  14. 14

    Ye W, Liu X, Lin S, Tan J, Pan J, Li D, Yang H. The vertical distribution of bacterial and archaeal communities in the water and sediment of Lake Taihu. FEMS Microbiol Ecol. 2009; 70(2):263–76. doi:10.1111/j.1574-6941.2009.00761.x.

    CAS  Article  Google Scholar 

  15. 15

    Vuillemin A, Ariztegui D. Geomicrobiological investigations in subsaline maar lake sediments over the last 1500 years. Quaternary Sci Rev. 2013; 71:119–30. doi:10.1016/j.quascirev.2012.04.011.

    Article  Google Scholar 

  16. 16

    Xiong W, Xie P, Wang S, Niu Y, Yang X, Chen W. Sources of organic matter affect depth-related microbial community composition in sediments of Lake Erhai, Southwest China. J Limnol. 2015; 74(2):310–23. doi:10.4081/jlimnol.2014.1106.

    Google Scholar 

  17. 17

    Spring S, Schulze R, Overmann J, Schleifer KH. Identification and characterization of ecologically significant prokaryotes in the sediment of freshwater lakes: molecular and cultivation studies. FEMS Microbiol Rev. 2000; 24(5):573–90. doi:10.1016/S0168-6445(00)00046-2.

    CAS  PubMed  Article  Google Scholar 

  18. 18

    Briée C, Moreira D, López-García P. Archaeal and bacterial community composition of sediment and plankton from a suboxic freshwater pond. Res Microbiol. 2007; 158(3):213–27. doi:10.1016/j.resmic.2006.12.012.

    PubMed  Article  CAS  Google Scholar 

  19. 19

    Lipp JS, Morono Y, Inagaki F, Hinrichs KU. Significant contribution of Archaea to extant biomass in marine subsurface sediments. Nature. 2008; 454(7207):991–4.

    CAS  PubMed  Article  Google Scholar 

  20. 20

    Lloyd KG, May MK, Kevorkian RT, Steen AD. Meta-analysis of quantification methods shows that archaea and bacteria have similar abundances in the subseafloor. Appl Environ Microbiol. 2013; 79(24):7790–9. doi:10.1128/AEM.02090-13.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. 21

    Borrel G, Lehours AC, Crouzet O, Jézéquel D, Rockne K, Kulczak A, Duffaud E, Joblin K, Fonty G. Stratification of Archaea in the deep sediments of a freshwater meromictic lake: vertical shift from methanogenic to uncultured archaeal lineages. PLoS ONE. 2012; 7(8):8. doi:10.1371/journal.pone.0043346.

    Article  CAS  Google Scholar 

  22. 22

    Zhang J, Yang Y, Zhao L, Li Y, Xie S, Liu Y. Distribution of sediment bacterial and archaeal communities in plateau freshwater lakes. Appl Microbiol Biotechnol. 2015; 99(7):3291–302. doi:10.1007/s00253-014-6262-x.

    CAS  PubMed  Article  Google Scholar 

  23. 23

    Rodrigues T, Catão E, Bustamante MMC, Quirino BF, Kruger RH, Kyaw CM. Seasonal effects in a lake sediment archaeal community of the Brazilian Savanna. Archaea. 2014; 2014:957145. doi:10.1155/2014/957145. 957145

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  24. 24

    Gächter R, Meyer JS, Mares A. Contribution of bacteria to release and fixation of phosphorus in lake sediments. Limnol Oceanography. 1988; 33(6):1542–58. doi:10.4319/lo.1988.33.6part2.1542.

    Google Scholar 

  25. 25

    Chan OC, Claus P, Casper P, Ulrich A, Lueders T, Conrad R. Vertical distribution of structure and function of the methanogenic archaeal community in Lake Dagow sediment. Environ Microbiol. 2005; 7(8):1139–49. doi:10.1111/j.1462-2920.2005.00790.x.

    CAS  PubMed  Article  Google Scholar 

  26. 26

    Polymenakou PN, Fragkioudaki G, Tselepides A. Bacterial and organic matter distribution in the sediments of the Thracian Sea (NE Aegean Sea). Cont Shelf Res. 2007; 27(17):2187–97. doi:10.1016/j.csr.2007.05.003.

    Article  Google Scholar 

  27. 27

    Lomstein BA, Langerhuus AT, D’Hondt S, Jørgensen BB, Spivack AJ. Endospore abundance, microbial growth and necromass turnover in deep sub-seafloor sediment. Nature. 2012; 484(7392):101–4. doi:10.1038/nature10905.

    CAS  PubMed  Article  Google Scholar 

  28. 28

    Xie S, Lipp JS, Wegener G, Ferdelman TG, Hinrichs K. -u.Turnover of microbial lipids in the deep biosphere and growth of benthic archaeal populations,. Proc Nat Acad Sci USA. 2013; 110(15):6010–4. doi:10.1073/pnas.1218569110.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  29. 29

    Rothfuss F, Bender M, Conrad R. Survival and activity of bacteria in a deep, aged lake sediment (Lake Constance). Microbial Ecol. 1997; 33(1):69–77. doi:10.1007/s002489900009.

    CAS  Article  Google Scholar 

  30. 30

    Hoehler TM, Jørgensen BB. Microbial life under extreme energy limitation. Nat Rev Microbiol. 2013; 11(2):83–94. doi:10.1038/nrmicro2939.

    CAS  PubMed  Article  Google Scholar 

  31. 31

    Baselga A. Partitioning the turnover and nestedness components of beta diversity. Global Ecol Biogeography. 2010; 19(1):134–43. doi:10.1111/j.1466-8238.2009.00490.x.

    Article  Google Scholar 

  32. 32

    Legendre P. Interpreting the replacement and richness difference components of beta diversity. Global Ecol Biogeography. 2014; 23(11):1324–34. doi:10.1111/geb.12207.

    Article  Google Scholar 

  33. 33

    Corinaldesi C, Beolchini F, Dell’Anno A. Damage and degradation rates of extracellular DNA in marine sediments: implications for the preservation of gene sequences. Mol Ecol. 2008; 17(17):3939–51. doi:10.1111/j.1365-294X.2008.03880.x.

    CAS  PubMed  Article  Google Scholar 

  34. 34

    Wunderlin T, Corella JP, Junier T, Bueche M, Loizeau JL, Girardclos S, Junier P. Endospore-forming bacteria as new proxies to assess impact of eutrophication in Lake Geneva (Switzerland-France). Aquat Sci. 2014; 76(S1):103–16. doi:10.1007/s00027-013-0329-0.

    CAS  Article  Google Scholar 

  35. 35

    Meyers PA, Ishiwatari R. Lacustrine organic geochemistry—an overview of indicators of organic-matter sources and diagenesis in lake-sediments. Org Geochem. 1993; 20(7):867–900. doi:10.1016/0146-6380(93)90100-p.

    CAS  Article  Google Scholar 

  36. 36

    Arndt S, Jørgensen BB, LaRowe DE, Middelburg JJ, Pancost RD, Regnier P. Quantifying the degradation of organic matter in marine sediments: a review and synthesis. Earth-science Rev. 2013; 123:53–86.

    CAS  Article  Google Scholar 

  37. 37

    Wurzbacher C, Attermeyer K, Kettner MT, Flintrop C, Warthmann N, Hilt S, Grossart HP, Monaghan MT. Metabarcoding of unfractionated water samples relates phyto-, zoo-and bacterioplankton dynamics and reveals a single-taxon bacterial bloom. bioRxiv. 2016:058628. doi:10.1101/058628.

  38. 38

    McKnight DM, Boyer EW, Westerhoff PK, Doran PT, Kulbe T, Andersen DT. Spectrofluorometric characterization of dissolved organic matter for indication of precursor organic material and aromaticity. Limnol Oceanography. 2001; 46(1):38–48. doi:10.4319/lo.2001.46.1.0038.

    CAS  Article  Google Scholar 

  39. 39

    Chao A, Gotelli NJ, Hsieh T, Sander EL, Ma K, Colwell RK, Ellison AM. Rarefaction and extrapolation with Hill numbers: a framework for sampling and estimation in species diversity studies. Ecol Monographs. 2014; 84(1):45–67. doi:10.1890/13-0133.1.

    Article  Google Scholar 

  40. 40

    Comte J, Lindström ES, Eiler A, Langenheder S. Can marine bacteria be recruited from freshwater sources and the air?ISME J. 2014; 8(12):1–8. doi:10.1038/ismej.2014.89.

    Article  Google Scholar 

  41. 41

    Capone DG, Kiene RP. Comparison of microbial dynamics in marine and freshwater sediments: contrasts in anaerobic carbon catabolism. Limnol Oceanography. 1988; 33(4):725–49.

    CAS  Google Scholar 

  42. 42

    Luna GM, Corinaldesi C, Rastelli E, Danovaro R. Patterns and drivers of bacterial alpha- and beta-diversity across vertical profiles from surface to subsurface sediments. Environ Microbiol Rep. 2013; 5(5):731–9. doi:10.1111/1758-2229.12075.

    PubMed  Google Scholar 

  43. 43

    Jørgensen SL, Hannisdal B, Lanzén A, Baumberger T, Flesland K, Fonseca R, Ovreås L, Steen IH, Thorseth IH, Pedersen RB, Schleper C. Correlating microbial community profiles with geochemical data in highly stratified sediments from the Arctic Mid-Ocean Ridge. Proc Nat Acad Sci USA. 2012; 109(42):2846–55. doi:10.1073/pnas.1207574109.

    Article  Google Scholar 

  44. 44

    Koizumi Y, Kojima H, Fukui M. Characterization of depth-related microbial community structure in lake sediment by denaturing gradient gel electrophoresis of amplified 16S rDNA and reversely transcribed 16S rRNA fragments. FEMS Microbiol Ecol. 2003; 46(2):147–57. doi:10.1016/S0168-6496(03)00212-5.

    CAS  PubMed  Article  Google Scholar 

  45. 45

    Sass H, Cypionka H, Babenzien HD. Vertical distribution of sulfate-reducing bacteria at the oxic-anoxic interface in sediments of the oligotrophic Lake Stechlin. FEMS Microbiol Ecol. 1997; 22(3):245–55.

    CAS  Article  Google Scholar 

  46. 46

    Deutzmann JS, Stief P, Brandes J, Schink B. Anaerobic methane oxidation coupled to denitrification is the dominant methane sink in a deep lake. Proc Nat Acad Sci USA. 2014; 111(51):18273–8. doi:10.1073/pnas.1411617111.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  47. 47

    Thevenon F, Graham ND, Chiaradia M, Arpagaus P, Wildi W, Poté J. Local to regional scale industrial heavy metal pollution recorded in sediments of large freshwater lakes in central Europe (lakes Geneva and Lucerne) over the last centuries. Sci Total Environ. 2011; 412–413:239–47. doi:10.1016/j.scitotenv.2011.09.025.

    PubMed  Article  CAS  Google Scholar 

  48. 48

    Treude T, Krause S, Maltby J, Dale AW, Coffin R, Hamdan LJ. Sulfate reduction and methane oxidation activity below the sulfate-methane transition zone in Alaskan Beaufort Sea continental margin sediments : implications for deep sulfur cycling. Geochimica et Cosmochimica Acta. 2014; 144:217–37. doi:10.1016/j.gca.2014.08.018.

    CAS  Article  Google Scholar 

  49. 49

    Conrad R, Chan OC, Claus P, Casper P. Characterization of methanogenic Archaea and stable isotope fractionation during methane production in the profundal sediment of an oligotrophic lake (Lake Stechlin, Germany). Limnol Oceanography. 2007; 52(4):1393–1406. doi:10.4319/lo.2007.52.4.1393.

    CAS  Article  Google Scholar 

  50. 50

    Lovley DR, Coates JD, Blunt-Harris EL, Phillips EJP, Woodward JC. Humic substances as electron acceptors for microbial respiration. Nature. 1996; 382(6590):445–8. doi:10.1038/382445a0.

    CAS  Article  Google Scholar 

  51. 51

    Torres IC, Inglett KS, Reddy KR. Heterotrophic microbial activity in lake sediments: effects of organic electron donors. Biogeochemistry. 2011; 104(1–3):165–81. doi:10.1007/s10533-010-9494-6.

    CAS  Article  Google Scholar 

  52. 52

    Glombitza C, Stockhecke M, Schubert CJ, Vetter A, Kallmeyer J. Sulfate reduction controlled by organic matter availability in deep sediment cores from the saline, alkaline Lake Van (Eastern Anatolia, Turkey). Front Microbiol. 2013; 4(JUL):209. doi:10.3389/fmicb.2013.00209.

    PubMed  PubMed Central  Google Scholar 

  53. 53

    Parkes RJ, Cragg B, Roussel E, Webster G, Weightman A, Sass H. A review of prokaryotic populations and processes in sub-seafloor sediments, including biosphere: geosphere interactions. Marine Geol. 2014; 352:409–25. doi:10.1016/j.margeo.2014.02.009.

    CAS  Article  Google Scholar 

  54. 54

    Luna GM, Manini E, Danovaro R. Large fraction of dead and inactive bacteria in coastal marine sediments: comparison of protocols for determination and ecological significance. Appl Environ Microbiol. 2002; 68(7):3509–13. doi:10.1128/AEM.68.7.3509-3513.2002.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  55. 55

    Danovaro R, Dell’Anno A, Corinaldesi C, Magagnini M, Noble R, Tamburini C, Weinbauer M. Major viral impact on the functioning of benthic deep-sea ecosystems,. Nature. 2008; 454(7208):1084–7. doi:10.1038/nature07268.

    CAS  PubMed  Article  Google Scholar 

  56. 56

    Engelhardt T, Kallmeyer J, Cypionka H, Engelen B. High virus-to-cell ratios indicate ongoing production of viruses in deep subsurface sediments. ISME J. 2014; 8(7):1–7. doi:10.1038/ismej.2013.245.

    Article  Google Scholar 

  57. 57

    Dell’Anno A, Corinaldesi C, Danovaro R. Virus decomposition provides an important contribution to benthic deep-sea ecosystem functioning. Proc Nat Acad Sci. 2015; 112(16):2014–9. doi:10.1073/pnas.1422234112.

    Article  CAS  Google Scholar 

  58. 58

    Kandel PP, Pasternak Z, van Rijn J, Nahum O, Jurkevitch E. Abundance, diversity and seasonal dynamics of predatory bacteria in aquaculture zero discharge systems. FEMS Microbiol Ecol. 2014; 89(1):149–61. doi:10.1111/1574-6941.12342.

    CAS  PubMed  Article  Google Scholar 

  59. 59

    Fuchs A, Selmeczy GB, Kasprzak P, Padisák J, Casper P. Coincidence of sedimentation peaks with diatom blooms, wind, and calcite precipitation measured in high resolution by a multi-trap. Hydrobiologia. 2016; 763(1):329–44. doi:10.1007/s10750-015-2388-9.

    Article  Google Scholar 

  60. 60

    Domaizon I, Savichtcheva O, Debroas D, Arnaud F, Villar C, Pignol C, Alric B, Perga ME. DNA from lake sediments reveals the long-term dynamics and diversity of Synechococcus assemblages. Biogeosci Discuss. 2013; 10(2):2515–64. doi:10.5194/bgd-10-2515-2013.

    Article  Google Scholar 

  61. 61

    Capo E, Debroas D, Arnaud F, Domaizon I. Is planktonic diversity well recorded in sedimentary DNA? Toward the reconstruction of past protistan diversity. Microb Ecol. 2015; 70(4):865–75. doi:10.1007/s00248-015-0627-2.

    CAS  PubMed  Article  Google Scholar 

  62. 62

    Gies EA, Konwar KM, Thomas Beatty J, Hallam SJ. Illuminating microbial dark matter in meromictic Sakinaw Lake. Appl Environ Microbiol. 2014; 80(21):6807–18. doi:10.1128/AEM.01774-14.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  63. 63

    Adrian L, Szewzyk U, Wecke J, Görisch H. Bacterial dehalorespiration with chlorinated benzenes. Nature. 2000; 408(6812):580–3. doi:10.1038/35046063.

    CAS  PubMed  Article  Google Scholar 

  64. 64

    Ferrer M, Guazzaroni ME, Richter M, García-Salamanca A, Yarza P, Suárez-Suárez A, Solano J, Alcaide M, van Dillewijn P, Molina-Henares MA, López-Cortés N, Al-Ramahi Y, Guerrero C, Acosta A, de Eugenio LI, Martínez V, Marques S, Rojo F, Santero E, Genilloud O, Pérez-Pérez J, Rosselló-Móra R, Ramos JL. Taxonomic and functional metagenomic profiling of the microbial community in the anoxic sediment of a sub-saline shallow lake (Laguna de Carrizo, Central Spain). Microb Ecol. 2011; 62(4):824–37. doi:10.1007/s00248-011-9903-y.

    CAS  PubMed  Article  Google Scholar 

  65. 65

    Hiraishi A. Biodiversity of dehalorespiring bacteria with special emphasis on polychlorinated biphenyl/dioxin dechlorinators. Microbes Environ. 2008; 23(1):1–12. doi:10.1264/jsme2.23.1.

    PubMed  Article  Google Scholar 

  66. 66

    Hug LA, Castelle CJ, Wrighton KC, Thomas BC, Sharon I, Frischkorn KR, Williams KH, Tringe SG, Banfield JF. Community genomic analyses constrain the distribution of metabolic traits across the Chloroflexi phylum and indicate roles in sediment carbon cycling. Microbiome. 2013; 1(1):22. doi:10.1186/2049-2618-1-22.

    PubMed  PubMed Central  Article  Google Scholar 

  67. 67

    Sørensen KB, Teske A. Stratified communities of active Archaea in deep marine subsurface sediments stratified communities of active Archaea in deep marine subsurface sediments. Appl Environ Microbiol. 2006; 72(7):4596–603. doi:10.1128/AEM.00562-06.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  68. 68

    Evans PN, Parks DH, Chadwick GL, Robbins SJ, Orphan VJ, Golding SD, Tyson GW. Methane metabolism in the archaeal phylum Bathyarchaeota revealed by genome-centric metagenomics. Science. 2015; 350(6259):434–8. doi:10.1126/science.aac7745.

    CAS  PubMed  Article  Google Scholar 

  69. 69

    Kadnikov VV, Mardanov AV, Beletsky AV, Shubenkova OV, Pogodaeva TV, Zemskaya TI, Ravin NV, Skryabin KG. Microbial community structure in methane hydrate-bearing sediments of freshwater Lake Baikal. FEMS Microbiol Ecol. 2012; 79(2):348–58. doi:10.1111/j.1574-6941.2011.01221.x.

    CAS  PubMed  Article  Google Scholar 

  70. 70

    Baker BJ, Comolli LR, Dick GJ, Hauser LJ, Hyatt D, Dill BD, Land ML, VerBerkmoes NC, Hettich RL, Banfield JF. Enigmatic, ultrasmall, uncultivated Archaea. Proc Nat Acad Sci. 2010; 107(19):8806–11. doi:10.1073/pnas.0914470107.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  71. 71

    Koizumi Y, Takii S, Nishino M, Nakajima T. Vertical distributions of sulfate-reducing bacteria and methane-producing archaea quantified by oligonucleotide probe hybridization in the profundal sediment of a mesotrophic lake. FEMS Microbiol Ecol. 2003; 44(1):101–8. doi:10.1016/S0168-6496(02)00463-4.

    CAS  PubMed  Article  Google Scholar 

  72. 72

    Schwarz JIK, Eckert W, Conrad R. Community structure of Archaea and Bacteria in a profundal lake sediment Lake Kinneret (Israel). Syst Appl Microbiol. 2007; 30(3):239–54. doi:10.1016/j.syapm.2006.05.004.

    CAS  PubMed  Article  Google Scholar 

  73. 73

    Koonin EV, Wolf YI. Genomics of bacteria and archaea: the emerging dynamic view of the prokaryotic world. Nucleic Acids Res. 2008; 36(21):6688–719. doi:10.1093/nar/gkn668.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  74. 74

    Angly FE, Dennis PG, Skarshewski A, Vanwonterghem I, Hugenholtz P, Tyson GW. CopyRighter: a rapid tool for improving the accuracy of microbial community profiles through lineage-specific gene copy number correction. Microbiome. 2014; 2(1):11. doi:10.1186/2049-2618-2-11.

    PubMed  PubMed Central  Article  Google Scholar 

  75. 75

    Spang A, Saw JH, Jørgensen SL, Zaremba-Niedzwiedzka K, Martijn J, Lind AE, van Eijk R, Schleper C, Guy L, Ettema TJ. Complex archaea that bridge the gap between prokaryotes and eukaryotes. Nature. 2015; 521(7551):173–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  76. 76

    Meng J, Xu J, Qin D, He Y, Xiao X, Wang F. Genetic and functional properties of uncultivated MCG archaea assessed by metagenome and gene expression analyses. ISME J. 2014; 8(3):650–9. doi:10.1038/ismej.2013.174.

    CAS  PubMed  Article  Google Scholar 

  77. 77

    Kubo K, Lloyd KG, F Biddle J, Amann R, Teske A, Knittel K. Archaea of the miscellaneous crenarchaeotal group are abundant, diverse and widespread in marine sediments. ISME J. 2012; 6(10):1949–65. doi:10.1038/ismej.2012.37.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  78. 78

    Carini P, Marsden PJ, Leff JW, Morgan EE, Strickland MS, Fierer N. Relic DNA is abundant in soil and obscures estimates of soil microbial diversity. Nat Microbiol. 2016; 2:16242. doi:10.1038/nmicrobiol.2016.242.

    PubMed  Article  CAS  Google Scholar 

  79. 79

    Corinaldesi C, Danovaro R, Dell’Anno A. Simultaneous recovery of extracellular and intracellular DNA suitable for molecular studies from marine sediments. Appl Environ Microbiol. 2005; 71(1):46–50. doi:10.1128/AEM.71.1.46-50.2005.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  80. 80

    Ogram A, Sayler GS, Barkay T. The extraction and purification of microbial DNA from sediments. J Microbiol Methods. 1987; 7(2–3):57–66. doi:10.1016/0167-7012(87)90025-X.

    CAS  Article  Google Scholar 

  81. 81

    Cai P, Huang QY, Zhang XW. Interactions of DNA with clay minerals and soil colloidal particles and protection against degradation by DNase. Environ Sci Technol. 2006; 40(9):2971–6. doi:10.1021/es0522985.

    CAS  PubMed  Article  Google Scholar 

  82. 82

    Tang CQ, Leasi F, Obertegger U, Kieneke A, Barraclough TG, Fontaneto D. The widely used small subunit 18S rDNA molecule greatly underestimates true diversity in biodiversity surveys of the meiofauna. Proc Nat Acad Sci. 2012; 109(40):16208–12. doi:10.1073/pnas.1209160109.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  83. 83

    Koschel R, Adams D. An approach to understanding a temperate oligotrophic lowland lake (Lake Stechlin, Germany). Archiv fur Hydrobiologie, Special Issues Adv Limnol. 2003; 58:1–9.

    Google Scholar 

  84. 84

    Hesslein RH. An in situ sampler for close interval pore water studies. Limnol Oceanography. 1976; 21(6):912–4. doi:10.4319/lo.1976.21.6.0912.

    CAS  Article  Google Scholar 

  85. 85

    Glud RN. Oxygen dynamics of marine sediments. Marine Biol Res. 2008; 4(4):243–89. doi:10.1080/17451000801888726.

    Article  Google Scholar 

  86. 86

    Casper P, Albino MF, Adams DD. Diffusive fluxes of CH4 and CO2 across the water-air interface in the eutrophic Lake Dagow, northeast Germany. Verh. Int. Verein. Limnol. 2009; 30(6):874–77.

    CAS  Google Scholar 

  87. 87

    Buesing N, Gessner MO. Incorporation of radiolabeled leucine into protein to estimate bacterial production in plant litter, sediment, epiphytic biofilms, and water samples. Microb Ecol. 2003; 45(3):291–301. doi:10.1007/s00248-002-2036-6. 3319

    CAS  PubMed  Article  Google Scholar 

  88. 88

    Attermeyer K, Premke K, Hornick T, Hilt S, Grossart HP. Ecosystem-level studies of terrestrial carbon reveal contrasting bacterial metabolism in different aquatic habitats. Ecology. 2013; 94(12):2754–66.

    PubMed  Article  Google Scholar 

  89. 89

    Azam F, Simon M. Protein content and protein synthesis rates of planktonic marine bacteria. Marine Ecol Prog Ser. 1989; 51(3):213. doi:10.3354/meps051201.

    Google Scholar 

  90. 90

    Kirchman DL. In: Kemp PF, Sherr BF, Sherr EB, Cole JJ, (eds).Leucine incorporation as a measure of biomass production by heterotrophic bacteria. Handbook of methods in aquatic microbial ecology. Boca Raton: Lewis Pub.; 1993. p. 509–12.

    Google Scholar 

  91. 91

    Nercessian O, Noyes E, Kalyuzhnaya MG, Lidstrom ME, Chistoserdova L. Bacterial populations active in metabolism of C 1 compounds in the sediment of Lake Washington, a freshwater lake. Appl Environ Microbiol. 2005; 71(11):6885–99. doi:10.1128/AEM.71.11.6885-6899.2005.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  92. 92

    Engelbrektson A, Kunin V, Wrighton KC, Zvenigorodsky N, Chen F, Ochman H, Hugenholtz P. Experimental factors affecting PCR-based estimates of microbial species richness and evenness. ISME J. 2010; 4(5):642–7. doi:10.1038/ismej.2009.153.

    CAS  PubMed  Article  Google Scholar 

  93. 93

    Parada AE, Needham DM, Fuhrman JA. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ Microbiol. 2016; 18(5):1403–14. doi:10.1111/1462-2920.13023.

    CAS  PubMed  Article  Google Scholar 

  94. 94

    Hölker F, Wurzbacher C, Weißenborn C, Monaghan MT, Holzhauer SIJ, Premke K. Microbial diversity and community respiration in freshwater sediments influenced by artificial light at night. Phil Trans R Soc London. Series B, Biol Sci. 2015; 370:20140130. doi:10.1098/rstb.2014.0130.

    Article  CAS  Google Scholar 

  95. 95

    Vasileiadis S, Puglisi E, Arena M, Cappa F, Cocconcelli PS, Trevisan M. Soil bacterial diversity screening using single 16S rRNA gene V regions coupled with multi-million read generating sequencing technologies. PloS one. 2012; 7(8):42671.

    Article  CAS  Google Scholar 

  96. 96

    Seedorf H, Kittelmann S, Henderson G, Janssen PH. Rim-db: a taxonomic framework for community structure analysis of methanogenic archaea from the rumen and other intestinal environments. PeerJ. 2014; 2:494.

    Article  Google Scholar 

  97. 97

    Hadziavdic K, Lekang K, Lanzen A, Jonassen I, Thompson EM, Troedsson C. Characterization of the 18S rrna gene for designing universal eukaryote specific primers. PLoS One. 2014; 9(2):87624.

    Article  CAS  Google Scholar 

  98. 98

    Pruesse E, Peplies J, Glöckner FO. SINA: accurate high throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics. 2012; 28(14):1823–9. doi:10.1093/bioinformatics/bts252.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  99. 99

    Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glöckner F. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013; 41(D1):590–6. doi:10.1093/Nar/Gks1219.

    Article  CAS  Google Scholar 

  100. 100

    Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’Hara RB, Simpson GL, Solymos P, Stevens HH, Wagner H. vegan: Community Ecology Package. R package version 2.3-5. 2016.

  101. 101

    Chen J, Bittinger K, Charlson ES, Hoffmann C, Lewis J, Wu GD, Collman RG, Bushman FD, Li H. Associating microbiome composition with environmental covariates using generalized unifrac distances. Bioinformatics. 2012; 28(16):2106–13. doi:10.1093/bioinformatics/bts342.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  102. 102

    Price MN, Dehal PS, Arkin AP. Fasttree 2—approximately maximum-likelihood trees for large alignments. PloS one. 2010; 5(3):9490. doi:10.1371/journal.pone.0009490.

    Article  CAS  Google Scholar 

  103. 103

    Roberts DW. Statistical analysis of multidimensional fuzzy set ordinations. Ecology. 2008; 89(5):1246–60.

    PubMed  Article  Google Scholar 

  104. 104

    Carvalho JC, Cardoso P, Gomes P. Determining the relative roles of species replacement and species richness differences in generating beta-diversity patterns. Global Ecol Biogeography. 2012; 21(7):760–71. doi:10.1111/j.1466-8238.2011.00694.x.

    Article  Google Scholar 

  105. 105

    Hsieh T, Ma K, Chao A. iNEXT: an R package for rarefaction and extrapolation of species diversity (Hill numbers). Methods Ecol Evol. 2016; 7(12):1451–6. doi:10.1111/2041-210X.12613.

    Article  Google Scholar 

  106. 106

    Chiu CH, Wang YT, Walther BA, Chao A. An improved nonparametric lower bound of species richness via a modified good-turing frequency formula. Biometrics. 2014; 70(3):671–82.

    PubMed  Article  Google Scholar 

  107. 107

    Ondov BD, Bergman NH, Phillippy AM. Interactive metagenomic visualization in a web browser. BMC Bioinformatics. 2011; 12(1):1.

    Article  Google Scholar 

Download references


We thank Michael Sachtleben and Federica Pinto for their help during sampling. We also thank Elke Mach, Uta Mallok, Christiane Herzog, Hans-Jürgen Exner, Grit Siegert, Susan Mbedi, Kirsten Richter, Camila Mazzoni, and Andreas Kleeberg for their assistance and measurements in the laboratory and Henrik R. Nilsson and two anonymous referees for the comments on earlier versions of this manuscript. We thank U.-K. Schkade (Federal Office for Radiation Protection, Berlin, Germany) for the dating of the Lake Stechlin sediments.


This research was funded by the Leibniz Association Pakt for Research and Innovation projects “Tembi” (SAW-2011-IGB-2) and “MycoLink” (SAW-2014-IGB-1) and by grants DFG GR1540/15-1 and 23-1. The funding bodies had no influence on the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials

The sequence data sets generated during and analyzed during the current study are available in the ENA repository under the accession number PRJEB14189, All environmental data generated and analyzed during this study are included in this published article [and its supplementary information files].

Authors’ contributions

All co-authors designed the experimental setup. CW, AF, KA, and KF performed the experiments and measurements. CW and AF analyzed the data. All co-authors compiled the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

The manuscript does not contain any individual person’s data in any form.

Ethics approval and consent to participate

The manuscript does not report data collected from humans or animals.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information



Corresponding author

Correspondence to Christian Wurzbacher.

Additional files

Additional file 1

Krona chart of recovered sediment taxa. Browsable Krona chart ([107], S1.html, please use an internet browser with network access to open the file) of all taxa based on the median occurrence of OTUs for each depth replicate and classified against the SILVA reference database (, version 111). (HTML 114 kb)

Additional file 2

Figure: detailed depth profiles of individual cores. Additional detailed depth profiles of individual core (A,B,C,D) variables of Lake Stechlin at 30-m depth. Units: w.c.(water content) [%]; C [%]; N [%]; S [%];H [%]; Ca [mg g −1 dry weight]; Mg [mg g −1 dry weight]; NO\(_{3}^{-}\) [mg l −1]; SO\(_{4}^{2-}\) [mg l −1]; Fe 2+/3+ [mg l −1]; Mn 2+ [mg l −1]; Al [mg g −1 dry weight]; Cd [mg g −1 dry weight]; Co [mg g −1 dry weight]; Cr [mg g −1 dry weight]; Cu [mg g −1 dry weight]; Mn [mg g −1 dry weight]; Ni [mg g −1 dry weight]; Ti [mg g −1 dry weight]; Zn [mg g −1 dry weight]. See [108] for comparison with previous data. (PDF 91 kb)

Additional file 3

Figure: richness component vs. depth. Increasing richness component with increasing depth. The first cm is an outlier of the observed linearity. (PDF 15 kb)

Additional file 4

Figure: taxonomic composition of the most structuring taxa. Hierarchical taxonomic presentation of the most structuring taxa (SCBD), i.e., all OTUs that account for more than 5 per mill of the total β-diversity (see inlet to the left). The pie chart is color coded according to the three domains: Bacteria (red), Archaea (green), and Eukaryota (blue). (PDF 197 kb)

Additional file 5

Figure: sediment DNA as a function of present taxonomic signals. Multiple linear regression on the sediment DNA content as a function of the occurrence of Eukaryota (75.6% of the variation) together with Bacteria (10.0% of the variation; model: R 2=0.856, p<0.001). (PDF 10 kb)

Additional file 6

Figure: UniFrac ordinations. Left panel - A nonmetric multidemsional scaling (analogous to Fig. 4 b) of all the samples based weighted UniFrac distances. This was also reflected in the distance between the surface and deep sediments on axis 1 (adonis: R 2=0.520,p<0.001). We were able to significantly recover the three depth zones (adonis: R 2=0.601,p<0.001). The overall community structure was correlated with both present (Mantel correlation: r=0.512,p<0.001) and past (r=0.333,p<0.001) parameters, which were nearly orthogonal in ordination. Right panel—a metric multidimensional scaling (principal coordinate analysis) of the UniFrac distance matrix that is displayed in Fig. 3 b, with the corresponding proportional eigenvalues for each axis. The curved shape may point to an ordination artifact. (PDF 27 kb)

Additional file 7

Text: sediment zonation according to taxonomic clustering, β-partitioning, and context data. (PDF 14 kb)

Additional file 8

Figure: DNA size distribution after extraction. DNA microgel electrophoresis (Experion, BioRad) from a random subset of samples from various sediment depths, showing the absence of small environmental DNA (<1 kb). (PDF 395 kb)

Additional file 9

Figure: sampling locations within Lake Stechlin, Germany. Depth map of Lake Stechlin (Germany) and the four replicate sampling sites (A,B,C,D) in the South-West bay with the corresponding oxygen penetration depth in cm (pink bars). (PDF 389 kb)

Additional file 10

Figure: principal component analysis of environmental parameters. Principal component analysis defining the “present” (left panel) and “past” (right panel) parameters. The samples are color coded according to the three depth clusters (a–c). (PDF 52 kb)

Additional file 11

Dataset: classified OTU table. Tab-separated text file that compiles the complete OTU table with all samples, read counts, representative FASTA sequences, and classifications. (CSV 1920 kb)

Additional file 12

Dataset: environmental parameters. Tab-separated text file that compiles all environmental parameters. (CSV 1570 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wurzbacher, C., Fuchs, A., Attermeyer, K. et al. Shifts among Eukaryota, Bacteria, and Archaea define the vertical organization of a lake sediment. Microbiome 5, 41 (2017).

Download citation


  • Archaea
  • Eukaryota
  • Bacteria
  • Community
  • Freshwater
  • Lake
  • DNA metabarcoding
  • Beta-diversity
  • Sediment
  • Turnover