Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 30 July 2025

Sec. Aquatic Microbiology

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1614055

Metagenomic insights into surface sediment microbial community and functional composition along a water-depth gradient in a subtropic deep lake

  • 1Department of Municipal Engineering, School of Civil Engineering, Southeast University, Nanjing, China
  • 2State Key Laboratory of Lake and Watershed Science for Water Security, Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences, Nanjing, China

Deep lakes play a critical role in global elemental cycling and serve as habitats for diverse microbial communities. However, studies on the effects of lake stratification on microbial composition and functional potential in surface sediments remain limited. Here, we investigated microbial community structure and functional composition using metagenomics of 38 surface sediments across a depth gradient of 0–90 m in Lugu Lake, China. Our results showed that Shannon diversity peaked at the thermocline for microbial communities, while a U-shaped pattern for functional genes. Microbial communities and functional genes in the surface sediments showed higher spatial heterogeneity at the shallow layer, whereas those at deeper layers tended toward more homogenized. Although water depth was the most important driver in explaining 29.9 and 26.5% of variance in microbial and functional gene composition, stochastic processes primarily governed the community assemblages, particularly dispersal limitation with the contribution of 43.7%. We further found the surface layer was enriched in genes mainly involved in aerobic metabolism and methanogenesis. In contrast, genes related to reduction reactions, including dissimilatory nitrate and sulfate reduction were more abundant in the thermocline and deep layer, reflecting lower redox potential in a deeper layer. Overall, our results provide evidence for microbial community stratification and functional partitioning in deep lakes.

Introduction

Deep lakes are critical ecosystems for global biogeochemical cycling, harboring diverse microbial communities that exhibit remarkable adaptability to stratified environmental conditions (Cabello-Yeves et al., 2019). Vertical gradients in temperature, oxygen, and nutrient availability profoundly influence microbial diversity and community assembly processes in the surface sediments across water layers (Wu et al., 2019; Zhao et al., 2019). Surface layers, characterized by dynamic hydrological cycling and higher oxygen levels, support microbes with high dispersal capacity (Lv et al., 2024). In contrast, deeper layers, with reduced oxygen and lower temperatures, impose strong environmental filtering, limiting microbial diversity but fostering specialized microbial assemblages adapted to extreme conditions (Cabello-Yeves et al., 2023). Recognizing microbial community composition and functional potential in surface sediment of deep lakes is essential for revealing their roles in biogeochemical cycles and their responses to environmental changes in the context of accelerating climate change and increasing anthropogenic influence (Falkowski et al., 2008).

Metabolic potential of microorganisms in the sediments also varies across water layers with distinct environmental conditions (Ayala-Muñoz et al., 2022; Haro-Moreno et al., 2018; Imhoff, 2016). Changes in redox status can drive metabolic transitions as depth increases (Bush et al., 2017; Louca et al., 2016; Peura et al., 2015). In the surface layer with abundant oxygen (Chang et al., 2022), microbes in the water columns rely on aerobic respiration and photosynthesis, with genes associated with carbon fixation, glycolysis, citrate cycling and organic nitrogen metabolism being more abundant (Peura et al., 2018; Yin et al., 2019). In contrast, anaerobic respiration and fermentation predominate in oxygen-depleted deep waters, utilizing alternative electron acceptors such as nitrate, sulfate, and carbon dioxide (Lipsewers et al., 2016; Song et al., 2022). Functional genes related to these reductive processes are expected to be more prevalent in the deeper layer, contributing significantly to the metabolic characteristics of the ecosystem (Voss et al., 2013). Although these microbial processes in water columns are well-studied, their variability and interactions with microbial groups in sediment along the depth gradient require further investigation, particularly in subtropical deep lakes.

Here, we collected 38 surface sediments along a depth gradient of 0–90 m in Lugu Lake, a deep plateau freshwater lake. Using metagenomic sequencing, we analyzed the structure and functional profiles of microbial communities and explored the influence of environmental factors including water depth, physical and chemical characteristics. This study aims to address the following three questions: (1) How do surface sediment microbial communities and functional genes vary along depth gradient in the deep lake? (2) What is the relative importance of stochastic and deterministic processes constraining microbial assembly in different layers? (3) Which functional genes in the surface sediment are enriched along the depth of the overlaying water column? (4) How are functional genes coupled with microbial taxa across layers? Our results elucidated the diversity patterns, assembly mechanisms, and metabolic characteristics of microbial communities in Lugu Lake across the depth gradient, providing evidence for stratification of microbial communities and their functional potential in deep lakes.

Materials and methods

Field sampling

We collected 38 surface sediments and water samples from 0 to 90 meters in August 2010, in Lugu Lake (27°41′–27°45′N, 100°45′–100°50′E), situated in Yunnan Province, China (Supplementary Figure S1) (Zhao et al., 2023). Lugu Lake is one of the deepest freshwater plateau lakes in the region, with a maximum depth of 93.5 meters and the water surface area of 50.5 km2 (Wu et al., 2020). In addition, the lake is situated at an elevation of 2,685 meters, with a catchment area of around 171.4 km2 (Zhao et al., 2023). Unlike many lakes in the surrounding area, Lugu Lake remains ice-free throughout the year, due to its warm temperate climate and semi-enclosed nature (Su et al., 2022). During the winter season, the lake exhibits a vertically uniform temperature profile, whereas other seasons feature thermal stratification (Liu et al., 2019). These conditions make Lugu Lake an ideal model for studying microbial community structure and functional enrichment in deep lakes, especially under high-altitude influences.

The detailed procedures for sample collection are described in a previous study (Wang et al., 2012). In brief, three sediment cores with a diameter of 6 cm were retrieved using a gravity core at each site along a water-depth gradient of 0–90 m. The surface sediments were then pooled together. All samples were freeze-dried using a vacuum freeze-dryer and stored at −20°C. Surface water samples were collected from the top 0.5 m of the water column corresponding to the surface sediments collection. Bottom water samples were collected from the sediment–water interface.

Characterization of environmental factors

We collected and measured a variety of physicochemical parameters of surface, bottom water and surface sediment to investigate the environmental factors influencing microbial communities and functional genes. For surface water, we measure temperature, pH, dissolved oxygen, conductivity, total nitrogen, total phosphorus, concentration of HCO3, chlorophyll-a, and silicon content (Wang et al., 2012). For bottom water, temperature, pH, dissolved oxygen, chlorophyll-a and conductivity were measured (Syarif Sukri et al., 2023). For surface sediment, we quantified water depth, total phosphorus, total nitrogen, loss on ignition, porosity, water content, 19 types of metal ions and particle size. The term “water depth” in this study specifically refers to the depth of surface sediments. Detailed methodologies for the measuring and calculating of these abiotic variables have been described in previous studies (Wang et al., 2007). To simplify the complexity of these metal ions, a principal component analysis (PCA) was performed. The first two principal components (PC1 and PC2) were extracted and incorporated as additional environmental parameters for further analysis (Wu et al., 2019; Zhang et al., 2024; Zhao et al., 2019). The details of environmental factors could be found in Supplementary Table S1.

DNA extraction and metagenome sequencing

Total DNA was extracted from approximately 0.4 g of freeze-dried sediment using the PowerSoil DNeasy Kit (QIAGEN, Germany), with the DNA quality for all samples was evaluated via UV spectrophotometry. Sequencing was performed using a 2 × 150 bp paired-end strategy on the Illumina NovaSeq6000 platform. Raw reads were subjected to quality control and adapter trimming using FastQC and Trimmomatic (Andrews et al., 2010). Reads with average Phred scores below 25 or shorter than 50 bp were discarded. These reads were assembled into contigs with MEGAHIT v1.2.9 using the ‘meta-sensitive’ mode, and the protein-coding genes were predicted by prodigal v2.6.3 from the contigs longer than 1,000 bp (Li et al., 2015). To reduce redundancy, the amino acid sequences of these genes were clustered at the global level using the Linclust algorithm in MMseqs2, with a minimum sequence identity of 0.3 and a fraction of aligned sequences of 0.5, respectively (Hyatt et al., 2010). The relative abundance of the clustered genes was quantified using Salmon v1.0.0, followed by the normalization of total sequencing reads for each gene (Supplementary Table S2). Finally, functional annotation was performed by mapping the clustered genes to the eggNOG 5.0 database with eggnog-mapper v2.1.11, using the DIAMOND v2.1.8 as the search engine (Cantalapiedra et al., 2021; Huerta-Cepas et al., 2019). The genes assigned to the KEGG Orthology (KOs) in the annotation results were aggregated and subjected to the downstream analyses (Andrews et al., 2010). The statistics of assembly and the KO-representing genes across samples were provided in the Supplementary Table S3.

Community analyses using the rpS3-based approach

