Vertical variations in microbial diversity, composition, and interactions in freshwater lake sediments on the Tibetan plateau

Microbial communities in freshwater lake sediments exhibit a distinct depth-dependent variability. Further exploration is required to understand their biodiversity pattern and microbial interactions in vertical sediments. In this study, sediment cores from two freshwater lakes, Mugecuo (MGC) and Cuopu (CP), on the Tibetan plateau were sampled and subsequently sliced into layers at a depth of every centimeter or half a centimeter. Amplicon sequencing was used to analyze the composition, diversity, and interaction of microbial communities. Results showed that sediment samples of both lakes could be clustered into two groups at a sediment depth of about 20 cm, with obvious shifts in microbial community compositions. In lake MGC, the richness component dominated β-diversity and increased with depth, indicating that the microbial communities in the deep layer of MGC was selected from the surface layer. Conversely, the replacement component dominated β-diversity in CP, implying a high turnover rate in the surface layer and inactive seed banks with a high variety in the deep layer. A co-occurrence network analysis showed that negative microbial interactions were prevalent in the surface layers with high nutrient concentrations, while positive microbial interactions were more common in the deep layers with low nutrient concentrations, suggesting that microbial interactions are influenced by nutrient conditions in the vertical sediments. Additionally, the results highlight the significant contributions of abundant and rare taxa to microbial interactions and vertical fluctuations of β-diversity, respectively. Overall, this work deepens our understanding of patterns of microbial interactions and vertical fluctuation in β-diversity in lake sediment columns, particularly in freshwater lake sediments from the Tibetan plateau.

Introduction

