Subscriber access provided by KEAN UNIV
Ecotoxicology and Human Environmental Health
Arsenic and Sulfamethoxazole Increase the Incidence of Antibiotic Resistance Genes in the Gut of Earthworm Hong-Tao Wang, Yong-Guan Zhu, Qiaoqiao Chi, Dong Zhu, Gang Li, Jing Ding, Xin-Li An, Fei Zheng, and Xi-Mei Xue Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.9b02277 • Publication Date (Web): 02 Aug 2019 Downloaded from pubs.acs.org on August 2, 2019
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 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 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.
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 37
Environmental Science & Technology
1
Arsenic and Sulfamethoxazole Increase the Incidence of
2
Antibiotic Resistance Genes in the Gut of Earthworm
3 4
Hong-Tao Wang,†,‡ Qiao-Qiao Chi,† Dong Zhu,†,‡ Gang Li,† Jing Ding,§ Xin-Li An,†
5
Fei Zheng,†,‡ Yong-Guan Zhu,†,‡,§ and Xi-Mei Xue*,†
6 7
†
8
Academy of Sciences, 1799 Jimei Road, Xiamen 361021, China
9
‡ University
Key Laboratory of Urban Environment and Health, Institute of Urban Environment, Chinese
of Chinese Academy of Sciences, 19A Yuquan Road, Beijing 100049, China
10
§
11
Sciences, Chinese Academy of Sciences, Beijing 100085, China
State Key Laboratory of Urban and Regional Ecology, Research Center for Eco-Environmental
12 13 14 15
*Address Correspondence to Xi-Mei Xue,
16
Address: Institute of Urban Environment, Chinese Academy of Sciences, 1799 Jimei
17
Road, Xiamen 361021, China
18
Phone number: +86(0)592 6190559
19
Fax number: +86(0)592 6190977
20
Email address:
[email protected] 21 22 23 24
ACS Paragon Plus Environment
Environmental Science & Technology
25
ABSTRACT
26
Combinations of metal(loid) contamination and antibiotics are considered to increase
27
the abundance of resistance genes in the environment, whereas the combined effect of
28
metal(loid)s and antibiotics on microbial communities and antibiotic resistance genes
29
(ARGs) in the gut of soil fauna remains unknown. We investigated herein the
30
alteration of ARGs and the gut microbial communities after the earthworm Metaphire
31
sieboldi was exposed to arsenate and/or sulfamethoxazole using high-throughput
32
quantitative PCR and Illumina sequencing analysis. Arsenic accumulation in the body
33
tissues of arsenic-exposed earthworms exerted a significant inhibition on growth and
34
survival. The synergistic interactions of arsenic and sulfamethoxazole increased
35
significantly the incidence of ARGs and mobile genetic elements in the earthworm
36
gut microbiota. In addition, co-exposure to arsenic and sulfamethoxazole altered the
37
structure of the gut microbial communities, and the changes correlated with ARG
38
profiles of the gut microbiota. Our results indicate that the gut of soil fauna is a
39
neglected hotspot of antibiotic resistance.
40 41 42 43 44 45 46 47
ACS Paragon Plus Environment
Page 2 of 37
Page 3 of 37
Environmental Science & Technology
48
TOC ART
49 50
The TOC was drawn by Dr. Hong-Tao Wang
51 52 53 54
ACS Paragon Plus Environment
Environmental Science & Technology
55
INTRODUCTION
56
Metal(loid)s generally coexist with antibiotics in animal manures and agricultural soil
57
amended with swine manure.1 Antibiotics, widely applied as human and veterinary
58
medicines to treat bacterial infections or to promote the growth of animals, may enter
59
the terrestrial environment when soil is fertilized with manure or sludge.2 Higher
60
concentrations of antibiotics were shown in organic vegetable fields than in traditional
61
fields due to the repeated manure application.3 Antibiotics in agricultural soil exhibit
62
various adsorption affinities on soil, and degradation rate of antibiotics can be
63
affected by their initial concentrations, soil microbial activities and oxygen status.4
64
Sulfamethoxazole belongs to sulphonamides, a broad-spectrum antibiotic against
65
most of gram-positive and many gram-negative bacteria by inhibiting their
66
multiplication, and is widely utilized to protect the health of animals or human.5
67
Sulfamethoxazole is a relatively recalcitrant antibiotic1 and can exert selective
68
pressure on soil microbial communities.3 Moreover, metalloid arsenic may be
69
introduced into managed agricultural soils due to the application of inorganic and
70
organic fertilisers, pesticides, sewage sludge, composts and livestock wastes where
71
arsenic compounds were utilized as food additives.6 Metal(loid)s can persist in soil for
72
long periods because these elements are non-degradable and difficult to be removed.7
73
Combined pollution of antibiotics and toxic metal(loid)s posed a major threat to soil
74
ecosystem and human health by altering the diversity and composition of soil
75
microbial communities7, 8 or exerting a co-selection pressure on antibiotic resistance
76
of microorganisms.1, 9
ACS Paragon Plus Environment
Page 4 of 37
Page 5 of 37
Environmental Science & Technology
77
Microbial antibiotic resistance is encoded by antibiotic resistance genes (ARGs),
78
which are regarded as emerging environmental contaminants.10 Unlike traditional
79
chemical contaminants of environment, ARGs can be spread vertically (bacterial
80
proliferation) and horizontally (exchange of genetic information between bacteria),
81
and can be transported globally.10-12 Multidrug-resistant pathogens acquiring
82
resistance from the environmental microbiome will be “superbugs”, one of the
83
greatest threats to human health.13 Antibiotic residues in the environment are
84
considered to be a selective pressure of ARGs emergence and transmission.14
85
Numerous antibiotics residues could induce selection of bacteria and promote the
86
horizontal gene transfer, increasing ARGs abundance.2,
87
ARGs abundance is sometimes independent of antibiotics and is affected by toxic
88
metals such as copper, nickel and zinc.16, 17 Arsenic contamination affects the growth
89
and development of organisms and poses a potential threat to human health via food
90
chain.18, 19 In addition, arsenic in soils is ingested and accumulated by non-target soil
91
organisms via ingestion and dermal contact, subsequently affecting the survival and
92
the microbial community structure of soil animals.20,
93
frequently detected together in agricultural fields after long-term application of
94
livestock manure due to the dietary additions of metalloids and antibiotic growth
95
promoters in the livestock industry.1 Metal(loid) contamination was considered to
96
increase the proliferation of ARGs via several co-selection mechanisms including
97
co-resistance (metal(loid) resistance genes and ARGs present on the same genetic
98
element, such as in plasmids, integrons or transposons), cross-resistance (the same
21
ACS Paragon Plus Environment
15
Moreover, alteration of
Arsenic and antibiotics are
Environmental Science & Technology
99
genetic element accounts for resistance to both of metal(loid)s and antibiotics), or
100
co-regulation (transcriptional or translational responses to metal(loid)s or antibiotics
101
form an interconnected response to either stress).9 Arsenic in swine manure has been
102
shown to increase the abundance of ARGs and mobile genetic elements (MGEs),
103
suggesting that co-selection of arsenic and antibiotic occurred in the environment.1
104
However, the combined effects of soil contamination with arsenicals and antibiotics
105
on ARGs in the gut microbiota of non-target soil fauna remains to be investigated.
106
Soil fauna, an important reservoir of biodiversity and soil quality bio-indicators,
107
is involved in the decomposition of organic matter and nutrient cycling in the soil
108
ecosystem service via their gut microbiota.22-24 The earthworm Metaphire sieboldi, a
109
worldwide distribution and highly productive species dwelling in the soil surface,
110
lives mainly on plant litter and can be applied in microcosm experiments of soil
111
fauna.25,
112
(As(V)) and sulfamethoxazole. The main aims of the study were 1) to detect the
113
impact of arsenic or/and antibiotic on the gut microbial communities by sequencing
114
bacterial 16S rRNA gene, 2) to characterize the diversity and abundance of ARGs in
115
the earthworm gut exposed to As(V) or/and sulfamethoxazole with high throughput
116
quantitative polymerase chain reaction (HT-qPCR), 3) to study the co-occurrence
117
patterns between the gut microbiota and ARGs.
118
MATERIALS AND METHODS
119
Organisms and reagents.
120
used in this study were purchased from Aohai company in Nanjing, China.
26
In this study, M. sieboldi was treated with combinations of Arsenate
Sexually mature earthworms of the species M. sieboldi
ACS Paragon Plus Environment
Page 6 of 37
Page 7 of 37
Environmental Science & Technology
121
Earthworms with the same age and approximately the same weight were cultured in
122
the collected field soil under the controlled conditions (20 oC, 12 h light/12 h dark
123
cycle and moderate relative air humidity (70%)) prior to test. Sodium arsenate
124
(Na3AsO₄·12H2O, CAS 15120-17-9, purity, >99%) was used in this study at 100 mg
125
kg−1 to monitor effect on microbiota and ARGs in the gut. Sulfamethoxazole (CAS
126
723-46-6, purity >99%) was applied at 10 mg kg-1.
127
Assay conditions. The earthworms were cultured in soil from an arsenic-free
128
vegetable field in Xiamen City, Southeast China, and oatmeal was used as the diet
129
(5-8 g oatmeal kg−1 dry weight soil) by thoroughly mixing the soil with the oatmeal
130
powder. Soil spiked with the contaminants (As(V) alone (at a final concentration of
131
100 mg kg-1 dry weight soil), sulfamethoxazole alone (at a final concentration of 10
132
mg kg−1 dry weight soil), combined contamination with As(V) (100 mg kg-1) and
133
sulfamethoxazole (10 mg kg-1), respectively) was obtained by combining soil with an
134
appropriate volume of As(V) or/and sulfamethoxazole in solution. The soil was added
135
by an equal volume of Milli-Q water (Millipore, 18.2 MΩ cm at 25 oC) as control.
136
Totally, there were four treatments including the control and three spiked soils (i.e.
137
As(V) alone, sulfamethoxazole alone, combined treatment), and each treatment has
138
three replicates. For each replicate, about 900 g of the soil was added to a
139
polyethylene plastic box (25×15×12 cm) which was covered with a vent lid to keep
140
ventilation. Ten adult earthworms similar in size were put in one box, and the
141
experiment was carried out at 70% relative humidity and with a 12:12 h light/dark
142
period in a 20 oC incubator which was supplemented regularly with Milli-Q water to
ACS Paragon Plus Environment
Environmental Science & Technology
143
keep relative humidity constant during the assays. After 28 days, the earthworms were
144
counted and weighed.
145
Sample collection and DNA extraction. All earthworms collected were killed with
146
liquid nitrogen and immediately transferred into absolute ethyl alcohol for a few
147
seconds and washed five times with sterilized water to reduce interference of
148
microorganisms from the body wall. Earthworm guts were removed under sterile
149
conditions with aseptic forceps, and transferred into a 2-mL Eppendorf tube. Total
150
DNA from 0.5 g of soil or gut was extracted by using a FastDNA® Spin Kit for Soil
151
(MP Biomedicals, CA) according to the manufacturer's directions. Agarose gel
152
electrophoresis (1.0%) and spectrophotometric analysis (Nanodrop ND-1000, Thermo
153
Fisher) were used to detect the quality and concentration of DNA. Soil and gut DNA
154
were stored at -20 oC until use.
155
16S rRNA gene amplification, sequencing and data analysis. The V4−V5 region of
156
the bacterial 16S rRNA gene was amplified with the forward primer 515F:
157
GTGCCAGCMGCCGCGG and reverse primer 907R with the unique barcode
158
CCGTCAATCMTTTRAGTTT. The PCR and gel purification of PCR fragments
159
were carried out as described previously.27 The thermal cycle of PCR consisted of an
160
initial enzyme activation at 95 °C for 5 min, 35 cycles of amplification: 95 °C for 30
161
s, 58 °C for 30 s, and 72 °C for 30 s. The amplicons reclaimed from agar gel were
162
sequenced on the Illumina HiSeq2500 platform (Novogene, Tianjin, China).
163
Quantitative Insights Into Microbial Ecology (QIIME, version 1.8.0)28 was used to
164
analyze the high-throughput sequencing data. The operational taxonomic units
ACS Paragon Plus Environment
Page 8 of 37
Page 9 of 37
Environmental Science & Technology
165
(OTUs) was defined at the ≥ 97% sequence similarity using UCLUST29 and OTUs
166
with only one representative sequence were removed prior to downstream analysis.
167
The taxonomic classification of OTUs was assigned against the reference sequences
168
in the Greengenes 13.8 16S rRNA gene database,30 and sequences was aligned via a
169
PyNAST aligner.31 The bacterial alpha-diversity was estimated with the metrics
170
observed species (OTUs) and alpha diversity indicators, such as Chao1, PD Whole
171
tree, and rarefaction curves. Principal coordinate analysis (PCoA) based on
172
Bray−Curtis distance and Adonis test was performed to compare microbial
173
communities from different samples.
174
HT-qPCR of ARGs. HT-qPCR was run on the Wafergen SmartChip Real-time PCR
175
System (Warfergen Inc., Fremont, CA) to estimate the diversity and abundance of
176
ARGs in the soil and earthworm gut. A total of 296 primer sets targeting 285 ARGs,
177
10 MGEs including 8 transposase genes and one clinical class 1 integron, one class 1
178
integron and one 16S rRNA gene (Table S1) in one ARGs run were used in this study.
179
For each run, a non-template control including three replicates in each primer pair was
180
amplified. The reaction system composed of nuclease-free PCR grade water, 1 ×
181
LightCycler 480 SYBR Green I Master mix, bovine serum albumin, primers, and
182
DNA template. The HT-qPCR program for ARGs analysis was carried out as
183
described previously.24 This program consisted of an initial 10 min denaturation step
184
at 95 °C, followed by 40 cycles of denaturation for 30 s at 95 °C and integrated
185
annealing for 30 s at 60 °C. The melting processes and qPCR results were all
186
generated by the corresponding software. Wells with amplification efficiencies
ACS Paragon Plus Environment
Environmental Science & Technology
187
(90%−110%) were accepted, and corresponding ARGs would to be discarded if only
188
one of three technical replicates were not amplified. In addition, the fold-change value
189
of ARGs (FC value) was used to indicate the enrichment of ARGs in the treated gut
190
compared to the control according to a previous study.27 A threshold cycle (CT) value
191
(31) was set as the detection limit, and relative copy number of ARGs and FC value
192
were calculated according to the formula:
193
R = Relative Gene Copy Number = 10((31−CT)/(10/3))
194
Normalized ARG Copy Number = (R /Relative16S rRNA Gene Copy Number) × 4.1
195
ΔCT = CT(ARG)−CT(16S)
196
ΔΔCT = CT(Target)−ΔCT(Ref)
197
FC = 2(−ΔΔCT)
198
CT was the threshold cycle, Target was the amended sample, Ref was the control
199
sample.27 To minimize potential variation in DNA extraction efficiencies, the relative
200
abundance of ARGs was normalized by 16S rRNA gene and converted to ARGs
201
copies per cell, of which 4.1 is the average number of 16S rRNA gene per
202
bacterium.32
203
Total and extractable arsenic. The freeze-dried earthworm body tissues and gut
204
contents were ground to fine powder in an agate mortar with liquid nitrogen prior to
205
analysis. The soil was dried at 25 oC and pulverized to 100 mesh prior to digestion.
206
The homogenized body tissues (30 mg) were precisely weighed into 50 mL
ACS Paragon Plus Environment
Page 10 of 37
Page 11 of 37
Environmental Science & Technology
207
polypropylene tubes containing 7 mL of HNO3: H2O2 mixture (5+2 v/v) and allowed
208
to stand at room temperature for 2 h. Approximately 30 mg gut contents and 200 mg
209
soil samples were accurately weighed into 50 mL polypropylene digestion tubes
210
consisting of 7 mL of HNO3: HF mixture (6+1, v/v). All the tubes were transferred to
211
the microwave-accelerated system (CEM Microwave Technology Ltd., Buckingham,
212
UK). The procedure of microwave system was carried out following three stages:
213
105 °C for 20 min, 160 °C for 20 min, 180 °C for 30 min. The extraction of
214
bioavailable arsenic was constructed as follows. Soil (100 mg) and gut contents (30
215
mg) were dissolved in 5 mL of 0.05 M aqueous ammonium sulfate,33 and mixed on a
216
rotary wheel at 450 rpm for overnight. The mixture was centrifuged (5000 g, 15 min)
217
to gather the supernatant. The digest and extraction solution of all samples were
218
diluted to 50 mL with Milli-Q water and filtered through 0.22 μm syringe filters
219
(PVDF, Millipore, USA) prior to analysis. Arsenic concentration was determined by
220
ICP-MS 7500cx (Agilent technologies, USA). Arsenic species was analysed using
221
HPLC (Agilent 1200, Agilent Technologies, USA)-ICP-MS (Agilent 7700, Agilent
222
Technologies, USA) based on the previous method.26 The method was validated
223
against certified reference material (CRM, GBW07403 and GBW10024) from the
224
National Institute of Metrology of China. The recovery rates of CRM ranged from
225
92.1% to 106.7%. The Bio-concentration Factors (BF) was calculated using the ratio
226
of the total arsenic concentration in the earthworm body tissues to that in the
227
corresponding soil. The extraction efficiency and arsenic species of all samples was
228
shown in Table S2. Other soil properties were listed in Table S3. The analysis method
ACS Paragon Plus Environment
Environmental Science & Technology
229
of soil properties was presented in Supporting Information.
230
Statistical Analysis. Mean values, standard errors (SE), and standard deviations (SD)
231
of all data were calculated using Microsoft Excel 2010. The abundance of bacterial
232
species and related genes in samples are shown as mean ± SE, and weight of
233
earthworm are presented as mean ± SD. PCoA based on the Bray−Curtis distance, the
234
diversity index (Shannon, PD Whole tree and Chao1), Adonis test, Procrustes test,
235
and mantel test were performed in R version 3.5.1. One-way analysis of variance,
236
t-tests and other significance tests were conducted with IBM SPSS V20 package.
237
RESULTS
238
Earthworm growth, mortality and arsenic bioaccumulation. The survival and the
239
growth rate of M. sieboldi under arsenic treatment, including As(V) alone and
240
combined with antibiotic, were decreased by 56.7-76.7% and 65-70% respectively,
241
compared to the control, in which lethality and growth inhibition were not found after
242
28-day culture (ANOVA, P < 0.01) (Figure 1a). Arsenic was not bio-accumulated by
243
the earthworm in the control and the treatment with sulfamethoxazole alone, of which
244
BF indexes were lower than 1. However, total arsenic concentrations in the body
245
tissues of M. sieboldi exposed to arsenic (As(V) alone and combined treatment), with
246
up to 190.1 and 145.5 mg kg-1, respectively, were much higher than those in the
247
corresponding soil, resulting in the increase of BF values (1.35-1.88) (Figure 1b).
248
Moreover, total arsenic concentrations of the body tissues were higher than those of
249
guts in the As(V) alone or the combination treatments. The extraction efficiency of
ACS Paragon Plus Environment
Page 12 of 37
Page 13 of 37
Environmental Science & Technology
250
arsenic in the gut was significantly higher than that in soil (28.6%) (t-test, P 1 (S was the standard deviation of the ΔΔCt value), was considered for
664
the enrichment. If 2−(ΔΔCt-2S) < 1, was considered for the decrease.
665 666
Figure 4. Donut charts of bacterial relative abundance at the phylum (a) and family
ACS Paragon Plus Environment
Environmental Science & Technology
667
(b) levels in the earthworm gut and surrounding soil.
668 669
Figure 5. Procrustes analysis and Mantel test of the relationship between the gut
670
microbial communities and ARG profiles based on Bray−Curtis dissimilarity metrics.
ACS Paragon Plus Environment
Page 32 of 37
Page 33 of 37
Environmental Science & Technology
ACS Paragon Plus Environment
Environmental Science & Technology
ACS Paragon Plus Environment
Page 34 of 37
Page 35 of 37
Environmental Science & Technology
ACS Paragon Plus Environment
Environmental Science & Technology
.9 %
3%
25
2% 4.
%
9.
4.3
%
2.7
3.9% %
3.8%
3.1% 3.0%
% 5 . 4 1
3
2.4%
2.7% ACS Paragon Plus Environment
. 2 %
10.5%
4 . 3 %
4 . 6 %
36%
4.6%
4.7% 4 . 7 %
14.2%
2.4%
%
9.1
4%
46.
.9%
4.1
5.1%
32
8.7%
6.7%
6% 8.
% 5 . 3
%
%
1
9.9
40
% 4 . 4
Page 36 of 37
Page 37 of 37
Environmental Science & Technology
ACS Paragon Plus Environment