Accepted Manuscript New Insights into Antibiotic Resistome in Drinking Water and Management Perspectives: A Metagenomic based Study of Small-sized Microbes Liping Ma, Bing Li, Tong Zhang PII:
S0043-1354(19)30030-2
DOI:
https://doi.org/10.1016/j.watres.2018.12.069
Reference:
WR 14368
To appear in:
Water Research
Received Date: 7 June 2018 Revised Date:
18 December 2018
Accepted Date: 22 December 2018
Please cite this article as: Ma, L., Li, B., Zhang, T., New Insights into Antibiotic Resistome in Drinking Water and Management Perspectives: A Metagenomic based Study of Small-sized Microbes, Water Research, https://doi.org/10.1016/j.watres.2018.12.069. This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
ACCEPTED MANUSCRIPT
RI PT
Graphic Abstract
Drinking Water
ARG Risks
Antimicrobial Resistome Profiling
SC
Pathogenic Bacteria Tracking
Pathogenic Host Identification Two-step Filtration
M AN U
Pore Size = 0.45 μm Pore Size = 0.2 μm
Uncultivated Small-sized Microbes
AC C
EP
TE D
(0.2–0.45 μm)
Public Health
ACCEPTED MANUSCRIPT
New Insights into Antibiotic Resistome in Drinking Water and
2
Management Perspectives: A Metagenomic based Study of
3
Small-sized Microbes
4
RI PT
1
Authors: Liping Ma1,2, Bing Li3 and Tong Zhang1,4,5*
6
Author affiliation:
7
1
Environmental Biotechnology Laboratory, The University of Hong Kong, Hong Kong
8
2
Shanghai Key Laboratory for Urban Ecological Processes and Eco-Restoration, School
9
of Ecological and Environmental Sciences, East China Normal University, Shanghai,
M AN U
SC
5
10
China
11
3
12
Environmental Safety, Graduate School at Shenzhen, Tsinghua University, China
13
4
School of Public Health, The University of Hong Kong, Hong Kong
14
5
International Center for Antibiotic Resistance in the Environment, School of
15
Environmental Science and Engineering, Southern University of Science and Technology,
16
Shenzhen, China
17
Corresponding author:
18
Prof. Tong Zhang
19
Environmental Biotechnology Lab, The University of Hong Kong, Pokfulam Road, Hong
20
Kong
21
Email:
[email protected];
[email protected]
22
Tel: +852-28578551
23
Fax: +852-28598987
AC C
EP
TE D
Guangdong Provincial Engineering Research Center for Urban Water Recycling and
ACCEPTED MANUSCRIPT
Abstract
25
The proliferation of antibiotic resistance genes (ARGs) in drinking water and their
26
potential horizontal transfer to pathogenic microbes may cause failure of antibiotics.
27
However, antimicrobial resistome monitoring in drinking water is not currently routine.
28
The bacterial hosts of ARGs, especially small-sized microbes in drinking water, may not
29
be effectively removed by membrane filtration disinfection and thus pose threats to
30
human health. In the present study, using metagenomic based approach, we investigated
31
antibiotic resistome of small-sized microbes (0.2–0.45 µm) in 20 household drinking
32
water samples from 12 cities in Mainland China, Hong Kong and Singapore. A total of
33
265 ARG subtypes belonging to 17 ARG types were detected at abundances ranging
34
from 4.0×10-2 to 1.0×100 copies/cell. Multidrug, bacitracin and aminoglycoside resistance
35
genes are dominant, and 43 ARG subtypes were specifically carried by small-size
36
microbes. Metagenomic assembly strategy revealed fragments of three opportunistic
37
pathogen, i.e. Pseudomonas alcaligenes, Pseudomonas aeruginosa and Mycobacterium
38
gordonae, carried mexW, aph(3')-I and aac(2')-I, respectively. Drinking water samples
39
were classified into three groups based on the presence of ARG, pathogen and ARG-
40
carrying pathogen. These new insights into the antibiotic resistome of small microbes in
41
drinking water over a broad scale indicate the need for more comprehensive ARGs
42
monitoring and surveillance of drinking water supplies. These findings, together with the
43
perspectives and strategies proposed in this study, could support initiatives to improve
44
drinking water safety.
45
Keywords: drinking water, public health, antibiotic resistome, pathogenic host,
46
metagenomic sequencing
AC C
EP
TE D
M AN U
SC
RI PT
24
ACCEPTED MANUSCRIPT
1. Introduction
48
Increases in the selective pressure to develop resistance to antibiotics used in human
49
medicine and livestock breeding have led to the wide distribution and exacerbation of
50
antibiotic-resistant bacteria (ARB) and antibiotic-resistance genes (ARGs) throughout
51
both natural and man-made environments (Nesme et al. 2014, Rodriguez-Mozaz et al.
52
2015). Antimicrobial resistance has thus become a global issue and is considered an
53
emerging environmental pollutant (UNEP 2017) that may pose risks to public health
54
(WHO 2018). Nevertheless, most relevant research efforts have focused on the sources of
55
ARGs, such as hospital environments and sewage (Ashbolt et al. 2013, Le Page et al.
56
2017), or management options (e.g., nutrient management, runoff control and
57
infrastructure upgrades) with the aim of reducing the emission of ARGs into the
58
environment (Pruden et al. 2013). To date, however, inadequate strategies have been
59
proposed to monitor and manage ARGs in household drinking water, which potentially
60
pose direct threats to human health (Su et al. 2018). Remarkably, although the presence
61
of ARGs is a health threat, a more critical issue related to ARGs is the potential
62
horizontal transfer of ARGs to pathogenic microbes as the statement implies (Ashbolt
63
2015). Accordingly, these implied potential risks make it imperative to better integrate
64
monitoring efforts and risk management strategies into public health decision-making to
65
enhance drinking water safety.
66
Previous studies of ARGs in drinking water have mainly used conventional methods,
67
such as bacterial isolation (Vaz-Moreira et al. 2017), PCR (Figueira et al. 2011) and
68
qPCR (Su et al. 2018). In the past decade, advances in genomic technologies have
69
yielded novel approaches and new visions for environmental health monitoring and ARG
AC C
EP
TE D
M AN U
SC
RI PT
47
ACCEPTED MANUSCRIPT
risk assessment (Cote et al. 2016). The sequencing of whole microbial communities has
71
revealed genomic information from the vast majority of unculturable bacteria (Port et al.
72
2014). Given the potential health risks associated with drinking water, researchers have
73
recently attempted to clarify the antibiotic resistomes (collection of all the ARGs) of
74
drinking water microbes using metagenomic sequencing technique (Chao et al. 2013, Shi
75
et al. 2013). However, current studies on this resistome have only focused on microbes
76
larger than 0.45 µm (Chao et al. 2013, Ma et al. 2017), and no report has addressed the
77
resistome of small microbes (50% of the ORFs on a contig were annotated
190
as the same taxon, this contig was then assigned to that taxon (Ishii et al. 2013). The
191
average coverage, indicating abundance (×/Gb), of each ARG-carrying contig was
192
quantified by mapping metagenomic sequencing reads to the gene cassette using CLC
193
Genomics Workbench with a minimum similarity of 95% and over 95% of the read
194
length. The species-level taxonomic list of ACCs was then further compared with the
195
summarized pathogen list proposed in a previous study to identify ARG-carrying
AC C
EP
TE D
M AN U
SC
RI PT
173
ACCEPTED MANUSCRIPT
pathogen (Li et al. 2015a). Furthermore, virulence factors carried by ARG-carrying
197
pathogen was screened by using VFDB (Virulence Factors of Pathogenic Bacteria) based
198
online annotation pipeline VFanalyzer (Chen et al. 2016).
199
RI PT
196
2.5. Pathogen tracking from small-sized microbes
201
To identify pathogen from small-sized microbes in drinking water, MetaPhlAn2
202
(Metagenomic Phylogenetic Analysis) (Truong et al. 2015) was used to conduct
203
taxonomic profiling of metagenomic sequences. The species taxa and corresponding
204
abundances were subsequently extracted from the obtained taxonomic profiles. The
205
species-level taxonomic list was then further compared with a previously summarized
206
pathogen list (Li et al. 2015a), to identify small-sized pathogen in the drinking water
207
samples.
TE D
208
M AN U
SC
200
2.6. qPCR
210
Abundance of important ARGs with prevalence in all samples were further quantified by
211
using quantitative real time PCR (qPCR) assay. The unit of ARGs were transformed into
212
relative abundance by normalizing to 16S rRNA gene, as copy of ARG per 16S rRNA
213
gene. Primers of ARGs and 16S rRNA gene were listed in supplementary materials Table
214
S9. Each 20 µL qPCR reaction consisted of 10 µL 2× SYBR® Green Master Mixes (Bio-
215
Rad, Hercules, California), 0.5 µM each primer, 5 ng DNA as template and 7 µL
216
nuclease-free PCR-grade water. The thermal cycle consisted of a 3-min initial enzyme
217
activation at 95 °C, followed by 45 cycles of denaturation at 95 °C for 30 s, annealing at
AC C
EP
209
ACCEPTED MANUSCRIPT
60 °C for 30 s and extension 72 °C for 60 s. Reactions were performed on a Bio-Rad
219
CFX96 Connect™ Real-Time PCR Detection System and quantified with Bio-Rad CFX
220
manager software (Bio-Rad, Hercules, California). All qPCRs were performed in
221
technical triplicates with negative controls. Specific amplification of qPCR was
222
confirmed by checking the melting curve of the target gene in each sample.
RI PT
218
SC
223
3. Results and discussion
225
3.1. Antibiotic resistome of small-sized microbes in drinking water
226
First, a metagenomic technique based on the structured SARG database was used to
227
profile the antibiotic resistome of small-sized microbes in drinking water and present a
228
comprehensive spectrum of the ARGs carried by these microbes that had
229
opportunistically passed through membrane disinfection and thus represent a direct risk
230
to human health via drinking water. A total of 265 ARG subtypes classified into 17 ARG
231
types were detected in small-sized microbes from the 20 drinking water samples collected
232
in 12 cities (Fig. 1, Table 1), with total abundances ranging from 4.0×10-2 to 1.0×100
233
copies/cell. Until now, very limited studies had been employed to show the abundance of
234
ARGs in household drinking water/tap water, because of its low biomass and limitation
235
of sample collection. Among the limited previous studies, Xu et al., 2016 used high-
236
throughput qPCR to quantify diverse ARGs in tap water, and the relative abundance of
237
ARGs normalizing to cell number ranged from 5.75×10-1 to 8.25×10-1 copies/cell. This is
238
in the range of our detection with 4.0×10-2 to 1.0×100 copies/cell.
AC C
EP
TE D
M AN U
224
ACCEPTED MANUSCRIPT
Genes encoding multidrug- (1.6×10-2 to 5.0×10-1 copies/cell), bacitracin- (1.4×10-2 to
240
2.6×10-1 copy/cell) and beta-lactam-resistance (2.6×10-4 to 1.2×10-1 copies/cell) were
241
detected at the highest levels of abundance in these microbes. Furthermore, these top 3
242
ARG types, together with genes conferring resistance to macrolide-lincosamide-
243
streptogramin (MLS), aminoglycoside, sulfonamide and fosmidomycin, were detected as
244
general ARG types in small-sized microbes from all samples. Among the 265 ARG
245
subtypes, bacA (bacitracin) and mexF (multidrug) were dominant, with respective
246
abundances of 1.4×10-2 to 2.4×10-1 and 1.5×10-3 to 6.5×10-2 copies/cell in the drinking
247
water samples.
248
Fig. 2 presents the profiles of 56 ARG subtypes carried by small-sized microbes in
249
drinking water samples with average abundances of >1.0×10-3 copies/cell. Notably, 27
250
ARG subtypes were prevalent in the drinking water samples from all cities. These
251
encoded resistance to bacitracin (bacA and bcrA), fosmidomycin (rosA and rosB), MLS
252
(macB and macA), multidrug (acrB, acrF, adeB, amrB, bpeF, ceoB, emrB, mdtB, mdtC,
253
mdtF, mexB, mexF, mexT, mexW, multidrug ABC transporter, multidrug transporter,
254
ompR, oprC), polymyxin (arnA), sulfonamide (sul1) and vancomycin (vanR). Among
255
them, 17 ARG subtypes were defined as generalist ARG subtypes present in all samples
256
(Fig. 2).
257
Moreover, primers of 4 generalist ARGs (bacA, acrB, mexF and sul1) were available in
258
previous studies by literature search, and seven samples were randomly selected to
259
further verify their prevalence by using qPCR method. All the examined generalist ARGs
260
were detected in these samples (supplementary materials Table S10). mexF and sul1 had
261
similar abundances comparing with using metagenomic based approach, while bacA and
AC C
EP
TE D
M AN U
SC
RI PT
239
ACCEPTED MANUSCRIPT
acrB were averagely ~100 times less. The metagenomic based ARG detection mainly
263
relies on ARG database, and therefore the most complete ARG database SARG v2.0
264
(containing 23 ARG types, 1227 ARG subtypes and 12307 reference sequences) was
265
used in the present study. In SARG v2.0 database, there are totally 241 non-redundant
266
bacA sequences and 117 acrB sequences. Thus, it covers more sequence diversity, not
267
only to annotate/amplify for specific sequence or region of target genes.The generalist
268
ARG types, namely, bacitracin, multidrug, sulfonamide, beta-lactam and aminoglycoside,
269
have previously been reported to be relatively highly abundant in drinking water
270
metagenomes (Li et al. 2015b, Xu et al. 2016). Bacterial strains isolated from drinking
271
water were also frequently reported to be resistant to these antibiotics (Figueira et al.
272
2011, Xi et al. 2009). At the ARG subtype level, bacA, macB, acrB and sul1 have
273
previously been detected in drinking water samples (Ma et al. 2017), whereas some of the
274
detected ARG subtypes, such as blaLCR-1, blaLRA-17 and blaOXA-142, were
275
identified in drinking water for the first time. As this is the first large-scale study to
276
determine the dominant and prevalent antibiotic resistome among small-sized microbes in
277
drinking water samples, we aimed to further clarify whether some of the detected ARGs
278
were specifically carried by specific small-sized microbes in drinking water, and how
279
these might pose risks to human health.
SC
M AN U
TE D
EP
AC C
280
RI PT
262
281
3.2. Specific ARGs in small-sized microbes in drinking water
282
To reveal specific ARG patterns in small-sized microbes collected from drinking water,
283
we compared the ARG profile with that of large-sized microbes (>0.45 µm) in the same
284
samples. Notably, the abundance of ARGs carried by small-sized microbes (4.0×10-2 to
ACCEPTED MANUSCRIPT
1.0×100 copies/cell) was, on average, 7.4% less than that carried by large-sized microbes
286
(4.9×10-2 to 1.2×100 copies/cell) (P