Prevalence of Opportunistic Pathogens and Diversity of Microbial Communities in the Water System of a Pulmonary Hospital＊

TANG Wei, MAO Yu, LI Qiu Yan, MENG Die, CHEN Ling,3, WANG Hong,3,

ZHU Ren2, and ZHANG Wei Xian1,3,#

1. State Key Laboratory for Pollution Control and Resource Reuse, College of Environmental Science and Engineering, Tongji University, Shanghai 200092, China; 2. Department of Infection Control, Shanghai Pulmonary Hospital, Tongji University, Shanghai 200433, China; 3. Shanghai Institute of Pollution Control and Ecological Security, Shanghai 200092, China

Abstract Objective Our objective was to investigate the occurrence of opportunistic pathogens and characterize the bacterial community structures in the water system of a pulmonary hospital.Methods The water samples were collected from automatic and manual faucets in the consulting room, treatment room, dressing room, respiratory ward, and other non-medical rooms in three buildings of the hospital. Quantitative polymerase chain reaction was used to quantify the load of several waterborne opportunistic pathogens and related microorganisms, including Legionella spp.,Mycobacterium spp., and M. avium. Illumina sequencing targeting 16S rRNA genes was performed to profile bacterial communities.Results The occurrence rates of Legionella spp., Mycobacterium spp., and M. avium were 100%, 100%,and 76%, respectively in all samples. Higher occurrence rates of M. avium were observed in the outpatient service building (building 1, 91.7%) and respiration department and wards (building 2, 80%)than in the office building (building 3), where no M. avium was found. M. avium were more abundant in automatic faucets (average 2.21 × 104 gene copies/L) than in manual faucets (average 1.03 × 104 gene copies/mL) (P ＜ 0.01). Proteobacteria, Actinobacteria, Bacteroidetes, Cyanobacteria, Firmicutes, and Acidobacteria were the dominant bacterial phyla. Disinfectant residuals, nitrate, and temperature were found to be the key environmental factors driving microbial community structure shifts in water systems.Conclusion This study revealed a high level of colonization of water faucets by opportunistic pathogens and provided insight into the characteristics of microbial communities in a hospital water system and approaches to reduce risks of microbial contamination.

Key words: Opportunistic pathogen; Microbial community; Hospital water system; Mycobacterium avium

INTRODUCTION

Drinking-water distribution systems including premise plumbing are complex microbial ecosystems that may harbor a diversity of bacteria, protozoa, virus, and fungi. As the most frequently studied group of microorganisms in water systems, bacteria play a key role in water quality and drinking water safety. In general,Proteobacteria,Bacteroidetes,Cyanobacteria,Verrucomicrobia, andActinobacteriaare the most frequently reported abundant bacterial phyla in the water systems[1]. Notably, the presence of pathogens or opportunistic pathogens in these phyla, such asLegionella pneumophila(L.pneumophila) andMycobacterium avium(M. avium),can pose potential health threats to humans.L.pneumophilaare Gram-negative obligate aerobic bacteria of the phylumProteobacteria, which can cause Legionnaires' disease (LD) or Pontiac fever.Similar to some species within theActinobacteriaphylum, non-tuberculous mycobacteria (NTM) can cause pulmonary diseases resembling tuberculosis,lymphadenitis, soft tissue infection, or disseminated infection in immunocompromised individuals such as acquired immune deficiency syndrome (AIDS) and transplanted patients[2]. So far, although over 150 different species of NTM have been reported, NTM infections are most commonly associated with species such asM. aviumcomplex,M. abscessus,M.fortuitum, andM. kansasii[3].

Owing to their ability to resist stressful environmental conditions (such as oxygen deprivation, nutrition deficiency, temperature shocks, pH changes, and chlorine resistance),Legionellaand NTM are often found in a wide range of natural and artificial aquatic environments[4,5].According to previous studies, these bacteria have been widely detected in municipal water treatment and distribution systems[6,7].Legionellaand NTM colonization in water distribution systems are influenced by many factors, including water temperature, organic matter, pipe material,corrosion, nutrient levels, hydraulic conditions, and biofilm formation in pipe networks[8-10]. The U.S.Environmental Protection Agency (EPA) has includedM. aviumandL. pneumophilain the contaminate candidate list 4 (CCL4), which is a list of contaminants that are currently not subject to any proposed or promulgated national primary drinking water regulations[11].

