Profile and Fate of Bacterial Pathogens in Sewage Treatment Plants

Aug 7, 2015 - *Phone: +852-28591968; fax: +852-25595337; e-mail: [email protected]. .... in a wastewater treatment plant and its receiving surface water...
0 downloads 0 Views 812KB Size
Subscriber access provided by UNIV OF PITTSBURGH

Article

Profile and fate of bacterial pathogens in sewage treatment plants revealed by high-throughput metagenomic approach Bing Li, Feng Ju, Lin Cai, and Tong Zhang Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.5b02345 • Publication Date (Web): 07 Aug 2015 Downloaded from http://pubs.acs.org on August 13, 2015

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 34

Environmental Science & Technology

1

Profile and fate of bacterial pathogens in sewage treatment plants

2

revealed by high-throughput metagenomic approach

3

Bing Li†#, Feng Ju‡#, Lin Cai‡, Tong Zhang‡*

4 5 6



#

Bing Li and Feng Ju contributed equally to this work

7

*

Corresponding author: Tong Zhang; email: [email protected]; phone: +852-28591968; fax: +852-25595337

8 9 10



Division of Energy & Environment, Graduate School at Shenzhen, Tsinghua University, China Environmental Biotechnology Lab, The University of Hong Kong, Hong Kong SAR, China

Bing Li: [email protected]; Feng Ju: [email protected]; Lin Cai: [email protected];Tong Zhang:[email protected]

Abi ot r op h ia d ef ec ti va Ac hr o m o b ac te r p ie c hau di i Ac h r om o ba ct er xy lo s ox i da ns Ac id a m i no co c cu s fer m en ta n s A ci neto ba c te r b au m an n ii Ac in et ob ac te r h a em ol y tic u s Ac ine tob ac te r j oh n so n ii A ci ne to b ac te r j un ii Aci ne to b ac te r l w of fii Ac in et ob ac ter r a d io r es is te n s A ct in om yc e s od o nt ol yt ic us Aer o c oc c us v ir i da ns Ae ro m on a s hy dr oph i la Ar c o ba c te r b ut zl er i B ac il lu s an th r a ci s B ac il lu s cer e us Baci llu s p um il us

MetaPhlAn

Ba ci ll u s s ub ti lis Bac il lu s th ur i ng i en si s Ba ct er o i de s ca c ca e Ba ct er o id e s eg g er t hi i B ac te r o id es f r ag il is

Metagenomic Phylogenetic Analysis

Sewage Treatment Plants (Influent, Effluent, AS, Biofilm, ADS) High-throughput Sequences

Metagenomic DNA

Illumina high-throughput sequencing (Hiseq 2000)

11 12 13

B ac te r o id es o v at us B ac te r oi de s p ec ti no p hi lu s Bac te r oi d es s te r c or i s B ac te r o id es th e ta io ta om i c r o n Bac te r o id es u ni fo r m is Bac te r oi d es v u lg at us B i fi do bac ter iu m de n ti u m B il oph B iloradw e ad te l slawaor v tiuhima Bo r de te ll a p a r ap er tu ss is Bo r d et el la p er t us s is C hr om obac te r iu m v i ol a ce u m C itr o b ac te r k o se r i C los tr i di um b otu lin u m C lo st r id iu m di ffi ci l e C l os tr id iu m p er f r in ge ns C oll in se ll a ae r ofa ci en s C om am o na s tes tost er o ni D e l f ti a a c id ov o r an s Eg ge r th e lla l en ta Ei ke n el la c or r o d en s En te r ob ac te r c a nc e r og e nu s En te r ob a ct er cl oac ae En te r o co cc u s ca s se li fl av us Ente r oc occ u s fa ec a l is Ente r oc occ u s fa ec i u m En te r oc o cc u s ga ll i n ar u m E sc h er i ch i a c ol i E ubac te r iu m li m os um Euba ct er i um r e c ta le Fi n ego ld ia m a gn a F u s ob ac ter iu m m or t ifer u m F u so b ac te r i um nu c l e at um F us ob a ct er i um u lc e r an s F us ob a ct er i um v a r iu m G ar dn er e l l a va g in al is G o r do ni a b r on c hi al is H ae m op hi lu s pa r a in flu e nz a e Kle bs ie ll a pne um o ni ae M a nnhe im ia h a em o ly tic a M e g am ona s hy p er m e ga le M y co b ac te r iu m ab s ce s su s M y co b ac te r i u m a vi um M y c ob act er i um bo vi s M y co bac ter iu m ka n sa s i i M yc ob a ct er i um l ep r ae M yc ob a ct er i um m a r in um M y c ob ac te r iu m s m e g m a ti s M yc ob a ct er i um u lc e r an s N ei ss e r ia el o ng a ta N e is se r ia f la ve s ce ns N eis s er i a gon or r h o ea e N ei ss e r ia m en in g iti di s N ei ss e r ia m uc o sa N e i s s er i a si cc a N e is s er i a su b fl a v a N o ca r d ia fa r c in i c a O c hr o b ac tr u m ant hr o pi P r op io ni ba c te r iu m ac n es P se u do m o n as f lu or e s ce n s P se ud o m o n as s tu tz e r i R a ls to n ia p ic ke tt i i R ho doc o cc u s eq ui R ho d oc oc c us e r y thr opo l is Ro th ia d e nt oc a r io sa Seba l del la te r m iti di s Sh i ge l la b o yd ii Sh ig el la d ys ent er i ae Shi g el la fl ex n er i Sh ig el l a s on n ei Sta ph y lo co c cu s a ur e u s Sta p hy lo c oc cus s ap r o ph y tic u s St en o tr o ph om o na s m al to phi l i a Str e p toc o cc u s ag al a ct ia e Str e p toc o cc u s an gi n os u s Str e pt oc o cc us bo v is S tr ep toc oc c us g o r do n ii Str e p to coc cu s m i ti s S tr ep to c oc c us m u ta ns Str e pt oc o cc u s sa li var iu s St r ep to c oc c us s u is Su tt er el la w a d sw o r th e ns is T suk am u re l l a pa u r om e ta bo la Ve il lo nel l a at yp ic a Ve il l o ne ll a di sp ar Ve ilVloi bnel a r iol ac pa ho rlev rulae V ib r io fu r n is si i Vi br i o m i m i cu s V i b r io pa r aha e m o ly ti cu s Vi br i o vu ln if ic us S ha t in - E ff . S h a ti n - In f. S h a tin - A S F o a m i ng

P a th o g e n sp e c i e s nu m be r

37

91

7

T o ta l a bu n da nc e (%)

11 . 8

2 3 .2

0 . 63

S h a t i n- A S S ta nl ey -A S St an le y - BF Sh at in - AD S SW H - AD S

4 0 . 33

TOC

14

1

ACS Paragon Plus Environment

27

24

6

3

0 . 64

0 .5 6

1 .2 3

0 .1 1

Environmental Science & Technology

15

Abstract: The broad-spectrum profile of bacterial pathogens and their fate in sewage treatment plants (STPs)

16

were investigated using high-throughput sequencing based metagenomic approach. This novel approach could

17

provide a united platform to standardize bacterial pathogen detection and realize direct comparison among

18

different samples. Totally, 113 bacterial pathogen species were detected in 8 samples including influent, effluent,

19

activated sludge (AS), biofilm, and anaerobic digestion sludge with the abundances ranging from 0.000095% to

20

4.89%. Among these 113 bacterial pathogens, 79 species were reported in STPs for the first time. Specially,

21

compared to AS in bulk mixed liquor, more pathogen species and higher total abundance were detected in upper

22

foaming layer of AS. This suggests that the foaming layer of AS might impose more threat to onsite workers

23

and citizens in the surrounding areas of STPs because pathogens in foaming layer are easily transferred into air

24

and cause possible infections. The high removal efficiency (98.0%) of total bacterial pathogens suggests that AS

25

treatment process is effective to remove most bacterial pathogens. Remarkable similarities of bacterial pathogen

26

compositions between influent and human gut indicated that bacterial pathogen profiles in influents could well

27

reflect the average bacterial pathogen communities of urban resident guts within the STP catchment area.

28 29

INTRODUCTION

30

Despite the remarkable improvement achieved in water and wastewater treatment technology as well as

31

management level, waterborne diseases which mainly result from the infection by pathogens are responsible for

32

4.0% of total death rate and 5.7% of the total burden of disease approximately. 1 Therefore, pathogens, including

33

bacteria, viruses, protozoa and helminths are considered to pose a major threat to human health around the 2

ACS Paragon Plus Environment

Page 2 of 34

Page 3 of 34

Environmental Science & Technology

34

world.2,3 Sewage treatment plant (STP) is one of the major sources for pathogens to enter into water

35

supplies/recreational waters and soil via water reuse and land application of sludge. 2-4,5 Considering the nature

36

of pathogens, that is, low infectious dose of some pathogens would associate with high level of community

37

risks,6 it is imperative to establish a rapid, accurate and universal detection method for the broad-spectrum

38

monitoring of pathogens in wastewater treatment plants.

39 40

Conventionally, simultaneous detection and identification of multiple kinds of pathogens using culture-based

41

methods is a time consuming and labor intensive undertaking. To get rid of such heavy work load, indicator

42

microorganisms, such as Escherichia coli and Clostridium perfringens, are widely used to estimate the overall

43

contamination level of pathogens in water and wastewater samples indirectly. However, due to the different

44

survival characteristics between indicator microorganisms and the majority of other pathogens, this method

45

could neither reflect the actual quantity of pathogens nor identify the pathogens explicitly. 4,7,8 In addition, the

46

ability of some pathogens to persist in a viable-but-non-culturable state will lead to a great underestimation of

47

pathogen contamination. 7 During the past 15 years, the development of molecular biology techniques provided

48

researchers some novel and sensitive approaches for pathogen detection. PCR, quantitative real-time PCR

