Centralized Drinking Water Treatment Operations ... - ACS Publications

control for microbial community structure in finished drinking water. .... site. The treatment process and the sampling locations are outlined in Figu...
0 downloads 4 Views 640KB Size
Subscriber access provided by Binghamton University | Libraries

Article

Centralized Drinking Water Treatment Operations Shape Bacterial and Fungal Community Structure Xiao Ma, Amit Vikram, Leonard W. Casson, and Kyle Bibby Environ. Sci. Technol., Just Accepted Manuscript • Publication Date (Web): 31 May 2017 Downloaded from http://pubs.acs.org on June 1, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 31

Environmental Science & Technology

1 2 3

Centralized Drinking Water Treatment Operations Shape Bacterial and Fungal Community Structure

4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39

Xiao Ma1, Amit Vikram1, Leonard Casson1,2, Kyle Bibby1,3*

1

Department of Civil and Environmental Engineering, University of Pittsburgh, Pittsburgh, PA 15261 2 Graduate School of Public Health, University of Pittsburgh 3 Department of Computational and Systems Biology, University of Pittsburgh Medical School, Pittsburgh, PA 15261, USA

*Corresponding Author: Kyle Bibby, 709 Benedum Hall, Pittsburgh, PA 15261 [email protected], 412-624-9207

Keywords Fungi, bacteria, drinking water, microbiome, filtration, disinfection

1 ACS Paragon Plus Environment

Environmental Science & Technology

Page 2 of 31

40

Abstract

41

Drinking water microbial communities impact opportunistic pathogen colonization and corrosion

42

of water distribution systems, and centralized drinking water treatment represents a potential

43

control for microbial community structure in finished drinking water. In this manuscript, we

44

examine bacterial and fungal abundance and diversity, as well as microbial community

45

taxonomic structure following each unit operation in a conventional surface water treatment

46

plant. Treatment operations drove the microbial composition more strongly than sampling time.

47

Both bacterial and fungal abundance and diversity decreased following sedimentation and

48

filtration; however, only bacterial abundance and diversity was significantly impacted by free

49

chlorine disinfection. Similarly, each treatment step was found to shift bacterial and fungal

50

community beta-diversity, with the exception of disinfection on fungal community structure. We

51

observed the enrichment of bacterial and fungal taxa commonly found in drinking water

52

distribution systems through the treatment process, e.g. Sphingomonas following filtration and

53

Leptospirillium and Penicillium following disinfection. Study results suggest that centralized

54

drinking water treatment processes shape the final drinking water microbial community via

55

selection of community members, and that the bacterial community is primarily driven by

56

disinfection while the eukaryotic community is primarily controlled by physical treatment

57

processes.

58 59 60 61 62

2 ACS Paragon Plus Environment

Page 3 of 31

Environmental Science & Technology

63

1. Introduction

64

Centralized drinking water treatment processes are designed as barriers to remove primary

65

pathogens from raw water and have greatly reduced the burden of waterborne infectious

66

disease.1,

67

treated drinking water is far from a sterile environment and hosts diverse microbial

68

communities.1, 3 Previous studies have estimated bacterial abundance in finished drinking water

69

to range from 103 to 105 cells/mL.3-5 Microbial eukaryotes, such as fungi and free-living amoeba,

70

have also frequently been isolated from drinking water.6-8 Within this post-treatment indigenous

71

drinking water microbiome, opportunistic pathogens have emerged as a significant public health

72

issue, posing the greatest infectious disease risk associated with drinking water in industrialized

73

countries.3,

74

negative impacts, such as nitrification and biocorrosion.1, 3

2

Although drinking water treatment processes greatly reduce microbial abundance,

9

Members of the drinking water microbiome also have potential to cause other

75 76

Complete sterilization of drinking water is infeasible. On the other hand, drinking water

77

treatment processes may selectively remove members of the microbial community,10,

78

providing a potentially beneficial opportunity to manage the drinking water microbial

79

community structure.3, 9 Previous studies of the microbial ecology of drinking water treatment

80

processes have primarily investigated the bacterial community structure.10-12 A notable previous

81

study demonstrated that bacteria colonizing filter media governed the downstream bacterial

82

community structure, suggesting this process to be a primary driver of downstream microbial

83

community composition.10 Additionally, although the drinking water microbiome also contains

84

diverse fungi and potential opportunistic pathogens such as Aspergillus spp. and Fusarium

85

spp.,13-15 fungal community dynamics within drinking water treatment processes have not been

11

3 ACS Paragon Plus Environment

Environmental Science & Technology

Page 4 of 31

86

holistically evaluated. Examining the fungal community through drinking water treatment

87

provides an opportunity to discern the role of treatment processes in shaping the composition of

88

microbial eukaryotes. A knowledge gap remains regarding the dynamics of bacterial and

89

eukaryote ecology along the drinking water treatment process and operational factors that shape

90

the community structure.3

91 92

In the present study we evaluated bacterial and fungal diversity, taxonomy, and abundance along

93

a conventional surface water treatment process train. We used high-throughput sequencing of

94

bacterial 16S rRNA and fungal ITS1 gene regions to characterize the microbial community

95

diversity and taxonomy. Additionally, we evaluated bacterial and fungal abundance, as well as

96

presence of Acanthamoeba spp. and Aspergillus fumigatus. The goal of this study was to

97

investigate the influence of drinking water treatment unit operations, including media filtration

98

and disinfection, on both bacterial and fungal community structure. An enhanced understanding

99

of the role of centralized drinking water treatment in driving microbial community composition

100

will contribute to efforts to engineer a desirable drinking water microbiome.

101 102

2. Materials and Methods

103

2.1. Sampling Site

104

An anonymous U.S. drinking water treatment plant fed with river water was selected as the study

105

site. The treatment process and the sampling locations are outlined in Figure S1. River water is

106

pumped into the plant as raw water, and ferric chloride and cationic polymer (approximately 25

107

mg/L and 1 mg/L, respectively) are added for coagulation. The water passes through a clarifier

108

for initial clarification, and is settled for an additional 24 hours in a secondary sedimentation

4 ACS Paragon Plus Environment

Page 5 of 31

Environmental Science & Technology

109

basin. The settled water passes through a media filter composed of coal, sand, and support

110

gravel. Free chlorine disinfection using sodium hypochlorite is the last step of treatment, and a

111

free chlorine residual of 0.4 to 1.0 mg/L was maintained during the study period.

112 113

2.2. Sample Collection

114

Bulk water samples were collected from four locations along the treatment process (Figure S1):

115

(1) raw water at the pump to the treatment plant; (2) post-sedimentation water, the effluent from

116

the secondary 24-hour sedimentation basin; (3) post-filtration water, the effluent from media

117

filtration; (4) post-disinfection water, the effluent from the disinfection clear well. Eight

118

duplicate monthly water samples were collected from these four locations on January, February,

119

March, April, May, June, July, and September 2014. Samples were collected in duplicate using a

120

1-L sterile Nalgene polypropylene sampling bottle without preservative.

121 122

Water samples were immediately transported on ice to the laboratory and processed within four

123

hours of collection. Upon arrival, 500 mL to 1900 mL of each water sample was filtered through

124

a 0.2 µm Supor® 200 Polyethersulfone membrane (Pall Corporation) housed in sterile Nalgene

125

analytical filter funnels (Thermo Scientific; Fisher) within a class II biological safety cabinet.

126

The filter membrane with biomass was then stored in a sterile 10 mL tube at -80◦C until DNA

127

extraction. The remaining subset of the water samples were measured for conductivity

128