Microbial community profiling was performed using a modified pipeline of ribosomal protein rpS3 as described previously (Diamond et al., 2019; Ren et al., 2024). Firstly, the prokaryotic species among the community was represented by the clustering of the conserved marker genes rpS3 genes from metagenomic assemblies. The marker rpS3 genes were initially identified from assembled contigs using hmmsearch v3.2.1 (Eddy, 2011) against a custom HMM database, followed by clustering at 99% sequence similarity with USEARCH v11.0.667 (Edgar, 2010) to define species-level operational taxonomic units. Secondly, the relative abundance of the rpS3-represented species was estimated through a metagenomic read mapping strategy. The clean reads were aligned to the longest rpS3-containing contigs using Bowtie2 v2.3.5 (Li et al., 2023). The hit reads with ≥99% similarity were counted using the ‘depth’ module of Samtools v1.15.1 (Li et al., 2009). The final abundance of each species in a sample was calculated as the total mapped bases on the representative sequence divided by the representative sequence length and normalized by the total sequencing bases in the sample. Thirdly, taxonomic classification was achieved using the sequence alignment followed by phylogenetic correction. A custom rpS3 reference database was constructed from the RefSeq prokaryotic genome collection (~27,000 genomes, July 2019). Initial taxonomic assignments were generated through BLASTP analysis with the e-value threshold 1 × 10−3, minimum 50% sequence identity against the reference database. Concurrently, a phylogenetic tree was established through multiple sequence alignment of rpS3 amino acid sequences using MAFFT v7.427 (Katoh et al., 2005), subsequent alignment refinement using trimAl v1.4.1 with automated parameters (Capella-Gutiérrez et al., 2009), and construction of an approximate maximum-likelihood tree using FastTree v2.1.11 (Price et al., 2010). Final taxonomic designations were determined through reconciliation of both complementary approaches, with phylogenetic evidence superseding BLAST-based assignments in cases of discordance or absence of significant database matches.

Statistical analysis

We conducted a series of statistical analyses to investigate the distribution of microbial communities and functional genes in the surface sediment across water depths. We quantified the alpha diversity of microbial communities and functional genes using Shannon diversity indices (Dixon, 2003) based on the rpS3-based species abundance table and the KEGG Orthology (KO) abundance table at gene level, respectively. To further explore the relationship between Shannon diversity of functional genes and water depth, we applied both linear and quadratic models, selecting the best model based on the lowest Akaike information criterion (AIC) value (Supplementary Table S4) (Spellerberg and Fedor, 2003). For spatial variability (beta diversity), we assessed community composition and functional gene diversity using Bray–Curtis dissimilarity (Gutiérrez-Cánovas et al., 2013). Water depth distances were quantified using Euclidean metrics. To examine the relationship between functional gene beta diversity and depth-related distances, we utilized a Gaussian generalized linear model to analyze distance decay patterns (Supplementary Table S5) (Morlon et al., 2008). The statistical significance of these patterns was evaluated using a Mantel test with 9,999 permutations. Additionally, we performed Non-metric Multidimensional Scaling (NMDS) based on Bray–Curtis dissimilarities to visualize the functional gene composition along the depth gradient.

To explore key drivers of microbial community structure and functional gene distribution, we first applied variable clustering to assess the collinearity among environmental variables, removing redundant variables (Spearman’s ρ2 > 0.7) to minimize potential confounding effects (Supplementary Figures S2, S3). The relationships between other environmental factors and water depth were presented in supplementary materials (Supplementary Figure S4). In addition, random forest model was employed to determine the optimal number of trees (2,000) using cross-validation and then assessed the importance of each variable (Liaw, 2002). The importance scores were normalized to reflect their relative contribution to the overall model, and we iterated this process, eliminating the least influential variables, until each remaining variable contributed more than 5% to the model (Zhao et al., 2019). Furthermore, we assessed the correlation between microbial community structure, functional gene composition, and environmental variables using Mantel tests with non-redundant variables (Aljohani et al., 2024). In this analysis, Bray–Curtis dissimilarities were used to quantify microbial composition, while Euclidean distances were used for environmental variables.

To further explore the ecological processes shaping microbial communities, we applied a phylogeny-based null model (iCAMP) framework (Ning et al., 2020). This framework utilizes the beta-net correlation index and Raup-Crick index to estimate phylogenetic beta diversity and taxonomic beta diversity (Ning et al., 2020; Stegen et al., 2013). We quantified the relative importance of five ecological processes, including heterogeneous selection, homogeneous selection, dispersal limitation, homogenizing dispersal, and drift. Their differences across water layers were evaluated using the Wilcoxon test (Divine et al., 2013).

To identify microbial functional profiles, we annotated metagenomic reads by comparing them with the Kyoto Encyclopedia of Genes and Genomes (KEGG) database.1 To better understand the functional profiles across the three water depth layers, we first identified the functional genes related to carbon, nitrogen, and sulfur cycling (Supplementary Table S3). A generalized linear model (GLM) with a negative binomial distribution was used to estimate differences in gene enrichment between the layers (Robinson et al., 2010). Statistical significance was determined using the Benjamini–Hochberg false discovery rate, with a threshold of p < 0.05. GLM modeling was performed using the glmFit function in the edgeR package (Chen et al., 2014). Additionally, to assess the importance of metabolic pathways in community-level metabolic potential and the functional contributions of microbes, we calculated the explained variance of each functional pathway (Supplementary Table S6) and the contribution of microbial genus to these pathways (Supplementary Table S7) (Anantharaman et al., 2016). All p-values generated from the Wilcoxon tests were adjusted using the Benjamini–Hochberg false discovery rate method (Chung et al., 2007). All statistical analyses were conducted using R version V4.3.1, with the packages such as vegan V2.6.4, randomForest V4.7.1.1, Hmisc V5.2-3, edgeR V3.30.3, dplyr V1.1.3, ggplot2 V3.4.3 and ggcor V0.9.8.

Results

The alpha and beta diversity of taxonomy and functional genes across layers

The compositions of taxonomy and functional genes in Lugu Lake were associated with the depth of the overlaying water column, with more similar profiles toward deeper layers (Figures 1af; Supplementary Figures S5, S6). We further evaluated the Shannon diversity index and β-diversity metrics such as Bray–Curtis dissimilarity of taxa (Figures 1b,c) and functional genes (Figures 1e,f) along depth. Our results showed that microbial Shannon diversity exhibited a hump-shaped pattern peaking at the thermocline (Figure 1b), while functional gene diversity showed a U-shaped pattern along depth (Figure 1e). In addition, our study revealed the three functional subgroups involved in C, N, and S cycling showed distinct patterns (Figures 2ac and Supplementary Table S4). Specifically, nitrogen cycling genes showed a significant hump-shaped trend along depth peaking around 50 m (R2adj = 0.68, p < 0.001; Figure 2b). In contrast, sulfur cycling genes decreased along with depth, showing a sharp decline beyond 50 m (R2adj = 0.56, p < 0.001; Figure 2c).

Figure 1
Graphs displaying taxonomic and functional gene data. (a) NMDS plot for taxonomy with groups SUR, THE, DEE showing different clustering and stress value 0.046. (b) Taxonomy shows Shannon diversity against water depth, noting an R² of 0.29, P<0.001. (c) Taxonomy Bray-Curtis dissimilarity violin plots for SUR, THE, DEE with statistical significance indicated. (d) NMDS plot for functional genes with groups SUR, THE, DEE and stress value 0.077. (e) Functional gene Shannon diversity versus water depth, with R² of 0.46, P<0.001. (f) Functional gene Bray-Curtis dissimilarity violin plots showing significant differences among groups.

Figure 1. Comparisons of microbial and functional diversity and composition among different water layers. Non-Metric Multidimensional Scaling (NMDS) plots of (a) taxonomy and (d) functional genes. Each point represents a sample, which was colored by water depth, from surface layer (SUR, 0–10 m) to thermocline (THE, 10–50 m) and then to deep layer (DEE, 50–90 m). The Shannon diversity of (b) taxonomy and (e) functional genes among three water layers. Differences in microbial beta diversity consisting of (c) taxonomic and (f) functional genes variation (determined by pair Bray–Curtis distance) among three water layers. Different asterisks in the violin plots denote significant differences in corresponding variables between layer (determined by a two-sided pairwise Wilcoxon test). *p < 0.05, **p < 0.01, and ***p < 0.001, and ns: non-significant. In boxplots, the lower and upper hinges of the box correspond to the first and third quartiles (the 25th and 75th percentiles); the upper and lower whiskers extend from the hinge to the largest and smallest values no further than 1.5 times the interquartile range (IQR), respectively; and the central lines represent the median.

Figure 2
Six scatter plots show relationships in cycling for carbon, nitrogen, and sulfur. Plots (a), (b), and (c) depict Shannon diversity against water depth, showing no correlation, a quadratic relationship, and a negative correlation, respectively. Adjusted R-squared values are 0.01, 0.68, and 0.56. Plots (d), (e), and (f) illustrate Bray-Curtis dissimilarity against differences in water depth, each with positive correlations. Correlation coefficients are 0.39, 0.38, and 0.57, all with significant p-values below 0.001. Different colors represent each cycling type.