49

(qPCR) and microarray assays have been successfully established to identify and quantify pathogens in water

50

and wastewater. 3,4,7,9 Nevertheless, PCR, qPCR approaches and microarray assays suffer from inherent low

51

throughput design; thereby numerous reactions are needed to detect diverse pathogens in environmental samples.

3

ACS Paragon Plus Environment

Environmental Science & Technology

52

On the other hand, since the limited available primers or probes usually could not cover all known pathogens,

53

these methods cannot reveal the broad profiles on the diversity and abundance of the overall pathogens.

54 55

Currently, the breakthrough of next-generation sequencing technology launched a new era for the

56

broad-spectrum monitoring and rapid quantification of pathogens. For example, 454 pyrosequencing of the 16S

57

rRNA gene was used to investigate the occurrence and quantity of more than 30 bacterial pathogens in biosolids,

58

activated sludge (AS), anaerobic digestion sludge (ADS), influent and effluent simultaneously.5,10-12

59

Undoubtedly, 454 pyrosequencing provides people new and important insights into the diversity and

60

contamination level of bacterial pathogens in various environmental samples. Nonetheless, this novel pathogen

61

detection approach still has the following limitations. Just like all the other PCR-based methods, the inherent

62

amplification biases during PCR step could not be avoided.13 Additionally, pyrosequencing of the 16S rRNA

63

gene will generate biased quantification due to the wide variability in copy numbers of the 16S gene even within

64

the same genus.

65

resolution at genus level which is neither accurate nor specific enough to investigate pathogens, 15 since not all

66

of the species under the pathogenic genera are pathogens.

14

Finally, the relatively short sequence length (~400 bp) can only provide the taxonomic

67

Metagenomic shotgun sequencing using high-throughput sequencing (HTS) technology, which typically

68

generates millions to billions of reads for the metagenomic DNA extracted from each environmental sample

69

directly, could provide a broad profile of microbial communities.15 Compared to PCR-based pyrosequencing of

70

the 16S rRNA gene, metagenomic approach is able to analyze the microbial community and estimate their 4

ACS Paragon Plus Environment

Page 4 of 34

Page 5 of 34

Environmental Science & Technology

71

abundances without PCR bias.14 Furthermore, metagenomic approach using unique clade-specific marker genes

72

can investigate the microbial composition at lower taxonomic level, i.e., species level. 15 This is especially

73

crucial for pathogen detection since the investigation of pathogens at genus level usually leads to some

74

inaccurate estimate.

75

approach to identify and evaluate the diversity of human viruses in freshwater lake, 18 reclaimed water

76

biosolid collected from STP. 20,21 However, so far, the detection of bacterial pathogens using metagenomic

77

approach is very limited although bacterial pathogens accounted for as high as 43% of the total pathogens

78

concerned by World Health Organization.22 Our previous study mainly focused on the establishment of bacterial

79

pathogen detection approach based on HTS-based metagenomic technique. 23 Compared to BLAST searching

80

against several databases using 16S rRNA genes and virulence factors, MetaPhlAn tool (Metagenomic

81

Phylogenetic Analysis) is the most reliable method because it relies on unique clade-specific marker genes

82

identified from 3,000 reference genomes for taxonomic classification (species level) and abundance

83

estimation.15,23 In addition, all the pathogen species detected by MetaPhlAn tool were also verified by the other

84

two approaches, further indicating the reliability of MetaPhlAn tool. 23 Therefore, MetaPhlAn tool was selected

85

in this study to investigate the profile of bacterial pathogens and their fate in STPs.

16.17

A few pioneer studies demonstrated the successful application of metagenomic 19

and

86

The objectives of this study were (1) to reveal the diversity and abundances of the bacterial pathogens in

87

various samples including influent, effluent, AS (foaming and non-foaming), biofilm, and anaerobic digestion

88

sludge (ADS); (2) to evaluate the fate of multiple bacterial pathogens simultaneously in secondary/biological

5

ACS Paragon Plus Environment

Environmental Science & Technology

89

treatment process and (3) to compare the overall bacterial pathogen community compositions in different

90

samples collected from STPs and human guts.

91 92 93

MATERIALS AND METHODS

Sample Collection

94

Two grab influent samples (~500 mL) and two grab effluent samples (~1000 mL) were collected from

95

Hong Kong Shatin STP (22.407° N, 114.214° E) in December of 2012 and January of 2013. Two DNA samples

96

of influents and effluents were pooled together respectively as composite samples named “Shatin-Inf.” and

97

“Shatin-Eff.”. AS samples (~25 mL) including “Shatin-AS Foaming” and “Shatin-AS” were collected from the

98

same aeration tank of Shatin STP in March of 2012 when the serious foaming phenomenon occurred.

99

“Shatin-AS Foaming” was gathered from the upper foaming layer while “Shatin-AS” was obtained from the

100

bulk mixed liquor which was about 30 centimeters below the wastewater surface. Biofilm (“Stanley-BF”,~25

101

mL) and the counterpart AS sample (“Stanley-AS”,~25 mL) were collected in March of 2012 from the hybrid

102

aeration tank of Hong Kong Stanley STP (22.219° N, 114.21° E). Biofilms grew on poly ethylene carriers with

103

25 mm diameter (Figure S1). Anaerobic digestion sludge samples (“SWH-ADS” and “Shatin-ADS”,~25 mL)

104

were obtained from the mesophilic anaerobic digestion tanks of Hong Kong Shek Wu Hui STP (22.51° N,

105

114.119° E) and Shatin STP in March of 2012, respectively. At sampling, AS, biofilm and ADS were kept in 50

106

ml sterilized centrifuge tubes and mixed with 100% ethanol at a 1:1 volume ratio for fixation. Influent and

107

effluent were kept in sterilized containers which were placed in an ice box and transported to the laboratory 6

ACS Paragon Plus Environment

Page 6 of 34

Page 7 of 34

Environmental Science & Technology

108

within 2 hours after sampling. When returned to the laboratory, influent samples were firstly centrifuged at 4000

109

rpm for 20 min to collect the cell pellets and then suspended and fixed using 50% ethanol. Effluent samples

110

were firstly filtered with 0.2 µm membrane to collect cell pellets and then fixed using 50% ethanol. At last, all

111

the fixed samples were stored at -20 °C before DNA extraction.23 The characteristics and treatment processes of

112

these three STPs were summarized in the Supporting Information (Table S1).

113 114

DNA Extraction and HTS

115

The DNA extraction was conducted using the FastDNA® Spin Kit for Soil (Q-biogene, CA, USA) and

116

DNA concentrations as well as purity were determined using microspectrophotometry (NanoDrop® ND-1000,

117

NanDrop Technologies, Willmington, DE). DNA extracted from two separate aliquots of each sample was

118

pooled to minimize the bias caused by variation of DNA extraction. Totally, eight prepared DNA samples (1

119

from influent, one from effluent, three from AS, one from biofilm and two from ADS) were sent out to Beijing

120

Genomics Institute (BGI) for shotgun library construction and Illumina HTS on HiSeq2000 platform. HTS was

121

performed using the strategy of index PE101+8+101 cycle (Paired-End sequencing, 101-bp reads and 8-bp

122

index sequence). The base-calling pipeline (Version Illumina Pipeline - 0.3) was used to process the raw

123

fluorescence images and call sequences. 24

124

Bioinformatics Analysis

125

To guarantee the quality of downstream analysis, raw reads containing three or more ambiguous

126

nucleotides, or with average quality score below 20, or with length less than 100 bp were removed using a 7

ACS Paragon Plus Environment

Environmental Science & Technology

127

custom-made script. On the other hand, the HiSeq2000 platform produces artificial replicates that are nearly

128

identical and failure to remove these replicated sequences could result in biased estimate.25-28 Therefore, the

129

reads after the above quality control processing were then filtered by another custom-made script to remove

130

replicated sequences.23,29 Subsequently, the remaining clean paired-end reads were paired-aligned by the third

131

custom-made script to merge 10-50 bp overlapped paired-end reads into 150-190 bp illumina Tags (iTags).23 The

132

numbers of iTags belonging to each sample before and after normalization was summarized in Table S2. After

133

that, MetaPhlAn (Version 1.6.0) was applied to conduct taxonomic classification at species level and quantify

134

abundance by mapping metagenomic reads against a catalog of clade-specific marker sequences currently

135

spanning the bacterial and archaeal phylogenies.15 Except for the threshold of the e-value of 1 × 10–15, all the

136

other parameters of MetaPhlAn utilized default settings. Finally, the output file of MetaPhlAn including

137

taxonomic results and the corresponding abundance of each identified species will be compared to our bacterial

138

pathogen database which contains 538 pathogenic species.30 The mapped pathogenic species names and their

139

abundances will be extracted as a csv file using the fourth custom-made script. The accuracy of all these four

140