(Accument conductivity meter, Fisher Scientific), pH (Russell RL 060P pH meter, Thermo

129

Scientific), and free chlorine (Hach colorimeter, Hach) following manufacturer’s protocols. Two

130

additional weekly sample sets from the described sampling locations were collected during June

131

27th to July 8th 2016. These additional samples were analyzed for heterotrophic plate count (HPC)

5 ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 31

132

and total filamentous fungal colony count to evaluate the abundance of culturable bacteria and

133

fungi through the studied water treatment processes train. Detailed methods for HPC and total

134

filamentous fungal colony count analysis are provided in the Supplementary Information.

135 136

2.3. DNA Extraction, PCR, and Sequencing

137

Frozen filter membranes with biomass were thawed at 4ºC and subjected to DNA extraction

138

using a RapidWater DNA Isolation Kit (MO BIO Laboratories) following the manufacturer’s

139

protocol. All duplicate samples were subjected to PCR targeting the fungal ITS1 gene region

140

using primers ITS1F and ITS2 modified for Illumina multiplex sequencing as previously

141

described.16 One replicate for each sample was subjected to 16S rRNA PCR using previously

142

described primers 515F and 806R.17 Detailed description of PCR conditions for ITS1 and 16S

143

rRNA PCR can be found in the Supplementary Information. DNA extraction negative controls

144

with a blank filter membrane were both measured by Qubit 2.0 DS DNA high sensitivity kit and

145

included in PCR for sequencing library preparation and no-template PCR negative controls were

146

included in all PCR runs to ensure reagents and equipment were not contaminated. All controls

147

were negative for contamination. Two pure culture fungal DNA samples, Penicillium

148

chrysogenum (CAES PC-1) and Aspergillus fumigatus (ATCC 34506), were included in the

149

ITS1 PCR as positive controls.

150 151

PCR products were purified with Agencourt AMPure XP magnetic beads (Beckman Coulter) and

152

confirmed by gel visualization. The purified PCR products were quantified with a Qubit 2.0

153

Fluorometer with a Qubit dsDNA HS Assay Kit (Invitrogen), and equimolar amounts were

154

pooled to construct fungal ITS amplicon and bacterial 16S rRNA amplicon libraries. One

6 ACS Paragon Plus Environment

Page 7 of 31

Environmental Science & Technology

155

additional purification step was performed for each pooled library with Agencourt AMPure XP

156

magnetic beads. The ITS and 16S rRNA libraries were then sequenced on an Illumina MiSeq

157

Sequencer (Illumina) in two separate runs using MiSeq Reagent Kit v2 (300-cycles) (Illumina).

158 159

2.4. Bioinformatics and Statistical Analyses

160

Raw 16S rRNA sequence data was demultiplexed and trimmed based on a quality score of 20

161

using QIIME 1.8.0.18 After trimming, sequences were grouped into operational taxonomic units

162

(OTUs) at 97% similarity and taxonomically classified using a closed reference OTU picking

163

strategy and Greengenes 13.5 as the reference database;19 17.6% of sequences did not match the

164

reference database were not included in downstream analysis (details in Table S1), and

165

sequences classified as chloroplast were removed from further analysis. As a previous study has

166

demonstrated that grouping fungal ITS sequencing reads into OTUs limits the taxonomic

167

classification accuracy for the fungal community,20 raw sequencing data of ITS1 library was

168

trimmed and demultiplexed based on a quality score of 20 using QIIME 1.8.018 and compared

169

directly against the reference ITS database21 using BLAST ver. 2.2.19.22 The BLASTn output

170

was sorted using FHiTINGs.23 Boxplots of relative abundances were generated using R packages

171

reshape2 and ggplot2 with default options; the outlier was defined as less than the first quartile -

172

1.5*interquartile range or greater than the third quartile + 1.5* interquartile range.24, 25 A two

173

sample t-test was conducted to compare mean relative abundance value of bacterial and fungal

174

taxa at each sampling site. Additional analyses using LEfSe (Linear discriminant analysis effect

175

size) tool was conducted to confirm enriched taxa at each sampling location (Kruskal-Wallis test

176

p ≤ 0.05; log10 LDA score threshold equals 2.0).26 The sequencing results of the two positive

177

ITS1 amplicon controls indicated accurate classification of fungal genera, as 99.9% of sequences

7 ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 31

178

of Aspergillus fumigatus were assigned to Aspergillus spp. and 99.9% of sequences of

179

Penicillium chrysogenum were accurately assigned to Penicillium spp. These results are

180

consistent with previous observations of ITS1 classification annotation accuracy to the genus

181

level.27

182 183

Unweighted UniFrac (membership-based) and weighted UniFrac (structure-based) matrices were

184

calculated to evaluate bacterial community similarity between samples.28 As UniFrac, which

185

incorporates phylogenetic distance information, is impractical for the fungal ITS gene region,29

186

Jaccard (membership-based) and Bray-Curtis (structure-based) matrices were calculated using R

187

package Ecodist based on fungal genera.30,

188

performed to evaluate variation between each sampling group and pairwise ADONIS analysis

189

(permutational multivariate analysis of variance)32 implemented within QIIME18, 33 was used to

190

further assess the statistical significance of differences in community structure among sampling

191

groups. The raw sequencing data are publically available on figshare under DOI:

192

10.6084/m9.figshare.5007743 and 10.6084/m9.figshare.5007737.

31

Principal coordinate analysis (PCoA) was

193 194

2.5. qPCR of Total Fungi, Total Bacteria, Acanthamoeba spp., and PCR of Aspergillus

195

fumigatus

196

Abundances of total fungi, total bacteria, and Acanthamoeba spp. were quantified by qPCR using

197

SYBR Green and fungal universal primers ITS1FI2 and ITS2,34, 35 universal bacterial primers

198

targeting the 16S rRNA gene region,36 and Acanthamoeba spp. primers AcantF900 and

199

AcantR1100.37 Detailed qPCR conditions can be found in the Supplementary Information. No-

200

template qPCR negative controls were included in all qPCR runs to ensure reagents and

8 ACS Paragon Plus Environment

Page 9 of 31

Environmental Science & Technology

201

equipment were not contaminated. Inhibition was evaluated by spiking sample DNA into a

202

known amount of standard; no qPCR inhibition was observed. Based on the standard curves and

203

volume of sample filtered, the lower detection limit of bacterial qPCR was 0.54 P. fluorescens

204

genome equivalents/mL; the lower detection limit of fungal qPCR was 0.22 A. fumigatus

205

genome equivalents/mL; and the lower detection limit of Acanthamoeba spp. qPCR was 0.002 A.

206

castellanii genome equivalents/mL (1.2 A. castellanii 18S rRNA gene copies/mL).

207 208

The presence of Aspergillus fumigatus was determined using previously described PCR primers

209

targeting the A. fumigatus alkaline protease gene.38,

210

(ATCC 34506) genomic DNA as template, and no-template negative controls were included. All

211

positive controls were positive and negative controls were negative for contamination. Detailed

212

PCR conditions can be found in Supplementary Information.

39

Positive controls using A. fumigatus

213 214

3. Results and Discussion

215

3.1. Sampling and Analysis Overview

216

In total, 64 samples (duplicate samples of four locations over eight months) were collected from

217

a single anonymous drinking water treatment plant. A process diagram of the drinking water

218

treatment plant with sampling locations indicated is shown in Figure S1. All duplicate samples

219

were subjected to universal fungal ITS1 PCR, with 61 out of 64 samples resulting in positive

220

amplification. ADONIS analysis based on Jaccard and Bray-Curtis matrices comparing

221

significance of community difference between the two replicates groups showed no significant

222

community variation was caused by duplicate samples (p-values > 0.05), which is consistent

223