Figure 2. Water-depth diversity patterns and distance-decay relationship for functional genes. We considered the Shannon diversity of the three subgroups of functional genes involved in carbon cycling (a), nitrogen cycling (b), and sulfur cycling (c) (Supplementary Table S3). The relationships between functional gene diversity and water depth were evaluated by linear and quadratic models. The better model was selected based on the lower value of the Akaike information criterion. The lower panels (d–f) show the relationships between water depth changes and Bray–Curtis dissimilarity of the three subgroups. Linear regressions of relationships based on a linear model are shown with a solid line. Mantel tests were used to examine correlations between differences in functional gene composition and differences in community composition using 9,999 permutations. The Mantel r-values are shown, with all p-values being less than 0.001. The term “water depth” in this study specifically refers to depth of surface sediments.

For beta diversity, the spatial variability of both microbes and functional genes significantly decreased with the depth of the overlaying water column in the deep lake (Figures 1c,f). This indicates distinct changes in the composition of microbial communities and functional genes across different depths and locations. This is supported by the fact that the compositions of microbial communities (Mantel r = 0.57, p < 0.001) and functional genes (r = 0.37, p < 0.001) showed a significant distance-decay relationship with water depth changes (Figures 2df; Supplementary Figures S7a,b and Supplementary Table S5). The relative abundance of different functional genes remained balanced across samples (Figures 1f, 3b,d; Supplementary Figure S8), compared with microbial communities (Figures 1c, 3c; Supplementary Figure S8).

Figure 3
(a) Heatmap displaying microbial genera abundance across three conditions: SUR, THE, DEE, with a color gradient indicating abundance. (b) Heatmap showing metabolic pathway activity under the same conditions with a color scale. (c) Bar chart of microbial genus relative abundance at different sample points, with depth indicated by color shades. (d) Bar chart illustrating relative abundance of metabolic pathways with varying depths.

Figure 3. The composition of taxonomic groups and functional traits across different water depths. (a) Microbial community and (b) functional gene profiles, with samples ordered by water depth. SUR, surface layer; THE, thermocline; DEE, deep layer. Darker colors correspond to higher relative abundances. The relative abundance of (c) microbial genus and (d) functional pathways in samples from different water depths. Only the top 15 microbial genera with high relative abundances are annotated in the figure.

Environmental drivers and ecological processes shaping microbial assemblages

We investigated the relative importance of environmental factors in shaping the diversity and composition of microbial communities and functional genes (Figure 4a; Supplementary Figure S9). Water depth had the highest contribution of 17.1% for taxonomic diversity and 26.4% for functional gene diversity. Microbial and functional gene compositions were mainly driven by similar environmental factors, including water depth with 29.9% for microbial composition and 26.5% for functional gene composition, followed by surface temperature and bottom conductivity (Figure 4a; Supplementary Figure S10). Depth was the most influential factor among all environmental factors, distinguishing samples from different depth layers in both taxonomic and functional compositions. Specifically, microbial community structures and functional gene composition in the surface layer were most strongly linked to sediment phosphorus, followed by porosity (Figure 4b). In the thermocline, they were significantly associated with depth, followed by bottom conductivity (Figure 4b). In the deep layer, dissolved oxygen and surface water pH showed the strongest correlations with them (Figure 4b).

Figure 4
(a) Bar charts display key environmental variables influencing taxonomic and functional gene diversity. (b) Network diagram shows interactions between environmental factors with color coding by correlation strength. (c) Violin plots illustrate relative importance of ecological processes for different sites. (d) Additional violin plots compare process importance across layers, highlighting dispersal limitation and other ecological processes.

Figure 4. Environmental factors and ecological processes shaping microbial community structure and functional genes. (a) Relative contribution of environmental factors to taxonomic and functional diversity. Random forest analysis identified and quantified significant predictors of Shannon diversity and composition. The first axis of NMDS was used to represent composition. We selected the explanatory variables with a relative contribution rate >5%. Details of variable abbreviations are provided in Supplementary Table S1. (b) Associations between microbial community structure and functional gene composition (determined by Bray–Curtis distance) with environmental factors (determined by Euclidean distance) using the partial Mantel test. Partial Mantel’s r values are indicated by the edge width, while the statistical significance is denoted by the edge color. Pairwise correlations of environmental variables are depicted with a color gradient reflecting Spearman’s correlation coefficient. (c) The relative contribution of each ecological process driving microbial community assembly within the layer based on null model analysis (n = 231). (d) Differences in the relative importance of ecological processes among three water layers (n = 231). Different lowercase letters in box plots indicate significant differences for the ecological processes with soil depth (determined by a two-sided Wilcoxon test, p < 0.05). SUR, surface layer; THE, thermocline layer; DEE, deep layer. Central line and whiskers in each box represent the median and 1.5 times the interquartile range, respectively. Boxes indicate the interquartile range between 25th and 75th percentiles. Single points are outliers.

Although environmental factors significantly influenced microbial community structure, stochastic ecological processes dominated microbial community assembly in the deep lake as revealed through null model analysis (Figures 4c,d; Supplementary Figure S11). The relative importance of stochastic processes markedly varied with water depth. Specifically, dispersal limitation with an average contribution of 43.7% was the most critical factor in shaping microbial assemblages in the surface layer, followed by homogeneous selection with 29.5% (Figure 4c). In the thermocline and deep layers, dispersal limitation remained the primary mechanism influencing microbial community assembly with 51.6 and 47.1%, followed by drift with 22.4 and 30.4%, respectively (Figure 4c).

Enrichment of layer-specific C, N, and S cycling gene

Our results revealed significant enrichment in the composition of functional genes involved in carbon, nitrogen, and sulfur cycles across different layers (Figures 5ac; Supplementary Figure S12 and Supplementary Table S3). Genes enriched in the surface layer were predominantly associated with rapid energy acquisition and organic matter degradation. We found that the genes involved in the tricarboxylic acid (TCA) cycle, glycolysis, and methanogenesis were enriched in the surface layer (Figure 5a). The abundance of methanogenesis-related genes (mcrA, mtrB) was higher in the surface layer (Figure 5a). Additionally, genes involved in nitrate assimilation (nrtABC) and organic nitrogen metabolism (CPS1) were more abundant in the surface layer, suggesting that surface microbes preferentially utilize inorganic nitrogen and organic nitrogen sources (Figure 5b). Genes related to thiosulfate assimilation were also enriched in the surface layer, suggesting that surface microbes assimilate sulfur to support higher primary productivity and organic sulfur demand (Figure 5c). Genes involved in anaerobic metabolic pathways or reduction reactions were more abundant in the deeper layers. The thermocline and deep layer were enriched in genes associated with nitrification (amoA, hao), denitrification (nirS, nirK), and dissimilatory nitrate reduction (nafH, narG) (Figure 5b). Genes related to sulfide oxidation (soxACDXY) were enriched in the thermocline. Additionally, phthalate degradation genes were more abundant in the thermocline. The deep layer was enriched in genes associated with dissimilatory sulfate reduction (aprB, sat).

Figure 5
Diagram depicting metabolic pathways and gene expression heatmaps in three panels, labeled a, b, and c. Panel a shows carbon metabolism pathways, with glycolysis and degradation pathways labeled. Panel b illustrates nitrogen metabolism, denitrification, nitrification, and amino acid synthesis. Panel c details sulfur metabolism, sulfate assimilation, and oxidation pathways. Heatmaps in each panel represent log fold changes in gene expression across three conditions: SUR, THE, and DEE, with a color gradient from red (positive) to blue (negative).

Figure 5. Differences in the abundance of functional genes involved in C, N, and S cycling across three layers. The heatmap shows the enrichment of functional genes involved in (a) carbon cycling, (b) nitrogen cycling, and (c) sulfur cycling among three water layers. Statistical significance of the changes in gene abundance was assessed by a generalized linear model with a negative binomial distribution using edgeR package. The p-values were obtained from two-sided likelihood ratio tests (LRTs) and adjusted for multiple comparisons via the Benjamini–Hochberg procedure. Genes with significant changes in abundance (p < 0.05) are indicated with an asterisk. LogFC, log2-fold change. The full names of the genes in this figure are listed in Supplementary Table S3. SUR, surface layer; THE, thermocline layer; DEE, deep layer.

Contribution of microbial taxa to C, N, and S cycling genes across the layer

