Antibiotic Resistome and Its Association with Bacterial Communities

in Swine Manure During Simulated Composting and Lagoon Treatments Microb. ...... variable and constant aeration regimes during sewage sludge biodr...
0 downloads 0 Views 1MB Size
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