with previous study.40 Previous research has demonstrated high reproducibility of 16S rRNA

9 ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 31

224

sequencing results between replicate drinking water samples and technical replicates.16, 40 Thus,

225

one replicate of each sample (32 samples) was subjected to universal 16S rRNA PCR, with all

226

samples resulting in positive amplification. All negative controls indicated no contamination of

227

reagents or equipment. In total, 446,902 ITS1 sequence reads (average 7326±7866 reads per

228

sample) and 548,098 16S rRNA sequence reads (average 17128±6181 reads per sample) were

229

used for taxonomic classification after demultiplexing and quality trimming.

230 231

3.2. qPCR of Total Bacteria, Total Fungi, Acanthamoeba spp., and PCR of Aspergillus

232

fumigatus

233

Quantitative PCR was utilized to determine microbial abundance through the treatment process.

234

The average total bacterial abundance was 5.5±0.4 log10 genome equivalents/mL in raw water,

235

3.4±0.4 log10 genome equivalents/mL post-sedimentation, 2.9±0.3 log10 genome equivalents/mL

236

post-filtration, and 1.2±0.3 log10 genome equivalents/mL post-disinfection (Table 1). Bacterial

237

abundance significantly decreased after each unit operation (all p-values < 0.05). The average

238

fungal abundance was 3.4±0.5 log10 genome equivalents/mL in raw water, 1.6±0.5 log10 genome

239

equivalents/mL post-sedimentation, 0.9±0.4 log10 genome equivalents/mL post-filtration, and

240

0.2±0.4 log10 genome equivalents/mL post-disinfection (Table 1). A significant decrease in

241

fungal abundance was observed post-sedimentation and post-filtration (p-values < 0.05). The

242

average post-disinfection fungal abundance was lower than post-filtration, but without a

243

statistically significant change (p > 0.05).

244 245

The inability of molecular methods to discriminate DNA from viable and non-viable cells may

246

inflate the bacterial and fungal abundance detected, especially in the post-disinfection samples.

247

To confirm observed results, we collected additional samples and analyzed culturable

10 ACS Paragon Plus Environment

Page 11 of 31

Environmental Science & Technology

248

heterotrophic bacteria and filamentous fungi. The culture analysis showed a similar decreasing

249

trend through the treatment process train for both heterotrophic plate count and filamentous fungi

250

colony counts (Table 1). Both the qPCR and subsequent culturing results suggest the treatment

251

process significantly reduced bacterial and fungal abundance, agreeing with previous culture-

252

based studies.41, 42

253 254

As the eukaryote Acanthamoeba spp. is an important human pathogen43 and may host other

255

opportunistic pathogens such as Legionella spp.,44,

256

presence and abundance of Acanthamoeba spp.. Six out of eight raw water samples and one out

257

of eight post-filtration samples had detectable level of Acanthamoeba spp., with average

258

abundances of -1.36±1.61 log10 genome equivalents/mL and -1.51 log10 genome equivalents/mL,

259

respectively. None of the post-sedimentation or post-disinfection samples had detectable

260

Acanthamoeba spp.. Previous studies demonstrated the isolation of Acanthamoeba spp. from

261

post-filtration water samples46 and drinking water distribution systems.8,

262

Acanthamoeba spp. in raw water of the current study indicate that Acanthamoeba spp. likely

263

originate from an environmental source and persist within the filtration process, and further study

264

investigating the diversity and association of Acanthamoeba spp. with filter media would be

265

helpful for understanding the role of media filtration in Acanthamoeba spp. persistence through

266

drinking water treatment. Similarly, as Aspergillus fumigatus is among the most notable human

267

fungal pathogens and has previously been isolated from drinking water,14 we conducted PCR

268

targeting A. fumigatus. No sample had detectable A. fumigatus based on gel visualization of PCR

269

products.

45

qPCR was conducted to evaluate the

47

The presence of

270

11 ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 31

271

3.3. Microbial Diversity Through The Treatment Process

272

We assessed the beta-diversity for both bacterial and fungal communities, which is a measure of

273

the degree of microbial community similarity between samples. For both bacterial and fungal

274

communities, the mean community dissimilarity distances between samples at different water

275

treatment unit operations were all significantly higher than the mean distances across samples of

276

different months at the same sampling location (p-values < 0.05) (Table S2). This result indicates

277

that water treatment unit operations more strongly drove both bacterial and fungal community

278

structure than sampling time.

279 280

Non-metric multidimensional scaling (NMDS) analysis of the bacterial community based on

281

unweighted UniFrac indicated that samples clustered within each sampling site and the bacterial

282

community structure shifted following each unit operation (Figure 1). The statistical significance

283

of the community shift following each unit operation was confirmed with pairwise ADONIS

284

analysis (p-values < 0.05, Table S3). NMDS plots and ADONIS analysis based on weighted

285

UniFrac demonstrated a similar trend of community variation, except that the variation between

286

post-sedimentation and post-filtration samples was not statistically significant (Figure 1, Table

287

S3).

288 289

NMDS plot based on the Jaccard distance indicated that the fungal community structure shifted

290

following sedimentation and media filtration. Post-disinfection samples clustered with post-

291

filtration samples, indicating that the fungal community structure in post-disinfection samples

292

was similar to post-filtration samples, and no further community structural variation occurred

293

following disinfection (Figure 1). Pairwise ADONIS based on the Jaccard distance further

12 ACS Paragon Plus Environment

Page 13 of 31

Environmental Science & Technology

294

indicated that the fungal community structure significantly shifted after sedimentation and media

295

filtration (p-values < 0.05, Table S4), but no significant change occurred post-disinfection (p >

296

0.05, Table S4). NMDS plots and ADONIS based on Bray-Curtis distance demonstrated the

297

same pattern of fungal community variation between unit operations (Figure 1, Table S4).

298 299

We also assessed the variation in community richness (total number of taxonomic units) through

300

the treatment process. The number of observed bacterial OTUs significantly decreased through

301

each step of the treatment process (p-values < 0.05, Table 1). The number of fungal genera

302

significantly decreased post-sedimentation and post-filtration (p-values < 0.05), but no

303

significant decrease occurred post-disinfection (p > 0.05, Table 1). The decrease in bacterial

304

richness through the water treatment process agrees with previous study.11 These results also

305

identify a similar trend to bacteria for fungal richness with the exception of no significant change

306

post-disinfection. Pearson correlation between water quality parameters (pH, conductivity, free

307

chlorine) and bacterial OTU number as well as fungal genera number were evaluated. However,

308

no significant correlation was observed (- 0.5 < r-values < 0.5; p-values > 0.05). The only

309

exception was that a significant negative correlation was observed between the number of

310

bacterial OTUs and free chlorine concentration (r = - 0.52; p < 0.01). Similarly, no significant

311

correlation between water quality parameters and dominant taxa relative abundances was found

312

(- 0.5 < r-values < 0.5; p-values > 0.05), consistent with previous studies.40, 48

313 314

Our results demonstrated that the bacterial community shifted following each treatment step, and

315

disinfection shaped the final bacterial community in treated drinking water. The final fungal

316

community structure in treated drinking water was primarily shaped by media filtration, as no

13 ACS Paragon Plus Environment

Environmental Science & Technology

Page 14 of 31

317

significant fungal community change was observed following disinfection. Previous research

318

identified a similar trend with the bacterial community structure changing in response to addition

319

of secondary disinfection in a premise plumbing system,49 with no change in fungal community

320

structure found after the initiation of secondary disinfection.50 Media filtration has been

321

previously shown to control the downstream bacterial community.10 The current result agrees

322

with previous research, and further suggests free chlorine disinfection shifts the bacterial

323

