Isolation, screening, and identification of functional bacteria
Eighteen psychrotolerant nitrogen-removing bacterial strains were isolated and preliminarily screened from activated sludge, river sediments, and frozen soil in northern China during winter when temperatures were consistently below 5°C. These strains exhibited capabilities for removing NH4+-N, NO3−-N, and TN. After re-screening, three strains each of nitrifying and aerobic denitrifying bacteria were obtained. The colony morphology and cellular characteristics of these six strains are detailed in Table 1 and Fig. S1.
Table 1
Colonial and cellular morphological characteristics of six isolated psychrotolerant bacteria
Function | Strain | Colony morphology | Gram staining | Cell morphology | GenBank Accession No. |
|---|
Nitrification | HN1 | Light yellow and transparent, circular morphology with irregular edges, smooth and moist surface with a center protrusion, and 2–3 mm in diameter. | G– | Straight long rod, (1.5–3.0) µm × (0.6–0.8) µm | PP913934 |
HN2 | Slightly yellow and transparent, circular morphology with regular edges, smooth and moist surface with a center protrusion, and 3–4 mm in diameter. | G– | Straight rod, (1.0–2.5) µm × (0.6–0.8) µm | PP911455 |
HN3 | Light yellow and semi-transparent, irregular in morphology, smooth, moist and glossy surface with a center protrusion, and 1–2 mm in diameter. | G– | Curved rod, (1.0–2.0) µm × (0.8–1.0) µm | PP911456 |
Denitrification | AD1 | White and transparent, circular morphology with regular edges, smooth and moist surface with a center high protrusion, and 2–4 mm in diameter. | G– | Short rod, (0.6–0.8) µm × (0.6–0.8) µm | PP913935 |
AD2 | Light yellow and transparent, irregular in morphology, smooth and moist surface with a center protrusion, and 2–3 mm in diameter. | G– | Straight long rod, (2.0–3.0) µm × (0.6–0.8) µm | PP911450 |
AD3 | White and opaque, circular morphology with regular edges, smooth and moist surface with a center high protrusion, and 2–3 mm in diameter. | G– | Curved rod, (1.5–2.5) µm × (0.8–1.0) µm | PP913933 |
Three screened nitrifying strains were taxonomically identified and designated as Pseudomonas veronii HN1, P. poae HN2, and P. peli HN3 based on the results of 16S rRNA gene sequence alignment against the GenBank database and phylogenetic tree analysis (Fig. S2A-C). After culture in NM at 10°C for 3 days with shaking at 150 rpm, these isolates exhibited significant removal of NH4+-N, achieving REs of 86.8 ± 1.7%, 85.7 ± 1.8%, and 86.4 ± 2.1%, with average removal rate of 1.27 ± 0.15, 1.25 ± 0.13, and 1.26 ± 0.16 mg/L/h, respectively, higher than that of B. simplex H-b (0.74 mg/L/h) at 10°C [10]. The Pseudomonas species have been extensively documented in previous studies for their significant nitrification capabilities under diverse environmental settings. Yang et al. [27] isolated a novel acid-resistant bacterium, P. citronellolis YN-21, which exhibited exceptional heterotrophic nitrification ability under acidic conditions with NH4+-N RR of 7.84 mg/L/h at pH 5.0. P. aeruginosa WS-03 could effectively remove NH4+-N at 30°C with an RR of 8.96 mg/L/h [28]. According to several reports, Pseudomonas species, characterized by exceptional growth rates, remarkable environmental adaptability, and superior NR efficiency, exhibit significant competitive advantages over conventional nitrifying bacteria [29]. Therefore, we selected these three isolated Pseudomonas species as candidates for CMA construction.
Isolated aerobic denitrifying strains exhibited the highest sequence identity (≥ 99%) with Aeromonas sp. strain J223, P. extremaustralis 14 − 3, and Serratia liquefaciens strain ATCC 27592 in the GenBank database (Fig. S2D-F). Thus, the three strains were identified and designated as Aeromonas sp. AD1, P. extremaustralis AD2, and Serratia liquefaciens AD3, which exhibited remarkable NO3−-N removal performance during 3-day cultivation in DM at 10°C and shaking at 150 rpm, with corresponding REs of 85.7 ± 1.6%, 86.3 ± 1.8%, and 85.8 ± 1.8%. Furthermore, the average RR of NO3−-N by these three strains were 1.25 ± 0.21, 1.26 ± 0.17, and 1.25 ± 0.15 mg/L/h, respectively, surpassing that of Acinetobacter tandoii MZ-5 (1.04 mg/L/h at 25°C) [30] and Arthrobacter arilaitensis Y-10 (0.35 mg/L/h at 15°C) [31]. During cultivation, all three bacterial strains accumulated detectable levels of NO2−-N, with concentrations of 2.49 ± 0.37 mg/L (Aeromonas sp. AD1), 2.25 ± 0.41 mg/L (P. extremaustralis AD2), and 2.73 ± 0.46 mg/L (S. liquefaciens AD3). Chen et al. [32] demonstrated that Aeromonas sp. HN-02, an HNAD bacterium, could maintain activity at 5°C, with RRs of 0.9 and 22.3 mg/L/h for ammonia and COD, respectively. P. extremaustralis, a psychrophilic bacterium native to the extreme environments in Antarctica, has been reported as an antibacterial agent [33] and a phenacetin-degrading strain [34]; however, its potential application in NR, especially under low-temperature conditions, has received limited attention. In addition, this study potentially represents the first documented evidence of the denitrification capability of S. liquefaciens, revealing a previously uncharacterized metabolic trait of this species.
Construction of psychrotolerant CMA
Pairwise antagonism assays conducted on the six screened bacterial isolates did not reveal any detectable antagonistic interactions, which indicated that these strains did not inhibit the growth of each other and can be cocultured for the preparation of CMA. Therefore, the CMA developed in this study comprised a psychrotolerant nitrifying consortium NC1 (including P. veronii HN1, P. poae HN2, and P. peli HN3) and an aerobic denitrifying consortium DC1 (including Aeromonas sp. AD1, P. extremaustralis AD2, and S. liquefaciens AD3). The results of comparative experiment illustrated in Fig. 1 indicate that the CMA achieved REs of 89.3 ± 1.8%, 88.3 ± 2.2%, 85.8 ± 1.3%, and 95.5 ± 2.7% for NH4+-N, NO3−-N, TN, and COD, respectively, in CBM, which were obviously improved compared with those of the individual strains and single-functional consortia cultivated under identical conditions. Furthermore, the CMA exhibited the highest OD600 of 1.99 ± 0.03 among different strain combinations. These results indicate that when phylogenetically distinct microorganisms coexist in the same environment but occupy different ecological niches, they can establish mutualistic, symbiotic relationships and leverage the synergistic metabolic interactions to improve microbial proliferation and contaminant removal efficacy [35]. Jing et al. [36] reported similar findings, demonstrating that the composite microbial consortium ECT, formed by combining three screened petroleum-degrading strains, Ralstonia sp. CP, Rhizobium sp. BX, and Acinetobacter sp. FL, achieved superior petroleum-degrading efficiency compared with individual strains. During CMA-mediated NR, the lowest concentration of accumulated NO2−-N was 0.63 ± 0.37 mg/L, indicating complete denitrification. Overall, the CMA exhibited excellent NR performance under low-temperature conditions, which indicated that the two functional bacterial consortia could be cocultured for synergistic elimination of nitrogen. Consequently, the cocultivation-derived CMA, integrating both nitrifying and denitrifying consortia, was selected for subsequent studies.
The biomass growth of CMA and its NR performance at varying inoculation ratios of nitrifying and denitrifying consortia after 3 days of cultivation at 10°C with shaking at 150 rpm are depicted in Fig. 2. At 1:1 inoculation ratio of NC1 and DC1 consortia, the CMA achieved the highest REs of 89.3 ± 1.5% for NH4+-N, 88.1 ± 1.8% for NO3−-N, 85.5 ± 2.1% for TN, and 95.3 ± 2.3% for COD. Concurrently, it exhibited maximal biomass growth (OD600, 1.98 ± 0.03) and minimal NO2−-N accumulation (0.68 ± 0.12 mg/L). The inoculation ratio significantly affected the growth of CMA and its NR efficiency. As reported by Guo et al. [37], the inoculation ratio determined the initial relative abundance of functionally different microorganisms, thereby influencing synergistic metabolic interactions within the cocultured microorganisms. Ebadi et al. [38] proposed that the decontamination performance of composite bacterial consortium is regulated by a dynamic equilibrium involving competitive relationships for nutrients (e.g., carbon, nitrogen, and phosphorus) and synergistic metabolic interactions among functionally distinct microorganisms. Notably, interspecies synergistic metabolism played a decisive role in the overall decontamination capacity of the system. Interspecies synergistic metabolism can be fully exploited only when functionally diverse bacterial strains are combined at their optimal inoculation ratios, thereby maximizing the RE of pollutants [36]. Therefore, the CMA cultured under the optimal inoculation ratio was designated NDC-6.
Effects of carbon sources on NR
The biomass growth and NR performance of NDC-6 using different carbon sources were evaluated, as shown in Fig. 3. When sodium succinate was used as the sole carbon source, with incubation at 10°C for 3 days and shaking at 150 rpm, NDC-6 achieved peak values for both biomass growth (OD600, 1.99 ± 0.01) and TN-removal efficiency (87.8 ± 1.1%). When cultured with sodium acetate or sodium citrate as sole carbon sources, NDC-6 maintained secondary biomass amount and NR efficacy relative to that achieved with the optimal carbon source, exhibiting OD600 values of 1.93 ± 0.01 and 1.97 ± 0.01, and TN-removal efficiencies of 79.0 ± 1.6% and 83.2 ± 1.4%, respectively. Relative to other carbon sources, NDC-6 exhibited the poorest growth (OD600: sucrose, 0.39 ± 0.01; glucose, 0.49 ± 0.01) and TN removal (sucrose 37.7 ± 1.3%; glucose 52.0 ± 2.1%) when cultured with sucrose or glucose as a carbon source. Numerous studies have indicated that sodium succinate and sodium citrate—intermediate metabolites in the tricarboxylic acid cycle during bacterial respiration—were more readily utilized by nitrifying or denitrifying bacteria and enhanced the activity of nitrate reductase [39]. Moreover, Wei et al. [40] demonstrated that sodium acetate, a small molecule organic compound with a simple metabolic pathway, induced higher efficiency of electron generation, transfer, and competition, thereby promoting rapid bacterial proliferation and complete denitrification. Collectively, sodium succinate, sodium citrate, and sodium acetate have been identified as optimal carbon sources for various nitrifying and denitrifying bacteria. For instance, the use of sodium succinate as the carbon source yielded the highest NH4+-N-removal efficiency in HNAD mixed bacteria HY-1 [41]. Sodium citrate served as an optimal carbon source for Acinetobacter calcoaceticus TY1 [42]. Wang et al. [43] reported significant enhancement of the denitrification efficiency of Paracoccus versutus JUST-3 by sodium acetate. In contrast, when macromolecular organic compounds, such as glucose or sucrose, were used as sole carbon source, they underwent hydrolysis into small organic acids prior to microbial utilization. This process reduced the carbon assimilation efficiency, inhibited bacterial growth and proliferation, and consequently compromised the NR performance. Similar findings were reported in studies on Rhizobium sp. WS7 [44] and Thauera linaloolentis [40] wherein sucrose was more unfavorable for denitrification. Therefore, based on its maximal promotion of NDC-6 growth and NR synthesis, sodium succinate was identified as the preferred carbon source and employed in subsequent experiments.
Optimization of culture conditions for NR by NDC-6
Twenty-nine experimental Box–Behnken design matrices for independent variables (C/N ratio, temperature, pH, and shaking speed) and corresponding response variable (Y) for TN-removal efficiency are presented in Table S4. Following second-order polynomial regression fitting, the quadratic polynomial regression equation in terms of coded independent variables was obtained (Eq. (3)).
$$\:\text{Y}\text{}\text{=}\text{}\text{84.4}\text{-1.59}\text{A}\text{+1.31}\text{B}\text{-5.05}\text{C}\text{-1.64}\text{D}\text{+0.72}\text{AB}\text{-0.32}\text{AC}\text{-0.18}\text{AD}\text{+0.18}\text{BC}\text{+0.83}\text{BD}$$
$$\:\text{-1.01CD-2.21}{\text{A}}^{\text{2}}\text{-7.28}{\text{B}}^{\text{2}}\text{-8.26}{\text{C}}^{\text{2}}\text{-3.65}{\text{D}}^{\text{2}}$$
3
As presented in Table 2, the significance test value (F-value, 208.26) and low probability value (p < 0.0001) for the established quadratic regression model in ANOVA indicated that the model could effectively describe the impact of various factors on the TN-removal efficiency. The determination coefficient R² and adjusted R2 were 0.987 and 0.983, respectively, which were close to the predicted R² of 0.983. The model exhibited nonsignificant lack of fit (p > 0.05), indicating negligible pure error and good agreement between the model-predicted and experimental values [45]. Collectively, these findings indicated that the fitted quadratic polynomial model for TN removal by NDC-6 had high accuracy and reliability in predicting the optimal values for the independent variables and analyzing the interactive effects of these variables on the TN-removal efficiency.
Table 2
ANOVA of the fitted quadratic polynomial model for TN removal by NDC-6
Source | Sum of squares | Mean square | F value | P value (prob > F) | Significance |
|---|
Model | 995.1 | 71.1 | 73.7 | < 0.0001 | ** |
A-C/N | 30.8 | 30.8 | 31.9 | < 0.0001 | ** |
B-Temperature | 26.3 | 26.3 | 23.3 | < 0.0001 | ** |
C-pH | 258 | 258 | 267 | < 0.0001 | ** |
D-Shaking speed | 23.0 | 23.0 | 27.9 | 0.0002 | ** |
AB | 2.06 | 2.06 | 2.13 | 0.166 | |
AC | 0.41 | 0.41 | 0.42 | 0.528 | |
AD | 0.13 | 0.13 | 0.13 | 0723 | |
BC | 0.12 | 0.12 | 0.13 | 0.727 | |
BD | 7.08 | 7.08 | 7.33 | 0.017 | * |
CD | 0.25 | 0.25 | 0.25 | 0.622 | |
A2 | 33.9 | 33.9 | 35.2 | < 0.0001 | ** |
B2 | 376 | 376 | 390 | < 0.0001 | ** |
C2 | 386. | 386 | 400 | < 0.0001 | ** |
D2 | 73.2 | 73.2 | 75.9 | < 0.0001 | ** |
Residual | 13.5 | 0.96 | | | |
Lack of Fit | 10.7 | 1.07 | 1.55 | 0.358 | not significant |
Pure Error | 2.77 | 0.69 | | | |
Cor Total | 1008.5 | | | | |
| Note: ** very significant, P value<0.01; * significant, P value<0.05 |
The ANOVA results revealed that all linear terms (A, B, C, and D) in the established model exhibited p-values less than 0.01, signifying that each of the four factors significantly affected the TN removal by NDC-6. This result was consistent with previous reports that C/N ratio, temperature, pH, and DO are critical environmental factors for BNR [46]. Using the F-value analysis, the factors influencing TN removal were ranked according to their influence degree as follows: pH > shaking speed > C/N ratio > temperature. pH had the most significant impact on TN removal. Similar findings were reported for the optimization of conditions for coculture of B. cereus G2 and B. pumilus G5 by Peng et al. [47], who demonstrated that excessively acidic or alkaline environments remarkably compromised the NR performance of CMA. Notably, temperature had a relatively weaker influence on TN removal compared with the other three factors. The reason might be that NDC-6 was composed of multiple psychrotolerant strains, which had a certain tolerance to low-temperature environments and a broad temperature-adaptation range. In addition, temperature and shaking speed significant influenced TN removal (p < 0.05), signifying that the NR performance of NDC-6 was susceptible to the interactions of these two environmental parameters.
The 3D response surface plots illustrating the interactive effects of the four independent variables (C/N ratio, temperature, pH, and shaking speed) on the removal of TN by NDC-6 are presented in Fig. 4. With increasing C/N ratio and temperature, the TN-removal efficiency by NDC-6 initially increased and subsequently decreased (Fig. 4A). The highest TN-removal efficiencies were achieved at a C/N ratio of approximately 5–6 and a temperature around 10°C. Carbon sources serve as electron donors essential for denitrification [48]. Carbon deficiency triggers competitive inhibition between heterotrophic nitrifying bacteria and aerobic denitrifying bacteria for limited carbon resources, causing concurrent decline in NH4+-N- and NO3−-N-removal efficiencies [49]. In contrast, excessive carbon availability stimulates the overproliferation of heterotrophic bacteria, which consumes substantial amounts of DO, resulting in the deterioration of water quality and carbon wastage. Therefore, effective enhancement of NR can only be achieved when the C/N ratio is maintained within an optimal range. Temperature stress induces microbial oxidative stress and cellular damage, thereby impairing the nutrient-removal capacity [43]. At lower temperatures, enzymatic activity is inhibited, reducing metabolic rates and substantially diminishing the NR performance [50]. On the contrary, excessively high temperatures induce enzyme denaturation [51]. Optimal metabolic rates can only be maintained within a defined temperature range to enhance NR.
As illustrated in Fig. 4F with the increase in initial pH of CBM and shaking speed of the oscillating incubator, the efficiency of TN removal by NDC-6 exhibited parabolic variation trends, characterized by an initial increase followed by a decrease. The highest efficiency of TN removal was achieved at initial pH in the range of 7.0–7.5 and a shaking speed of 150–160 rpm. Similar to other NR-functional strains, such as Paracoccus versutus JUST-3 [43] and Glutamicibacter halophytocola MD1 [52], NDC-6 preferred a neutral or weakly alkaline environment (pH 6–9). Exposure to strongly acidic or alkaline conditions induces aggregation of surface charge on bacterial cells, disrupting the integrity of cell structure and inhibiting microbial proliferation capacity and enzymatic activity, thereby impairing nitrogen transformation [53]. Consequently, precise pH regulation is a critical determinant in optimizing microbial NR.
In this study, DO level in the culture medium was adjusted by regulating the shaking speed of the oscillating incubator. Jin et al. [54] found that the DO content increased linearly with shaking speed within a certain range. Preliminary experiments in this study indicated that the DO content in the culture medium in 250 mL conical flask was approximately 3.3 mg/L at a shaking speed of 100 rpm. Each 10 rpm increment in shaking speed increased the DO content by approximately 0.4 mg/L. The highest TN-removal efficiency was observed for the shaking speed of 150–160 rpm, corresponding to a DO content of 5.3–5.7 mg/L. Li et al. [55] reported that excessively low DO content inhibited the activity of ammonia monooxygenase in nitrifying bacteria, reducing the nitrification performance. Moreover, low DO concentrations suppress the catalytic efficiency of denitrifying enzymes, thereby slowing down mass transfer and substrate utilization [56]. Conversely, an excessive high DO content not only inhibits the activity of nitrite reductase in the denitrification process but also induces competition for electrons between O2 and NO3−-N, leading to incomplete nitrate reduction [57]. This highlights the importance of optimizing the shaking speed to ensure efficient removal of nitrogen by NDC-6.
The cultivation conditions optimized for maximizing TN removal were as follows: A = 0, B = 0.075, C = − 0.3, and D = − 0.2. Accordingly, the model predicted that maximum TN-removal efficiencies of 90.2% would be achieved at a C/N ratio of 6, a temperature of 10.2°C, pH 7.2, and a shaking speed of 156 rpm. Under the optimal conditions, a 3-day cultivation in the verification experiment yielded a TN-removal efficiency of 89.8 ± 1.8%%, which closely approached the theoretically predicted maximum RE for TN. This confirmed the capability of the model to accurately optimize cultivation conditions.
Nitrogen metabolic pathways in NDC-6
Nitrogen balance analysis
The single functional strains, heterotrophic nitrifying consortium, aerobic denitrifying consortium, and NDC-6 were inoculated into optimized CBM (containing 52.5 ± 0.2 mg/L NH4+-N and 52.6 ± 0.2 mg/L NO3−-N) at optimal inoculation ratios. Following 3-day incubation at 10.2°C with shaking at 156 rpm, different forms of nitrogen were determined to establish nitrogen mass balance (Table 3). Notably, a distinct accumulation of NO2−-N (> 5 mg/L) was consistently observed during NR by P. veronii HN1, P. poae HN2, P. peli HN3, and the nitrifying consortium NC1, suggestive of the absence of significant denitrification capability of these strains. Furthermore, only the individual denitrifier, denitrifying consortium, and NDC-6 were capable of converting TN into gaseous-N. Of the initial TN content, 40.2 ± 1.6% was converted to gaseous-N by NDC-6, which was significantly higher than that in the case of Aeromonas sp. AD1 (25.5 ± 1.4%), P. extremaustralis AD2 (28.2 ± 1.8%), S. liquefaciens AD3 (24.9 ± 2.1%), and the denitrifying consortium DC1 (28.4 ± 1.2%). Gas chromatography revealed N2O concentrations of 0.26 ± 0.02 mg/L for Aeromonas sp. AD1, 0.33 ± 0.02 mg/L for P. extremaustralis AD2, 0.29 ± 0.02 mg/L for S. liquefaciens AD3, 0.15 ± 0.01 mg/L for DC1, and 0.07 ± 0.01 mg/L for NDC-6, indicating that the gaseous-N produced during denitrification by these strains was predominantly N2 instead of N2O. These observations indicate enhanced completeness of denitrification and metabolic efficiency by NDC-6, wherein a higher proportion of NO3−-N or NO2−-N was fully reduced to the terminal product N2 rather than accumulating as intermediate gases (e.g., N2O) or residual ions. Similar results were reported by Fang et al. [58], who found that coculturing the yeast Kazachstania exigua T14-1 with the bacterium Methylobacterium sp. T5-6 substantially improved the nitrogen metabolism capacity compared with monocultures of either strain. Furthermore, following 3 days of cultivation at 10.2°C, > 30 mg/L of intracell-N was detected in all the tested microbial combinations, indicating that > 30% of the initial TN content was assimilated into cellular nitrogen via biosynthetic assimilation. Among these, NDC-6 exhibited the highest intracell-N conversion rate at 49.5 ± 1.4%. In contrast to the reported microbial NR pathways primarily relying on assimilation [59], NDC-6 achieved NR under low-temperature aerobic conditions primarily via dual pathways of bacterial assimilation and dissimilation converting inorganic N (mainly NH4+-N and NO3−-N) into intracell-N and N2. These findings indicated that the composite consortium NDC-6 exhibited exceptional NR performance under low-temperature conditions, attributable to synergistic interactions among constituent strains that significantly enhance nitrogen metabolism efficiency, highlighting its greater potential for application in wastewater treatment.
Table 3
Nitrogen balance analysis during nitrogen removal by different strain combinations
Strain combinations | Nitrogen concentration (mg/L) | RETN (%) |
|---|
NH4+-N | NO3―-N | NO2―-N | Organic-N | Intracell-N | Gaseous-N | TN |
|---|
Pseudomonas veronii HN1 | 7.21 ± 0.34 | 50.67 ± 1.65 | 7.64 ± 0.52 | 1.53 ± 0.31 | 37.79 ± 1.05 | – | 67.29 ± 1.53 | 36.19 ± 1.28 |
Pseudomonas poae HN2 | 7.67 ± 0.42 | 51.12 ± 1.78 | 5.85 ± 1.03 | 1.78 ± 0.13 | 38.66 ± 0.92 | – | 66.42 ± 1.24 | 37.01 ± 1.16 |
Pseudomonas peli HN3 | 7.13 ± 0.38 | 50.08 ± 2.14 | 7.35 ± 1.05 | 1.94 ± 0.21 | 38.58 ± 1.07 | – | 66.50 ± 1.56 | 36.94 ± 1.03 |
Aeromonas sp. AD1 | 31.5 ± 1.03 | 7.51 ± 0.36 | 2.24 ± 0.87 | 2.51 ± 0.22 | 34.43 ± 1.21 | 26.89 ± 1.59 | 43.76 ± 1.25 | 58.52 ± 2.05 |
Pseudomonas extremaustralis AD2 | 29.81 ± 1.22 | 7.26 ± 0.13 | 2.06 ± 0.73 | 2.27 ± 0.46 | 33.91 ± 1.06 | 29.77 ± 2.31 | 41.40 ± 1.37 | 60.73 ± 2.14 |
Serratia liquefaciens AD3 | 31.64 ± 1.74 | 7.68 ± 0.29 | 2.78 ± 0.46 | 2.64 ± 0.38 | 34.07 ± 1.23 | 26.27 ± 1.69 | 44.74 ± 1.59 | 57.57 ± 1.98 |
Nitrifying consortium NC1 | 6.77 ± 0.09 | 50.48 ± 1.71 | 5.99 ± 1.22 | 2.08 ± 0.18 | 40.31 ± 1.02 | – | 65.32 ± 2.11 | 38.06 ± 1.56 |
Denitrifying consortium DC1 | 30.22 ± 1.39 | 7.12 ± 0.48 | 1.21 ± 0.05 | 2.08 ± 0.24 | 34.47 ± 1.25 | 29.98 ± 1.47 | 40.63 ± 1.72 | 61.47 ± 1.78 |
NDC-6 | 4.61 ± 0.55 | 4.13 ± 0.38 | 0.48 ± 0.17 | 1.64 ± 0.38 | 52.15 ± 1.69 | 42.37 ± 2.04 | 10.88 ± 0.57 | 89.68 ± 1.35 |
Functional genes related to NR
The genes hao, napA, nirS, nirK, cnorB, and nosZ are generally considered as the functional genes encoding key enzymes involved in microbial nitrogen metabolism. In this study, PCR amplification of these six functional genes was conducted across all bacterial strains to further validate their nitrification and denitrification capabilities, thereby elucidating the potential metabolic pathways of NR under low-temperature conditions. The amplification results are presented in Fig. 5 and summarized in Table S5.
The hao gene, encoding hydroxylamine oxidase (HAO), is a critical biomarker for nitrification. HAO catalyzes the oxidation of hydroxylamine (NH2OH), a key intermediate in nitrification, to NO2−-N [60]. In this study, the successful amplification and expression of the hao gene in P. veronii HN1 (282 bp), P. poae HN2 (265 bp), and P. peli HN3 (1752 bp) confirmed the involvement of hao in nitrification in these three strains. The hao gene was likewise amplified in other nitrification-capable Pseudomonas strains, such as P. mendocina SCZ-2 [61] and P. citronellolis YN-21 [27]. Furthermore, PCR amplification revealed the absence of functional genes encoding nitrite reductases (NIR) in these three strains. As NIR catalyzes the reduction of NO2−-N to gaseous NO or N2O [62], this genetic characteristic explained distinct accumulation of NO2−-N observed during the NR process. Notably, the nosZ gene, encoding nitrous oxide reductase (NOS), was successfully amplified from strains HN1 (396 bp), HN2 (352 bp), and HN3 (298 bp). Given the well-established catalytic role of NOS in reducing N2O to N2, as demonstrated by Guo et al. [63], this finding signified that these three strains can convert N2O generated from ancillary denitrification processes into N2, thereby mitigating the emission of greenhouse gases.
As illustrated in Fig. 5 and Table S5, typical denitrification genes (napA, nirS, nirK, cnorB, and nosZ) were PCR-amplified from Aeromonas sp. AD1, P. extremaustralis AD2, and S. liquefaciens AD3. The napA gene, a biomarker involved in aerobic denitrification, encodes the periplasmic nitrate reductase (NAP), which catalyzes the reduction of NO3−-N to NO2−-N as demonstrated by Zheng et al. [64]. The napA gene has also been amplified from other psychrotolerant denitrifiers, such as Priestia aryabhattai KX-3 [65] and Psychrobacter cryohalolentis strain F5-6 [66]. Concurrently, nirS and nirK were identified as functional genes encoding NIR. These two types of functional genes in this study could simultaneously exist in the same strain and promoted the removal of NO2−-N. This genetic configuration correlated with the observed low NO2−-N accumulation during the removal of NO3−-N by strains AD1, AD2, and AD3, and their composite consortium, indicating synergistic functionality of multiple NIR isoforms in nitrite metabolism. The cnorB gene encoding nitric oxide reductase (NOR) for the conversion of NO to N2O and the nosZ gene responsible for the reduction of N2O to N2 were detected in Aeromonas sp. AD1, P. extremaustralis AD2, and S. liquefaciens AD3, respectively, indicating that these three strains prevented NO and N2O emission by directly reducing them to N2 during denitrification. These findings aligned with the previously observed low concentrations of N2O produced by these three strains during denitrification. In summary, successful amplification of napA, nirS, nirK, cnorB, and nosZ in strains AD1, AD2, and AD3 indicated that NAP, NIR, NOR, and NOS are involved in different steps of the dissimilatory nitrate reduction pathway. This genetic evidence confirms their capability to achieve complete denitrification without intermediate accumulation. This denitrification pathway was consistent with the nitrogen-removal mechanism in P. aeruginosa WS-03 reported by Wei et al. [28].
Pathways for nitrogen metabolism
Based on the results of nitrogen balance analysis and PCR amplification of functional genes related to NR, the proposed pathways for nitrogen metabolism in NDC-6 are depicted in Fig. 6. The composite consortium NDC-6, constructed from heterotrophic nitrifiers and aerobic denitrifiers, was speculated to possess a complete pathway for nitrogen metabolism and to achieve synchronous and efficient removal of NH4+-N and NO3−-N through dual metabolic pathways: assimilatory conversion into intracell-N and dissimilatory transformation (heterotrophic nitrification and aerobic denitrification) to N2. The initial step of heterotrophic nitrification in the composite consortium NDC-6 involved the oxidation of NH4+-N to NH2OH using sodium succinate as the carbon source and electron donor, mediated by constituent strains P. veronii HN1, P. poae HN2, and P. peli HN3. Subsequently, HAO from these nitrifying strains catalyzed further oxidation of NH2OH into NO2−-N. Concurrently, during aerobic denitrification, NAP expressed by Aeromonas sp. AD1, P. extremaustralis AD2, and S. liquefaciens AD3 reduced NO3−-N to NO2−-N. In the following steps, the NO2−-N generated through nitrification and initial denitrification steps underwent sequential enzymatic conversions: NIR catalyzed the reduction of NO2−-N to NO, followed by NOR-mediated transformation of NO to N2O, and final reduction to N2 via NOS. The functional enzymes involved in these steps were collectively expressed by Aeromonas sp. AD1, P. extremaustralis AD2, and S. liquefaciens AD3.
Overall, the composite consortium NDC-6 achieved efficient removal of NH4+-N, NO3−-N, and NO2−-N through synergistic metabolism and complementary functions among its constituent strains, which facilitated a broad ecological niche for this CMA and enhanced its environmental adaptability, resulting in superior metabolic capabilities for diverse nitrogen substrates in wastewater treatment.