To evaluate the relative importance of metabolic pathways in the deep lake and the contribution of microbial taxa to metabolism, we quantified the explained variance of each metabolic pathway and assessed the relative contributions of various microbial phyla (Figure 6a and Supplementary Tables S6, S7). Pathways with higher explained variance were considered more critical to community-level metabolism. Among the metabolic pathways, nitrification emerged as the most critical process with explained variation of 91.07%, followed by denitrification with 89.56%, methanogenesis with 81.72%, and dissimilatory nitrate reduction with 77.45% (Figure 6a and Supplementary Table S6). This demonstrates the important roles of carbon, nitrogen, and sulfur cycles in regulating microbial energy metabolism in the lake. In the surface layer, metabolic activity was predominantly driven by Nitrospira, Gloeomargarita and Methanothrix, (Figure 6b and Supplementary Table S7). The thermocline exhibited more diverse microbial communities, with significant contributions from Nitrospira, Haliea, Methanothrix, Pedosphaera and Anaeromyxobacter (Figure 6b and Supplementary Table S7). In the deep layer, Aromatoleum, Rubrivivax and Candidatus Methylomirabilis were the primary contributors to functional cycling (Figure 6b and Supplementary Table S7).

Figure 6
A two-part figure showing microbial data. Part (a) is a heatmap with explained variation percentages and microbial taxa. Columns represent processes, and rows represent taxa, with colors indicating the strength and direction of relationships. Circle sizes denote importance percentages. Part (b) is a Sankey diagram linking study areas (SUR, THE, DEE) to microbial taxa and processes, demonstrating flow and relationships between them. Different colors represent different pathways.

Figure 6. Contribution of microbial communities to biogeochemical processes across different water layers. (a) The contribution of microbial genus to each metabolic pathway through random forest analysis is represented by circles of different sizes. The color gradient reflects the strength of the Spearman correlation coefficient, where dark blue indicates a strong positive correlation, and dark red represents a strong negative correlation. Statistical significance is denoted by asterisks: ***p < 0.001, **p < 0.01, and *p < 0.05. (b) The Sankey diagram illustrates the difference in the contributions of microbial groups to specific biogeochemical processes across three water layers, with the taxonomic classification of microbial groups and their associated category of functional pathways. The three columns represent, from left to right, water layers, taxonomic groups, and metabolic pathways, respectively. SUR, surface layer; THE, thermocline layer; DEE, deep layer.

Discussion

Using metagenome sequencing technology, we identified four main findings: (i) The diversity of microbial communities and functional genes exhibited distinct depth patterns, with microbial Shannon diversity exhibiting a hump-shaped pattern and functional gene diversity showing a U-shaped pattern; (ii) Despite water depth significantly driving microbial and functional gene composition, stochastic ecological processes such as dispersal limitation dominated community assembly in the deep lake; (iii) Functional gene enrichment analysis revealed that the surface layer was primarily linked to aerobic respiration and methanogenesis, while the thermocline and deep layers were enriched in genes associated with anaerobic metabolism; (iv) The metabolism at the surface layer was contributed by taxa that thrive on light and oxygen for rapid energy acquisition and decomposition, while deeper layers shifted to anaerobic metabolism.

Taxonomic and functional diversity across layers

The thermocline effectively separated the mixed or surface waters from the colder bottom layers (Chang et al., 2022). Oxygen levels fluctuate in the thermocline region of Lugu Lake, and there are steep gradients of temperature and redox potential (Ren and Wang, 2022; Tran et al., 2021). Frequent sediment resuspension events potentially altering redox gradients in the surface sediments (Dadi et al., 2017; Nevers et al., 2020). This unique environment creates diverse microhabitats, providing opportunities for different microbial communities to thrive (Xing et al., 2019). In the deep layer, the challenging conditions of sediments such as low temperatures, hypoxia, and nutrient scarcity restricted microbial growth, potentially leading to a decrease in microbial Shannon diversity (Cabello-Yeves et al., 2019). In the thermocline region, steep oxygen gradients made nitrogen cycling genes diverse and active (Chang et al., 2022). The dynamic environment of this layer facilitated the coexistence and functional differentiation of nitrogen-cycling microbes, leading to peak functional gene diversity (Martens-Habbena et al., 2009). In the deep and oxygen-deprived layer, sulfur cycling shifted toward sulfate reduction, a process typically dominated by specific taxa like sulfate-reducing bacteria (Chang et al., 2022; Zhang et al., 2024). The harsh conditions imposed strong selective pressure, allowing well-adapted microbes carrying specific functional genes (e.g., dsrAB, aprB, sat) to persist, thus limiting sulfur cycling functional gene diversity (Müller et al., 2015; Wasmund et al., 2017).

Microbial communities of sediments in the surface layer exhibited higher spatial heterogeneity, while those in the deeper layers tended to be more homogenized. This trend aligned with previous research showing that the total beta diversity of microbial communities in the sediments showed significant decreasing trends toward deeper water in deep lake sediments (Wu et al., 2020). The homogenization in deep layers can be attributed to the unique physicochemical characteristics of deep lakes, such as depth-specific environmental gradients and limited nutrient availability, which restricts species exchange and diffusion between microbial communities in the water columns at different depths (Lear et al., 2014). Microbial communities in the surface layers are more influenced by nutrient exchange, fostering a more varied ecological environment (Wang et al., 2022).

Environmental drivers and ecological processes shaping microbial assemblages

The responses of microbial communities and functional genes to environmental factors were strongly coupled, contrasting with previous studies in soil and ocean where environmental factors strongly influenced functional groups but had a weaker effect on community composition (Louca et al., 2016). Depth and temperature were identified as robust predictors of microbial communities and functional genes of the sediments in lake ecosystems (Zhang et al., 2024; Zhao et al., 2019). As a multifaceted proxy, depth encapsulated variations in temperature, dissolved oxygen, light availability, and nutrient levels, influencing microbes directly or indirectly (Kitazawa et al., 2018; Nevalainen, 2012; Rojas-Jimenez et al., 2021). These environmental gradients modulated resource availability, energy-related processes, and stressor prevalence, exerting selective pressure on microbial communities and functional composition (Gurung et al., 2001; Wu et al., 2020).

Despite the relatively open nature of the surface environment, spatial resource heterogeneity imposed significant limitations on microbial dispersal (Martiny et al., 2006). Homogeneous selection was likely linked to the strong selective pressures exerted by key environmental factors such as light availability and dissolved oxygen concentration (Gotelli, 2000; Liu et al., 2020; Zhao et al., 2019). The thermocline, characterized by the steepest temperature changes within the water column, may constrain microbial migration and dispersal by altering water density and viscosity (Morrison et al., 2017). This pronounced these temperature changes created a strong physical stratification, impeding material exchange between upper and lower water layers and restricting microbial migration across different depths (Peura et al., 2015; Zhang et al., 2024). The deep layer was characterized by stable but harsh physicochemical conditions, limiting the ability of different microbial species to colonize and establish diverse communities (Wang et al., 2013). Furthermore, the importance of drift increased significantly with depth (Figure 4d), likely due to the scarcity of resources and smaller community sizes in deeper layers (Ofiţeru et al., 2010). These factors amplified the impact of random population fluctuations on community structure. We also found homogenizing selection decreased with depth (Figure 4d), reflecting the diminishing influence of environmental factors such as light and oxygen (Wang et al., 2013; Zhao et al., 2021). The effects of environmental filtering weakened, and stochastic processes became more dominant in shaping community structure as depth increases.

Functional enrichment of C, N, and S cycling genes across layers

The enrichment of genes associated with energy acquisition and organic matter degradation in the surface layer suggests that microbial communities in this zone are adapted to the oxygen-rich, light-rich conditions where rapid energy cycling is critical for supporting primary productivity (Cunliffe et al., 2008; Peura et al., 2018). Moreover, the high abundance of methanogenesis-related genes in the surface layer may be attributed to the presence of localized anaerobic microenvironments, such as particle deposition zones, which facilitate methane production processes (Bridgham et al., 2013). The preferential utilization of inorganic and organic nitrogen sources by surface microbes highlights the competitive pressures in this layer, where primary producers like phytoplankton dominate and microbial communities must efficiently assimilate available nitrogen to sustain their growth (Chrost et al., 2009). Similarly, the enrichment of sulfur-related genes suggests that surface microbes also play a role in supporting primary productivity by assimilating sulfur, an essential nutrient for microbial growth and metabolic processes.