community entering the water distribution system, but does not pose a significant effect on the

324

post-filtration fungal community.

325 326

3.4. Temporal Beta-Diversity Within Water Treatment Processes

327

Temporal beta-diversity following each treatment process unit operations was evaluated by

328

comparing the average dissimilarity distance within each sampling location. Mean temporal

329

bacterial community dissimilarity distances became apparently higher post-filtration and post-

330

disinfection; and significantly higher post-disinfection regarding membership based dissimilarity

331

(Unweighted UniFrac, p < 0.05) (Table S5). For the fungal community, membership based mean

332

temporal community dissimilarity distance significantly increased post-filtration (Jaccard, p
0.05).

were the most

361

15 ACS Paragon Plus Environment

Environmental Science & Technology

Page 16 of 31

362

Within Proteobacteria, Alpha-, Beta-, and Gammaproteobacteria were the most abundant classes

363

and were each detected in all samples. The average relative abundance of Alphaproteobacteria

364

was 9.5±2.7% in raw water, 10.4±7.1% post-sedimentation, 30.5±30.1% post-filtration, and

365

25.6±25.4% post-disinfection. A trend of increasing relative abundance post filtration was

366

observed, but was not statistically significant (p > 0.05). Alphaproteobacteria was previously

367

found to be a dominant bacterial class within filter media, suggesting media filtration may

368

provide a survival niche for Alphaproteobacteria, or conversely that Alphaproteobacteria are less

369

efficiently removed by the filtration process.10, 12 In this study, Sphingomonas spp. was the most

370

abundant genus within Alphaproteobacteria, and LEfSe analysis indicates that Sphingomonas

371

spp. was the most enriched taxa in post-filtration samples (LDA score = 5.3, Table S6).

372

Sphingomonas spp. is a common member of drinking water system biofilms and has an

373

important role in biofilm formation,54-56 suggesting that Sphingomonas spp. may form a biofilm

374

within filtration media and subsequently influence the post filtration bacterial community

375

structure. Betaproteobacteria had an average relative abundance of 41.1±10.1% in raw water,

376

35.5±11.9% post-sedimentation, 31.3±15.8% post-filtration, and 19.6±15.0% post-disinfection,

377

with a significant decreasing trend comparing raw water to post-disinfection water (p < 0.05)

378

(Figure 2). The average relative abundance of Gammaproteobacteria was 5.2±2.4% in raw water,

379

8.5±11.2% post-sedimentation, 4.1±4.2% post-filtration, and significantly increased to

380

17.1±11.3% post-disinfection (p < 0.05). This increasing trend of Gammaproteobacteria relative

381

abundance was confirmed by LEfSe analysis, showing Gammaproteobacteria was enriched post-

382

disinfection (LDA score > 2.0, Table S6). Acinetobacter spp. was found to be the most abundant

383

Gammaproteobacteria genus to increase in relative abundance post-disinfection (LDA score =

384

4.8, Table S6). Several strains of Acinetobacter spp. were previously found to survive without

16 ACS Paragon Plus Environment

Page 17 of 31

Environmental Science & Technology

385

significant inactivation when exposed to 4 mg/L free chlorine for 2 minutes.57 Previous study

386

also

387

Betaproteobacteria to the same level of free chlorine.58 This differential resistance to disinfection

388

likely drove the shift in the composition of Proteobacteria phylum during disinfection.

demonstrated

higher

resistance

by

Gammaproteobacteria

than

Alpha-

and

389 390

The Actinobacteria class was the dominant class within the Actinobacteria phylum and was

391

detected in all samples with an average relative abundance of 17.0±9.4% in raw water,

392

12.3±8.6% post-sedimentation, 6.5±4.8% post-filtration, and 6.5±5.4% post-disinfection. A

393

significant trend of decreasing relative abundance was observed comparing raw water to post-

394

disinfection water (p-values < 0.05) (Figure 2). This finding is consistent with a previous study

395

which reported decreased relative abundance of Actinobacteria comparing treated drinking water

396

to raw water.10 Actinobacteria relative abundance was not decreased comparing post-filtration to

397

post-disinfection samples, suggesting the potential disinfection resistance of this bacterial group.

398 399

The Bacteroidetes phylum had an average relative abundance of 15.4±2.7% in raw water,

400

24.7±8.9% post-sedimentation, 23.8±15.8% post-filtration, and 2.0±2.7% post-disinfection.

401

Bacteroidetes had significantly increased relative abundance post-sedimentation (p < 0.05), no

402

significant shift following filtration, and a significantly decreased relative abundance post-

403

disinfection (p < 0.05). Cytophagia and Saprospirae were the most abundant bacterial classes

404

within the Bacteroidetes phylum and reflected the same trend of relative abundance as the

405

Bacteroidetes phylum. The decreasing trend of relative abundance for Bacteroidetes post-

406

disinfection suggested this phylum might be more sensitive to chlorination, and further

17 ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 31

407

evaluation of disinfection kinetics would help to explain the mechanism for the observed

408

negative selection.

409 410

The Cyanobacteria and Firmicutes phyla shared similar trends of low relative abundance without

411

a statistically significant change (p-values > 0.05) before disinfection, but significantly increased

412

relative abundance after disinfection (p-values < 0.05). LEfSe analysis also indicated that

413

Cyanobacteria and Firmicutes were significantly more abundant post-disinfection (LDA scores >

414

2.0, Table S6). The average relative abundance of Cyanobacteria was 1.0±1.3% in raw water,

415

3.3±4.9% post-sedimentation, 0.7±1.1% post-filtration, and statistically significantly increased to

416

10.7±10.4% post disinfection (p < 0.05). The increase of Cyanobacteria was associated with the

417

4C0d-2 class, which had a significantly increased post-disinfection relative abundance (p < 0.05)

418

and became the dominant Cyanobacterial class post-disinfection, suggesting potential resistance

419

to chlorination by this non-photosynthetic Cyanobacteria class.59 This class has previously been

420

identified to be abundant in water systems, including drinking water treatment plants60 and a

421

simulated drinking water distribution system.61 Firmicutes had an average relative abundance of

422

1.1±0.9% in raw water, 1.0±1.3% post-sedimentation, 0.4±0.7% post-filtration, and statistically

423

significantly increased to 10.8±7.5% post-disinfection (p < 0.05). Within the Firmicutes phylum,

424

Bacilli and Clostridia were the most abundant classes and a statistically significant increase in

425

relative abundance was observed for Bacilli and Clostridia following disinfection (p-values
0.05) (Figure 2). Leptospirillum spp.

434

was identified in all eight post-disinfection samples but was only detected in 10 out of 24

435

(41.7%) of pre-disinfection samples. LEfSe analysis output also indicated that Leptospirillum

436

spp. were enriched in post-disinfection samples (LDA score = 4.7, Table S6), suggesting

437

Leptospirillum spp. to be potentially disinfection resistant. The Leptospirillum genus, which

438

contains iron-oxidizing species,64 was previously found in a cast iron water pipe biofilm and has

439

been suggested to contribute to iron pipe corrosion.65 Further evaluation of the resistance of

440

Leptospirillum spp. to disinfectant may explain its presence and inform its control in drinking

441

water distribution systems.

442 443

3.6. Fungal Taxonomy

444

Fungal ITS1 sequences were classified to 1043 fungal genera across all samples. Three phyla

445

were detected in all samples: Ascomycota, Basidiomycota, and Zygomycota. These three phyla

446

represented an average relative abundance of 98.4±1.4%. The relative abundance of fungal

447

genera detected in all samples is shown in Figure 3.

448 449

The average relative abundance of the Ascomycota phylum was 66.2±19.3% in raw water,