Lake sediment is an essential global organic carbon sink (Dean and Gorham, 1998), formed by the interaction of deposition and decomposition of organic matter (OM) in the overlying water (Ambrosetti et al., 2003). Like marine sediments, as the depth increases, the contribution of refractory organic matter also increases while the availability of electron acceptors decrease (Ambrosetti et al., 2003; Wurzbacher et al., 2017; D'Hondt et al., 2019), which may lead to changes in the microbial communities in the lake sediments (Rissanen et al., 2019). Environmental factors such as nutrients, salinity, and pH, which vary with depth, also impact the microbial community structure in sediment environments (Edlund et al., 2006; Nelson et al., 2007; Zhao et al., 2011). Microorganisms in sediment play a crucial role in OM mineralization and in nutrient cycling and energy transfer processes (Mason et al., 2009). Understanding the microbial communities present in vertical lake sediments enables us to investigate the depth-related changes in ecosystem processes.

The β-diversity, as the extent of change in community composition, was first proposed in 1960 (Whittaker, 1960). Previous studies have found that β-diversity could reflect nestedness and spatial turnover (Harrison et al., 1992; Baselga, 2007, 2010). Nestedness, also named richness difference, reflects the change in the relative abundance of ever-present taxa (Legendre, 2014). Spatial turnover, also called replacement, means the extent of replacement for different species without the influence of richness difference (Baselga, 2010). Partitioning β-diversity into nestedness and turnover components could reveal ecological processes for species composition variation in the sediment core. For example, a nested microbial community selected from the lake surface layer was found in the lake’s deep layer (Vuillemin et al., 2018), indicating that microbial community variation is mainly driven by the environment filter effect. However, a high replacement effect was also found to be dominant in the upper sediments of lakes, indicating the impact of cellular turnover and random appearances on microbial communities (Wurzbacher et al., 2017).

Microorganisms interact with each other in different ways. Some strains may cooperate and rely heavily on their allies, while others may compete with each other for basic finite resources like nutrients, light, or territory (Bauer et al., 2018). The nutrient concentration has been found to affect bacterial interactions in pure culture studies and in situ habitats (Khan et al., 2018; Velez et al., 2018; Bergk Pinto et al., 2019; Hester et al., 2019). Those studies suggested that positive interactions predominate in environments with limited nutrients, while negative interactions predominate in nutrient-rich environments. Previous research on marine and lake surface sediment has shown that negative microbial interactions increase with nutrient supplementation (Dai et al., 2022; Zhang et al., 2022). However, the extent to which microbial interactions change with depth in vertical lake sediment with varied nutrient concentrations remains unclear.

Microbial communities generally consist of a few abundant and many rare species (Logares et al., 2014). Abundant and rare species are defined based on their relative abundance with a specific cutoff level. For instance, species with a relative abundance higher than 0.1% are defined as abundant (Jiao and Lu, 2020). Besides, the frequency of species occurrence is also taken into consideration (Dai et al., 2022). Abundant and rare taxa show different characteristics and ecological roles (Pedros-Alio, 2012; Logares et al., 2014). Abundant taxa play a crucial role in biomass, carbon flow, and nutrient cycling, while rare taxa are essential for biodiversity and redundant functions (Pedros-Alio, 2012; Debroas et al., 2015). Previous research has revealed differences in the distribution patterns and functional traits of abundant and rare subcommunities in various environments (Zhang X. et al., 2021; Ren et al., 2022). Therefore, incorporating both abundant and rare subcommunities can provide a more comprehensive evaluation of the impact of microbial communities on vertical lake sediments.

High-altitude lakes, such as those on the Tibetan plateau, are more sensitive to the impacts of human activities and climatic change due to their limited water exchange compared with low-altitude lakes (Sala et al., 2000; Messerli et al., 2004). Freshwater lakes make up a significant portion of the total lake area on the Tibetan plateau (Zheng et al., 1993) and play important roles in ecosystem balance (Liu et al., 2009). Most freshwater sediments have many different environmental niches, even on a millimeter-scale (Spring et al., 2000). The sediment microbial community in these lakes can serve as an indicator of ecosystem health and environmental conditions (Dai et al., 2016; Zhang et al., 2019). Further studies on the vertical β-diversity and microbial interactions could deepen our understanding of the potential ecological mechanism governing microbial community in freshwater lake sediments. In this study, we obtained sediment cores from two freshwater lakes, Mugecuo Lake and Cuopu Lake. Vertical variations in microbial communities were analyzed through quantitative PCR (qPCR) and high-throughput 16S rRNA gene sequencing techniques. Our objective is to examine the following: (1) The relationship between microbial community composition and potential function with depth in freshwater lake sediments, (2) The change in vertical beta diversity with depth, and (3) The shifts in microbial interactions with sediment depth and nutrient concentration.

Materials and methods

Study site and sampling

Mugecuo Lake (MGC, 101.8414 E–101.8708 E, 30.1626 N–30.1483 N) and Cuopu Lake (CP, 99.3344E–99.6458E, 30.3731 N–30.6353 N) are freshwater lakes located on the eastern Tibetan plateau (Supplementary Figure S1). MGC and CP have altitudes of 3,793 and 4,230 m, water area of 4 km² and 0.9 m², and maximum water depth of 70 and 15 m, respectively, (Ma et al., 2019; Jiang et al., 2020). In October 2017, lake sediment cores from MGC lake (30.1492 N, 101.8531E) and CP lake (30.4917 N, 99.5442E) were sampled using a Gravity sampler. The water depth at the sampling location at MGC was ~35 m and at LC was ~14 m. The collected sediment cores of were 38.5 and 45 cm long for MGC and CP, respectively. The cores were sliced into 76 and 45 samples at 0.5 or 1 cm intervals for MGC and CP, respectively. Each sample was divided into two subsamples—one for analyzing soil physicochemical index and the other for molecular biological analysis.

Sediment properties measurement

The sediment samples were analyzed for various properties, including total organic carbon (TOC), total nitrogen (TN), ammonia nitrogen (NH₄⁺-N), nitrate nitrogen (NO₃⁻-N), and inorganic phosphorous (IP). The TOC was detected using a TOC analyzer (TOC-L&SSM, Shimadzu Corp, Japan) with the dry combustion method. The TN was measured using an ultraviolet spectrophotometer (UV 3600, Shimadzu Corp, Japan) with the sulfate digestion method. The NH₄⁺-N and nitrate nitrogen NO₃⁻-N were measured using a continuous flow analyzer (Skalar, Breda, and Netherlands). The IP was measured following the Standards Measurements and Testing (SMT) protocol (Wu et al., 2018).

In addition, the specific activity of ²²⁶Ra and ²¹⁰Pb was measured using a high-resolution high purity Germanium (HPGe) γ-spectrometer (EG&GORTEC, GWL-120-15, United States). The age of the sediment deposition was calculated using the constant rate of supply (CRS) of the ²¹⁰Pb model (Sanchez-Cabeza and Ruiz-Fernández, 2012).

$t=\frac{1}{\lambda }\times \ln \left(\frac{A_h}{A_0}\right)$

Where λ is the decay constant of ²¹⁰Pb (λ = 0.03114/a), A_h is the ²¹⁰Pb_ex activity of each sedimentary layer, and A₀ is the total activity of ²¹⁰Pb in the sedimentary core.

DNA extraction, quantitative PCR, and Illumina sequencing

To minimize the risk of DNA contamination, all the molecular experiments were conducted in a sterile environment using 75% ethanol-cleaned and UV sterilized laboratory tools. The total DNA of lake sediments was extracted using the FastDNA SPIN Kit (MP Biomedicals, United States), and its concentration was quantified using a Nanodrop ultraviolet spectrophotometer (Thermo Scientific, Wilmington, United States). The extracted DNA was stored at −20°C until further analysis.

The abundance of 16S rRNA gene in the sediment samples was detected by quantitative PCR (qPCR) using 515F (5′-GTGCCAG CMGCCGCGGTAA-3′) and 907R (5′-CCGYCAATTYMTTTRA GTTT-3′) primer pairs (Kolb et al., 2003). The qPCR reactions were performed on a CFX96 Real-Time System (BioRad, Laboratories Inc., United States), with an initial denaturation (95°C, 3 min) followed by 40 cycles of denaturation at 95°C for 10 s, annealing at 55°C for 20 s, and extension at 72°C for 20 s. The qPCR reaction volume included 2 μL of template DNA, 0.4 μL of each primer (10 μM), 10 μL of SYBR Green and ROX mixture (2×, Accurate Biotechnology Co., Hunan, China), and 7.2 μL of sterile water. Melting curve analysis was used to test the specificity of PCR products. The R² and amplification efficiency of the standard curve were 0.996 and 99%, respectively.

PCR amplicon was carried out targeting 16S rRNA gene V4–V5 region using primer pairs 515F and 907R (Kolb et al., 2003). The PCR was performed in a 50 μL reaction volume, including 2 μL of the template DNA, 1 μL of each primer (10 μM), 25 μL of EasyTaq ® PCR SuperMix (2×), and 22 μL of sterile water. The PCR cycling was performed on a S100™ Thermal Cycler (BioRad, Laboratories Inc., United States), with an initial denaturation (95°C, 1 min) followed by 28 cycles (denaturation at 95°C for 20 s, annealing at 55°C for 20 s, and extension at 72°C for 30 s), followed by a final extension step (72°C for 10 min). The PCR products were purified with Omega Bio-Tek E.Z.N.A.® Cycle Pure Kit and quantified using Nanodrop ultraviolet spectrophotometer (Thermo Scientific, Wilmington, United States). Finally, the purified DNA was sequenced on an Illumina HiSeq2500 platform at Magigene Biotechnology Co. Ltd. (Guangzhou, China).

Sequencing analysis

The 16S rRNA raw data were processed using USEARCH (v10.0) and Mothur (v1.45.2) to perform the sequencing analysis. Firstly, the paired reads of raw sequences were merged using USEARCH’s-fastq_mergepairs command, and the forward and reverse primers were removed. After filtering quality, sequences with a maximum expected error greater than one were excluded. The Operational taxonomic units (OTUs) were generated by clustering the sequence at 97% sequence similarity using USEARCH’s-cluster_otus command. The obtained OTU table was then classified using the command Mothur’s classify. Seqs command, using the SILVA 138 reference database. The alpha diversity indices, including Chao1, Shannon, Simpson, and Richness, were calculated in USEARCH using the-alpha_div command based on the normalized OTU table. All the 16S rRNA sequences have been deposited in the SRA database of the national center for biotechnology information (NCBI) with the accession number PRJNA876099.

Statistical analyses

To assess the microbial community by considering subcommunities with different relative abundance, we defined the “abundant” OTUs as the OTUs with relative abundance >0.1% in more than 50% of the samples and present in more than 80% of the samples. “Rare” OTUs refer to OTUs with an average relative abundance of <0.1% in all samples. “Moderate” OTUs refer to other OTUs except for rare and abundant OTUs (Galand et al., 2009; Wu et al., 2019; Dai et al., 2022). All further statistical analyses were performed in R (v 4.1.2, http://www.r-project.org).

The β-diversity and its distribution into replacement (β_repl) and richness difference (β_rich) components were calculated based on the Jaccard dissimilarity via the “adespatial” package in R (Shen et al., 2020). Vertical β-diversity was calculated as the β-diversity between the sum of OTUs in the surface (1 cm) layer and each deeper layer below 1 cm (Wurzbacher et al., 2017; Zhang Y. et al., 2021).

The representative OTUs-based heatmap was to show the similarity of the microbial community across different lake layers and to identify the depth trends. Hellinger transformation of the OTU table and Principal components analysis (PCA) was conducted. To construct the heatmap, representative OTUs were selected with the highest loadings of PC1, PC2, and PC3. The heatmap was created based on Manhattan distance and the clustering algorithm of ward.D using the function “heatmap.2” of the “gplots” package in R. Nonmetric multidimensional scaling (NMDS) was constructed to display differences in microbial communities among different lakes and depths and was plotted using the “phyloseq” package based on Bray–Curtis distance.

The normalized stochasticity ratio (NST) was used to estimate the determinacy and stochasticity of the microbial assembly processes in two vertical lake sediments (Ning et al., 2019). The NST index was calculated based on the null-model analysis, with 50% as the boundary point between more deterministic (<50%) and more stochastic (>50%) assembly processes (Ning et al., 2019). A canonical analysis of principal coordinates (CAP), also known as distance-based redundancy analysis (db-RDA), was used to explain further the relationship between environmental factors and microbial communities in the two lakes. Besides, the relationship between environmental factors and vertical β-diversity and its components (replacement and richness) was analyzed by CAP. CAP was plotted based on Bray–Curtis distance using the “capscale” function of the “vegan” package. Spearman correlation analysis was performed to show the relationship between environmental factors and alpha diversity of microbial communities using the function “corr.test” of the “psych” package in R.

The function annotation of the soil bacterial taxonomy was evaluated based on the Functional Annotation of Prokaryotic Taxa Database (Louca et al., 2016). The top 20 functional groups with the highest proportions of OTUs were visualized using the function “heatmap” in the software Tbtools (Chen et al., 2020). To analyze the depth variation of microbial co-occurrence patterns, networks were constructed for two depth groups or all lake samples based on OTUs present in more than 30% of samples (Dai et al., 2022). The correlations of OTUs were calculated using Spearman’s rank-based correlations, and only highly significant correlations between OTUs were retained (correlation coefficient ρ > 0.7 and ρ < −0.7, p-value <0.05) in the network. The network was performed by the “igraph” package in R (Hartvigsen, 2011) and visualized using the interactive platform Gephi (Bastian et al., 2009).

Results

Environmental factors of sediments in two lakes

Depth distributions of the physicochemical properties in the sediment of two lakes are shown in Figure 1. The results showed that the IP contents in lake MGC decreased with depth, ranging from 836.5 to 49.5 mg kg⁻¹, while TN, TOC, NO-3-N, and NH+ 4-N contents did not directly correlate with depth. In lake CP, TN and TOC contents decreased with depth and dropped sharply at 25–30 cm, from 9.97 to 0.94 mg g⁻¹ and 52.50 to 4.60 mg g⁻¹, respectively. IP and NO-3-N contents were also decreased with depth, while NH+ 4-N content was not significantly correlated with depth. The ²¹⁰Pb dating indicated that the sedimentary cores in both lakes were deposited over the past ca. 167 and 176 years (Figure 1).

FIGURE 1

Figure 1. Depth profiles of key physicochemical characteristics (TOC, TN, IP, NO₃⁻-N, NH₄⁺-N, and depositional age) and microbial characteristics (copy numbers of 16S rRNA genes), the relative abundance of Cyanobacteria, and alpha diversity indices (Chao1, Richness, Shannon, and Simpson).

Alpha-diversity and abundance of bacteria varied with depth

The abundance of bacteria in the sediments of both lakes decreased with depth (Figure 1), as quantified using the qPCR method. The copy numbers of 16S rRNA genes in lake CP ranged from 2.65 × 10⁸ to 2.21 × 10¹¹ g⁻¹ dry soil, while in lake MGC ranged from 1.61 × 10⁹ to 8.96 × 10¹⁰ g⁻¹ dry soil.

Alpha diversity indices, including Chao1, Richness, Simpson, and Shannon, were calculated based on the normalized OTU table (Figure 1). Lake CP had higher Chao1, Richness, and Shannon indices than lake MGC. Simpson indices were the opposite. The Chao1, Richness, and Shannon indices all negatively correlated with depth in lake MGC (p < 0.001, Supplementary Table S1). Conversely, the Simpson indices increased with depth (p < 0.001, Supplementary Table S1). The Chao1, Richness, and Shannon indices in lake CP remained constant or very slightly increased with depth between 1 and 35 cm and declined below 35 cm.

Moreover, Spearman correlation analysis showed the association between alpha diversity and environmental factors (Supplementary Table S1), with Chao1 and Shannon indices negatively correlated with depositional age and NO− 3-N (p < 0.001) In lake MGC, and positively correlated with TN and IP (p < 0.05). The Simpson index in lake MGC showed an opposite relationship with environmental factors. In lake CP, the correlation between the four alpha diversity indices and environmental factors was not significant (p > 0.05).

Bacterial community composition in phylum- and class-level

Operational taxonomic units were classified into 65 bacterial phyla in sediments. The phyla with relative abundance higher than 1% in the sediments of both lakes were regarded as the abundant phyla (Figure 2). Proteobacteria was the most abundant phylum in both lakes, accounting for 62.00 and 38.42% of bacteria in the MGC and CP lakes, respectively. Gammaproteobacteria was the largest bacterial class in both lakes, accounting for 57.99 and 30.07% in MGC and CP, respectively. In MGC, Chloroflexi (9.74%) was the second most abundant phylum, followed by Planctomycetes (4.35%) and Bacteroidetes (3.10%). In contrast, Cyanobacteria (15.79%) was the second most abundant phylum in lake CP, followed by Chloroflexi (7.23%) and Bacteroidetes (6.49%).

FIGURE 2

Figure 2. The relative abundance of dominant phyla in sediment cores. Only phyla with an average relative abundance >1% were shown.

In lake MGC, the relative abundance of Gammaproteobacteria increased with depth (Supplementary Figure S2, p < 0.001). Cyanobacteria, Bacteroidetes, and Deltaproteobacteria were enriched in shallow layers, and their relative abundance declined with depth (Supplementary Figure S2, p < 0.001). In lake CP, the relative abundance of Cyanobacteria was much higher than that in lake MGC, even accounting for 8.30% (average relative abundance) below 25 cm in CP lake.

We separated the OTUs into abundant, moderate, and rare OTUs. The abundant OTUs mainly belonged to Proteobacteria in lake MGC (37.74%) and CP (33.44%). The relative abundance of the abundant OTUs in total reads was high, exceeding 60% in each lake (Supplementary Figure S3). However, the number of these abundant OTUs accounts for less than 1% of the total OTUs (Supplementary Figure S3).

The β-diversity of vertical sediments

The total β-diversity increased with depth, between 0.62–0.90 and 0.57–0.85 in MGC and CP lakes, respectively, (Figure 3). In lake MGC, species replacement (β_repl) generally decreased with depth, and richness differences (β_rich) increased with depth (Figure 3). Conversely, in lake CP, species replacement (β_repl) was consistently high and increased with depth (35–45 cm) in the deep layer, and richness differences (β_rich) were the opposite, which remained at a lower level and decreased with depth at 35–45 cm (Figure 3). The trend in β-diversity and its components (β_repl and β_rich) of the rare OTUs was consistent with the entire community in each lake. In contrast, the β-diversity and its components (β_repl and β_rich) of abundant OTUs remained stable and were consistently low in both lakes.

FIGURE 3

Figure 3. Vertical β-diversity and its components (replacement and richness) based on the total OTUs and partial OTUs (abundant OTUs, moderate OTUs, and rare OTUs).

The heatmap analysis based on representative OTUs showed that the microbial communities in each lake had a depth-related trend (Supplementary Figure S4), with samples from deeper layer being separated from shallow layers. In the MGC lake, samples from various layers were divided into shallow layers (1–17.5 cm) and deep layers (18–38.5 cm). The heatmap of lake CP showed that the samples of shallow layers (1–22 cm) were separated from samples of deep layers (23–45 cm). The OTUs in class Gammaproteobacteria were abundant in both shallow and deep layers of MGC. However, the different OTUs in Gammaproteobacteria have different depth-related distributions. OTU4, OTU469, and OTU14 were only abundant in the deep layer, while OTU2, OTU17, and OTU18 were enriched in both shallow and deep layers. The depth distribution of Gammaproteobacteria OTUs in lake CP was consistent with that in lake MGC. In addition, OTUs belonging to Cyanobacteria, which include OTU7 and OTU43, have a higher relative abundance in the shallow layers (1–22 cm) of lake CP.

Moreover, NMDS revealed the difference in microbial community structure among samples with different lakes and depths (Supplementary Figure S5). There was a clear separation of microbial community composition between lake MGC and lake CP by NMDS (R = 0.72, p = 0.001). Furthermore, depth-related trends of microbial communities were observed in these two lakes (R = 0.15, p = 0.001).

Potential function prediction of microbial community

Potential function prediction was calculated by FAPROTAX analysis. A total of 1,184 OTUs (9.29% of all OTUs) were assigned to 59 functional groups. The results of the predictions showed that the potential functions varied and changed with the sediment depth (Figure 4).

FIGURE 4

Figure 4. Predicted functional groups of soil microbial community. Only the top 20 functional groups were shown in the plot. The relative abundances of the functional groups were indicated by bubble size.

The functions involved in the carbon cycle were the primary and abundant functions. In lake MGC, photoautotrophy was abundant in surface sediments (1–5 cm), while chloroplast was most abundant in shallow sediments (1–25 cm) of lake CP. Chemoheterotrophy, aerobic chemoheterotrophy, and ureolysis were the main functions in deeper sediments of lake MGC (15–38.5 cm) and lake CP (30–45 cm). In the sediments below 10 cm, the functional groups involved in methanogenesis were more abundant. In addition, for the functions involved in the nitrogen cycle, nitrification was active in the surface sediments (1–5 cm). In comparison, nitrate reduction was more abundant in the deep sediments below 35 cm in lake MGC. In the sediments above 25 cm in lake CP, the functional groups involved in the sulfur cycle, such as the respiration of sulfur compounds and sulfite respiration, were more abundant (Figure 4).

Correlation between microbial community and vertical β-diversity with environmental factors

The CAP was performed to illustrate the influence of environmental factors on microbial communities and vertical β-diversity (Supplementary Figure S6). All environmental factors, including TN, TC, IP, NO− 3-N, and NH+ 4-N, depth and depositional age, accounted for 18.92 and 62.42% of the total variations for microbial communities in lake MGC and CP, respectively, in Supplementary Figure S6A. CAP1 and CAP2 are the first two constrained axes in CAP coordinates, which account for 52.73 and 18.22% of the variation for lake MGC and 78.54 and 15.56% of the variation for lake CP, respectively. This result was consistent with NST, which indicated that the community assembly of lake CP (NST = 46.41%) was more deterministic than that of lake MGC (NST = 75.57%). For the single environmental factor, depth was found to be the most critical factor in both lakes, accounting for a significant portion of the total variation.

The CAP results showed that environmental factors accounted for 68.05 and 81.93% of the total variations in β-diversity in MGC and CP, respectively, (Supplementary Figure S6B). Depth is the most important factor, contributing 60.74 and 48.80% of the variation in β-diversity in MGC and CP (p < 0.001). In MGC, the second crucial variable was IP contents, which explained 53.83% of the variation (p < 0.001), while in CP, depositional age and TOC were the second and third crucial variables, explaining 46.17 and 36.92% of the variation, respectively. In terms of the relationship between replacement and richness with environmental factors, the results differ between MGC and CP. In MGC, replacement increased with IP but decreased with depth and NO− 3-N, while richness showed the opposite trend. In CP, replacement increased with depth but decreased with TOC, TN, NO− 3-N, and IP, and richness had a weak correlation with these environmental factors.

Co-occurrence network analysis

Networks were built based on all samples in each lake to explore co-occurrence patterns among the microbial community (Supplementary Figure S7). The network of lake MGC consisted of 522 nodes and 10,293 links, with main modules represented in different colors. Modules 1 and 2 had high relative abundance in shallow sediments (1–17.5 cm) and mainly consisted of Chloroflexi and Nanoarchaeaeota. The microbial community of module 3 was enriched in deep sediments (18–38.5 cm) and mainly consisted of Chloroflexi and Acidobacteria. Besides, in modules 1 and 4, the abundant class Gammaproteobacteria was enriched in deep sediments (18–38.5 cm). In the network of lake CP, there are 1,173 nodes and 47,402 links. Modules 1 and 4 had higher abundance in deep sediments (23–45 cm). Gammaproteobacteria was abundant in module 1, while Euryarchaea were enriched in module 4. Module 2 mainly consisted of the phylum Cyanobacteria, which was enriched in shallow sediments (1–22 cm).

The networks of the microbial community in surface and deep sediments were established for each lake to analyze the co-occurrence patterns (Figure 5). The network topology indices are shown in Supplementary Table S2. The networks from different depth layers showed that both rare and moderate OTUs accounted for a high proportion in both lakes (Figure 5A). The results showed that the degree and betweenness centrality of abundant OTUs were significantly higher compared to rare and moderate OTUs (Figure 5B, p < 0.001, Kruskal–Wallis test). The positive and negative interactions among OTUs were also calculated, with a decrease in the proportion of negative interactions observed with increasing depth. Furthermore, negative interactions mainly existed in moderate/rare OTUs or between rare and moderate OTUs (Figure 5B).

FIGURE 5

Figure 5. The co-occurrence network and its topological properties based on OTUs of each depth group in lake sediment. (A) The networks for different depth groups. The network name, e.g., MGC-S and MGC-D, S and D represent the shallow and deep depth group, respectively. Network node represents a bacterial taxon at the OTU level, and the size of the node represents the “degree” (the number of connections). Nodes are denoted by different colors representing their affiliated taxa (abundant taxa, moderate taxa, and rare taxa). The significant spearman correlations between OTUs are retained (correlation coefficient ρ > 0.7 and <−0.7, p < 0.05). The red and green edges signify positive and negative correlations, respectively. (B) The topological properties of the network, including nodes’ degree, betweenness centrality, and proportion of network links. The difference in nodes’ degree and betweenness centrality level within different taxa were tested by Kruskal–Wallis test. The proportion of network links was shown by the affiliated taxa of the associated OTUs. AT, abundant taxa; MT, moderate taxa; RT, rare taxa.

Discussion

Vertical variation of microbial community composition in two lakes

Microbial communities in two lake sediments varied with the depth due to the depth gradient of biogeochemical properties, which provided specific niches for microorganisms with diverse metabolic capabilities (Wilms et al., 2006). The microbial community can be grouped into different clusters based on depth (Supplementary Figures S4, S7), a phenomenon that has been observed in other lake sediments (Wurzbacher et al., 2017). The shallow layer clusters mainly consist of the phyla Proteobacteria, Cyanobacteria and Bacteroidetes. Proteobacteria are commonly abundant in lake sediments, and the composition of their significant classes and their proportions varies with lakes (Zhang et al., 2015). In lake CP, Cyanobacteria was consistently abundant at depth of 1–35 cm. Previous studies have noted the predominance of Cyanobacteria in surface sediment, which may be related to cyanobacterial blooms (Latour et al., 2007; Liu et al., 2010). Lake CP was oligotrophic before the 1880s, but its nutrient conditions improved over time, which may have contributed to the increased growth of algae (Jiang et al., 2020). Cyanobacteria can also adapt to the cold, dark, and anoxic deep sediment, possibly due to its cold-shock DEAD-box protein A, which has been shown to play a role in cyanobacterial adaptation to low temperatures (Chamot et al., 1999; Redder et al., 2015). Bacteroidetes, which are chemoorganotrophic bacteria with phototrophic capabilities, are also commonly found in surface lake sediments (Newton et al., 2011) and positively correlated with DOC concentration (Eiler et al., 2003), likely explaining their predominance in shallow layers.

Vertical variation of β-diversity and its components in lake sediments

Exploring the β-diversity and its components (replacement and richness) can provide insights into the processes driving the vertical variation of microbial communities. In this study, the richness effect constantly increased with depth while the replacement effect decreased with depth in lake MGC (Figure 3). This trend reflected that the microbial community in the deep layer was a subgroup of the surface microbial community, with the community shifting from a complicated and rich community at the surface to a more atrophic community in the deeper sediment. This transition was mainly driven by the filter effect, which might result from abiotic limitations, such as nutrient reductions and space shrinking (Hu et al., 2012). In contrast, the trend of replacement and richness effect in lake CP was the opposite. The richness effect maintained a low level, while the replacement effect was constantly abundant and increased with depth (Figure 3). This trend of increased replacement effects in the deep layer has also been found in marine sediment (Zhang Y. et al., 2021), which may indicate the existence of an inactive seed bank with a great variety in the deep layer. The trend of β-diversity components in lake CP was consistent with the microbial community structural model that proposed niche specialists adapted to the specific environmental conditions at different depths of layers (Wurzbacher et al., 2017). During the burial process for some lake sediments, the different patterns of environment parameters with depth may result in different β-diversity portioning patterns in vertical sediments. Our results suggested that most environmental factors in MGC were not depth-related, and only IP contents decreased with increasing depth. IP was found to play a critical role in shaping the β-diversity pattern in MGC. While in CP, most depth-related environmental factors such as TOC and TN had a significant impact on the β-diversity portioning pattern. Furthermore, spatial and biotic factors also influenced the β-diversity portioning patterns in lake sediments (Wu et al., 2020). Further research can be conducted in the future by incorporating more parameters to further explore these patterns.

The co-occurrence network of microbial interactions at various depths of lake sediments

In this study, the proportion of negative correlation among OTUs in the network decreased with sediment depth in both MGC and CP (Figure 5). Resource competition and metabolic cross-feeding are the main drivers of the microbial community assembly (Machado et al., 2021). The interactions between microbial communities are significantly influenced by the availability of nutrients (Gulis and Suberkropp, 2003; West et al., 2007). Most microorganisms are auxotrophic and require external nutrients for growth. In a nutrient-rich environment, active microbial growth may result in competition for necessary nutrients and the inhibition of the growth of competitor growth in the community (Ratzke et al., 2020). In comparison, collaboration may encourage microbial colonization and growth in low-nutrient conditions, facilitating adaptation to oligotrophic environments (Velez et al., 2018). Lake sediments can provide microbes with different nutrient environments at various depths, such as decreasing availability of electron acceptors and substrates with depth. Consequently, negative microbial interactions decreased with depth due to the reduction of major nutrients including IP in MGC and TOC, TN, and IP in CP. This phenomenon was also found in marine and other lake sediments (Dai et al., 2022; Zhang et al., 2022), where negative microbial interactions were observed to be more prevalent in rich-nutrient conditions, and positive interactions were more predominant in low-nutrient conditions.

Our analysis also revealed that the number of rare OTUs accounted for a higher proportion of total nodes, and that most negative correlations existed among moderate or rare OTUs. This suggested that rare and moderate OTUs have more overlapping ecological niches. Besides, abundant OTUs had higher degree and betweenness centrality levels, which aligns with previous findings (Zhang et al., 2020, Zhang Y. et al, 2021), indicating that abundant OTUs play a critical role in maintaining the structure of microbial interactions in the network.

Conclusion

In this study, we conducted a comprehensive analysis of the microbial community composition, diversity, and interaction in sediment samples collected from two lakes on the Tibetan plateau. Our findings demonstrate the significance of both abundant and rare taxa in shaping the β-diversity and microbial interactions. Our results reveal similarities in microbial community compositions at the phylum level between the two lakes, with Cyanobacteria being more prevalent in lake CP. Additionally, we observed varying β-diversity patterns correlated with environmental factors. Our results indicate that the pattern of microbial interactions varied with depth, with evidence suggesting that higher nutrient levels may lead to negative interactions among microbes, while positive interactions are more prevalent in low-nutrient environments. These results broaden our understanding of patterns of β-diversity partitioning and microbial interactions in vertical sediment samples on the Tibetan plateau and provide new insights into mechanisms driving microbial community variation.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, PRJNA876099.

Author contributions

YD, TH, and CheH designed the study. TH and ZZ sampled the sediments. TH and CheH were responsible for microbial lab work. XZ and YD analyzed the data, prepared figures, and wrote the initial manuscript. LC, KL, YL, and ChaH reviewed and edited the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 41971077) and the Second Tibetan Plateau Scientific Expedition and Research (STEP) program (Grant No. 2019QZKK0503).

Conflict of interest

The authors declare that the research was conducted without any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1118892/full#supplementary-material

References

Ambrosetti, W., Barbanti, L., and Sala, N. (2003). Residence time and physical processes in lakes. J. Limnol. 62, 1–15. doi: 10.4081/jlimnol.2003.s1.1

CrossRef Full Text | Google Scholar

Baselga, A. S. (2007). Disentangling distance decay of similarity from richness gradients: response to Soininen. Ecography 30, 838–841. doi: 10.1111/j.2007.0906-7590.05387.x

CrossRef Full Text | Google Scholar

Baselga, A. (2010). Partitioning the turnover and nestedness components of beta diversity. Glob. Ecol. Biogeogr. 19, 134–143. doi: 10.1111/j.1466-8238.2009.00490.x

CrossRef Full Text | Google Scholar

Bastian, M., Heymann, S., and Jacomy, M., Gephi: an open source software for exploring and manipulating networks, Proceedings of the International AAAI Conference on Web and Social Media; 2009. 3. p. 361–362

Google Scholar

Bauer, M. A., Kainz, K., Carmona-Gutierrez, D., and Madeo, F. (2018). Microbial wars: competition in ecological niches and within the microbiome. Microb. Cell. 5, 215–219. doi: 10.15698/mic2018.05.628

PubMed Abstract | CrossRef Full Text | Google Scholar

Bergk Pinto, B., Maccario, L., Dommergue, A., Vogel, T. M., and Larose, C. (2019). Do organic substrates drive microbial community interactions in Arctic snow? Front. Microbiol. 10:2492. doi: 10.3389/fmicb.2019.02492

PubMed Abstract | CrossRef Full Text | Google Scholar

Chamot, D., Magee, W. C., Yu, E., and Owttrim, G. W. (1999). A cold shock-induced cyanobacterial RNA helicase. J. Bacteriol. 181, 1728–1732. doi: 10.1128/JB.181.6.1728-1732.1999

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C., Chen, H., Zhang, Y., Thomas, H. R., Frank, M. H., He, Y., et al. (2020). TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 13, 1194–1202. doi: 10.1016/j.molp.2020.06.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, T., Wen, D., Bates, C. T., Wu, L., Guo, X., Liu, S., et al. (2022). Nutrient supply controls the linkage between species abundance and ecological interactions in marine bacterial communities. Nat. Commun. 13:175. doi: 10.1038/s41467-021-27857-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, Y., Yang, Y., Wu, Z., Feng, Q., Xie, S., and Liu, Y. (2016). Spatiotemporal variation of planktonic and sediment bacterial assemblages in two plateau freshwater lakes at different trophic status. Appl. Microbiol. Biotechnol. 100, 4161–4175. doi: 10.1007/s00253-015-7253-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Dean, W. E., and Gorham, E. (1998). Magnitude and significance of carbon burial in lakes, reservoirs, and peatlands. Geology 26, 535–538. doi: 10.1130/0091-7613(1998)026<0535:MASOCB>2.3.CO;2

CrossRef Full Text | Google Scholar

Debroas, D., Hugoni, M., and Domaizon, I. (2015). Evidence for an active rare biosphere within freshwater protists community. Mol. Ecol. 24, 1236–1247. doi: 10.1111/mec.13116

PubMed Abstract | CrossRef Full Text | Google Scholar

D'Hondt, S., Pockalny, R., Fulfer, V. M., and Spivack, A. J. (2019). Subseafloor life and its biogeochemical impacts. Nat. Commun. 10, 3519–3513. doi: 10.1038/s41467-019-11450-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Edlund, A., Soule, T., Sjöling, S., and Jansson, J. K. (2006). Microbial community structure in polluted Baltic Sea sediments. Environ. Microbiol. 8, 223–232. doi: 10.1111/j.1462-2920.2005.00887.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Eiler, A., Langenheder, S., Bertilsson, S., and Tranvik, L. J. (2003). Heterotrophic bacterial growth efficiency and community structure at different natural organic carbon concentrations. Appl. Environ. Microbiol. 69, 3701–3709. doi: 10.1128/AEM.69.7.3701-3709.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Galand, P. E., Casamayor, E. O., Kirchman, D. L., and Lovejoy, C. (2009). Ecology of the rare microbial biosphere of the Arctic Ocean. Proc. Natl. Acad. Sci. U. S. A. 106, 22427–22432. doi: 10.1073/pnas.0908284106

PubMed Abstract | CrossRef Full Text | Google Scholar

Gulis, V., and Suberkropp, K. (2003). Interactions between stream fungi and bacteria associated with decomposing leaf litter at different levels of nutrient availability. Aquat. Microb. Ecol. 30, 149–157. doi: 10.3354/ame030149

CrossRef Full Text | Google Scholar

Harrison, S., Ross, S. J., and Lawton, J. H. (1992). Beta diversity on geographic gradients in Britain. J. Anim. Ecol. 61:151. doi: 10.2307/5518

CrossRef Full Text | Google Scholar

Hartvigsen, G. (2011). Using R to build and assess network models in biology. Math. Model Nat. Phenom. 6, 61–75. doi: 10.1051/mmnp/20116604

CrossRef Full Text | Google Scholar

Hester, E. R., Jetten, M. S. M., Welte, C. U., and Lücker, S. (2019). Metabolic overlap in environmentally diverse microbial communities. Front. Genet. 10:989. doi: 10.1101/653881

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, G., Wu, J., Feeley, K. J., Xu, G., and Yu, M. (2012). The effects of landscape variables on the species-area relationship during late-stage habitat fragmentation. PLoS One 7:e43894. doi: 10.1371/journal.pone.0043894

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Q., Li, S., Chen, Z., Huang, C., Wu, W., Wan, H., et al. (2020). Disturbance mechanisms of lacustrine organic carbon burial: case study of Cuopu Lake, Southwest China. Sci. Total Environ. 746:140615. doi: 10.1016/j.scitotenv.2020.140615

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiao, S., and Lu, Y. (2020). Soil pH and temperature regulate assembly processes of abundant and rare bacterial communities in agricultural ecosystems. Environ. Microbiol. 22, 1052–1065. doi: 10.1111/1462-2920.14815

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan, N., Maezato, Y., McClure, R. S., Brislawn, C. J., Mobberley, J. M., Isern, N., et al. (2018). Phenotypic responses to interspecies competition and commensalism in a naturally-derived microbial co-culture. Sci. Rep. 8, 297–2322. doi: 10.1038/s41598-017-18630-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kolb, S., Knief, C., Stubner, S., and Conrad, R. (2003). Quantitative detection of methanotrophs in soil by novel pmoA-targeted real-time PCR assays. Appl. Environ. Microbiol. 69, 2423–2429. doi: 10.1128/AEM.69.5.2423-2429.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Latour, D., Salençon, M.-J., Reyss, J.-L., and Giraudet, H. (2007). Sedimentary imprint Ofmicrocystis Aeruginosa(cyanobacteria) blooms in Grangent reservoir (Loire, France). J. Phycol. 43, 417–425. doi: 10.1111/j.1529-8817.2007.00343.x

CrossRef Full Text | Google Scholar

Legendre, P. (2014). Interpreting the replacement and richness difference components of beta diversity. Glob. Ecol. Biogeogr. 23, 1324–1334. doi: 10.1111/geb.12207

CrossRef Full Text | Google Scholar

Liu, L., Peng, Y., Zheng, X., Xiao, L., and Yang, L. (2010). Vertical structure of bacterial and Archaeal communities within the sediment of a eutrophic Lake as revealed by culture-independent methods. J. Freshw. Ecol. 25, 565–573. doi: 10.1080/02705060.2010.9664406

CrossRef Full Text | Google Scholar

Liu, Y., Yao, T., Zhu, L., Jiao, N., Liu, X., Zeng, Y., et al. (2009). Bacterial diversity of freshwater alpine Lake puma Yumco on the Tibetan plateau. Geomicrobiol J. 26, 131–145. doi: 10.1080/01490450802660201

CrossRef Full Text | Google Scholar

Logares, R., Audic, S., Bass, D., Bittner, L., Boutte, C., Christen, R., et al. (2014). Patterns of rare and abundant marine microbial eukaryotes. Curr. Biol. 24, 813–821. doi: 10.1016/j.cub.2014.02.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Louca, S., Parfrey, L. W., and Doebeli, M. (2016). Decoupling function and taxonomy in the global ocean microbiome. Science 353, 1272–1277. doi: 10.1126/science.aaf4507

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, X., Hu, W., Smilauer, P., Yin, M., and Wolinska, J. (2019). Daphnia galeata and D. dentifera are geographically and ecologically separated whereas their hybrids occur in intermediate habitats: a survey of 44 Chinese lakes. Mol. Ecol. 28, 785–802. doi: 10.1111/mec.14991

PubMed Abstract | CrossRef Full Text | Google Scholar

Machado, D., Maistrenko, O. M., Andrejev, S., Kim, Y., Bork, P., Patil, K. R., et al. (2021). Polarization of microbial communities between competitive and cooperative metabolism. Nat. Ecol. Evol. 5, 195–203. doi: 10.1038/s41559-020-01353-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Mason, O. U., Meo-Savoie, D., Carol, A., Van Nostrand, J. D., Zhou, J., Fisk, M. R., et al. (2009). Prokaryotic diversity, distribution, and insights into their role in biogeochemical cycling in marine basalts. ISME J. 3, 231–242. doi: 10.1038/ismej.2008.92

PubMed Abstract | CrossRef Full Text | Google Scholar

Messerli, B., Viviroli, D., and Weingartner, R. (2004). Mountains of the world: vulnerable water towers for the 21st century. AMBIO: a journal of the human. Environment 33, 29–34. doi: 10.1007/0044-7447-33.sp13.29

CrossRef Full Text | Google Scholar

Nelson, D. M., Ohene-Adjei, S., Hu, F. S., Cann, I. K. O., and Mackie, R. I. (2007). Bacterial diversity and distribution in the Holocene sediments of a northern temperate lake. Microb. Ecol. 54, 252–263. doi: 10.1007/s00248-006-9195-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Newton, R. J., Jones, S. E., Eiler, A., McMahon, K. D., and Bertilsson, S. (2011). A guide to the natural history of freshwater Lake bacteria. Microbiol. Mol. Biol. Rev. 75, 14–49. doi: 10.1128/MMBR.00028-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Ning, D., Deng, Y., Tiedje, J. M., and Zhou, J. (2019). A general framework for quantitatively assessing ecological stochasticity. Proc. Natl. Acad. Sci. U. S. A. 116, 16892–16898. doi: 10.1073/pnas.1904623116

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedros-Alio, C. (2012). The rare bacterial biosphere. Annu. Rev. Mar. Sci. 4, 449–466. doi: 10.1146/annurev-marine-120710-100948

CrossRef Full Text | Google Scholar

Ratzke, C., Barrere, J., and Gore, J. (2020). Strength of species interactions determines biodiversity and stability in microbial communities. Nat. Ecol. Evol. 4, 376–383. doi: 10.1101/671008

PubMed Abstract | CrossRef Full Text | Google Scholar

Redder, P., Hausmann, S., Khemici, V., Yasrebi, H., and Linder, P. (2015). Bacterial versatility requires DEAD-box RNA helicases. FEMS Microbiol. Rev. 39, 392–412. doi: 10.1093/femsre/fuv011

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, Z., Zhang, C., Li, X., Ma, K., and Cui, B. (2022). Abundant and rare bacterial taxa structuring differently in sediment and water in Thermokarst Lakes in the Yellow River source area, Qinghai-Tibet plateau. Front. Microbiol. 13:774514. doi: 10.3389/fmicb.2022.774514

PubMed Abstract | CrossRef Full Text | Google Scholar

Rissanen, A. J., Peura, S., Mpamah, P. A., Taipale, S., Tiirola, M., Biasi, C., et al. (2019). Vertical stratification of bacteria and archaea in sediments of a small boreal humic lake. FEMS Microbiol. Lett. 366, 1–11. doi: 10.1093/femsle/fnz044

PubMed Abstract | CrossRef Full Text | Google Scholar

Sala, O. E., Stuart Chapin, F., Armesto, J. J., Berlow, E., Bloomfield, J., Dirzo, R., et al. (2000). Global biodiversity scenarios for the year 2100. Science 287, 1770–1774. doi: 10.1126/science.287.5459.1770

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanchez-Cabeza, J. A., and Ruiz-Fernández, A. C. (2012). ²¹⁰Pb sediment radiochronology: an integrated formulation and classification of dating models. Geochim. Cosmochim. Acta 82, 183–200. doi: 10.1016/j.gca.2010.12.024

CrossRef Full Text | Google Scholar

Shen, C., Gunina, A., Luo, Y., Wang, J., He, J. Z., Kuzyakov, Y., et al. (2020). Contrasting patterns and drivers of soil bacterial and fungal diversity across a mountain gradient. Environ. Microbiol. 22, 3287–3301. doi: 10.1111/1462-2920.15090

PubMed Abstract | CrossRef Full Text | Google Scholar

Spring, S., Schulze, R., Overmann, J.R., and Schleifer, K.-H.. (2000). Identification and characterization of ecologically significant prokaryotes in the sediment of freshwater lakes: molecular and cultivation studies. FEMS Microbiol Rev, 24, 573–590. doi: 10.1111/j.1574-6976.2000.tb00559.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Velez, P., Espinosa-Asuar, L., Figueroa, M., Gasca-Pineda, J., Aguirre-von-Wobeser, E., Eguiarte, L. E., et al. (2018). Nutrient dependent cross-kingdom interactions: fungi and bacteria from an oligotrophic desert oasis. Front. Microbiol. 9:1755. doi: 10.3389/fmicb.2018.01755

PubMed Abstract | CrossRef Full Text | Google Scholar

Vuillemin, A., Ariztegui, D., Horn, F., Kallmeyer, J., Orsi, W. D., and Team, P. S. (2018). Microbial community composition along a 50,000-year lacustrine sediment sequence. FEMS Microbiol. Ecol. 94, 1–14. doi: 10.1093/femsec/fiy029

PubMed Abstract | CrossRef Full Text | Google Scholar

West, S. A., Diggle, S. P., Buckling, A., Gardner, A., and Griffin, A. S. (2007). The social lives of microbes. Annu. Rev. Ecol. Evol. Syst. 38, 53–77. doi: 10.1146/annurev.ecolsys.38.091206.095740

CrossRef Full Text | Google Scholar

Whittaker, R. H. (1960). Vegetation of the Siskiyou mountains. Oregon Calif. Ecol. Monogr. 30, 279–338. doi: 10.2307/1943563

CrossRef Full Text | Google Scholar

Wilms, R., Sass, H., Koجêpke, B., Koجêster, J., Cypionka, H., and Engelen, B. (2006). Specific bacterial, archaeal, and eukaryotic communities in tidal-flat sediments along a vertical profile of several meters. Appl. Environ. Microbiol. 72, 2756–2764. doi: 10.1128/AEM.72.4.2756-2764.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Y., Huang, T., Huang, C., Shen, Y., Luo, Y., Yang, H., et al. (2018). Internal loads and bioavailability of phosphorus and nitrogen in Dianchi Lake, China. Chin. Geogr. Sci. 28, 851–862. doi: 10.1007/s11769-018-0994-y

CrossRef Full Text | Google Scholar

Wu, L., Ning, D., Zhang, B., Li, Y., and Sanabria, J. (2019). Global diversity and biogeography of bacterial communities in wastewater treatment plants. Nat. Microbiol. 4, 1183–1195. doi: 10.1038/s41564-019-0426-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, K., Zhao, W., Li, M., Picazo, F., Soininen, J., Shen, J., et al. (2020). Taxonomic dependency of beta diversity components in benthic communities of bacteria, diatoms and chironomids along a water-depth gradient. Sci. Total Environ. 741:140462. doi: 10.1016/j.scitotenv.2020.140462

PubMed Abstract | CrossRef Full Text | Google Scholar

Wurzbacher, C., Fuchs, A., Attermeyer, K., Frindte, K., Grossart, H. P., Hupfer, M., et al. (2017). Shifts among Eukaryota, bacteria, and Archaea define the vertical organization of a lake sediment. Microbiome 5, 41–18. doi: 10.1186/s40168-017-0255-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Hou, F., Xie, W., Wang, K., Zhou, X., Zhang, D., et al. (2020). Interaction and assembly processes of abundant and rare microbial communities during a diatom bloom process. Environ. Microbiol. 22, 1707–1719. doi: 10.1111/1462-2920.14820

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Qi, L., Li, W., Hu, B. X., and Dai, Z. (2021). Bacterial community variations with salinity in the saltwater-intruded estuarine aquifer. Sci. Total Environ. 755:142423. doi: 10.1016/j.scitotenv.2020.142423

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Yang, L., Li, Y., Wang, C., Zhang, W., Wang, L., et al. (2022). Pollution gradients shape the co-occurrence networks and interactions of sedimentary bacterial communities in Taihu Lake, a shallow eutrophic lake. J. Environ. Manag. 305:114380. doi: 10.1016/j.jenvman.2021.114380

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Yao, P., Sun, C., Li, S., Shi, X., Zhang, X. H., et al. (2021). Vertical diversity and association pattern of total, abundant and rare microbial communities in deep-sea sediments. Mol. Ecol. 30, 2800–2816. doi: 10.1111/mec.15937

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L., Zhao, T., Wang, Q., Li, L., Shen, T., and Gao, G. (2019). Bacterial community composition in aquatic and sediment samples with spatiotemporal dynamics in large, shallow, eutrophic Lake Chaohu, China. J. Freshw. Ecol. 34, 575–589. doi: 10.1080/02705060.2019.1635536

CrossRef Full Text | Google Scholar

Zhao, D. Y., Ma, T., Zeng, J., Yan, W. M., Jiang, C. L., Feng, J. W., et al. (2011). Phospholipid fatty acids analysis of the vertical distribution of microbial communities in eutrophic lake sediments. Int. J. Environ. Sci. Technol. 8, 571–580. doi: 10.1007/BF03326243

CrossRef Full Text | Google Scholar

Zheng, M., Tang, J., Liu, J., and Zhang, F. (1993). Chinese saline lakes. Hydrobiologia 267, 23–36. doi: 10.1007/BF00018789

CrossRef Full Text | Google Scholar