Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES
Article
Antibiotic resistome and its association with bacterial communities during sewage sludge composting Jianqiang Su, Weiying Ouyang, Bei Wei, Fuyi Huang, Yi Zhao, Huijuan Xu, and Yong-Guan Zhu Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.5b01012 • Publication Date (Web): 27 May 2015 Downloaded from http://pubs.acs.org on June 1, 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 30
Environmental Science & Technology
1
TITLE
2
Antibiotic resistome and its association with bacterial communities during sewage
3
sludge composting
4
Jian-Qiang Su1, Bei Wei1,2, Wei-Ying Ou-Yang1,2, Fu-Yi Huang1, Yi Zhao1, Hui-Juan, Xu1,2,
5
Yong-Guan Zhu1*
6
1. Key Lab of Urban Environment and Health, Institute of Urban Environment,
7
Chinese Academy of Sciences, 1799 Jimei Road, Xiamen 361021, China.
8
2. University of Chinese Academy of Sciences, Beijing 100049, China.
9 10
*Correspondence: Yong-Guan Zhu, 1799 Jimei Road, Xiamen, 361021, China
11
E-mail:
[email protected] 12
Phone: +86-592-6190997
13
Fax: +86-6190977
14
ACS Paragon Plus Environment
Environmental Science & Technology
15
ABSTRACT
16
Composting is widely used for recycling of urban sewage sludge to improve soil
17
properties, which represents a potential pathway of spreading antibiotic resistant
18
bacteria and genes to soils. However, the dynamics of antibiotic resistance genes
19
(ARGs) and the underlying mechanisms during sewage sludge composting was not
20
fully explored. Here, we used high-throughput quantitative PCR and 16S rRNA gene
21
based illumina sequencing to investigate the dynamics of ARGs and bacterial
22
communities during a lab-scale in-vessel composting of sewage sludge. A total of 156
23
unique ARGs and mobile genetic elements (MGEs) were detected encoding
24
resistance to almost all major classes of antibiotics. ARGs were detected with
25
significantly increased abundance and diversity, distinct patterns, and were enriched
26
during composting. Marked shift in bacterial community structures and compositions
27
were observed during composting, with Actinobacteria being the dominant phylum
28
at the late phase of composting. The large proportion of Actinobacteria may partially
29
explain the increase of ARGs during composting. ARGs patterns were significantly
30
correlated with bacterial community structures, suggesting that the dynamic of ARGs
31
was strongly affected by bacterial phylogenetic compositions during composting.
32
These results imply that direct application of sewage sludge compost on field may
33
lead to the spread of abundant ARGs in soils.
34
Key words: urban sewage sludge, composting, antibiotic resistance genes, bacterial
35
community
ACS Paragon Plus Environment
Page 2 of 30
Page 3 of 30
Environmental Science & Technology
36
TOC art
37 38
ACS Paragon Plus Environment
Environmental Science & Technology
Page 4 of 30
39 40
INTRODUCTION
41
The increasing emergence and spread of antibiotic resistant bacteria (ARB) and
42
resistance genes (ARGs) have caused intensive concern worldwide, threatening the
43
achievements of modern medicines and posing risks to human health
44
incomplete metabolism in human bodies and disposal of unused antibiotics have
45
resulted in the release of large amounts of antibiotics into municipal wastewater 3.
46
The presence of antibiotics and biological treatment processes in wastewater
47
treatment plants (WWTPs) may potentially stimulate the development of antibiotic
48
resistance 4. While, traditional WWTPs are not designed for the removal of
49
antibiotics or ARGs, and ARB were frequently detected in WWTPs effluent, even after
50
disinfection or mixed-media filtration
51
hotspots for continually spreading antibiotics and antibiotic resistance into
52
environments 4, 7.
5, 6
1, 2
. The
. WWTPs are considered as one of the
53 54
Sewage sludge is one of the richest reservoirs of antibiotic resistance owing to the
55
input of ARB and ARGs from human and veterinary waste water 8. It is characterized
56
with high density and diversity of microbial flora, which may facilitate horizontal gene
57
transfer (HGT) via mobile genetic elements (MGEs), such as plasmids
58
integrons 11. By using culture-based method 12, quantitative PCR 13 and metagenomic
59
investigation
60
sludge, and that these antibiotic resistance factors were associated with clinical
61
settings, raising a public health issue when disseminated in downstream
62
environment of the sewage plant 15, 16.
14
9, 10
and
, diverse and abundant ARB and ARGs had been detected in sewage
63 64
The high antibiotic resistance and large amount of sewage sludge produced by
65
WWTPs pose great economic and environmental challenges for safe treatment and
66
recycling. It is estimated that more than 600 wastewater treatment plants have been
ACS Paragon Plus Environment
Page 5 of 30
Environmental Science & Technology
67
established by the year 2007 in China, producing around 5,000 million kg of sewage
68
sludge every year 17. Composting is aerobic digestion of sewage sludge to create an
69
end product that can be applied as soil conditioner and fertilizer to improve soil
70
physical prosperities, especially texture, water holding capacity and soil fertility 18-20.
71
About 120,000 tonnes of class B-equivalent biosolids are applied to soil annually in
72
Ontario, Canada
73
continue to grow because of constraints and environmental concerns of land-filling
74
and incineration
75
metals, bacterial pathogens and organic pollutants other than pharmaceutical
76
residues and ARGS. Land application of biosoilds and compost could lead to the input
77
of antibiotics into the soil
78
tetracyclines
79
soils and further to vegetables
80
resistance genes and integrons
81
integrons were much slower than previous reports 25. These findings suggested that
82
mitigating the spread of ARGs should pay more attention to the shift of ARGs during
83
sewage sludge treatment.
24
21
22
. The interest of land application of treated sewage sludge
. However, recycling of sewage sludge mainly focused on heavy
23
, including fluoroquinolones, macrolides, and
. The emergence of ARGs was also detected in biosolids-amended 21
, exhibiting enrichment of tetracycline, macrolides
11
, and that the decay rates for ARGs and class 1
84 85
Degradation of antibiotics and decline of ARGs during manure composting have been
86
reported previously
87
digestion of sewage sludge on specific types of ARGs, such as tetracyclines,
88
sulfonamides, marcrolides resistance genes 29. Although sewage sludge is recognized
89
as a rich reservoir of ARGs
90
compost
91
sludge is not well characterized, and the mechanisms underlying the shift of ARGs
92
has not been fully explored. By combining high throughput quantitative PCR (295
93
primer sets targeting almost all major classes of ARGs and MGEs) and illumina
94
sequencing of bacterial 16S rRNA gene, this study aimed: 1) to characterize the shift
21, 30
26-28
. However, only a few studies focus on the effect of various
8
and diverse ARGs are frequently detected in sludge
, the temporal change of ARGs during aerobic composting of sewage
ACS Paragon Plus Environment
Environmental Science & Technology
95
of the abundance and patterns of ARGs during a lab-scale composting of urban
96
sewage sludge; and 2) to investigate the dynamics of bacterial communities to
97
address the factors affecting the ARGs variation.
98 99
EXPERIMENTAL SECTION
100
Composting experiment set-up
101
A lab-scale of urban sewage sludge composting was carried out with dewatered
102
sewage sludge collected from the Tong’an Sewage Treatment Plant, Xiamen, China.
103
Raw materials for composting, sawdust and rice straw (Table S1) were crushed in a
104
mill (2 mm mesh size), autoclaved and were thoroughly mixed with sewage sludge
105
using a blender at the ratio of 1:2:2 (v/v). The composting experiment was conducted
106
for 50 days in quadruplicate, each containing ca. 35 kg of mixture in a PVC box (65 cm
107
in length, 50 cm in width and 40 cm in height). Moisture was maintained between 50%
108
and 65% by adding ddH2O irregularly during the composting 31, 32.
109 110
Sample collection and DNA extraction
111
Samples were collected on days 0, 1, 2, 8, 20 and 50, including the following stages: (i)
112
municipal sewage sludge before composting was collected as control sample, day 0;
113
(ii) mesophilic phase (calefactive phase) at 45 0C, day 1; (iii) thermophilic phase
114
above 55 0C, day 2 (60 0C), this phase maintained over 55 0C for 2 days; (iv) cooling
115
phase at 45 0C, day 8; (v) maturation phase at 35 0C and 30 0C, day 20 and 50,
116
respectively. Samples were collected by mixing subsamples from the upper, central
117
and lower portion of the compost uniformly to achieve high representativeness 31. All
118
the samples were stored at -80 0C for further analysis.
119 120
DNA was extracted from 0.5 g fresh samples using FastDNA Spin Kit for soil (MP
121
Biomedical, France) according to the manufacturer’s instructions. The total DNA
ACS Paragon Plus Environment
Page 6 of 30
Page 7 of 30
Environmental Science & Technology
122
were eluted with 100 μL of provided DES solution and stored at -20 0C until use. The
123
concentration and quality of extracted DNA were checked by spectrophotometric
124
analysis using NanoDrop ND-1000 (Nanodrop, USA).
125 126
High-throughput quantitative PCR (HT-qPCR) and data analysis
127
To evaluate the abundance of ARGs in samples, high-throughput quantitative PCR
128
(HT-qPCR) of ARGs were performed using the SmartChip Real-time PCR (Warfergen
129
Inc. USA) as described previously with a slight modification 7. A total of 296 primer
130
sets (Table S2) were used including 293 validated and used primer sets targeting
131
284 ARGs conferring resistance to major classes of antibiotics, 8 transposase and 1
132
16S rRNA gene
133
(intI) 35 and 1 for clinical class 1 integron-integrase gene (cintI) 36. HT-qPCR data was
134
preprocessed as described previously, for each primer set, amplifications efficiency
135
beyond the range (90% - 110%) were discarded, and amplification was confirmed
136
with more than two positive replicates
137
figure out ARGs’ fold change (FC value) of compost samples compared to the control
138
33, 38
139
absolute copy numbers by normalizing to 16S rRNA gene copy numbers which was
140
quantified separately from the Wafergen platform 37, 39. The average number of 16S
141
rRNA-encoding genes per bacterium is currently estimated at 4.1 based on the
142
Ribosomal RNA Operon Copy Number Database (rrnDB version 4.3.3) 40. Bacterial cell
143
numbers was then estimated by dividing 16S rRNA gene copy numbers by this value
144
and normalized copy number of ARGs per bacterial cell was calculated 41.
33
, 1 for blaNDM-1 34, 1 for universal class I integron-integrase gene
7, 37
. Comparative CT method was used to
. Relative copy number of ARGs and MGEs was calculated and transformed to
145 146
16S rRNA gene amplication, sequencing and data processing
147
To investigate the bacterial community structures and compositions, the V3 region of
148
bacterial 16S rRNA gene were amplified, purified, quantified, pooled and sequenced
149
on an Illumina Hiseq2000 platform at BGI, Shenzhen, China
ACS Paragon Plus Environment
42
. Raw pair-end reads
Environmental Science & Technology
Page 8 of 30
150
were assembled after filtering adaptor, low-quality reads, ambiguous N and barcode
151
to generate clean joined reads capturing the complete V3 region of the 16S rRNA
152
gene by BGI. The generated high quality sequences were processed and analyzed
153
using Quantitative Insights Into Microbial Ecology (QIIME)
154
operational taxonomic unit (OTU) picking was performed following the online
155
instruction of QIIME. OTUs were defined at the 97% similarity level using UCLUST
156
clustering 44. Representative sequence of each OTU was selected by default method,
157
which was assigned to taxonomy using RDP classifier with a confidence threshold of
158
0.80 (80%)
159
Sequences were aligned using PyNAST aligner
160
using the Fasttree algorithm 47 for downstream analysis. The levels of alpha diversity
161
were estimated using the metrics observed species (OTU), chao1 and PD Whole tree,
162
and rarefaction curves were generated to compare the level of bacterial OTU
163
diversity. The difference between microbial communities was compared using the
164
weighted and unweighted Unifrac metric and Bray-Curtis distances followed by
165
principal coordinate analysis (PCoA) and Adonis test. All sequences have been
166
deposited in the National Center for Biotechnology Information Sequence Read
167
Archive under the accession number SRP052220.
43
. The open-reference
45
. Global singletons were removed prior to downstream analysis. 46
and a phylogenetic tree was built
168 169
Statistical analysis
170
Averages, standard deviations, and fold change values of ARGs were determined
171
using Excel 2013 (Microsoft Office 2013, Microsoft, USA). ARGs were considered
172
statistically enriched or decreased if the range created by two standard deviations of
173
the mean fold change was entirely >1 or