450

55.9±16.8% post-sedimentation, 56.5±19.8% post-filtration, and 70.7±23.9% post-disinfection,

451

with no statistically significant change through the water treatment process (p-values > 0.05).

452

Within the Ascomycota phylum, the most abundant genera were Penicillium, Aspergillus,

19 ACS Paragon Plus Environment

Environmental Science & Technology

Page 20 of 31

453

Alatospora, Taphrina, and Tuber, each detected in all samples. Penicillium spp. had an average

454

relative abundance of 1.5±1.8% in raw water, 9.3±13.5% post-sedimentation, 25.2±32.5% post-

455

filtration, and 41.3±38.3% post-disinfection, with a significant increasing trend comparing raw

456

water to post-disinfection water (p < 0.05) (Figure 3). The average relative abundance of

457

Aspergillus spp. was 0.6±0.8% in raw water, 1.2±1.1% post-sedimentation, 1.2±1.2% post-

458

filtration, and 3.6±2.1% post-disinfection. Similar to Penicillium spp., the relative abundance of

459

Aspergillus spp. significantly increased from raw water to post-disinfection water (p < 0.05);

460

including a significant increase in relative abundance from post-filtration to post-disinfection

461

samples (p < 0.05) (Figure 3). LEfSe analysis also shows Penicillium spp. and Aspergillus spp.

462

were significantly enriched post-disinfection (LDA scores > 2.0, Table S7). Alatospora spp. and

463

Taphrina spp., with overall average relative abundances of 2.1±3.9% and 2.8±5.8% respectively,

464

did not demonstrate a significant change in relative abundance through the treatment process (p-

465

values > 0.05). A decreasing trend of Tuber spp. relative abundance was observed from

466

19.0±24.5% in raw water to 8.0±7.0% post-disinfection, yet was not statistically significant (p-

467

values > 0.05).

468 469

The average relative abundance of the Basidiomycota phylum was 23.1±19.2% in raw water,

470

33.9±19.1% post-sedimentation, 35.7±20.1% post-filtration, and 22.0±17.1% post-disinfection,

471

with no significant change through the treatment process (p-values > 0.05). Mrakia and

472

Peniophora were the two genera within the Basidiomycota phylum detected in all samples.

473

Peniophora spp., a soil fungus belonging to Basidiomycota,66 was previously detected in tap

474

water samples.50,

67

No significant change in relative abundance across the water treatment

20 ACS Paragon Plus Environment

Page 21 of 31

Environmental Science & Technology

475

process was observed for either Mrakia or Peniophora (p-values > 0.05), mirroring the trend on

476

the phylum level (Figure 3).

477 478

The average relative abundance of the Zygomycota phylum was 8.4±7.1% in raw water,

479

8.6±8.7% post-sedimentation, 6.7±10.5% post-filtration, and 6.1±7.8% post-disinfection. No

480

statistically significant change was observed through the treatment process (p-values > 0.05).

481

The Schizangiella and Basidiobolus genera belonging to Zygomycota were detected in all

482

samples, and relative abundance significantly decreased through the treatment process (p-values

483

< 0.05).

484 485

Isolation of the fungal phyla Ascomycota, Basidiomycota, and Zygomycota from untreated

486

surface water has previously been reported.68-70 The effect of centralized water treatment process

487

on the fungal community structure in drinking water has not been previously evaluated; however,

488

Penicillium spp. and Aspergillus spp. have been widely isolated from drinking water.13, 71, 72 Our

489

results demonstrated that the relative abundance of Penicillium spp. and Aspergillus spp.

490

significantly increased through the water treatment process, especially post-disinfection,

491

suggesting that these two genera were less efficiently removed by the conventional water

492

treatment process or were more resistant to the selection pressure posed by water treatment

493

processes. Penicillium and Aspergillus produce unicellular spores with relatively small size (1.8-

494

6 µm × 2-4.5 µm).73-75 The small spore size may facilitate the passage of Penicillium and

495

Aspergillus through sedimentation and media filtration, while other fungal groups with larger

496

spore sizes were removed. Further research is needed to investigate the role of cell size in the

497

community structure change along the drinking water treatment process. Stronger resistance to

21 ACS Paragon Plus Environment

Environmental Science & Technology

Page 22 of 31

498

oxidative disinfectants by waterborne fungi, especially Penicillium and Aspergillus, comparing

499

to common waterborne bacteria has been demonstrated by previous studies.51,

500

insignificant post-disinfection fungal community variation is likely related to the resistance to

501

chlorination. Further research focusing on fungal disinfection and ecology in municipal drinking

502

water could benefit the control of waterborne fungi, especially for Aspergillus as this genus

503

contains potential human pathogens such as A. fumigatus and A. flavus.77

52, 76

The

504 505

3.7. Implications

506

In the current study we applied molecular approaches to investigate bacterial and fungal

507

abundance, diversity, and taxonomy through centralized surface water treatment. Bacterial and

508

fungal abundance and richness were observed to significantly decrease following sedimentation

509

and media filtration, but only the bacterial community was affected by disinfection. Assessing

510

the bacterial and fungal community beta-diversity, each treatment step was found to shift

511

bacterial and fungal community structure, with the exception of disinfection on the fungal

512

community structure.

513 514

The current study shows disinfection by free chlorine shaped the bacterial community entering

515

the distribution system. Differential sensitivity to disinfection likely drives the selection of

516

bacterial taxa in final product water. The fungal community was found to be more influenced by

517

physical processes (sedimentation and media filtration), highlighting the difference in

518

community shifts between bacteria and microbial eukaryotes during drinking water treatment.

519

The observed difference in community drivers between bacteria and fungi suggests the necessity

520

of considering a broader microbiome including both bacteria and eukaryotes when investigating

22 ACS Paragon Plus Environment

Page 23 of 31

Environmental Science & Technology

521

drivers of the drinking water microbiome. Ultimately, this study demonstrates that centralized

522

drinking water treatment shapes both the bacterial and eukaryotic microbiome of finished

523

drinking water, suggesting the future potential to harness these processes to select for a more

524

desirable drinking water microbial community.

525 526 527 528 529 530 531 532 533 534

Supporting Information Additional descriptions of experimental methods, statistics of 16S rRNA sequencing reads not matching reference database (Table S1), comparison of mean community dissimilarity distance across water treatment processes and time (Table S2 & S5), results of ADONIS analysis (Table S3 & S4), a water treatment processes diagram with sampling locations (Figure S1), principal coordinate analysis plots (Figure S2 & S3), detailed LEfSe analysis output (Table S6 & S7), and additional bacterial beta diversity analysis based on Jaccard and Bray-Curtis matrices (Table S8 – S10, Figure S4).

535

Acknowledgement

536 537 538 539 540 541

The authors wish to thank the staff at the anonymous drinking water treatment plant for assistance with sample collection.

TOC Art

542

23 ACS Paragon Plus Environment

Environmental Science & Technology

543 544 545 546 547

548 549 550 551 552 553 554 555 556

Page 24 of 31

Tables

Table 1. Bacterial and fungal diversity and abundance and water quality parameters. Detection equivalent/mL for bacteria qPCR and -0.66 log10 genome equivalent/mL for fungi. PostRaw water sedimentation Bacterial community Total number of observed OTUs 6152 3736 Average number of observed OTUs per sample 1744.5±482.6 968.3±364.9 Average bacterial abundance (log10 genome equivalent/mL) 5.5±0.4 3.4±0.4 Average heterotrophic plate count (CFU/mL) 501.7±120.2 0.5±0.7 Fungal community Total number of observed genera 980 585 Average number of observed genera per sample 396.5±173.4 188.9±60.8 Average fungal abundance (log10 genome equivalent/mL) 3.4±0.5 1.6±0.5 Average total filamentous fungal colony (CFU/L) 17750±13081.5 52.5±31.8 pH 7.5±0.4 7.5±0.5 Conductivity (µS/cm) 279.7±79.7 304.5±97.8 Free chlorine residual (mg/L) ND ND ND: not detected

