Next Article in Journal
Population Attributable Fraction of Early Age of Onset of Alcohol Use in Alcohol Abuse and Dependence: A 3-Year Follow-Up Study in University Students
Next Article in Special Issue
Impact of Freeze–Thaw Cycles on Die-Off of E. coli and Intestinal Enterococci in Deer and Dairy Faeces: Implications for Landscape Contamination of Watercourses
Previous Article in Journal
Spatio-Temporal Variations of Satellite-Based PM2.5 Concentrations and Its Determinants in Xinjiang, Northwest of China
Previous Article in Special Issue
Legionella Detection in Water Networks as per ISO 11731:2017: Can Different Filter Pore Sizes and Direct Placement on Culture Media Influence Laboratory Results?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterisation of Microbial Community Associated with Different Disinfection Treatments in Hospital hot Water Networks

1
Department of Biomedical, Metabolic and Neural Sciences, Section of Public Health, University of Modena and Reggio Emilia, 41125 Modena, Italy
2
University Hospital Policlinico of Modena, 41124 Modena, Italy
3
Department of Movement, Human and Health Sciences, Public Health Unit, University of Rome ‘Foro Italico’, 00135 Rome, Italy
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(6), 2158; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17062158
Submission received: 19 February 2020 / Revised: 20 March 2020 / Accepted: 22 March 2020 / Published: 24 March 2020

Abstract

:
Many disinfection treatments can be adopted for controlling opportunistic pathogens in hospital water networks in order to reduce infection risk for immunocompromised patients. Each method has limits and strengths and it could determine modifications on bacterial community. The aim of our investigation was to study under real-life conditions the microbial community associated with different chemical (monochloramine, hydrogen peroxide, chlorine dioxide) and non-chemical (hyperthermia) treatments, continuously applied since many years in four hot water networks of the same hospital. Municipal cold water, untreated secondary, and treated hot water were analysed for microbiome characterization by 16S amplicon sequencing. Cold waters had a common microbial profile at genera level. The hot water bacterial profiles differed according to treatment. Our results confirm the effectiveness of disinfection strategies in our hospital for controlling potential pathogens such as Legionella, as the investigated genera containing opportunistic pathogens were absent or had relative abundances ≤1%, except for non-tuberculous mycobacteria, Sphingomonas, Ochrobactrum and Brevundimonas. Monitoring the microbial complexity of healthcare water networks through 16S amplicon sequencing is an innovative and effective approach useful for Public Health purpose in order to verify possible modifications of microbiota associated with disinfection treatments.

1. Introduction

The complexity of healthcare facilities and the increasing proportion of immunologically vulnerable patients make prevention of healthcare-associated waterborne infections a priority. Several factors make hospital buildings suitable for colonization with bacteria and moulds: large and complex water systems with areas of low flow predispose to stagnation and biofilm formation and water temperatures optimal for healthcare use may also be ideal for bacterial growth [1]. Moreover, water characteristics, age and corrosion of the pipes, or metabolic activity of colonizing bacteria can influence the microbial community [2,3,4]. Therefore, the microbial ecology of water networks varies in the water distribution systems [5] and can serve as reservoirs for waterborne pathogens such as Legionella spp., Pseudomonas aeruginosa, Stenotrophomonas maltophilia, Acinetobacter spp., and non-tuberculous mycobacteria [6,7,8,9,10,11].
A range of on-site physical and chemical treatments for secondary disinfection have been proposed with the aim of controlling microbial water contamination, differing in effectiveness, application methods (in continuous and shock), costs and management. Each of these methods has benefits and disadvantages and to date, the most effective procedure has not been defined [12,13,14]. Among physical methods, circulating water at 60 °C, such that the temperature at each outlet reaches at least 50 °C and preferably 55 °C, is the method most commonly used to control legionellae in hot water system. It is relatively easy to implement and to monitor continuously, but it has the possible disadvantage of increasing energy consumption and the risk of scalding [15]. Moreover, the temperature control regimen may fail, especially in large complex water systems with poor engineering control, leading to the necessary implementation of additional control methods, such as biocides [13,16]. Among chemical disinfection treatments, the most popular are continuous systems such chlorine-based biocides and hydrogen peroxide [17]. All these methods have proven effective against Legionella and other waterborne pathogens, but active concentrations of biocide need to be continuously monitored since no one eliminates the bacteria once the water network is contaminated [12,13,18,19]. The main disadvantages associated with chemical treatments are corrosiveness and, limited to chlorine-based biocides, formation of toxic disinfection by-products (DBPs) [3,12,13,20].
Wang et al. [4] suggested that disinfection treatments can modify the composition of microbial community of drinking water distribution systems both in terms of abundance and diversity at the genus level. Additionally, the effectiveness of disinfection in removing pathogens from drinking water is mediated by the microbial ecology of drinking water system [5]. For instance, the growth of nitrifying bacteria in chloraminated building plumbing systems can contribute to accelerated chloramine residuals decay [21]. Moreover, it is widely accepted that bacteria like Legionella spp. and Mycobacterium spp. embedded in biofilm are more resistant to disinfection than individual planktonic cells [22,23].
It is difficult to cultivate most environmental microorganisms under laboratory conditions, thus limiting our knowledge of the whole bacterial community [24,25]. In the last decade, the use of next generation sequencing (NGS) and bioinformatics tools has offered a more extensive approach for examining the culturable and unculturable microorganisms, allowing to characterize the entire microbial community in an ecological niche [26,27,28,29,30].
To date, there is still a notable lack of studies concerning the effect of disinfection on the microbial community present in the healthcare water networks where secondary disinfection is frequently adopted in order to reduce the risk of waterborne infections. Since drinking water distribution systems are designed and maintained to transport chemically and biologically safe potable water to consumers [31] and treatments can determine important modification in the microbiota structure present in the water networks [29], the aim of our investigation was to study under real-life conditions the microbial community associated with different chemical (monochloramine, hydrogen peroxide, chlorine dioxide) and non-chemical (hyperthermia) treatments, continuously applied since many years in four hot water networks of the same hospital, by using 16S amplicon sequencing approach.

2. Materials and Methods

2.1. Water Distribution Systems

The study was carried out in the 621-bedded University Hospital Policlinico of Modena, Italy, consisting of separate buildings constructed between 1970s and 1990s. The incoming cold groundwater, disinfected with chlorine dioxide, is provided by the municipality. A single main cold water line branches out into several supply lines, each of which enters the water station of each hospital building. Once inside, cold water is heated using heat exchangers in order to produce hot water which is moved in a water recirculation system. No secondary disinfection treatment and no storage are applied on the cold water. Continuous hydrogen peroxide, monochloramine and chlorine dioxide systems are used since many years as secondary disinfection systems on hot water lines of three different buildings defined as building 1, 2, and 3, respectively. Hydrogen peroxide device is operating since 2012, monochloramine and chlorine dioxide devices since 2009 [13]. A 48% hydrogen peroxide solution (O2 s.r.l., Bergamo, Italy) is injected continuously into the recirculating hot water by a dosing pump in order to ensure concentrations of 15 to 20 mg/L at distal outlets. Monochloramine is generated in situ by reaction between a stabilized chlorine-based precursor and an ammonium salt (Sanipur s.r.l., Brescia, Italy). The monochloramine generator is set to keep a concentration of biocide in the recirculation loop between 1.5 and 3.0 mg/L. Monochloramine residual levels are in line with the guideline value of 3 mg/L and the maximum contaminant level of 4 mg/L, established by WHO and the U.S. Environmental Protection Agency, respectively. Chlorine dioxide is produced in situ by injecting hydrochloric acid and sodium chlorite (Sanipur s.r.l., Brescia, Italy). A dosing system is set up to add the disinfectant to the recirculating hot water, ensuring concentrations of at least 0.3 mg/L at distal points. This dose was reported as effective in controlling Legionella contamination in European hospitals [32].
Hot water in a fourth building (building 4) is not chemically treated but is produced at a temperature of at least 60 °C and distributed continuously at >50 °C (hyperthermia). This hot water network supplies medical studies and laboratories, while the hospital rooms have instantaneous water heaters which are directly fed with cold water and maintained at 60 °C, for Legionella contamination control.

