Multiple Transcriptional Mechanisms Collectively Mediate Copper

Mar 28, 2019 - We investigated the mechanism of Cu resistance in Cupriavidus gilardii CR3, a Cu-resistant bacterium with a fully sequenced, annotated ...
0 downloads 0 Views 821KB Size
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 g1 to 49.81 mg g1 and

193

intracellular bioaccumulation increased from 42.46 mg g1 to 194.29 mg g1 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 dm3 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