- 1College of Geography and Environmental Science, Zhejiang Normal University, Jinhua, China
- 2The Management Center of Wuyanling National Natural Reserve in Zhejiang, Wenzhou, China
Water reservoir sediments represent a distinct habitat that harbors diverse microbial resources crucial for nitrogen cycling processes. The discovery of resuscitation promoting factor (Rpf) has been recognized as a crucial development in understanding the potential of microbial populations. However, our understanding of the relationship between microorganisms containing rpf-like genes and nitrogen-cycling functional populations remains limited. The present study explored the distribution patterns of rpf-like genes and nitrogen-cycling genes in various water reservoir sediments, along with their correlation with environmental factors. Additionally, the co-occurrence of rpf-like genes with genes associated with the nitrogen cycle and viable but non-culturable (VBNC) formation was investigated. The findings indicated the ubiquitous occurrence of Rpf-like domains and their related genes in the examined reservoir sediments. Notably, rpf-like genes were predominantly associated with Bradyrhizobium, Nitrospira, and Anaeromyxobacter, with pH emerging as the primary influencing factor for their distribution. Genera such as Nitrospira, Bradyrhizobium, Anaeromyxobacter, and Dechloromonas harbor the majority of nitrogen-cycling functional genes, particularly denitrification genes. The distribution of nitrogen-cycling microbial communities in the reservoir sediments was mainly influenced by pH and NH4+. Notably, correlation network analysis revealed close connections between microorganisms containing rpf-like genes and nitrogen-cycling functional populations, as well as VBNC bacteria. These findings offer new insights into the prevalence of rpf-like genes in the water reservoir sediments and their correlation with nitrogen-cycling microbial communities, enhancing our understanding of the significant potential of microbial nitrogen cycling.
1 Introduction
Reservoirs, formed by the construction of dams for purposes such as hydroelectricity production, flood control, and recreation activities, offer a distinctive ecological habitat (Zhang et al., 2015; Yang et al., 2021). The significance of reservoirs has grown in meeting the increasing energy and water demands. The importance of maintaining favorable chemical and biological water quality in reservoirs cannot be overstated, as it not only helps in reducing treatment costs but also ensures the supply of safe drinking water. Acting as both sinks and sources in nutrient cycles within water reservoirs, sediments serve as vital components of element migration and transformation (Zhang et al., 2015; Li et al., 2019). Microorganisms in water reservoir sediments, which exhibit high diversity and variability, are crucial for pollutant degradation and nutrient biogeochemical cycling, contributing significantly to the overall health and functioning of reservoir ecosystems (Gao et al., 2021; Wu et al., 2021).
Considerable attention has been directed towards examining the distribution and dynamics of functional microbial populations in reservoir sediments, particularly those involved in nitrogen cycling (Yang et al., 2021; Zhang et al., 2021; Lv et al., 2022; Shi et al., 2022). Microbe-driven nitrogen cycling encompasses various distinct processes, such as nitrogen fixation, nitrification, denitrification, anaerobic ammonium oxidation (anammox), dissimilatory/assimilatory nitrate reduction to ammonium (DNRA/ANRA), and ammonia assimilation/ammonification (Kuypers et al., 2018; Wang et al., 2022a; Nie et al., 2023; Zhang et al., 2023c). Previous research has detected microorganisms involved in these processes in reservoir ecosystems using molecular methods. However, most of these studies have focused on specific nitrogen cycling processes or functional populations by targeting limited genes through high-throughput amplicon sequencing and/or quantitative real-time PCR methods (Wang et al., 2021; Sun et al., 2023; Yuan et al., 2023). Comprehensive investigations of various nitrogen-cycling functional genes in water reservoir sediments using metagenomic approaches are still scarce. Thus, it is essential to conduct a comprehensive investigation into the functional microbial communities that participate in nitrogen cycles within reservoir sediments.
Additionally, water reservoir sediments serve as a unique habitat harboring diverse microbial resources. However, the identification of functional strains within these sediments remains limited. It should be noted that the majority of microorganisms on Earth are inaccessible and cannot be cultivated in pure cultures, with less than 1% being successfully cultured (Lin et al., 2023; Zhou et al., 2023). A significant challenge in obtaining functional strains arises from the presence of microorganisms in the viable but non-culturable (VBNC) state, which are incapable of forming colonies on culture plates (Su et al., 2018; Han et al., 2023). The VBNC state, similar to other non-growing states like persistence and dormancy, is a common occurrence for long-term survival under stressful conditions (Xie et al., 2021; Lin et al., 2023). The conditions that trigger the VBNC state, such as extreme temperatures, osmotic pressure, pH fluctuations, starvation, and high concentrations of pollutants, are commonly found in natural environments (Han et al., 2023; Yang et al., 2024). To date, more than 100 microbial species have been documented to enter the VBNC state (Su et al., 2023). VBNC microorganisms have been recognized in diverse environments, and the resuscitation of these microorganisms plays a vital role in environmental bioremediation and biogeochemical cycles (Lennon and Jones, 2011; Shi et al., 2024a).
The discovery of resuscitation promoting factor (Rpf) from Micrococcus luteus has been a significant breakthrough in the recovery of VBNC microorganisms, due to its potent anti-dormancy effects at picomolar concentrations (Su et al., 2023; Zhou et al., 2023). Studies have demonstrated its ability to resuscitate and enhance the activity of functional populations engaged in the degradation of diverse xenobiotics, including phenol (Su et al., 2018), polychlorinated biphenyls (PCBs) (Su et al., 2019c), and anthraquinone dyes (Cai et al., 2021). Notably, a strain Pseudomonas sp. SSPR1, resuscitated through Rpf supplementation, exhibited excellent performance in simultaneous nitrification and denitrification (Su et al., 2019b). The Rpf-responsive bacterial communities are distributed across phyla Actinobacteria, Proteobacteria, Bacteroidetes, and Firmicutes, encompassing numerous functional populations associated with biogeochemical cycles (Su et al., 2019a, 2023; Zhou et al., 2023). Furthermore, rpf-like genes have been identified in various genera, such as Rhodococcus, Mycobacterium, Streptomyces, Achromobacter, Nocardia, Corynebacterium, and Saccharopolyspora, exhibiting similar muralytic activity (Mukamolova et al., 2006; Kana et al., 2008; Hoang et al., 2021; Shi et al., 2024b). These genes are widely distributed across various environments and are closely linked to the transition between the VBNC and active states in microbial communities (Lennon and Jones, 2011; Shi et al., 2024b). Although previous studies have explored the use of Rpf homologues in vaccine development and bioremediation enhancements (Kana et al., 2008; Su et al., 2018; Zhou et al., 2023), their presence in natural environments and their potential correlation with nitrogen-cycling genes have yet to be investigated. Therefore, exploring the distribution and co-occurrence of rpf-like genes with the nitrogen cycle and VBNC formation in water reservoir sediments is essential.
In this study, three water reservoirs located in Wuyanling, a national-level nature reserve in China nearest to the East China Sea, were chosen for investigating sediment-associated functional microbial communities through metagenomic sequencing. The study aimed to: (1) demonstrate the structure and composition of sediment microbial communities and their correlation with environmental factors in the three water reservoirs; (2) analyze the distribution patterns of rpf-like genes and nitrogen-cycling genes and their association with environmental factors in diverse water reservoir sediments; (3) explore the co-occurrence of rpf-like genes with genes involved in the nitrogen cycle and VBNC formation. The present study is the first to report on the presence of rpf-like genes in water reservoir sediments and their association with functional bacterial populations participating in nitrogen cycles. These findings shed light on the crucial role of sediment microbes in nitrogen cycles within water reservoirs and highlight the potential of Rpf-like proteins in stimulating functional populations.
2 Materials and methods
2.1 Study site and sampling
During the optimal season for nitrogen cycling processes, sediment sampling was conducted in July 2022 at three water reservoirs, namely Sanchaxi (SCX), Zhangnenzi (ZNZ), and Wenyang (WY), located in the Wuyanling National Nature Reserve in the Zhejiang Province, China. These reservoirs are crucial for hydroelectricity production, with SCX being the largest among them. A total of five sampling sites (SCX1–SCX5) were designated for SCX, while three sampling sites were designed for both ZNZ (ZNZ1–ZNZ3) and WY (WY1–WY3). The sampling points within each reservoir were sequentially numbered from upstream to downstream. Detailed information regarding the sampling sites can be found in Supplementary Figure S1. At each sampling site, surface sediments were collected from the 0–30 cm depth using a Petersen stainless steel grab sampler. To ensure sample representativeness, a composite surface sediment sample was created by combining three sediment samples collected in close proximity to the sampling site (Zhang et al., 2015; Gao et al., 2021). Additionally, corresponding water samples were collected using a vertical water sampler. The sediment samples were sealed in sample bags, and the water samples were stored in polyethylene bottles. All samples were then placed in ice boxes for transportation to the laboratory. Subsequently, the water samples underwent immediate analysis of their physicochemical characteristics, while the sediment samples were stored separately in sampling bags (at 4°C) and polypropylene tubes (at −80°C) for further analysis of physicochemical characteristics and microbial communities, respectively.
2.2 Physicochemical characteristics analysis
The sediment pH was measured using a PHSJ-4F pH meter (Leici, China) with a sediment suspension ratio of 1:2.5 (sediment to deionized water). Prior to measurement, the deionized water was purged with nitrogen gas for 15 min. The moisture content (MC) of sediments was assessed by oven-drying fresh sediment at 105°C until a constant weight was achieved (Wang et al., 2021). To extract exchangeable dissolved inorganic nitrogen, including ammonium (NH4+), nitrite (NO2−), and nitrate (NO3−), from fresh sediments, a 2 M KCl solution pre-purged with high-purity nitrogen gas was utilized at a ratio of 1:5 (sediment to KCl solution, w/v). The resulting solution was filtered through a 0.45 μm membrane filter and then measured using a TU-1810 UV/VIS spectrophotometer, following the standard method (Li et al., 2019). The total organic matter (TOM) content was estimated by measuring the weight loss after combustion at 550°C using the LOI550 (loss on ignition at 550°C) method (Wang et al., 2021). Additionally, water quality parameters including pH, NH4+, NO2−, NO3−, total phosphorus, total organic carbon, and total nitrogen were determined using standard methods (Rice et al., 2012). All the physicochemical characteristics were measured in triplicate. The summaries of the physicochemical characteristics of both the sediment and water samples are presented in Supplementary Tables S1, S2, respectively.
2.3 Metagenomic sequencing
The total metagenomic DNA of the sediment samples was extracted using FastDNA Spin Kit for Soil (MP Biomedicals, United States). The concentration of the extracted DNA was determined using the TBS-380, while the purity was assessed using the NanoDrop2000. Equal masses of DNA from each of the three replicate samples were consolidated to reduce potential variability. The DNA extract was then fragmented to an average size of approximately 400 bp using the Covaris M220 to facilitate the construction of paired-end libraries. For library construction, the NEXTflex rapid DNA-seq kit (Bioo Scientific, United States) was utilized, with adapters that included the full set of sequencing primer hybridization sites ligated to the blunt ends of the DNA fragments. Paired-end sequencing was conducted on the Illumina Novaseq 6,000 (Illumina Inc., San Diego, CA, United States) at Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China). The raw sequence data have been deposited in the NCBI Sequence Read Archive under the project accession number PRJNA1061319.
2.4 Metagenome assembly and annotation
To ensure the reliability of the metagenomic data, adaptors and low-quality reads were trimmed. The remaining high-quality reads were then assembled using MEGAHIT, with contigs ≥300 bp chosen as the ultimate assembly. Gene annotation and prediction were performed on the contigs. Gene abundance quantification was conducted through the reads per kilobase per million mapped reads (RPKM) method. For each contig, open reading frames (ORFs) were predicted using MetaGene (Noguchi et al., 2006), and those with a length of ≥100 bp were extracted for translation into amino acid sequences using the NCBI translation table. Subsequently, a non-redundant gene catalog was built using CD-HIT, applying a 95% sequence identity and 90% coverage threshold. The high-quality reads were then aligned to this non-redundant gene catalog using SOAPaligner for determining gene abundance with 95% identity. For taxonomic annotations, representative sequences from the non-redundant gene catalog were aligned to the NR database using Diamond (Buchfink et al., 2015). The alpha diversity of the microbial community was assessed by calculating the number of reads at the genus level. The Conserved Domain Database (CDD) was utilized to annotate the structural domains of non-redundant genes (Jia et al., 2022; Shi et al., 2024b). The functions of non-redundant genes were annotated against the Kyoto Encyclopedia of Genes and Genomes (KEGG) database using Diamond (Buchfink et al., 2015).
2.5 Identification of rpf-like genes and genes associated with the nitrogen cycle and VBNC formation
The presence of rpf-like genes in the water reservoir sediments was investigated by analyzing the detected Rpf-like domains (Supplementary Table S3), including LT_TF-like, chitinase_GH19, lyz_endolysin_autolysin, chitosanase_GH46, pesticin_lyz-like, N-acetylmuramidase_GH108, lambda_lys-like, LT-like, Slt35-like, LT_IagB-like, Slt70-like, MLTF-like, LT_VirB1-like, LT_MltC_MltE, LT_Slt70-like, MltD-like, and Lyz-like.1 Moreover, an investigation was conducted on the distribution of nitrogen-cycling genes (Supplementary Table S4) encompassing nitrogen fixation (nifH, nifD, nifK), nitrification (amoA/pmoA, amoB/pmoB, amoC/pmoC, hao, nxrB), denitrification (narG, nirS, nirK, nosZ, norB), DNRA (nrfA, nirB, nirD, nrfH), ANRA (nasA, nirA, narB), anammox (hzsA_1/hzsA_2, hdh/hzo, hzsB, hzsC), and ammonification (ureC, gdhA) in the water reservoir sediments (Kuypers et al., 2018; Zhang et al., 2022; Wang et al., 2022a). Additionally, the genes involved in VBNC formation (Supplementary Table S5) were also examined in the water reservoir sediments to explore their potential co-occurrence with rpf-like genes and nitrogen-cycling functional genes.
2.6 Statistical analysis
The physicochemical properties data obtained from the experiment in triplicate were presented as mean values and standard deviations. Principal coordinates analysis (PCoA) with Bray–Curtis dissimilarity was utilized to assess the dissimilarity of sediment microbial communities from three water reservoirs (Zhou et al., 2023). The significance of these differences was evaluated using the analysis of similarities (ANOSIM) test. To compare the microbial community variations among the three water reservoir sediments, the Kruskal–Wallis H test was performed, followed by the Tukey–Kramer test at a 95% confidence level. Redundancy analysis (RDA) was employed to examine the association between environmental variables and functional genes, as well as their associated microbial communities. The environmental variables considered in the RDA included pH, NH4+, NO2−, NO3−, and TOM, which were selected to avoid collinearity and assessed using variance inflation factors. The significance of the RDA analysis was determined using Permutest (Shi et al., 2022). The relationships between functional genes and environmental factors, as well as dominant taxa and environmental factors, were evaluated using Spearman’s correlation coefficient (Wang et al., 2022a). Mantel test was employed to determine the correlation between different microbial communities and environmental factors. Additionally, to mitigate collinearity effects among variables, only functional genes or genera present in more than half (>50%) of all 11 sediment samples were included for network construction (Liu et al., 2023). Co-occurrence networks were generated using Spearman rank correlations, setting a threshold of |r| ≥ 0.7 and p < 0.01 to investigate the co-occurrence relationships among rpf-like genes, and genes related to nitrogen cycle and VBNC formation, as well as among the gene-associated genera. Spearman correlation coefficients were computed using the R package “psych,” and the resulting network diagrams were visualized using the Gephi 0.9.2 platform (Liu et al., 2023).
3 Results and discussion
3.1 Microbial diversity and community composition
The physiochemical properties (Supplementary Table S1) of the sediment samples exhibited notable spatial heterogeneity at a statistically significant level (p < 0.05). SCX sediments were marginally alkaline, while WY sediments were slightly acidic. ZNZ sediments had varying pH levels (6.80 to 8.51), with ZNZ1 upstream being alkaline, potentially due to proximity to a village and impact of domestic wastewater. WY sediments had the highest MC values (58.26–71.39%), followed by SCX (54.45–63.94%) and ZNZ (33.54–53.93%). The trend in ammonia concentrations was WY > SCX > ZNZ, with a decrease from upstream to downstream. Nitrite concentrations followed a pattern of ZNZ > SCX > WY. There were no significant differences in NO3− and TOM concentrations among the reservoir sediments. High-quality metagenomic sequencing data was obtained for sediment samples from the three reservoirs (Supplementary Table S6), with alpha diversity indices of microbial communities detailed in Supplementary Table S7. The Shannon diversity of sediment samples showed a significant positive correlation with NH4+ (r = 0.76, p < 0.05). Previous studies have indicated that an appropriate range of ammonium nitrogen concentration can supply necessary nutrients for microbial growth, facilitating microbial reproduction in the absence of other nitrogen sources (Kuypers et al., 2018).
PCoA revealed significant differences in the microbial community compositions within sediment samples from the three reservoirs, as evidenced by ANOSIM test (r = 0.38, p = 0.014) (Supplementary Figure S2). Microbial communities at phylum and genus levels, with notable variations across the reservoirs, are illustrated in Figure 1. Proteobacteria (45.01–50.69%) was the predominant phylum in all sediment samples except ZNZ1, where Bacteroidetes (40.63%) dominated (Figure 1A). The prevalence of Proteobacteria in reservoir ecosystems is extensively documented (Li et al., 2019; Wu et al., 2021). Furthermore, significant differences (p < 0.05) in the relative abundance of 38 phyla were observed across the reservoirs (Figure 1B). In addition, pairwise comparisons indicated notable variations in Euryarchaeota between SCX and ZNZ (p < 0.01). Candidatus Omnitrophica exhibited significant differences in relative abundance between WY and SCX (p < 0.001), as well as between WY and ZNZ (p < 0.001). Lentisphaerae also demonstrated significant differences between WY and ZNZ (p < 0.05) (Figure 1B). Notably, Euryarchaeota, Candidatus Omnitrophica, and Lentisphaerae are known to play roles in geomicrobiological nitrogen cycling (Speth et al., 2016).
Figure 1. Microbial community composition and abundance in sediments from three reservoirs. (A) Relative abundance of microbial community at the phylum level. (B) Comparative analysis of the relative abundance of the top 15 most abundant microbial communities at the phylum level. (C) Relative abundance of microbial community at the genus level. (D) Comparative analysis of the relative abundance of the top 15 most abundant microbial communities at the genus level. The “Unclassified” column indicates the combined relative sequence abundances of taxonomically unclassified sequences. The “Other phyla”/"Other genera” column indicates combined relative sequence abundances of all rare phyla/genera, defined as relative sequence abundances of <1% for all samples. Significant differences are denoted by asterisks: *p < 0.05, **p < 0.01, and ***p < 0.001.
At the genus level, Bradyrhizobium (1.70–9.44%) and Anaeromyxobacter (0.52–7.03%) were the predominant genera in the three reservoirs, except for ZNZ1 (Figure 1C). These genera are known for their involvement in nitrogen fixation and the DNRA process in the nitrogen cycle (Nelson et al., 2016; Sun et al., 2023). In ZNZ1, the most predominant genus was the methane-producing bacterium Alkaliflexus (6.12%), favored by the alkaline conditions of ZNZ1, providing an optimal environment for the growth of Alkaliflexus (Wang et al., 2022b). Furthermore, Dechloromonas (2.95%), known as an organohalide-respiring bacterium, was the second most abundant genus in ZNZ1, associated with the degradation of aromatic compounds (Yuan et al., 2023). Additionally, based on multiple comparisons, the relative abundance of 261 genera, including Dechloromonas, exhibited statistically significant differences among the three reservoirs (p < 0.05) (Figure 1D). Pairwise comparisons revealed significant differences in Methanothrix between SCX and ZNZ (p < 0.05), and in Methanoregula between ZNZ and SCX (p < 0.01), as well as between ZNZ and WY (p < 0.01). Both Methanothrix and Methanoregula, members of the class Methanomicrobia, are known to play a crucial role in methane production (Mand and Metcalf, 2019; Wu et al., 2021). Significant variations in the relative abundance of the genus Opitutus, involved in the DNRA process (Nelson et al., 2016), were observed between all pairs of reservoirs (Figure 1D). It is noteworthy that a substantial proportion of the genera in the three reservoirs were unclassified (Figure 1C), indicating limited accessibility to most of them. Restoring the activity and culturability of VBNC microorganisms in the reservoirs is essential to bridge the significant gaps between functional strains obtained and those existing in the environments (Lennon and Jones, 2011; Zhou et al., 2023).
3.2 Distribution of rpf-like genes and their correlation with environmental factors
The presence of Rpf-like domains and associated genes in sediment samples from the three reservoirs is depicted in Figure 2A and detailed in Supplementary Table S3. Specifically, SCX, ZNZ, and WY contained 577, 1,246, and 327 genes respectively, spanning across 14, 17, and 14 Rpf-like domains. The samples SCX1 and ZN1, located closest to the upstream (near the village) in SCX and ZN, exhibited the greatest occurrence of rpf-like genes, while sample WY3 in the WY reservoir displayed the greatest occurrence of rpf-like genes. Notably, six Rpf-like domains, including MltD-like, Slt70-like, LT-like, MLTF-like, Slt35-like, and LT_Slt70-like, were particularly prevalent across the three reservoirs, collectively representing 87.58% of the total rpf-like genes. Moreover, a systematic assessment was conducted to investigate the relationship between rpf-like genes and environmental variables. RDA analysis showed that 45.93% of the total variation in rpf-like genes could be attributed to five environmental factors, including pH, NH4+, NO2−, NO3− and TOM. Importantly, pH and NH4+ were identified as exerting significant influences (p < 0.05) on the distribution of rpf-like genes in the sediment samples (Figure 2B). Specifically, the relative abundance of lyz_endolysin_autolysin exhibited a significant positive correlation (r = 0.63, p < 0.05) with NH4+ concentration, whereas the relative abundance of LT-like showed a significant negative correlation (r = −0.76, p < 0.01) with NH4+ (Figure 2C). Additionally, pesticin_lyz-like and LT_VirB1-like were significantly associated with TOM, while LT_Slt70-like displayed a significant positive correlation (r = 0.63, p < 0.05) with NO2− concentration. These findings suggested the significant influence of NH4+ concentration on shaping the distribution of rpf-like genes in sediment samples from water reservoirs.
Figure 2. Distribution of Rpf-like domains and associated genes in sediment samples, and their correlation with environmental factors. (A) Percentage of abundance for Rpf-like domains and associated genes at each sampling site. (B) Redundancy analysis (RDA) illustrating the relationship between rpf-like genes and physiochemical parameters of sediments. The variables significantly influencing rpf-like genes, as selected by the forward selection procedure, are highlighted in red on the ordinations. (C) Spearman correlation analysis between rpf-like genes and environmental factors. The red and blue ellipses indicate positive and negative correlations, respectively. MC, moisture content; TOM, total organic matter. Significant correlations are denoted by asterisks: *p < 0.05 and **p < 0.01.
In order to examine the distribution patterns of microbial communities harboring rpf-like genes in the reservoir sediments, an investigation into the microbial composition linked to the presence of rpf-like genes was carried out. At the phylum level (Figure 3A), Proteobacteria (31.56–78.04%) was the primary phylum harboring rpf-like genes in the reservoir sediments. Bacteroidetes, the second most abundant phylum containing rpf-like genes, exhibited an increased abundance from the upstream to downstream sections in the SCX and WY reservoirs. Remarkably, the presence of Bacteroidetes with rpf-like genes was significantly higher in the ZNZ1 (31.36%) compared to other sediment samples (0.51–10.81%). A range of 0.68 to 3.91% of Actinobacteria were also found to contain rpf-like genes, with the highest proportion detected in SCX1. Extensive documentation exists regarding the presence of rpf-like genes in high G + C Gram-positive Actinobacteria (Mukamolova et al., 2006; Shi et al., 2024b). At the class level (Figure 3B), Alphaproteobacteria (6.14–46.37%), a subgroup of the Proteobacteria, was the predominant class hosting rpf-like genes in the reservoir sediments, followed by Betaproteobacteria (6.47–32.62%). Additionally, the abundance of Bacteroidia (25.68%) containing rpf-like genes was notably elevated in the ZNZ1 in comparison to the other sediment samples, consistent with the findings from the phylum analysis. These results revealed that rpf-like genes in the three reservoirs were predominantly found in Proteobacteria, specifically in the Alphaproteobacteria and Betaproteobacteria.
Figure 3. Taxonomic distribution of microbial communities harboring rpf-like genes, and their associated environmental factors. Distribution of phyla (A) and classes (B) with abundance greater than 1%. (C) Circos figures illustrating relationships between sediments and genera with abundance greater than 1%. (D) Spearman correlation analysis between environmental factors and dominant genera (with abundance greater than 1%) as well as Streptomyces and Mycobacterium. The red and blue ellipses indicate positive and negative correlations, respectively. Significant correlations are denoted by asterisks: *p < 0.05 and **p < 0.01. (E) Redundancy analysis (RDA) demonstrating the link between microbial communities containing rpf-like genes and physiochemical parameters of sediments. The variables significantly influencing microbial communities, selected by the forward selection procedure, are highlighted in red on the ordinations. (F) Mantel test based on Spearman rank correlation exhibiting the correlations of environmental factors with microbial communities and the correlations among different physiochemical properties of sediments. MC, moisture content; TOM, total organic matter. Significance is denoted by asterisk: *p < 0.05.
At the genus level (Figure 3C), most of the genera containing rpf-like genes in the reservoir sediments were categorized as unclassified. In the SCX reservoir, Bradyrhizobium (5.24–10.20%) was the predominant genus, except for SCX1 where Anaeromyxobacter was dominant. Significant differences were observed in the distribution of genera containing rpf-like genes from the upstream to downstream sections of the ZNZ reservoir. The predominant genera were Thermotoga, Nitrospira, and Anaeromyxobacter in ZNZ1, ZNZ2, and ZNZ3, respectively. In the WY reservoir, Nitrospira (4.20–5.17%) was the dominant genus, except for WY3 where Dechloromonas was prevalent. It is noteworthy that genera harboring rpf-like genes, such as Bradyrhizobium, Anaeromyxobacter, and Nitrospira, have been extensively reported for their involvement in nitrogen cycling (Wang et al., 2021; Sun et al., 2023). Although several studies have reported the presence of rpf-like genes in Streptomyces and Mycobacterium (Mukamolova et al., 2006; Kana et al., 2008), Streptomyces with rpf-like genes were only found in SCX1 (0.08%) and SCX2 (0.33%), and Mycobacterium with rpf-like genes were only detected in SCX5 (0.08%) and ZNZ3 (0.13%).
The association between genera harboring rpf-like genes and environmental factors was examined through heat map analysis (Figure 3D). The relative abundance of Aromatoleum exhibited a significant (p < 0.05) and negative correlation with pH, while Dechloromonas showed a significant (p < 0.05) and negative correlation with NO3− concentration. Similarly, the relative abundances of Hyphomicrobium and Mycobacterium displayed significant (p < 0.05) and negative correlations with NH4+ concentration. The relative abundance of Permianibacter was significantly (p < 0.05) and negatively correlated with NO2− concentration. In contrast, the relative abundance of Smithella was significantly and positively correlated with MC (p < 0.05) and TOM (p < 0.01). The relative abundance of Plasticicumulans exhibited a significant (p < 0.05) and positive correlation with NO3−, while the relative abundance of Serpentinimonas showed a significant (p < 0.05) and positive correlation with pH. The heat map analysis revealed that environmental factors had diverse effects on species harboring rpf-like genes, thus influencing the microbial community structure. Furthermore, by integrating local environmental variables with community structures, the correlation between environmental factors and community compositions containing rpf-like genes was identified. As illustrated in Figure 3E, the first two components explained 41.18% of the microbial diversity containing rpf-like genes. The RDA results supported the significant influence of pH (p < 0.05) on the distribution of microbial communities with rpf-like genes in the reservoir sediments. Additionally, pairwise Spearman’s correlation and partial Mantel tests were further employed to reveal the environmental drivers of communities harboring rpf-like genes in the three reservoirs (Figure 3F). The microbial communities containing rpf-like genes in SCX and WY reservoirs were significantly (r = 0.58, p < 0.05; r = 0.52, p < 0.05) correlated with pH, indicating that pH exerted a primary influence on the species compositions containing rpf-like genes.
3.3 Microbial community and functional genes involved in nitrogen cycling and their correlation with environmental factors
Based on the KEGG annotation, nitrogen-cycling functional genes associated with nitrogen fixation, nitrification, denitrification, DNRA, ANRA, anammox, and ammonification were annotated in all sediment samples (Supplementary Figure S3 and Supplementary Table S4). Denitrification was found to be the predominant function, constituting 34.67% of the total nitrogen transformation functions (Supplementary Figure S4). Ammonification and DNRA genes accounted for 18.16 and 16.73% of the total nitrogen transformation functions, respectively, while nitrogen fixation genes exhibited the lowest proportion at 7.23%. Compared to denitrification, the DNRA process prefers highly reduced environments (Zhang et al., 2023c). The dominance of denitrification over DNRA in the reservoir sediments may be attributed to the less reduced environments. Notably, the key genes, hzsA, hzsB, and hdh, known to be involved in anammox (Kuypers et al., 2018; Wang et al., 2022a), were absent in all samples. This absence could be explained by the low tolerance of anammox bacteria to oxygen, which inhibits their activity even under low oxygen levels. Molecular oxygen can directly poison redox enzymes, particularly Fe-S proteins with low redox potentials (such as ferredoxins and hydrogenases), essential for electron transfer. Furthermore, during the oxygen reduction process, intracellular reactive oxygen species are generated due to electron leakage, leading to oxidative stress that further hampers the activity of anammox bacteria (Zhang et al., 2023b).
The heatmap illustrating the correlations between nitrogen-cycling genes and environmental factors is presented in Figure 4. Nitrogen fixation genes were predominantly located downstream of the WY reservoir, particularly in WY3, with a minor distribution observed in the downstream sections of SCX4 and SCX5 in the SCX reservoir. The abundance of nifH, involved in nitrogen fixation, demonstrated a significant positive correlation with the TOM. The relatively higher abundance of nifH in the downstream section may be attributed to reduced flow, resulting in TOM deposition (Wakelin et al., 2008). Nitrification genes were predominantly present in the WY reservoir, with the amoB gene, responsible for ammonia oxidation, exhibiting a significant positive correlation with NH4+ concentration, a relationship previously observed by Hallam et al. (2006). Denitrification and DNRA genes were mainly distributed in the ZNZ2, with a high distribution also observed in the WY reservoir. The abundance of nirS and norB genes was significantly associated with pH, while the fluctuation in nirD corresponded to that of NO2− concentration. Furthermore, ANRA genes were primarily distributed downstream of the ZNZ reservoir, with significant negative correlations observed between the abundance of nasA and nirA genes and NH4+ concentration. The narB gene exhibited significant correlations with NO3− and NO2− concentrations. Additionally, the abundance of ammonification genes in the three reservoirs displayed a decreasing trend from the upstream to downstream sections.
Figure 4. Heatmap depicting the abundance of functional genes related to the nitrogen cycle in reservoir sediment samples (left plot), with high abundance represented in red and low abundance in blue. The abundance of functional genes was normalized using row scaling. The right plot illustrates the correlations (positive in red, negative in blue) between nitrogen-cycling genes and environmental factors. MC, moisture content; TOM, total organic matter. Significance levels are indicated by asterisks: *p < 0.05 and **p < 0.01.
Following taxonomic assignment, a total of 617 genera within 61 phyla containing nitrogen-cycling genes were identified in all sediment samples. At the phylum level (Figure 5A), the majority of nitrogen-cycling functional genes were associated with Proteobacteria (47.40–63.66%), followed by Acidobacteria (1.64–15.67%), Actinobacteria (0.84–11.49%), Nitrospirae (0.62–8.58%), Bacteroidetes (0.40–16.72%), and Chloroflexi (1.40–5.84%). Notably, the abundance of Bacteroidetes in ZNZ1 was significantly higher compared to the other sediments. At the genus level (Figure 5B), apart from unclassified genera, 23 genera with relative abundance greater than 1% were detected, with Nitrospira (0.07–7.07%), Bradyrhizobium (0.34–7.32%), Anaeromyxobacter (0.03–13.93%), and Dechloromonas (0.30–12.53%) being the most dominant. To elucidate the primary chemical characteristics of sediment supporting the nitrogen-cycling microbial community, RDA was conducted to assess the sediment chemical properties exerting the most significant influence. RDA1 and RDA2 could explain 44.70% of the total variation (Figure 5C), with SCX and WY samples exhibiting greater clustering, while ZNZ samples appeared more dispersed. The composition of sediment nitrogen-cycling microbial communities in the SCX, ZNZ, and WY exhibited significant correlations with NH4+ (p < 0.01) and pH (p < 0.05). Furthermore, partial Mantel tests were employed to further elucidate the environmental drivers of nitrogen-cycling microbial communities in the sediments (Figure 5D). The community structure of denitrifiers was significantly influenced by pH (r = 0.41, p < 0.05) and MC (r = 0.42, p < 0.05), while the nitrogen-fixing microbial community was primarily affected by TOM (r = 0.61, p < 0.01). Figure 5D also illustrated the correlations among different sediment parameters, with the exclusion of high collinearity (|r| > 0.8) (Ding et al., 2022). Specifically, TOM was found to be significantly correlated with MC (r = 0.75, p < 0.05). The results suggested that TOM and MC may jointly or independently impact the composition of nitrogen-fixing microbial communities (Zhang et al., 2023a).
Figure 5. Taxonomic distribution of nitrogen-cycling functional microorganisms in reservoir sediments, and their associated environmental factors. Distribution of phyla (A) and genera (B) with abundance greater than 1%. (C) Redundancy analysis (RDA) illustrating the relationship between microbial communities hosting nitrogen-cycling functional genes and sediment physiochemical parameters. Variables significantly influencing microbial communities, identified through a forward selection procedure, are highlighted in red. (D) Mantel test, based on Spearman rank correlation, demonstrating the correlations of environmental factors with microbial communities involved in the six major nitrogen cycling processes in sediments, as well as the correlations among different physiochemical properties of the sediments. Significance is denoted by asterisk: *p < 0.05.
Additionally, the relative contributions of different taxa to nitrogen-cycling genes across various nitrogen transformation processes were evaluated (Figure 6). The majority of nitrogen fixation genes were associated with Anaeromyxobacter, Geobacter, and Methanoregula. Notably, Methanoregula was found to harbor the nifH gene, which encodes iron-containing electron transfer proteins (Kuypers et al., 2018), and exhibited significant enrichment in downstream sites, regulating the downstream nitrogen-fixing process. Regarding the nitrification process, amoA, amoB, amoC, hao, and nxrB were predominantly assigned to the genus Nitrospira, known for its ability to oxidize ammonia in ammonia-limited environments (Kits et al., 2017). The prevalence of Nitrospira in the SCX reservoir suggested its potential involvement in the complete ammonia oxidation process. Furthermore, the main denitrifiers, namely, Bradyrhizobium and Nitrospira, exhibited higher abundance downstream of the reservoirs compared to upstream. Bradyrhizobium, also identified as the predominant genus in the ANRA process, was mainly associated with the abundance of nasA, encoding the catalytic subunit of assimilatory nitrate reductase (Hu et al., 2022), and demonstrated increased prevalence at downstream sites. Hidalgo-García et al. (2019) demonstrated the generation of nitrous oxide through the interplay of the ANRA and denitrification processes, indicating the promotion of the initial step of the ANRA process at downstream sites, subsequently providing more substrate for denitrification. Meanwhile, most of the ammonification genes were also assigned to the genus Bradyrhizobium, as well as Anaeromyxobacter. In the DNRA process, nrfA and nrfH were predominantly attributed to the genus Anaeromyxobacter, and unclassified genera in Deltaproteobacteria and Anaeromyxobacteraceae, while nirB was mainly assigned to unclassified genera in Betaproteobacteria. Previous studies have highlighted the role of DNRA in generating NH4+ to enhance primary productivity in oligotrophic marine and soil systems (Lam et al., 2009; Shan et al., 2016), a process that may also occur in water reservoir sediments.
Figure 6. Taxonomic contributions to nitrogen-cycling functional genes across various nitrogen transformation processes in diverse water reservoir sediment samples.
3.4 Co-occurrence of rpf-like genes with genes associated with nitrogen cycle and VBNC formation
The network analysis based on Spearman correlation revealed co-occurrence patterns among rpf-like genes, nitrogen-cycling functional genes, and VBNC-related genes with strong correlations (Figure 7A). The network consisting of 81 nodes connected by 638 edges, with 88.1% positive and 11.9% negative interactions. Notably, the rpf-like genes demonstrated 11 positive connections and 4 negative connections with nitrogen-cycling functional genes. Among these, the MLTF-like gene, which showed high abundance in all three reservoirs, displayed the highest number of connections with nitrogen-cycling functional genes and primarily positive correlations, indicating its crucial role in the positive correlation between rpf-like genes and nitrogen-cycling functional genes. Furthermore, the rpf-like genes demonstrated 33 positive links and 13 negative links with VBNC-related genes, with the LT_Slt70-like gene specifically linked to VBNC-related genes through 28 positive connections. Additionally, VBNC-related genes showed the highest self-co-occurrence rate, comprising around 74.5% of the occurrences.
Figure 7. Co-occurrence patterns among rpf-like genes, nitrogen-cycling functional genes, and VBNC-related genes (A), as well as among associated genera (B). Links represent statistically significant associations between nodes (|r| ≥ 0.7, p < 0.01). The size of each node is proportional to the number of connections. The red link indicates a positive correlation, and the blue link indicates a negative correlation. The thickness of each connection between two nodes (i.e., edge) is proportional to the strength of the correlation.
To unravel the intricate microbial co-occurrence network, correlation network analysis was conducted on the gene-associated genera. As depicted in Figure 7B, the network consisted of 364 nodes connected by 1874 edges, with 82.0% representing positive correlations and 18.0% reflecting negative correlations among these genera. A modularity index exceeding 0.4 is widely acknowledged as indicative of a modular structure (Jiao et al., 2016), and the observed modularity index of 0.58 suggested a modular construction of functional microbial populations in the water reservoir sediments. Within the network, the key genera Pseudolabrys, Bradyrhizobium, Nitrospira, Anaeromyxobacter, and Rubrivivax, which all contained nitrogen-cycling functional genes, rpf-like genes, and VBNC-related genes, collectively constituted 5.4, 4.9, 3.3, 2.2, and 1.7% of the network, respectively. Previous studies have documented the participation of Pseudolabrys in denitrification processes (Kämpfer et al., 2006; Li et al., 2022). Notably, Pseudolabrys appeared to be a significant genus within the network, possessing all three types of genes, with a higher degree (degree: 105) and betweenness centrality. Dechloromonas was noted for harboring both nitrogen-cycling functional genes and rpf-like genes. Furthermore, 11 genera were found to contain nitrogen-cycling functional genes and VBNC-related genes (Supplementary Table S8), suggesting the possibility of microorganisms involved in the nitrogen cycle entering the VBNC state. The results indicated that these nitrogen-cycling functional microbial populations included microorganisms with rpf-like genes and VBNC-related genes. Overall, the interaction between rpf-like genes and genes involved in nitrogen cycle and VBNC formation, along with their associated genera, revealed a strong correlation of Rpf homologs with functional microbial communities in water reservoir sediments.
4 Conclusion
The present study demonstrated the distinct microbial communities and functional populations present in the water reservoir sediments. Importantly, it revealed the widespread occurrence of Rpf-like domains and related genes across the three reservoir sediments. The rpf-like genes were primarily associated with Bradyrhizobium, Nitrospira, and Anaeromyxobacter. The distribution of rpf-like genes and their related genera in the sediments showed correlations with various environmental factors, with pH emerging as the primary influential factor. Furthermore, denitrification genes were identified as the primary nitrogen-cycling functional genes in the reservoir sediments, followed by ammonification and DNRA genes, with nitrogen fixation genes were comparatively less abundant. The majority of nitrogen-cycling functional genes were found in Bradyrhizobium and Dechloromonas within Proteobacteria, as well as the genera Nitrospira and Anaeromyxobacter. The presence of nitrogen-cycling microbial communities in the reservoir sediments was predominantly influenced by pH and NH4+. Notably, a strong correlation was observed between rpf-like genes and genes associated with the nitrogen cycle and VBNC formation, indicating the potential of Rpf homologs in resuscitating and stimulating functional microbial populations in the water reservoir sediments. This study provided novel insights into the prevalence of rpf-like genes in water reservoir sediments and advanced our understanding of the relationship between microorganisms containing rpf-like genes and nitrogen-cycling functional populations, offering new perspectives for the huge potential of microbial nitrogen cycling in water reservoirs.
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 in the article/Supplementary material.
Author contributions
AH: Data curation, Investigation, Writing – original draft. HF: Data curation, Investigation, Writing – original draft. LL: Investigation, Resources, Writing – review & editing. XS: Conceptualization, Project administration, Supervision, Writing – original draft, Writing – review & editing. SZ: Project administration, Supervision, Writing – original draft, Writing – review & editing. JL: Resources, Software, Writing – review & editing. FS: Resources, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The authors appreciate the funding by the National Natural Science Foundation of China (No. 42277025), the Natural Science Foundation of Zhejiang Province, China (No. LY21D010006), and the Jinhua Science and Technology Bureau, China (No. 2022-4-033).
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.
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.2024.1433046/full#supplementary-material
Footnotes
References
Buchfink, B., Xie, C., and Huson, D. H. (2015). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60. doi: 10.1038/nmeth.3176
Cai, J. F., Pan, A. D., Li, Y. L., Xiao, Y. Y., Zhou, Y., Chen, C. J., et al. (2021). A novel strategy for enhancing anaerobic biodegradation of an anthraquinone dye reactive blue 19 with resuscitation-promoting factors. Chemosphere 263:127922. doi: 10.1016/j.chemosphere.2020.127922
Ding, Y. T., Pan, B. Z., Zhao, X. H., Zhao, G. N., Han, X., and Li, M. (2022). Will a heavy sediment load affect responses of phytoplankton functional groups to aquatic environmental changes in different water body types? Sci. Total Environ. 837:155863. doi: 10.1016/j.scitotenv.2022.155863
Gao, Y., Zhang, W. L., and Li, Y. (2021). Microbial community coalescence: does it matter in the three gorges reservoir? Water Res. 205:117638. doi: 10.1016/j.watres.2021.117638
Hallam, S. J., Mincer, T. J., Schleper, C., Preston, C. M., Roberts, K., Richardson, P. M., et al. (2006). Pathways of carbon assimilation and ammonia oxidation suggested by environmental genomic analyses of marine Crenarchaeota. PLoS Biol. 4:e95. doi: 10.1371/journal.pbio.0040095
Han, Z., Lin, Q. H., Zhang, S. S., Zhou, X. R., Li, S., Sun, F. Q., et al. (2023). High PCBs mineralization capability of a resuscitated strain Bacillus sp. LS1 and its survival in PCB-contaminated soil. Sci. Total Environ. 856:159224. doi: 10.1016/j.scitotenv.2022.159224
Hidalgo-García, A., Torres, M. J., Salas, A., Bedmar, E. J., Girard, L., and Delgado, M. J. (2019). Rhizobium etli produces nitrous oxide by coupling the assimilatory and denitrification pathways. Front. Microbiol. 10:980. doi: 10.3389/fmicb.2019.00980
Hoang, H. H., Wang, P. C., and Chen, S. C. (2021). Recombinant resuscitation-promoting factor protein of Nocardia seriolae, a promissing vaccine candidate for largemouth bass (Micropterus salmoides). Fish Shellfish Immunol. 111, 127–139. doi: 10.1016/j.fsi.2021.01.015
Hu, R. W., Liu, S. F., Huang, W. M., Nan, Q., Strong, P. J., Saleem, M., et al. (2022). Evidence for assimilatory nitrate reduction as a previously overlooked pathway of reactive nitrogen transformation in estuarine suspended particulate matter. Environ. Sci. Technol. 56, 14852–14866. doi: 10.1021/acs.est.2c04390
Jia, Y. Y., Li, X. Y., Xu, F. J., Liu, Z. F., Fu, Y. L., Xu, X., et al. (2022). Single-cell-level microfluidics assisted with resuscitation-promoting factor technology (SMART) to isolate novel biphenyl-degrading bacteria from typical soils in eastern China. Environ. Pollut. 311:119864. doi: 10.1016/j.envpol.2022.119864
Jiao, S., Liu, Z. S., Lin, Y. B., Yang, J., Chen, W. M., and Wei, G. H. (2016). Bacterial communities in oil contaminated soils: biogeography and co-occurrence patterns. Soil Biol. Biochem. 98, 64–73. doi: 10.1016/j.soilbio.2016.04.005
Kämpfer, P., Young, C. C., Arun, A. B., Shen, F. T., Jäckel, U., Rosselló-Mora, R., et al. (2006). Pseudolabrys taiwanensis gen. Nov., sp. nov., an alphaproteobacterium isolated from soil. Int. J. Syst. Evol. Microbiol. 56, 2469–2472. doi: 10.1099/ijs.0.64124-0
Kana, B. D., Gordhan, B. G., Downing, K. J., Sung, N., Vostroktunova, G., Machowski, E. E., et al. (2008). The resuscitation-promoting factors of Mycobacterium tuberculosis are required for virulence and resuscitation from dormancy but are collectively dispensable for growth in vitro. Mol. Microbiol. 67, 672–684. doi: 10.1111/j.1365-2958.2007.06078.x
Kits, K. D., Sedlacek, C. J., Lebedeva, E. V., Han, P., Bulaev, A., Pjevac, P., et al. (2017). Kinetic analysis of a complete nitrifier reveals an oligotrophic lifestyle. Nature 549, 269–272. doi: 10.1038/nature23679
Kuypers, M. M., Marchant, H. K., and Kartal, B. (2018). The microbial nitrogen-cycling network. Nat. Rev. Microbiol. 16, 263–276. doi: 10.1038/nrmicro.2018.9
Lam, P., Lavik, G., Jensen, M. M., van de Vossenberg, J., Schmid, M., Woebken, D., et al. (2009). Revising the nitrogen cycle in the Peruvian oxygen minimum zone. PNAS 106, 4752–4757. doi: 10.1073/pnas.0812444106
Lennon, J. T., and Jones, S. E. (2011). Microbial seed banks: the ecological and evolutionary implications of dormancy. Nat. Rev. Microbiol. 9, 119–130. doi: 10.1038/nrmicro2504
Li, Y., Gao, Y., Zhang, W. L., Wang, C., Wang, P. F., Niu, L. H., et al. (2019). Homogeneous selection dominates the microbial community assembly in the sediment of the three gorges reservoir. Sci. Total Environ. 690, 50–60. doi: 10.1016/j.scitotenv.2019.07.014
Li, S. J., Liao, Y. H., Pang, Y. M., Dong, X. L., Strous, M., and Ji, G. D. (2022). Denitrification and dissimilatory nitrate reduction to ammonia in long-term lake sediment microcosms with iron(II). Sci. Total Environ. 807:150835. doi: 10.1016/j.scitotenv.2021.150835
Lin, Q. H., Ding, J. Y., Yang, Y. Y., Sun, F. Q., Shen, C. F., Lin, H. J., et al. (2023). Simultaneous adsorption and biodegradation of polychlorinated biphenyls using resuscitated strain Streptococcus sp. SPC0 immobilized in polyvinyl alcohol-sodium alginate. Sci. Total Environ. 868:161620. doi: 10.1016/j.scitotenv.2023.161620
Liu, J. J., Guo, Y. P., Gu, H. D., Liu, Z. X., Hu, X. J., Yu, Z. H., et al. (2023). Conversion of steppe to cropland increases spatial heterogeneity of soil functional genes. ISME J. 17, 1872–1883. doi: 10.1038/s41396-023-01496-9
Lv, K., Guo, X. J., Wang, C. F., Su, Q. Q., Liu, D. F., Xiao, S. B., et al. (2022). Sediment nitrogen contents controlled by microbial community in a eutrophic tributary in three gorges reservoir, China. Environ. Pollut. 314:120312. doi: 10.1016/j.envpol.2022.120312
Mand, T. D., and Metcalf, W. W. (2019). Energy conservation and hydrogenase function in methanogenic archaea, in particular the genus Methanosarcina. Microbiol. Mol. Biol. Rev. 83:e00020. doi: 10.1128/MMBR.00020-19
Mukamolova, G. V., Murzin, A. G., Salina, E. G., Demina, G. R., Kell, D. B., Kaprelyants, A. S., et al. (2006). Muralytic activity of Micrococcus luteus Rpf and its relationship to physiological activity in promoting bacterial growth and resuscitation. Mol. Microbiol. 59, 84–98. doi: 10.1111/j.1365-2958.2005.04930.x
Nelson, M. B., Martiny, A. C., and Martiny, J. B. H. (2016). Global biogeography of microbial nitrogen-cycling traits in soil. PNAS 113, 8033–8040. doi: 10.1073/pnas.1601070113
Nie, S. Q., Mo, S. M., Gao, T. W., Yan, B., Shen, P. H., Kashif, M., et al. (2023). Coupling effects of nitrate reduction and sulfur oxidation in a subtropical marine mangrove ecosystem with Spartina alterniflora invasion. Sci. Total Environ. 862:160930. doi: 10.1016/j.scitotenv.2022.160930
Noguchi, H., Park, J., and Takagi, T. (2006). MetaGene: prokaryotic gene finding from environmental genome shotgun sequences. Nucleic Acids Res. 34, 5623–5630. doi: 10.1093/nar/gkl723
Rice, E. W., and Bridgewater, L.Association, A.P.H (2012). Standard methods for the examination of water and wastewater, vol. 10. Washington, DC: American Public Health Association.
Shan, J., Zhao, X., Sheng, R., Xia, Y. Q., Ti, C. P., Quan, X. F., et al. (2016). Dissimilatory nitrate reduction processes in typical Chinese paddy soils: rates, relative contributions, and influencing factors. Environ. Sci. Technol. 50, 9972–9980. doi: 10.1021/acs.est.6b01765
Shi, J., Yang, Y. Y., Zhang, S. S., Lin, Q. H., Sun, F. Q., Lin, H. J., et al. (2024a). New insights into survival strategies and PCB bioremediation potential of resuscitated strain Achromobacter sp. HR2 under combined stress conditions. J. Hazard. Mater. 465:133242. doi: 10.1016/j.jhazmat.2023.133242
Shi, J. X., Zhang, B. G., Liu, J., Fang, Y., and Wang, A. J. (2022). Spatiotemporal dynamics in microbial communities mediating biogeochemical cycling of nutrients across the Xiaowan reservoir in Lancang River. Sci. Total Environ. 813:151862. doi: 10.1016/j.scitotenv.2021.151862
Shi, J., Zhou, X. R., Zhang, S. S., Sun, F. Q., Shen, C. F., and Su, X. M. (2024b). Unveiling the distribution characteristics of rpf-like genes and indigenous resuscitation promoting factor production in PCB-contaminated soils. J. Environ. Manag. 357:120803. doi: 10.1016/j.jenvman.2024.120803
Speth, D. R., IN’T Zandt, M. H., Guerrero-Cruz, S., Dutilh, B. E., and Jetten, M. S. M. (2016). Genome-based microbial ecology of anammox granules in a full-scale wastewater treatment system. Nat. Commun. 7:11172. doi: 10.1038/ncomms11172
Su, X. M., Wang, Y. Y., Xue, B. B., Hashmi, M. Z., Lin, H. J., Chen, J. R., et al. (2019a). Impact of resuscitation promoting factor (Rpf) in membrane bioreactor treating high-saline phenolic wastewater: performance robustness and Rpf-responsive bacterial populations. Chem. Eng. J. 357, 715–723. doi: 10.1016/j.cej.2018.09.197
Su, X. M., Wang, Y. Y., Xue, B. B., Zhang, Y. G., Mei, R. W., Zhang, Y., et al. (2018). Resuscitation of functional bacterial community for enhancing biodegradation of phenol under high salinity conditions based on Rpf. Bioresour. Technol. 261, 394–402. doi: 10.1016/j.biortech.2018.04.048
Su, X. M., Xie, M. Q., Han, Z., Xiao, Y. Y., Wang, R., Shen, C. F., et al. (2023). Resuscitation-promoting factor accelerates enrichment of highly active tetrachloroethene/polychlorinated biphenyl-dechlorinating cultures. Appl. Environ. Microbiol. 89, e01951–e01922. doi: 10.1128/aem.01951-22
Su, X. M., Xue, B. B., Wang, Y. Y., Hashmi, M. Z., Lin, H. J., Chen, J. R., et al. (2019b). Bacterial community shifts evaluation in the sediments of Puyang River and its nitrogen removal capabilities exploration by resuscitation promoting factor. Ecotoxicol. Environ. Saf. 179, 188–197. doi: 10.1016/j.ecoenv.2019.04.067
Su, X. M., Zhou, M., Hu, P., Xiao, Y. Y., Wang, Z., Mei, R. W., et al. (2019c). Whole-genome sequencing of an acidophilic Rhodotorula sp. ZM1 and its phenol-degrading capability under acidic conditions. Chemosphere 232, 76–86. doi: 10.1016/j.chemosphere.2019.05.195
Sun, X., Tao, R. D., Xu, D. Q., Qu, M. J., Zheng, M. M., Zhang, M., et al. (2023). Role of polyamide microplastic in altering microbial consortium and carbon and nitrogen cycles in a simulated agricultural soil microcosm. Chemosphere 312:137155. doi: 10.1016/j.chemosphere.2022.137155
Wakelin, S. A., Colloff, M. J., and Kookana, R. S. (2008). Effect of wastewater treatment plant effluent on microbial function and community structure in the sediment of a freshwater stream with variable seasonal flow. Appl. Environ. Microbiol. 74, 2659–2668. doi: 10.1128/AEM.02348-07
Wang, P. D., Li, J. L., Luo, X. Q., Ahmad, M., Duan, L., Yin, L. Z., et al. (2022a). Biogeographical distributions of nitrogen-cycling functional genes in a subtropical estuary. Funct. Ecol. 36, 187–201. doi: 10.1111/1365-2435.13949
Wang, S. Y., Wang, X. M., Jiang, Y. Y., Han, C., Jetten, M. S. M., Schwark, L., et al. (2021). Abundance and functional importance of complete ammonia oxidizers and other nitrifiers in a riparian ecosystem. Environ. Sci. Technol. 55, 4573–4584. doi: 10.1021/acs.est.0c00915
Wang, X. T., Zhang, Y. F., Wang, B., Wang, S., Xing, X., Xu, X. J., et al. (2022b). Enhancement of methane production from waste activated sludge using hybrid microbial electrolysis cells-anaerobic digestion (MEC-AD) process—a review. Bioresour. Technol. 346:126641. doi: 10.1016/j.biortech.2021.126641
Wu, D. B., Zhao, Y., Cheng, L., Zhou, Z., Wu, Q. S., Wang, Q., et al. (2021). Activity and structure of methanogenic microbial communities in sediments of cascade hydropower reservoirs, Southwest China. Sci. Total Environ. 786:147515. doi: 10.1016/j.scitotenv.2021.147515
Xie, M. Q., Xu, L. N., Zhang, R., Zhou, Y., Xiao, Y. Y., Su, X. M., et al. (2021). Viable but nonculturable state of yeast Candida sp. strain LN1 induced by high phenol concentrations. Appl. Environ. Microbiol. 87, e01110–e01121. doi: 10.1128/AEM.01110-21
Yang, Y. Y., Zhang, Q., Lin, Q. H., Sun, F. Q., Shen, C. F., Lin, H. J., et al. (2024). Unveiling the PCB biodegradation potential and stress survival strategies of resuscitated strain Pseudomonas sp. HR1. Environ. Pollut. 344:123320. doi: 10.1016/j.envpol.2024.123320
Yang, N., Zhang, C., Wang, L. Q., Li, Y., Zhang, W. L., Niu, L. H., et al. (2021). Nitrogen cycling processes and the role of multi-trophic microbiota in dam-induced river-reservoir systems. Water Res. 206:117730. doi: 10.1016/j.watres.2021.117730
Yuan, S. P., Han, X. K., Yin, X. Y., Su, P. X., Zhang, Y. Y., Liu, Y. F., et al. (2023). Nitrogen transformation promotes the anaerobic degradation of PAHs in water level fluctuation zone of the three gorges reservoir in Yangtze River, China: evidences derived from in-situ experiment. Sci. Total Environ. 864:161034. doi: 10.1016/j.scitotenv.2022.161034
Zhang, L., Dong, T. J., Yang, J. C., Hao, S. W., Sun, Z. C., and Peng, Y. Z. (2023b). Anammox coupled with photocatalyst for enhanced nitrogen removal and the activated aerobic respiration of anammox bacteria based on cbb3-type cytochrome c oxidase. Environ. Sci. Technol. 57, 17910–17919. doi: 10.1021/acs.est.3c02435
Zhang, H. H., Huang, T. L., Chen, S. N., Yang, X., Lv, K., and Sekar, R. (2015). Abundance and diversity of bacteria in oxygen minimum drinking water reservoir sediments studied by quantitative PCR and pyrosequencing. Microb. Ecol. 69, 618–629. doi: 10.1007/s00248-014-0539-6
Zhang, Z. F., Nie, S. Q., Sang, Y. M., Mo, S. M., Li, J. H., Kashif, M., et al. (2022). Effects of Spartina alterniflora invasion on nitrogen fixation and phosphorus solubilization in a subtropical marine mangrove ecosystem. Microbiol. Spectrum 10:e00682. doi: 10.1128/spectrum.00682-21
Zhang, H. J., Sun, L. W., Li, Y., Zhang, W. L., Niu, L. H., and Wang, L. F. (2021). The bacterial community structure and N-cycling gene abundance in response to dam construction in a riparian zone. Environ. Res. 194:110717. doi: 10.1016/j.envres.2021.110717
Zhang, H. H., Yang, Y. S., Liu, X., Huang, T. L., Ma, B., Li, N., et al. (2023a). Novel insights in seasonal dynamics and co-existence patterns of phytoplankton and micro-eukaryotes in drinking water reservoir, Northwest China: DNA data and ecological model. Sci. Total Environ. 857:159160. doi: 10.1016/j.scitotenv.2022.159160
Zhang, X. L., Yao, C., Zhang, B. S., Tan, W., Gong, J., Wang, G. Y., et al. (2023c). Dynamics of benthic nitrate reduction pathways and associated microbial communities responding to the development of seasonal deoxygenation in a castal mariculture zone. Environ. Sci. Technol. 57, 15014–15025. doi: 10.1021/acs.est.3c03994
Zhou, X. R., Zhang, S. S., Wang, R., An, Z. J., Sun, F. Q., Shen, C. F., et al. (2023). A novel strategy for enhancing bioremediation of polychlorinated biphenyl-contaminated soil with resuscitation promoting factor and resuscitated strain. J. Hazard. Mater. 447:130781. doi: 10.1016/j.jhazmat.2023.130781
Keywords: reservoir sediments, rpf-like genes, nitrogen-cycling genes, microbial community, co-occurrence network
Citation: Hou A, Fu H, Liu L, Su X, Zhang S, Lai J and Sun F (2024) Exploring the distribution and co-occurrence of rpf-like genes and nitrogen-cycling genes in water reservoir sediments. Front. Microbiol. 15:1433046. doi: 10.3389/fmicb.2024.1433046
Edited by:
Ren-Cun Jin, Hangzhou Normal University, ChinaReviewed by:
Qi Han, Beijing Forestry University, ChinaXianbiao Lin, Ocean University of China, China
Cheng Wang, Sun Yat-sen University, China
Caiyun Yang, Southwest University, China
Copyright © 2024 Hou, Fu, Liu, Su, Zhang, Lai and Sun. 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: Xiaomei Su, cHVycGxlQHpqbnUuY24=; Shusheng Zhang, enloeDI2NThAMTYzLmNvbQ==
†These authors have contributed equally to this work