Subscriber access provided by UNIV OF LOUISIANA
Ecotoxicology and Human Environmental Health
Multiple transcriptional mechanisms collectively mediate copper resistance in Cupriavidus gilardii CR3 Ning Huang, Juan Mao, Yan Zhao, Mingzhong Hu, and Xiaoyu Wang Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b06787 • Publication Date (Web): 28 Mar 2019 Downloaded from http://pubs.acs.org on April 1, 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 30
Environmental Science & Technology
1
Multiple transcriptional mechanisms collectively
2
mediate copper resistance in Cupriavidus gilardii CR3
3
Ning Huang,† Juan Mao,‡ Yan Zhao,§ Mingzhong Hu,║ XiaoYu Wang, *, †,‡
4
†
5
University, Changchun 130117, PR China
6
‡
7
Environmental Protection, Northeast Normal University, Changchun 130117, PR
8
China
9
§
Engineering Lab for Water Pollution Control and Resources, Northeast Normal
Key Laboratory of Wetland Ecology and Vegetation Restoration of National
School of Chemistry and Environmental Engineering, Changchun University of
10
Science and Technology, Changchun 130022, PR China
11
║
12
130012, PR China
13
Abstract Bacteria resist copper (Cu) stress by implementing several metabolic
14
mechanisms. However, these mechanisms are not fully understood. We investigated
15
the mechanism of Cu resistance in Cupriavidus gilardii CR3, a Cu-resistant bacterium
16
with a fully sequenced, annotated genome. The growth of CR3 was inhibited by
17
higher Cu concentrations (≥1.0 mM) but not by lower ones (≤0.5 mM). CR3
18
accumulated Cu intracellularly (ratios of intercellular to extracellular Cu were 11.6,
School of Chemical Engineering, Changchun University of Technology, Changchun
ACS Paragon Plus Environment
1
Environmental Science & Technology
Page 2 of 30
19
4.24, and 3.9 in 0.1, 0.5, and 1.5 mM Cu treatments, respectively). A comparative
20
transcriptome analysis of CR3 respectively revealed 310 and 413 differentially
21
expressed genes under 0.5 mM and 1.5 mM Cu stress, most of which were
22
up-regulated under Cu treatment. Gene Ontology and Kyoto Encyclopedia of Genes
23
and Genomes functional enrichment analyses uncovered several genotype-specific
24
biological processes related to Cu stress. Besides revealing known Cu
25
resistance-related genes, our global transcriptomics approach indicated that sulfur
26
metabolism, iron-sulfur cluster, and cell secretion systems are involved in mediating
27
Cu resistance in strain CR3. These results suggest that bacteria collectively use
28
multiple systems to cope with Cu stress. Our findings concerning the global
29
transcriptome response to Cu stress in CR3 provide new information for
30
understanding the intricate regulatory network of Cu homeostasis in prokaryotes.
31
32
Introduction
33
Copper (Cu) is the third most used metal in the world. Because it functions as a
34
cofactor for enzymes involved in respiration in both prokaryotic and eukaryotic cells,
35
Cu is an essential micronutrient. However, excess Cu is toxic to cells.1, 2 Cu can cause
ACS Paragon Plus Environment
2
Page 3 of 30
Environmental Science & Technology
36
toxic effects on algae and fish or amphibians at as low as 5–10 ppb.3 Excessive Cu in
37
the environment is due to industrial processes and agricultural activities, such as
38
application of fertilizers, pesticides and herbicides.4 To survive in the presence of Cu,
39
many bacteria have evolved an extensive network for Cu biosorption, transport, and
40
intracellular distribution, allowing them to tolerate various extreme environments. The
41
elucidation of bacterial Cu-resistance mechanisms is important for understanding how
42
bacteria can be exploited for efficient bioremediation of Cu-contaminated sites.5, 6
43
Previous studies have identified several resistance mechanisms in bacteria that
44
can survive in Cu-rich environments, Such as extracellular sequestration, efflux,
45
bioaccumulation, and transport mechanisms.7, 8 Under the extracellular sequestration
46
scenario, extracellular polymeric substances (EPS) are secreted by cells, and
47
functional groups in the cell wall can bind Cu ions to reduce the transportation of Cu
48
inside the cell.9-14 Intracellular Cu-resistance mechanisms are more complicated
49
processes that participate in efflux, bioaccumulation, and transport.8 From a molecular
50
prospective, various copper gene clusters including cue, cus, and cop, have been
51
identified from several different bacteria, e.g. Escherichia coli, Pseudomonas
52
syringae, and Cupriavidus metallidurans CH34.15-17 The functions of some resistance
53
genes have been elucidated. For example, cusCFBA in E. coli and copRSABCD in P.
54
syringae encode Cu efflux operons,18,
55
involved in multi-Cu oxidase activity to convert periplasmic Cu+ to less toxic Cu2+ in
56
vivo.20 This multiplicity of resistance mechanisms serves as the basis for the use of
57
bacteria in microbial bioremediation.21 Even though some mechanisms of Cu
19
and cueO-controlled genes are primarily
ACS Paragon Plus Environment
3
Environmental Science & Technology
Page 4 of 30
58
resistance have been revealed, however, the bacterial response to Cu stress is a very
59
complex set of integrated processes and the details are largely unexplored.
60
The genus Cupriavidus is known for its heavy metal resistance, especially against
61
Cu.22, 23 This genus currently comprises 11 species.23-25 In particular, C. metallidurans
62
CH34 exhibits periplasmic and cytoplasmic detoxification of Cu ions by
63
copSRABCD, copI, and copF, as determined by both microarray and quantitative
64
PCR methods.17 In our previous study, we reported the complete genome sequence of
65
Cupriavidus gilardii CR326 (here after, CR3) and characterized its Cu response genes
66
under 0.5 mM Cu stress27 and its application in Cu removal.28 Our study revealed that
67
CR3 has distinctive genomic features compared with other heavy metal-resistant
68
bacteria. In most bacteria, metal resistance genes are generally encoded in plasmid(s)
69
or co-located in plasmid(s) and chromosomal DNA. CR3 has no plasmids, however,
70
and all of its heavy metal resistance genes are encoded on its two chromosomes. This
71
distinctive genomic feature of strain CR3 suggested that it may have unique heavy
72
metal resistance mechanisms. Thus, CR3 was identified as a good candidate
73
bacterium for studies on Cu resistance mechanisms.
74
Detailed transcriptome analyses can provide insights into bacterial response to
75
heavy metal stress.29 Current knowledge of bacterial Cu molecular resistance
76
mechanisms is based on the results of studies using methods such as microarrays,
77
real-time polymerase chain reaction, and gene mutation.17,
78
methods can only be used to analyze known genes and pathways involved in bacterial
79
response. The emerging technique of transcriptome sequencing (RNA-seq) is a more
30-32
However, these
ACS Paragon Plus Environment
4
Page 5 of 30
Environmental Science & Technology
80
powerful tool to determine how an organism responds transcriptionally to a particular
81
abiotic condition relative to transcriptional levels under normal physiological
82
conditions. Thus, RNA-seq can extend the application of transcriptome studies to
83
analyses of previously unidentified genes.33 To date, most RNA-seq studies on Cu
84
stress responses have focused on eukaryotic organisms, for example Arabidopsis,
85
grapevine, and hard-shelled mussel (Mytilus coruscus).34-36 Only one reported study
86
has involved the use of RNA-seq to assess the stress response of prokaryotic cells
87
(Pseudomonas aeruginosa),37 possibly because of the lack or incompleteness of
88
reference genomes for Cu-resistant bacteria. A large number of differentially
89
expressed genes (DEGs) were recently identified when C. gilardii CR3 was exposed
90
to 0.5 mM Cu,27 but the Cu-resistance mechanisms and physiological changes of this
91
strain were not explored.
92
In this study, we used C. gilardii CR3 as a representative bacterium and
93
examined its physiological changes under Cu stress by analysis of growth curves,
94
measurement of extracellular biosorption and intracellular bioaccumulation, scanning
95
electron microscopy (SEM) and confocal laser scanning microscopy (CLSM). We
96
then used RNA-seq to elucidate the Cu-resistance mechanisms of CR3. Our aims were
97
to explore: (a) the transcriptional responses of CR3 to Cu stress, (b) the pathways
98
involved in Cu resistance, and (c) the connection between extracellular responses and
99
the intercellular efflux system.
100
Materials and Methods
101
Microorganisms, media, and growth curves
ACS Paragon Plus Environment
5
Environmental Science & Technology
Page 6 of 30
102
The CR3 strain of C. gilardii used in this study was isolated from a natural
103
asphalt deposit, as reported in our previous study.26 The strain was cultivated at 30°C
104
with shaking at 150 rpm in Luria Bertani (LB) medium (10.0 g tryptone, 5.0 g yeast
105
extract, and 10.0 g NaCl per liter, pH 6.7).
106
Growth curves for CR3 were obtained by growing cells in LB medium
107
supplemented with different amounts of CuCl2 (0.1, 0.5, 1.0, and 1.5 mM). CR3
108
grown in LB medium without Cu served as the control. Cells were incubated at 30°C
109
on a shaker at 150 rpm for 144 h. OD600 values were recorded every 12 hours using a
110
UV-visible spectrophotometer (Pgeneral T6, Beijing, China) and bacterial growth
111
curves were constructed. Bacterial biomass can be correlated with optical density by
112
reference to a standard curve.38 The linear correlation between dry cell weight (DCW)
113
(g l−1) and the OD600 of strain CR3 could be described by the following equation:
114
DCW (g l−1) = 0.45848 OD600 0.00093 (with R2 = 0.99926). The OD600 was
115
considered to represent bacterial biomass.
116
Data were reported as mean values, with error bars representing the 95%
117
confidence intervals around the average value. To identify the differences in growth
118
curves under different Cu concentrations, a one-way analysis of variance with
119
interaction was conducted using SPSS software (SPSS Inc., Chicago, IL, USA).
120
Extracellular biosorption and intracellular bioaccumulation
121
To determine the extracellular biosorption and intracellular bioaccumulation of
122
Cu, CR3 was grown in LB mediaium (100 ml) containing 0 mM, 0.1 mM, 0.5 mM, or
123
1.5 mM Cu at 30°C for 96 h. The culture was centrifuged at 12,000 g for 5 min. After
ACS Paragon Plus Environment
6
Page 7 of 30
Environmental Science & Technology
124
washing three times with deionized water, the cell pellets were resuspended in
125
5 mM sterilized EDTA at 25 °C and agitated at 150 rpm for 30 min to remove Cu ions
126
adsorbed to the cell surface.39, 40 The supernatant was centrifuged at 10,000 g for
127
5 min and then used to determine Cu concentration of extracellular biosorption by
128
atomic absorption spectrometry (Shimadzu AA-6300, Kyoto, Japan). The cell pellets
129
were dried at 85°C for 24 h, and their dry weights were determined.
130
After dry weight measurements, the biomass was subjected to acid digestion
131
using 70% trace metal-grade nitric acid (Fluka, Buchs, Switzerland) at 80°C for 30
132
min and then cooled to room temperature. The concentration of Cu ions in the
133
digested cell pellet represented the fraction of Cu bioaccumulated by the cells. The
134
amounts of Cu by extracellular biosorption (q), or intracellular bioaccumulation (q),
135
are calculated on dry cell weight basis, respectively.41, 42 The formula is expressed as
136
followed: q (mg g−1) = V·C/ m,
137
concentrations (mg l−1), and m is dry weight biomass (g) respectively.
138
RNA Extraction, library preparation, and RNA-Seq
where V is
the
sample
volume, C is
the
Cu
139
Transcriptome studies were carried out by selecting two concentrations based on
140
the CR3 growth rate at different Cu concentrations: one that affected bacterial growth
141
(1.5 mM) and one that did not (0.5 mM). Free Cu served as the control. Three
142
biological repeats were performed for each treatment; and nine samples in total were
143
sequenced to RNA-Seq. CR3 cells were cultured in 100 ml LB medium to an OD600
144
value of 0.2 and then exposed to 0, 0.5, and 1.5 mM Cu. The cells were allowed to
145
continue growing to an OD600 value of 1.0, and 1 ml of bacterial culture was then
ACS Paragon Plus Environment
7
Environmental Science & Technology
Page 8 of 30
146
collected and centrifuged. RNA extraction, library preparation, and RNA-Seq were
147
performed as described in Supporting Information. Gene expression data were
148
obtained and quantified as reads per kilobase of coding sequence per million reads.43
149
Raw reads were submitted to the NCBI Sequence Read Archive under the accession
150
number SRP154974.
151
Transcriptome annotation and functional enrichment analyses
152
Gene functional annotation was presented in Supporting Information.
153
We used GOATOOLS to perform GO enrichment analysis of the DEGs44 and GO
terms
for
gene
functional
categorization.
KOBAS
154
obtained
155
(http://kobas.cbi.pku.edu.cn/home.do) was used for pathway enrichment analysis
156
based on the KEGG database, where ‘pathway’ was set as the basic unit. Terms with a
157
P-adj (corrected p-value) < 0.05 in the GO and KEGG enrichment analyses were
158
considered to be enriched.
159
Quantitative
160
(qRT-PCR)
reverse
transcriptase-polymerase
chain
reaction
analyses
161
Eight DEGs related to Cu resistance (cysD, cusA, copS, copR, copA, copB, iscA,
162
and gshA) were randomly selected for validation of the RNA-seq results by qRT-PCR.
163
A detailed description of qRT-PCR, and selection of reference genes can be found in
164
the Supporting Information (Table S1).45, 46
165
Biofilm EPS composition analysis by SEM and CLSM
166
According to the transcriptomic data, type III secretion system genes were
167
up-regulated in C. gilardii CR3. To verify the occurrence of type III secretion on the
ACS Paragon Plus Environment
8
Page 9 of 30
Environmental Science & Technology
168
surface of C. gilardii CR3 cells under Cu stress, SEM and CLSM were used to
169
examine C. gilardii CR3 in the absence (control) and presence of different Cu
170
concentrations (0.1, 0.5, and 1.5 mM) for 96 h. SEM was used to observe cellular
171
surface changes under Cu stress, while CLSM was used to investigate changes in
172
proteins and polysaccharides in biofilm EPS under different Cu conditions.47 SEM
173
and CLSM were performed as detailed in Supporting Information.
174
All experiments were performed in triplicate.
175
RESULTS
176
Effect of initial Cu concentration on growth
177
The growth curves of CR3 at lower Cu concentrations (0.1 and 0.5 mM) were
178
similar to the curve of CR3 in the control (0 mM Cu), with no significant differences
179
observed in OD600 values (p > 0.05) (Figure 1). When the Cu concentration was
180
increased to 1.0 or 1.5 mM, the bacterial biomass sharply decreased and growth was
181
significantly inhibited (p < 0.05). The maximum biomass obtained after 96 h of
182
growth in the 1.0 and 1.5 mM Cu treatments was decreased by 41.2% and 45.2%,
183
respectively, compared with that of the control. According to these results, CR3
184
growth was therefore inhibited by Cu at high concentrations (1.0 and 1.5 mM) but not
185
at low Cu concentrations (0.1 and 0.5 mM) (Figure 1).
ACS Paragon Plus Environment
9
Environmental Science & Technology
186 187 188 189
Page 10 of 30
Figure 1. Growth of Cupriavidus gilardii CR3 in culture medium supplemented with different concentrations of Cu.
Extracellular biosorption and intracellular bioaccumulation of Cu
190
To understand the extra- and intracellular bioaccumulation of Cu, the
191
bioaccumulation of Cu in C. gilardii CR3 at different Cu concentrations was
192
determined. Extracellular biosorption increased from 3.66 mg g1 to 49.81 mg g1 and
193
intracellular bioaccumulation increased from 42.46 mg g1 to 194.29 mg g1 as the Cu
194
concentration increased from 0.1 mM to 1.5 mM (Figure 2). The analysis revealed
195
that intracellular Cu was more abundant, as the ratio of intercellular to extracellular
196
Cu was 11.6, 4.24, and 3.9 in the 0.1, 0.5, and 1.5 mM Cu treatments, respectively.
197
These results demonstrate the significantly higher capacity of CR3 for intracellular Cu
198
bioaccumulation than for extracellular biosorption (p < 0.05) (Figure 2).
ACS Paragon Plus Environment
10
Page 11 of 30
199 200 201 202 203
Environmental Science & Technology
Figure 2. Amount of copper (Cu) removed due to biosorption and bioaccumulation after culturing Cupriavidus gilardii CR3 in media with different Cu concentrations. Initial Cu concentrations were 0, 0.1, 0.5, and 1.0 mM. NA: not available.
Differential gene expression under Cu exposure as determined by RNA-Seq
204
To investigate genes associated with Cu stress responses in CR3, RNA-Seq was
205
used for sequencing analysis. Details of the Illumina sequencing data for the nine
206
tested samples are listed in Table S2. In total, 220 million raw reads (2.2 Gbps) were
207
generated. After removing residual rRNA and low-quality reads, 216 million clean
208
reads (2.16 Gbps) were obtained for further alignment and analysis. Out of all the
209
RNA sequence reads, 94.2%–99.35% were mapped to the genome of CR3 (Table S2),
210
indicating that the sequencing data were completely reliable. The principal component
211
analysis showed good parallel correlations and obvious differences for each set
212
(Figure S1A).
213
Compared with CR3 in the control group, CR3 treated with 0.5 mM or 1.5 mM
214
Cu showed strong responses at the transcriptome level (Figure S1B). In total, 310
215
genes (302 up- and 8 down-regulated) and 413 genes (380 up- and 33 down-regulated)
ACS Paragon Plus Environment
11
Environmental Science & Technology
Page 12 of 30
216
were significantly differentially expressed (FDR < 0.05, |log2FC| 1) by 0.5 mM and
217
1.5 mM Cu, respectively (Figure S1B). The overall expression profile of the DEGs
218
was analyzed by hierarchical clustering (Figure S1C). The resulting heatmaps were
219
largely indistinguishable between libraries from the same treatment, but large
220
differences were observed in transcript expressions between Cu stress treatments and
221
the control (Figure S1C). The differential expression between 0.5 mM and 1.5 mM Cu
222
stress was not significant.
223
To validate the RNA-seq results, eight genes differentially expressed in response
224
to Cu treatment were randomly selected for qRT-PCR analysis. In seven of these eight
225
genes, the expression fold changes in strain CR3 were in good agreement between
226
RNA-seq and qRT-PCR (Table S3). The exception was the iscA gene, where the
227
log2FC was 2.92 according to RNA-seq vs. 0.295 by qRT-PCR (Table S3). A
228
significant correlation was observed between the RNA-Seq and qRT-PCR data (R2 =
229
0.858, p < 0.001), thus indicating that the RNA-Seq data were reliable (Figure S2).
230
GO and KEGG functional enrichment analyses of DEGs
231
To determine the possible metabolic pathways affected by Cu stress in CR3, GO
232
and KEGG pathway enrichment analyses were performed on the set of identified
233
DEGs (P-adj < 0.05). The enrichment analysis of GO terms classified the genes into
234
two functional categories: biological process and molecular function (Figure 3). In
235
total, 14 and 13 GO terms were significantly overrepresented at 0 vs. 0.5 mM and 0 vs.
236
1.5 mM Cu, respectively (Figure 3). Thirteen enriched GO terms were similar
237
between the two levels of Cu stress. A little exception was oxidoreductase activity
ACS Paragon Plus Environment
12
Page 13 of 30
Environmental Science & Technology
238
(GO:0016667), which was enriched by 0.5 mM Cu but not 1.5 mM Cu. The 13 GO
239
terms included four Cu ion binding or transporter homeostasis terms (GO:0005507,
240
0006878, 0055070, and 0006825); five sulfur compound transport or transporter
241
activity terms (GO: 1902358, 0015116, 0008272, 001541, and 1901682), two
242
metallo-sulfur or iron-sulfur (Fe-S) cluster assembly terms (GO: 0031163 and
243
0016226), one inorganic anion transmembrane transport term (GO:0098661), and one
244
metal ion binding term (GO: 0046914, 0005507) (Figure 3).
245 246 247 248 249 250 251 252
Figure 3. GO functional enrichment scatter diagram of Cupriavidus gilardii CR3 after exposure to different concentrations of copper (Cu). GO functional enrichment was represented by an enrichment factor, the P-value, and the number of unigenes enriched (enrichment factor is ratio of numbers of differentially expressed genes annotated in a particular pathway term to numbers of all genes annotated in this pathway term). Names of GO functional enrichment classes are listed along the y-axis. P-adj ranged from 0 to 0.05, with a P-adj closer to 0 indicateing greater enrichment.
253
According to the KEGG enrichment results, four pathways were significantly
254
enriched after 0.5 mM Cu treatment (Table 1): sulfur metabolism, bacterial secretion
255
system, seleno-compound metabolism, and the two-component system (P-adj < 0.05).
256
The sulfur metabolism pathway was significantly enriched after 1.5 mM Cu treatment
257
(P-adj < 0.001), and the bacterial secretion system and two-component system were
ACS Paragon Plus Environment
13
Environmental Science & Technology
Page 14 of 30
258
also enriched, with significant uncorrected p-values (p < 0.05) (Table 1). Thus, Cu at
259
different concentrations induced similar transcriptional responses in CR3. Pathways
260
associated with Cu resistance, sulfur metabolism, Fe-S cluster assembly, and cell
261
secretion system were significantly enriched under Cu stress, which suggests that they
262
play important roles in Cu resistance.
263 264 265
Table 1. Enrichment of KEGG pathways in Cupriavidus gilardii CR3 exposed to copper (Cu) (0.5 mM and 1.5 mM) compared with a non-Cu control. P-adj < 0.05 indicates significant enrichment. pathway
KEGG Id
0.5 mM
1.5 mM
Adjusted p-value
Adjusted p-value
Sulfur metabolism
ko00920
6.56E-05
2.21E-07
Bacterial secretion system
ko03070
0.006905
0.188
Two-component system
ko02020
0.032918
0.146
Selenocompound metabolism
ko00450
0.021674
0.682
266
Cu resistance genes exhibited the biggest response to both concentrations of Cu.
267
Most genes (31 for 0.5 mM Cu and 35 for 1.5 mM Cu) showed significant
268
up-regulation under Cu stress (FDR < 0.05). These DEGs included five putative Cu
269
resistance
270
copKHFIDCBARS, cusFABgene1401, and merRcopAcopZ (Table S4). Cu resistance
271
genes that were induced under Cu stress in CR3 encoded proteins that are involved in
272
the two-component signal system (copRS), copper efflux system (cusFABgene1401),
273
copper binding (copK, copI, copZ, and copH), copper ion transport (copB and copD),
274
and multicopper oxidase (copA).
275
gene
clusters:
copSRABCD,
copQLFGJIDCBARS_ompC_copKBA,
GO and KEGG analyses indicated that sulfur metabolism was enriched under Cu
ACS Paragon Plus Environment
14
Page 15 of 30
Environmental Science & Technology
276
stress (24 and 20 genes were up-regulated by 0.5 mM and 1.5 mM Cu, respectively;
277
Table 1 and Figure 3). At both concentrations, Cu induced significant up-regulation of
278
genes encoding proteins in sulfur assimilation pathways, including cysPUWA, cysND,
279
cysH, cysI, ssuP, and tauD. Additionally, genes encoding proteins involved in
280
glutathione29 biosynthesis (gshA, gshB, and E1.11.1.9), tyrosine metabolism (fahA and
281
hmgA) and the biosynthesis of phenylalanine (phhA and paaA), glutamate (putA), and
282
glycine (betA) were significantly up-regulated in CR3 treated with Cu at 0.5 mM and
283
1.5 mM (Table S4), supporting their involvement in the Cu detoxification response.
284
The GO analyses also revealed that some of the DEGs in CR3 treated with Cu at
285
0.5 mM and 1.5 mM were enriched in pathways related to Fe-S clusters (GO:0016226)
286
(Figure 3). Further analyses of our RNA-seq data revealed that Cu supplementation
287
resulted in significant up-regulation of the isc operon (iscRSU-iscA-hscBA), which
288
drives the synthesis of Fe-S clusters (9 and 10 genes up-regulated by treatment with
289
0.5 mM and 1.5 mM Cu, respectively) (Table S4).
290
Two cell secretion systems, type II secretion system (T2SS) and type III
291
secretion system (T3SS), were detected in strain CR3. Interestingly, the T3SS in CR3
292
was significantly induced under 0.5mM Cu stress (ko03070) (p < 0.01) (Table 1),
293
whereas the T2SS was not. These results suggest that Cu may affect EPS secretion
294
and the cell membrane via the T3SS.
295
Characterization of extracellular Cu in CR3
296
The SEM analyses revealed the mucous nature of the EPS surrounding bacterial
297
cells and biofilm. The EPS on the cell membrane surface was significantly thicker in
ACS Paragon Plus Environment
15
Environmental Science & Technology
Page 16 of 30
298
the 0.1mM Cu treatment than in the control (Figure 4a, b); in the 0.5 mM Cu
299
treatment, however, the EPS was very similar to that in the control (Figure 4c). In the
300
1.5 mM Cu treatment, cell surface secretion was significantly reduced and the biofilm
301
was affected (Figure 4d), which suggests that Cu affects cellular activity. The EDS
302
analysis revealed two additional Cu signals in the EDX spectrum following Cu stress
303
(Figure 4e, f), indicating that the cell surface bioaccumulated Cu.
304 305 306 307 308
Figure 4. Scanning electron micrographs and energy dispersive X-ray spectrographs of Cupriavidus gilardii CR3 exposed to different concentrations of copper (Cu). (a) 0 mM Cu; (b) 0.1 mM Cu; (c) 0.5 mM Cu; (d) 1.5 mM Cu; (e) 0 mM Cu (control); (f) 0.1 mM Cu (treated).
309
The ratio of protein and polysaccharide in biofilm EPS was determined by
310
CLSM after 96 h exposure to different Cu concentrations (Figure S3). In the biofilm
311
without any Cu exposure, the ratio of protein to polysaccharide was 0.15 ± 0.04, and
312
the composition of the biofilm EPS was mainly polysaccharides (Figure S3A). In
313
biofilm under exposure to various Cu concentrations, ratios of protein to
314
polysaccharide were 1.52 ± 0.49, 0.38 ± 0.14 and 0.93 ± 0.38 for 0.1, 0.5, and 1.5 mM
ACS Paragon Plus Environment
16
Page 17 of 30
Environmental Science & Technology
315
Cu treatments, respectively (Figure S3B, C, D), which were significantly higher than
316
the ratio in the absence of Cu treatment (t-test, p < 0.05). These results demonstrate
317
that Cu stimulated the more secretion of extracellular protein.
318
DISCUSSION
319
The goal of this work was to obtain a global description of bacterial Cu
320
resistance response mechanisms under Cu stress. The RNA-seq analyses showed that
321
Cu at both tested concentrations led to transcriptional changes in CR3 (302 up- and 8
322
down-regulated and 380 up- and 33 down-regulated genes in 0.5 mM and 1.5 mM Cu
323
Cu2+ treatments, respectively). In a previous study, 293 up- and 126 down-regulated
324
genes were detected in response to 0.5 mM Cu in P. aeruginosa.37 Although there
325
were a similar number of DEGs under Cu stress in the two bacteria, the transcriptional
326
results for P. aeruginosa did not exhibit significant enrichment of particular metabolic
327
pathways. However, analyses of CR3 revealed significant enrichment of several GO
328
and KEGG pathways under Cu stress, including Cu resistance genes, sulfur
329
metabolism, Fe-S cluster assembly, and bacterial secretion systems. The results
330
suggest that these metabolic pathways respond to Cu and affect how the bacterium
331
tolerates Cu stress.
332
Cu resistance-related genes
333
Cu resistance genes are important for bacterial persistence in Cu-contaminated
334
areas. In our study, many virulence genes were up-regulated in CR3 under Cu stress.
335
In particular, the expression levels of genes encoding copABCD, copGDCBARSKB,
ACS Paragon Plus Environment
17
Environmental Science & Technology
Page 18 of 30
336
and copKDCBARS were increased 2.95- to 12.98-fold in this study (Table S4), much
337
higher than the up-regulation of copRS, copA1, copZ1, and copZ2 reported in P.
338
aeruginosa under 0.5 mM CuSO4 stress (1.79- to 3.95-fold).37 This difference
339
indicates that Cu resistance genes have a stronger transcriptional response in CR3
340
than in P. aeruginosa. Two-component regulatory systems (copRS) can detect elevated
341
levels of Cu ions outside the cytoplasm and are involved in the regulation of Cu
342
homoeostasis genes.48 Such a system was first described in P. syringae49 and has now
343
been identified and characterized in a variety of bacteria. The Cu efflux system is
344
responsible for acquiring Cu(I) ions from the cytoplasm and facilitating their efflux
345
into the periplasmic space.50 A Cu-transporting ATPase transports a significant
346
proportion of internal Cu across cell membranes. Cu-binding proteins can bind Cu(I)
347
or Cu(II) at different sites,51-54 and multicopper oxidase can oxidize Cu+ to the less
348
toxic Cu2+.16 These pathways were enriched by Cu in CR3, suggesting that enhanced
349
Cu transport led to an increase in Cu efflux out of the cell and maintenance of Cu
350
homeostasis to resist Cu toxicity. This multi-pronged approach would allow the
351
bacterium to survive in areas contaminated with high concentrations of Cu and helps
352
to explain why CR3 can grow normally in the presence of 0.5 mM Cu2+. In contrast,
353
non-resistant bacteria show decreased microbial growth in response to Cu at lower
354
concentrations; for example, 16 μM in E. coli, 9 ppm in Desulfovibrio sp., and 100
355
mg dm3 in Rhizopus arrhizus.55-57
356
Sulfur metabolism
357
As a vital element for all living organisms, sulfur is required for the synthesis of
ACS Paragon Plus Environment
18
Page 19 of 30
Environmental Science & Technology
358
proteins and essential cofactors.58 In our study, the enrichment of genes involved in
359
sulfur assimilation suggested the important role of this process in Cu detoxification.
360
Sulfate assimilation pathway genes that were highly expressed in the presence of Cu
361
encode a key enzyme in cysteine biosynthesis (PAPS),59,
362
(CysHIJ),59 and enzymes involved in GSH precursor synthesis (Figure 5, Table S4).
363
Sulfur and cysteine have been previously reported to participate in Cu detoxification
364
in microorganisms.32, 61 For instance, Camila et al. conducted microarray analyses and
365
found that cysD, cysI, and cysN were up-regulated under 325 mM Cu stress in
366
Acidithiobacillus ferrooxidans strain D2.61 cysDN encode ATP sulfurylase which
367
converts sulfate to adenosine 5-phosphosulfate (APS),62 while cysI encodes an
368
iron-sulfur protein that uses a siroheme cofactor to reduce sulfite to sulfide.63 Garrett
369
et al. detected significant up-regulation of the M. sedula cysNHI operon in the
370
presence of 4 or 8 mM Cu, but this operon was not up-regulated by other metal
371
stimuli. 32
60
sulfite reductase
372
Sulfate assimilation pathways also provide sulfate for the synthesis of GSH,
373
glycine, and tyrosine.64 GSH is composed of a tripeptide of L-cysteine, L-glutamate
374
and glycine,65 which is able to chelate metal ions and thus increase organism
375
resistance to heavy metals.66, 67 Previous studies have extensively revealed the roles of
376
GSH in enhancing resistance to heavy metals in yeast and plant cells.64, 68-70 However,
377
there are fewer studies investigating the roles of GSH in binding Cu in bacteria. Miras
378
et al. have posited that GSH increases intracellular accumulation of Cu and proposed
379
that Cu binds to GSH and various thiol-rich proteins in E. coli.71 In this study, we
ACS Paragon Plus Environment
19
Environmental Science & Technology
Page 20 of 30
380
found that genes encoding proteins involved in GSH biosynthesis (gshA, gshB, and
381
E1.11.1.9) were significantly up-regulated in CR3 under Cu stress (Figure 5, Table
382
S4), which provides transcriptional evidence that GSH is involved in the Cu resistance
383
of C. gilardii CR3. Chelation with Cu during the process of GSH detoxification might
384
be responsible for the higher intracellular bioaccumulation of Cu in CR3. The
385
up-regulation of genes, including those related to tyrosine metabolism (fahA and
386
hmgA), methionine (msrA), and the biosynthesis of phenylalanine (phhA and paaA),
387
glutamate (putA), and glycine (betA), also contribute to the bioaccumulation of Cu in
388
CR3. These genes are all involved in the biosynthesis of Cu-resistance genes (copCD,
389
azurin) associated with binding to Cu2+ or Cu+ in cells.48 According to our study
390
results, transcription under Cu stress can be linked to the higher intercellular
391
bioaccumulation of Cu.
392
Iron-sulfur clusters
393
In organisms, Iron-sulfur clusters (Fe-S) play an important role in the growth,
394
respiration, and metabolism, such as DNA repair, gene regulation, biosynthetic
395
pathways, and aerobic and anaerobic respiration.72,
396
generally bind Fe, but may bind Cu if Cu replaces Fe in the Fe-S cluster by occupying
397
the cysteine-thiolate donor positions.55 This binding promotes the accumulation of Cu
398
in the cell. Genes encoding Fe-S clusters and cysteine biosynthesis proteins are
399
up-regulated in Bacillus subtilis under Cu stress; the higher Fe-S cluster and cysteine
400
contents would help to complex intracellular Cu, hence decreasing its toxic effect.73, 74
401
In present study, GO and KEGG analyses revealed that genes in the Fe-S cluster
73
Fe-S assembly systems
ACS Paragon Plus Environment
20
Page 21 of 30
Environmental Science & Technology
402
assembly system were significantly enriched under Cu stress (Figure 3, Table 1), thus
403
indicating that our results are consistent with those of the previous study. The
404
up-regulation of Fe-S cluster genes therefore partly contributed to the higher observed
405
intracellular accumulation of Cu in CR3.
406
Cell secretion system
407
To our surprise, we observed significant enrichment of genes encoding
408
components of the T3SS in CR3 cells grown in the presence of Cu. The T3SS can
409
secrete substrate proteins into the extracellular environment and translocate effector
410
proteins to the inside of eukaryotic cells via a syringe-like needle apparatus that
411
extends from the bacterial surface.75 Previous studies on Xanthomonas citri have
412
shown that the T3SS is essential for bacterial biofilm formation.76 The lack of the
413
T3SS led to altered expression of proteins involved in metabolism, energy, EPS, and
414
cell motility.76 In our SEM analyses, we observed a similar morphological distribution
415
of EPS between the 0.5 mM Cu treatment and the control (Figure 4), which suggests
416
that the cellular secretion system might stimulate a compensatory mechanism to
417
counteract the toxicity induced by 0.5 mM Cu. At the same time, the CLMS analysis
418
showed that protein secretion in CR3 biofilm EPS increased significantly after Cu
419
treatment (Figure S3). Some studies have demonstrated that Cu directly affects the
420
cellular activity and extracellular secretions of bacteria,77,
421
proteins may be directly involved in heavy metal resistance by enhancing
422
biosorption79-85.
423
Gram-negative bacteria.
However,
surface
display
systems
78
and that cell surface
vary
among
different
ACS Paragon Plus Environment
21
Environmental Science & Technology
Page 22 of 30
424
A role for the secretory system to counter heavy metal stress has been rarely
425
reported. In a recent study of the pathogenic bacterium Edwardsiella tarda, Cu stress
426
inhibited the expression of genes encoding components of the T3SS, resulting in
427
reduced biofilm growth and host cell adhesion.86 Our results are inconsistent with
428
those for E. tarda, genes encoding components of the T3SS were up-regulated in CR3
429
under Cu stress. This opposite results might be due to the fact that strain CR3 is a
430
typical Cu-resistant bacterium, suggesting that the T3SS probably contributes to Cu
431
resistance. An elevated Cu concentration in the external environment may stimulate
432
the T3SS in CR3, resulting in increased secretion of EPS and the adsorption of Cu
433
ions to the cell surface to prevent entry of Cu into the cell. To our knowledge, T3SS
434
involvement has not been reported in any Cu-resistant bacterium. Hence, more studies
435
are needed to verify whether the T3SS in Cu-resistant bacteria plays a role in Cu
436
resistance processes in the future.
437
Multifaceted response of CR3 to Cu exposure
438
Overall, our transcriptome analysis linked the observed physiological effects
439
under Cu stress and revealed that the resistance of CR3 to Cu is a multi-system
440
collaborative process (Figure 5). We found that CR3 exhibits four significant
441
mechanistic responses to Cu stress: Cu resistance gene expression, as well as
442
activation of the Fe-S assembly system, sulfur metabolism, and the cell secretion
443
system. Together, these responses allow CR3 to achieve a relatively stable cellular
444
environment for survival. Under excess Cu stress, genes in the cell secretion system
445
are up-regulated, and this stimulates the production of proteins on the cell surface.
ACS Paragon Plus Environment
22
Page 23 of 30
Environmental Science & Technology
446
Although extracellular proteins can enhance Cu adsorption and prevent Cu from
447
entering the cell, excessive Cu can still be transported into cells and trigger
448
up-regulated expression of Cu resistance genes. The up-regulation of Cu resistance
449
genes leads to an increase in the efflux of intercellular Cu into the periplasmic space.
450
The sulfur metabolism system is also highly expressed. The sulfur metabolism
451
pathway is closely related to the biosynthesis of cysteine and GSH, which are related
452
to the chelation and detoxification of Cu ions.32, 61, 71 We therefore infer that sulfur
453
metabolism plays an vital role in the Cu resistance process. Strong expression of the
454
Fe-S assembly system under Cu stress possibly means that Cu ions can bind to Fe
455
binding sites to achieve Cu chelation. Collectively, the results of our combined
456
transcriptome and physiological analyses provide a better understanding of the Cu
457
resistance and response mechanisms of bacteria.
458
ACS Paragon Plus Environment
23
Environmental Science & Technology
459 460 461 462 463
Page 24 of 30
Figure 5. Proposed model for resistance response mechanisms of Cupriavidus gilardii CR3 to copper (Cu) stress and C. gilardii survival strategies. Colored text indicates genes that are highly up-regulated under Cu stress: Cu-resistance genes (mustard green), sulfate assimilation pathway genes (purple), iron-sulfur cluster genes (green), Type III secretion system genes (orange), and genes encoding GSH and related enzymes (red).
464
AUTHOR INFORMATION
465
Corresponding Author
466
* Phone: +86 13504403055; E-mail:
[email protected];
467
Notes
468
The authors declare no competing financial interest.
469
ACKNOWLEDGMENTS
470
This work was funded by the National Natural Science Foundation of China (no.
471
51678122) and the Jilin Provincial Science and Technology Department Project (no.
472
20150414046GH). We also thank Dr. Yanhong Xiao, Experimental Center of the
473
School of the Environment, Northeast Normal University, for her assistance with
474
primer design.
475
Supporting Information Available
476
RNA Extraction, library preparation, and RNA-Seq; Transcriptome annotation;
477
Quantitative reverse transcriptase-polymerase chain reaction analyses (qRT-PCR);
478
Scanning electron microscope-energy dispersive spectroscopy analyses; Biofilm
479
composition determination by confocal laser scanning microscopy (CLMS);
480
Transcriptome analysis of Cupriavidus gilardii CR3 in response to different
481
concentrations of copper (Cu) (0, 0.5, and 1.5 mM) (Figure S1); Correlation between
ACS Paragon Plus Environment
24
Page 25 of 30
Environmental Science & Technology
482
RNA-Seq and qRT-PCR for the eight selected genes (Figure S2); CLSM images of
483
Cupriavidus gilardii CR3 exposed to different concentrations of copper (Cu) (Figure
484
S3) (PDF)
485
The list of primers used for qRT-PCR (Table S1), sequencing and assembly statistics
486
for the nine transcriptome data of copper treatment and control group(Table S2), the
487
data of RNA-Seq and qRT-PCR for the selected eight genes (Table S3) and list of
488
DEGs belonging to copper resistance mechanism (Table S4) were provided in Excel,
489
and separately (XLS)
490
This information is available free of charge via the Internet at http://pubs.acs.org
491
REFERENCES
492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514
(1) Amachawadi, R. G.; Scott, H. M.; Alvarado, C. A.; Mainini, T. R.; Vinasco, J.; Drouillard, J. S.; Nagaraja, T. G. Occurrence of the transferable copper resistance gene tcrB among fecal enterococci of US feedlot cattle fed copper-supplemented diets. Appl. Environ. Microbiol. 2013, 79 (14), 4369-4375. (2) Gerstein, A. C.; Ono, J.; Lo, D. S.; Campbell, M. L.; Kuzmin, A.; Otto, S. P. Too much of a good thing: the unique and repeated paths toward copper adaptation. Genetics 2015, 199 (2), 555-571. (3) Chua, H.; Yu, P. H.; Sin, S. N.; Cheung, M. W. Sub-lethal effects of heavy metals on activated sludge microorganisms. Chemosphere 1999, 39 (15), 2681-2692. (4) Flemming, C.; Trevors, J. Copper toxicity and chemistry in the environment: a review. Water, Air, & Soil Pollution. 1989, 44 (1-2), 143-158. (5) Merroun, M. L. Interactions between metals and bacteria: Fundamental and applied research. Commun. Curr Res & Educ Top & Trends. Appl Microbiol. 2007, 4, 108-119. (6) Méndez-Vilas, A., Communicating current research and educational topics and trends in applied microbiology. Formatex: Badajoz, 2007; Vol. 2. (7) Lemire, J. A.; Harrison, J. J.; Turner, R. J. Antimicrobial activity of metals: mechanisms, molecular targets and applications. Nat. Rev. Microbiol. 2013, 11 (6), 371-384. (8) Gillan, D. C. Metal resistance systems in cultivated bacteria: are they found in complex communities? Curr. Opin. Biotechnol. 2016, 38, 123-130. (9) Sklodowska, A.; Matlakowska, R. Relative surface charge, hydrohobicity of bacterial cells and their affinity to substrate during copper bioleaching from post-flotation wastes. Biotechnol. Lett. 1998, 20 (3), 229-233. (10) Sheng, G. P.; Yu, H. Q.; Li, X. Y. Extracellular polymeric substances (EPS) of microbial
ACS Paragon Plus Environment
25
Environmental Science & Technology
515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558
Page 26 of 30
aggregates in biological wastewater treatment systems: A review. Biotechnol. Adv. 2010, 28 (6), 882-894. (11) González, A. G.; Shirokova, L. S.; Pokrovsky, O. S.; Emnova, E. E.; Martínez, R. E.; Santana-Casiano, J. M.; González-Dávila, M.; Pokrovski, G. S. Adsorption of copper on Pseudomonas aureofaciens: Protective role of surface exopolysaccharides. J. Colloid Interface Sci. 2010, 350 (1), 305-314. (12) Fang, L.; Wei, X.; Cai, P.; Huang, Q.; Chen, H.; Liang, W.; Rong, X. Role of extracellular polymeric substances in Cu(II) adsorption on Bacillus subtilis and Pseudomonas putida. Bioresour. Technol. 2011, 102 (2), 1137-1141. (13) Bai, J.; Chao, Y.; Chen, Y.; Wang, S.; Qiu, R. Immobilization of Cu by Bacillus subtilis DBM and the role of extracellular polymeric substances. Water, Air, Soil Pollut. 2017, 228 (3), 86. (14) Pan, J. h.; Liu, R. x.; Tang, H. x. Surface reaction of Bacillus cereus biomass and its biosorption for lead and copper ions. J. Environ. Sci. 2007, 19 (4), 403-408. (15) Outten, F. W.; Huffman, D. L.; Hale, J. A.; O'Halloran, T. V. The independent cue and cus systems confer copper tolerance during aerobic and anaerobic growth in Escherichia coli. J. Biol. Chem. 2001, 276 (33), 30670-30677. (16) Puig, S.; Rees, E. M.; Thiele, D. J. The ABCDs of periplasmic copper trafficking. Structure. 2002, 10 (10), 1292-1295. (17) Monchy, S.; Benotmane, M. A.; Wattiez, R.; Van, S. A.; Auquier, V.; Borremans, B.; Mergeay, M.; Taghavi, S.; Van, D. D. L.; Vallaeys, T. Transcriptomic and proteomic analyses of the pMOL30-encoded copper resistance in Cupriavidus metallidurans strain CH34. Microbiology. 2006, 152 (6), 1765-1776. (18) Cha, J. S.; Cooksey, D. A. Copper resistance in Pseudomonas syringae mediated by periplasmic and outer membrane proteins. Proc. Natl. Acad. Sci. U. S. A. 1991, 88 (20), 8915-8919. (19) Franke, S.; Grass, G.; Rensing, C.; Nies, D. H. Molecular analysis of the copper-transporting efflux system CusCFBA of Escherichia coli. J. Bacteriol. 2003, 185 (13), 3804-3812. (20) Singh, S. K.; Grass, G.; Rensing, C.; Montfort, W. R. Cuprous oxidase activity of CueO from Escherichia coli. J. Bacteriol. 2004, 186 (22), 7815-7817. (21) Haferburg, G.; Kothe, E. Microbes and metals: interactions in the environment. J. Basic Microbiol. 2007, 47 (6), 453-467. (22) Slaveykova, V. I.; Pinheiro, J. P.; Floriani, M.; Garcia, M. Interactions of core-shell quantum dots with metal resistant bacterium Cupriavidus metallidurans: consequences for Cu and Pb removal. J. Hazard. Mater. 2013, 261, 123-129. (23) Chen, W. M.; Wu, C. H.; James, E. K.; Chang, J. S. Metal biosorption capability of Cupriavidus taiwanensis and its effects on heavy metal removal by nodulated Mimosa pudica. J. Hazard. Mater. 2008, 151 (2–3), 364-371. (24) Ferreira, P. A. A.; Bomfeti, C. A.; Soares, C. R. F. S.; Soares, B. L.; Moreira, F. M. S. Cupriavidus necator strains: Zinc and cadmium tolerance and bioaccumulation. Sci. Agric. 2018, 75 (6), 452-460. (25) Li, L. G.; Cai, L.; Zhang, T. Genome of Cupriavidus sp. HMR-1, a heavy metal-resistant bacterium. Genome Announce. 2013, 1 (1), e00202-12. (26) Janssen, P. J.; Wang, X.; Chen, M.; Xiao, J.; Hao, L.; Crowley, D. E.; Zhang, Z.; Yu, J.; Huang, N.; Huo, M.; Wu, J. Genome sequence analysis of the naphthenic acid degrading and
ACS Paragon Plus Environment
26
Page 27 of 30
559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602
Environmental Science & Technology
metal resistant bacterium Cupriavidus gilardii CR3. Plos One 2015, 10 (8), e0132881. (27) Huang, N.; Mao, J.; Hu, M.; Wang, X.; Huo, M. Responses to copper stress in the metal-resistant bacterium Cupriavidus gilardii CR3: a whole-transcriptome analysis. J. Basic Microbiol. 2019; DOI: 10.1002/jobm.201800693. (28) Yang, Y.; Hu, M.; Zhou, D.; Fan, W.; Wang, X.; Huo, M. Bioremoval of Cu2+ from CMP wastewater by a novel copper-resistant bacterium Cupriavidus gilardii CR3: characteristics and mechanisms. RSC Adv. 2017, 7 (30), 18793-18802. (29) Zhang, Y.; Jiang, R.; Wu, H.; Liu, P.; Xie, J.; He, Y.; Pang, H. Next-generation sequencing-based transcriptome analysis of Cryptolaemus montrouzieri under insecticide stress reveals resistance-relevant genes in ladybirds. Genomics 2012, 100 (1), 35-41. (30) Yamamoto, K.; Ishihama, A. Transcriptional response of Escherichia coli to external copper. Mol. Microbiol. 2005, 56 (1), 215-227. (31) Zhang, X.-X.; Rainey, P. B. Regulation of copper homeostasis in Pseudomonas fluorescens SBW25. Environ. Microbiol. 2008, 10 (12), 3284-3294. (32) Wheaton, G. H.; Mukherjee, A.; Kelly, R. M.; Kostka, J. E. Transcriptomes of the extremely thermoacidophilic archaeon Metallosphaera sedula exposed to metal “Shock” reveal generic and specific metal responses. Appl. Environ. Microbiol. 2016, 82 (15), 4613-4627. (33) Chen, M.; Zhang, M.; Borlak, J.; Tong, W. A decade of toxicogenomic research and its contribution to toxicological science. Toxicol. Sci. 2012, 130 (2), 217-228. (34) Bernal, M.; Casero, D.; Singh, V.; Wilson, G. T.; Grande, A.; Yang, H.; Dodani, S. C.; Pellegrini, M.; Huijser, P.; Connolly, E. L.; Merchant, S. S.; Kramer, U. Transcriptome sequencing identifies SPL7-regulated copper acquisition genes FRO4/FRO5 and the copper dependence of iron homeostasis in Arabidopsis. Plant Cell 2012, 24 (2), 738-761. (35) Leng, X.; Jia, H.; Sun, X.; Shangguan, L.; Mu, Q.; Wang, B.; Fang, J. Comparative transcriptome analysis of grapevine in response to copper stress. Sci Rep. 2015, 5, 17749. (36) Xu, M.; Jiang, L.; Shen, K. N.; Wu, C.; He, G.; Hsiao, C. D. Transcriptome response to copper heavy metal stress in hard-shelled mussel (Mytilus coruscus). Genomics Data 2016, 7, 152-154. (37) Quintana, J.; Novoa-Aponte, L.; Argüello, J. M. Copper homeostasis networks in the bacterium Pseudomonas aeruginosa. J. Biol. Chem. 2017, 292 (38), 15691-15704. (38) Myers, J. A.; Curtis, B. S.; Curtis, W. R. Improving accuracy of cell and chromophore concentration measurements using optical density. BMC Biophys. 2013, 6 (4), 1-15. (39) Sedlakova-Kadukova, J.; Kopcakova, A.; Gresakova, L.; Godany, A.; Pristas, P. Bioaccumulation and biosorption of zinc by a novel Streptomyces K11 strain isolated from highly alkaline aluminium brown mud disposal site. Ecotoxicol. Environ. Saf. 2019, 167, 204-211. (40) Hassler, C. S.; Slaveykova, V. I.; Wilkinson, K. J. Discriminating between intra- and extracellular metals using chemical extractions. Limnol. Oceanogr.: Methods 2004, 2 (4), 237-247. (41) Olguín, E. J.; Sánchezgalván, G. Heavy metal removal in phytofiltration and phycoremediation: the need to differentiate between bioadsorption and bioaccumulation. New Biotechnol. 2012, 30 (1), 3-8. (42) Velásquez, L.; Dussan, J. Biosorption and bioaccumulation of heavy metals on dead and living biomass of Bacillus sphaericus. J. Hazard. Mater. 2009, 167 (1), 713-716. (43) Cock, P. J. A.; Fields, C. J.; Goto, N.; Heuer, M. L.; Rice, P. M. The Sanger FASTQ file format for sequences with quality scores, and the Solexa/Illumina FASTQ variants. Nucleic Acids
ACS Paragon Plus Environment
27
Environmental Science & Technology
603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646
Page 28 of 30
Res. 2009, 38 (6), 1767-1771. (44) Young, M. D.; Wakefield, M. J.; Smyth, G. K.; Oshlack, A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010, 11 (2), R14. (45) Kenneth, J.; Livak, T. D. Analysis of relative gene expression data using real—time quantitative PCR and the 2(-Delta Delta C(T)) method. Method 2001, 25 (4), 402-408. (46) Vandesompele, J.; Preter, K. D.; Pattyn, F.; Poppe, B.; Roy, N. V.; Paepe, A. D.; Speleman, F. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3 (7), research0034.1. (47) Shen, Y.; Huang, C.; Monroy, G. L.; Janjaroen, D.; Derlon, N.; Lin, J.; Espinosa-Marzal, R.; Morgenroth, E.; Boppart, S. A.; Ashbolt, N. J. Response of simulated drinking water biofilm mechanical and structural properties to long-term disinfectant exposure. Environ. Sci. Technol. 2016, 50 (4), 1779-1787. (48) Rensing, C.; Mcdevitt, S. F., The copper metallome in prokaryotic cells. In Metallomics and the Cell; Banci, L., Ed. Springer Netherlands: Dordrecht, 2013; pp 417-450. (49) Mills, S. D.; Jasalavich, C. A.; Cooksey, D. A. A two-component regulatory system required for copper-inducible expression of the copper resistance operon of Pseudomonas syringae. J. Bacteriol. 1993, 175 (6), 1656-1664. (50) Su, C. C.; Long, F.; Yu, E. W. The Cus efflux system removes toxic ions via a methionine shuttle. Protein Sci. 2011, 20 (1), 6-18. (51) Bondarczuk, K.; Piotrowska-Seget, Z. Molecular basis of active copper resistance mechanisms in Gram-negative bacteria. Cell Biol. Toxicol. 2013, 29 (6), 397-405. (52) Kim, E. H.; Nies, D. H.; McEvoy, M. M.; Rensing, C. Switch or funnel: how RND-type transport systems control periplasmic metal homeostasis. J. Bacteriol. 2011, 193 (10), 2381-2387. (53) Ash, M. R.; Chong, L. X.; Maher, M. J.; Hinds, M. G.; Xiao, Z.; Wedd, A. G. Molecular Basis of the Cooperative Binding of Cu(I) and Cu(II) to the CopK Protein from Cupriavidus metallidurans CH34. Biochemistry 2011, 50 (43), 9237-9247. (54) Sendra, V.; Gambarelli, S.; Bersch, B.; Coves, J. Site-directed mutagenesis reveals a conservation of the copper-binding site and the crucial role of His24 in CopH from Cupriavidus metallidurans CH34. J. Inorg. Biochem. 2009, 103 (12), 1721-1728. (55) Macomber, L.; Imlay, J. A. The iron-sulfur clusters of dehydratases are primary intracellular targets of copper toxicity. Proc. Natl. Acad. Sci. U. S. A. 2009, 106 (20), 8344-8349. (56) Cabrera, G.; Perez, R.; Gomez, J. M.; Abalos, A.; Cantero, D. Toxic effects of dissolved heavy metals on Desulfovibrio vulgaris and Desulfovibrio sp. strains. J. Hazard. Mater. 2006, 135 (1-3), 40-46. (57) Uslu, G.; Dursun, A. Y.; Ekiz, H. İ.; Aksu, Z. The effect of Cd (II), Pb (II) and Cu (II) ions on the growth and bioaccumulation properties of Rhizopus arrhizus. Process Biochem. 2003, 39 (1), 105-110. (58) Albanesi, D.; Mansilla MCSchujman, G. E.; De, M. D. Bacillus subtilis cysteine synthetase is a global regulator of the expression of genes involved in sulfur assimilation. J. Bacteriol. 2005, 187 (22), 7631-7638. (59) Wolschendorf, F.; Ackart, D.; Shrestha, T. B.; Hascall-Dove, L.; Nolan, S.; Lamichhane, G.; Wang, Y.; Bossmann, S. H.; Basaraba, R. J.; Niederweis, M. Copper resistance is essential for virulence of Mycobacterium tuberculosis. Proc. Natl. Acad. Sci. 2011, 108 (4), 1621-1626. (60) Tetaz, T. J.; Luke, R. K. Plasmid-controlled resistance to copper in Escherichia coli. J.
ACS Paragon Plus Environment
28
Page 29 of 30
647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690
Environmental Science & Technology
Bacteriol. 1983, 154 (3), 1263-1268. (61) Salazar, C. N.; Acosta, M.; Galleguillos, P. A.; Shmaryahu, A.; Quatrini, R.; Holmes, D. S.; Demergasso, C. S. Analysis of gene expression in response to copper stress in Acidithiobacillus ferrooxidans strain D2, isolated from a copper bioleaching operation. Adv. Mater. Res. 2013, 825, 157-161. (62) Kredich, N. M. The molecular basis for positive regulation of cys promoters in Salmonella typhimurium and Escherichia coli. Mol. Microbiol. 2010, 6 (19), 2747-2753. (63) Guédon, E.; Martin-Verstraete, I., Cysteine Metabolism and Its Regulation in Bacteria. In Amino Acid Biosynthesis ~ Pathways, Regulation and Metabolic Engineering, Wendisch, V. F., Ed. Springer Berlin Heidelberg: Berlin, Heidelberg, 2007; pp 195-218. (64) Mendoza-Cózatl, D.; Loza-Tavera, H.; Hernández-Navarro, A.; Moreno-Sánchez, R. Sulfur assimilation and glutathione metabolism under cadmium stress in yeast, protists and plants. FEMS Microbiol. Rev. 2005, 29 (4), 653-671. (65) Ezraty, B.; Gennaris, A.; Barras, F.; Collet, J.-F. Oxidative stress, protein damage and repair in bacteria. Nat. Rev. Microbiol. 2017, 15 (7), 385-396. (66) Zhou, Q.; Guo, J. J.; He, C. T.; Shen, C.; Huang, Y. Y.; Chen, J. X.; Guo, J.; Yuan, J. G.; Yang, Z. Y. Comparative transcriptome analysis between low- and high-cadmium-accumulating genotypes of pakchoi (Brassica chinensis L.) in response to cadmium stress. Environ. Sci. Technol. 2016, 50 (12), 6485-6494. (67) Rubino, F. M. Toxicity of Glutathione-Binding Metals: A Review of Targets and Mechanisms. Toxics 2015, 3 (1), 20-62. (68) Wang, K.; Tang, S. F.; Hou, X. Molecular mechanism investigation on the interactions of copper (II) ions with glutathione peroxidase 6 from Arabidopsis thaliana. Spectrochim. Acta, Part A. 2018, 203, 428-433. (69) Liang, Z. Y.; Pilonsmits, E. A.; Jouanin, L.; Terry, N. Overexpression of glutathione synthetase in indian mustard enhances cadmium accumulation and tolerance. Plant Physiol. 1999, 119 (1), 73-79. (70) Zhu, Y. L.; Pilonsmits, E. A. H.; Tarun, A. S.; Weber, S. U.; Jouanin, L.; Terry, N. Cadmium tolerance and accumulation in indian mustard is enhanced by overexpressing γ-Glutamylcysteine synthetase. Plant Physiol. 1999, 121 (4), 1169-1178. (71) Miras, R.; Morin, I.; Jacquin, O.; Cuillel, M.; Guillain, F.; Mintz, E. Interplay between glutathione, Atx1 and copper. 1. Copper(I) glutathionate induced dimerization of Atx1. JBIC, J. Biol. Inorg. 2008, 13 (2), 195-205. (72) Romsang, A.; Duang-Nkern, J.; Leesukon, P.; Saninjuk, K.; Vattanaviboon, P.; Mongkolsuk, S. The iron-sulphur cluster biosynthesis regulator IscR contributes to iron homeostasis and resistance to oxidants in Pseudomonas aeruginosa. Plos One 2014, 9 (1), e86763. (73) Chillappagari, S.; Seubert, A.; Trip, H.; Kuipers, O. P.; Marahiel, M. A.; Miethke, M. Copper stress affects iron homeostasis by destabilizing iron-sulfur cluster formation in Bacillus subtilis. J Bacteriol. 2010, 192 (10), 2512-2524. (74) Fang, Z.; Dos Santos, P. C. Protective role of bacillithiol in superoxide stress and Fe-S metabolism in Bacillus subtilis. MicrobiologyOpen 2015, 4 (4), 616-631. (75) Galán, J. E.; Wolfwatz, H. Protein delivery into eukaryotic cells by type III secretion machines. Nature 2006, 444 (7119), 567-573. (76) Zimaro, T.; Thomas, L.; Marondedze, C.; Sgro, G. G.; Garofalo, C. G.; Ficarra, F. A.; Gehring,
ACS Paragon Plus Environment
29
Environmental Science & Technology
691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718
Page 30 of 30
C.; Ottado, J.; Gottig, N. The type III protein secretion system contributes to Xanthomonas citri subsp. citri biofilm formation. BMC Microbiol. 2014, 14 (1), 96. (77) Deschatre, M.; Ghillebaert, F.; Guezennec, J.; Colin, C. S. Sorption of copper(II) and silver(I) by four bacterial exopolysaccharides. Appl. Biochem. Biotechnol. 2013, 171 (6), 1313-1327. (78) Shameer, S. Biosorption of lead, copper and cadmium using the extracellular polysaccharides (EPS) of Bacillus sp., from solar salterns. Biotech 2016, 6 (2), 194. (79) Valls, M.; Atrian, S.; Lorenzo, V. D.; Fernández, L. A. Engineering a mouse metallothionein on the cell surface of Ralstonia eutropha CH34 for immobilization of heavy metals in soil. Nat. Biotechnol. 2000, 18 (6), 661-665. (80) Ueki, T.; Sakamoto, Y.; Yamaguchi, N.; Michibata, H. Bioaccumulation of Copper Ions by Escherichia coli Expressing Vanabin Genes from the Vanadium-Rich Ascidian Ascidia sydneiensis samea. Appl. Environ. Microbiol. 2003, 69 (11), 6442-6446. (81) Bae, W.; Mulchandani, A.; Chen, W. Cell surface display of synthetic phytochelatins using ice nucleation protein for enhanced heavy metal bioaccumulation. J. Inorg. Biochem. 2002, 88 (2), 223-227. (82) Sousa, C.; Kotrba, P.; Ruml, T.; Cebolla, A.; De, L. V. Metalloadsorption by Escherichia coli cells displaying yeast and mammalian metallothioneins anchored to the outer membrane protein LamB. J. Bacteriol. 1998, 180 (9), 2280-2284. (83) Mejàre, M.; Ljung, S.; Bülow, L. Selection of cadmium specific hexapeptides and their expression as OmpA fusion proteins in Escherichia coli. Protein Eng. 1998, 11 (6), 489-494. (84) Ravikumar, S.; Yoo, I. K.; Lee, S. Y.; Hong, S. H. Construction of copper removing bacteria through the integration of two-component system and cell surface display. Appl. Biochem. Biotechnol. 2011, 165 (7-8), 1674-1681. (85) Patel, J.; Zhang, Q.; Mckay, R. M.; Vincent, R.; Xu, Z. Genetic Engineering of Caulobacter crescentus for Removal of Cadmium from Water. Appl. Biochem. Biotechnol. 2010, 160 (1), 232-243. (86) Hu, Y. H.; Dang, W.; Liu, C. S.; Sun, L. Analysis of the effect of copper on the virulence of a pathogenic Edwardsiella tarda strain. Lett. Appl. Microbiol. 2010, 50 (1), 97-103.
719
ACS Paragon Plus Environment
30