limits were -0.27 log10 genome Postfiltration

Postdisinfection

2207 546.4±238.1

995 192.1±68.7

2.9±0.3 ND

1.2±0.3 ND

315 91.1±51.5 0.9±0.4 7.8±1.1 7.6±0.4 311.9±66.0 ND

294 91.5±28.6 0.2±0.4 3.3±1.1 7.6±0.2 303.1±75.8 0.5±0.3

24 ACS Paragon Plus Environment

Page 25 of 31

557 558 559 560 561 562

Environmental Science & Technology

Figures

1.2

Unweighted UniFrac (Bacterial)

Weighted UniFrac (Bacterial)

2

1

Raw water

1.5

0.8 1 0.4

NMDS2

NMDS2

0.6

0.2 0

Post-sedimentation

0.5 Raw water 0

Post-filtration

-0.5

-0.2 -1 -0.4

Post-disinfection -1.5

-0.6 -0.8

stress = 0.07

-2

stress = 0.11 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3 NMDS1

-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3 NMDS1 2.5

2

2

1.5

1.5

1

1

0.5

NMDS2

NMDS2

Jaccard (Fungal) 2.5

0 -0.5

0

-0.5

-1

-1

-1.5

-1.5

-2

563 564 565 566 567 568

0.5

stress = 0.14 -2 -2 -1.5 -1 -0.5 0 0.5 NMDS1

1

1.5

2

2.5

Bray-Curtis (Fungal)

Ra w wa ter

stress = 0.18 -2 -1.5 -1 -0.5 0 0.5 NMDS1

1

1.5

2

Figure 1. Non-metric Multidimensional Scaling (NMDS) plot of the bacterial community (based on unweighted and weighted UniFrac distances) and fungal community (based on Jaccard and Bray-Curtis distances).

25 ACS Paragon Plus Environment

Environmental Science & Technology

Page 26 of 31

Relative abundance (%)

80

60

Proteobacteria_Alphaproteobacteria Proteobacteria_Betaproteobacteria Proteobacteria_Gammaproteobacteria Actinobacteria_Actinobacteria Bacteroidetes_Cytophagia Bacteroidetes_[Saprospirae] Cyanobacteria_4C0d-2 Firmicutes_Bacilli Firmicutes_Clostridia Nitrospirae_Nitrospira

40

20

0 RW

569 570 571 572 573 574 575 576

P_S

P_F

P_D

Figure 2. Relative abundance of dominant bacterial classes through the water treatment processes. Boxplot range defined as the first and third quartile; outliers defined as 1.5 times the interquartile range greater than the third quartile or less than the first quartile.

26 ACS Paragon Plus Environment

Environmental Science & Technology

Relative abundance (%)

Page 27 of 31

75 Ascomycota_Penicillium Ascomycota_Aspergillus Ascomycota_Alatospora Ascomycota_Taphrina Ascomycota_Tuber Basidiomycota_Mrakia Basidiomycota_Peniophora Zygomycota_Schizangiella Zygomycota_Basidiobolus

50

25

0 R_W

577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601

P_S

P_F

P_D

Figure 3. Relative abundance of fungal genera detected in all samples through the water treatment process. Boxplot range defined as the first and third quartile; outliers defined as 1.5 times the interquartile range greater than the third quartile or less than the first quartile.

References: 1. Berry, D.; Xi, C.; Raskin, L., Microbial ecology of drinking water distribution systems. Curr. Opin. Biotechnol. 2006, 17, (3), 297-302. 2. McGuire, M. J., Eight revolutions in the history of US drinking water disinfection. J. - Am. Water Works Assoc. 2006, 98, (3), 123-149. 3. Proctor, C. R.; Hammes, F., Drinking water microbiology — from measurement to management. Curr. Opin. Biotechnol. 2015, 33, (0), 87-94. 4. Lautenschlager, K.; Hwang, C.; Liu, W.-T.; Boon, N.; Köster, O.; Vrouwenvelder, H.; Egli, T.; Hammes, F., A microbiology-based multi-parametric approach towards assessing biological stability in drinking water distribution networks. Water Res. 2013, 47, (9), 3015-3025. 5. Hammes, F.; Berney, M.; Wang, Y.; Vital, M.; Köster, O.; Egli, T., Flow-cytometric total bacterial cell counts as a descriptive microbiological parameter for drinking water treatment processes. Water Res. 2008, 42, (1–2), 269-277. 6. Doggett, M. S., Characterization of Fungal Biofilms within a Municipal Water Distribution System. Appl. Environ. Microbiol. 2000, 66, (3), 1249-1251. 7. Hageskal, G.; Lima, N.; Skaar, I., The study of fungi in drinking water. Mycol. Res. 2009, 113, 165172.

27 ACS Paragon Plus Environment

Environmental Science & Technology

602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649

Page 28 of 31

8. Delafont, V.; Brouke, A.; Bouchon, D.; Moulin, L.; Héchard, Y., Microbiome of free-living amoebae isolated from drinking water. Water Res. 2013, 47, (19), 6958-6965. 9. Wang, H.; Edwards, M. A.; Falkinham, I., J.O.; Pruden, A., Probiotic Approach to Pathogen Control in Premise Plumbing Systems? A Review. Environ. Sci. Technol. 2013, 47, 10117-10128. 10. Pinto, A. J.; Xi, C.; Raskin, L., Bacterial Community Structure in the Drinking Water Microbiome Is Governed by Filtration Processes. Environ. Sci. Technol. 2012, 46, 8851-8859. 11. Lin, W.; Yu, Z.; Zhang, H.; Thompson, I. P., Diversity and dynamics of microbial communities at each step of treatment plant for potable water generation. Water Res. 2014, 52, (0), 218-230. 12. Albers, C. N.; Ellegaard-Jensen, L.; Harder, C. B.; Rosendahl, S.; Knudsen, B. E.; Ekelund, F.; Aamand, J., Groundwater Chemistry Determines the Prokaryotic Community Structure of Waterworks Sand Filters. Environ. Sci. Technol. 2015, 49, (2), 839-846. 13. Hageskal, G.; Knutsen, A. K.; Gaustad, P.; Hoog, G. S.; Skaar, I., Diversity and Significance of Mold Species in Norwegian Drinking Water. Appl. Environ. Microbiol. 2006, 72, (12), 7586-7593. 14. Anaissie, E. J.; Stratton, S. L.; Dignani, M. C.; Lee, C. K.; Summerbell, R. C.; Rex, J. H.; Monson, T. P.; Walsh, T. J., Pathogenic molds (including Aspergillus species) in hospital water distribution systems: a 3-year prospective study and clinical implications for patients with hematologic malignancies. Blood 2003, 101, 2542-2546. 15. Anaissie, E. J.; Kuchar, R. T.; Rex, J. H.; Francesconi, A.; Kasai, M.; Müller, F. M.; Lozano-Chiu, M.; Summerbell, R. C.; Dignani, M. C.; Chanock, S. J.; Walsh, T. J., Fusariosis associated with pathogenic fusarium species colonization of a hospital water system: a new paradigm for the epidemiology of opportunistic mold infections. Clin. Infect. Dis. 2001, 33, (11), 1871-1878. 16. Smith, D. P.; Peay, K. G., Sequence Depth, Not PCR Replication, Improves Ecological Inference from Next Generation DNA Sequencing. PLoS ONE 2014, 9: e90234. doi: 10.1371/journal.pone.0090234. 17. Caporaso, J. G.; Lauber, C. L.; Walters, W. A.; Berg-Lyons, D.; Huntley, J.; Fierer, N.; Owens, S. M.; Betley, J.; Fraser, L.; Bauer, M.; Gormley, N.; Gilbert, J. A.; Smith, G.; Knight, R., Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012, 6, (8), 16211624. 18. Caporaso, J. G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F. D.; Goodrich, J.; al., e., QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 2010, 7, 335-336. 19. McDonald, D.; Price, M. N.; Goodrich, J.; Nawrocki, E. P.; DeSantis, T. Z.; Probst, A.; Andersen, G. L.; Knight, R.; Hugenholtz, P., An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012, 6, (3), 610-618. 20. Yamamoto, N.; Bibby, K., Clustering of Fungal Community InternalTranscribed Spacer (ITS) Sequence Data Obscures Taxonomic Diversity. Environ. Microbiol. 2014, 16, (18), 2491-2500. 21. Nilsson, R. H.; Bok, G.; Ryberg, M.; Kristiansson, E.; Hallenberg, N., A software pipeline for processing and identification of fungal ITS sequences. Source Code Biol. Med. 2009, 4, (1). 22. Altschul, S. F.; Gish, W.; Miller, W.; Myers, E. W.; Lipman, D. J., Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403-410. 23. Dannemiller, K. C.; Reeves, D.; Bibby, K.; Yamamoto, N.; Peccia, J., Fungal High-throughput Taxonomic Identification tool for use with Next-Generation Sequencing (FHiTINGS). J. Basic Microbiol. 2013, 54, (4), 315-321. 24. Wickham, H., Reshaping Data with the reshape Package. J. Stat. Softw. 2007, 21, (12), 1-20. 25. Wickham, H., ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag: New York, 2009. 26. Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W. S.; Huttenhower, C., Metagenomic biomarker discovery and explanation. Genome Biol. 2011, 12, (6), R60-R60. 27. Yamamoto, N.; Dannemiller, K.; Bibby, K.; Peccia, J., Identification accuracy and diversity reproducibility associated with ITS-based fungal taxonomic library preparation. Environ. Microbiol. 2014, 16, (9), 2764-2776. 28 ACS Paragon Plus Environment