Many modes have been identified for the transmission of these opportunistic pathogens to humans. Inhalation or micro aspiration of aerosols generated by contaminated water, such as from nebulizers, showers, faucets, fountains, hot water systems, and cooling towers, is the most common route of exposure. The hospital environment is of importance in the cross-transmission of opportunistic pathogens. For example, pathogen contamination on the surfaces of medical devices or other healthcare facilities is a potential exposure route to opportunistic pathogens for susceptible people, which increases the potential risk of healthcare-associated infections. A review of 27 LD outbreaks investigated by the U.S. CDC shown that 33% LD cases are linked to healthcare[12,13].Airborne/aerosols mode of infection are also important migration routes in hospital.Legionella spp. and NTM can survive in water systems and spread in the form of aerosols during showering,bathing, or through faucet water[14]. Like other bacteria, these opportunistic pathogens mostly grow in biofilms, which provide an ecological niche for microorganisms on the surfaces of plumbing materials. Once established within the water system,opportunistic pathogens are difficult to be eradicated.

There is an extensive concern that nosocomial infections are associated with waterborne opportunistic pathogens. Some studies indicated that the water system is one of the main sources for LD outbreaks in hospitals[15,16].Legionellaspp. in hospital drinking water samples have been associated with patient isolates in Italy[17]. The aerators/rectifiers connected to the faucets of the water supply in a Japanese hospital were also suspected to be related to the increased frequency of NTM-positive disease cases[18]. However, there are gaps in the knowledge about the existence of opportunistic pathogens in the water systems of respiratory disease hospitals that hospitalize many patients with pulmonary diseases susceptible to infection with opportunistic pathogens. In addition,little is known about the composition of microbial community, especially that of opportunistic pathogens, associated with bacterial populations in hospital water systems. Considering the complexity of hospital water systems and the vulnerability of hospitalized patients, it is important to conduct studies to survey the presence of opportunistic pathogens and investigate the composition of microbial community in water systems of respiratory disease hospitals.

In this study, representative opportunistic pathogens such asLegionellaspp.,Mycobacteriumspp., andM. aviumwere targeted in a pulmonary hospital water system. This hospital is one of the largest hospitals specializing in respiratory disease treatment and research in China. In addition,amplicon sequencing of 16S rRNA genes was performed. The aims of this study are as follows: (1)to reveal the occurrence of opportunistic pathogens in tap water; (2) to evaluate the correlations of water quality parameters like temperature, pH, total chlorine, nitrate, nitrite, and ammonia with opportunistic pathogens; (3) to characterize microbial communities and to explore the influence of drinking water systems on the drinking water microbiome. The information from this study might be useful for improving hospital water systems to prevent nosocomial infections caused by opportunistic pathogens.

MATERIALS AND METHODS

Sample Collection and Processing

The work was conducted in a pulmonary hospital in Shanghai, China, which is one of the largest hospitals specialized in respiratory disease treatment, research and education in China. This hospital was founded in 1933, covering an area of 103,000 square meters. There are 1,348 employees,13 clinical departments (such as thoracic surgery,oncology department, respiration department,tuberculosis department, occupational disease department), and 11 medical technique departments. It reports 1,044,704 outpatient visits and discharges 88,768 patients on a yearly average.