The enrichment of genes involved in anaerobic metabolic processes in the deeper layers suggests a shift toward energy pathways that do not rely on oxygen, as microbes adapt to low-oxygen conditions in these layers. The presence of genes related to nitrification, denitrification, and nitrate reduction in the thermocline and deep layers indicates that microbes in these layers rely on anaerobic nitrogen transformations, such as nitrate and nitrite reduction, to adapt to the limited oxygen availability in deeper waters (Robert Hamersley et al., 2009). The enrichment of genes related to sulfide oxidation in the thermocline was probably due to the elevated sulfide concentrations at the redox interface (Chen et al., 2024). Microbes in this layer may oxidize sulfide to generate intermediate sulfur compounds, to adapt to the oxidative conditions of this transitional environment (Zhang et al., 2021). Phthalate degradation genes were more abundant in the thermocline. Phthalates and their derivatives, common pollutants, may originate from sediment release or external inputs (Tuan Tran et al., 2022). The sharp temperature gradient and stable chemical stratification of the thermocline created a unique ecological niche supporting the degradation of complex organic compounds (Kurt, 2019). Microbes in this layer may leverage adaptive metabolic capabilities to selectively degrade structurally complex organic carbon, contributing significantly to the carbon cycle in this stratified environment. The enrichment of genes associated with dissimilatory sulfate reduction in deep layer reflected microbial reliance on sulfate as a terminal electron acceptor for anaerobic respiration (Zhu et al., 2018). This metabolic trait corresponded to the anoxic conditions of the deep layer, where sulfate reduction serves as a critical energy-yielding process (Watanabe et al., 2013).

Layer-specific contribution of microbial taxa to functional genes

The surface layer, characterized by high organic carbon and nutrient availability, supports robust microbial activity (Ertefai et al., 2008). Gloeomargarita, an oxygenic photoautotroph, contributed to both photosynthesis and organic carbon degradation (Bacchetta et al., 2022; Moreira et al., 2017). Its distribution aligns with previous observations in soil (Nelson et al., 2016) and marine systems (Zehr and Capone, 2020), suggesting that Gloeomargarita distribution is strongly influenced by factors such as light intensity, temperature, and nutrient availability. Its predominance highlights the importance of oxygenic metabolism and primary production in surface waters, where light and oxygen are readily available.

Nitrospira, as a key nitrifying genus, indicates high nitrogen turnover in the thermocline, suggesting active nitrification processes (Ren and Wang, 2022; Winter et al., 2009). The thermocline with its steep temperature gradient and variable oxygen levels, supports both aerobic and anaerobic metabolic processes (Gorham and Boyce, 1989). Haliea commonly exhibited photoheterotrophic traits, indicating the persistence of light-driven energy capture in the thermocline (Yamamoto et al., 2020; Yang et al., 2020). These metabolic activities likely provide additional energy for microbial processes in the thermocline. The deep layer, characterized by low oxygen and nutrient availability, exhibited significant anaerobic and heterotrophic metabolic activities. Aromatoleum remained dominant in the deep layer, utilizing its ability to metabolize organic substrates under anaerobic conditions (Becker et al., 2022). Its metabolic strategies likely reflect adaptations to limited energy resources, emphasizing the importance of efficient resource utilization in deep lake environments (Vagts et al., 2018).

Conclusion

Our study provides metagenomic analyses of microbial community structure and functional potential along a water depth gradient in Lugu Lake, a deep lake in subtropic zone. We revealed significant shifts in microbial diversity and functional gene composition across the lake sediment in response to the depth of the overlaying water column and its redox state, revealing tight metabolic coupling between sediment and water column ecosystems. Microbial communities in the shallower layers exhibited higher spatial heterogeneity, while those in the deep layers were more homogenized. The thermocline, with steep gradients of temperature and redox potential, created diverse microhabitats that support different microbial lineages to thrive. We also observed that microbial communities involved in alternative electron accepting processes were more diverse in the thermocline, likely due to the lower redox potential and complex nutrient strategies in this layer.

Functional genes involved in carbon, nitrogen, and sulfur cycling showed layer-specific enrichment. Surface waters were dominated by taxa that exploit abundant light and oxygen, favoring rapid energy acquisition and organic matter decomposition. In contrast, the thermocline and deep layers shifted toward anaerobic and specialized metabolic pathways, reflecting adaptations to oxygen-limited and nutrient-poor conditions. The diverse metabolic strategies observed across depth gradients underscore the critical role of microbial communities in regulating biogeochemical cycles in deep lakes. These findings emphasize that the thermocline significantly affects the shaping of microbial community and functional gene distributions. Future research could focus on multi-season sampling to capture annual variability in microbial dynamics and ecosystem processes.

Data availability statement