Page 29 of 31

650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697

Environmental Science & Technology

28. Lozupone, C.; Knight, R., UniFrac: a New Phylogenetic Method for Comparing Microbial Communities. Appl. Environ. Microbiol. 2005, 71, (12), 8228-8235. 29. Lindahl, B. D.; Nilsson, R. H.; Tedersoo, L.; Abarenkov, K.; Carlsen, T.; Kjøller, R.; Kõljalg, U.; Pennanen, T.; Rosendahl, S.; Stenlid, J.; Kauserud, H., Fungal community analysis by high-throughput sequencing of amplified markers – a user's guide. New Phytol. 2013, 199, (1), 288-299. 30. Goslee, S. C.; Urban, D. L., The ecodist package for dissimilarity-based analysis of ecological data. J. Stat. Softw. 2007, 22, (7), 1-19. 31. R Development Core Team R: A language and environment for statistical computing, R Foundation for Statistical Computing: Vienna, Austria, 2010. 32. Anderson, M. J., A new method for non-parametric multivariate analysis of variance. Aust. Ecol. 2001, 26, (1), 32-46. 33. Oksanen, J.; Kindt, R.; Legendre, P.; O’Hara, B.; Stevens, M. H. H.; Oksanen, M. J.; Suggests, M., The vegan package. Community ecology package 2007, 10. 34. White, T.; Bruns, T.; Lee, S.; Taylor, J., Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In PCR protocols: a guide to methods and applications, Academic Press, San Diego: 1990; pp 315-322. 35. Schmidt, P. A.; Balint, M.; Greshake, B.; Bandow, C.; Rombke, J.; Schmitt, I., Illumina metabarcoding of a soil fungal community. Soil Biol. Biochem. 2013, 65, 128-132. 36. Maeda, H.; Fujimoto, C.; Haruki, Y.; Maeda, T.; Kokeguchi, S.; Petelin, M.; Arai, H.; Tanimoto, I.; Nishimura, F.; Takashiba, S., Quantitative real-time PCR using TaqMan and SYBR Green for Actinobacillus actinomycetemcomitans, Porphyromonas gingivalis, Prevotella intermedia, tetQ gene and total bacteria. FEMS Immunol. Med. Microbiol. 2003, 39, (1), 81-86. 37. Qvarnstrom, Y.; Visvesvara, G. S.; Sriram, R.; da Silva, A. J., Multiplex Real-Time PCR Assay for Simultaneous Detection of Acanthamoeba spp., Balamuthia mandrillaris, and Naegleria fowleri. J. Clin. Microbiol. 2006, 44, (10), 3589-3595. 38. Gabal, M. A.; El-Sherif, A. M.; Enany, M. S.; Soliman, S. S., A polymerase chain reaction ‘PCR’ for a quick diagnosis of aspergillosis. Mycoses 1999, 42, (9-10), 515-520. 39. Wu, Z.; Wang, X.-R.; Blomquist, G., Evaluation of PCR primers and PCR conditions for specific detection of common airborne fungi. J. Environ. Monit. 2002, 4, 377-382. 40. 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. 41. Charnock, C.; Otterholt, E.; Nordlie, A.-L., Treatment stage associated changes in cellular and molecular microbial markers during the production of drinking water at the Vansjø water works. Water Res. 2015, 81, 240-249. 42. Hoefel, D.; Monis, P. T.; Grooby, W. L.; Andrews, S.; Saint, C. P., Profiling bacterial survival through a water treatment process and subsequent distribution system. J. Appl. Microbiol. 2005, 99, (1), 175-186. 43. Khan, N. A., Acanthamoeba: biology and increasing importance in human health. FEMS Microbiol. Rev. 2006, 30, (4), 564-595. 44. Buse, H. Y.; Ashbolt, N. J., Differential growth of Legionella pneumophila strains within a range of amoebae at various temperatures associated with in-premise plumbing. Lett. Appl. Microbiol. 2011, 53, (2), 217-224. 45. Berry, D.; Horn, M.; Xi, C.; Raskin, L., Mycobacterium avium Infections of Acanthamoeba Strains: Host Strain Variability, Grazing-Acquired Infections, and Altered Dynamics of Inactivation with Monochloramine. Appl. Environ. Microbiol. 2010, 76, (19), 6685-6688. 46. Hoffmann, R.; Michel, R., Distribution of free-living amoebae (FLA) during preparation and supply of drinking water. Int. J. Hyg. Environ. Health 2001, 203, (3), 215-219. 29 ACS Paragon Plus Environment

Environmental Science & Technology

698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745

Page 30 of 31

