Yiqiang Chen1,2
Xu Jiang3
Juanjuan Zhao4Maosheng Yang3Yong Chen1*Hong Ling1Yang Liu1Feng Deng4
Zhu Wang2*- 1Key Laboratory of Environmental Remediation and Ecological Health, Ministry of Industry and Information Technology, Jiangsu Environmental Engineering Technology Co., Ltd, Jiangsu Environmental Protection Group Co., Ltd., Nanjing, China
- 2Institute of Environmental Research at Greater Bay, Key Laboratory for Water Quality and Conservation of the Pearl River Delta, Ministry of Education, Guangzhou University, Guangzhou, China
- 3Key Laboratory of Integrated Regulation and Resources Development on Shallow Lakes, Ministry of Education, College of Environment, Hohai University, Nanjing, China
- 4College of Environmental Science and Engineering, Nanjing Tech University, Nanjing, China
This study investigated the responses of the bacterial community structure and metabolic pathways in a sulfur-based autotrophic denitrification filter (SADF) system to fast elevated sulfate salinity, from 0.04 to 1.2% in 30 days. Results showed that the SADF system exhibited robust sulfate salinity stress tolerance at low nitrate concentrations. In the context of sulfate scenarios, the genus Thiobacillus significantly proliferated and was identified as the dominant sulfur-oxidizing player in the SADF system, achieving a relative abundance of 63.79% under 1.2% sulfate salinity. Cooperative and competitive interactions were found in the SADF-related microorganisms, promoting stable denitrification performance under high salinity. Surprisingly, with a low hydraulic retention time (HRT) of 60 min, metagenomic sequencing revealed a upregulated abundance of functional genes encoding for enzymes associated with nitrogen and sulfur metabolism, while positive correlations were observed between these two pathways in response to sulfate salinity. Furthermore, global wastewater treatment plants were thoroughly explored for the distribution of the SADF-related microorganisms identified in this study. Interestingly, one-way ANOVA analysis showed that the SADF-related microorganisms were widely distributed globally, demonstrating their universality in potential engineering applications worldwide.
1 Introduction
Nitrate, a primary contaminant in wastewater effluents, poses significant ecological and health risks when discharged in excessive quantities (Abascal et al., 2022; Xu et al., 2025). Overloading of nitrate into aquatic systems can trigger severe eutrophication, characterized by algal blooms and hypoxic conditions, ultimately disrupting aquatic ecosystems (Wang L. et al., 2023; Wang T. et al., 2023; Yue et al., 2023). Beyond environmental concerns, elevated nitrate levels in drinking water sources directly threaten human health, contributing to potential health issues such as methemoglobinemia and increased risk of certain cancers (Picetti et al., 2022; Liu et al., 2024). The effective removal of nitrate from wastewater is therefore critical to mitigate environmental degradation and safeguard public health.
Biological denitrification, a microbial mediated process that functions under oxygen-limited conditions (DO <0.5 mg/L), could achieve complete reduction of nitrate to nitrogen gas (N₂) through two distinct metabolic pathways: (1) heterotrophic denitrification (Shao et al., 2025) utilizing organic carbon sources (e.g., methanol, acetate) and (2) autotrophic denitrification (Ye et al., 2024) driven by inorganic electron donors (e.g., H₂, SO). While this process provides an effective nitrogen removal strategy for wastewater treatment, conventional heterotrophic denitrification poses significant limitations, including (1) elevated operational costs from continuous organic carbon supplementation (Xi et al., 2022), (2) excessive sludge production (Wang et al., 2021), and (3) risks of secondary pollution due to residual chemical oxygen demand (COD) in effluents (Sun et al., 2017).
In alignment with global carbon neutrality goals, autotrophic denitrification has emerged as a promising alternative due to its low operational costs, reduced sludge yield, and minimal organic carbon residues. Among autotrophic denitrification technologies, SADF system has attracted increasing attention due to its pronounced advantages over conventional heterotrophic denitrification, including a 55% reduction in sludge generation, over 87% reduction in treatment costs, and significant decreases in greenhouse gas emissions (e.g., carbon dioxide and nitrous oxide) by 80–100% (Zheng et al., 2024; Wang L. et al., 2023; Wang T. et al., 2023; Pang and Wang, 2021).
Recent studies have demonstrated the feasibility of the sulfur autotrophic denitrification (SADN) process for treating nitrate-laden wastewater, achieving satisfactory nitrogen removal performance with total nitrogen removal efficiencies consistently exceeding 90%. For example, Zhao et al. (2024) invented a special filler based on combination with sulfur and siderite, achieving robust removal rates of 91.2% at an extreme cold temperature of 7°C. Li R. et al. (2021) and Li W. et al. (2021) successfully removed both arsenite and nitrate in the SADF system, with the maximum removal efficiency of nitrate at 98.42%. Bai et al. (2023) fused FeS and S with 2:1 as a filler, finding that the optimal hydraulic retention time (HRT) to be 3 h with a removal rate of 90–100%. These findings underscore the effectiveness of SADN as an eco-friendly technology for achieving sustainable wastewater treatment, particularly in scenarios with extreme low C/N ratio and expectations of low carbon emission.
Industrial wastewater from petroleum refining, pharmaceutical manufacturing, and textile dyeing processes typically contains high salinity, primarily composed of sodium chloride (NaCl) and sodium sulfate (Na₂SO₄). The performance of the SADN process varies significantly under different NaCl salinities when treating saline wastewater. Notably, the performance of the SADN process could be improved within a NaCl salinity range from 0.1 to 1.5%. For instance, previous study (Shen et al., 2023) reported a significant improvement in nitrate removal efficiency from 81 to 88%, while the NaCl salinity increased from 0.16 to 1.5%. These findings (Wang T. et al., 2024; Wang X. et al., 2024; Li Y. Y. et al., 2023; Li Z. R. et al., 2023; Chen B. et al., 2022; Chen H. et al., 2022) highlights the adaptability of the SADN process to moderate NaCl salinity levels, making it a promising approach for treating saline industrial wastewater. However, while chloride salts (e.g., NaCl) primarily exert osmotic stress on microbial cells, sulfate salts (e.g., Na₂SO₄) not only impose similar osmotic pressure but may also cause more complex and severe inhibitory effects due to the toxicity of their metabolic by-products (e.g., H₂S, SO₄2−) (Chatla et al., 2023). These differences suggest that sulfate-induced salinity may interfere with microbial communities through stronger electron acceptor competition, metabolic pathway disruption, and dramatic shifts in microbial structure.
Previous studies have primarily focused on the effects of NaCl salinity on SADN processes in wastewater treatment systems. Sulfate generated as a byproduct of SADN processes can significantly disrupt the chemical equilibrium under excessive conditions. Therefore, investigating the impact of sulfate-mediated salinity on the SADN process is of great importance, especially for optimizing nitrogen removal performance in high-sulfate industrial wastewater. However, there is limited understanding of the responses of the bacterial community structure and metabolic pathways of SADN process when applied to treat industrial wastewater characterized by high sulfate salinity.
In this study, a SADF system packed with elemental sulfur (S0) as the sole electron donor was specifically designed to thoroughly assess the effects of various Na₂SO₄ salinities (0.04–1.2%) on the performance of the SADN process and its underlying metabolic mechanisms. To our knowledge, no prior studies have specifically examined the effects of Na₂SO₄ salinity, unlike NaCl salinity, on the SADF system, especially with a fast elevated salinity in 30 days and a low HRT of 60 min. The key research aspects include: (1) assessing the robust tolerance to high salt concentrations and resistance mechanisms within the SADF system, (2) investigating the shift in microbial community structure with increasing Na₂SO₄ salinity, (3) analyzing the competitive and cooperative interactions among different microorganisms in the SADF system, (4) examining the upregulated fluctuations in the abundance of functional genes involved in nitrogen and sulfur metabolism under Na₂SO₄ salinity stress, (5) evaluating the global distribution of SADN-related microbial communities across wastewater treatment plants. This comprehensive investigation aims to provide a deeper understanding of the interplay among Na₂SO₄ salinity, microbial dynamics, and metabolic response in SADF systems.
2 Materials and methods
2.1 Experimental setup and bioreactor operation
A lab-scale upflow packed-bed SADF bioreactor was constructed with a working volume of 2 L (Supplementary Figure 1). The bioreactor was filled with spherical filler, which was prepared by mixing elemental sulfur and calcium carbonate into a uniform granule with a particle size of 3–5 mm at a ratio of sulfur to alkali source of 3:1 (wt./wt.). The reactor was maintained at a constant temperature of 25°C via a water bath, and the hydraulic retention time (HRT) was controlled at 60 min using a peristaltic pump to ensure consistent flow rates throughout the experiment. The seed sludge used was sourced from the anoxic process of a municipal wastewater treatment plant in Jiangsu Province. Initially, the sludge was subjected to intermittent feeding, followed by a gradual transition to continuous feeding over a period of 2 weeks. During this time, the concentration and organic loading rate progressively increased to acclimate the microbial community. Prior to the formal experimental phase, the system had stably operated for 3 months, ensuring that the seed sludge was well-adapted.
To investigate the effects of Na₂SO₄ salinity stress on denitrification performance, the reactor was operated continuously for 30 days, comprising five stages (S0, S1, S2, S3, and S4), each lasting 6 days. During these stages, the salinity (expressed as total dissolved solids, TDS) was set at 400, 3,000, 6,000, 9,000, and 12,000 mg/L, respectively. Reagent-grade Na₂SO₄ was added to the influent to adjust the salinity levels.
The influent NO3−-N concentration was set to be 20 mg/L. Each liter of synthetic wastewater contained the following components: 145 mg of KNO3, 151.2 mg of NaHCO3, 3.821 mg of NH4Cl, 6.6 mg of K2HPO4 and 12 mg/L of MgCl2. The influent pH was tightly regulated within the range of 7.9 to 8.4 through the combined buffering action of calcium carbonate packing and external supplementation with sodium bicarbonate (NaHCO3).
2.2 Chemical analysis
Water samples were routinely collected from the SADF system for chemical analysis. Subsequently, they were filtered through 0.45 μm membrane filters. A water quality analyzer (LH-3BAV12, Lianhua, China) was employed to determine the total nitrogen, nitrate nitrogen, and nitrite nitrogen. A multi - function portable meter (DZS-708 L, Leici, China) was used to directly measure the pH, dissolved oxygen, conductivity, and total dissolved solids (TDS).
Regarding the description of the filler’s surface morphology, samples were taken from S0 to S4 and fixed with 2.5% glutaraldehyde (pH = 6.8) at 4°C for 1.5 h. After that, the filler samples were rinsed three times with 50 mmol·L−1 phosphate buffer solutions. Then, a series of ethanol washes were conducted at concentrations of 50, 70, 80, and 90%, each wash lasting 15 min. Subsequently, three additional 15-min washes with 100% ethanol were carried out. Next, dehydration was performed using a 1:1 mixture of ethanol and isoamyl acetate for 15 min. Finally, the samples were dehydrated with isoamyl acetate for 15 min. After freeze-drying and gold sputtering, the treated samples were attached to a scanning electron microscope (SEM, FEI Quanta FEG250, USA) for observation.
2.3 Metagenomic sequencing and data analysis
Sludge samples were, respectively, collected from stages S0 to S4 with salinities of 400, 3,000, 6,000, 9,000, and 12,000 mg/L. These samples were immediately frozen using dry ice and then shipped to Novogene (Tianjin, China) for library construction and subsequent metagenomic sequencing. DNA sequencing was carried out using the Illumina Novaseq 6,000 PE150 platform. Metagenomic sequencing for each sample yielded 10 Gb raw data. Clean reads from all samples were pooled, co-assembled and binned using metaWrap v1.3 (Uritskiy et al., 2018), which ultimately resulted in 68 bins. Taxonomy annotation was performed using GTDB-tk v2.2.3 software (Chaumeil et al., 2020), and the relative abundance calculation was calculated with coverm v0.6.1 software.1 The assembled contigs were aligned with the KO (KEGG Orthology) database (Kanehisa et al., 2016) to annotate the functional genes.
Taxonomy abundance plots and heatmaps were created using R v4.3.02 in conjunction with the ggplot2 package. The nitrogen and sulfur metabolic pathways were depicted based on the KEGG pathway. A co-occurrence network was constructed using R v4.3.0 to illustrate the interaction relationships among the major microbial components in the SADF system, and data visualization was accomplished using Gephi v.0.10.1 software.
2.4 Global mapping for SADF-related microorganisms
High-throughput DNA sequencing data from 472 global activated sludge samples were downloaded from the National Center for Biotechnology Information (NCBI) database.3 These samples were sourced from 57 wastewater treatment plants (WWTPs) across 39 cities in 12 countries spanning 4 continents. The assembled bins in this study were aligned with the 472 downloaded global samples using the coverm v0.6.1 software. To analyze the geographic differences, the mapped percentages among the 57 WWTPs were calculated, and a one-way ANOVA test was conducted using SPSS software (v26.0).
3 Results and discussion
3.1 SADF system exhibited robust Na₂SO₄ salinity tolerance
In conventional saline wastewater treatment systems, influent undergoes preliminary desalination via physicochemical pretreatment (e.g., coagulation-flocculation or evaporation) to maintain salinity below the 1.5% (15,000 mg/L TDS) operational threshold for biological processes (Yadai and Suzuki, 2023; Kavitha et al., 2019). To systematically evaluate Na₂SO₄ salinity impacts, we therefore designed a controlled gradient series spanning 400 mg/L (baseline), 3,000 mg/L, 6,000 mg/L, 9,000 mg/L, and 12,000 mg/L, strategically aligned with typical industrial discharge concentrations to bridge laboratory findings and field applications.
As illustrated in Figure 1a, during the initial acclimation phase (S0, Days 1–6), the SADF system achieved rapid start-up under low-salinity conditions (TDS = 400 mg/L), attaining an effluent total nitrogen (TN) concentration as low as 3.4 mg/L by the end of S0. Subsequent incremental salinity increases to 6,000 mg/L in phases S1 and S2 exhibited minimal impact on nitrate removal performance, maintaining average effluent NO₃−-N concentrations of 3.2 mg/L (S1) and 2.4 mg/L (S2). Notably, nitrite accumulation showed progressive elevation from 0.8 mg/L (S0) to 2.3 mg/L (S2), as depicted in Figure 1b. At the maximum salinity condition (S4: TDS = 12,000 mg/L), the system experienced a marked performance decline, characterized by a 35% reduction in total nitrogen (TN) removal efficiency (from 77% at S2 to 50% at S4) (Figure 1a) and a concurrent 17% decrease in nitrate removal efficiency (from 88 to 71%) (Figure 1b). Moreover, nitrite accumulation demonstrated a progressive increase, culminating in an average concentration of 4.5 mg/L (Figure 1b). The sustained accumulation of nitrite may exert inhibitory effects on microbial metabolism, particularly under high salinity conditions, where intensified redox stress could impair the activity of key metabolic enzymes and disrupt the denitrification pathway, ultimately compromising process efficiency (Bian et al., 2024). Moreover, as an intermediate metabolite, elevated nitrite levels may drive shifts in microbial community composition, rendering the system more susceptible to environmental perturbations (Wang T. et al., 2024; Wang X. et al., 2024).
Figure 1. SADF system performance under sulfate salinity stress. (a) Effluent TN concentrations and TN removal rate. (b) Effluent NO3−-N, NO2−-N and NO3−-N removal rate. (c) Influent and effluent pH and salinity.
This observation aligns with findings from previous studies on sulfate-induced inhibition of sulfur-based autotrophic denitrifying bacteria. A recent study (Zhu et al., 2023) demonstrated that denitrification was suppressed and nitrite accumulation occurred when sulfate concentrations exceeded 6,000 mg/L. Notably, another study (Campos et al., 2008) reported an earlier inhibition threshold, observing that denitrification activity began to decline at 1,500 mg/L SO₄2− and was reduced by 85% at 15,000 mg/L SO₄2−. Interestingly, Chung et al. (2014) observed no apparent inhibitory effect on autotrophic denitrification even when sulfate concentrations were over 10,000 mg/L, with complete nitrate removal still achieved under these conditions. The observed discrepancies may arise from variations in experimental design, including the use of elemental sulfur (S0) as the electron donor in our study compared to thiosulfate (S₂O₃2−) in other’s systems and our continuous-flow operation with a HRT of only 60 min, contrasting with their batch-mode operation utilizing flask bottles with a prolonged incubation period of up to 350 h.
Throughout the experimental process, the influent pH remained consistently at an average of 8.18, while the effluent pH remained relatively stable at an average of 6.44 (as shown in Figure 1c). This stability in effluent pH was achieved by supplementing the system with sodium bicarbonate (NaHCO₃) to neutralize the acidity generated during sulfur oxidation.
In general, the SADF system demonstrated exceptional operational resilience under two challenging conditions: (1) rapid Na₂SO₄ salinity escalation (400 to 12,000 mg/L TDS within 30 days) and (2) a comparatively low HRT of 60 min. Despite these stressors, the system maintained a 71–88% nitrate removal efficiency alongside sustained high volumetric nitrate loading rate (0.26–0.32 kg NO₃−-N/m3·d). Given that HRT and prolonged salt-tolerant sludge acclimation are critical operational parameters for bioreactor performance in hypersaline wastewater treatment, we hypothesize that optimizing HRT and gradual microbial adaptation could further enhance TN removal in scenarios involving nitrate-laden wastewater with elevated Na₂SO₄ salinity. These findings not only highlight the significant potential of this approach but also provide a theoretical foundation for optimizing treatment strategies in saline wastewater management, thereby warranting comprehensive investigation through future experimental validation studies. However, despite the good denitrification performance and tolerance of the SADF system under high sulfate concentrations, the elevated sulfate levels may pose potential impacts on downstream treatment units. These impacts include changes in microbial community structure and corresponding alterations in metabolic pathways. Additionally, high sulfate concentrations can accelerate the dissolution of calcium carbonate fillers, damage their structural stability, and weaken their buffering capacity. Therefore, while the SADF system shows promise for treating high-sulfate wastewater, these potential downstream effects should be carefully considered in the design and operation of integrated treatment processes.
3.2 The fluctuation in denitrification performance was attributed to disturbance in the microbial community structure
Through SEM analysis, microbial colonization and biofilm formation on the surface of the filler were observed under different salinity stress conditions. At 500x magnification, the surface of the new, unused filler was uneven, featuring numerous grooves and irregularities (Supplementary Figure 2). After start-up at a salinity of 400 mg/L, the filler surface exhibited depressed cavities and numerous irregular cracks. This phenomenon indicated that after reacting with nitrate, the filler provided a larger specific surface area for biofilm growth, enhancing denitrification and significantly accelerating the biofilm formation rate in wastewater treatment (Zhou et al., 2021).
Upon further magnification of the S0 filler surface to 10,000x, numerous short rod-shaped, long rod-shaped, and spherical bacteria were observed on the filler surface, suggesting a gradually stabilizing biofilm structure. Additionally, some filamentous bacteria were detected on the filler biofilm. These filamentous bacteria could act as a bridge connecting various bacteria, facilitating the colonization of diverse microorganisms on the membrane surface (Secor et al., 2015). With fast elevated salinity, the microbial density on the filler surface also rose. Interestingly, SEM observations revealed that the microbial composition on the surface was similar across different salinities, mainly consisting of rod-shaped and spherical bacteria. However, the number of filamentous bacteria gradually decreased, indicating that high salinity had a certain inhibitory effect on the growth of filamentous bacteria (Quartaroli et al., 2019). For the S3 and S4 filler samples, the biofilm on the surface became much denser, with a further increase in the number of short rod-shaped and spherical bacteria. These results suggested that the sulfur-based filler supported a relatively rich microbial population after continuous operation. Nevertheless, the microbial structure seemed to change significantly with increasing salinity. Meanwhile, a more detailed investigation of the taxonomic changes in the main functional microorganisms and the microbial community structure at each stage of the reactor was required. Subsequently, metagenomic sequencing was employed for downstream microbial analysis at the genomic and genetic levels.
At the phylum level (Figure 2a), Proteobacteria was the most abundant phylum, accounting for 88.2–93.0% of the entire microbial community. This was not unexpected, as many denitrification microorganisms belong to Proteobacteria (Yang et al., 2021). Bacteroidota, the second most abundant phylum, was also found to be salinity-resistant, with its relative abundance increasing from 2.5 to 8.1%. These two phyla have been previously identified as typical phyla involved in autotrophic denitrification processes and were also shown to exhibit selectivity to NaCl salinity (Li Y. Y. et al., 2023; Li Z. R. et al., 2023).
Figure 2. Microbial community shift at five salinity stages at the (a) phylum level and the (b) genus level. (c) Heatmaps of the top 30 abundant bins with standardized Z score corresponding to the relative abundance.
At the genus level (Figure 2b), it was obvious that salinity stress remarkably shaped microbial community structure. At stage S0, the top three dominant genera were Thiobacillus (48.37%), Sulfurimicrobium (30.90%), and Dokdonella_A (6.30%), all of which were denitrifiers. Interestingly, as salinity increased, the relative abundances of SADF-related bacteria showed opposite trends. For example, in our SADF system, there were three main genera that utilized elemental sulfur as an electron donor for sulfur autotrophic denitrification: Thiobacillus, Sulfurimicrobium, and Sulfurimonas. The abundance of Thiobacillus increased from 48.37% (at S0) to 63.79% (at S4), indicating that this genus was highly resistant to salinity stress and played a key role in nitrogen removal in the SADF system. This might be attributed to its salt-tolerance genes (Li R. et al., 2021; Li W. et al., 2021). Notably, regardless of the salinity stress, the dominant genus Thiobacillus maintained the highest abundance, and its proportion was positively correlated with salinity (R2 = 0.84, Supplementary Figure 3). Similar results were also observed in the previous study (Xia Y. et al., 2019; Xia Z. et al., 2019), where both studies noted the changes in denitrifying bacteria under different salinity conditions. Under high salinity conditions, the relative abundance of certain denitrifying bacteria increased, showing strong salt tolerance.
In contrast, the relative abundances of Sulfurimicrobium and Sulfurimonas decreased from 30.90 to 3.23% and from 3.35 to 0.95%, respectively. It could be inferred that these two genera might be more sensitive to salinity stress. Overall, the total abundances of these three genera reached 81.43% at a salinity of 400 mg/L and then rapidly dropped to 66.65% at a salinity of 3,000 mg/L. Surprisingly, their abundances remained relatively stable in the later stages, even showing a slight increase to 67.97% at a salinity of 12,000 mg/L. A high proportion of sulfur-oxidizing bacteria might significantly contribute to the fluctuating yet relatively high nitrogen removal efficiency observed in the SADF system. Additionally, Figures 2a,b indicated that the top 10 most abundant bacteria in the system were all rod-shaped bacteria, while cocci or filamentous bacteria were not commonly detected in any of the five stages. This was slightly different from the SEM results, suggesting that the spherical morphology observed in the SEM images might not be cocci but could be extracellular polymeric substances (EPS) (Zhang et al., 2021).
As shown in Figure 2c, it was clear that the standardized Z score of almost all the top 30 abundant bins underwent significant changes under the influence of salinity. Moreover, a distinct pattern change was observed at stage S3, where salinity-resistant bacteria began to predominate, indicating that the increase in salinity to 9,000 mg/L had a crucial impact on the system (Chen B. et al., 2022; Chen H. et al., 2022). In this study, under low nitrogen loads, the shift in the microbial community structure in response to salinity stress influenced the denitrification performance in the SADF system.
3.3 Competition and cooperation among different microorganisms enhance salinity resilience in the SADF system
To gain deeper insights into the symbiotic relationships among SADF-related microorganisms, co-occurrence network was specifically constructed for 4 abundant autotrophic and heterotrophic denitrifiers, namely Thiobacillus, Sulfurimicrobium, Thermomonas and Dokdonella_A, to illustrate their interactions with other bins. Bins connected to these genera were mainly classified into six phyla, among which Proteobacteria and Bacteroidota were most widely distributed, accounting for over 77% of all nodes.
As the most abundant denitrifier, Thiobacillus was found to have 13 positive and 9 negative correlations with other bins (Figure 3). Notably, Thiobacillus exhibited complex relationships with other denitrifiers. For instance, it was positively correlated with Ferruginibacter and unclassified genera belonging to Ignavibacteriaceae, which could be explained by the salinity resistance of these bins (Figure 2c). Previously, Ferruginibacter was reported to be related to biofilm formation (Xia Y. et al., 2019; Xia Z. et al., 2019) and Ignavibacteriaceae was found to be capable for dissimilatory nitrate reduction to ammonia (DNRA) through nitrite (Lawson et al., 2017). Surprisingly, the other 3 denitrifiers, Dokdonella_A, Sulfurimicrobium and Thermomonas were all negatively correlated with Thiobacillus, possibly due to the competition for the same substrate and ecological niche (Pastore et al., 2021).
Figure 3. Co-occurrence network for 4 major autotrophic and heterotrophic denitrifiers. (Node size represents the abundance of symbiotic bins. Red colour indicates the positive correlation, while the green color for the negative correlation).
Sulfurimicrobium was the second most abundant bacterium at S0, but its abundance rapidly declined under sulfate salinity stress, which was previously identified as an autotroph utilizing thiosulfate, tetrathionate, and elemental sulfur as an electron donor (Kojima et al., 2021). In the co-occurrence network, Sulfurimicrobium was positively correlated with another sulfur-oxidizing denitrifier Sulfurimonas, whose abundance also decreased significantly as salinity increased. This is consistent with the results reported in studies on salinity stress that salinity inhibits sulfur metabolic activity, indicating that these microbial communities generally face the problem of limited electron donor utilization in high-salinity environments (Zhou et al., 2024). Thermomonas and Dokdonella_A were likely to be involved in heterotrophic denitrification in WWTPs. However, Thermomonas was also identified to possess an autotrophic pathway in the mixotrophic denitrification system (Xing et al., 2018). These two heterotrophic genera had a positive relationship in our SADF system, which was consistent with the previous study that Thermomonas and Dokdonella_A were, respectively, responsible for nitrate and nitrite reduction (Wang et al., 2019).
The microbial co-occurrence network for the 4 major denitrifiers in the SADF system reveals intricate patterns of competition and cooperation, which may contribute to the system’s enhanced denitrification stability under salinity stress. The predominance of positive over negative correlations suggests that cooperative interactions among microorganisms may outweigh competitive ones within the SADF system. This tendency may reflect a collective microbial strategy to cope with salinity stress, in accordance with interaction patterns observed in response to temperature fluctuations (Zhang et al., 2022). Such correlations may indicate similar ecological responses to environmental changes. However, these inferences remain speculative and require further validation through experimentally verified gene expression and metabolic dynamics.
3.4 Stimulation of nitrogen and sulfur metabolic pathways under sulfate stress
To explore the microbial similarities among samples with varying salinities, the distribution of gene counts and the similarity in correlation analysis based on genetic composition were investigated accordingly. As depicted in Figure 4, the number of genes increased with rising salinity, reaching a peak at a salinity of 9,000 mg/L. Nevertheless, when the salinity reached 12,000 mg/L, the gene count decreased, suggesting that extremely high salinity may impede gene expression in microorganisms. Surprisingly, a total of 171,972 genes were shared among all samples, accounting for 59.75, 53.84, 49.33, 40.25, and 46.49% of the total genes in each respective sample. This indicated that despite significant variations in salinity, the SADF system largely shared similar functional genes. Regarding unique genes, the five samples had 7,476, 6,764, 5,098, 12,694, and 780 unique genes, respectively. It was observed that the number of unique genes in S4 was significantly lower than that in samples with other salinities, suggesting that high salinity might suppress gene expression. The varying percentages of unique genes could be attributed to shift in the microbial community composition, which may potentially cause fluctuations in the denitrification performance of the SADF system (Wang et al., 2018).
Figure 4. Microbial similarity based on gene patterns at 5 different salinity stages.
In pairwise comparisons, it was notable that S3 and S4 shared a relatively large number of 73,671 genes, while S1 and S2 shared 21,728 genes. Both values were significantly higher than the number of shared genes in other pairs. This result was consistent with the findings in Supplementary Figure 4, where S1 and S2, as well as S3 and S4, exhibited strong correlations. Additionally, the correlation between S0 and other samples was relatively weak, indicating that sulfate salinity had a significant impact on the genetic composition in the SADF system.
As mentioned above, the number of genes escalated with the salinity increase. However, it remains to be resolved whether the abundances of related enzymes in the nitrogen and sulfur metabolic pathways changed in the same direction. Based on the KEGG database, the nitrogen and sulfur metabolic pathways in the SADF system were reconstructed. Among all the genes shown in Figure 5, genes encoding nitrate reductase 1 (EC 1.7.99.-) exhibited the highest relative abundance (1.35–2.16%) in both the nitrogen and sulfur metabolic pathways.
Figure 5. Nitrogen (a) and sulfur (b) metabolic pathways of the SADF system. The scale bar represents the relative abundance. Numbers near the box represent EC numbers for each enzyme.
Regarding the nitrogen metabolic pathway, 12 out of 16 N-related genes identified in S0 showed varying degrees of increase under high salinity stress. The total N-related gene abundance increased from 4.92% in S0 to 7.28% in S4. Interestingly, among various nitrogen metabolic steps, the relative abundance of nitrate reductase genes increased most notably (from 2.94 to 4.25%), suggesting that the conversion of nitrate to nitrite may be particularly sensitive to rapidly increasing salinity (Jeong et al., 2018).
Simultaneously, 7 out of 9 genes related to sulfur metabolic pathway showed an increasing trend in relative abundance in response to salinity stress, with the total abundance rising from 3.03% at S0 to 4.87% at S4. There were mainly two pathways for the transformation of elemental sulfur to sulfate, one generating APS while the other directly produced sulfate from sulfite (Wu et al., 2021). Results showed that at different salinity levels, the total enzyme abundance for the former pathway (1.24–2.38%) was consistently higher than the latter one (0.97–1.55%), suggesting that APS was an important intermediate in the SADF system. Surprisingly, despite the high external sulfate concentrations, the relative abundance of genes involved in the sulfite-to-sulfate conversion pathway also increased, suggesting a potential metabolic adaptation under hypersaline conditions.
Overall, the total relative abundance of genes related to nitrogen and sulfur metabolic pathways exhibited similar upregulated patterns with the fast elevated salinity, with positive correlations identified (R2 = 0.82), indicating that they both played an important role in the SADF system.
3.5 SADF-related microbial community were widely distributed across global wastewater treatment plants
DNA sequencing data were collected from 472 global activated sludge samples sourced from 57 WWTPs located in 39 cities across 12 countries and 4 continents (Figure 6). Raw reads were successively mapped with the assembled SADF-related genomic bins to assess whether the microbial community identified in this study displayed regional characteristics.
Figure 6. Global distribution of SADF-related bins identified in this study across wastewater treatment plants in 39 cities worldwide, calculated with mapped percentage.
Notably, the average mapped percentage for 39 cities was 0.99%, with the highest value reaching 3.52% and the lowest being 0.06%. The highest value of 3.52% were from Zarate, Argentina, and Xiamen, China, respectively. These two locations are geographically distant, and the WWTP in Zarate is an industrial facility for petroleum refinery, while the one in Xiamen is a full-scale municipal WWTP. It is worth noting that a sample from Foshan, China, which is geographically close to Xiamen, had a mapped percentage of only 0.16%. These results strongly suggest that SADF-related microorganisms are widely distributed in global WWTPs.
This finding was further corroborated by statistical analysis. In this study, the 472 samples were classified into 5 groups according to their geographical locations: China, Asia (excluding China), Europe, USA, and Argentina. Intriguingly, one-way ANOVA analysis indicated no statistically significant differences between these groups (p-value > 0.05, see Supplementary Table 1). This implies that geographical location does not significantly affect the relative abundance of the SADF-related microbial community in WWTPs.
Previous studies have identified regional factors such as annual average temperature, humidity, and altitude as major determinants influencing the structure of activated sludge communities (Tian and Wang, 2020; Xia Y. et al., 2019; Xia Z. et al., 2019). However, this study revealed that for a specific subset of functional microorganisms in the activated sludge system, such as the SADF-related microbial community, their distribution remains relatively stable despite regional disparities. In other words, the widespread presence of the SADF-related microbial community in global WWTPs suggests its potential for indigenous enrichment and application of the SADF system on a global scale (Molari et al., 2023; Sriaporn et al., 2021). Previous studies have identified regional factors such as average annual temperature, humidity, and elevation as major determinants of microbial community composition in activated sludge systems (Tian and Wang, 2020; Xia Y. et al., 2019; Xia Z. et al., 2019). However, our findings highlight that certain functionally specialized microbial groups—such as those associated with SADF system—can exhibit relatively consistent global distribution patterns despite environmental and operational heterogeneity. This widespread presence underscores the potential for in situ enrichment and implementation of SADF-based processes across diverse geographic settings (Molari et al., 2023; Sriaporn et al., 2021). Nonetheless, it is important to acknowledge potential limitations that may affect the practical applicability of SADF systems across regions. Variations in wastewater composition, operational parameters (e.g., hydraulic retention time, oxidation–reduction potential), and indigenous microbial communities could influence system performance. Therefore, while the SADF-associated microbiota appears globally prevalent, their functional expression and ecological competitiveness may still be context-dependent, warranting further investigation under site-specific conditions.
4 Conclusion
This study conducted in-depth genomic and genetic analyses to explore the denitrification performance and the shift in the microbial community of the SADF system under different levels of sulfate salinity. The dominant genus, Thiobacillus, was found to increase in abundance as salinity rose. Co-occurrence network analysis indicated a larger number of positive correlations, which contribute to promoting stable denitrification performance. The upregulated stimulations in the nitrogen and sulfur metabolic pathways suggested that they exhibited strong adaptability patterns with the fast elevated salinity. Additionally, the SADF-related microorganisms were found to be widely distributed in global WWTPs, with no statistically significant differences detected among different groups. These findings could provide valuable insights for applications of the SADF system in treating wastewater with high sulfate salinity. Nevertheless, several limitations remain in the present study. First, the findings of our study should be interpreted with caution regarding their representativeness of long-term denitrification processes. Second, it should be noted that the study did not account for the potential influence of complex pollutants (e.g., heavy metals, organic compounds) typically present in real industrial wastewater. Third, the results primarily reflect taxonomic composition and inferred functional potential, rather than offering direct evidence of gene expression or metabolic activity. Future investigations should aim to address these limitations by extending the exposure duration to more accurately simulate long-term denitrification processes, conducting experiments using real industrial wastewater containing complex pollutants to enhance the applicability and reliability of the findings, and incorporating advanced analytical techniques (e.g., transcriptomics, proteomics, or metabolomics) to provide direct evidence of gene expression and metabolic dynamics, thereby offering a more comprehensive understanding of the underlying mechanisms.
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov/sra, accession number: PRJNA1265719.
Author contributions
YiC: Funding acquisition, Conceptualization, Writing – original draft, Methodology. XJ: Visualization, Writing – review & editing, Formal analysis, Investigation. JZ: Formal analysis, Writing – original draft, Investigation. MY: Writing – review & editing. YoC: Funding acquisition, Writing – review & editing, Supervision. HL: Writing – review & editing. YL: Writing – review & editing. FD: Writing – review & editing, Supervision. ZW: Funding acquisition, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. The authors acknowledge the generous support from the National Natural Science Foundation of China (42407190), the National Key R&D Program of China (2023YFE0100900), the Natural Science Foundation of Jiangsu Province (BK20230763) and the Young Elite Scientists Sponsorship Program by JSAST (JSTJ-2023-XH022).
Conflict of interest
YiC, YoC, HL, and YL were employed by Jiangsu Environmental Engineering Technology Co., Ltd. and Jiangsu Environmental Protection Group Co., Ltd.
The remaining 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.
Generative AI statement
The author(s) 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.1615317/full#supplementary-material
Footnotes
References
Abascal, E., Gómez-Coma, L., Ortiz, I., and Ortiz, A. (2022). Global diagnosis of nitrate pollution in groundwater and review of removal technologies. Sci. Total Environ. 810:152233. doi: 10.1016/j.scitotenv.2021.152233
Bai, Y., Wang, S., Zhussupbekova, A., Shvets, I. V., Lee, P.-H., and Zhan, X. (2023). High-rate iron sulfide and sulfur-coupled autotrophic denitrification system: nutrients removal performance and microbial characterization. Water Res. 231:119619. doi: 10.1016/j.watres.2023.119619
Bian, D.-D., Shi, Y.-X., Zhu, X.-R., Sun, X.-L., Xu, X., Ding, L., et al. (2024). Effects of nitrite exposure on the oxidative stress, immune response and intestine microbiota of Procambarus clarkii. Aquaculture 593:741331. doi: 10.1016/j.aquaculture.2024.741331
Campos, J. L., Carvalho, S., Portela, R., Mosquera-Corral, A., and Méndez, R. (2008). Kinetics of denitrification using Sulphur compounds: effects of S/N ratio, endogenous and exogenous compounds. Bioresour. Technol. 99, 1293–1299. doi: 10.1016/j.biortech.2007.02.007
Chatla, A., Almanassra, I. W., Abushawish, A., Laoui, T., Alawadhi, H., Atieh, M. A., et al. (2023). Sulphate removal from aqueous solutions: state-of-the-art technologies and future research trends. Desalination 558:116615. doi: 10.1016/j.desal.2023.116615
Chaumeil, P.-A., Mussig, A. J., Hugenholtz, P., and Parks, D. H. (2020). GTDB-Tk: a toolkit to classify genomes with the genome taxonomy database. Bioinformatics 36, 1925–1927. doi: 10.1093/bioinformatics/btz848
Chen, H., Ma, K., Huang, Y., Fu, Q., Qiu, Y., and Yao, Z. (2022). Significant response of microbial community to increased salinity across wetland ecosystems. Geoderma 415:115778. doi: 10.1016/j.geoderma.2022.115778
Chen, B., Qaisar, M., Wang, K., Li, W., and Cai, J. (2022). Response of simultaneous sulfide and nitrate removal process on acute toxicity of substrate concentration and salinity: single toxicity and combined toxicity. Sci. Total Environ. 836:155639. doi: 10.1016/j.scitotenv.2022.155639
Chung, J., Amin, K., Kim, S., Yoon, S., Kwon, K., and Bae, W. (2014). Autotrophic denitrification of nitrate and nitrite using thiosulfate as an electron donor. Water Res. 58, 169–178. doi: 10.1016/j.watres.2014.03.071
Jeong, D., Cho, K., Lee, C.-H., Lee, S., and Bae, H. (2018). Effects of salinity on nitrification efficiency and bacterial community structure in a nitrifying osmotic membrane bioreactor. Process Biochem. 73, 132–141. doi: 10.1016/j.procbio.2018.08.008
Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M., and Tanabe, M. (2016). KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 44, D457–D462. doi: 10.1093/nar/gkv1070
Kavitha, J., Rajalakshmi, M., Phani, A. R., and Padaki, M. (2019). Pretreatment processes for seawater reverse osmosis desalination systems—a review. J. Water Process Eng. 32:100926. doi: 10.1016/j.jwpe.2019.100926
Kojima, H., Kanda, M., Umezawa, K., and Fukui, M. (2021). Sulfurimicrobium lacus gen. nov., sp. nov., a sulfur oxidizer isolated from lake water, and review of the family Sulfuricellaceae to show that it is not a later synonym of Gallionellaceae. Arch. Microbiol. 203, 317–323. doi: 10.1007/s00203-020-02029-0
Lawson, C. E., Wu, S., Bhattacharjee, A. S., Hamilton, J. J., McMahon, K. D., Goel, R., et al. (2017). Metabolic network analysis reveals microbial community interactions in anammox granules. Nat. Commun. 8:15416. doi: 10.1038/ncomms15416
Li, R., Guan, M., and Wang, W. (2021). Simultaneous arsenite and nitrate removal from simulated groundwater based on pyrrhotite autotrophic denitrification. Water Res. 189:116662. doi: 10.1016/j.watres.2020.116662
Li, Y. Y., Lin, L., and Huang, X. (2023). Partial nitritation-anammox for treatment of saline wastewater: hydrazine-assisted salinity adaptation and nitrate control. Chem. Eng. J. 470:144268. doi: 10.1016/j.cej.2023.144268
Li, W., Liu, J., Zhen, Y., Lin, M., Sui, X., Zhao, W., et al. (2021). Simultaneous removal of nitrite and organics in a biofilm-enhanced high-salt wastewater treatment system via mixotrophic denitrification coupled with sulfate reduction. J. Water Process Eng. 40:101976. doi: 10.1016/j.jwpe.2021.101976
Li, Z. R., Zhang, X.-N., Wang, H.-C., Cheng, H.-Y., Wang, A.-J., Zhang, Y.-Q., et al. (2023). Effects of salinity on sulfur-dominated autotrophic denitrification microorganisms: microbial community succession, key microorganisms and response mechanisms. Chem. Eng. J. 478:147308. doi: 10.1016/j.cej.2023.147308
Liu, X., Beusen, A. H. W., van Grinsven, H. J. M., Wang, J., van Hoek, W. J., Ran, X., et al. (2024). Impact of groundwater nitrogen legacy on water quality. Nat. Sustain. 7, 891–900. doi: 10.1038/s41893-024-01369-9
Molari, M., Hassenrueck, C., Laso-Pérez, R., Wegener, G., Offre, P., Scilipoti, S., et al. (2023). A hydrogenotrophic Sulfurimonas is globally abundant in deep-sea oxygen-saturated hydrothermal plumes. Nat. Microbiol. 8, 651–665. doi: 10.1038/s41564-023-01342-w
Pang, Y., and Wang, J. (2021). Various electron donors for biological nitrate removal: a review. Sci. Total Environ. 794:148699. doi: 10.1016/j.scitotenv.2021.148699
Pastore, A. I., Barabás, G., Bimler, M. D., Mayfield, M. M., and Miller, T. E. (2021). The evolution of niche overlap and competitive differences. Nat. Ecol. Evol. 5, 330–337. doi: 10.1038/s41559-020-01383-y
Picetti, R., Deeney, M., Pastorino, S., Miller, M. R., Shah, A., Leon, D. A., et al. (2022). Nitrate and nitrite contamination in drinking water and cancer risk: a systematic review with meta-analysis. Environ. Res. 210:112988. doi: 10.1016/j.envres.2022.112988
Quartaroli, L., Silva, C. M., Silva, L. C. F., Lima, H. S., de Paula, S. O., Dias, R. S., et al. (2019). Effect of the gradual increase of salt on stability and microbial diversity of granular sludge and ammonia removal. J. Environ. Manag. 248:109273. doi: 10.1016/j.jenvman.2019.109273
Secor, P. R., Sweere, J. M., Michaels, L. A., Malkovskiy, A. V., Lazzareschi, D., Katznelson, E., et al. (2015). Filamentous bacteriophage promote biofilm assembly and function. Cell Host Microbe 18, 549–559. doi: 10.1016/j.chom.2015.10.013
Shao, B., Xie, Y.-G., Zhang, L., Ruan, Y., Liang, B., Zhang, R., et al. (2025). Versatile nitrate-respiring heterotrophs are previously concealed contributors to sulfur cycle. Nat. Commun. 16:1202. doi: 10.1038/s41467-025-56588-1
Shen, Z., Xie, L., Lyu, C., Xu, P., Yuan, Y., Li, X., et al. (2023). Effects of salinity on nitrite and elemental sulfur accumulation in a double short-cut sulfur autotrophic denitrification process. Bioresour. Technol. 369:128432. doi: 10.1016/j.biortech.2022.128432
Sriaporn, C., Campbell, K. A., Van Kranendonk, M. J., and Handley, K. M. (2021). Genomic adaptations enabling Acidithiobacillus distribution across wide-ranging hot spring temperatures and pHs. Microbiome 9:135. doi: 10.1186/s40168-021-01090-1
Sun, H., Wu, Q., Yu, P., Zhang, L., Ye, L., Zhang, X.-X., et al. (2017). Denitrification using excess activated sludge as carbon source: performance and the microbial community dynamics. Bioresour. Technol. 238, 624–632. doi: 10.1016/j.biortech.2017.04.105
Tian, L., and Wang, L. (2020). A meta-analysis of microbial community structures and associated metabolic potential of municipal wastewater treatment plants in global scope. Environ. Pollut. 263:114598. doi: 10.1016/j.envpol.2020.114598
Uritskiy, G. V., DiRuggiero, J., and Taylor, J. (2018). MetaWRAP—a flexible pipeline for genome-resolved metagenomic data analysis. Microbiome 6:158. doi: 10.1186/s40168-018-0541-1
Wang, X., Cheng, L., Tan, Q., Xing, Y., Wu, H., Tian, Q., et al. (2024). Lag in pollutant removal efficiency relative to microbial community dynamics in a denitrifying phosphorus removal bioreactor. Chem. Eng. J. 497:154498. doi: 10.1016/j.cej.2024.154498
Wang, S.-S., Cheng, H.-Y., Zhang, H., Su, S.-G., Sun, Y.-L., Wang, H.-C., et al. (2021). Sulfur autotrophic denitrification filter and heterotrophic denitrification filter: comparison on denitrification performance, hydrodynamic characteristics and operating cost. Environ. Res. 197:111029. doi: 10.1016/j.envres.2021.111029
Wang, D., Li, T., Huang, K., He, X., and Zhang, X.-X. (2019). Roles and correlations of functional bacteria and genes in the start-up of simultaneous anammox and denitrification system for enhanced nitrogen removal. Sci. Total Environ. 655, 1355–1363. doi: 10.1016/j.scitotenv.2018.11.321
Wang, T., Li, X., Wang, H., Xue, G., Zhou, M., Ran, X., et al. (2023). Sulfur autotrophic denitrification as an efficient nitrogen removals method for wastewater treatment towards lower organic requirement: a review. Water Res. 245:120569. doi: 10.1016/j.watres.2023.120569
Wang, L., Shang, S., Liu, W., She, D., Hu, W., and Liu, Y. (2023). Hydrodynamic controls on nitrogen distribution and removal in aquatic ecosystems. Water Res. 242:120257. doi: 10.1016/j.watres.2023.120257
Wang, T., Wang, H., Ran, X., and Wang, Y. (2024). Salt stimulates sulfide−driven autotrophic denitrification: microbial network and metagenomics analyses. Water Res. 257:121742. doi: 10.1016/j.watres.2024.121742
Wang, L., Zhang, J., Li, H., Yang, H., Peng, C., Peng, Z., et al. (2018). Shift in the microbial community composition of surface water and sediment along an urban river. Sci. Total Environ. 627, 600–612. doi: 10.1016/j.scitotenv.2018.01.203
Wu, B., Liu, F., Fang, W., Yang, T., Chen, G.-H., He, Z., et al. (2021). Microbial sulfur metabolism and environmental implications. Sci. Total Environ. 778:146085. doi: 10.1016/j.scitotenv.2021.146085
Xi, H., Zhou, X., Arslan, M., Luo, Z., Wei, J., Wu, Z., et al. (2022). Heterotrophic nitrification and aerobic denitrification process: promising but a long way to go in the wastewater treatment. Sci. Total Environ. 805:150212. doi: 10.1016/j.scitotenv.2021.150212
Xia, Y., He, X., Feng, Z., Zhang, Q., and Yang, H. (2019). A comprehensive analysis of the microbial diversity in natural and engineered ecosystems based on high-throughput sequencing of 16S rRNA gene. Int. Biodeterior. Biodegradation 140, 160–168. doi: 10.1016/j.ibiod.2019.03.018
Xia, Z., Wang, Q., She, Z., Gao, M., Zhao, Y., Guo, L., et al. (2019). Nitrogen removal pathway and dynamics of microbial community with the increase of salinity in simultaneous nitrification and denitrification process. Sci. Total Environ. 697:134047. doi: 10.1016/j.scitotenv.2019.134047
Xing, W., Li, J., Li, D., Hu, J., Deng, S., Cui, Y., et al. (2018). Stable-isotope probing reveals the activity and function of autotrophic and heterotrophic denitrifiers in nitrate removal from organic-limited wastewater. Environ. Sci. Technol. 52, 7867–7875. doi: 10.1021/acs.est.8b01993
Xu, X., Zou, Y., Pan, H., Zhang, R., and Gu, B. (2025). Safeguarding groundwater nitrate within regional boundaries in China. Environ. Sci. Technol. 59, 467–477. doi: 10.1021/acs.est.4c08197
Yadai, T., and Suzuki, Y. (2023). Development of softening and ballasted flocculation as a pretreatment process for seawater desalination through a reverse osmosis membrane. NPJ Clean Water 6:7. doi: 10.1038/s41545-023-00226-0
Yang, Z., Zhou, Q., Sun, H., Jia, L., Zhao, L., and Wu, W. (2021). Metagenomic analyses of microbial structure and metabolic pathway in solid-phase denitrification systems for advanced nitrogen removal of wastewater treatment plant effluent: a pilot-scale study. Water Res. 196:117067. doi: 10.1016/j.watres.2021.117067
Ye, J., Ren, G., Liu, L., Zhang, D., Zeng, R. J., van Loosdrecht, M. C. M., et al. (2024). Wastewater denitrification driven by mechanical energy through cellular piezo-sensitization. Nat. Water 2, 531–540. doi: 10.1038/s44221-024-00253-2
Yue, F.-J., Li, S.-L., Waldron, S., Oliver, D. M., Chen, X., Li, P., et al. (2023). Source availability and hydrological connectivity determined nitrate-discharge relationships during rainfall events in karst catchment as revealed by high-frequency nitrate sensing. Water Res. 231:119616. doi: 10.1016/j.watres.2023.119616
Zhang, N., Sun, Y. L., Yao, B. M., Zhang, B., and Cheng, H. Y. (2022). Insight into the shaping of microbial communities in element sulfur-based denitrification at different temperatures. Environ. Res. 215:114348. doi: 10.1016/j.envres.2022.114348
Zhang, C., Yin, L., Ou, Y., Yang, G., Huang, L., and Li, F. (2021). Contribution of selective bacterial extracellular polymeric substances to the polymorphism and morphologies of formed ca/mg carbonates. Int. Biodeterior. Biodegradation 160:105213. doi: 10.1016/j.ibiod.2021.105213
Zhao, Q., Wang, L., Jia, T., Li, X., Zhang, Q., and Peng, Y. (2024). Elemental sulfur–siderite composite filler empowers sustainable tertiary treatment of municipal wastewater even at an ultra-low temperature of 7.3 °C. Nat. Water 2, 782–792. doi: 10.1038/s44221-024-00285-8
Zheng, R., Zhang, K., Kong, L., and Liu, S. (2024). Research progress and prospect of low-carbon biological technology for nitrate removal in wastewater treatment. Front. Environ. Sci. Eng. 18:80. doi: 10.1007/s11783-024-1840-3
Zhou, Y., Chen, F., Chen, N., Peng, T., Dong, S., and Feng, C. (2021). Denitrification performance and mechanism of biofilter constructed with sulfur autotrophic denitrification composite filler in engineering application. Bioresour. Technol. 340:125699. doi: 10.1016/j.biortech.2021.125699
Zhou, Y., Chen, Z., Lu, M., Wan, X., Nong, R., Liu, X., et al. (2024). Effect of salinity on the operational efficiency, microbial community, and metabolism of a bioelectrochemical-granular sludge coupled system for treating sulfate wastewater. J. Water Process Eng. 57:104737. doi: 10.1016/j.jwpe.2023.104737
Keywords: sulfur-based autotrophic denitrification filter, sulfate salinity, Thiobacillus, salinity stress, microbial interactions
Citation: Chen Y, Jiang X, Zhao J, Yang M, Chen Y, Ling H, Liu Y, Deng F and Wang Z (2025) Microbial response under sulfate stress in a sulfur-based autotrophic denitrification system. Front. Microbiol. 16:1615317. doi: 10.3389/fmicb.2025.1615317
Edited by:
Bo Li, Guangdong University of Technology, ChinaReviewed by:
Wei Li, East China University of Science and Technology, ChinaMingyi Xu, South China University of Technology, China
Bin Ji, Wuhan University of Science and Technology, China
Jun Li, Nanjing University, China
Copyright © 2025 Chen, Jiang, Zhao, Yang, Chen, Ling, Liu, Deng 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: Zhu Wang, d2FuZ3podUBnemh1LmVkdS5jbg==; Yong Chen, eW9uZy5jaGVuQGhvdG1haWwuY29t