The metagenomic data presented in this study have been deposited in the National Omics Data Encyclopedia (NODE, https://www.biosino.org/node/index) under the accession number OEX011297 (experiment ID).

Author contributions

PZ: Data curation, Writing – review & editing, Formal analysis, Writing – original draft. MR: Methodology, Data curation, Resources, Writing – review & editing. YX: Writing – review & editing. JW: Conceptualization, Writing – review & editing, Formal analysis.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was financially supported by the National Natural Science Foundation of China (U24A20578, 42372353, 42225708, 92251304, 92351303, 42002304) and the Science and Technology Planning Project of NIGLAS (NIGLAS2022GS09).

Conflict of interest

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

The reviewer DH declared a shared parent affiliation with the authors MR and JW to the handling editor at the time of review.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

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.2025.1614055/full#supplementary-material

Footnotes

References

Aljohani, A., Clarke, D., Byrne, M., and Fleming, G. (2024). The bacterial microbiome and resistome of house dust mites in Irish homes. Sci. Rep. 14:19621. doi: 10.1038/s41598-024-70686-y

PubMed Abstract | Crossref Full Text | Google Scholar

Anantharaman, K., Brown, C. T., Hug, L. A., Sharon, I., Castelle, C. J., Probst, A. J., et al. (2016). Thousands of microbial genomes shed light on interconnected biogeochemical processes in an aquifer system. Nat. Commun. 7:13219. doi: 10.1038/ncomms13219

Crossref Full Text | Google Scholar

Andrews, S., Krueger, F., Segonds-Pichon, A., Biggins, L., Krueger, C., and Wingett, S. (2010). FastQC. A quality control tool for high throughput sequence data. 370, Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc. (Accessed July 22, 2025).

Google Scholar

Ayala-Muñoz, D., Macalady, J. L., Sánchez-España, J., Falagán, C., Couradeau, E., and Burgos, W. D. (2022). Microbial carbon, sulfur, iron, and nitrogen cycling linked to the potential remediation of a meromictic acidic pit lake. ISME J. 16, 2666–2679. doi: 10.1038/s41396-022-01320-w

PubMed Abstract | Crossref Full Text | Google Scholar

Bacchetta, T., López-García, P., Gutiérrez-Preciado, A., Mehta, N., Skouri-Panet, F., Benzerara, K., et al. (2022). Description of Gloeomargarita ahousahtiae sp. nov., a thermophilic member of the order Gloeomargaritales with intracellular carbonate inclusions. Eur. J. Phycol. 59, 11–20. doi: 10.1080/09670262.2023.2216257

Crossref Full Text | Google Scholar

Becker, P., Kirstein, S., Wünsch, D., Koblitz, J., Buschen, R., Wöhlbrand, L., et al. (2022). Systems biology of aromatic compound catabolism in facultative anaerobic Aromatoleum aromaticum EbN1T. mSystems 7:e0068522. doi: 10.1128/msystems.00685-22

Crossref Full Text | Google Scholar

Bridgham, S. D., Cadillo-Quiroz, H., Keller, J. K., and Zhuang, Q. (2013). Methane emissions from wetlands: biogeochemical, microbial, and modeling perspectives from local to global scales. Glob. Change Biol. 19, 1325–1346. doi: 10.1111/gcb.12131

PubMed Abstract | Crossref Full Text | Google Scholar

Bush, T., Diao, M., Allen, R. J., Sinnige, R., Muyzer, G., and Huisman, J. (2017). Oxic-anoxic regime shifts mediated by feedbacks between biogeochemical processes and microbial community dynamics. Nat. Commun. 8:789. doi: 10.1038/s41467-017-00912-x

PubMed Abstract | Crossref Full Text | Google Scholar

Cabello-Yeves, P. J., Picazo, A., Roda-Garcia, J. J., Rodriguez-Valera, F., and Camacho, A. (2023). Vertical niche occupation and potential metabolic interplay of microbial consortia in a deeply stratified meromictic model lake. Limnol. Oceanogr. 68, 2492–2511. doi: 10.1002/lno.12437

Crossref Full Text | Google Scholar

Cabello-Yeves, P. J., Zemskaya, T. I., Zakharenko, A. S., Sakirko, M. V., Ivanov, V. G., Ghai, R., et al. (2019). Microbiome of the deep Lake Baikal, a unique oxic bathypelagic habitat. Limnol. Oceanogr. 65, 1471–1488. doi: 10.1002/lno.11401

Crossref Full Text | Google Scholar

Cantalapiedra, C. P., Hernández-Plaza, A., Letunic, I., Bork, P., and Huerta-Cepas, J. (2021). eggNOG-mapper v2: functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Mol. Biol. Evol. 38, 5825–5829. doi: 10.1093/molbev/msab293

PubMed Abstract | Crossref Full Text | Google Scholar

Capella-Gutiérrez, S., Silla-Martínez, J. M., and Gabaldón, T. (2009). trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25, 1972–1973. doi: 10.1093/bioinformatics/btp348

PubMed Abstract | Crossref Full Text | Google Scholar

Chang, F., Hou, P., Wen, X., Duan, L., Zhang, Y., and Zhang, H. (2022). Seasonal stratification characteristics of vertical profiles and water quality of Lake Lugu in Southwest China. Water 14:2554. doi: 10.3390/w14162554

Crossref Full Text | Google Scholar

Chen, Y., McCarthy, D., Robinson, M., and Smyth, G. K. (2014) edgeR: differential expression analysis of digital gene expression data user’s guide. Bioconductor user’s guide. Available online at: https://www.bioconductor.org/packages/devel/bioc/vignettes/edgeR/inst/doc/edgeRUsersGuide.pdf (Accessed July 22, 2025).

Google Scholar

Chen, X., Sheng, Y., Wang, G., Zhou, P., Liao, F., Mao, H., et al. (2024). Spatiotemporal successions of N, S, C, Fe, and As cycling genes in groundwater of a wetland ecosystem: enhanced heterogeneity in wet season. Water Res. 251:121105. doi: 10.1016/j.watres.2024.121105

PubMed Abstract | Crossref Full Text | Google Scholar

Chrost, R. J., Adamczewski, T., Kalinowska, K., and Skowronska, A. (2009). Inorganic phosphorus and nitrogen modify composition and diversity of microbial communities in water of mesotrophic lake. Pol. J. Microbiol. 58, 77–90. Available at: https://pubmed.ncbi.nlm.nih.gov/19469290/

Google Scholar

Chung, P.-J., Bohme, J. F., Mecklenbrauker, C. F., and Hero, A. O. (2007). Detection of the number of signals using the Benjamini–Hochberg procedure. IEEE Trans. Signal Process. 55, 2497–2508. doi: 10.1109/TSP.2007.893749

Crossref Full Text | Google Scholar

Cunliffe, M., Schäfer, H., Harrison, E., Cleave, S., Upstill-Goddard, R., and Murrell, J. C. (2008). Phylogenetic and functional gene analysis of the bacterial and archaeal communities associated with the surface microlayer of an estuary. ISME J. 2, 776–789. doi: 10.1038/ismej.2008.28

PubMed Abstract | Crossref Full Text | Google Scholar

Dadi, T., Wendt-Potthoff, K., and Koschorreck, M. (2017). Sediment resuspension effects on dissolved organic carbon fluxes and microbial metabolic potentials in reservoirs. Aquat. Sci. 79, 749–764. doi: 10.1007/s00027-017-0533-4

Crossref Full Text | Google Scholar

Diamond, S., Andeer, P. F., Li, Z., Crits-Christoph, A., Burstein, D., Anantharaman, K., et al. (2019). Mediterranean grassland soil C–N compound turnover is dependent on rainfall and depth, and is mediated by genomically divergent microorganisms. Nat. Microbiol. 4, 1356–1367. doi: 10.1038/s41564-019-0449-y

PubMed Abstract | Crossref Full Text | Google Scholar

Divine, G., Norton, H. J., Hunt, R., and Dienemann, J. (2013). A review of analysis and sample size calculation considerations for Wilcoxon tests. Anesth. Analg. 117, 699–710. doi: 10.1213/ANE.0b013e31827f53d7

PubMed Abstract | Crossref Full Text | Google Scholar

Dixon, P. (2003). VEGAN, a package of R functions for community ecology. J. Veg. Sci. 14, 927–930. doi: 10.1111/j.1654-1103.2003.tb02228.x

Crossref Full Text | Google Scholar

Eddy, S. R. (2011). Accelerated profile HMM searches. PLoS Comput. Biol. 7:e1002195. doi: 10.1371/journal.pcbi.1002195

PubMed Abstract | Crossref Full Text | Google Scholar

Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461

PubMed Abstract | Crossref Full Text | Google Scholar

Ertefai, T. F., Fisher, M. C., Fredricks, H. F., Lipp, J. S., Pearson, A., Birgel, D., et al. (2008). Vertical distribution of microbial lipids and functional genes in chemically distinct layers of a highly polluted meromictic lake. Org. Geochem. 39, 1572–1588. doi: 10.1016/j.orggeochem.2008.07.009

Crossref Full Text | Google Scholar

Falkowski, P. G., Fenchel, T., and Delong, E. F. (2008). The microbial engines that drive earth’s biogeochemical cycles. Science 320, 1034–1039. doi: 10.1126/science.1153213

PubMed Abstract | Crossref Full Text | Google Scholar

Gorham, E., and Boyce, F. M. (1989). Influence of lake surface area and depth upon thermal stratification and the depth of the summer thermocline. J. Great Lakes Res. 15, 233–245. doi: 10.1016/S0380-1330(89)71479-9

Crossref Full Text | Google Scholar

Gotelli, N. J. (2000). Null model analysis of species co-occurrence patterns. Ecology 81, 2606–2621. doi: 10.1890/0012-9658(2000)081[2606:Nmaosc]2.0.Co;2

Crossref Full Text | Google Scholar

Gurung, T. B., Kagami, M., Yoshida, T., and Urabe, J. (2001). Relative importance of biotic and abiotic factors affecting bacterial abundance in Lake Biwa: an empirical analysis. Limnology 2, 19–28. doi: 10.1007/s102010170012

Crossref Full Text | Google Scholar

Gutiérrez-Cánovas, C., Millán, A., Velasco, J., Vaughan, I. P., and Ormerod, S. J. (2013). Contrasting effects of natural and anthropogenic stressors on beta diversity in river organisms. Glob. Ecol. Biogeogr. 22, 796–805. doi: 10.1111/geb.12060

Crossref Full Text | Google Scholar

Haro-Moreno, J. M., López-Pérez, M., De La Torre, J. R., Picazo, A., Camacho, A., and Rodriguez-Valera, F. (2018). Fine metagenomic profile of the Mediterranean stratified and mixed water columns revealed by assembly and recruitment. Microbiome 6:128. doi: 10.1186/s40168-018-0513-5

PubMed Abstract | Crossref Full Text | Google Scholar

Huerta-Cepas, J., Szklarczyk, D., Heller, D., Hernández-Plaza, A., Forslund, S. K., Cook, H., et al. (2019). eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 47, D309–D314. doi: 10.1093/nar/gky1085

PubMed Abstract | Crossref Full Text | Google Scholar

Hyatt, D., Chen, G.-L., LoCascio, P. F., Land, M. L., Larimer, F. W., and Hauser, L. J. (2010). Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 11:119. doi: 10.1186/1471-2105-11-119

PubMed Abstract | Crossref Full Text | Google Scholar

Imhoff, J. (2016). New dimensions in microbial ecology—functional genes in studies to unravel the biodiversity and role of functional microbial groups in the environment. Microorganisms 4:19. doi: 10.3390/microorganisms4020019

PubMed Abstract | Crossref Full Text | Google Scholar

Katoh, K., Kuma, K.-i., Toh, H., and Miyata, T. (2005). MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 33, 511–518. doi: 10.1093/nar/gki198

PubMed Abstract | Crossref Full Text | Google Scholar

Kitazawa, D., Yoshida, T., Zhou, J., and Park, S. (2018). Comparative study on vertical circulation in deep lakes: Lake Biwa and Lake Ikeda. 2018 OCEANS-MTS/IEEE Kobe Techno-Oceans (OTO)

Google Scholar

Kurt, H. (2019). “Depth distribution of microbial diversity in lakes” in Freshwater microbiology (London: Academic Press), 225–262.

Google Scholar

Lear, G., Bellamy, J., Case, B. S., Lee, J. E., and Buckley, H. L. (2014). Fine-scale spatial patterns in bacterial community composition and function within freshwater ponds. ISME J. 8, 1715–1726. doi: 10.1038/ismej.2014.21

PubMed Abstract | Crossref Full Text | Google Scholar

Li, X., Han, Q., Li, M., Luo, Q., Zhu, S., Zheng, Y., et al. (2023). Complete mitochondrial genome sequence, characteristics, and phylogenetic analysis of Oenanthe javanica. Agronomy 13:2103. doi: 10.3390/agronomy13082103

Crossref Full Text | Google Scholar

Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079. doi: 10.1093/bioinformatics/btp352

PubMed Abstract | Crossref Full Text | Google Scholar

Li, D., Liu, C.-M., Luo, R., Sadakane, K., and Lam, T.-W. (2015). MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 31, 1674–1676. doi: 10.1093/bioinformatics/btv033

PubMed Abstract | Crossref Full Text | Google Scholar

Liaw, A. (2002). Classification and regression by randomForest. R News 2, 18–22.

Google Scholar

Lipsewers, Y. A., Hopmans, E. C., Meysman, F. J., Sinninghe Damsté, J. S., and Villanueva, L. (2016). Abundance and diversity of denitrifying and anammox bacteria in seasonally hypoxic and sulfidic sediments of the saline lake Grevelingen. Front. Microbiol. 7:1661. doi: 10.3389/fmicb.2016.01661

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, Y., Chen, C., and Yang, S. (2019). Assessment of anthropogenic impact versus climate change on the succession of the diatom community in Lugu Lake (Yunnan-Guizhou Plateau, China) using the sedimentary record of geochemical elements. Water 11:655. doi: 10.3390/w11040655

Crossref Full Text | Google Scholar

Liu, Y., Ren, Z., Qu, X., Zhang, M., Yu, Y., Zhang, Y., et al. (2020). Microbial community structure and functional properties in permanently and seasonally flooded areas in Poyang Lake. Sci. Rep. 10:4819. doi: 10.1038/s41598-020-61569-z

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

Lv, H., Yang, M., Cheng, Y., Li, K., Ji, G., Huang, T., et al. (2024). Disentangling the assembly patterns and drivers of microbial communities during thermal stratification and mixed periods in a deep-water reservoir. Sci. Total Environ. 946:174398. doi: 10.1016/j.scitotenv.2024.174398

PubMed Abstract | Crossref Full Text | Google Scholar

Martens-Habbena, W., Berube, P. M., Urakawa, H., de la Torre, J. R., and Stahl, D. A. (2009). Ammonia oxidation kinetics determine niche separation of nitrifying Archaea and Bacteria. Nature 461, 976–979. doi: 10.1038/nature08465

PubMed Abstract | Crossref Full Text | Google Scholar

Martiny, J. B., Bohannan, B. J., Brown, J. H., Colwell, R. K., Fuhrman, J. A., Green, J. L., et al. (2006). Microbial biogeography: putting microorganisms on the map. Nat. Rev. Microbiol. 4, 102–112. doi: 10.1038/nrmicro1341

PubMed Abstract | Crossref Full Text | Google Scholar

Moreira, D., Tavera, R., Benzerara, K., Skouri-Panet, F., Couradeau, E., Gérard, E., et al. (2017). Description of Gloeomargarita lithophora gen. nov., sp. nov., a thylakoid-bearing, basal-branching cyanobacterium with intracellular carbonates, and proposal for Gloeomargaritales ord. nov. Int. J. Syst. Evol. Microbiol. 67, 653–658. doi: 10.1099/ijsem.0.001679

PubMed Abstract | Crossref Full Text | Google Scholar

Morlon, H., Chuyong, G., Condit, R., Hubbell, S., Kenfack, D., Thomas, D., et al. (2008). A general framework for the distance-decay of similarity in ecological communities. Ecol. Lett. 11, 904–917. doi: 10.1111/j.1461-0248.2008.01202.x

PubMed Abstract | Crossref Full Text | Google Scholar

Morrison, J. M., Baker, K. D., Zamor, R. M., Nikolai, S., Elshahed, M. S., and Youssef, N. H. (2017). Spatiotemporal analysis of microbial community dynamics during seasonal stratification events in a freshwater lake (Grand Lake, OK, USA). PLoS One 12:e0177488. doi: 10.1371/journal.pone.0177488

PubMed Abstract | Crossref Full Text | Google Scholar

Müller, A. L., Kjeldsen, K. U., Rattei, T., Pester, M., and Loy, A. (2015). Phylogenetic and environmental diversity of DsrAB-type dissimilatory (bi)sulfite reductases. ISME J. 9, 1152–1165. doi: 10.1038/ismej.2014.208

PubMed Abstract | Crossref Full Text | Google Scholar

Nelson, M. B., Martiny, A. C., and Martiny, J. B. (2016). Global biogeography of microbial nitrogen-cycling traits in soil. Proc. Natl. Acad. Sci. U.S.A. 113, 8033–8040. doi: 10.1073/pnas.1601070113

PubMed Abstract | Crossref Full Text | Google Scholar

Nevalainen, L. (2012). Distribution of benthic microcrustaceans along a water depth gradient in an Austrian Alpine lake—sedimentary evidence for niche separation. Limnologica 42, 65–71. doi: 10.1016/j.limno.2011.08.003

Crossref Full Text | Google Scholar

Nevers, M. B., Byappanahalli, M. N., Nakatsu, C. H., Kinzelman, J. L., Phanikumar, M. S., Shively, D. A., et al. (2020). Interaction of bacterial communities and indicators of water quality in shoreline sand, sediment, and water of Lake Michigan. Water Res. 178:115671. doi: 10.1016/j.watres.2020.115671

PubMed Abstract | Crossref Full Text | Google Scholar

Ning, D., Yuan, M., Wu, L., Zhang, Y., Guo, X., Zhou, X., et al. (2020). A quantitative framework reveals ecological drivers of grassland microbial community assembly in response to warming. Nat. Commun. 11:4717. doi: 10.1038/s41467-020-18560-z

PubMed Abstract | Crossref Full Text | Google Scholar

Ofiţeru, I. D., Lunn, M., Curtis, T. P., Wells, G. F., Criddle, C. S., Francis, C. A., et al. (2010). Combined niche and neutral effects in a microbial wastewater treatment community. Proc. Natl. Acad. Sci. U.S.A. 107, 15345–15350. doi: 10.1073/pnas.1000604107

PubMed Abstract | Crossref Full Text | Google Scholar

Peura, S., Buck, M., Aalto, S. L., Morales, S. E., Nykänen, H., and Eiler, A. (2018). Novel autotrophic organisms contribute significantly to the internal carbon cycling potential of a boreal lake. mBio 9:e00916-18. doi: 10.1128/mbio.00916-18

Crossref Full Text | Google Scholar

Peura, S., Sinclair, L., Bertilsson, S., and Eiler, A. (2015). Metagenomic insights into strategies of aerobic and anaerobic carbon and nitrogen transformation in boreal lakes. Sci. Rep. 5:12102. doi: 10.1038/srep12102

PubMed Abstract | Crossref Full Text | Google Scholar

Price, M. N., Dehal, P. S., and Arkin, A. P. (2010). FastTree 2—approximately maximum-likelihood trees for large alignments. PLoS One 5:e9490. doi: 10.1371/journal.pone.0009490

PubMed Abstract | Crossref Full Text | Google Scholar

Ren, M., Hu, A., Zhang, L., Yao, X., Zhao, Z., Kimirei, I. A., et al. (2024). Acidic proteomes are linked to microbial alkaline preference in African lakes. Water Res. 266:122393. doi: 10.1016/j.watres.2024.122393

PubMed Abstract | Crossref Full Text | Google Scholar

Ren, M., and Wang, J. (2022). Phylogenetic divergence and adaptation of Nitrososphaeria across lake depths and freshwater ecosystems. ISME J. 16, 1491–1501. doi: 10.1038/s41396-022-01199-7

PubMed Abstract | Crossref Full Text | Google Scholar

Robert Hamersley, M., Woebken, D., Boehrer, B., Schultze, M., Lavik, G., and Kuypers, M. M. M. (2009). Water column anammox and denitrification in a temperate permanently stratified lake (Lake Rassnitzer, Germany). Syst. Appl. Microbiol. 32, 571–582. doi: 10.1016/j.syapm.2009.07.009

PubMed Abstract | Crossref Full Text | Google Scholar

Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616

PubMed Abstract | Crossref Full Text | Google Scholar

Rojas-Jimenez, K., Araya-Lobo, A., Quesada-Perez, F., Akerman-Sanchez, J., Delgado-Duran, B., Ganzert, L., et al. (2021). Variation of bacterial communities along the vertical gradient in Lake Issyk Kul, Kyrgyzstan. Environ. Microbiol. Rep. 13, 337–347. doi: 10.1111/1758-2229.12935

PubMed Abstract | Crossref Full Text | Google Scholar

Song, W., Liu, J., Qin, W., Huang, J., Yu, X., Xu, M., et al. (2022). Functional traits resolve mechanisms governing the assembly and distribution of nitrogen-cycling microbial communities in the Global Ocean. mBio 13:e0383221. doi: 10.1128/mbio.03832-21

PubMed Abstract | Crossref Full Text | Google Scholar

Spellerberg, I. F., and Fedor, P. J. (2003). A tribute to Claude Shannon (1916–2001) and a plea for more rigorous use of species richness, species diversity and the ‘Shannon–Wiener’ Index. Glob. Ecol. Biogeogr. 12, 177–179. doi: 10.1046/j.1466-822X.2003.00015.x

Crossref Full Text | Google Scholar

Stegen, J. C., Lin, X., Fredrickson, J. K., Chen, X., Kennedy, D. W., Murray, C. J., et al. (2013). Quantifying community assembly processes and identifying features that impose them. ISME J. 7, 2069–2079. doi: 10.1038/ismej.2013.93

PubMed Abstract | Crossref Full Text | Google Scholar

Su, K., Wang, Q., Li, L., Cao, R., and Xi, Y. (2022). Water quality assessment of Lugu Lake based on Nemerow pollution index method. Sci. Rep. 12:13613. doi: 10.1038/s41598-022-17874-w

PubMed Abstract | Crossref Full Text | Google Scholar

Syarif Sukri, A., Saripuddin, M., Nasrul,, Karama, R., Talanipa, R., Kadir, A., et al. (2023). Utilization management to ensure clean water sources in coastal areas. J. Hum. Earth Future 4, 23–35. doi: 10.28991/hef-2023-04-01-03

Crossref Full Text | Google Scholar

Tran, P. Q., Bachand, S. C., McIntyre, P. B., Kraemer, B. M., Vadeboncoeur, Y., Kimirei, I. A., et al. (2021). Depth-discrete metagenomics reveals the roles of microbes in biogeochemical cycling in the tropical freshwater Lake Tanganyika. ISME J. 15, 1971–1986. doi: 10.1038/s41396-021-00898-x

PubMed Abstract | Crossref Full Text | Google Scholar

Tuan Tran, H., Lin, C., Bui, X. T., Ky Nguyen, M., Dan Thanh Cao, N., Mukhtar, H., et al. (2022). Phthalates in the environment: characteristics, fate and transport, and advanced wastewater treatment technologies. Bioresour. Technol. 344:126249. doi: 10.1016/j.biortech.2021.126249

PubMed Abstract | Crossref Full Text | Google Scholar

Vagts, J., Scheve, S., Kant, M., Wöhlbrand, L., and Rabus, R. (2018). Towards the response threshold for p-hydroxyacetophenone in the denitrifying bacterium “Aromatoleum aromaticum” EbN1. Appl. Environ. Microbiol. 84:e01018. doi: 10.1128/AEM.01018-18

PubMed Abstract | Crossref Full Text | Google Scholar

Voss, M., Bange, H. W., Dippner, J. W., Voss, M., Bange, H. W., Dippner, J. W., et al. (2013). The marine nitrogen cycle: recent discoveries, uncertainties and the potential relevance of climate change. Philos. Trans. R. Soc. B 368:20130121. doi: 10.1098/rstb.2013.0121

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/aem.00062-07

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Y., Guo, M., Li, X., Liu, G., Hua, Y., Zhao, J., et al. (2022). Shifts in microbial communities in shallow lakes depending on trophic states: feasibility as an evaluation index for eutrophication. Ecol. Indic. 136:108691. doi: 10.1016/j.ecolind.2022.108691

Crossref Full Text | Google Scholar

Wang, J., Shen, J., Wu, Y., Tu, C., Soininen, J., Stegen, J. C., et al. (2013). Phylogenetic beta diversity in bacterial assemblages across ecosystems: deterministic versus stochastic processes. ISME J. 7, 1310–1321. doi: 10.1038/ismej.2013.30

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Q., Yang, X., Hamilton, P. B., and Zhang, E. (2012). Linking spatial distributions of sediment diatom assemblages with hydrological depth profiles in a plateau deep-water lake system of subtropical China. Fottea 12, 59–73. doi: 10.5507/fot.2012.005

PubMed Abstract | Crossref Full Text | Google Scholar

Wasmund, K., Mussmann, M., and Loy, A. (2017). The life sulfuric: microbial ecology of sulfur cycling in marine sediments. Environ. Microbiol. Rep. 9, 323–344. doi: 10.1111/1758-2229.12538

PubMed Abstract | Crossref Full Text | Google Scholar

Watanabe, T., Kojima, H., Takano, Y., and Fukui, M. (2013). Diversity of sulfur-cycle prokaryotes in freshwater lake sediments investigated using aprA as the functional marker gene. Syst. Appl. Microbiol. 36, 436–443. doi: 10.1016/j.syapm.2013.04.009

Crossref Full Text | Google Scholar

Winter, C., Kerros, M.-E., and Weinbauer, M. G. (2009). Seasonal changes of bacterial and archaeal communities in the dark ocean: evidence from the Mediterranean Sea. Limnol. Oceanogr. 54, 160–170. doi: 10.4319/lo.2009.54.1.0160

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

Wu, K., Zhao, W., Wang, Q., Yang, X., Zhu, L., Shen, J., et al. (2019). The relative abundance of benthic bacterial phyla along a water-depth gradient in a plateau Lake: physical, chemical, and biotic drivers. Front. Microbiol. 10:1521. doi: 10.3389/fmicb.2019.01521

PubMed Abstract | Crossref Full Text | Google Scholar

Xing, P., Tao, Y., Luo, J., Wang, L., Li, B., Li, H., et al. (2019). Stratification of microbiomes during the holomictic period of Lake Fuxian, an alpine monomictic lake. Limnol. Oceanogr. 65, S134–S148. doi: 10.1002/lno.11346

Crossref Full Text | Google Scholar

Yamamoto, K., Matsutani, M., Shiwa, Y., Ishige, T., Sakamoto, H., Saitoh, H., et al. (2020). Comparative analysis of bacterial diversity and community structure in the rhizosphere and root endosphere of two halophytes, Salicornia europaea and Glaux maritima, collected from two brackish lakes in Japan. Microbes Environ. 35:ME20072. doi: 10.1264/jsme2.ME20072

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, Q., Jiang, Z., Zhou, X., Zhang, R., Xie, Z., Zhang, S., et al. (2020). Haliea alexandrii sp. nov., isolated from phycosphere microbiota of the toxin-producing dinoflagellate Alexandrium catenella. Int. J. Syst. Evol. Microbiol. 70, 1133–1138. doi: 10.1099/ijsem.0.003890

PubMed Abstract | Crossref Full Text | Google Scholar

Yin, X., Chen, L., Tang, D., Zhang, Y., Liu, G., Hua, Y., et al. (2019). Seasonal and vertical variations in the characteristics of the nitrogen-related functional genes in sediments from urban eutrophic lakes. Appl. Soil Ecol. 143, 80–88. doi: 10.1016/j.apsoil.2019.05.027

Crossref Full Text | Google Scholar

Zehr, J. P., and Capone, D. G. (2020). Changing perspectives in marine nitrogen fixation. Science 368:eaay9514. doi: 10.1126/science.aay9514

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, B., Liu, J., Sheng, Y., Shi, J., and Dong, H. (2021). Disentangling microbial syntrophic mechanisms for hexavalent chromium reduction in autotrophic biosystems. Environ. Sci. Technol. 55, 6340–6351. doi: 10.1021/acs.est.1c00383

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, P., Ren, M., Zhang, W., Xu, Y., and Wang, J. (2024). Quantifying the functional genes of C, N, P, and S cycling in a deep lake: depth patterns and drivers. Ecol. Indic. 166:112532. doi: 10.1016/j.ecolind.2024.112532

Crossref Full Text | Google Scholar

Zhao, W., Hu, A., Ni, Z., Wang, Q., Zhang, E., Yang, X., et al. (2019). Biodiversity patterns across taxonomic groups along a lake water-depth gradient: effects of abiotic and biotic drivers. Sci. Total Environ. 686, 1262–1271. doi: 10.1016/j.scitotenv.2019.05.381

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, B., Hu, Y., Yu, H., Chen, S., Xing, T., Guo, S., et al. (2023). A method for researching the eutrophication and N/P loads of plateau lakes: Lugu Lake as a case. Sci. Total Environ. 876:162747. doi: 10.1016/j.scitotenv.2023.162747

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, W., Pan, F., Meng, F., Soininen, J., Tupola, V., Ren, M., et al. (2021). Cross-taxon congruence of aquatic microbial communities across geological ages in Iceland: stochastic and deterministic processes. Sci. Total Environ. 774:145103. doi: 10.1016/j.scitotenv.2021.145103

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, J., He, Y., Zhu, Y., Huang, M., and Zhang, Y. (2018). Biogeochemical sulfur cycling coupling with dissimilatory nitrate reduction processes in freshwater sediments. Environ. Rev. 26, 121–132. doi: 10.1139/er-2017-0047

Crossref Full Text | Google Scholar

Keywords: microbial community, functional gene, metagenomic, water depth, deep lake

Citation: Zhang P, Ren M, Xu Y and Wang J (2025) Metagenomic insights into surface sediment microbial community and functional composition along a water-depth gradient in a subtropic deep lake. Front. Microbiol. 16:1614055. doi: 10.3389/fmicb.2025.1614055

Received: 18 April 2025; Accepted: 10 July 2025;
Published: 30 July 2025.

Edited by:

Hideyuki Tamaki, National Institute of Advanced Industrial Science and Technology (AIST), Japan

Reviewed by:

Da Huo, Chinese Academy of Sciences (CAS), China
Xianzhe Gong, Shandong University, China

Copyright © 2025 Zhang, Ren, Xu and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Minglei Ren, bWxyZW5AbmlnbGFzLmFjLmNu; Yan Xu, eHV4dWNhbG1tQHNldS5lZHUuY24=

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