47. Cateau, E.; Delafont, V.; Hechard, Y.; Rodier, M., Free-living amoebae: what part do they play in healthcare-associated infections? J. Hosp. Infect. 2014, 87, (3), 131-140. 48. Pinto, A. J.; Schroeder, J.; Lunn, M.; Sloan, W.; Raskin, L., Spatial-Temporal Survey and Occupancy-Abundance Modeling To Predict Bacterial Community Dynamics in the Drinking Water Microbiome. mBio 2014, 5, (3), e01135-14. 49. Baron, J. L.; Vikram, A.; Duda, S.; Stout, J. E.; Bibby, K. J., Shift in Microbial Ecology of a Hospital Hot Water System Following the Introduction of an On-Site Monochloramine Disinfection System. PLoS ONE 2014, 9, (7). 50. Ma, X.; Baron, J. L.; Vikram, A.; Stout, J. E.; Bibby, K., Fungal Diversity and Presence of Potentially Pathogenic Fungi in a Hospital Hot Water System Treated with On-Site Monochloramine. Water Res. 2015, 71, (15), 197-206. 51. Ma, X.; Bibby, K., Free chlorine and monochloramine inactivation kinetics of Aspergillus and Penicillium in drinking water. Water Res. 2017, 120, 265-271. 52. Pereira, V. J.; Fernandes, D.; Carvalho, G.; Benoliel, M. J.; San Romao, M. V.; Barreto Crespo, M. T., Free chlorine inactivation of fungi in drinking water sources. Water Res. 2013, 47, 517-523. 53. Ling, F.; Hwang, C.; LeChevallier, M. W.; Andersen, G. L.; Liu, W.-T., Core-satellite populations and seasonality of water meter biofilms in a metropolitan drinking water distribution system. ISME J 2015. 54. Hong, P.-Y.; Hwang, C.; Ling, F.; Andersen, G. L.; LeChevallier, M. W.; Liu, W.-T., Pyrosequencing Analysis of Bacterial Biofilm Communities in Water Meters of a Drinking Water Distribution System. Appl. Environ. Microbiol. 2010, 76, (16), 5631-5635. 55. Zhang, M.; Liu, W.; Nie, X.; Li, C.; Gu, J.; Zhang, C., Molecular Analysis of Bacterial Communities in Biofilms of a Drinking Water Clearwell. Microbes Environ. 2012, 27, (4), 443-448. 56. Bereschenko, L. A.; Stams, A. J. M.; Euverink, G. J. W.; van Loosdrecht, M. C. M., Biofilm Formation on Reverse Osmosis Membranes Is Initiated and Dominated by Sphingomonas spp. Appl. Environ. Microbiol. 2010, 76, (8), 2623-2632. 57. Karumathil, D. P.; Yin, H.-B.; Kollanoor-Johny, A.; Venkitanarayanan, K., Effect of Chlorine Exposure on the Survival and Antibiotic Gene Expression of Multidrug Resistant Acinetobacter baumannii in Water. Int. J. Environ. Res. Public Health 2014, 11, (2), 1844-1854. 58. Mathieu, L.; Bouteleux, C.; Fass, S.; Angel, E.; Block, J. C., Reversible shift in the α-, β- and γproteobacteria populations of drinking water biofilms during discontinuous chlorination. Water Res. 2009, 43, (14), 3375-3386. 59. Soo, R. M. In search of non-photosynthetic Cyanobacteria. University of Queensland, 2015. 60. Ji, P.; Parks, J.; Edwards, M. A.; Pruden, A., Impact of water chemistry, pipe material and stagnation on the building plumbing microbiome. PloS one 2015, 10, (10), e0141087. 61. Williams, M.; Domingo, J.; Meckes, M.; Kelty, C.; Rochon, H., Phylogenetic diversity of drinking water bacteria in a distribution system simulator. J. Appl. Microbiol. 2004, 96, (5), 954-964. 62. Ridgway, H. F.; Olson, B. H., Chlorine resistance patterns of bacteria from two drinking water distribution systems. Appl. Environ. Microbiol. 1982, 44, (4), 972-987. 63. Venczel, L. V.; Arrowood, M.; Hurd, M.; Sobsey, M. D., Inactivation of Cryptosporidium parvum oocysts and Clostridium perfringens spores by a mixed-oxidant disinfectant and by free chlorine. Appl. Environ. Microbiol. 1997, 63, (4), 1598-1601. 64. Rawlings, D.; Tributsch, H.; Hansford, G., Reasons why ‘Leptospirillum’-like species rather than Thiobacillus ferrooxidans are the dominant iron-oxidizing bacteria in many commercial processes for the biooxidation of pyrite and related ores. Microbiology 1999, 145, (1), 5-13. 65. Teng, F.; Guan, Y. T.; Zhu, W. P., Effect of biofilm on cast iron pipe corrosion in drinking water distribution system: Corrosion scales characterization and microbial community structure investigation. Corros. Sci. 2008, 50, (10), 2816-2823. 30 ACS Paragon Plus Environment

Page 31 of 31

746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776

Environmental Science & Technology

66. George, T. S.; Simpson, R. J.; Gregory, P. J.; Richardson, A. E., Differential interaction of Aspergillus niger and Peniophora lycii phytases with soil particles affects the hydrolysis of inositol phosphates. Soil Biol. Biochem. 2007, 39, (3), 793-803. 67. Buse, H. Y.; Lu, J.; Struewing, I. T.; Ashbolt, N. J., Eukaryotic diversity in premise drinking water using 18S rDNA sequencing: implications for health risks. Environ. Sci. Pollut. Res. 2013, 20, (9), 63516366. 68. Pereira, V. J.; Basílio, M. C.; Fernandes, D.; Domingues, M.; Paiva, J. M.; Benoliel, M. J.; Crespo, M. T.; San Romão, M. V., Occurrence of filamentous fungi and yeasts in three different drinking water sources. Water Res. 2009, 43, (15), 3813-3819. 69. Jasrotia, P.; Green, S. J.; Canion, A.; Overholt, W. A.; Prakash, O.; Wafula, D.; Hubbard, D.; Watson, D. B.; Schadt, C. W.; Brooks, S. C.; Kostka, J. E., Watershed-Scale Fungal Community Characterization along a pH Gradient in a Subsurface Environment Cocontaminated with Uranium and Nitrate. Appl. Environ. Microbiol. 2014, 80, (6), 1810-1820. 70. López-Archilla, A. I.; González, A. E.; Terrón, M. C.; Amils, R., Ecological study of the fungal populations of the acidic Tinto River in southwestern Spain. Can. J. Microbiol. 2004, 50, (11), 923-934. 71. Gonçalves, A. B.; Paterson, R. R. M.; Lima, N., Survey and significance of filamentous fungi from tap water. Int. J. Hyg. Environ. Health 2006, 209, 257-264. 72. Anaissie, E. J.; Stratton, S. L.; Dignani, M. C.; Summerbell, R. C.; Rex, J. H.; Monson, T. P.; Spencer, T.; Kasai, M.; Francesconi, A.; Walsh, T., Pathogenic aspergillus species recovered from a hospital water system: a 3-year prospective study. Clin. Infect. Dis. 2002, 34, (6), 780-789. 73. Yamamoto, N.; Bibby, K.; Qian, J.; Hospodsky, D.; Rismani-Yazdi, H.; Nazaroff, W. W.; Peccia, J., Particle-size distributions and seasonal diversity of allergenic and pathogenic fungi in outdoor air. ISME J. 2012, 6, 1801-1811. 74. Cole, G. T.; Samson, R. A., The conidia. In Mould allergy, Al-Doory, Y.; Domson, J. F., Eds. Philadelphia, 1984; pp 66-104. 75. Kwon-Chung, K. J.; Sugui, J. A., Aspergillus fumigatus—What Makes the Species a Ubiquitous Human Fungal Pathogen? PLoS Pathog. 2013, 9, (12), e1003743. 76. Hageskal, G.; Tryland, I.; Liltved, H.; Skaar, I., No simple solution to waterborne fungi: various responses to water disinfection methods. Water Sci. Technol.: Water Supply 2012, 12, (2), 220-226. 77. Alangaden, G. J., Nosocomial fungal infections: epidemiology, infection control, and prevention. Infect. Dis. Clin. North Am. 2011, 25, (1), 201-225.

777

31 ACS Paragon Plus Environment