2.2. Samples Collection

One inlet cold water and six hot water samples were collected from every building. The cold water sample was collected within the water station at the entry of each building, while hot water was taken from different sites throughout each building (hot water heater, recirculation loop and taps). In detail, hot water samples were collected from the water network treated with hydrogen peroxide (HP) within the building 1, from the network treated with monochloramine (M) within the building 2, from the network treated with chlorine dioxide (CD) within the building 3, and from the not chemically treated network (NT) within the building 4.
Water samples were collected during two sampling sessions (February and May 2018). In February, the hot water samples were collected in duplicate (called “a” and “b”) in order to verify reproducibility of methodology. Due to the positive results, there were no replicates in May.
Two litres of hot water and five litres of inlet cold water were collected in sterile plastic bottles without flaming and after 1 min flushing. Sodium thiosulphate was added in order to neutralize residual free chlorine or other oxidizing biocides. Water temperature and residual disinfectant concentration were determined during samples collection as explained: temperature (digital thermometer), free chlorine and chlorine dioxide (DPD method, Microquant, Merck, Darmstadt, Germany), monochloramine (DR900 colorimeter, indophenol method, Hach Lange, Milan, Italy), hydrogen peroxide (reflectometer RQflex 2, Merck, Darmstadt, Germany). Collection bottles were returned to the laboratory immediately after sampling.
The University Hospital gave approval and permission to collect water samples for this study.

2.3. Characterization of Microbial Community

Water samples were concentrated by membrane filtration (0.40 μm polycarbonate membrane, extraction DNA kit, Minerva Biolabs, Germany). The filter membrane was stored in a Petri plate at −20 °C. The filter membrane was turned upside down onto an incubation dish filled with 2 mL of Lysis buffer (Minerva BioLabs, Germany) and then incubated at 37 °C for 30 min. Subsequently, the lysis solution was transferred into an incubation tube with 0.1 mg of glass beads (Sigma Aldrich, USA) and incubated at 56 °C for 15 min after vortexing for 1 min. DNA purification was carried out using Aqua screen Fast Extract Kit according to the manufacturer’s protocol (Minerva BioLabs, Germany).

Bacterial Community 16S Profiling

