Diversity of microbial community structure
By using the Megahit (version 1.12) assembly software, a total of 14254368 contigs were obtained, and the total sequence length predicted by metaGene was 9282883724bp. This study was conducted at the level of phylum and genus according to the experimental conditions.
At phylum level, the abundance of Proteobacteria was more than 59.87%, Bacteroidetes was more than 5.77%, Actinobacteria was more than 4.88%, Acidobacteria was more than 3.39%, and planktomycetes was more than 3.25%, see Fig. 2(a). With the dry-wet alternation time (DAWT 4h/8h/12h, see Fig. 2(b)), Proteobacteria showed a "V" pattern, while other bacteria, include Nitrospirae showed an inverted "V" pattern.
At genus level, the abundance of Mesorhizobium, Hyphomicrobium, Bradyrhizobium, Hydrogenophaga and Pseudoxanthomonas were more than 3.08%, 2.70%, 1.92%, 1.70% and 1.15%, respectively, see Fig. 3(a). With the dry-wet alternation time (DAWT 4h/8h/12h,Fig. 3(b)), the main microbial species showed a decreasing trend at genus level: Mesorhizobium, Hyphomicrobium, Hydrogenophaga, Bradyrhizobium and Genophaga showed an inverted "V" decreasing trend, the number of denitrifying bacteria such as Sphingomonas and Sphingobium increased in a "V" pattern, while the number of ammonia-oxidizing bacteria such as Nitrotomonas and Nitrosospira increased in a "V" pattern. It can be seen that the abundance of microorganisms varied greatly at DAWT 8h, which indicated that the microbes were adjusting the structure of their own genus to adapt to the change of dry-wet alternate environment.
Community distribution characteristics of Denitrogen bacteria
The genus level distribution of denitrogen bacteria is shown in Table 2 and Fig. 4, the total abundance of denitrifying bacteria is more than 10.0%. Among them, Nitrosomonas, Nitrosospira and Nitrosococcus are the main ammonia-oxidizing bacteria, the total abundance is 0.18%, the proportion of denitrogen bacteria is about 1.8%. Nitrospira, Nitrospira and Nitrococcus are the main nitrifying bacteria, the total abundance is 0.72%, the proportion of denitrogen bacteria is about 7.1%. The main denitrifying bacteria, such as Pseudomonas, Pseudomonas, chloomonas, Flavobacterium, Rubrivivax, Bacillus, Pedobacter, Thauera, Hyphomicrobium, Hydrogenophaga, Rhizobium, Azoira and Acidovorax, the total abundance was 9.04%, accounting for about 90.4% of the total denitrogen bacteria. The abundance of Candidatus, Candidatus, Candidatus and Candidatus was 0.07%, which accounted for about 0.7% of the total Anammox. Among all the denitrogen bacteria, denitrifying bacteria were the most abundant, nitrifying bacteria the second, and Anammox the least.
Table 2
The relative abundance of nitrifying bacteria in Genus cell(%)
Taxonomy | H-121814 | H-121815 | H-121816 | H120521 | H120522 | H120523 | H-112331 | H-112335 | H-112336 |
Nitrosomonas | 0.003571 | 0.004048 | 0.006582 | 0.004355 | 0.004333 | 0.006356 | 0.002755 | 0.007518 | 0.007767 |
Nitrosococcus | 0.004818 | 0.003665 | 0.006676 | 0.004224 | 0.004126 | 0.005666 | 0.002613 | 0.006286 | 0.006902 |
Nitrosospira | 0.004216 | 0.004019 | 0.005335 | 0.003475 | 0.004108 | 0.00517 | 0.001665 | 0.005384 | 0.005637 |
Nitrospira | 0.036161 | 0.154328 | 0.07075 | 0.033127 | 0.070847 | 0.060423 | 0.042152 | 0.060783 | 0.025703 |
Nitrospina | 0.001368 | 0.001577 | 0.002308 | 0.00126 | 0.001281 | 0.001649 | 0.000671 | 0.002243 | 0.00215 |
Nitrococcus | 0.00154 | 0.001276 | 0.002043 | 0.001223 | 0.001207 | 0.00156 | 0.000495 | 0.001904 | 0.00178 |
Bacillus | 0.008502 | 0.01072 | 0.012414 | 0.006831 | 0.008149 | 0.009092 | 0.003667 | 0.013039 | 0.011989 |
Terrimonas | 0.031969 | 0.04859 | 0.04786 | 0.029913 | 0.050711 | 0.030891 | 0.020916 | 0.11891 | 0.068656 |
Flavobacterium | 0.030308 | 0.022656 | 0.048869 | 0.02005 | 0.023516 | 0.017929 | 0.011496 | 0.040032 | 0.036726 |
Pedobacter | 0.01716 | 0.015249 | 0.017477 | 0.010259 | 0.014051 | 0.009979 | 0.007036 | 0.026693 | 0.024509 |
Pseudomonas | 0.066751 | 0.042267 | 0.065602 | 0.038607 | 0.04459 | 0.052721 | 0.021127 | 0.063141 | 0.067214 |
Rubrivivax | 0.021594 | 0.026745 | 0.029076 | 0.017428 | 0.014441 | 0.017231 | 0.010591 | 0.026823 | 0.027149 |
Dechloromonas | 0.009804 | 0.007698 | 0.014551 | 0.013937 | 0.019105 | 0.019209 | 0.007687 | 0.013236 | 0.024752 |
Thauera | 0.013404 | 0.01308 | 0.016417 | 0.011293 | 0.012353 | 0.017332 | 0.007325 | 0.014217 | 0.015694 |
Hyphomicrobium | 0.255264 | 0.221496 | 0.235097 | 0.262102 | 0.23177 | 0.238969 | 0.085411 | 0.230893 | 0.243098 |
Hydrogenophaga | 0.194108 | 0.067588 | 0.164097 | 0.113839 | 0.136871 | 0.10198 | 0.116131 | 0.213241 | 0.146844 |
Rhizobium | 0.085387 | 0.072496 | 0.05181 | 0.047517 | 0.035739 | 0.044031 | 0.01544 | 0.059142 | 0.086671 |
Azospira | 0.005949 | 0.005127 | 0.008731 | 0.005874 | 0.007195 | 0.013896 | 0.003384 | 0.006305 | 0.01382 |
Acidovorax | 0.214717 | 0.102938 | 0.131994 | 0.071753 | 0.069273 | 0.075746 | 0.109424 | 0.182268 | 0.120183 |
Candidatus_Brocadia | 0.002418 | 0.002779 | 0.003385 | 0.001546 | 0.00213 | 0.002533 | 0.000781 | 0.003412 | 0.003058 |
Candidatus_Jettenia | 0.0009 | 0.001128 | 0.00145 | 0.000578 | 0.000938 | 0.001272 | 0.00044 | 0.001554 | 0.001138 |
Candidatus_Kuenenia | 0.000877 | 0.001065 | 0.001572 | 0.000752 | 0.000917 | 0.00115 | 0.000412 | 0.001434 | 0.001158 |
Candidatus_Scalindua | 0.000701 | 0.000911 | 0.001137 | 0.000539 | 0.000665 | 0.000873 | 0.000303 | 0.001156 | 0.000877 |
The composition of denitrogen bacteria in constructed wetland reactor was different under different dry-wet alternate conditions:
(1) At 10℃, compared the abundance of ammonia-oxidizing bacteria, DWAT12h (H-121816) > DWAT4h (H-121814) > DWAT8h (H-121815); Compared with nitrifying bacteria, DWAT8h(H-121815) > DWAT4h(H-121814) > DWAT12h(H-121816); Compared with denitrifying bacteria, DWAT4h(H-121814) > DWAT12h (H-121816) > DWAT8h(H-121815); Compared with Anammox, DWAT12h (H-121816) > DWAT8h (H-121815) > DWAT4h (H-121814).
(2) At 20℃, compared the abundance of ammonia-oxidizing bacteria, DWAT12h (H-121816) > DWAT8h (H-121815) > DWAT4h (H-121814); Compared with nitrifying bacteria, DWAT4h(H-121814) > DWAT12h(H-121816) > DWAT8h(H-121815); Compared with denitrifying bacteria, DWAT8h(H-121815) > DWAT4h(H-121814) > DWAT12h(H-121816); Compared with Anammox, DWAT12h (H-121816) > DWAT8h (H-121815) > DWAT4h (H-121814).
(3) At 30℃, compared the abundance of ammonia-oxidizing bacteria, DWAT12h (H-121816) > DWAT8h (H-121815) > DWAT4h (H-121814);; Compared with nitrifying bacteria, DWAT8h (H-121815) > DWAT4h(H-121814) > DWAT12h(H-121816); Compared with denitrifying bacteria, DWAT8h(H-121815) > DWAT12h (H-121816) > DWAT4h (H-121814); Compared with Anammox, DWAT12h (H-121816) > DWAT8h (H-121815) > DWAT4h (H-121814).
The results showed that the abundance of DWAT4h and DWAT8h of ammonia-oxidizing bacteria was higher at 12h than at 8h of nitrifying bacteria, and higher at 20 °c than at 4h of DWAT12h of nitrifying bacteria Compared with DWAT4h and DWAT12h, the abundance of DWAT4h and DWAT12h of denitrifying bacteria was highest at 10°C, and that of Anammox was highest at 12h.
Response Analysis of denitrogen bacteria and dry-wet alternation
(1) Response Analysis at 10℃
Table 3 shows that the main ammonia-oxidizing bacteria (Nitrosomonas, Nitrospira and Nitrosococcus) and Nitrospira (Nitrospira, Nitrospina and Nitrococcus) are significantly correlated with NO2-N, DO and DWAT at 10℃(P < 0.05), the results showed that the concentration of NO2-N and DWAT affected ammonia oxidation and nitrification to some extent. The main denitrifying bacteria such as Terrimonas, Pseudomonas, Dechloromonas and Flavobacterium were significantly correlated with NO3-N, ON and TOC (P < 0.05), which indicated that the concentration of NO3-N affected the growth of Denitrification cells to some extent. The main Anammox bacteria such as Candidatus_Brocadia, Candidatus_Jettenia, Candidatus_kuenenia and Candidatus_Scalindua were significantly correlated with NH3-N, NO2-N, DO and DWAT (P < 0.05). The results showed that NO2-N, NO3-N, DO and DWAT were important factors affecting the distribution of denitrogen bacteria in the substrate at 10°C.
Table 3
Analysis of correlation between denitrogen bacteria and environmental factors (10℃)
Taxonomy | TN | NH3-N | NO3-N | NO2-N | ON | DO | TOC | DWAT |
Nitrosomonas/Nitrosospira/Nitrosococcus | -0.346 | -0.512 | -0.03 | -0.912* | -0.173 | -0.870* | 0.430 | 0.612* |
Nitrospira/Nitrospinav/ Nitrococcus | 0.419 | 0.578 | -0.049 | 0.877* | 0.215 | 0.828* | -0.357 | -0.587 |
Terrimonas/Pseudomonas/Dechloromonas/Flavobacterium /… | -0.698 | -0.555 | 0.915* | -0.763 | -0.835* | 0.820* | -1.000* | 0.555 |
Candidatus_Brocadia/Candidatus_Jettenia/Candidatus_kuenenia/Candidatus_Scalindua | 0.527 | 0.674* | -0.171 | -0.811* | 0.333 | 0.753* | -0.240 | -0.674* |
Note: N = 3;*, P < 0.05༛**, P < 0.01༛ |
(2) Response Analysis at 20℃
Table 4 shows that the main ammonia-oxidizing bacteria (Nitrosomonas, Nitrospira and Nitrosococcus) is significantly correlated with TN, NO3-N, DO and DWAT at 20℃(P < 0.05). The Nitrospira (Nitrospira, Nitrospina and Nitrococcus) is significantly correlated with NO3-N, NO2-N, DO and DWAT (P < 0.05). The main denitrifying bacteria such as Terrimonas, Pseudomonas, Dechloromonas and Flavobacterium were significantly correlated with NO3-N, NO2-N, DO and DWAT (P < 0.05). The main Anammox bacteria such as Candidatus_Brocadia, Candidatus_Jettenia, Candidatus_kuenenia and Candidatus_Scalindua were significantly correlated with NH3-N, NO2-N and TOC (P < 0.05). The results showed that NO2-N, NO3-N, DO and DWAT were important factors affecting the distribution of denitrogen bacteria in the substrate at 20°C.
Table 4
Correlation between denitrogen bacteria and environmental factors (20℃)
Taxonomy | TN | NH3-N | NO3-N | NO2-N | ON | DO | TOC | DWAT |
Nitrosomonas/Nitrosospira/Nitrosococcus | -0.772* | -0.460 | -0.979* | 0.608 | -0.410 | -1.000** | 0.5385 | 0.824* |
Nitrospira/Nitrospinav/ Nitrococcus | -0.608 | -0.248 | -0.908* | 0.999* | -0.193 | -0.974* | -0.333 | -0.675* |
Terrimonas/Pseudomonas/Dechloromonas/Flavobacterium /… | 0.333 | -0.064 | 0.734* | -0.931* | -0.120 | 0.857* | 0.025 | 0.414 |
Candidatus_Brocadia/Candidatus_Jettenia/Candidatus_kuenenia/Candidatus_Scalindua | 0.524 | 0.816* | 0.067 | 0.870* | 0.847* | -0.137 | -0.992* | 0.447 |
Note: N = 3;*, P < 0.05༛**, P < 0.01༛ |
(3) Response Analysis at 30℃
Table 5 shows that the main ammonia-oxidizing bacteria (Nitrosomonas, Nitrospira and Nitrosococcus) is significantly correlated with TN, NO3-N, NO2-N, ON, DO and DWAT at 30℃(P < 0.05). The Nitrospira (Nitrospira, Nitrospina and Nitrococcus) is significantly correlated with NO3-N, ON and TOC (P < 0.05). The main denitrifying bacteria such as Terrimonas, Pseudomonas, Dechloromonas and Flavobacterium were significantly correlated with NO3-N and TOC (P < 0.05). The main Anammox bacteria such as Candidatus_Brocadia, Candidatus_Jettenia, Candidatus_kuenenia and Candidatus_Scalindua were significantly correlated with TN, NO3-N, NO2-N ON, DO and TOC (P < 0.05). The results showed that NO2-N, NO3-N, DO and DWAT were important factors affecting the distribution of denitrogen bacteria in the substrate at 30°C.
Table 5
Correlation between denitrogen bacteria and environmental factors (30℃)
Taxonomy | TN | NH3-N | NO3-N | NO2-N | ON | DO | TOC | DWAT |
Nitrosomonas/Nitrosospira/Nitrosococcus | 0.953* | -0.482 | 0.871* | -0.984* | 0.999* | 0.700* | -0.362 | -0.827* |
Nitrospira/Nitrospinav/ Nitrococcus | 0.293 | 0.447 | -0.899* | -0.410 | 0.609* | -0.191 | -0.973* | -0.006 |
Terrimonas/Pseudomonas/Dechloromonas/Flavobacterium /… | 0.026 | -0.709 | -0.713* | 0.098 | -0.325 | 0.494 | 0.996* | -0.312 |
Candidatus_Brocadia/Candidatus_Jettenia/Candidatus_kuenenia/Candidatus_Scalindua | 0.823* | -0.204 | 0.977* | -0.887* | 0.970* | 0.459 | -0.620 | -0.625* |
Note: N = 3;*, P < 0.05༛**, P < 0.01༛ |
Characteristics of denitrogen bacteria and environmental factors
The results of Redundancy Analysis (RDA) between denitrogen bacteria and environmental factors are shown in Fig. 5. The main environmental factors, NH3-N (VIF value 2.32597), NO3-N (VIF value 1.90636), TN (VIF value 1.931025), were screened by analysis of Variance inflation factor (VIF). Then, the main environmental variables were analyzed, and the results showed that NH3-N (AN) and TN (TN) had a significant effect (Table 6, PAN0.019 < 0.05, PTN0.001 < 0.05), the correlation coefficients between denitrogen bacteria and environmental factors were 0.63 and 0.84, respectively. At the same time, different DAWT4h, 8h and 12h groups were distributed in four quadrants, and the distances between the points were larger, which indicated that the functional composition of the samples were different. The analysis shows that NH3-N (AN) and TN are the main factors that influence the difference of denitrogen microbial community, that is, the difference of environmental factors caused by dry-wet alternation affects the structure and distribution of denitrogen microbial community [8, 21–22].
Table 6
Redundancy analysis (RDA) of denitrogen bacteria
| RDA1 | RDA2 | r2 | P_values |
A_N | 0.996577 | 0.082664 | 0.632249 | 0.019 |
X_N | 0.956327 | -0.2923 | 0.476429 | 0.061 |
T_N | -0.56045 | -0.82819 | 0.837563 | 0.001 |