Drinking water samples were collected from three buildings during March to May, 2018 in this hospital, where chloramine was used as the secondary disinfectant. The first building is used for outpatient service, the second building are wards of respiration department, and the third building is an office building free of medical services, serving as a control site in this study. The detailed information for sampling is shown in Supplementary Table S1(available in www.besjournal.com). Samples(1,000 mL) were collected after flushing for 3 min,and then the bulk water sample was filtered through a sterile 0.22 μm nitrocellulose membrane(Millipore, MA) using a peristaltic pump. A filter membrane with collected biomass was transferred into a 2 mL tube (the lysing matrix A tube) from FastDNA SPIN Kit (MP Biomedicals, OH). The sample was further processed according to the manufacturer’s protocol. The extracted DNA was purified, precipitated, and dissolved in 50 μL of sterile nuclease-free water. The concentration and purity of the total DNA were measured using NanoDrop 2000 (Thermo Fisher Scientific).

Water Quality Analysis

Temperature, pH, total chlorine, nitrate, nitrite,and ammonia were measured at the time of collection. pH and temperature were monitored using a portable pH meter (Bante Instrument,Bante903P). Total chlorine, nitrate, nitrite, and ammonia were analyzed using HACH reagents, on a DR900 colorimeter (HACH, Loveland, CO).

Quantitative PCR (qPCR)

Legionellaspp.,Mycobacteriumspp.,M. avium,and total bacterial 16S rRNA genes were quantified by qPCR assays using previously published primers and probes[19]. All qPCR reactions were performed in triplicate in a 10 μL system containing TaqMan Probes or SybrGreen Premix (TaKaRa). For SYBR Green assays, the 10 μL reaction mixture contained 5 μL of SYBR Green Premix, 0.8 μL each of 5 μmol/L forward and reverse primers, 0.16 μL ROX, 2.24 μL DNase/RNase-free water, and 1 μL template DNA.For TaqMan assays, the 10 μL reaction mixture contained 5 μL of Probes Premix, 0.4 μL each of 5 μmol/L primers, 0.35 μL of 10 μmol/L probe, 0.1 μL ROX, 2.75 μL DNase/RNase-free water, and 1 μL template DNA. The qPCR conditions were as follows:30 s at 95 °C, 40 cycles of 5 s at 95 °C, 34 s at 55 °C,and 40 s at 72 °C (for 16S rRNA genes); 1 min at 95 °C, 40 cycles of 5 s at 95 °C, 34 s at 64 °C, and 40 s at 72 °C (forM. avium); 30 s at 95 °C, 40 cycles of 5 s at 95 °C, 34 s at 58.5 °C, and 40 s at 72 °C (forLegionellaspp.); 30 s at 95 °C, 45 cycles of 5 s at 95 °C, 34 s at 55 °C, and 40 s at 72 °C (forMycobacteriumspp.). Samples were diluted at 1:10 ratios to minimize PCR inhibition. Standards (positive controls) and sterile water (negative controls) were included in each qPCR run. qPCR standards were analyzed to determine the amplification efficiencies(range: between 90% and 110%) and the correlation coefficient (R2＞ 0.99 for all the standard curves).Samples yielding a detectable threshold cycle in at least two experiments out of three were scored as positive.

Illumina MiSeq Amplicon Sequencing

The V3-V4 region of 16S rRNA genes were amplified for all DNA extracts with primers 338F(forward primer sequence: 5’-ACTCCTACGGG AGGCAGCA-3’) and 806R (reverse primer sequence:5’-GGACTACHVGGGTWTCTAAT-3’). PCRs were performed with a 20 μL system using TransStart Fastpfu DNA Polymerase (Beijing TransGen Biotech).The PCR products from triplicate amplifications were examined using a 2% agarose gel, and then purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, USA) according to the manufacturer’s instructions. Purified amplicons were pooled in equimolar amounts prior to amplicons sequencing on an Illumina MiSeq platform using paired-end 250 bp kits according to the standard protocols at the Majorbio Bio-Pharm Technology Corporation in Shanghai, China.

Paired-end reads were stitched and screened using Trimmomatic to remove short sequences,singletons, and noisy reads, prior to alignment using FLASH[20,21]. Chimeras were identified and removed from further analysis using UCHIME[22]. Sequences were clustered into operational taxonomic units(OTUs) using the average neighbor algorithm with a similarity cut-off of 97%[23]. The taxonomic information of sequences was aligned against the Greengenes database, with a threshold confidence level of 70%[24].

The software Mothur was used to analyze the alpha diversity indices, including the coverage,observed OTUs richness (Sobs), Shannon index,abundance-based coverage estimator (ACE) and Chao1[25]. The intersection of OTUs lists for different samples was analyzed using the vegan packages in R to reveal overlapped and unique OTUs in different buildings. Non-metric multidimensional scaling(NMDS) analysis and analysis of similarity (ANOSIM)were performed based on generated Bray-Curtis distance matrices to evaluate the variation of the microbial community structure among different sample groups. Redundancy analysis (RDA) was used to evaluate the relationship between the abundance of microbial species and environmental variables,using the vegan and ggplot2 packages in R v3.5.0 software.

Statistical Analysis

qPCR and physiochemical data were analyzed with a significance level ＜ 0.05. The Shapiro-Wilk test was used to test the normality of data. Oneway analysis of variance (ANOVA) or Kruskal-WallisHtest was performed to compare target microbes in different buildings depending on the normality and variance of the data. Wilcoxon rank-sum test was used to compare physiochemical and/or biological data between the two different types of faucets. Spearman correlation analysis was used to identify associations between qPCR data and physiochemical parameters using the IBM SPSS statistics v24.0 software. Maximal Informationbased Nonparametric Exploration (MINE) analysis,a tool for finding the pairwise associations (linear or non-linear association) in a dataset, with a false detection rate correction at 0.05 was used to characterize the associations betweenLegionellaspp.,Mycobacteriumspp., and other bacterial genera[26,27].

RESULTS

Water Quality Characteristics in the Hospital Water System

The water quality parameters in the hospital water system are shown in Table 1 and Supplementary Table S2, available in www.besjournal.com. The temperature, pH, total chlorine,NO3-N, NO2-N, and NH4-N varied in the range of 12.7-27.5 °C, 7.17-7.91, 0.06-0.93 mg/L, 0.3-2.9 mg/L, ＜ 0.001-0.046 mg/L, and 0.07-0.20 mg/L among all water samples, respectively. The levels of total chlorine in building 2 were significantly higher than in building 1 and building 3 (P＜ 0.05). NO3-N shared a similar trend: building 3 ＜ building 1 ＜building 2 (P＜ 0.05). However, the variation of temperature was opposite, with the average temperature in building 3 about 7.0 °C higher than that of building 2 (P＜ 0.05). The concentrations of NH4-N in building 1 and building 2 were approximately two times higher than building 3 (P＜0.05). No significant difference was observed for NO2-N among the three buildings.

Table 1. Physiochemical parameters of water samples collected from the hospital water system

Occurrence of Total Bacteria, Legionella spp.,Mycobacterium spp., and Its Species M. avium

The frequencies of detection (FOD) ofLegionellaspp.,Mycobacteriumspp., andM.aviumin the three hospital buildings were 100%,100%, and 76%, respectively (Figure 1 and Supplementary Table S3, available in www.besjournal.com). The total bacteria,Legionellaspp., andMycobacteriumspp. were detected in all water samples taken from 3 different buildings,ranging from (0.70-8.03) × 106, (0.11-1.30) × 105,and (0.56-12.09) × 104gene copies/L, respectively.M. aviumwas only detected in water samples of building 1 and building 2, while noM. aviumwas found in building 3. The FOD ofM. avium(91.7%)in building 1 was higher than that in building 2(80.0%), and the FOD in all samples was 76.0%(19 out of 25). There was a significant difference in the number ofM. aviumgene copies between the automatic faucet and manual faucet (P＜ 0.01).Different levels of total bacteria andLegionellagene copies were observed among the three buildings (P＜ 0.05).

Correlations between Microorganisms and Water Quality Parameters

The correlations between the water quality parameters and the gene copies of total bacteria and opportunistic pathogens are shown in Table 2. In all samples, significant negative correlations were found between the number of total bacterial 16S rRNA gene copies and total chlorine (r= -0.650,P＜0.01), and between the number ofM. aviumgene copies and total chlorine (r= -0.705,P＜ 0.01).Positive correlations were observed between the number of total bacterial 16S rRNA gene copies and temperature (r= 0.508,P＜ 0.01), pH (r= 0.535,P＜0.01), and nitrite (r= 0.475,P＜ 0.05). Correlation analysis was also performed between target organisms. The number of gene copies ofMycobacterium(r= 0.554),M. avium(r= 0.689)showed positive correlations with 16S rRNA gene numbers, respectively (P＜ 0.01).

Microbial Community Composition in the Hospital Water System

Illumina sequencing of 16S rRNA gene amplicons was performed to characterize microbial community compositions in the hospital water system (Figure 2).A total of 1,427,942 sequences were obtained from 25 water samples. These sequences were clustered into 3,182 OTUs at a 97% similarity level. The relative abundances of the most dominant phyla across the 25 samples were Proteobacteria (26.20%-90.22%),Actinobacteria (0.84%-42.63%), Bacteroidetes(0.85%-29.33%), Cyanobacteria (1.30%-14.24%),Firmicutes (0.04%-26.34%), Acidobacteria (0.32%-2.38%), Planctomycetes (0.15%-3.58%), Nitrospirae(0.04%-5.34%), Chlamydiae (0.04%-1.12%), and Gemmatimonadetes (0.05%-1.23%). At the class level, Alphaproteobacteria, Betaproteobacteria,Actinobacteria, Gammaproteobacteria, Chloroplast,Clostridia, Saprospirae, Acidimicrobiia, Cytophagia,Deltaproteobacteria, Flavobacteriia, Nitrospira, and Bacilli were the dominant classes in these samples(the sum of relative abundance ＞ 1%).Novosphingobium, Mycobacterium, and Sphingomonas were the most prevailing genera in all the water samples, the combined relative abundances of which were 11.21%, 8.36%, and 7.17%, respectively. Other major genera in the samples were Schlegelella, Pseudomonas,Sediminibacterium, Hyphomicrobium,Bradyrhizobium, and Nitrospira (the relative abundance ＞ 1%). MINE analysis identified 22 statistically significant associations betweenMycobacteriumspp. and other microbial genera under a false detection rate of 0.05 (Supplementary Table S4 available in www.besjournal.com); whereas no significant association was found betweenLegionellaspp. and other microbial genera. The highest maximal information coefficient (MIC)betweenMycobacterium spp. and other microbial genera was withMethylosinus(MIC = 0.9427, TIC =0.7910), followed bySediminibacterium(MIC =0.8411, TIC = 0.3093),Polaromonas(MIC = 0.8325,TIC = 0.2773) andLimnohabitans(MIC = 0.7983,TIC = 0.4796).

The main differences between the microbial community in phylum and class (top 20) in different groups are shown in Supplementary Tables S5-S6,available in www.besjournal.com. For the three buildings, the significant differences in abundance(top 20) were observed in Verrucomicrobia (P＜0.01), Gemmatimonadetes (P＜ 0.01), Synergistetes(P＜ 0.01) at the phylum level, and in Deltaproteobacteria (P＜ 0.01), Planctomycetia (P＜0.01), Holophagae (P＜ 0.05) and Sphingobacteriia(P＜ 0.01) at the class level. For the two types of faucets, the relative abundances of microorganisms(top 10) were different. The relative abundances of Firmicutes (phylum,P＜ 0.01), Verrucomicrobia(phylum,P＜ 0.01), Clostridia (class,P＜ 0.01), and Gammaproteobacteria (class,P＜ 0.05) in samples taken from manual faucets were higher than thosein samples taken from automatic faucets.

Table 2. Correlations between the water quality parameters and the gene copies of total bacteria and opportunistic pathogens

To evaluate the shared OTUs in different building taps in this hospital, Venn diagrams were constructed (Supplementary Figure S1, available in www.besjournal.com). There were 122, 347, and 271 shared OTUs present in water samples collected from building 1, building 2, and building 3,respectively. Five hundred and forty shared OTUs were found in the three buildings, with the 10 most abundant genera identified asNovosphingobium,Mycobacterium,Sphingomonas,Pseudomonas,Sediminibacterium,Hyphomicrobium,Bradyrhizobium,Nitrospira,Sphingopyxis, andErythromicrobium.

Microbial community richness and diversity indices in all the water samples are shown in Supplementary Table S7, available in www.besjournal.com. The coverage of all samples was above 99.7%, showing that sequencing depth was sufficient in this work. There were significant differences in the Sobs, ACE and Chao1 indexes in three buildings (P＜ 0.05), indicating the variation of the richness and diversity of the microbial community in different buildings.

NMDS analysis of the samples based on Bray-Curtis distance matrices is illustrated in Figure 3. The result showed that the samples were clustered into three groups by different buildings. However, the samples of washroom and drinking room water collected from building 1 clustered closer to those of building 2, suggesting higher community similarity to building 2 water samples than the rest in building 1.ANOSIM analysis also demonstrated the difference in microbial community structures between different buildings (r= 0.44,P= 0.001).

Relationships between Water Quality Properties and Microbial Community Composition

The relationships between water quality properties and microbial community composition were assessed with RDA (Figure 4). RDA1 axis showed a significant correlation with temperature,chlorine, and nitrate (P＜ 0.05), and the RDA2 axis showed a significant correlation with temperature (P＜ 0.05). Over 41.6% of the variation in bacterial community composition can be explained by the RDA1 axis and RDA2 axis. Similar to NMDS (Figure 3),the RDA results showed that the samples from building 3 tightly clustered into a group, and the samples from building 2 had the similar trend except that two samples were taken from the respiratory ward (i.e., points located in the upper of the figure);however, the samples collected from building 1 dispersed without obvious clustering. Total chlorine,nitrate, and temperature were the major environmental driver factors, showing a significant correlation with the bacterial community composition (P＜ 0.05).

DISCUSSION

This research revealed a high prevalence ofLegionellaspp. andMycobacteriumspp. in water faucets and provided insight into these opportunistic pathogens and the broader bacterial community composition in a respiratory disease hospital water system, which may have implications for reducing microbial contamination risks of hospital water supplies.

Legionellaspp. andMycobacteriumspp. were ubiquitously found in water supply systems, such as drinking water treatment trains, premise plumbing,water reservoirs, water storage tanks, water heaters and showerheads[28-33]. Because of the complexity of hospital water systems and the vulnerability of hospitalized patients, many of the infections associated with opportunistic pathogens have been reported in hospital environments[34]. However, to our knowledge, there have been few reports on a large-scale specialized respiratory system hospital in China using culture-independent methods.

In this study,Legionellaspp. andMycobacteriumspp. were ubiquitously present (100%) in the hospital water system, whileM. aviumwas only found in buildings for medical services. The finding was similar to one previous study on secondary water supply systems in Shanghai, demonstrating the ubiquity ofLegionellaand mycobacterial genera and the absence ofM. aviumin public and residential buildings[35]. Another one-year survey of opportunistic pathogens in the tap-water (including hospital, hotel, school, public office, and private dining room) of one northern city of China also showed high detection rates ofLegionellaspp. andMycobacteriumspp. (i.e., 100% FOD in 44 samples)by qPCR[36]. However, lower FOD ofLegionellawas also reported in private hospitals and university hospitals in Italy, which used chlorine dioxide as disinfection rather than chloramine[37,38].

NTM are environmental mycobacteria other thanMycobacterium tuberculosisandMycobacterium leprae. Species such asM. avium,M. intracellulare,M. mucogenicum,M. kansasii,M. gordonae,M.fortuitum, andM. chelonaare frequently reported opportunistic pathogens for humans. In this study,the high FOD ofM. aviumin tap water was in contrast to another study isolating NTM from patients in the same hospital during 2008-2012,demonstratingM. kansasii(45.0%),M. intracellulare(20.8%),M. chelonae/abscessus(14.9%) as the most frequently identified clones, withM. aviumonly accounting for 3.6% of clinical isolates[39]. Although we did not investigate other mycobacterial species in water samples and biofilm samples in the present study, it is likely that the existence of biofilms in the water systems can harbor a high level of opportunistic pathogens and potentially release them into the water. Different analytic methods employed in these two studies might also account for the differences. Our work applied molecular techniques with advantages of lower detection limits and higher specificity, while the study employed conventional culture methods, which was considered as the gold standard for identifying live pathogens[39]. It is also possible that the composition ofMycobacteriumspp. in this hospital may demonstrate temporal changes.

Ecological interactions among bacteria are ubiquitous and complex in water, the pattern of which is determined by the bacterial population and water quality parameters[40]. As shown in Supplementary Table S4, MINE analysis identified 22 statistically significant associations betweenMycobacteriumspp. and other microbial genera,such asMethylosinus,Sediminibacterium,Novosphingobium,Bdellovibrio,Staphylococcus,Flavobacterium, andRalstonia(P＜ 0.05). Similar to our previous study, the associations betweenMycobacteriumspp. and four bacteria genera includingNovosphingobium,Flavobacterium,

Ralstonia, andArcobacterwere found, suggesting their antagonistic or synergistic relationships in water systems[41]. In addition,Flavobacteriumcan produce some metabolites and components, such as dimethyl polysulfides, which result in the odor in tap water[42]. Notably, a positive significant association was identified betweenMycobacteriumspp. and pathogenic staphylococcus, which may increase virulence potential in water systems. On the contrary, a negative significant association was found betweenMycobacteriumspp. andBdellovibrio.Bdellovibriois a parasite or predator of other bacteria including some pathogens, which may act as natural inhibitors of bacterial growth in manmade water systems[43].

Chloramine has replaced free chlorine as a disinfectant in municipal water systems in shanghai,because it provides many benefits, including higher stability, lower production trihalomethanes, and haloacetic acids, and improves maintenance of disinfectant residual in distribution pipes[44,45].Chloramine has been known to penetrate biofilm more strongly compared to chlorine, while it less effectively inactivated microorganisms near the biofilm surface[46]. It is reported that mycobacterial FODs in a chloraminated water potable system was higher than other chlorinated systems[47]. Donohue et al.[48]also found that the levels of NTM in chloraminated water were statistically higher than those in chlorinated water.M. aviumwas reported as one of the most frequently detected NTM in many water systems. Mycobacteria can be categorized into fast-growing species (colony formation 3-7 d,includingM. abscessus,M. chelonae, andM.fortuitum) and slow-growing species (colony formation ＞ 7 d, includingM. avium,M.intracellulare,M. xenopi, andM. kansasii)[49,50]. As a slow-growing mycobacterium,M. aviumwas more disinfectant-resistant than fast-growing species[51].

According to another survey by Coniglio et al.[52],monochloramine was associated with a lower prevalence ofLegionellacolonization in complex piping systems. However, in our study, the occurrence ofLegionellawas 100%. HigherLegionellaprevalence is likely associated with relatively low disinfectant levels in the present work.About 30% (7/25, 28.0%) of samples had total chorine concentrations lower than 0.5 mg/L, which is the disinfectant concentration recommended for the prevention ofLegionellosisby the American Society of Heating, Refrigeration, and Air Conditioning Engineering[53]. In addition, there are many other factors influencing the growth and survival ofLegionellaspecies. For example, the free-living amoeba can serve as hosts forLegionellaspp.,providing an intracellular environment to protect from disinfectants in water systems[54].

Our study showed that the number ofM. aviumgene copies in water samples collected from the automatic faucet was higher than that collected from the manual faucet (P＜ 0.01). Faucets were likely to be an important source of healthcareassociated infections. One outbreak ofM.mucogenicumbacteremia was reported to be associated with contaminated electronic faucets in a pediatric hematology-oncology department[55].Studies also demonstrated thatLegionellaspecies andPseudomonas aeruginosahad a similar trend with higher FODs in automatic/sensor (notouch)/electronic faucets[56,57]. Some possible reasons included: low water flow, materials of electronic faucets internal components, bacterial colonization of internal magnetic valves and retrograde contamination from the faucet outlet[58].

Hospital water systems are complex ecosystems that colonize trillions of microorganisms interacting with each other and with their environment. Nitrate was one of the major drivers in microbial composition across all samples in this study(Figure 4). High relative abundances ofNitrospirahas also been noted in our study. It has been reported that nitrification can occur in a chloraminated system, which can promote the growth of nitrifying bacteria[59]. Nitrification can also reduce water quality, leading to inadequate disinfectant residual and shift of microbial community structures[28].

Proteobacteria, one of the bacterial phyla with highest abundance in environmental water samples,was the most abundant phylum of gram-negative bacteria, mainly including five classes, viz Alphaproteobacteria, Betaproteobacteria,Gammaproteobacteria, Deltaproteobacteria, and Epsilonproteobacteria[60-62]. Proteobacteria widely existed in water systems, which is consistent with previous studies[35,63]. For example, Vaz-Moreira et al.[63]reported that Alphaproteobacteria predominated in tap water, while Gamma dominated in source water. The ubiquitous tap water genera comprised Proteobacteria of the class Alpha- (Blastomonas, Brevundimonas,Methylobacterium, Sphingobium, Sphingomonas),Beta- (Acidovorax, Ralstonia), and Gamma-(Acinetobacter and Pseudomonas). At phylum level,the relative abundances of Verrucomicrobia (P＜0.01) and Gemmatimonadetes (P＜ 0.01) were noted among the three hospital buildings. Verrucomicrobia and Acidobacteria classified as ‘sensitive taxa’ could act as biomarkers to distinguish samples between source and tap water according to a study by Han et al.[1], which assessed the impact of source water on tap water bacterial communities in 46 drinking water supply systems in China.

Several bacterial genera with a relatively high abundance were known to contain pathogenic species.Mycobacteriumwas (8.36%) the second most prevailing genus in all the water samples.Mycobacteriumcontained a series of pathogenic species, such asM. abscessus,M. chelonae,M.avium,M. gordonae, andM. kansasii. According to some studies, NTM infection in hospitals should consider the possibility of patient colonization or environmental contamination. In particular,M.gordonaederived from contamination of water supply systems or other environmental sources have been associated with the pseudo-outbreaks of hospital infection[18,64]. In thePseudomonasgenus,some species, such asP. aeruginosa, can cause serious infections in immunocompromised patients with cancer, severe burns, and cystic fibrosis[65]. 30%to 50% of nosocomial infections byP. aeruginosawere caused by contaminated water in ICU wards,and the types of faucets and connecting pipes materials are important factors influencing the prevalence ofP. aeruginosa[66].

CONCLUSION

Our results showed a ubiquitous presence of

Legionellaspp.,Mycobacteriumspp., andM. aviumin the pulmonary hospital drinking water system,implying the hospital water system is an important reservoir for opportunistic pathogens. High-throughput sequencing indicated that a shared community of bacteria, mainly includingNovosphingobium,

Mycobacterium, andSphingomonas, was present constantly in the three hospital buildings. The relative abundances of Firmicutes (phylum), Verrucomicrobia(phylum), Clostridia (class), and Gammaproteobacteria(class) in samples taken from manual faucets were higher than that those of samples taken from automatic faucets. Disinfectant residuals, nitrate, and temperature could significantly affect opportunistic pathogens and microbial diversity in a hospital water system. Our study also showed different pathogen contamination levels and variations of microbial community structures in different hospital buildings,implying their significant role in shaping the drinking water microbiome and determining the pathogen occurrence rates. Further research needs to focus on the survival mechanisms of opportunistic pathogens,such asM. avium, in hospital water systems and approaches to reduce risks of microbial contamination.

CONFLICTS OF INTEREST STATEMENT

The authors declare that they have no competing interests.

Received: June 5, 2019;

Accepted: February 20, 2020