16S rRNA paired-end sequencing was performed according to the “16S Metagenomic Sequencing Library Preparation” protocol (Part# 15044223 rev. A; Illumina, USA). Briefly, the primers containing Illumina adapter and linker sequence and targeting the V1-V2 regions of bacterial 16S rRNA genes were used [30,33,34]. Three libraries with unique tags were generated for each sample as technical replicates. Each amplification reaction had a total volume of 25 μL containing 12.5 μL of KAPA HiFi HotStart ReadyMix (Roche, Pleasanton, CA, USA), 5 μL of each primer (1 μM), and 2 μL template DNA. Reactions were carried out on a Techne®TC-PLUS thermalcycler (VWR International, LLC, Radnor, USA). Following amplification, 5 μL of PCR product from each reaction was used for agarose gel (1%) electrophoresis to confirm amplification. The final concentration of cleaned DNA amplicon was determined using the Qubit PicoGreen dsDNA BR assay kit (Invitrogen, Grand Island, NY, USA) and validated on Bioanalyzer DNA 1000 chip (Agilent, Santa Clara, CA, USA).
Due to failed DNA amplification, 10 out of 80 samples collected were excluded. Failed DNA amplification may have been due to multiple factors: low DNA concentrations, the presence of PCR inhibitors (i.e., humic substances, phenolic compounds) and/or extensive DNA damage caused by high levels of disinfectant (i.e., chlorine) [31,35,36].
Libraries were prepared using the MiSeq Reagent Kit Preparation Guide (Illumina, San Diego, CA, USA). Raw sequence data was processed using an in-house pipeline which was built on the Galaxy platform and incorporated various software tools to evaluate the quality of the raw sequence data (e.g., FastQC, http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). All data sets were rigorously screened to remove low quality reads (short reads <200 nt, zero-ambiguous sequences). Demultiplexing was performed to remove PhiX sequences and sort sequences; moreover, to minimize sequencing errors and ensure sequence quality, the reads were trimmed based on the sequence quality score using Btrim [37]. OTUs (operational taxonomic units) were clustered at a 97% similarity level and final OTUs were generated based on the clustering results. Taxonomic annotation of individual OTUs was based on representative sequences using RDP’s 16S Classifier 2.5. Observed OTUs were defined as observed species. The sequence reads were analysed, also, in the cloud environment BaseSpace through the 16S Metagenomics app (version 1.0.1 (Illumina®, San Diego, CA, USA)): the taxonomic database used was the Illumina-curated version (May 2013 release of the Greengenes Consortium Database [38]).
The raw sequencing data have been submitted to NCBI Sequence Read Archive (http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra/) with the project accession number of PRJNA574025.

2.4. Data Analysis

Relative abundances of community members were determined with rarefied data and summarized at each taxonomic level. Alpha and Beta diversity were calculated using EstimateS software at a level of 97% sequence similarity. Regarding Alpha diversity, microbial richness was computed based on number of OTUs observed and biodiversity was estimated through Shannon index and equitability index at species level [39]. Morisita Horn index (Beta diversity) was calculated in order to evaluate similarity/dissimilarity between two samples or two pools [40,41,42]. The index ranges from 0 (no similarity) to 1 (complete similarity). Principal components analysis (PCA) was performed using METAGENassist platform [43], in order to investigate the dissimilarity between groups. Moreover, Beta diversity clustering was analysed using ANOSIM (ANalysis Of SIMilarities) and it was run to assess significant differences in relative abundance of OTUs on different samples [44]. In order to assess sequencing depth, alpha rarefaction plots were done in mothur (version 1.31.1, www.mothur.org) and R (version 3.1.3, www.R-project.org) using packages ‘ggplot2’ and ‘vegan’ (R Core team 2013).
All statistical analyses were performed with software package IBM SPSS statistics, version 22 (IBM Corp., Armonk, NY, USA) Comparison for normal variables was performed by t-test and one-way analysis of variance (ANOVA) with the Bonferroni test and for non-normal variables (OTUs) by non-parametric test (Mann–Whitney test and Kruskal–Wallis test). The Pearson correlation coefficient and the Spearman’s rank correlation coefficient were calculated to analyse the correlation between normal variables and between non-normal variables, respectively.

3. Results

A total of 70 samples were analysed (8 cold water and 62 hot water samples). Forty-six samples (65.7%) were collected in February during the first sampling session and the remaining in May (second session). Rarefaction curves were calculated for each sample (Figure S1). All the curves reached a stable plateau. Separately amplified and barcoded replicates of hot water samples from the first sampling session (n. 19 couples) were sequenced to verify reproducibility of procedure. Morisita Horn analyses of replicates demonstrated high levels of community similarity at genera level, with an average equal to 0.904 ± 0.130 and values higher than 0.850 for 89% of replicates (n. 17 couples) and higher than 0.540 for the remaining 11% of the replicates (n. 2 couples).

3.1. Cold Water Samples

Eight cold water samples were analysed, 4 from the first and 4 from the second sampling session. Table 1 shows the values of temperature, OTUs, Shannon and equitability indices for each sample. No significant differences were observed for waters collected in February compared to those collected in May. The Principal Components Analysis showed common profiles for all the cold waters (Figure 1). The similarity of cold waters from four different buildings was supported by ANOSIM (R value = 0.208, p = 0.112). In cold water samples, the predominant phyla were Proteobacteria (mean 76.30%), Bacteroidetes (mean 11.50%) and Actinobacteria (mean 7.75%). Figure 2 shows the predominant genera in inlet cold water to the four buildings. Burkholderia (mean 15.82%), Bradyrhizobium (mean 10.11%) and Sediminibacterium (mean 9.24%) were the predominant and common genera for all cold water samples.

3.2. Hot Water Samples

A total of 62 hot water samples were analysed, 14 treated with hydrogen peroxide (HP), 17 with monochloramine (M), 15 with chlorine dioxide (CD), and 16 not chemically treated (NT). The mean temperature of the samples collected during the first sampling session (n. 42) did not differed from the mean of those collected during the second one (n. 20) (46.0 ± 5.6 vs 43.2 ± 7.0 °C).
Table 2 presents water temperatures, biocide concentrations, Shannon and equitability indices and number of OTUs. Samples from network NT had significantly higher mean temperature than the other networks (HP, M, CD). Kruskal–Wallis test showed a global statistically significant difference in terms of OTUs (p = 0.008). Median values of OTUs resulted significantly higher in not chemically treated network than in networks treated with hydrogen peroxide and monochloramine (p = 0.001 NT vs HP; p = 0.038 NT vs M; p = 0.286 NT vs CD, Mann–Whitney test), while no significant differences in OTUs median were detected between the three chemically treated hot networks (p = 0.115, Kruskal–Wallis test). No significant differences were observed between the four hot networks in terms of alpha diversity expressed as Shannon (p = 0.481) and Equitability (p = 0.254) indices (ANOVA). Also limiting the analysis to the three networks chemically treated, ANOVA tests showed no significant differences in terms of Shannon index (p = 0.385) and Equitability index (p = 0.228). Moreover, we observed a negative, although not significant, correlation between the number of OTUs and biocide concentration, with Spearman’s Rho equal to −0.225 (p = 0.440) for hydrogen peroxide, −0.153 (p = 0.557) for monochloramine and −0.318 (p = 0.248) for chlorine dioxide.
Figure 3 shows the distribution of the relative abundance of phyla for each hot water sample. The predominant phylum was Proteobacteria (69.31%), followed by Actinobacteria (8.83%), Bacteroidetes (8.77%), Firmicutes (5.22%), Thermi (1.37%), and Acidobacteria (1.36%). The mean relative abundance of Proteobacteria was higher in the samples from network M (82.13%) and network CD (83.58%) compared to samples from network HP (57.85%) or network NT (52.34%). Conversely, the phylum Actinobacteria was prevalent in the samples from network HP (29.87%) compared to those from network M (3.03%) or network CD (4.28%) or network NT (0.89%). The phyla Bacteroidetes and Firmicutes had a higher mean relative abundance in not chemically treated hot water samples (12.87% and 14.52%, respectively) compared to the chemically treated ones. The mean relative abundance of Thermi was higher in the samples from network HP (3.23%) compared to those from other networks. Relative abundance of Acidobacteria was higher in network CD (2.37%) and network NT (2.55%) than other networks.
The genera distribution of samples pooled for the four hot water networks is illustrated in Figure 4 and the three predominant genera for each network are shown in Table 3. Some genera resistant to unfavourable environmental conditions were predominant in the presence of biocides, such as Bradyrhizobium (networks treated with hydrogen peroxide and chlorine dioxide), Mycobacterium (hydrogen peroxide), Gallionella (monochloramine) and Blastomonas (chlorine dioxide). The thermophilic genera (Thermobaculum) predominated in network NT. The Morisita Horn indices were calculated by comparing the four pools in order to evaluate the similarity/dissimilarity of microbiological community. This analysis revealed different microbiological profiles at genera level according to the water networks with the average value of the Morisita Horn indices equal to 0.221 (range 0.118−0.294). The microbiological similarity of the hot water samples from the same network emerged from the Principal Components Analysis, as shown in Figure 5. Moreover, ANOSIM was calculated by comparing the four hot waters pools. The result indicated a higher dissimilarity between hot waters from different networks than those within the same network (R value = 0.621, p = 0.001).
In networks HP and CD there was no association between the biocide concentration and genera relative abundances, while Sediminibacterium negatively correlated with the monochloramine concentration in network M (Pearson = −0.580, p = 0.015).

3.3. Genera Containing Opportunistic Pathogens

Sequence data was further analysed in order to observe the differences in genera containing opportunistic pathogens among the different hot water networks. Escherichia, Shigella, Salmonella, Klebsiella and Enterobacter genera were investigated but not detected in any samples. Genera detected were Burkholderia, Brevundimonas, Sphingomonas, Mycobacterium, Legionella, Acinetobacter, Chryseobacterium, Ochrobactrum, Pseudomonas, and Stenotrophomonas. All these genera had relative abundances ≤1%, except for Mycobacterium in the chemically treated networks, Sphingomonas in all networks, Ochrobactrum in network M (mean 4.03%) and Brevundimonas in network CD (mean 1.07%). Mycobacterium had higher relative abundance in hot water samples treated with hydrogen peroxide (mean 36.36%) followed by chlorine dioxide (4.16%) and monochloramine (1.80%), and lower levels in NT hot water network (mean 0.07%). Mean relative abundance of Sphingomonas was <2% in all networks (range 1.05%–1.71%).
We also investigated species for each genus containing opportunistic pathogen with mean relative abundance above 1% in hot water (Table S1). We identified 46 species of non-tuberculous Mycobacterium, 28 species of Sphingomonas, none of which associated with human infection, 6 species of Ochrobactrum, among which O. anthropi, and 5 species of Brevundimonas, included B. diminuta.

4. Discussion

In this investigation, we evaluated under real-life conditions the microbial community associated with different chemical and non-chemical (hyperthermia) treatments, by using 16S amplicon sequencing approach. We analysed both cold and hot waters collected in four buildings of our University hospital, all receiving the same incoming cold water but differing in secondary disinfection treatments continuously applied to hot water since many years in order to reduce the risk of waterborne infections, as reported in previous papers [13,17,45].
Similar to other studies [31,46], the inlet cold water provided by the municipality presented a common profile in all hospital buildings, characterized by the same predominant phyla and the same predominant environmental ubiquitous genera. The first three predominant phyla in hot water were the same found in cold water samples, which are common constituents of water microbial community [4,47]. Interestingly, hot water networks showed a different microbiological profile at genera level in relation to treatment. In not chemically treated water network, heat-resistant genera were among the predominant, confirming that temperature influences the water microbial community [48,49]. Indeed, the water temperatures of this network were above 50 °C compared to chemically treated ones, which were kept at ≤45 °C in order to avoid the degradation of biocides. In the chemically treated networks, the genera resistant to unfavourable environmental conditions were predominant. Bradyrhizobium was mostly present in the networks treated with hydrogen peroxide and chlorine dioxide. Mycobacterium, Gallionella and Blastomonas were found among the prevalent genera in network treated with hydrogen peroxide, monochloramine and chlorine dioxide, respectively.
Bradyrhizobium belongs to the Rhizobiaceae family [50]. Rhizobiaceae includes nitrogen-fixing bacteria producing nitrogenase enzyme, which reduces molecular nitrogen (N2) in ammonia (NH3). Of note, our results showed a lower relative abundance of Bradyrhizobium in water network treated with monochloramine than the other chemically treated ones, in accordance with Gomez-Alvarez et al. [51]. The authors studied the biofilm community succession in bench-scale annular reactors simulating arrestment of chloraminated drinking water nitrification. In more detail, Bradyrhizobium dominated in early stages of biofilm formation when reduced nutrient levels found in water habitats together with minimal monochloramine residuals (<0.1 mg/l) provided the metabolically versatile Bradyrhizobium a competitive advantage to outcompete other biofilm groups. In last periods of biofilm formation, the increase of monochloramine residual from ≈0.4 to 3.0 mg/l led to a decline in Bradyrhizobium relative abundance.
Mycobacterium genus is ubiquitous in the environment, especially in water [4,46,52]. This genus had higher relative abundance in the chemically treated networks compared to the not chemically treated one, in line with other investigations [46,53,54]. These studies reported higher Mycobacterium relative abundances in the samples treated with monochloramine, while in our hospital the highest relative abundance (36.4%) have been achieved in network treated with hydrogen peroxide. The dissimilarity between chemically and not chemically treated hot waters highlights Mycobacterium resistance to biocides, which eliminate or reduce the bacterial genera sensitive to their action, allowing the proliferation of the resistant bacteria. Regarding the possible survival mechanisms, the self-produced extracellular polymeric substances (EPS) are responsible for surface adhesion and allow the biofilm formation [55], that can protect Mycobacterium from direct exposure to disinfectants [56]. Although the composition of EPS varies significantly among different bacterial species, exopolysaccharides, proteins and lipids remain the key components of EPS [57]. The mycobacterial pellicles contain large quantities of free mycolic acids [58,59], that contribute to increase the cell-surface hydrophobicity to facilitate the cell-to-cell interaction for biofilm formation. Biofilm is difficult to penetrate by the disinfectants; thus, bacteria are protected inside it. Since monochloramine is able to penetrate biofilm [23], this mechanism can explain the dissimilarity of Mycobacterium relative abundance between hydrogen peroxide and monochloramine treated networks in our hospital.
In network with monochloramine treatment, Gallionella was the predominant genus. It is an iron oxidizing bacterium, aerobic microorganism that derive energy from the oxidation of iron from the ferrous state (Fe2+) to the ferric state (Fe3+) by accelerating the natural reaction between oxygen and ferrous ions in the water or on the surface of metal pipes such as galvanized iron, cast iron and steel pipes [60,61,62]. The presence of Gallionella has a potential negative effect on the pipes, due to the precipitation of its metabolism product (hydrated ferric hydroxide) reducing the pipes section [60]. Members of Gallionellaceae family requires also ammonia in the environment for autotrophic growth, as reported by other authors [61,63]. A previous study reported high Gallionella concentration in groundwater containing high levels of iron, manganese and ammonia [64]. Moreover, Stanish et al. reported that seasonally chloramine-treated systems had significantly higher abundances of several metal-oxidizing taxa including Gallionella spp. than chlorinated systems [65]. Since the incomplete reaction of monochloramine precursors can lead the ammonia release [66], it is not surprising that Gallionella was the predominant genus in waters treated with this biocide. This information can be relevant in monitoring of this disinfection treatment. In case of increased ammonia concentrations, corrective actions have to be taken in order to ensure that the proportions of chlorine and ammonia precursors during monochloramine formation are correct.
Regarding the Blastomonas genus, the relative abundance was higher in hot water treated with chlorine dioxide, suggesting a resistance to this treatment. Other studies found the Blastomonas resistance to chlorination [67,68]. Soto-Giron et al. detected Blastomonas as the second most abundant genus, after Mycobacterium, in hospital shower hoses [69]. They have also reported that Blastomonas and Mycobacterium expressed genes associated with resistance to disinfectants, including genes encoding proteins, such as catalase and superoxide dismutase, able to confer protection against oxidative stress [69,70].
A lower microbial richness (OTUs number) was observed in chemically treated hot water than in not-chemically treated ones. It is possible that disinfectants, killing some selected vulnerable genera, allowed an increase of dominance of the disinfectant-resistant microbial groups, thereby decreasing OTUs counts. Indeed, the negative, although not significant, correlation between the OTUs number and each biocide concentration suggests a decrease number of bacterial species with the increasing disinfectant concentration. A recent study about monochloramine impact on the microbiota of a hospital water network confirmed a lower bacterial richness in treated samples compared to control [71].
The overall microbial community composition was similar within each water network during our study as confirmed by PCA analysis, thus suggesting that the observed differences may be influenced by the treatment. The selective pressure of disinfectants on the microbial community has been observed also in other studies [4,54,72,73].
Several waterborne pathogen-containing genera were examined in order to observe the differences in relative abundance due to the disinfectant in use. The concentrations of biocides were at levels effective in reducing Legionella contamination [17,74]. Our results confirmed the effectiveness of biocides used in our hospital for controlling the presence of potential pathogens such as Legionella, as most of the investigated genera containing opportunistic pathogens were absent or present at very low level of relative abundance (≤1%), even if it cannot be excluded that low levels of pathogenic bacteria may pose a potential risk to public health. The increased relative abundance of non-tuberculous mycobacteria, especially in waters treated with hydrogen peroxide underlines that there exists the potential for unwanted consequences of supplemental disinfectant addition. Understanding the impact of secondary disinfection on water system microbial community is necessary to maximize biocide effectiveness and to ensure that supplemental disinfectant does not select for alternative opportunistic pathogens [53], also considering that other system factors such as temperature and pipe materials may impact microbial ecology [4]. In addition, although the 16S amplicon sequencing method has the limitation of not providing information about live/dead status, it can help identify these “hot spots” that could be acted upon with culture-based methods which are likely to remain the gold standard. In a context of risk prevention and rapid screening, molecular methods allow to support more extensive, timely, and economical evaluation of problematic sites that could be followed-up upon with culturing [75].
A limitation of this study is the absence of a pre-post treatment approach, not applicable as the programme to control Legionella contamination in the hospital’s hot water distribution systems was implemented since 2000 in accordance with mandatory National Guidelines. At that time, the hospital buildings were highly contaminated and over years the most promising disinfection treatments have been selected to better control the contamination, as reported in a previous paper [14]. In addition, other variables would be needed to better understand building-to-building variations, but this study aims to be one of the first approaches under real-life conditions. Most of the prior research in this area investigated bench or pilot-scale simulated water distribution systems [4,31,51,54], where a rigorous control of numerous variables such as pipe materials, degree of stagnation, average ambient temperature is possible. However, important strengths of our study are to investigate simultaneously microbiota of four different hospital hot water networks and to look at the incoming cold water for microbial profile as well as longitudinally at the system over different seasons in order to observe how this affects microbial communities. The high level of similarity among replicates is an evidence of reproducibility of our procedure.
Our research provides novel knowledge about the influence of disinfection on microbial community within water distribution systems under real-life conditions, an issue still scarcely investigated that requires further studies to better comprehend the treatments’ selective pressure on microbiota of water networks in healthcare settings. As Wang et al. [76] stated, research is needed to better understand how engineering controls may individually or in combination select a desirable microbiome, and how the microbiome itself may mediate proliferation of opportunistic pathogens, suggesting a probiotic approach to pathogens control in premise plumbing system.

5. Conclusions

We found different microbial profiles associated with four hot water networks. These findings may be related to the selective pressures exerted by disinfection systems on microbial communities. Monitoring the microbial complexity of healthcare water networks through 16S amplicon sequencing is an innovative and effective approach in its application in Public Health in order to verify possible modifications of microbiota associated with disinfection.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1660-4601/17/6/2158/s1. Figure S1: Rarefaction curves of the bacterial 16S rRNA gene sequences from the 70 water samples based on OTUs calculated at 97% identity. The total number of sequences analysed is plotted against the number of OTUs observed in the same community. For hot water samples: HP: network with hydrogen peroxide, M: network with monochloramine, CD: network with chlorine dioxide and NT: not chemically treated network. For cold water samples: CW: Cold Water, Building is indicated as 1, 2, 3 and 4, F: February and M: May, Table S1: Identified species of investigated genera containing opportunistic pathogens in hot waters.

Author Contributions

Conceptualization, A.B., S.P., P.B. and I.M.; methodology, I.M. and S.P.; formal analysis, F.V., M.E.C. and S.P.; investigation, I.M. and G.F. and S.P.; resources, I.M., E.V. and L.S.; data curation, M.E.C., F.V. and V.R.S.; writing—original draft preparation, S.P., M.E.C. and A.B.; writing—review and editing, A.B., S.P. and I.M.; supervision, A.B., P.B. and V.R.S.; project administration, A.B.; funding acquisition, A.B. and P.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Decker, B.K.; Palmore, T.N. Hospital water and opportunities for infection prevention. Curr. Infect. Dis. Rep. 2014, 16, 432. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Bargellini, A.; Marchesi, I.; Righi, E.; Ferrari, A.; Cencetti, S.; Borella, P.; Rovesti, S. Parameters predictive of Legionella contamination in hot water systems: Association with trace elements and heterotrophic plate counts. Water Res. 2011, 45, 2315–2321. [Google Scholar] [CrossRef] [PubMed]
  3. Giovanardi, R.; Bononi, M.; Messori, M.; Bargellini, A.; Paduano, S.; Borella, P.; Marchesi, I. Corrosion resistance of commonly used plumbing materials for water distribution systems exposed to disinfection treatments. Corros. Eng. Sci. Technol. 2020, 1–8. [Google Scholar] [CrossRef]
  4. Wang, H.; Masters, S.; Edwards, M.A.; Falkinham, J.O.; Pruden, A. Effect of disinfectant, water age, and pipe materials on bacterial and eukaryotic community structure in drinking water biofilm. Environ. Sci. Technol. 2014, 48, 1426–1435. [Google Scholar] [CrossRef]
  5. Berry, D.; Xi, C.; Raskin, L. Microbial ecology of drinking water distribution systems. Curr. Opin. Biotechnol. 2006, 17, 297–302. [Google Scholar] [CrossRef]
  6. Anaissie, E.J.; Penzak, S.R.; Dignani, M.C. The hospital water supply as a source of nosocomial infections: A plea for action. Arch. Intern. Med. 2002, 162, 1483–1492. [Google Scholar] [CrossRef] [Green Version]
  7. Cervia, J.S.; Ortolano, G.A.; Canonica, F.P. Hospital tap water as a source of Stenotrophomonas maltophilia infection. Clin. Infect. Dis. 2008, 46, 1485–1487. [Google Scholar] [CrossRef]
  8. Falkinham, J.O., 3rd. Common features of opportunistic premise plumbing pathogens. Int. J. Environ. Res. Public Health 2015, 12, 4533–4545. [Google Scholar] [CrossRef]
  9. Kusnetsov, J.; Torvinen, E.; Perola, O.; Nousiainen, T.; Katila, M.L. Colonization of hospital water systems by legionellae, mycobacteria and other heterotrophic bacteria potentially hazardous to risk group patients. APMIS 2003, 111, 546–556. [Google Scholar] [CrossRef]
  10. Merlani, G.M.; Francioli, P. Established and emerging waterborne nosocomial infections. Curr. Opin. Infect. Dis. 2003, 16, 343–347. [Google Scholar] [CrossRef]
  11. Wang, J.L.; Chen, M.L.; Lin, Y.E.; Chang, S.C.; Chen, Y.C. Association between contaminated faucets and colonization or infection by nonfermenting gram-negative bacteria in intensive care units in Taiwan. J. Clin. Microbiol. 2009, 47, 3226–3230. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Lin, Y.E.; Stout, J.E.; Yu, V.L. Controlling Legionella in Hospital Drinking Water: An Evidence-Based Review of Disinfection Methods. Infect. Control Hosp. Epidemiol. 2011, 32, 166–173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Marchesi, I.; Ferranti, G.; Mansi, A.; Marcelloni, A.M.; Proietto, A.R.; Saini, N.; Borella, P.; Bargellini, A. Control of Legionella Contamination and Risk of Corrosion in Hospital Water Networks following Various Disinfection Procedures. Appl. Environ. Microbiol. 2016, 82, 2959–2965. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Marchesi, I.; Marchegiano, P.; Bargellini, A.; Cencetti, S.; Frezza, G.; Miselli, M.; Borella, P. Effectiveness of different methods to control legionella in the water supply: Ten-year experience in an Italian university hospital. J. Hosp. Infect. 2011, 77, 47–51. [Google Scholar] [CrossRef] [PubMed]
  15. The European Centre for Disease Prevention and Control (ECDC). European Technical Guidelines for the Prevention, Control and Investigation of Infections Caused by Legionella Species; ECDC: Solna, Sweden, 2017. [Google Scholar]
  16. Darelid, J.; Lofgren, S.; Malmvall, B.E. Control of nosocomial Legionnaires’ disease by keeping the circulating hot water temperature above 55 degrees C: Experience from a 10-year surveillance programme in a district general hospital. J. Hosp. Infect. 2002, 50, 213–219. [Google Scholar] [CrossRef]
  17. Borella, P.; Bargellini, A.; Marchegiano, P.; Vecchi, E.; Marchesi, I. Hospital-acquired Legionella infections: An update on the procedures for controlling environmental contamination. Ann. Ig. 2016, 28, 98–108. [Google Scholar] [CrossRef]
  18. Coniglio, M.A.; Ferrante, M.; Yassin, M.H. Preventing Healthcare-Associated Legionellosis: Results after 3 Years of Continuous Disinfection of Hot Water with Monochloramine and an Effective Water Safety Plan. Int. J. Environ. Res. Public Health 2018, 15, 1594. [Google Scholar] [CrossRef] [Green Version]
  19. Casini, B.; Buzzigoli, A.; Cristina, M.L.; Spagnolo, A.M.; Del Giudice, P.; Brusaferro, S.; Poscia, A.; Moscato, U.; Valentini, P.; Baggiani, A.; et al. Long-term effects of hospital water network disinfection on Legionella and other waterborne bacteria in an Italian university hospital. Infect. Control Hosp. Epidemiol. 2014, 35, 293–299. [Google Scholar] [CrossRef]
  20. Hua, G.; Reckhow, D.A. Comparison of disinfection byproduct formation from chlorine and alternative disinfectants. Water Res. 2007, 41, 1667–1678. [Google Scholar] [CrossRef]
  21. Rhoads, W.J.; Pruden, A.; Edwards, M.A. Interactive Effects of Corrosion, Copper, and Chloramines on Legionella and Mycobacteria in Hot Water Plumbing. Environ. Sci. Technol. 2017, 51, 7065–7075. [Google Scholar] [CrossRef]
  22. Abdel-Nour, M.; Duncan, C.; Low, D.E.; Guyard, C. Biofilms: The stronghold of Legionella pneumophila. Int. J. Mol. Sci. 2013, 14, 21660–21675. [Google Scholar] [CrossRef] [Green Version]
  23. Lee, W.H.; Pressman, J.G.; Wahman, D.G. Three-Dimensional Free Chlorine and Monochloramine Biofilm Penetration: Correlating Penetration with Biofilm Activity and Viability. Environ. Sci. Technol. 2018, 52, 1889–1898. [Google Scholar] [CrossRef]
  24. Badhai, J.; Ghosh, T.S.; Das, S.K. Taxonomic and functional characteristics of microbial communities and their correlation with physicochemical properties of four geothermal springs in Odisha, India. Front. Microbiol. 2015, 6, 1166. [Google Scholar] [CrossRef] [PubMed]
  25. Paduano, S.; Valeriani, F.; Romano-Spica, V.; Bargellini, A.; Borella, P.; Marchesi, I. Microbial biodiversity of thermal water and mud in an Italian spa by metagenomics: A pilot study. Water Sci. Technol. Water Supply 2018, 18, 1456–1465. [Google Scholar] [CrossRef]
  26. Amin, A.; Ahmed, I.; Salam, N.; Kim, B.-Y.; Singh, D.; Zhi, X.-Y.; Xiao, M.; Li, W.-J. Diversity and Distribution of Thermophilic Bacteria in Hot Springs of Pakistan. Microb. Ecol. 2017, 74, 116–127. [Google Scholar] [CrossRef] [PubMed]
  27. Bautista-de Los Santos, Q.M.; Schroeder, J.L.; Blakemore, O.; Moses, J.; Haffey, M.; Sloan, W.; Pinto, A.J. The impact of sampling, PCR, and sequencing replication on discerning changes in drinking water bacterial community over diurnal time-scales. Water Res. 2016, 90, 216–224. [Google Scholar] [CrossRef] [PubMed]
  28. Li, Y.; Zhang, J.; Zhang, J.; Xu, W.; Mou, Z. Microbial Community Structure in the Sediments and Its Relation to Environmental Factors in Eutrophicated Sancha Lake. Int. J. Environ. Res. Public Health 2019, 16, 1931. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Proctor, C.R.; Hammes, F. Drinking water microbiology--from measurement to management. Curr. Opin. Biotechnol. 2015, 33, 87–94. [Google Scholar] [CrossRef]
  30. Valeriani, F.; Crognale, S.; Protano, C.; Gianfranceschi, G.; Orsini, M.; Vitali, M.; Spica, V.R. Metagenomic analysis of bacterial community in a travertine depositing hot spring. New Microbiol. 2018, 41, 126–135. [Google Scholar]
  31. Potgieter, S.; Pinto, A.; Sigudu, M.; du Preez, H.; Ncube, E.; Venter, S. Long-term spatial and temporal microbial community dynamics in a large-scale drinking water distribution system with multiple disinfectant regimes. Water Res. 2018, 139, 406–419. [Google Scholar] [CrossRef] [Green Version]
  32. Sidari, F.P.; Stout, J.E.; Vanbriesen, J.M.; Bowman, A.M.; Grubb, D.; Neuner, A.; Wagener, M.M.; Yu, V.L. Keeping Legionella out of water systems. J. Am. Water Works Assoc. 2004, 96, 111–119. [Google Scholar] [CrossRef]
  33. Kittelmann, S.; Seedorf, H.; Walters, W.A.; Clemente, J.C.; Knight, R.; Gordon, J.I.; Janssen, P.H. Simultaneous amplicon sequencing to explore co-occurrence patterns of bacterial, archaeal and eukaryotic microorganisms in rumen microbial communities. PLoS ONE 2013, 8, e47879. [Google Scholar] [CrossRef] [PubMed]
  34. Valeriani, F.; Protano, C.; Gianfranceschi, G.; Leoni, E.; Galasso, V.; Mucci, N.; Vitali, M.; Spica, V.R. Microflora Thermarum Atlas project: Biodiversity in thermal spring waters and natural SPA pools. Water Sci. Technol. Water Supply 2018, 18, 1472–1483. [Google Scholar] [CrossRef]
  35. Schrader, C.; Schielke, A.; Ellerbroek, L.; Johne, R. PCR inhibitors—Occurrence, properties and removal. J. Appl. Microbiol. 2012, 113, 1014–1026. [Google Scholar] [CrossRef]
  36. Van Aken, B.; Lin, L.S. Effect of the disinfection agents chlorine, UV irradiation, silver ions, and TiO2 nanoparticles/near-UV on DNA molecules. Water Sci. Technol. 2011, 64, 1226–1232. [Google Scholar] [CrossRef]
  37. Kong, Y. Btrim: A fast, lightweight adapter and quality trimming program for next-generation sequencing technologies. Genomics 2011, 98, 152–153. [Google Scholar] [CrossRef] [Green Version]
  38. Wang, S.; Hou, W.; Dong, H.; Jiang, H.; Huang, L.; Wu, G.; Zhang, C.; Song, Z.; Zhang, Y.; Ren, H.; et al. Control of temperature on microbial community structure in hot springs of the Tibetan Plateau. PLoS ONE 2013, 8, e62901. [Google Scholar] [CrossRef] [Green Version]
  39. Daly, A.J.; Baetens, J.M.; De Baets, B. Ecological Diversity: Measuring the Unmeasurable. Mathematics 2018, 6, 119. [Google Scholar] [CrossRef] [Green Version]
  40. Chao, A.; Chazdon, R.L.; Colwell, R.K.; Shen, T.-J. Abundance-based similarity indices and their estimation when there are unseen species in samples. Biometrics 2006, 62, 361–371. [Google Scholar] [CrossRef]
  41. Horn, H.S. Measurement of "Overlap" in Comparative Ecological Studies. Am. Nat. 1966, 100, 419–424. [Google Scholar] [CrossRef]
  42. Wolda, H. Similarity indices, sample size and diversity. Oecologia 1981, 50, 296–302. [Google Scholar] [CrossRef] [PubMed]
  43. Arndt, D.; Xia, J.; Liu, Y.; Zhou, Y.; Guo, A.C.; Cruz, J.A.; Sinelnikov, I.; Budwill, K.; Nesbø, C.L.; Wishart, D.S. METAGENassist: A comprehensive web server for comparative metagenomics. Nucleic Acids Res. 2012, 40, W88–W95. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Clarke, K.R.; Warwick, R.M. Similarity-based testing for community pattern: The two-way layout with no replication. Mar. Biol. 1994, 118, 167–176. [Google Scholar] [CrossRef]
  45. Marchesi, I.; Ferranti, G.; Bargellini, A.; Marchegiano, P.; Predieri, G.; Stout, J.E.; Borella, P. Monochloramine and chlorine dioxide for controlling Legionella pneumophila contamination: Biocide levels and disinfection by-product formation in hospital water networks. J. Water Health 2013, 11, 738–747. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Gomez-Alvarez, V.; Revetta, R.P.; Santo Domingo, J.W. Metagenomic analyses of drinking water receiving different disinfection treatments. Appl. Environ. Microbiol. 2012, 78, 6095–6102. [Google Scholar] [CrossRef] [Green Version]
  47. Van Assche, A.; Crauwels, S.; De Brabanter, J.; Willems, K.A.; Lievens, B. Characterization of the bacterial community composition in water of drinking water production and distribution systems in Flanders, Belgium. MicrobiologyOpen 2018, 8, e726. [Google Scholar] [CrossRef]
  48. Botero, L.M.; Brown, K.B.; Brumefield, S.; Burr, M.; Castenholz, R.W.; Young, M.; McDermott, T.R. Thermobaculum terrenum gen. nov., sp. nov.: A non-phototrophic gram-positive thermophile representing an environmental clone group related to the Chloroflexi (green non-sulfur bacteria) and Thermomicrobia. Arch. Microbiol. 2004, 181, 269–277. [Google Scholar] [CrossRef]
  49. Buse, H.Y.; Ji, P.; Gomez-Alvarez, V.; Pruden, A.; Edwards, M.A.; Ashbolt, N.J. Effect of temperature and colonization of Legionella pneumophila and Vermamoeba vermiformis on bacterial community composition of copper drinking water biofilms. Microb. Biotechnol. 2017, 10, 773–788. [Google Scholar] [CrossRef]
  50. Berrada, H.; Fikri-Benbrahim, K. Taxonomy of the Rhizobia: Current Perspectives. Br. Microbiol. Res. J. 2014, 4, 616–639. [Google Scholar] [CrossRef] [Green Version]
  51. Gomez-Alvarez, V.; Schrantz, K.A.; Pressman, J.G.; Wahman, D.G. Biofilm community dynamics in bench-scale annular reactors simulating arrestment of chloraminated drinking water nitrification. Environ. Sci. Technol. 2014, 48, 5448–5457. [Google Scholar] [CrossRef]
  52. September, S.M.; Brözel, V.S.; Venter, S.N. Diversity of nontuberculoid Mycobacterium species in biofilms of urban and semiurban drinking water distribution systems. Appl. Environ. Microbiol. 2004, 70, 7571–7573. [Google Scholar] [CrossRef] [Green Version]
  53. Baron, J.L.; Vikram, A.; Duda, S.; Stout, J.E.; Bibby, K. Shift in the microbial ecology of a hospital hot water system following the introduction of an on-site monochloramine disinfection system. PLoS ONE 2014, 9, e102679. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Williams, M.M.; Santo Domingo, J.W.; Meckes, M.C. Population diversity in model potable water biofilms receiving chlorine or chloramine residual. Biofouling 2005, 21, 279–288. [Google Scholar] [CrossRef] [PubMed]
  55. Flemming, H.-C.; Wingender, J. The biofilm matrix. Nat. Revi. Microbiol. 2010, 8, 623–633. [Google Scholar] [CrossRef] [PubMed]
  56. Faria, S.; Joao, I.; Jordao, L. General Overview on Nontuberculous Mycobacteria, Biofilms, and Human Infection. J. Pathog. 2015, 2015, 809014. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Sutherland, I.W. The biofilm matrix--an immobilized but dynamic microbial environment. Trends Microbiol. 2001, 9, 222–227. [Google Scholar] [CrossRef]
  58. Ojha, A.; Anand, M.; Bhatt, A.; Kremer, L.; Jacobs, W.R., Jr.; Hatfull, G.F. GroEL1: A dedicated chaperone involved in mycolic acid biosynthesis during biofilm formation in mycobacteria. Cell 2005, 123, 861–873. [Google Scholar] [CrossRef] [Green Version]
  59. Ojha, A.K.; Baughn, A.D.; Sambandan, D.; Hsu, T.; Trivelli, X.; Guerardel, Y.; Alahari, A.; Kremer, L.; Jacobs, W.R., Jr.; Hatfull, G.F. Growth of Mycobacterium tuberculosis biofilms containing free mycolic acids and harbouring drug-tolerant bacteria. Mol. Microbiol. 2008, 69, 164–174. [Google Scholar] [CrossRef] [Green Version]
  60. Li, D.; Li, Z.; Yu, J.; Cao, N.; Liu, R.; Yang, M. Characterization of bacterial community structure in a drinking water distribution system during an occurrence of red water. Appl. Environ. Microbiol. 2010, 76, 7171–7180. [Google Scholar] [CrossRef] [Green Version]
  61. Ridgway, H.F.; Means, E.G.; Olson, B.H. Iron bacteria in drinking-water distribution systems: Elemental analysis of Gallionella stalks, using x-ray energy-dispersive microanalysis. Appl. Environ. Microbiol. 1981, 41, 288–297. [Google Scholar] [CrossRef] [Green Version]
  62. Kaestli, M.; O’Donnell, M.; Rose, A.; Webb, J.R.; Mayo, M.; Currie, B.J.; Gibb, K. Opportunistic pathogens and large microbial diversity detected in source-to-distribution drinking water of three remote communities in Northern Australia. PLoS Negl. Trop. Dis. 2019, 13, e0007672. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. He, S.; Tominski, C.; Kappler, A.; Behrens, S.; Roden, E.E. Metagenomic Analyses of the Autotrophic Fe(II)-Oxidizing, Nitrate-Reducing Enrichment Culture KS. Appl. Environ. Microbiol. 2016, 82, 2656–2668. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Li, X.-K.; Chu, Z.-R.; Liu, Y.-J.; Zhu, M.-T.; Yang, L.; Zhang, J. Molecular characterization of microbial populations in full-scale biofilters treating iron, manganese and ammonia containing groundwater in Harbin, China. Bioresour. Technol. 2013, 147, 234–239. [Google Scholar] [CrossRef] [PubMed]
  65. Stanish, L.F.; Hull, N.M.; Robertson, C.E.; Harris, J.K.; Stevens, M.J.; Spear, J.R.; Pace, N.R. Factors Influencing Bacterial Diversity and Community Composition in Municipal Drinking Waters in the Ohio River Basin, USA. PLoS ONE 2016, 11, e0157966. [Google Scholar] [CrossRef]
  66. Pressman, J.G.; Lee, W.H.; Bishop, P.L.; Wahman, D.G. Effect of free ammonia concentration on monochloramine penetration within a nitrifying biofilm and its effect on activity, viability, and recovery. Water Res. 2012, 46, 882–894. [Google Scholar] [CrossRef]
  67. Buffet-Bataillon, S.; Bonnaure-Mallet, M.; de la Pintiere, A.; Defawe, G.; Gautier-Lerestif, A.L.; Fauveau, S.; Minet, J. Heterotrophic bacterial growth on hoses in a neonatal water distribution system. J. Microbiol. Biotechnol. 2010, 20, 779–781. [Google Scholar]
  68. Liu, R.; Yu, Z.; Guo, H.; Liu, M.; Zhang, H.; Yang, M. Pyrosequencing analysis of eukaryotic and bacterial communities in faucet biofilms. Sci. Total Environ. 2012, 435–436, 124–131. [Google Scholar] [CrossRef]
  69. Soto-Giron, M.J.; Rodriguez-R, L.M.; Luo, C.; Elk, M.; Ryu, H.; Hoelle, J.; Santo Domingo, J.W.; Konstantinidis, K.T. Biofilms on Hospital Shower Hoses: Characterization and Implications for Nosocomial Infections. Appl. Environ. Microbiol. 2016, 82, 2872–2883. [Google Scholar] [CrossRef] [Green Version]
  70. Dhandayuthapani, S.; Mudd, M.; Deretic, V. Interactions of OxyR with the promoter region of the oxyR and ahpC genes from Mycobacterium leprae and Mycobacterium tuberculosis. J. Bacteriol. 1997, 179, 2401–2409. [Google Scholar] [CrossRef] [Green Version]
  71. Baron, J.L.; Harris, J.K.; Holinger, E.P.; Duda, S.; Stevens, M.J.; Robertson, C.E.; Ross, K.A.; Pace, N.R.; Stout, J.E. Effect of monochloramine treatment on the microbial ecology of Legionella and associated bacterial populations in a hospital hot water system. Syst. Appl. Microbiol. 2015, 38, 198–205. [Google Scholar] [CrossRef]
  72. Hwang, C.; Ling, F.; Andersen, G.L.; LeChevallier, M.W.; Liu, W.-T. Microbial community dynamics of an urban drinking water distribution system subjected to phases of chloramination and chlorination treatments. Appl. Environ. Microbiol. 2012, 78, 7856–7865. [Google Scholar] [CrossRef] [Green Version]
  73. Yang, Y.; Cheng, D.; Li, Y.; Yu, L.; Gin, K.Y.-H.; Chen, J.P.; Reinhard, M. Effects of monochloramine and hydrogen peroxide on the bacterial community shifts in biologically treated wastewater. Chemosphere 2017, 189, 399–406. [Google Scholar] [CrossRef] [PubMed]
  74. Regione Emilia Romagna. Approvazione delle linee guida regionali per la sorveglianza e il controllo della legionellosi—DGR n. 828/2017. 2017. Available online: http://salute.regione.emilia-romagna.it/documentazione/leggi/regionali/dgr-2127-2016/dgr-n-828-2017-linee-guida-regionali-per-la-sorveglianza-e-il-controllo-della-legionellosi/view (accessed on 1 November 2018).
  75. Wang, H.; Bedard, E.; Prevost, M.; Camper, A.K.; Hill, V.R.; Pruden, A. Methodological approaches for monitoring opportunistic pathogens in premise plumbing: A review. Water Res. 2017, 117, 68–86. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Wang, H.; Edwards, M.A.; Falkinham, J.O., 3rd; Pruden, A. Probiotic approach to pathogen control in premise plumbing systems? A review. Environ. Sci. Technol. 2013, 47, 10117–10128. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Principal components analysis (PCA) of the normalized relative abundance of 16S rRNA sequences in all samples. Data are plotted following the genus-level classification. The PCA demonstrates the clustering of 16S rRNA sequences from cold water (Red) and hot water (Green).
Figure 1. Principal components analysis (PCA) of the normalized relative abundance of 16S rRNA sequences in all samples. Data are plotted following the genus-level classification. The PCA demonstrates the clustering of 16S rRNA sequences from cold water (Red) and hot water (Green).
Ijerph 17 02158 g001
Figure 2. Distribution of genera in inlet cold water pools according to building (cut off 2%). Genera with relative abundance ≤ 2% are listed as “other”.
Figure 2. Distribution of genera in inlet cold water pools according to building (cut off 2%). Genera with relative abundance ≤ 2% are listed as “other”.
Ijerph 17 02158 g002
Figure 3. Distribution of phyla in hot water samples (cut off 1%). Phyla with relative abundance ≤1% are included in "other". HP: network with hydrogen peroxide, M: network with monochloramine, CD: network with chlorine dioxide, and NT: not chemically treated network.
Figure 3. Distribution of phyla in hot water samples (cut off 1%). Phyla with relative abundance ≤1% are included in "other". HP: network with hydrogen peroxide, M: network with monochloramine, CD: network with chlorine dioxide, and NT: not chemically treated network.
Ijerph 17 02158 g003
Figure 4. Distribution of genera in hot water samples (cut off 2%), clustering for hot water network (pools). Genera with relative abundance ≤2% are included in "other". HP: network with hydrogen peroxide, M: network with monochloramine, CD: network with chlorine dioxide, and NT: not chemically treated network.
Figure 4. Distribution of genera in hot water samples (cut off 2%), clustering for hot water network (pools). Genera with relative abundance ≤2% are included in "other". HP: network with hydrogen peroxide, M: network with monochloramine, CD: network with chlorine dioxide, and NT: not chemically treated network.
Ijerph 17 02158 g004
Figure 5. Principal components analysis (PCA) of the normalized relative abundance of 16S rRNA sequences in all samples. Data are plotted following the genus-level classification. Light blue dots for network with hydrogen peroxide (HP), blue for monochloramine (M), red for chlorine dioxide (CD), and green for not chemically treated (NT). Black arrows indicate inlet cold water samples for each network.
Figure 5. Principal components analysis (PCA) of the normalized relative abundance of 16S rRNA sequences in all samples. Data are plotted following the genus-level classification. Light blue dots for network with hydrogen peroxide (HP), blue for monochloramine (M), red for chlorine dioxide (CD), and green for not chemically treated (NT). Black arrows indicate inlet cold water samples for each network.
Ijerph 17 02158 g005
Table 1. Temperature (t), operational taxonomic units (OTUs), Shannon and equitability index for each cold water sample.
Table 1. Temperature (t), operational taxonomic units (OTUs), Shannon and equitability index for each cold water sample.
Inlet Cold Water
MonthSamplet (°C)OTUsShannon IndexEquitability Index
FebruaryBuilding 113.05181.9670.340
Building 212.96642.3050.327
Building 312.93501.7720.363
Building 412.84841.5030.344
meanmedianmeanmean
12.95011.8870.344
MayBuilding 113.82912.5270.267
Building 213.43832.5070.332
Building 312.52801.9190.366
Building 416.73762.0620.393
meanmedianmeanmean
14.13342.2540.340
Table 2. Comparison between hot water according to the networks. HP: network with hydrogen peroxide, M: network with monochloramine, CD: network with chlorine dioxide, and NT: not chemically treated network.
Table 2. Comparison between hot water according to the networks. HP: network with hydrogen peroxide, M: network with monochloramine, CD: network with chlorine dioxide, and NT: not chemically treated network.
Hot Water ParametersNetwork HPNetwork MNetwork CDNetwork NTp-Value
Temperature
Mean (°C)
(Range)
45.2
(33.4–49.2)
40.0
(31.0–46.3)
42.1
(36.2–44.9)
53.2
(50.4–56.6)
<0.0001
Biocide concentration
Mean (ppm)
(Range)
21.3
(18.1–25.0)
3.0
(2.6–3.6)
0.2
(0.1–0.4)
--
OTUs
Median
(Range)
341
(274–414)
384
(252–816)
410
(271–956)
509
(312–981)
0.008
Shannon index
Mean
(Range)
2.316
(0.874–3.549)
2.110
(1.488–2.674)
1.998
(0.926–2.769)
2.135
(1.544–2.532)
0.481
Equitability index
Mean
(Range)
0.395
(0.155–0.629)
0.353
(0.249–0.433)
0.327
(0.165–0.487)
0.343
(0.259–0.425)
0.254
Table 3. The three predominant genera for each hot water network. HP: network with hydrogen peroxide, M: network with monochloramine, CD: network with chlorine dioxide and NT: not chemically treated network.
Table 3. The three predominant genera for each hot water network. HP: network with hydrogen peroxide, M: network with monochloramine, CD: network with chlorine dioxide and NT: not chemically treated network.
NetworkPredominant GeneraRelative Abundance (%)
MeanMinimumMaximum
HPMycobacterium36.361.5984.50
Bradyrhizobium8.582.2322.77
Thiocystis3.450.4311.10
MGallionella17.710.9660.42
Ectothiorhodospira10.420.4427.11
Sediminibacterium4.241.459.87
CDHalorhodospira10.350.0651.93
Blastomonas9.721.4837.28
Bradyrhizobium8.913.3725.82
NTChitinophaga11.704.0421.06
Methylonatrum6.900.0250.20
Thermobaculum6.192.4212.30

Share and Cite

MDPI and ACS Style

Paduano, S.; Marchesi, I.; Casali, M.E.; Valeriani, F.; Frezza, G.; Vecchi, E.; Sircana, L.; Romano Spica, V.; Borella, P.; Bargellini, A. Characterisation of Microbial Community Associated with Different Disinfection Treatments in Hospital hot Water Networks. Int. J. Environ. Res. Public Health 2020, 17, 2158. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17062158

AMA Style

Paduano S, Marchesi I, Casali ME, Valeriani F, Frezza G, Vecchi E, Sircana L, Romano Spica V, Borella P, Bargellini A. Characterisation of Microbial Community Associated with Different Disinfection Treatments in Hospital hot Water Networks. International Journal of Environmental Research and Public Health. 2020; 17(6):2158. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17062158

Chicago/Turabian Style

Paduano, Stefania, Isabella Marchesi, Maria Elisabetta Casali, Federica Valeriani, Giuseppina Frezza, Elena Vecchi, Luca Sircana, Vincenzo Romano Spica, Paola Borella, and Annalisa Bargellini. 2020. "Characterisation of Microbial Community Associated with Different Disinfection Treatments in Hospital hot Water Networks" International Journal of Environmental Research and Public Health 17, no. 6: 2158. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17062158

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop