2.1 Variation of water quality parameters
As shown in Table 1, the temperature range of the stagnant water was 11.3 - 17.9℃ and the fresh water was 8.2 - 11.4℃. The temperature of water samples increased significantly after stagnation (P < 0.001). The pH value of stagnant water was not significantly different (P > 0.05) from that of fresh water. After stagnation, the pH value was also within the normal range specified in the official Chinese standards for drinking water quality (GB 5749-2006). The results obtained in this study was similar to the observation made by Zhang et al. [19], in which the temperature of the stagnant overnight water conducted range from 15 to 17℃, and the temperature of the fresh water dropped to 4 - 6℃ after 10-minute flush in winter season.
Table 1
Water quality parameters of stagnant and fresh water in winter season.
Water quality
|
Stagnation
|
Fresh
|
t-test
|
T (℃)
|
13.84 ± 1.50
|
10.26 ± 0.76
|
***
|
pH
|
7.79 ± 0.24
|
7.75 ± 0.19
|
NS
|
Free residual chlorine (mg/L)
|
ND
|
0.18 ± 0.05
|
NS
|
Total residual chlorine (mg/L)
|
0.10 ± 0.01
|
0.26 ± 0.02
|
***
|
NO3−-N (mg/L)
|
1.71 ± 0.16
|
1.86 ± 0.16
|
**
|
NH4+-N (mg/L)
|
0.06 ± 0.02
|
0.07 ± 0.02
|
NS
|
TP (mg/L)
|
0.02 ± 0.01
|
0.02 ± 0.01
|
NS
|
TOC (mg/L)
|
1.45 ± 0.02
|
1.51 ± 0.03
|
*
|
Fe (mg/L)
|
0.11 ± 0.05
|
0.02 ± 0.01
|
**
|
Values showed as means and standard deviations ( n = 15). * P < 0.05, ** P < 0.01, and *** P < 0.001 represent statistical significance using t-test. NS represents no significance. ND represents no detected. |
Residual chlorine is one of the most common indicators that affects water quality and bacterial community structure in drinking water. As shown in Table 1, the concentration of total and free residual chlorine decreased significantly after stagnation (P < 0.001). The total residual chlorine concentrations decreased by 48.22 - 65.07% after stagnation and the free residual chlorine concentrations were below the detection limit in the stagnant water samples. These results were consistent with a previous study [40], which indicated that lower pH and higher temperature accelerated the chlorine decay rate in copper pipes. Under the condition of stagnation, the residual chlorine decayed over time due to the consumption of microorganisms, organic matter and invasive material [41].
The t-tests showed that there were no significant differences in the NH4+-N and TP concentrations before and after flushing (P > 0.05) (Table 1). In contrast, the concentrations of NO3−-N, Fe and TOC differed significantly (P < 0.05). Briefly, the average contents of NO3−-N was reduced by 6.1% after stagnation. This result was contrary to the findings presented by Masters et al. [42], which found the concentration of NO3−-N increased when the stagnant time was less than 1 day in their study. Masters et al. [42] also found that facultative denitrifying bacteria (measured by the abundance of denitrifying gene) were ubiquitous in drinking water pipes and the gene abundance increased with water age. Autotrophic denitrifying bacteria utilized H2 produced by iron as electron donor and nitrate as electron acceptor to convert NO3−-N to N2 under a water environment with metal iron, which might be responsible for the reduction of NO3−-N content. At the same time, this reaction might cause the rise in pH value and induce iron corrosion.
As shown in Table 1, the average concentration of Fe in the stagnant water samples (0.11 mg/L) was 4.43-fold higher than that of fresh water samples (0.02 mg/L). Iron was released from the pipe into the water during stagnation. This conclusion was similar to the results observed by Tian et al. [4], in which the stagnant time was found to be an essential influencing factor for the release of Fe. However, only long periods of stagnation (>8 h) may cause considerable risk and risk associated with overnight stagnation is considerably low. According to previous researches, water quality factors such as chloride, dissolved organic matter and pH would affect the stability of Fe in the pipeline [43–45]. Their studies disclosed that chloride was deemed to increase Fe concentration in DWDS. On the contrary, the increasing of pH was proved to enhance the Fe stability and reduce the Fe release in actual distribution systems.
As an indispensable element for bacterial regrowth in DWDS, trace amount of organic carbon (1 µg/L) could support up to 103 - 104 cell/mL microbes to reproduce [17, 46]. The average TOC concentration of fresh water (1.51 mg/L) was generally higher than that of stagnation (1.45 mg/L) in this study. A previous study also observed downtrends in TOC content with increased stagnation time, and the 30% reduction was examined after 168-hour of stagnation in the TOC concentrations [17]. Moreover, organic matter was known to serve as precursors for disinfection byproducts, and the presence of organic matters supplement increased the level of free chlorine residual required to control the microorganisms [47].
2.2 Changes in biomass and bio-activity
As “the energy currency” of biological cells, ATP was used to evaluate the bacterial activity. Fig. 1a, b displayed a higher concentration of ATP before flushing (P < 0.05). The concentration of total ATP and intracellular ATP increased by 1.65 and 1.69 times after stagnation, respectively. Series of studies also found that ATP increased after overnight stagnation [31, 34], which probably due to the suitability of indoor temperature for bacterial regeneration in the pipelines. In addition, the enhancement in metabolic activity and reproductive capacity of the cells in the stagnant samples were expressed by increasing intracellular ATP which is consistent with conclusions made in a recent study [8]. However, the intracellular ATP concentration per cell in stagnant water was lower than that measured by Zhang et al. [8], which might be due to the affection of temperature and the proportion of LNA cells to total cell count. Wang et al. [30] compared the growth rate of LNA bacteria at 12, 20, 25 and 30℃, and the results showed that the growth rate of LNA bacteria was the highest at 30℃. Meanwhile, ATP concentration of LNA bacteria increased continuously during the stagnation process.
As shown in Fig. 1c and Supplementary Figure 1, ICC was significantly higher after stagnation (P < 0.001). The average intact cell count in stagnant water was 6.99 × 104 cells/mL (range = 3.96 - 12.31 × 104 cells/mL), which was 1.53-fold as much as that in fresh water (2.44 × 104 cells/mL). These results were consistent with those conducted in previous studies [17, 48, 49]. The increasing ICC might be due to the variations of in the water environment in taps during indoor heating, including higher temperature and lower chlorine content [48]. Moreover, the free chlorine concentration was also an important factor limiting the number of bacteria, and the low chlorine concentration at the end of the faucet created favorable conditions for bacterial regeneration [17]. Furthermore, the release of biofilms was also an important source of bacteria. Ling et al. [49] found that the communities of water sampled in taps were more similar to that in biofilms than that in pipeline water. At the same time, the scour effect caused by the change of hydraulic conditions during the sampling process triggered the transportation of bacteria from the biofilm and pipe sediments to the water [1]. In general, the combination of the water environment, biofilm release, and scour effect caused the higher ICC in stagnation samples.
LNA bacteria could adopt dormancy strategy to deal with adverse environmental conditions [28]. LNA bacteria have higher survival competitiveness than HNA bacteria in oligotrophic environment, and become the dominant bacteria in the environment. Supplementary Figure 2 indicated the FCM graph of LNA and HNA cells. The percentage of LNA cells in the Intact cell counts (LNA%) is shown in Fig. 1d. The number of LNA cells in stagnant water samples (range = 1.53 - 4.38 × 104 cells/mL) increased significantly compared with that in fresh water samples (range = 1.24 - 2.11 × 104 cells/mL). However, LNA% in stagnant water (range = 24.29% - 48.48%) samples presented a downward trend compared with that in fresh water samples (range = 41.08% - 67.88%). The increase in LNA and HNA bacteria during stagnation was agreed with the results observed by Zlatanović et al. [17]. Gabrielli et al. [34] also found that LNA% generally decreased during night (due to a higher utilization of nutrients was conducive to the growth of HNA bacteria and/or the passage from LNA to HNA).
In addition to the relationship between HNA/LNA bacteria and ATP concentration, Wang et al. [25] indicated that increased counts of HNA cells might be associated with a 4-fold increase in ATP concentration. The transfer of LNA cells to HNA cells could account for a slight increase of ATP concentrations of stagnation. The conclusion was consistent with a previous study by Liu et al. [26], which reported that ATP-per-cell content of HNA cells were 10-fold higher than the LNA cells.
2.3 Bacterial community structure
Through the sequencing, 205,760 reads were detected from the stagnant water samples (11,165 to 45,878 reads per sample) and the fresh water samples (10,924 to 40,003 reads per sample). As shown in Supplementary Table 1, the mean length of clean reads was 1374 bp, and the length of clean reads were obtained from the stagnant water ranged from 1313 to 1400.71 bp per sample and in fresh water samples from1302.64 to 1401.16 bp per sample). It is worth noting that the sequence number mainly focused on sequencing ranged from 1401 to 1500 bp.
In total, 34,960 and 34,955 OTUs were produced in the fresh and stagnant water samples, respectively (Table 2). The average OTUs number in stagnant water was 227 (range = 183 - 280), which was slightly higher than OTUs observed for freshwater samples (218; range =199 - 236). A previous reported found that the average OTUs number was 449 per sample in stagnant water, and 323 per sample in fresh water (8). The low number of OTUs observed in the present study, possibly due to the full-length sequencing with clearer bacterial information and made the clustering more explicit.
Table 2
Bacterial community richness and diversity indices in each stagnant and fresh water sample.
Simple ID
|
0.97 level
|
OTUs
|
Chao 1
|
Coverage
|
Shannon
|
Simpson
|
HG-S
|
280
|
319 (298, 364)
|
0.999
|
3.74 (3.72, 3.76)
|
0.0582 (0.0568, 0.0596)
|
HG-F
|
214
|
240 (226, 274)
|
0.997
|
3.70 (3.67, 3.73)
|
0.0459 (0.0472, 0.0506)
|
YF-S
|
228
|
265 (244, 313)
|
0.996
|
3.83 (3.80, 3.86)
|
0.0435 (0.0422, 0.0448)
|
YF-F
|
233
|
264 (248, 298)
|
0.994
|
3.64 (3.60, 3.67)
|
0.0621 (0.0597, 0.0645)
|
SS-S
|
183
|
294 (273, 336)
|
0.998
|
3.04 (3.05, 3.06)
|
0.1330 (0.1304, 0.1357)
|
SS-F
|
199
|
321 (287, 364)
|
0.998
|
2.90 (2.88, 2.92)
|
0.1511 (0.1484, 0.1538)
|
YGC-S
|
218
|
260 (239, 303)
|
0.993
|
3.51 (3.47, 3.55)
|
0.0716 (0.0687, 0.0745)
|
YGC-F
|
210
|
258 (236, 304)
|
0.995
|
3.06 (3.02, 3.10)
|
0.1353 (0.1305, 0.1401)
|
BG-S
|
225
|
262 (246, 296)
|
0.994
|
3.87 (3.84, 3.90)
|
0.0427 (0.0410, 0.0444)
|
BG-F
|
236
|
293 (265, 350)
|
0.993
|
3.81 (3.77, 3.84)
|
0.0442 (0.0427, 0.0506)
|
“S” represents the sample of stagnant samples, and “F” represents the sample of fresh samples. |
Table 2 exhibits Chao 1, Shannon and Simpson indices obtained for the stagnant and fresh water samples. The OTUs values of the bacterial communities were all lower than the Chao 1 richness index, indicating that there were many unknown bacterial sequences in DWDS (50). Moreover, greater richness and diversity was discovered after stagnation due to higher alpha diversity [18, 51].
As shown in Fig. 2a, the dominant phylum in all water samples was Proteobacteria which accounted for 87.21 ± 6.54% of the classified sequences, followed by Actinobacteria (8.25 ± 4.55%). Higher abundance of Proteobacteria in the stagnant water samples (87.94 ± 5.80%) was observed as compared to the fresh water samples (86.47 ± 6.16%). Approximate results were obtained in DWDS by Zhang et al. [8] and Ling et al. [49], which presented that the OTUs related to Proteobacteria, Actinobacteria, and Bacteroidetes became more abundant after stagnation.
The dominant bacterial genera observed in the stagnant and fresh water samples are shown in Fig. 2b. Obvious changes in the bacterial community structure were observed after stagnation. In fresh water samples, the number of OTUs assigned to Curvibacter sp., Mesorhizobium sp. and Methylovirgula sp. were higher i.e., 25.42 ± 13.63% vs 36.68 ± 26.38%, 2.38 ± 1.33% vs 4.19 ± 2.37% and 1.84 ± 0.94% vs 2.89 ± 0.81%, respectively. Ralstonia sp., Acinetobacter sp. and Phreatobacter sp. became more abundant after stagnation, i.e., 17.09 ± 13.69% vs 11.55 ± 8.34%, 4.10 ± 3.76% vs 0.03 ± 0.00% and 3.55 ± 5.03% vs 2.43 ± 1.76%, respectively (Supplementary Figure 3). The OTUs number with increased abundance at the species level is shown in Supplementary Figure 4. Furuhata et al. [52] found that Ralstonia sp. was the most dominant bacteria in cold water samples from water dispensers, which also presented in healthy human bodies and widely distributed in the upper respiratory tract, such as the oral cavity, pharynx and trachea. Meanwhile, earlier research has referred that it was very rare for a healthy person to develop an infection caused by these bacteria, but for those people with an immune deficiency or reduced resistance might develop respiratory infections or blood poisoning due to Ralstonia sp. [53]. In addition, Acinetobacter sp. was a common species in tap water and it became the dominant bacteria when the stagnant time reached 12 h [54]. With its increasing antibiotic tolerance and its potential risk leading to hospital infections, natural disaster-related infections and community-acquired infections, Acinetobacter sp. had become an issue of major public health [55]. Furthermore, Perrin et al. [56] found that the relative abundance of Phreatobacter sp. was significantly higher during warmer period. A relevant report indicated that the abundance of Phreatobacter sp. was negatively associated with the content of chlorine in water bodies [57]. The abundance of Phreatobacter sp. in stagnant water might be caused by biofilm release [58].
As shown in Supplementary Figure 3, the abundance of Pseudomonas sp. and Mycobacterium sp. increased after overnight stagnation, i.e., 120.08% and 585.71% higher, respectively. This result was consistent with the previous reports [5, 8, 59] in which they indicated that the abilities of Mycobacteria sp. and other OPPPs to persist and survive in the pipeline environment were attributed to the ability of biofilms formation, the relative resistance to disinfectant, and the strong adaptability to oligotrophic conditions. More importantly, the pathogen communities could cause waterborne diseases such as pneumonia, sepsis and infection [55, 60]. Thus, more attention should be paid to these microbes.
AIC index was used to explain the fitting of the abundance distribution models of several theoretical taxa to stagnant and fresh water data [37]. As shown in Fig. 3, the lognormal model was both fit for stagnant water (AIC = 6213.95) and fresh water (AIC = 5794.42). A recent study confirmed the conclusion, which indicated that the global distribution of abundance in microbial communities was lognormal [61]. The lognormal model represents a distribution in which both extremely abundant and extremely rare species have small proportions [62]. The community dynamics presented by the lognormal model is a random zero-sum process, which indicate that the death or migration of one individual in a community is immediately followed by the emergence of another random individual to fill the vacancy [63]. Since only a few species dominate, the species abundance distribution of drinking water might be determined by ecological niche.
As a neo-analytical method, network analysis is generally used to study the interactions of microbe-microbe and microbe-environment interactions in complicated ecological systems [10, 64]. The different networks containing bacterial communities at genus were analyzed and presented in Fig. 4 (|r| > 0.6, P < 0.05). The topological properties of networks are presented in Supplementary Table 2. As shown in Fig. 4, 108 edges were calculated (76.11% positive and 23.89% negative) for bacteria in stagnant water samples, and 112 edges (73.15% positive and 26.85% negative) for bacteria in fresh water samples. This conclusion revealed that the interaction with respect to symbiosis, coaggregation, and cocolonization among microorganisms was more abundant in the drinking water pipeline ecosystem [38]. Compared with fresh water, the average path length of stagnant water was lower, which illustrated microorganisms in stagnant water possessed the stronger ability of material exchange and energy transfer [65]. In addition, compared with fresh water (6 modules), the number of modules increased after stagnation (8 modules), which indicated that the survival modes of bacterial communities were more complicate after stagnation [10].
Node microorganisms with a large number of connections were critical in influencing the structure and function of the microbial community in water system. The bacterial interactions varied obviously after stagnation. The nodes with the greatest connections in fresh water were Phreatobacter sp., Variovorax sp. and Tsukamurella sp., but the key nods were Chitinophaga sp. and Mesorhizobium sp. in stagnant water. Previous studies informed that Chitinophaga sp. and Mesorhizobium sp. were commonly isolated from soil and plants which were proved to affect the survival of other microorganisms. The chitinolytic activity of Chitinophaga sp. might antagonize a range of aquatic microorganisms [66]. Krick et al. [67] investigated a marine Mesorhizobium sp. and reported the presence of several N-acyl homoserine lactones (NAHLs). The bacteria reported to use AHLs for signal transmission between genera with certain antibacterial and cytotoxic activities. Chitinophaga sp. and Mesorhizobium sp. were not most abundant but played a crucial role in interspecific regulation of bacteria. However, the interspecies influence of bacteria in the stagnant water remains to be further studied, although some supporting results have been achieved in this study.
2.4 Correlation analysis
SEM was obtained to identify the interrelation between water quality and microorganisms (Fig. 5). Temperature had an effect on bacterial community structure (std. coeff = 0.14, P < 0.05), biomass (std. coeff = 0.17), and environmental factors (std. coeff = -0.54, P < 0.001) in all water samples. This result indicated that temperature might be the potential factor for the variations of bacterial community and ICC observed. This is in consistent with the result of a previous study by Zlatanović et al. [17], in which temperature was reported as one of the crucial controlling factors on growth and reproduction had been proved to occur an effective interrelation with bacterial counts and bacterial community structure in different ecosystems.
Other environmental factors were significant negatively correlated with the bacterial community structure (std. coeff = -0.16) and bacterial cell count (std. coeff = -0.40, P < 0.05). Nevertheless, the bacterial cell count was positively correlated with the bacterial community structure (std. coeff = 0.08) and this is in agreement with the finding of Zhang et al. [18], which reportrd that stagnation was not only induced the number of bacteria cells, but also enriched the diversity of bacterial community structure.
RDA analysis was conducted to determine the characteristic environmental elements affecting the bacterial community at the genus level. RDA1 and RDA2 accounted for 41.2% and 15.6% of the total variation in Fig. 6, respectively (P < 0.05). Moreover, temperature was considered as the most significant explanatory variable. In addition, RDA results showed the obvious clustering between stagnant and fresh water, indicating the significant differences in bacterial community structures between stagnant and fresh water.
The RDA results exhibited that the bacterial community structure was significantly influenced by temperature, which was consistent with the result observed by Zhang et al [18] and it was similar to the SEM results presented. Furthermore, chlorine and Fe concentrations also strongly influenced the bacterial community and this is in agreement with Baron et al. [6], which indicated that the abundance of OPPPs and denitrifying bacteria increased during the decay of chlorine. Similarly, Ji et al. [7] suggested that the total chlorine concentration and metal were associated with the variation in pipeline water microbiome composition. Recent research has proposed that the system of chlorine disinfection exhibited a positive contribution in the aspect of controlling bacteria counts, but played a limited role in controlling chlorine-resistant genes in bacteria [68]. Chlorination changed the structure of chlorine-resistant genes by affecting the relative abundance of chlorine-resistant genes, which enriched chlorine-resistant genes and their carriers. In the chlorination system, bacteria that could produce high extracellular polymeric substances (EPS) revealed stronger ability of survival, because EPS had a positive effect on maintaining the stability of cell adhesion structure and disinfection resistance [69]. The pipe material was found to be an important factor affecting the abundance of specific microorganisms and the bacterial community structure in DWDS [70]. Metal pipe materials and corrosion products could select bacteria that carries metal resistance genes through bacteria co-growth processes [13]. These results demonstrated that the variation in the environmental factors including temperature and chlorine attenuation caused significant variations in the bacterial community structure before and after stagnation.