custom-made scripts (https://github.com/RichieJu520/Metagenomics.git) has been validated before use to avoid

141

both false-positive and false-negative errors (S1, Supporting Information). One data set named “Shatin AS” was

142

previously used to evaluate different bacterial pathogen identification methods in our previous study.23 Based on

143

the optimized identification method, this data set was compared with seven new data sets using the most

144

updated and comprehensive bacterial pathogen database in the present study to investigate the broad-spectrum

145

profile and fate of bacterial pathogens in STPs. 8

ACS Paragon Plus Environment

Page 8 of 34

Page 9 of 34

146

Environmental Science & Technology

Bacterial Pathogen List

147

The bacterial pathogen database used in the present study includes 538 pathogenic species (Table

148

S3),30much more than the 35 pathogenic species contained in a previous database, 23 and thus could guarantee

149

broad-spectrum monitoring of bacterial pathogens in various environmental samples.

150 151

RESULTS AND DISCUSSION

152

The Occurrence and Abundance of Bacterial Pathogens in Sewage Treatment Plants

153

To distinguish the bacterial pathogen community composition in various types of sludge and wastewater

154

samples via metagenomic approach, it is essential to evaluate the reproducibility of the whole methodology,

155

including DNA extraction, shotgun library construction, Illumina HTS, taxonomic classification using

156

MetaPhlAn. The results shown in Figure S2 and S2 (Supporting Information, SI) revealed that the established

157

metagenomic approach was reliable enough to compare the abundance and variation of bacterial pathogen

158

community among different types of samples. To guarantee the subsequent comparison results, the abundance

159

ratio of 4.0 was selected as the threshold for the significant difference in the following analysis. That is, the

160

abundances for a specific bacterial pathogen occurred in different samples are considered to be significantly

161

different when the corresponding ratio is greater than 4.0 (Figure S2 and S2).

162

Figure 1 demonstrates the diversity and abundance of the overall bacterial pathogens in various samples

163

including influent, effluent, AS (foaming and non-foaming), biofilm and anaerobic digestion sludge. Totally,

164

113 bacterial pathogen species were detected in the above samples. Only three bacterial pathogens with total 9

ACS Paragon Plus Environment

Environmental Science & Technology

165

abundance of 0.11% were detected in SWH ADS which had the lowest bacterial pathogen diversity and

166

abundance. Not surprisingly, influent samples possessed the highest bacterial pathogen diversity, i.e. 91 species,

167

with total abundances of 23.2%. The bacterial pathogen composition in influent could reflect the overall

168

bacterial pathogen profiles in human gut of the population within the STP catchment area to the greatest extent.

169

This will be further discussed in detail in Subsection “PCoA and CA Analysis of Bacterial Pathogen in

170

Sewage Treatment Plants”.

171

Shatin Influent

172

Among the 91 species of bacterial pathogens existed in the influent sample, nine dominant bacterial

173

pathogens, including Megamonas hypermegale, Eubacterium rectale, Escherichia coli, Collinsella aerofaciens,

174

Streptococcus salivarius, Klebsiella pneumonia, Bacteroides vulgatus, Streptococcus suis and Streptococcus

175

bovis were detected with high abundances greater than 0.5%. Megamonas hypermegale, Eubacterium rectale,

176

and Escherichia coli were the three top abundant pathogens, i.e., 4.89%, 4.05% and 2.65% within the total

177

prokaryotic community of influent. To the best of our knowledge, 79 bacterial pathogen species were first found

178

in influents,4,7,9,23,31-33 mainly benefiting from the innovative metagenomic approach which could provide the

179

broad-spectrum scan of bacterial pathogens simultaneously (Table S4). Lee et al. 4 and Shannon et al.9 utilized

180

qPCR/oligonucleotide microarray for bacterial pathogen detection in STPs. Compared to the present study,

181

much less bacterial pathogens, i.e., only 7 and 9 species were detected in raw sewage of a Canada STP,

182

respectively. One of the major reasons for such great difference might be due to that the final detection results in

183

the previous study were limited by the primers or probes which only covered 17 target bacterial pathogen 10

ACS Paragon Plus Environment

Page 10 of 34

Page 11 of 34

Environmental Science & Technology

184

species. Escherichia coli population was considered to dominate the bacterial pathogen community in influent,

185

that is, more abundant than other pathogens by 2~4 orders of magnitude.4 However, except for Escherichia coli,

186

the other 8 species, i.e., Megamonas hypermegale, Eubacterium rectale, Collinsella aerofaciens, Streptococcus

187

salivarius, Klebsiella pneumonia, Bacteroides vulgatus, Streptococcus suis and Streptococcus bovis were also

188

reported as the dominant bacterial pathogen species with comparable abundances in this study. Therefore, the

189

HTS-based metagenomic approach could provide us new insights into broad-spectrum profiles of bacterial

190

pathogens in influents. On the other hand, our previous study using 454 pyrosequencing approach revealed that

191

Arcobacter and Streptococcus with the abundance of 2.88% and 3.30% were the most abundant bacterial

192

pathogen genera in the influent collected from the same STP.10 The abundance of Streptococcus in this study

193

was 3.69%, very similar to that obtained via 454 pyrosequencing method. Nevertheless, Arcobacter only

194

accounted for 0.22% in the present study, one order of magnitude lower than the 454 pyrosequencing result.

195

One of the possible reasons might be due to the inherent amplification biases during PCR step conducted before

196

454 pyrosequencing process.13 In addition, pyrosequencing of the 16S rRNA gene will also generate bias

197

because of the highly variable copy number of the 16S rRNA gene, ranging from 1 to 15. 14

198

Shatin Effluent

199

37 bacterial pathogen species were detected in effluent with the total abundances of 11.8%. Compared with

200

previous studies, 4,9,23,31-33 32 bacterial pathogen species were discovered in effluents of STPs for the first time

201

(Table S5). For example, except for Escherichia coli, other 8 dominant bacterial pathogens, Bacillus cereus,

202

Bacillus thuringiensis, Megamonas hypermegale, Collinsella aerofaciens, Bacillus subtilis, Streptococcus 11

ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 34

203

salivarius, Eubacterium rectale and Bacteroides vulgatus were also found with high abundance more than

204

0.10% in STP effluent. In addition, the number of disease types caused by Escherichia coli only accounts for

205

~10% of the number of total disease types caused by all the pathogens existed in effluents (Table S6). Therefore,

206

the profile of bacterial pathogens in STP effluent obtained in the present study could facilitate us to conduct a

207

comprehensive risk assessment for wastewater quality as well as public health. Undoubtedly, this is the essential

208

anchor for setting appropriate regulatory guidelines of wastewater reclamation and reuse. 2

209 210

Among the 37 bacterial pathogen species totally detected in STP effluent, 3 species, including Escherichia

211

coli, Mycobacterium avium and Shigella boydii belong to bacterial pathogens in water environment and are

212

highly concerned by WHO.22 Additionally, Mycobacterium avium which could cause pulmonary tuberculosis

213

was also selected as the microbial contaminant candidate in Contaminant Candidate List 3 (CCL3) compiled by

214

United

215

http://water.epa.gov/scitech/drinkingwater/dws/ccl/ccl3.cfm). These 3 bacterial pathogen species deserve to be

216

paid more attention in effluents because more and more wastewater reuse projects have been carried out around

217

the world and these pathogens are highly concerned in water quality control.34-36

States

Environmental

Protection

Agency

(US

EPA,

218

Shatin AS

219

Compared to influent and effluent, much less bacterial pathogen species and total abundances were found

220

in “Shatin-AS” and

221

(Arcobacter butzleri, Collinsella aerofaciens, Eubacterium rectale, Gordonia bronchialis, Streptococcus

“Shatin-AS Foaming” (Figure 1). To the best of our knowledge, all the pathogen species

12

ACS Paragon Plus Environment

Page 13 of 34

Environmental Science & Technology

222

salivarius, Vibrio furnissii and Vibrio mimicus) detected in “Shatin-AS” or “Shatin-AS Foaming” were first

223

reported in AS system of STPs. 9,10, 23, 32 Shannon et al.

224

and Escherichia coli as well as Klebsiella pneumonia were found in AS of a Canada STP. Although Escherichia

225

coli and Klebsiella pneumonia were also involved in our pathogen database used in this study, they were not

226

detected in the activated sludge samples. The possible reason might be due to the lower detection limit of qPCR

227

compared to the metagenomic approach at current sequencing depth. That is, the metagenomic approach is less

228

sensitive than qPCR approach and some rare pathogens could not be covered at the sequencing depth applied in

229

this study. On the other hand, the previous study found no pathogens at species level in Shatin AS although

230

metagenomic approach was also utilized.23 The major reason should be due to the different sizes of pathogen

231

databases, i.e., 538 pathogenic species in the database used in this study, much more comprehensive than the 35

232

pathogenic species in the database used in the previous study.23

9

applied qPCR to detect 13 bacterial pathogens in AS

233

Specially, three pathogens, i.e., Arcobacter butzleri, Gordonia bronchialis and Vibrio furnissii, only existed

234

in the “Shatin-AS Foaming” which was collected from the upper foaming layer of AS. This is consistent with

235

previous studies which indicated that Gordonia spp. was responsible for AS foaming.37,38 Nevertheless, up to

236

now, there are no reports about the existence of Arcobacter butzleri or Vibrio furnissii in foaming AS and thus

237

their roles during AS foaming were not well understood currently. More pathogen species and higher total

238

abundance of detected pathogens were found in “Shatin-AS Foaming” than “Shatin AS”, which was obtained

239

from the bulk mixed liquor. This implies the foaming layer of AS might impose more threat to workers in STP

240

and the people in the surrounding areas of STPs because pathogens existed in foaming layer of AS are easily 13

ACS Paragon Plus Environment

Environmental Science & Technology

241

transferred into air and cause infections accordingly.39

242

Stanley AS and Biofilm

243

27 and 24 bacterial pathogen species were detected in “Stanley AS” and “Stanley BF” with the total

244

abundances of 0.64% and 0.56%, respectively. Although the diversity and total abundance of the overall

245

bacterial pathogens in these two samples were comparable with each other, the occurrence of some specific

246

pathogens exhibited significant difference between AS and BF samples (Table S7). Among the 9 bacterial

247

pathogen species which occurred with higher abundances in “Stanley BF”, Mycobacterium abscessus,

248

Mycobacterium avium, Mycobacterium bovis, Mycobacterium kansasii, Mycobacterium marinum, Nocardia

249

farcinica and Ralstonia pickettii were commonly found to be able to form biofilm under various

250

environments.40-46 It should be noted that most studies on Mycobacterium pathogens are specific to drinking

251

water or hospital water distribution systems due to their close relationship with the public health.42,44,48 However,

252

the related studies for activated sludge or biofilm of wastewater system are very limited currently.47 Therefore,

253

the results of this study expanded the horizons of Mycobacterium pathogens into different environments. It

254

should be emphasized that, compared to planktonic bacterial pathogens, the pathogens existed in the biofilm

255

might impose more threat to human health because they could demonstrate stronger and more persistent

256

resistance to disinfectant when they disperse into the effluent via the detached biofilm.48-52 Therefore, these

257

pathogens should deserve special attention for reclamation and reuse of wastewater which are treated using

258

attached growth biological treatment processes.

259 14

ACS Paragon Plus Environment

Page 14 of 34

Page 15 of 34

Environmental Science & Technology

260

Acinetobacter johnsonii, one of the ten pathogens which occurred with much higher abundance in “Stanley

261

AS” was regarded as an important bridging microorganism for sludge floc formation in AS.53-55 This pathogen

262

showed a high affinity for many of settled sewage bacteria and activated sludge flocs.55 The extent and pattern

263

of coaggregation as well as the aggregate size exhibited good correlation with cellular characteristics of the

264

coaggregating partners.53 FISH result indicates that the abundance of Acinetobacter spp. in AS of a Japan STP is

265

2.2%,55 much higher than that of Acinetobacter johnsonii (0.19%) in the present study. The possible reasons for

266

such differences might be due to different samples and different detection methods (FISH vs. Metagenomic

267

approach) applied in these two studies. In addition, the abundance of 2.2% represents the total abundance for all

268

the Acinetobacter species, rather than Acinetobacter johnsonii only.55

269

Shatin ADS and SWH ADS

270

Six and three bacterial pathogen species were detected in “Shatin-ADS” and “SWH-ADS” with the total

271

abundances of 1.23% and 0.11%, respectively (Figure 1). This is comparable with the previous study which

272

reported that the bacterial pathogen abundance ranged between 0.55% and 1.57% for the mesophilic ADS.

273

Bibby et al.5 used 454 pyrosequencing of the 16S rRNA gene to reveal bacterial pathogen diversity and

274

abundances in mesophilic and thermophilic anaerobic digestion sludge. Due to the limitation of PCR based 16S

275

pyrosequencing and RDP taxonomic assignment approach, bacterial pathogens existed in digestion sludge could

276

only be identified at genus level rather than species level. The dominant pathogen genera were Clostridium and

277

Mycobacterium and the total abundances of all the pathogens in mesophilic and thermophilic anaerobic

278

digestion sludge were 0.08% and 0.022%, respectively.5 Very similar pathogen species composition and total 15

ACS Paragon Plus Environment

12

Environmental Science & Technology

279

abundances were found between “SWH-ADS” and the mesophilic anaerobic digestion sludge reported by Bibby

280

et al..5 However, the total abundance of pathogen in “Shatin-ADS” exhibited one order of magnitude higher than

281

that of “SWH-ADS”. Additionally, the pathogen compositions of “Shatin-ADS” and “SWH-ADS” also

282

exhibited remarkable difference (Table S8). The possible reason is due to that “Shatin-ADS” mainly derived

283

from AS treating saline sewage (salinity of ~1.2%) while “SWH-ADS” primarily originated from AS treating

284

freshwater sewage. These total different environmental conditions would probably lead to the distinct difference

285

of pathogen species and abundance in these ADS samples.

286 287

Among the eight bacterial pathogen species occurred in “Shatin-ADS” and “SWH-ADS”, Bacteroides vulgatus,

288

Collinsella aerofaciens, Eubacterium rectale, Streptococcus salivarius, and Streptococcus suis were discovered

289

in ADS of STP for the first time.5,56 Arcobacter butzleri was frequently found in ADS collected from an Italy

290

STP in different seasons, and the mean abundance is 7649 MPN/g dry ADS.5 Collinsella aerofaciens and

291

Eubacterium rectale are commonly found in the human gut.57-59 It was reported that Collinsella aerofaciens, a

292

high G+C, Gram-positive member of the Actinobacteria exhibited salt tolerance due to the presence of galE,

293

murB and mazG genes. This is consistent with the phenomenon that higher abundance of Collinsella

294

aerofaciens was observed in ADS collected from Shatin STP which treats saline sewage. Some indicator

295

microorganisms, such as Escherichia coli and Clostridium perfringens were not detected in ADS. However, it

296

does not mean the absence of potential risk for the ADS due to the occurrence of other pathogens (Table S8).

297

Instead, this further indicates the limitation of indicator microorganism method which might underestimate, 16

ACS Paragon Plus Environment

Page 16 of 34

Page 17 of 34

298

Environmental Science & Technology

sometimes even misestimate the potentially infectious risk for human health.

299

Although the final disposal for ADS in Hong Kong is landfill, application of ADS to agriculture is

300

widespread in the USA and Europe. Considering the strong survival ability of these bacterial pathogens56 and

301

the various diseases caused by them (Table S9), it is believed that the land application of ADS could result in

302

serious pollution in soil, farm produce and underground water. Even for landfill treatment, there is still potential

303

risk to transfer these pathogens into underground water and thus cause the spread of multiple diseases due to the

304

leachate leakage. Therefore, the full understanding of the bacterial pathogens profile in ADS is the crucial step

305

for the safe disposal and reuse of ADS.

306 307

Emerging/Re-emerging Bacterial Pathogens in Sewage Treatment Plants

308

Emerging/re-emerging pathogen species are notable for their taxonomic, biological, epidemiological and

309

geographical diversity. Emerging bacterial pathogen is defined as the bacterial infectious agents whose

310

incidences are increasing following their first introduction into the new host population. Re-emerging bacterial

311

pathogens are bacterial infectious agents whose incidences are increasing in the existing host population as a

312

result of long-term changes in their underlying epidemiology.30 It was reported that a significant number of

313

emerging and re-emerging bacterial pathogens had been recognized and water transmission had caused

314

outbreaks including Escherichia coli O157:H7, Helicobacter pylori, Legionella pneumophila, Campylobacter

315

jejuni as well as Mycobacterium avium complex.3 17 and 8 emerging/re-emerging bacterial pathogens were

316

found in influent/effluent and AS/ADS, respectively (Figure 2). The abundances of these emerging/re-emerging 17

ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 34

317

pathogens were 3.84% in STP influent and 0.33% in STP effluent. Emerging/re-emerging pathogens accounted

318

for 16.6% and 2.79% of the total bacterial pathogens in influent and effluent, respectively. The decrease of

319

percentages of the emerging/re-emerging pathogens in the total bacterial pathogens after the secondary

320

biological treatment implies that, compared to other non-emerging/re-emerging pathogens, some

321

emerging/re-emerging pathogens are easier to be removed in AS treatment process. The detailed removal

322

efficiency for the specific pathogen will be discussed in Subsection “Removal of Bacterial Pathogen in

323

Activated Sludge Process”.

324 325

As shown in Figure 2a, Escherichia coli and Klebsiella pneumoniae were most abundant, accounting for

326

88.9~92.0% of the total emerging/re-emerging pathogens in effluent and influent. Except for Escherichia coli

327

and

328

emerging/re-emerging pathogens) in STP effluent. Escherichia coli is usually used as an indicator

329

microorganism to estimate the presence of pathogenic microorganisms in wastewater. Klebsiella pneumoniae

330

was previously detected in the influents of a constructed soil filter in India and the raw sewage of a STP with

331

bioreactor “BIO-PAK” in Poland.31,39 As shown in Figure 2b, emerging/re-emerging pathogens accounted for

332

14.8%, 56.5% and 82.7 % of the total bacterial pathogens in “Stanley-AS”, “Stanley-BF” and “SWH-ADS”,

333

respectively. The abundances of these emerging/re-emerging pathogens ranged from 0.092% to 0.32%. No

334

emerging/re-emerging bacterial pathogens were detected in “Shatin-AS foaming”, “Shatin-AS” and

335

“Shatin-ADS”. Mycobacterium spp. were widely found in AS, BF and ADS samples. In “Stanley-BF”,

Klebsiella

pneumoniae,

Mycobacterium

ulcerans

was

also

18

ACS Paragon Plus Environment

abundant

(7.6%

of

the

total

Page 19 of 34

Environmental Science & Technology

336

Mycobacterium bovis is the most abundant emerging/re-emerging pathogen followed by Bordetella pertussis.

337

Although Mycobacterium bovis is able to form biofilm,44 there was no report on its occurrence in STP

338

previously. Bordetella pertussis, Mycobacterium ulcerans and Aeromonas hydrophila which were first reported

339

in AS in this study are the dominant emerging/re-emerging pathogen species in “Stanley-AS”. Mycobacterium

340

ulcerans and Clostridium botulinum made up the total emerging/re-emerging pathogens in “SWH-ADS”.

341

Clostridium botulinum, an anaerobic spore-former, was detected in sludge samples from biogas plants before

342

digestion, but disappeared after digestion.58 It implies that digestion is an effective process to remove

343

Clostridium botulinum. However, we could not calculate the removal efficiency of Clostridium botulinum in the

344

present study since we have no corresponding abundance data in “SWH-AS”.

345 346 347 348

Removal of Bacterial Pathogen in Activated Sludge Process

Removal efficiency of each specific bacterial pathogen was calculated based on Equation (1):

Re moval efficiency (%) = (1 −

VSSeffluent × Abundance pathogen,eff VSSinfluent × Abundance pathogen,inf

) × 100

(1)

349

Where VSSinfluent and VSSeffluent represent the volatile suspended solids (VSS) in influent and effluent,

350

respectively. Abundance pathogen,eff and Abundance pathogen,inf are the abundance of one specific bacterial pathogen

351

in effluent and influent, respectively. From the point view of engineering, we used Equation (1) to calculate the

352

removal efficiencies due to the fact that increasing volatile suspended solids removal also increases bacterial

353

removal as bacteria are often solid-associated. Additionally, we assumed that the pathogenic bacteria represent

354

only some percentage of total bacteria in VSS. VSS of influent and effluent during the sampling period were 19

ACS Paragon Plus Environment

Environmental Science & Technology

355

summarized in Table S10.

356 357

The removal efficiency of total bacterial pathogens was 98.0%. This suggests that AS treatment process is

358

efficient and effective to reduce the potential health risk to public health caused by bacterial pathogens. Based

359

on the removal efficiency results, bacterial pathogens in influents were divided into five groups (Figure 3).

360

Group I referred to that these 3 pathogen species could not be removed efficiently and their remove efficiencies

361

were lower than 1-log. Group II indicated that the 12 pathogen species could be removed with efficiencies

362

ranging from 1-log to 2-log. For Group III, the 14 pathogens were removed by 2-log to 3-log. Group IV

363

included 4 pathogens with removal efficiencies of 3-log ~ 4-log. For Group V, the 58 pathogens could be

364

reduced by >4-log during AS treatment process. As mentioned in Subsection “Emerging/Re-emerging

365

Bacterial Pathogens in Sewage Treatment Plants”, the decrease of percentages of the emerging/re-emerging

366

pathogens in the total bacterial pathogens after the secondary biological treatment implies that some

367

emerging/re-emerging pathogens are easier to be removed in AS treatment process. As shown in Figure 3, 13

368

out of 15 emerging/re-emerging pathogen species belong to Group III ~ Group V with high removal efficiencies

369

ranging from >99% to >99.99%. Additionally, the detected bacterial pathogen listed in CCL3 (Shigella sonnei)

370

and transmitted through drinking-water concerned by WHO (Shigella boydii, Escherichia coli, Vibrio cholera,

371

Shigella sonnei, Shigella flexneri and Shigella dysenteriae) were well removed by the activated sludge process

372

(>99% ~ >99.99%). Since only limited bacterial pathogens have been detected in previous studies, the

373

comparison of removal efficiencies between the present study and previous references could not cover all the 20

ACS Paragon Plus Environment

Page 20 of 34

Page 21 of 34

Environmental Science & Technology

374

pathogens illustrated in Figure 3. Shannon et al.9 used qPCR to detect selected bacterial pathogens during

375

municipal wastewater treatment and found that Escherichia coli, Klebsiella pneumoniae, Clostridium

376

perfringens as well as Enterococcus faecalis were all efficiently removed with high removal efficiencies

377

(99.93~99.99%) calculated based on the gene copy numbers per 100 mL wastewater.9 In this study, Escherichia

378

coli and Klebsiella pneumoniae exhibited similar removal pattern, i.e., 99.34%~99.42% of them could be

379

removed during AS treatment process. For Clostridium perfringens and Enterococcus faecalis, the removal

380

efficiencies were as high as >99.9% and >99.99%, respectively in the present study. Wen et al. 2 once detected

381

Clostridium perfringens in Bolivar STP utilizing secondary AS treatment process in Australia and the

382

corresponding removal efficiency was 94.22%, much lower than those reported by Shannon et al.9 and our study.

383

The possible reasons might be as follows. Clostridium perfringens is an anaerobe but is able to survive under

384

aerobic stage as spores. 6 h of anoxic stirring for nitrogen removal would cause the detachment of Clostridium

385

perfringens spores from the sludge flocs and consequently reducing the removal efficiency.2 Aeromonas

386

hydrophila was only detected in raw wastewater and was undetectable thereafter in Shannon et al.’s study.9

387

Similarly, the removal efficiency of Aeromonas hydrophila in this study was 99.999%, also consistent with the

388

results of 99.98% reported by Lee et al..4

389 390

PCoA and CA Analysis of Bacterial Pathogen in Sewage Treatment Plants

391

Principal coordinate analysis (PCoA) and cluster analysis (CA) were conducted to compare the overall

392

bacterial pathogen community compositions for different types of samples collected from STPs. PCoA results 21

ACS Paragon Plus Environment

Environmental Science & Technology

393

revealed that bacterial pathogen communities in these various samples could be clustered into five groups

394

(Figure 4): (1) Group I influent sample obtained from Shatin STP as well as two human gut samples published

395

in MG-RAST (http://metagenomics.anl.gov/, Table S11); (2) Group II contains effluent sample collected from

396

Shatin STP; (3) Group III is ADS from Shek Wu Hui SPT which treats freshwater sewage containing

397

slaughterhouse wastewater; (4) Group IV contains AS and BF samples from Stanley STP which is located inside

398

a cave and treats freshwater sewage; (5) Group V contains AS, foaming AS and ADS samples from Shatin STP

399

which treats saline sewage (salinity of ~1.2%). This grouping pattern was reasonable and in consistent with the

400

CA results shown in Figure S3. That is, the clustering has strong relatedness with both the sample types

401

(influent vs. effluent vs. AS/ADS samples) and the STP characteristics (freshwater STP or saline STP). It should

402

be pointed out that SWH ADS and Shatin ADS were clustered into two separate groups although they belong to

403

the same sample type, i.e., anaerobic digestion sludge. In addition, “SWH-ADS” clustered more closely to

404

“Stanley-AS” and “Stanley-BF” while “Shatin-ADS” clustered more tightly with “Shatin-AS” and “Shatin-AS

405

foaming”. Similarly, “Stanley-AS” and “Shatin-AS” were also separated into Group IV and Group V,

406

respectively. As mentioned above, both Stanley and Shek Wu Hui STPs treat freshwater sewage while Shatin

407

STP treats saline sewage. Therefore, it might imply that STP characteristic (salinity effect) is the main

408

determinant for this grouping pattern.

409 410

We hypothesized that there would be remarkable similarities of bacterial pathogen compositions between

411

influent samples and human gut samples. In other words, the bacterial pathogen profiles in influents could 22

ACS Paragon Plus Environment

Page 22 of 34

Page 23 of 34

Environmental Science & Technology

412

reflect the average bacterial pathogen communities of the urban resident guts within the STP catchment area to a

413

great extent. The grouping pattern demonstrated by the PCoA and CA results effectively confirmed this

414

assumption. Compared to other samples, the influent sample indeed clustered much more tightly with two

415

human gut samples (Figure 4). This result might be driven primarily by the shared presence of 14 dominant

416

pathogens in both influent and human gut samples, i.e., Eubacterium rectale, Bacteroides vulgatus, Collinsella

417

aerofaciens, Bacteroides stercoris, Escherichia coli, Streptococcus salivarius, Bacteroides ovatus, Sutterella

418

wadsworthensis, Bacteroides caccae, Bacteroides uniformis, Bilophila wadsworthia, Eggerthella lenta,

419

Bacteroides thetaiotaomicron and Bacteroides fragilis (Table S12). Assessing average bacterial community

420

structure in urban residents from catchment areas has also been studied in previous works using 16S rRNA.

421

The results indicated that influent sewage bacterial profile reflected the human microbiome and human gut

422

bacterial community was the dominant force shaping influent sewage bacterial profile. These findings were

423

consistent with the results obtained in our study.

11

424 425

Advantage and Disadvantage of Metagenomic Approach

426

The breakthrough of next-generation sequencing technology launched a new era for the broad-spectrum

427

monitoring and rapid detection of pathogens. To the best of our knowledge, this is the first study to conduct the

428

broad-spectrum monitoring of bacterial pathogens in various environmental samples simultaneously using

429

metagenomic approach combined with HTS technology. The full understanding of the bacterial pathogens

430

profile in effluent and ADS could facilitate us to conduct a comprehensive risk assessment for public health 23

ACS Paragon Plus Environment

Environmental Science & Technology

431

during the wastewater reclamation/reuse as well as application of ADS to agriculture. Compared to other

432

molecular biology detection method of bacterial pathogen (PCR/qPCR or 454 pyrosequencing of the 16S rRNA

433

gene), the metagenomic approach used in this study not only provides us a broad picture of the occurrence of

434

bacterial pathogens, but also shows the following advantages. Firstly, it is able to accurately analyze the

435

bacterial pathogen community and quantify the abundance without PCR bias and copy number variability.11,31

436

Secondly, it can investigate the bacterial pathogen composition at a lower taxonomic resolution, i.e., species

437

level. This is especially crucial to characterize and quantify pathogens since the investigation of pathogens at

438

genus level usually leads to some inaccurate estimate.15 Finally but more importantly, it could standardize the

439

bacterial pathogen detection method in different laboratories and realize the direct comparison of results among

440

different samples and studies on a united platform.

441 442

However, the metagenomic approach used in this study still has its limitations. Firstly, MetaPhlAn tool which

443

is used for pathogen annotation only applies to bacterial pathogens and thus could not cover the other pathogen

444

types including viruses/prions, fungi, protozoa and helminthes.15 These four pathogen types account for as high

445

as 61.8% of the total reported human pathogen species.26 In addition, some novel bacterial pathogens which was

446

excluded in the database of MetaPhlAn could not be annotated even they actually existed in the samples (Table

447

S13). Secondly, metagenomic approach might be less sensitive than qPCR approach at the current sequencing

448

depth and rare pathogens might not be detected, especially for complicated environmental samples, such as AS

449

and ADS. Thirdly, similar to the other molecular biology techniques, this metagenomic approach might 24

ACS Paragon Plus Environment

Page 24 of 34

Page 25 of 34

Environmental Science & Technology

450

overestimate the abundance of metabolically active bacterial pathogens since it is based on the total DNA

451

extracted from both metabolically active bacteria and dead bacterial cells. Moreover, it should be pointed out

452

that the strains within the identified pathogen species vary greatly in pathogenicity and therefore the bacterial

453

species identified in this study should be “potential pathogens”. Finally, HTS for metagenomics is still costly at

454

present and cannot be applied to routine analysis in STPs, although the cost has been greatly decreased in recent

455

years.61

456 457

Acknowledgment

458

This work was supported by the Research Grants Council of Hong Kong (HKU7201/11E). Dr. Lin Cai thank

459

The University of Hong Kong for the postdoctoral fellowships. Mr. Feng Ju thanks The University of Hong

460

Kong for the postgraduate studentship. The authors also wish to thank Dr. Feng Guo for his valuable

461

suggestions and discussion.

462

Supporting Information Available

463

Supplementary data includes S1, S2, Fig. S1- Fig. S3 and Table S1-Table S13. This information is available

464

free of charge via the Internet at http://pubs.acs.org.

465

References

466

(1) Prüss, A.; Kay, D.; Fewtrell, L.; Bartram, J. Estimating the burden of disease from water, sanitation, and

467

hygiene at a global level. Environ. Health Persp. 2002,110 (5), 537-542.

25

ACS Paragon Plus Environment

Environmental Science & Technology

468

(2) Wen, Q. X.; Tutuka, C.; Keegan, A.; Jin, B. Fate of pathogenic microorganisms and indicators in secondary

469

activated sludge wastewater treatment plants. J. Environ. Manage. 2009, 90 (3), 1442-1447.

470

(3) Aw, T. G.; Rose, J. B. Detection of pathogens in water: from phylochips to qPCR to pyrosequencing. Curr.

471

Opin. Biotech. 2012, 23 (3), 422-430.

472

(4) Lee, D. Y.; Shannon, K.; Beaudette, L. A. Detection of bacterial pathogens oligonucleotide microarray in

473

municipal wastewater using and real-time quantitative PCR. J. Microbiol. Meth. 2006, 65 (3), 453-467.

474

(5) Bibby, K.; Viau, E.; Peccia, J. Pyrosequencing of the 16S rRNA gene to reveal bacterial pathogen diversity

475

in biosolids. Water Res. 2010, 44 (14), 4252-4260.

476

(6) Rose, J. B.; Farrah, S. R. Microbial health risks and water quality. In Urban Water Security: Managing Risks.

477

Edited by Jimenez B, Rose J. UNESCO, 2008, 21-52.

478 479

(7) Toze, S. PCR and the detection of microbial pathogens in water and wastewater. Water Res. 1999, 33 (17),

480

(8) Ulrich, H.; Klaus, D.; Irmgard, F.; Annette, H.; Juan, L. P.; Regine, S. Microbiological investigations for

481

sanitary assessment of wastewater treated in constructed wetlands. Water Res. 2005, 39, (20), 4849-4858.

482

(9) Shannon, K. E.; Lee, D. Y.; Trevors, J. T.; Beaudette, L. A. Application of real-time quantitative PCR for the

483 484

detection of selected bacterial pathogens during municipal wastewater treatment. Sci. Total Environ. 2007, 382

485

(10) Ye, L.; Zhang, T. Pathogenic bacteria in sewage treatment plants as revealed by 454 pyrosequencing.

486 487

Environ. Sci. Technol. 2011, 45 (17), 7173-7179.

488

Appl Microbiol Biotechnol 2014, 98,3317-3326.

489

(12) Stiborova1, H.; Wolfram, J.; Demnerova K.; Macek T.; Uhlik, O. Bacterial community structure in treated

490 491

sewage sludge with mesophilic and thermophilic anaerobic digestion. Folia Microbiol 2015, DOI

492 493

(13) Haft, D. H.; Tovchigrechko, A. High-speed microbial community profiling. Nat. Methods 2012, 9 (8),

494

(14) Liu, B.; Gibbons, T.; Ghodsi, M.; Treangen, T.; Pop, M. Accurate and fast estimation of taxonomic profiles

495

from metagenomic shotgun sequences. BMC Genomics 2011,12:S4

496

(15) Segata, N.; Waldron, L.; Ballarini, A.; Narasimhan, V.; Jousson, O.; Huttenhower, C. Metagenomic

497

microbial community profiling using unique clade-specific marker genes. Nat. Methods 2012, 9 (8), 811-814.

3545-3556.

(1), 121-129.

(11) Cai, L.; Ju, F.; Zhang, T. Tracking human sewage microbiome in a municipal wastewater treatment plant.

10.1007/s12223-015-0396-9

793-794.

498 26

ACS Paragon Plus Environment

Page 26 of 34

Page 27 of 34

Environmental Science & Technology

499 500

phenotypic and genotypic methods for detection of enterotoxigenic Escherichia coli toxins and colonization

501

factors. J. Clin. Microbiol. 2007, 45 (10), 3295−3301.

502

(17) Cao, C.; Jiang, W. J.; Wang, B. Y.; Fang, J. H.; Lang, J. D.; Tian, G.; Jiang, J. K.; Zhu, T. F. Inhalable

503 504

microorganisms in Beijing’s PM2.5 and PM10 pollutants during a severe smog event. Environ. Sci. Technol. 2014,

505

(18) Djikeng, A.; Kuzmickas, R.; Anderson, N. G.; Spiro, D. J. Metagenomic analysis of RNA viruses in a fresh

506

water lake. Plos One 2009. 4 (9).

507

(19) Rosario, K.; Nilsson, C.; Lim, Y. W.; Ruan, Y. J.; Breitbart, M. Metagenomic analysis of viruses in

508

reclaimed water. Environ. Microbiol. 2009, 11, 2806-2820.

509

(20) Bibby, K.; Viau, E.; Peccia, J. Viral metagenome analysis to guide human pathogen monitoring in

510

environmental samples. Lett. Appl. Microbiol. 2011, 52 (4), 386-392.

511

(21) Bibby, K.; Peccia, J. Identification of viral pathogen diversity in sewage Sludge by metagenome analysis.

512

Environ. Sci. Technol. 2013, 47 (4), 1945-1951.

513 514

(22) WHO. Recommendations. Geneva. Guidelines for Drinking-water Quality: Microbial aspects, fourth ed.;

515

(23) Cai, L.; Zhang, T. Detecting human bacterial pathogens in wastewater treatment plants by a

516

high-throughput shotgun sequencing technique. Environ. Sci. Technol. 2013, 47 (10), 5433–5441.

517 518

(24) Qin, J. J.; Li, R. Q.; Raes, J.; Arumugam, M.; Burgdorf, K. S.; Manichanh, C.; Nielsen, T.; Pons, N.;

519 520 521 522

Q.; Zheng, H. S.; Xie, Y. L.; Tap, J.; Lepage, P.; Bertalan, M.; Batto, J. M.; Hansen, T.; Le Paslier, D.;

523 524

Meta, H. I. T. C. A human gut microbial gene catalogue established by metagenomic sequencing. Nature 2010,

525

(25) Gomez-Alvarez, V.; Teal, T. K.; Schmidt, T. M. Systematic artifacts in metagenomes from complex

526

microbial communities. ISME J. 2009, 3 (11), 1314-1317.

527

(26) Brown, J. M.; LaBarre, B. A.; Hewson, I. Characterization of Trichodesmium-associated viral communities

528

in the eastern Gulf of Mexico. FEMS Microbiol. Ecol. 2013, 84 (3), 603-613.

529

(27) Lipson, D. A.; Haggerty, J. M.; Srinivas, A.; Raab, T. K.; Sathe, S.; Dinsdale, E. A. Metagenomic insights

530

into anaerobic metabolism along an Arctic peat soil profile. Plos One 2013, 8 (5).

(16) Sjöling, Å.; Wiklund, G.; Savarino, S. J.; Cohen, D. I.; Svennerholm, A. M. Comparative analyses of

48 (3), 1499–1507.

http://whqlibdoc.who.int/publications/2011/9789241548151_eng.pdf (accessed 19.08.13.).

Levenez, F.; Yamada, T.; Mende, D. R.; Li, J. H.; Xu, J. M.; Li, S. C.; Li, D. F.; Cao, J. J.; Wang, B.; Liang, H.

Linneberg, A.; Nielsen, H. B.; Pelletier, E.; Renault, P.; Sicheritz-Ponten, T.; Turner, K.; Zhu, H. M.; Yu, C.; Li, S. T.; Jian, M.; Zhou, Y.; Li, Y. R.; Zhang, X. Q.; Li, S. G.; Qin, N.; Yang, H. M.; Wang, J.; Brunak, S.; Dore, J.; Guarner, F.; Kristiansen, K.; Pedersen, O.; Parkhill, J.; Weissenbach, J.; Bork, P.; Ehrlich, S. D.; Wang, J.;

464 (7285), 59-70.

27

ACS Paragon Plus Environment

Environmental Science & Technology

531 532

(28) Solonenko, S. A.; Ignacio-Espinoza, J. C.; Alberti, A.; Cruaud, C.; Hallam, S.; Konstantinidis, K.; Tyson,

533

BMC Genomics 2013,14:320

534

(29) Yang, Y.; Li, B.; Ju, F.; Zhang, T. Exploring variation of antibiotic resistance genes in activated sludge

535

over a four-year period through a metagenomic approach. Environ. Sci. Technol. 2013, 47(18), 10197–10205.

536 537 538 539

(30) Woolhouse, M.; Gowtage-Sequeria, S.; Evans, B. T16: Quantitative analysis of the characteristics of

540

(31) Kadam, A. M.; Oza, G. H.; Nemade, P. D.; Shankar, H. S. Pathogen removal from municipal wastewater in

541

constructed soil filter. Ecol. Eng. 2008, 33 (1), 37-44.

542

(32) Wéry, N.; Lhoutellier, C.; Ducray, F.; Delgenès, J. P.; Godon, J. J. Behaviour of pathogenic and indicator

543

bacteria during urban wastewater treatment and sludge composting, as revealed by quantitative PCR. Water Res.

544

2008, 42 (1-2), 53-62.

545 546

(33) Eddabra, R.; Moussaoui, W.; Prevost, G.; Delalande, F.; Van Dorsselaer, A.; Meunier, O.; Scheftel, J. M.;

547

World J. Microb. Biot. 2011, 27 (5), 1099-1108.

548

(34) Kajenthira, A.; Siddiqi, A.; Anadon, L. D. A new case for promoting wastewater reuse in Saudi Arabia:

549

Bringing energy into the water equation. J. Environ. Manage. 2012, 102, 184-192.

550

(35) Xi, X.; Poh, K. L. Using system dynamics for sustainable water resources management in Singapore. In

551 552

2013 Conference on Systems Engineering Research, Paredis, C. J. J.; Bishop, C.; Bodner, D.; Eds. 2013, Vol.

553 554

(36) Woods, G. J.; Kang, D.; Quintanar, D. R.; Curley, E. F.; Davis, S. E.; Lansey, K. E.; Arnold, R. G.

555

Pl -ASCE 2013,139 (3), 313-324.

556

(37) Soddell, J. A.; Stainsby, F. M.; Eales, K. L.; Seviour, R. J.; Goodfellow, M. Gordonia defluvii sp nov.; an

557

actinomycete isolated from activated sludge foam. Int. J. Syst. Evol. Micr. 2006, 56, 2265-2269.

558

(38) De los Reyes, F.L. Foaming. In: Microbial Ecology of Activated Sludge. IWA Pulishing, London, UK,

559

2010, pp. 215-258.

560 561

(39) Korzeniewska, E.; Filipkowska, Z.; Gotkowska-Plachta, A.; Janczukowicz, W.; Dixon, B.; Czulowska, M.

562

2009, 43 (11), 2841-2851.

G.; Wincker, P.; Sullivan, M. B. Sequencing platform and library preparation choices impact viral metagenomes.

emerging and re-emerging human pathogens. Centre for Infectious Diseases, University of Edinburgh. http://webarchive.nationalarchives.gov.uk/20121212135622/http://www.bis.gov.uk/assets/foresight/docs/infecti ous-diseases/t16.pdf, (accessed 25.07.15.).

Mimouni, R. Occurrence of Vibrio cholerae non-O1 in three wastewater treatment plants in Agadir (Morocco).

16, pp 157-166.

Centralized versus decentralized wastewater reclamation in the Houghton area of Tucson, Arizona. J.Water Res.

Determination of emitted airborne microorganisms from a BIO-PAK wastewater treatment plant. Water Res.

28

ACS Paragon Plus Environment

Page 28 of 34

Page 29 of 34

Environmental Science & Technology

563

(40) Vaerewijck, M. J. M.; Huys, G.; Palomino, J. C.; Swings, J.; Portaels, F. Mycobacteria in drinking water

564

distribution systems: ecology and significance for human health. FEMS Microbiol. Rev. 2005, 29 (5), 911-934.

565

(41) Alavi, M. R.; Shukla, H. D.; Whitaker, B.; Arnold, J.; Shahamat, M. Attachment and biofilm formation of

566 567

Mycobacterium marinum on a hydrophobic surface at the air interface. World J. Microb. Biot. 2007, 23 (1),

568 569

(42) Williams, M. M.; Yakrus, M. A.; Arduino, M. J.; Cooksey, R. C.; Crane, C. B.; Banerjee, S. N.; Hilborn, E.

570

Mycobacteria. Appl. Environ. Microb. 2009, 75 (7), 2091-2098.

571

(43) Cheng, W. C.; Kuo, M. C. T. The effect of toluene-substrate concentration on biological clogging and

572

sloughing in porous media. Ground Water Monit. R. 2010, 30, (4), 73-80.

573

(44) Ishida, S.; Arai, M.; Niikawa, H.; Kobayashi, M. Inhibitory effect of cyclic trihydroxamate siderophore,

574 575

desferrioxamine E, on the biofilm formation of Mycobacterium species. Biol. Pharm. Bull. 2011, 34 (6),

576

(45) Mijnendonckx, K.; Provoost, A.; Ott, C. M.; Venkateswaran, K.; Mahillon, J.; Leys, N.; Van Houdt, R.

577

Characterization of the survival ability of Cupriavidus metallidurans and Ralstonia pickettii from space-related

578

environments. Microb. Ecol. 2013, 65 (2), 347-360.

579 580 581

(46) Ovrutsky, A. R.; Chan, E. D.; Kartalija, M.; Bai, X. Y.; Jackson, M.; Gibbs, S.; Falkinham, J. O.; Iseman,

Mycobacteria in hospital water networks, and preferential growth of Mycobacterium avium in Acanthamoeba

582

lenticulata. Appl. Environ. Microb. 2013, 79 (10), 3185-3192.

583 584

(47) Radomski, N.; Betelli, L.; Moilleron, R.; Haenn, S.; Moulin, L.; Cambau, E.; Rocher, V.; Goncalves, A.;

585

coil and Enterococci. Environ. Sci. Technol. 2011, 45 (12), 5380-5386.

586

(48) Costerton, J. W.; Stewart, P. S.; Greenberg, E. P. Bacterial biofilms: A common cause of persistent

587

infections. Science 1999, 284 (5418), 1318-1322.

588

(49) Tsai, Y. P. Interaction of chlorine concentration and shear stress on chlorine consumption, biofilm growth

589

rate and particle number. Bioresource Technol. 2006, 97 (15), 1912-1919.

590

(50) Chowdhury, N.; Nakhla, G.; Zhu, J. Load maximization of a liquid-solid circulating fluidized bed

591

bioreactor for nitrogen removal from synthetic municipal wastewater. Chemosphere 2008, 71, (5), 807-815.

592

(51) Lackner, S.; Holmberg, M.; Terada, A.; Kingshott, P.; Smets, B. F. Enhancing the formation and shear

593

resistance of nitrifying biofilms on membranes by surface modification. Water Res. 2009, 43 (14), 3469-3478.

93-101.

D.; Donlan, R. M. Structural analysis of biofilm formation by rapidly and slowly growing nontuberculous

917-920.

M. D.; Reynolds, P. R.; McDonnell, G.; Thomas, V. Cooccurrence of free-living amoebae and nontuberculous

Lucas, F. S. Mycobacterium behavior in wastewater treatment plant, a bacterial model distinct from Escherichia

29

ACS Paragon Plus Environment

Environmental Science & Technology

594

(52) Epstein, A. K.; Pokroy, B.; Seminara, A.; Aizenberg, J. Bacterial biofilm shows persistent resistance to

595

liquid wetting and gas penetration. P. Natl. Acad. Sci. USA. 2011, 108 (3), 995-1000.

596

(53) Malik, A.; Sakamoto, M.; Hanazaki, S.; Osawa, M.; Suzuki, T.; Tochigi, M.; Kakii, K. Coaggregation

597 598

among nonflocculating bacteria isolated from activated sludge. Appl. Environ. Microb. 2003, 69 (10),

599

(54) Malik, A.; Kakii, K. Intergeneric coaggregations among Oligotropha carboxidovorans and Acinetobacter

600

species present in activated sludge. FEMS Microbiol. Lett. 2003, 224 (1), 23-28.

601

(55) Phuong, K.; Kakii, K.; Nikata, T. Role of Acinetobacter johnsonii S35 isolate in floc-formation in activated

602

sludge process. J. Biotechnol. 2010, 150, S44-S44.

603

(56) Stampi, S.; De Luca, G.; Varoli, O.; Zanetti, F. Occurrence, removal and seasonal variation of thermophilic

604

campylobacters and Arcobacter in sewage sludge. ZBL. Hyg. Umweltmed. 1999, 202 (1), 19-27.

605 606

(57) Aminov, R. I.; Walker, A. W.; Duncan, S. H.; Harmsen, H. J. M.; Welling, G. W.; Flint, H. J. Molecular

607

gut bacteria related to Roseburia spp. or Eubacterium rectale. Appl. Environ. Microb. 2006,72 (9), 6371-6376.

608 609

(58) Walker, A. W.; Ince, J.; Duncan, S. H.; Webster, L. M.; Holtrop, G.; Ze, X. L.; Brown, D.; Stares, M. D.;

610 611

Dominant and diet-responsive groups of bacteria within the human colonic microbiota. ISME J. 2011, 5 (2),

612

(59) Culligan, E. P.; Sleator, R. D.; Marchesi, J. R.; Hill, C. Functional metagenomics reveals novel salt

613

tolerance loci from the human gut microbiome. ISME J. 2012, 6 (10), 1916-1925.

614

(60) Bagge, E.; Persson, M.; Johansson, K. E. Diversity of spore-forming bacteria in cattle manure,

615

slaughterhouse waste and samples from biogas plants. J. Appl. Microbiol. 2010, 109 (5), 1549-1565.

616

(61) Mardis, E. R. A decade's perspective on DNA sequencing technology. Nature 2011, 470 (7333), 198-203.

6056-6063.

diversity, cultivation, and improved detection by fluorescent in situ hybridization of a dominant group of human

Scott, P.; Bergerat, A.; Louis, P.; McIntosh, F.; Johnstone, A. M.; Lobley, G. E.; Parkhill, J.; Flint, H. J.

220-230.

617

30

ACS Paragon Plus Environment

Page 30 of 34

Environmental Science & Technology

Abiotrophia defectiva Achromobacter piechaudii Achromobacter xylosoxidans Acidaminococcus fermentans Acinetobacter baumannii Acinetobacter haemolyticus Acinetobacter johnsonii Acinetobacter junii Acinetobacter lwoffii Acinetobacter radioresistens Actinomyces odontolyticus Aerococcus viridans Aeromonas hydrophila Arcobacter butzleri Bacillus anthracis Bacillus cereus Bacillus pumilus Bacillus subtilis Bacillus thuringiensis Bacteroides caccae Bacteroides eggerthii Bacteroides fragilis Bacteroides ovatus Bacteroides pectinophilus Bacteroides stercoris Bacteroides thetaiotaomicron Bacteroides uniformis Bacteroides vulgatus Bifidobacterium dentium Bilophila wadsworthia Bordetella avium Bordetella parapertussis Bordetella pertussis Chromobacterium violaceum Citrobacter koseri Clostridium botulinum Clostridium difficile Clostridium perfringens Collinsella aerofaciens Comamonas testosteroni Delftia acidovorans Eggerthella lenta Eikenella corrodens Enterobacter cancerogenus Enterobacter cloacae Enterococcus casseliflavus Enterococcus faecalis Enterococcus faecium Enterococcus gallinarum Escherichia coli Eubacterium limosum Eubacterium rectale Finegoldia magna Fusobacterium mortiferum Fusobacterium nucleatum Fusobacterium ulcerans Fusobacterium varium Gardnerella vaginalis Gordonia bronchialis Haemophilus parainfluenzae Klebsiella pneumoniae Mannheimia haemolytica Megamonas hypermegale Mycobacterium abscessus Mycobacterium avium Mycobacterium bovis Mycobacterium kansasii Mycobacterium leprae Mycobacterium marinum Mycobacterium smegmatis Mycobacterium ulcerans Neisseria elongata Neisseria flavescens Neisseria gonorrhoeae Neisseria meningitidis Neisseria mucosa Neisseria sicca Neisseria subflava Nocardia farcinica Ochrobactrum anthropi Propionibacterium acnes Pseudomonas fluorescens Pseudomonas stutzeri Ralstonia pickettii Rhodococcus equi Rhodococcus erythropolis Rothia dentocariosa Sebaldella termitidis Shigella boydii Shigella dysenteriae Shigella flexneri Shigella sonnei Staphylococcus aureus Staphylococcus saprophyticus Stenotrophomonas maltophilia Streptococcus agalactiae Streptococcus anginosus Streptococcus bovis Streptococcus gordonii Streptococcus mitis Streptococcus mutans Streptococcus salivarius Streptococcus suis Sutterella wadsworthensis Tsukamurella paurometabola Veillonella atypica Veillonella dispar Veillonella parvula Vibrio cholerae Vibrio furnissii Vibrio mimicus Vibrio parahaemolyticus Vibrio vulnificus

Pathogen species number Total abundance (%)

Abiotrophia_defectiva Achromobacter_piechaudii Achromobacter_xylosoxidans Acidaminococcus_fermentans Acinetobacter_baumannii Acinetobacter_haemolyticus Acinetobacter_johnsonii Acinetobacter_junii Acinetobacter_lwoffii Acinetobacter_radioresistens Actinomyces_odontolyticus Aerococcus_viridans Aeromonas_hydrophila Arcobacter_butzleri Bacillus_anthracis Bacillus_cereus Bacillus_licheniformis Bacillus_mycoides Bacillus_pumilus Bacillus_subtilis Bacillus_thuringiensis Bacteroides_caccae Bacteroides_eggerthii Bacteroides_fragilis Bacteroides_ovatus Bacteroides_pectinophilus Bacteroides_stercoris Bacteroides_thetaiotaomicron Bacteroides_uniformis Bacteroides_vulgatus Bifidobacterium_dentium Bilophila_wadsworthia Bordetella_avium Bordetella_parapertussis Bordetella_pertussis Brevibacillus_brevis Chromobacterium_violaceum Citrobacter_koseri Clostridium_botulinum Clostridium_difficile Clostridium_perfringens Collinsella_aerofaciens Comamonas_testosteroni Corynebacterium_amycolatum Delftia_acidovorans Eggerthella_lenta Eikenella_corrodens Enterobacter_cancerogenus Enterobacter_cloacae Enterococcus_casseliflavus Enterococcus_faecalis Enterococcus_faecium Enterococcus_gallinarum Escherichia_coli Eubacterium_limosum Eubacterium_rectale Finegoldia_magna Fusobacterium_mortiferum Fusobacterium_nucleatum Fusobacterium_ulcerans Fusobacterium_varium Gardnerella_vaginalis Gordonia_bronchialis Haemophilus_parainfluenzae Kingella_denitrificans Klebsiella_pneumoniae Legionella_pneumophila Mannheimia_haemolytica Megamonas_hypermegale Mycobacterium_abscessus Mycobacterium_avium Mycobacterium_bovis Mycobacterium_kansasii Mycobacterium_leprae Mycobacterium_marinum Mycobacterium_smegmatis Mycobacterium_ulcerans Neisseria_elongata Neisseria_flavescens Neisseria_gonorrhoeae Neisseria_meningitidis Neisseria_mucosa Neisseria_sicca Neisseria_subflava Nocardia_farcinica Ochrobactrum_anthropi Propionibacterium_acnes Pseudomonas_fluorescens Pseudomonas_stutzeri Pseudoramibacter_alactolyticus Ralstonia_pickettii Rhodococcus_equi Rhodococcus_erythropolis Rothia_dentocariosa Sebaldella_termitidis Shigella_boydii Shigella_dysenteriae Shigella_flexneri SWH-ADS Shatin-Eff. Shatin-Inf. Shatin-AS Shatin-AS Stanley-AS Stanley-BFShatin-ADS Shigella_sonnei Foaming Staphylococcus_aureus Staphylococcus_epidermidis Staphylococcus_saprophyticus 91 7 4 37 27 24 6 3 Stenotrophomonas_maltophilia Streptococcus_agalactiae Streptococcus_anginosus 11.8 23.2 0.63 0.33 0.64 0.56 1.23 0.11 Streptococcus_bovis Streptococcus_gordonii Streptococcus_mitis Streptococcus_mutans Streptococcus_salivarius Streptococcus_suis Sutterella_wadsworthensis T sukamurella_paurometabola Veillonella_atypica Veillonella_dispar Veillonella_parvula Vibrio_cholerae Vibrio_furnissii Vibrio_mimicus Vibrio_parahaemolyticus Vibrio_vulnificus

5.0

1

1.0

0

0.1

-1

0.01

-2

0.001

-3

Abundance (%)

Page 31 of 34

0.0001

-4

ND

-5

Figure 1 Species and relat ive abundance of bacterial pathogen in influent, effluent, activated sludge, biofilm and anaerobic digestion sludge ( listed in CCL3;

represents emerg ing/re-emerg ing bacterial pathogen;

represents bacterial pathogen

represents bacterial pathogen transmitted through drinking-water concerned by WHO). Blue color

in the heat map indicates the relative abundance of zero.

SWH.ADS

Shatin.ADS

Stanley.BF

Stanley.AS

Shatin.AS

Shatin.AS.Foaming

Shatin.Inf.Winter

Shatin.Inf.Summer

Shatin.Eff.Winter

Shatin.Eff.Summer

ACS Paragon Plus Environment

Environmental Science & Technology

100

(a)

4

Non-emerging/re-emerging pathogens Vibrio vulnificus Vibrio parahaemolyticus Vibrio cholerae Staphylococcus aureus Shigella dysenteriae Neisseria meningitidis Neisseria gonorrhoeae Mycobacterium ulcerans Mycobacterium marinum Mycobacterium avium Klebsiella pneumoniae Escherichia coli Enterococcus faecium Enterococcus faecalis Clostridium difficile Bacillus anthracis Aeromonas hydrophila

80

2 40

Abundance (%)

Percentage (%)

3 60

1 20

0

0

Abundance

Shatin-Eff.

Shatin-Inf.

100

0.4

(b)

80

60 0.2 40

Abundance (%)

0.3 Percentage (%)

Page 32 of 34

Non-emerging/re-emerging pathogens Mycobacterium ulcerans Mycobacterium marinum Mycobacterium leprae Mycobacterium bovis Mycobacterium avium Clostridium botulinum Bordetella pertussis Aeromonas hydrophila

0.1 Abundance

SWH-ADS

Shatin-ADS

Stanley-BF

Stanley-AS

Shatin-AS

0

Shatin-AS Foaming

20

0.0

Figure 2 Emerging/re-emerging bacterial pathogens in (a) influent/effluent and (b) AS/BF/ADS. Outer ring: Percentage of emerging/re-emerg ing bacterial pathogens in total bacterial pathogen. Inner ring: Relat ive abundance of emerg ing/re-emerging bacterial pathogens in the corresponding samples . The “abundance” obtained via MetaPhlAn was estimated by weighting read counts assigned using the direct method with the total nucleotide size of all the markers in the clade and normalizing by the sum of all direct ly estimated weighted read counts. The “percentage” represents the ratio of abundance of each individual pathogen divided by the total pathogen abundance in one specific sample.

ACS Paragon Plus Environment

Page 33 of 34

Environmental Science & Technology

Figure 3 Removal of bacterial pathogens in activated sludge process . pathogens.

represents bacterial pathogens listed in CCL3.

represents emerg ing/re-emerg ing bacterial

represents bacterial pathogens transmitted through

drinking-water concerned by WHO.

ACS Paragon Plus Environment

Environmental Science & Technology

Page 34 of 34

0.5

IV

P2-Percent variation explained 22.0%

0.4

Stanley-BF Stanley-AS

III

0.3

SWH-ADS

0.2 0.1

I

II

Shatin-Eff.

0.0

Human gut-1 Shatin-Inf.

-0.1 Human gut-2 -0.2

V Shatin-AS Shatin-ADS Shatin-AS Foaming

-0.3 -0.4 -0.6

-0.5

-0.4

-0.3

-0.2

-0.1

0.0

0.1

0.2

0.3

0.4

0.5

P1- Percent variation explained 31.8% Figure 4 Principal coordinate analysis (PCoA) of bacterial pathogen profiles in influent, effluent, activated sludge, biofilm, anaerobic digestion sludge and human gut samples. PCoA was conducted based on the Bray-Curtis distance calculated from the matrix of bacterial pathogen relative abundance using PAST software.

ACS Paragon Plus Environment