Subscriber access provided by NEW YORK MED COLL
Article
Identification of putative genes involved in bisphenol A degradation using differential protein abundance analysis of Sphingobium sp. BiD32 Nicolette Angela Zhou, Henrik Kjeldal, Heidi Lois Gough, and Jeppe Nielsen Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.5b02987 • Publication Date (Web): 21 Sep 2015 Downloaded from http://pubs.acs.org on September 26, 2015
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 31
Environmental Science & Technology
1
Identification of putative genes involved in
2
bisphenol A degradation using differential protein
3
abundance analysis of Sphingobium sp. BiD32
4 5
Nicolette A. Zhou1,2, Henrik Kjeldal1, Heidi L. Gough2, Jeppe L. Nielsen1*
6 7
1
8
Aalborg, Denmark
9
2
Aalborg University, Department of Chemistry and Bioscience, Fredrik Bajers Vej 7H, DK-9220
University of Washington, Department of Civil and Environmental Engineering; More Hall 201
10
Box 352700, Seattle, Washington 98195-2700, USA.
11
*corresponding author, (+45) 99 40 85 06; fax (+45) 98 14 18 08;
[email protected] 12 13
Keywords: label-free quantitative proteomics, metabolomics, microbial degradation, wastewater
14
treatment, bisphenol A, PPCP, BPA, TOrC
ACS Paragon Plus Environment
Environmental Science & Technology
15
Abstract
16
Discharge of the endocrine disrupting compound bisphenol A (BPA) with wastewater treatment
17
plant (WWTP) effluents into surface waters result in deleterious effects on aquatic life.
18
Sphingobium sp. BiD32 was previously isolated from activated sludge based on its ability to
19
degrade BPA. This study investigated BPA metabolism by Sphingobium sp. BiD32 using label-
20
free quantitative proteomics. The genome of Sphingobium sp. BiD32 was sequenced to provide a
21
species-specific platform for optimal protein identification. The bacterial proteomes of
22
Sphingobium sp. BiD32 in the presence and absence of BPA were identified and quantified. A
23
total of 2155 proteins were identified, 1174 of these proteins were quantified, and 184 of these
24
proteins had a statistically significant change in abundance in response to the presence/absence
25
of BPA (p ≤ 0.05). Proteins encoded for by genes previously identified to be responsible for
26
protocatechuate degradation were upregulated in the presence of BPA. The analysis of the
27
metabolites from BPA degradation by Sphingobium sp. BiD32 detected a hydroxylated
28
metabolite. A novel p-hydroxybenzoate hydroxylase enzyme detected by proteomics was
29
implicated in the metabolic pathway associated with the detected metabolite. This enzyme is
30
hypothesized to be involved in BPA degradation by Sphingobium sp. BiD32, and may serve as a
31
future genetic marker for BPA degradation.
2 ACS Paragon Plus Environment
Page 2 of 31
Page 3 of 31
Environmental Science & Technology
32
1.
Introduction
33
Wastewater treatment plants (WWTP) discharge endocrine disrupting compounds (EDC) to
34
aquatic environments, resulting in deleterious effects on aquatic life. A common EDC detected in
35
WWTP effluents and surface waters is bisphenol A (BPA). BPA is a widely used plasticizer for
36
polycarbonate plastics and epoxy resins manufacturing 1 with 3.2 million metric tons produced
37
globally in 2003 2. BPA is an estrogen and androgen receptor agonist 3 and has negative effects
38
on aquatic life at concentrations measured in the environment 4–6. While BPA is partially
39
removed during wastewater treatment, residual concentrations enter the receiving aquatic
40
environments 3,7 with BPA detected in 41.2% of US streams tested and a median and maximum
41
concentration of 140 ng/L and 12,000 ng/L, respectively 8. Up to 29 ng/L BPA was detected in
42
groundwater impacted by WWTP effluent 9. New technologies are needed to improve the
43
removal of EDCs, such as BPA, during wastewater treatment.
44 45
Bioaugmentation with EDC-degrading bacteria is a possible solution to improve the removal of
46
EDCs during wastewater treatment 10. Previous studies have isolated many bacteria capable of
47
degrading various contaminants 11–15. Among these is Sphingobium sp. BiD32, which was
48
isolated from activated sludge and is capable of degrading BPA to 5 ng/L 16. As average effluent
49
concentrations range from 4.2 to 1600 ng/L 17,18, bioaugmentation with these EDC-degrading
50
bacteria can likely improve the removal of EDCs in WWTPs. Understanding how the bacteria
51
degrade these contaminants will improve bioaugmentation techniques by potentially allowing for
52
easier monitoring of EDC degradation.
53
3 ACS Paragon Plus Environment
Environmental Science & Technology
54
Identifying genes that could be used as genetic indicators of BPA degradation would enable
55
tracking of biodegradation potential in WWTPs. Proteomic studies have been effective at
56
identifying degradation genes and pathways of other EDCs and pharmaceuticals, as
57
demonstrated by recent work with ibuprofen 19 and 17β-estradiol (E2) 20. Previous proteomic
58
studies have also led to the identification of gene markers for E2 exposure using label-free
59
proteomics 21. However, this approach has not yet been applied for BPA degradation.
60 61
The objective of this study was to identify enzymes involved in BPA degradation by
62
Sphingobium sp. BiD32 and to elucidate the degradation pathway using a combination of
63
genomics, proteomics, and metabolite analysis. Label-free quantitative proteomics was used to
64
investigate the effect of BPA on the proteome of Sphingobium sp. BiD32 and the results were
65
used to identify candidate BPA catabolic genes in the organism’s genome. Analysis of the
66
metabolites from BPA then assisted in confirming potential genes involved and the pathway.
67 68
2.
Material and methods
69
2.1
Chemicals, reagents, and bacterial strain
70
BPA (≥99%) was obtained from Sigma Aldrich (St. Louis, MO, USA). Organic solvents and
71
water used for BPA sample extraction and analysis were high performance liquid
72
chromatography (HPLC) grade. A 35 mg/L BPA aqueous stock was prepared by adding BPA
73
dissolved in acetonitrile (1000 mg/L), evaporating the acetonitrile at 60ºC, adding water, and
74
autoclave sterilizing (121ºC, 20 min). A control of acetonitrile was prepared the same way as the
75
BPA aqueous stock. Sphingobium sp. BiD32 (GenBank accession number JX87937) was
76
previously isolated from activated sludge 16.
4 ACS Paragon Plus Environment
Page 4 of 31
Page 5 of 31
Environmental Science & Technology
77 78
2.2
Library preparation, genome sequencing, assembly, and analysis
79
The whole genome sequence of bacterial strain Sphingobium sp. BiD32 was prepared and
80
completed (details in Supporting Information (SI)). The sequences were trimmed and assembled
81
using de novo assembly by the CLC Genomics Workbench version 5.5.1 (CLC bio, Aarhus,
82
Denmark). During trimming, the quality score limit was set at 0.05 and reads < 50 bp were
83
discarded. Default parameters for assembly were used with the exception of the word size (62),
84
the minimum contig length (1000), and the similarity fraction (0.95). Assembled contigs were
85
curated by CodonCode Aligner v 3.7 (CodonCode Corp.) and analyzed with the Rapid
86
Annotation using Subsystem Technology (RAST) annotation server for subsystem classification
87
and functional annotation 24,23. This whole-genome shotgun project was deposited at
88
DDBJ/EMBL/GenBank under the accession numbers CAVK010000001 to CAVK010000262.
89 90
2.3
Sample preparation for proteomic analysis
91
Samples were prepared for proteomic analysis by growing Sphingobium sp. BiD32 in the
92
presence and absence of BPA (Figure 1). Sphingobium sp. BiD32 was inoculated by transferring
93
a single colony from solid media into R2B media (pre-inoculum culture, PIC), a carbon-rich
94
liquid media prepared based on the composition of the solid media R2A 23. After reaching an
95
optical density of 0.603 at 600 nm (OD600), 5 mL of Sphingobium sp. BiD32 was transferred into
96
a second flask of R2B (inoculum culture, IC). When the IC reached an OD600 of 0.302, 8 mL was
97
inoculated into experimental flasks containing R2B for an estimated initial volatile suspended
98
solids (VSS) concentration of 5 mg/L. Five replicate cultures of Sphingobium sp. BiD32 were
99
grown (100 rpm, 28ºC) with an initial concentration of 1.5 mg/L BPA and five replicates were
5 ACS Paragon Plus Environment
Environmental Science & Technology
101
carbon in the media. OD600 was measured and samples for BPA quantification were collected
102
every one to two hours. An HPLC instrument with ultraviolet detector (HPLC-UV) was used to
103
quantify total BPA concentrations (sorbed and dissolved fraction) using previously reported
104
methods 16. BPA was re-spiked into the cultures at 7.6 hours to ensure that BPA was present and
105
actively being degraded when the cells were harvested. Sphingobium sp. BiD32 cultures (with
106
and without BPA) were harvested when reaching an OD600 between 0.385 and 0.412. At the
107
point of harvest, the cultures were still in exponential growth and approximately 60% of the BPA
108
had been degraded. Cells were harvested by washing the cell pellet twice with 20 mL phosphate-
109
buffered saline (PBS) and collected by centrifugation (10,000 xg, 4ºC, 15 min). The proteins
110
were extracted and digested in-solution (details in SI). 2.5
0.5
2
0.4
1.5
0.3
1
0.2
0.5
0.1
0
0 0
111
Growth (OD600)
grown without BPA. The BPA added to five of the replicates contributed less than 0.6% of the
Bisphenol A (mg/L)
100
5
10
15
Time (hours)
112
Figure 1. Growth and harvesting of Sphingobium sp. BiD32 for differential proteomics analysis.
113
Degradation of BPA is shown as squares. Growth of Sphingobium sp. BiD32 both in the
114
presence and absence of BPA is shown as diamonds; no differences were seen between the two
115
conditions (p=0.45). BPA was re-spiked at 7.6 hours. Cells were harvested for protein extraction
116
at 11.5 hours. Error bars represent the range of measurements from replicates (n=5).
117 6 ACS Paragon Plus Environment
Page 6 of 31
Page 7 of 31
Environmental Science & Technology
118
2.4
LC-MS/MS analyses for proteomics
119
Tryptic digests from the cytosolic and membrane fractions were analyzed separately by an
120
automated liquid chromatography electrospray ionization tandem mass spectrometer (LC-ESI-
121
MS/MS) with an UltiMate 3000 RSLCnano system on-line coupled to a Q Exactive mass
122
spectrometer via a Nanospray Flex ion source (Thermo Fisher Scientific, Hvidovre, Denmark) as
123
described elsewhere 24 with the following modifications: survey scans were detected at 350-1850
124
m/z and a dynamic exclusion of 45 sec was used for minimizing repetitive selection of the same
125
ions and normalized collision energy was set to 35. Peptides were eluted using 148 min linear
126
gradient, ranging from 13-50% (v/v) of solvent consisting of 0.1% (v/v) formic acid, 0.005%
127
(v/v) heptafluorobutyric acid, and 90% (v/v) acetonitrile.
128 129
2.5
Label-free quantification
130
Raw data files from the Q Exactive instrument were processed with MaxQuant v.1.5.1.2 (Max
131
Planck Institute of Biochemistry, Martinsried, Germany; Cox and Mann 2008).
132
Carboxymethylation was set as a fixed modification, oxidation of methionine was set as a
133
variable modification, and a false discovery rate of 1% was utilized. Label-free quantification
134
was determined using the label-free quantitation (LFQ) feature of MaxQuant. All other settings
135
were default. The raw data was searched against a sequence database retrieved from UniProt
136
(downloaded October 21, 2014) and comprised of the predicted open reading frames (ORFs) of
137
Sphingobium sp. BiD32. Differences in protein abundances between control and treated (with
138
BPA) were statistically evaluated using a Student's t-test with the statistical significance set at p
139
≤ 0.05. Protein abundances are represented as log2 transformed LFQ values. Identified proteins
7 ACS Paragon Plus Environment
Environmental Science & Technology
140
include proteins with a single peptide in at least one LC-MS run. Quantified proteins include
141
proteins with at leasttwo quantifiable unique peptides in at least two LC−MS runs.
142 143
2.6
Functional annotation of proteins
144
Quantified proteins more abundant in the presence of BPA were compared to a database
145
containing all proteins with a Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology
146
(KO) identifier (downloaded from UniProt 2nd of November 2013) using CLC Main Workbench
147
v 6.6.2 (CLCbio, Denmark) and The Basic Local Alignment Search Tool for proteins (BLASTp)
148
with default settings. Quantified proteins more abundant in the presence of BPA were also
149
compared to a database containing all proteins in the Patric database (downloaded from Patric
150
16th of August 2015) 26 using CLC Main Workbench v 6.6.2 (CLCbio, Denmark) and The Basic
151
Local Alignment Search Tool for proteins (BLASTp) with default settings. Matches with the
152
highest similarity scores and an E-value ≤ 1E-05 were recorded for each protein. Proteins were
153
mapped using the KO identifier from the BLASTp search and the Pathway mapping feature at
154
the KEGG (v. 68.0) website (http://www.genome.jp/KEGG/).
155 156
2.7
Pathway prediction of BPA
157
The BPA SMILES (simplified molecular-input line-entry system) line notation
158
(Oc1ccc(cc1)C(c2ccc(O)cc2)(C)C) was submitted to the EAWAG biocatalysis/biodegradation
159
database pathway prediction system (EAWAG-BBD PPS, http://eawag-bbd.ethz.ch/predict/) 27.
160 161
2.8
Sample preparation for metabolite analysis
8 ACS Paragon Plus Environment
Page 8 of 31
Page 9 of 31
Environmental Science & Technology
162
Samples for metabolite analysis were collected from log-phase growth Sphingobium sp. BiD32
163
cultures that were actively degrading BPA. Sphingobium sp. BiD32 was inoculated into R2B
164
containing either 0.5 or 2 mg/L BPA by transferring a single colony from solid media. The
165
bacteria were grown at 100 rpm and 28ºC. Samples were preserved by adding an equal volume
166
of acetonitrile, vortexing, and centrifuging (10 min, 10,000 ×g) to remove particles. Decanted
167
samples were stored at 4ºC in glass vials until concentrated for use. Samples were prepared as
168
described elsewhere with minor modifications 28. To concentrate, 1 mL of the sample was
169
evaporated after centrifugation at 10,000×g, under a gentle stream of nitrogen while at 60ºC. The
170
residue was resuspended in 50 µL of mobile phase A/B (1:1;v:v (A) 10 mM ammonium formate
171
in acetonitrile with 0.1% (v/v) formic acid and (B) acetonitrile with 0.1% (v/v) formic acid.
172 173
2.9
LC/MS analyses for metabolite identification
174
An Ultra High Definition (UHD) Accurate-Mass Quadrupole Time-of-Flight (Q-TOF) coupled
175
with an ultra-HPLC system (Agilent Technologies, Inc.; Santa Clara, CA, USA) was used to
176
determine the presence of BPA and the corresponding metabolites (details in SI).
177 178
3.
Results
179
3.1
Whole genome sequencing
180
Whole genome sequencing of Sphingobium sp. BiD32 was completed to provide an organism-
181
specific genomic template for the MS-based label-free quantitative proteomics. Information
182
obtained from the assembly of the Sphingobium sp. BiD32 genome is summarized in Table 1.
183
169 out of 206 (82%) bacterial essential genes proposed by Gil et al. were identified in the
184
genome, suggesting that the genome is mostly complete (Table S1) 26.
9 ACS Paragon Plus Environment
Environmental Science & Technology
185 186
Table 1. Sphingobium sp. BiD32 whole genome assembly. Parameter
Value
Contigs
262
Reads
2,632,020
Fold coverage
65
Mean read length (bp)
118.31
N50 contig length (Kb)
33.3
Draft genome size (b)
4,788,021
Predicted coding sequences
4,670
GC content (%)
31.3
187 188
3.2
Label-free quantitative MS-based proteomics of Sphingobium sp. BiD32
189
Proteins extracted from five biological replicates of Sphingobium sp. BiD32 grown in the
190
presence and absence of BPA were subjected to shotgun-proteomics and label-free
191
quantification. A total of 2180 proteins were identified (full proteins lists available in
192
SI_proteinlist.xlsx), corresponding to approximately 47% of the predicted CDS of Sphingobium
193
sp. BiD32. The number of proteins quantified, differentially abundant in response to the
194
presence/absence of BPA (p ≤ 0.05), and differentially abundant by more than twofold (-1.0 ≤
195
log2 (fold abundance) ≥ 1.0, p ≤ 0.05) in response to BPA are summarized in Table 2. Proteins
196
differentially abundant by more than twofold are shown by the shaded regions of Figure 2 for
197
both the cytosolic and membrane proteins; 43 proteins were upregulated in the cells exposed to
10 ACS Paragon Plus Environment
Page 10 of 31
Page 11 of 31
Environmental Science & Technology
198
BPA, while 48 proteins were downregulated in the presence of BPA. The upregulated proteins
199
were further evaluated to identify those potentially involved in BPA degradation.
200 201
Table 2. Sphingobium sp. BiD32 proteomic analysis. Proteins
Value
Identified
2155
Quantified
1174
Differentially abundant
184
in response to BPA Differentially abundant
90
by more than twofold 202
identified, protein was measured in 1 of 5 replicate incubations
203
quantified, protein was measured in 2 or more of 5 replication incubations
204 b) 4
3.5
3.5
3
3 -Log10(p-value)
-Log10(p-value)
a) 4
2.5 2 1.5
2.5 2 1.5
1
1
0.5
0.5 0
0 -5
-4
-3
-2 -1 0 1 2 log2(fold abundance)
3
4
-5
5
-4
-3
-2
-1
0
1
2
log2(fold abundance)
11 ACS Paragon Plus Environment
3
4
5
Environmental Science & Technology
205
Figure 2. Volcano plot of a) cytosolic and b) membrane quantified proteins. Statistically
206
significant (p ≤ 0.05) differentially abundant proteins (-1 ≤ log2 (fold abundance) ≥ 1) are shaded
207
gray.
208 209
3.3
Functional annotation of proteins
210
Proteins were functionally annotated by comparing them to those with KO identifiers and those
211
in the Patric database. Twenty-seven (27) of the 43 upregulated proteins matched to proteins
212
with KO identifiers. These proteins were annotated to the following categories: metabolism,
213
genetic information processing, environmental information processing, cellular processes, and
214
organismal systems. Metabolism included proteins related to xenobiotics biodegradation and
215
metabolism (N1MT47, N1MMX0, N1MNM9, N1MWA7, N1MSB0, and N1MS32), energy
216
metabolism (N1ML71, N1MN00, N1MTK7, N1MRW4, N1MPN2, and N1MLU6), lipid
217
metabolism (N1MTK7), amino acid metabolism (N1MF55, N1MFT4, N1MFX4, N1MK88,
218
N1ML71, N1MMX0, N1MPU3, N1MQ53, N1MS39, N1MWX7, and N1MRW4), metabolism
219
of other amino acids (N1MPU3), metabolism of cofactors and vitamins (N1MH82 and
220
N1MPU3), glycan biosynthesis and metabolism (N1MQ53), biosynthesis of other secondary
221
metabolites (N1MTJ4), metabolism of terpenoids and polyketides (N1MTK7), and carbohydrate
222
metabolism (N1MJQ9, N1MKJ6, N1ML71, N1MN00, N1MNR5, N1MTJ4, N1MTK7, and
223
N1MWX7). Genetic information processing included proteins related to translation (N1MH82).
224
Environmental information processing included proteins related to membrane transport
225
(N1MGE8) and signal transduction (N1MTJ4). Cellular processes included transport and
226
catabolism (N1ML71) and cell growth and death (N1MPU5). Organismal systems included
227
endocrine system (N1MN00). Functional information for four additional proteins was
228
determined by matching upregulated proteins to proteins in the Patric database. These proteins 12 ACS Paragon Plus Environment
Page 12 of 31
Page 13 of 31
Environmental Science & Technology
229
had Enzyme Class (EC) assignments of glycosidases (N1MKW1), enoyl-(acyl-carrier-protein)
230
reductase (NADH) (N1MSL9), and glutamate dehydrogenase (N1MUY1). A protein without an
231
EC assignment, had a Gene Ontology (GO) assignment of translation release factor activity and
232
translation termination (N1MIL6).
233 234
The xenobiotic biodegradation and metabolism orthology were given particular focus as these
235
were most likely to be involved with BPA degradation. These accounted for more of the
236
functionally annotated upregulated proteins than any other orthology. The level of homology
237
between the proteins annotated to xenobiotic biodegradation and metabolism and the entries in
238
the KO database varied, ranging from moderate (44% sequence identity) to high (88% sequence
239
identity) (Table 3). These proteins were characterized as dehydrogenases, a dioxygenase, a
240
hydratase, a hydroxylase, and a cycloisomerase and were related functionally to the degradation
241
of aromatic contaminants (Table 3). The presence of BPA resulted in a three-to-six-fold
242
increased abundance (log2 of 1.6 to 2.7) of these proteins. While two of the proteins were
243
exclusively associated with ‘benzoate degradation’ (N1MT47 and N1MWA7), the remaining
244
four were associated with degradation of multiple aromatic xenobiotics, including polycyclic
245
aromatic hydrocarbons, toluene, and aminobenzoate. Additionally, the hydroxymuconate-6-
246
semialdehyde (N1MMX0), which was almost four-fold upregulated in samples exposed to BPA
247
relative to the control, was found to be highly homologous (98% sequence identity, BLAST
248
against the UniProtKB database) to a predicted oxidoreductase of the aromatic hydrocarbon
249
degrading bacterium Sphingobium sp. Ant17.
250
13 ACS Paragon Plus Environment
Environmental Science & Technology
Page 14 of 31
K10220
■
1.8
Q2G4H6
86
567
K10219
■
■
2.0
Q2G4H5
85
497
K04101
■
■
■
2.7
B8GZG3
78
184
K00481
■
1.6
Q2G4H4
78
225
K04100
■
■
■
1.6
Q39HW9
44
179
K01856
■
■
proJ
hydratase
N1MMX0
4-carboxy-2-
proC
hydroxymuconate-
Polycyclic aromatic hydrocarbon degradation
638
Toluene degradation
88
Chlorocyclohexane and chlorobenzene degradation
Q2G4H3
Aminobenzoate degradation
KEGG Orthology
2.5
Benzoate degradation
Greatest bit score
4-oxalomesaconate
N1MT47
Accession (Hit)
Greatest identity %
metabolism proteins. Log2 (Fold Abundance)
252
Protein name assigned in genome annotation
Table 3. Proteins upregulated in response to BPA related to xenobiotic degradation and
UniProtKB accession number and gene
251
6-semialdehyde dehydrogenase N1MNM9
Protocatechuate
proB
4,5-dioxygenase beta chain
N1MWA7
p-hydroxybenzoate hydroxylase
N1MSB0
Protocatechuate
proA
4,5-dioxygenase alpha chain
N1MS32
Muconate
■
■
cycloisomerase
253 254
Four (4) of the proteins categorized as xenobiotic biodegradation and metabolism (N1MMX0,
255
N1MNM9, N1MSB0, and N1MT47) were homologous respectively to the genes proC, proB,
256
proA, and proJ, which are part of a proposed protocatechuate degradation pathway for BPA
14 ACS Paragon Plus Environment
Page 15 of 31
Environmental Science & Technology
257
(Table 4). Previous descriptions of the protocatechuate degradation pathway have suggested two
258
major structural organizations of the involved genes: the Comamonas-type, in which the genes of
259
the pathway are confined to a single operon, and the Sphingobium-type, in which the genes
260
constitute several transcriptional units 29. The structural organization of the protocatechuate
261
degradation genes of Sphingobium sp. BiD32 (Figure 3) appeared to be analogous to the
262
Sphingobium-type (SI_operon.xlsx) and the genes were confined to possibly three operons:
263
proFNTEL, proCBAJXR, and proKUI (Figure 3).
264 265
The proteins related to the proCBAJ operon were all highly upregulated (~four-fold) in response
266
to BPA exposure. Genes proB and proA encode for the small and large subunit of the
267
protocatechuate 4,5-dioxygenase, respectively. This enzyme mediates the first reaction in the
268
protocatechuate 4,5-cleavage pathway, catalyzing the oxygenolytic fission of protocatechuate
269
and its conversion to 4-carboxy-2-hydroxymuconate-6-semialdehyde dehydrogenase (CHMS).
270
The proC gene of Sphingobium sp. BiD32 encoded for a protein that is highly homologous (~78-
271
85% sequence identity) to the CHMS genes of other protocatechuate degraders (2811, ligC, and
272
pmdC) (SI_operon.xlsx) 29. CHMS catalyzes the oxidation of the hemiacetal form of CHMS
273
resulting in 2-pyrone-4,6-dicarboxylate (PDC) 29. PDC is then converted to the kenol-enol
274
tautomers of 4-oxalomesaconate (OMA) by the enzyme 2-pyrone-4,6-dicarboxylate hydrolase
275
(PDC hydrolase), which is encoded for by proI in Sphingobium sp. BiD32. proU is believed to
276
encode for an OMA tautomerase that catalyzes the tautomerization of the enol form of OMA to
277
the keto form. Previous studies performed on OMA tautomerase mutants indicates that this gene
278
is essential for OMA conversion 30. OMA is then converted to 4-carboxy-4-hydroxy-2-
279
oxoadipate (CHA) by the enzyme 4-oxalomesaconate hydratase (OMA hydratase), which is
15 ACS Paragon Plus Environment
Environmental Science & Technology
Page 16 of 31
280
encoded for by proJ in Sphingobium sp. BiD32. proJ is highly homologous to ligJ (~85%
281
sequence identity) from Sphingobium sp. strain SYK-6, an organism with the ability to grow
282
using a wide variety of lignin-derived mono- and diaromatic compounds 31. CHA is converted to
283
pyruvate and oxaloacetate catalyzed by the enzyme 4-carboxy-4-hydroxy-2-oxoadipate aldolase
284
(CHA aldolase), which is encoded for by proK in Sphingobium sp. BiD32. proK shows a high
285
degree of sequence similarity (~91% sequence identity) with the CHA aldolase of Sphingomonas
286
sp. LB126 (fldZ), an organism with the ability to grow on the polycyclic aromatic hydrocarbon
287
fluorene as the sole source of energy and carbon 32. Pyruvate and oxaloacetate are then readily
288
assimilated and used in the Krebs cycle. Upstream of proK in the genome of Sphingobium sp.
289
BiD32 is a LysR-type transcriptional regulator (proR), which likely modulates transcription of
290
the proRKUI operon through positive regulation 33.
291 292
Table 4. Gene cluster containing the protocatechuate transformation genes encoding for proteins
293
upregulated by more than twofold in response to bisphenol A and close proximity genes. Gene
UniProtKB
Protein name assigned in genome
Length
log2
accession
annotation
(bp)
(fold
number EBBID32_29090
p-value
abundance)
N1MMW4
p-hydroxybenzoate hydroxylase
404
NQ
N1MNM4
Uncharacterized protein
186
NQ
N1MSA4
Major facilitator superfamily
435
ND
290
ND
(proF) EBBID32_29100 (proN) EBBID32_29110 (proT) EBBID32_29120 (proE)
MFS_1 N1MT38
Transcriptional regulator, AraC family
16 ACS Paragon Plus Environment
Page 17 of 31
Environmental Science & Technology
EBBID32_29130
N1MP91
Uncharacterized protein
75
ND
N1MMX0
4-carboxy-2-hydroxymuconate-6-
312
1.8
3.0E-04
299
2.0
2.7E-03
140
1.6
2.0E-04
2.5
1.2E-02
(proL) EBBID32_29140 (proC) EBBID32_29150
semialdehyde dehydrogenase N1MNM9
(proB) EBBID32_29160
beta chain N1MSB0
(proA) EBBID32_29170
Protocatechuate 4,5-dioxygenase
Protocatechuate 4,5-dioxygenase alpha chain
N1MT47
4-oxalomesaconate hydratase
341
N1MP96
3-hydroxyisobutyrate
315
NQ
383
NQ
224
NQ
(proJ) EBBID32_29180 (proX)
dehydrogenase and related betahydroxyacid dehydrogenases
EBBID32_29190
N1MMX5
(proR) EBBID32_29200
Transcriptional regulator, LysR family
N1MNN2
(proK)
Putative siderophore biosynthesis protein,related to 2demethylmenaquinone methyltransferase
EBBID32_29210
N1MSB6
2-methylaconitate isomerase
351
NQ
N1MT55
2-pyrone-4,6-dicarboxylic acid
302
NQ
(proU) EBBID32_29220 (proI)
hydrolase
EBBID32_29230
N1MPA3
Carboxylic ester hydrolase
487
ND
EBBID32_29240
N1MMY1
Uncharacterized protein
117
ND
EBBID32_29250
N1MNN8
Uncharacterized protein
37
ND
294
NQ, not quantified; proteins were not quantified, but were identified
295
ND, not detected; protein was not identified or quantified 17 ACS Paragon Plus Environment
Environmental Science & Technology
296
297 298
Figure 3. Gene organization of the proposed protocatechuate operons (proCBAJ and proRKUI)
299
and close proximity genes on contig 143 (CAVK010000143.1) of Sphingobium sp. BiD32.
300
Proposed gene names for Sphingobium sp. BiD32 were based upon their similarity to genes of
301
the protocatechuate pathway from other organisms (SI_operon.xlsx). A black arrow indicates
302
that proteins for the corresponding gene were more abundant by more than twofold in response
303
to BPA exposure (Table 4), a grey arrow indicates that the protein was detected, but not
304
quantified, and a white arrow indicates that the gene was detected but the protein was not present
305
in the proteome.
306 307
Located outside the protocatechuate gene cluster, another protein (N1MN05) was functionally
308
annotated (K00100) to bisphenol degradation (SI_proteinlist.xlsx). This protein was found to be
309
more abundant in cells exposed to BPA relative to the control, though by less than twofold (log2
310
= 0.72). It was functionally annotated to the transformation of several bisphenol metabolites,
311
including the transformation of bis(4-hydroxyphenyl)-methanol to 4,4’-dihydroxybenzophenone
312
and transformation of 1-(4’hydroxyphenyl)-ethanol to 4’-hydroxyacetophenone.
313 314
Thirty (30) out of 48 of the downregulated proteins matched to proteins with a KO identifier.
315
These proteins were functionally annotated to the following categories: metabolism, genetic
316
information processing, environmental information processing, cellular processes, and human 18 ACS Paragon Plus Environment
Page 18 of 31
Page 19 of 31
Environmental Science & Technology
317
diseases. Metabolism included proteins related to carbohydrate metabolism (N1MKW5,
318
N1MJA6, N1MGN1, N1MG01, N1MIP5, N1MKD6, N1MHY2, N1MSX0, and N1MLE2),
319
nucleotide metabolism (N1MHY2, N1MJA6, and N1MMB5), energy metabolism (N1MKW5,
320
N1MSC3, N1MG01, N1MIP5, N1MSX0, and N1MLE2), lipid metabolism (N1MLC6 and
321
N1MKW5), amino acid metabolism (N1MWN4, N1MIR7, N1MKD6, N1MQT8, N1MN31,
322
N1MS44), metabolism of other amino acids (N1MLE2), glycan biosynthesis and metabolism
323
(N1MI43), metabolism of terpenoids and polyketides (N1MKW5), metabolism of cofactors and
324
vitamins (N1MRQ2, N1MT40, N1MNZ6, N1MIR7, and N1MY92) and biosynthesis of other
325
secondary metabolites (N1MPK8). Genetic information processing included proteins related to
326
translation (N1MMK2, N1MMZ3, N1MNQ3 and N1MQ71) and folding, sorting, and
327
degradation (N1MTE1). Environmental information processing included proteins related to
328
membrane transport (N1MSV8), cellular processes included transport and catabolism
329
(N1MLE2), and human diseases included drug resistance (N1MPK8 and N1MI43) and infectious
330
diseases (N1MJ91).
331 332
3.4
Metabolite identification
333
The parent compound (BPA) and one metabolite (BPA-M; 1,2-Bis(4-hydroxyphenyl)-2-
334
propanol) were detected during metabolite analysis (Table 5). BPA-M was detected in the
335
samples with an initial BPA concentration of 2 mg/L, but not in the samples with an initial BPA
336
concentration of 0.5 mg/L. This suggests that BPA-M did not reach concentrations above the
337
method detection limits at the lower initial BPA concentration. Pure chemicals of the
338
intermediate were not available so that standards for quantification could not be prepared.
339
19 ACS Paragon Plus Environment
Environmental Science & Technology
340
Page 20 of 31
Table 5. BPA metabolites identified during BPA degradation by Sphingobium sp. BiD32. Identified BPA metabolites
Retention time (min)
m/z
BPA
20.6
227.10775
15
243.10267
BPA-M (1,2-Bis(4-hydroxyphenyl)-2propanol)
341 342
3.5
Catabolic pathway of bisphenol A
343
The catabolic pathway of BPA by Sphingobium sp. BiD32 was predicted using the online
344
prediction tool EAWAG-BBD PPS (Figure 4, BPA to BPA-M1). It was predicted that BPA-M
345
was a potential product of the first step in BPA degradation (Figure 4, top panel). More than one
346
pathway was predicted by EAWAG-BBD PPS, however metabolites and enzymes were observed
347
only for the pathway shown in Figure 4. Smaller compounds predicted as intermediates are
348
unlikely to be detected even if present as they have molecular masses too small for detection by
349
the HPLC-TOF method used for metabolite analysis, and these metabolites are likely to be
350
readily transformed.
351 352
Enzymes predicted to be involved in the transformation of BPA to BPA-M were compared to the
353
proteins upregulated in Sphingobium sp. BiD32. This transformation involves substitution by a
354
hydroxyl group. The top two matches to enzymes predicted by the EAWAG-BBD PPS were
355
N1MMX0 (83% identity similarity to Q9KWL3) and N1MWA7 (68% identity similarity to
356
Q03298). N1MWA7, identified as a p-hydroxybenzoate hydroxylase by genome annotation
357
(Table 3), was identified as a protein likely involved in the degradation of BPA by Sphingobium
20 ACS Paragon Plus Environment
Page 21 of 31
Environmental Science & Technology
358
sp. BiD32 (Figure 4). Following hydroxylation of BPA to BPA-M, it was predicted that BPA-M
359
was dehydrated forming 4,4'-dihydroxy-alpha-methylstilbene (BPA-M1). No enzyme catalyzing
360
this step is known, and hence no candidate was found for BiD32. The middle steps (BPA-M2 to
361
BPA-M5), suggested that BPA-M1 is broken into two aromatic constituents (BPA-M2 and BPA-
362
M4) based upon proteomics data. BPA-M2 is oxidized to BPA-M3 and BPA-M4 is processed
363
via the protocatechuate pathway. It remains unclear if BPA-M3 is processed by Sphingobium sp.
364
BiD32, though a previous study proposed the demethylation of BPA-M3 resulting in 4-
365
hydroxybenzaldehyde, which was further converted to benzoic acid and eventually to carbon
366
dioxide and water 34. The bottom part of the pathway (BPA-M5 to BPA-M12) represents the
367
protocatechuate pathway described from other organisms, and the Sphingobium sp. BiD32
368
enzymes catalyzing each step.
369
21 ACS Paragon Plus Environment
Environmental Science & Technology
370 371
Figure 4. Proposed pathway for BPA degradation by Sphingobium sp. BiD32 based on measured
372
metabolites (bold) and metabolites and proteins predicted to be involved in BPA degradation by
373
the EAWAG Biocatalysis/Biodegradation Pathway Prediction System (top panel). The
22 ACS Paragon Plus Environment
Page 22 of 31
Page 23 of 31
Environmental Science & Technology
374
degradation pathway of BPA metabolites (bottom panels) was proposed based upon previous
375
findings35. Solid arrows indicate identification of a candidate enzyme upregulated in
376
Sphingobium sp. BiD32 by more than twofold in response to BPA relative to the control. Grey
377
arrows indicate identification of a candidate enzyme only quantified in the proteome of BPA-
378
exposed cells. Thin arrows indicate that the step is not catalyzed enzymatically. Metabolites are
379
as follows: BPA, bisphenol A; BPA-M, 1,2-bis(4-hydroxyphenyl)-2-propanol; BPA-M1, 4,4’-
380
(1E)-prop-1-ene-1,2-diyldiphenol; BPA-M2, 1-(4’-hydroxyphenyl)-ethanol; BPA-M3, 4’-
381
hydroxyacetophenone; BPA-M4, 4-hydroxybenzoate; BPA-M5, protocatechuate (PCA); BPA-
382
M6, 4-carboxy-2-hydroxymuconate-6-semialdehyde (CHMS); BPA-M7, 4-carboxy-2-
383
hydroxymuconate-6-semialdehyde (CHMS hemiacetal form); BPA-M8, 2-pyrone-4,6-
384
dicarboxylate (PDC); BPA-M9, 4-oxalomesacconate (OMA) enol form; BPA-M10, OMA
385
ketoform; BPA-M11, OMA ketoform; BPA-M12, 4-carboxy-4-hydroxy-2-oxoadipate (CHA).
386 387
4.
Discussion
388
In this study, the metabolism of BPA by Sphingobium sp. BiD32 was investigated using three
389
different tiers of biochemistry (genome, proteome, and metabolome) and was successfully
390
employed to identify previously undocumented proteins putatively involved in bacterial BPA
391
degradation. The overarching goal was to identify genes involved in BPA degradation for use as
392
potential biomarkers in future work. The gene involved in the initial transformation of BPA may
393
also enable identification of other new BPA-degrading bacteria.
394 395
A p-hydroxybenzoate hydroxylase (N1MWA7) was identified as involved in BPA degradation
396
based on proteomics results and metabolite analysis. This protein was significantly more
397
abundant in the presence of BPA (Table 3), was functionally annotated to proteins involved in 23 ACS Paragon Plus Environment
Environmental Science & Technology
398
xenobiotics biodegradation and metabolism (Table 3), and closely matched enzymes predicted to
399
be involved in the transformation of BPA to BPA-M (Figure 4). Upregulated Sphingobium sp.
400
BiD32 proteins that were functionally annotated to xenobiotics biodegradation and metabolism
401
were distantly related to each other.
402 403
Proteins involved in protocatechuate transformation are also likely to be involved in the BPA
404
degradation pathway (N1MMX0, N1MNM9, N1MSB0, and N1MT47). These proteins were
405
more abundant in the presence of BPA (Table 4), were functionally annotated to proteins
406
involved in xenobiotics biodegradation and metabolism (Table 3), and matched genes previously
407
identified to be involved in protocatechuate degradation in a variety of organisms including
408
Sphingobium sp. strain SYK-6 36, Sphingomonas sp. strain LB126 37, Novosphingobium
409
aromaticivorans DSM 12444 37, Comamonas sp. strain E6 35, Comamonas testosterone BR6020
410
38
411
protocatechuate degradation genes in Sphingobium sp. BiD32 most notably resembles that of
412
Sphingobium sp. strain SYK-6, a bacterium capable of degrading a wide variety of lignin-derived
413
biaryls and monoaryls 36.
, Pseudomonas straminea NGJ1 39, and Arthrobacter keyseri 12B 40. The organization of the
414 415
While several studies have documented microbial degradation of BPA and predicted BPA
416
degradation pathways, few proteins involved in BPA degradation have been reported 41. Included
417
amongst these was a cytochrome P450 monooxygenase reported to catalyze the transformation
418
of BPA to BPA-M by Sphingomonas sp. AO1 42. This cytochrome P450 monooxygenase was
419
found to be involved in BPA degradation by using a cytochrome P450 inhibitor 43 and then
420
confirmed by purifying the components of the cytochrome P450 monooxygenase system and
24 ACS Paragon Plus Environment
Page 24 of 31
Page 25 of 31
Environmental Science & Technology
421
testing their ability to degrade BPA 42. The Sphingomonas sp. AO1 cytochrome P450 (B1Q2N5)
422
was related to a Sphingobium sp. BiD32 putative cytochrome P450 hydroxylase (N1MI62) that
423
was identified in the genome, but not detected based on the proteomics portion of this study
424
(Figure S1; 30% identity match; E-value = 5.1e-45). Also, a Sphingomonas sp. AO1 ferredoxin
425
that is part of the cytochrome P450 monooxygenase system (B1Q2N4) was closely related to a
426
Sphingobium sp. BiD32 ferredoxin gene (N1MLS9) detected in the genome and proteome
427
(Figure S1; 47% identity match; E-value = 4.4e-23), but not differentially expressed. Lack of
428
increased proteome expression of these enzymes by Spingobium sp. BiD32 suggested that
429
differing pathways may be involved in BPA degradation by these two organisms.
430 431
Other enzymes previously reported to be involved in BPA degradation included a second
432
cytochrome P450, an ammonia monooxygenase, and an extracelluar laccase. The cytochrome
433
P450 was reported to catalyze the ipso substitution as the initial step in BPA degradation by
434
Sphingomonas sp. TTNP3 44. An ammonia monooxygenase in Nitrosomonas europaea 45 and an
435
extracellular laccase in Pseudomonas sp. LBC1 46 were also reported to be involved in BPA
436
degradation, though specific degradation steps mediated were not identified. No proteins similar
437
to these three were identified in the proteome of Sphingobium sp. BiD32.
438 439
Investigation of the metabolome resulted in the identification of BPA-M (Table 5). BPA-M has
440
also been detected during BPA degradation by other isolated bacteria 41, generally resulting in
441
complete mineralization. For example, BPA degradation by Sphingomonas sp. MV1 resulted in
442
the transformation of 85% of the BPA to BPA-M47. BPA-M was further degraded and finally
443
oxidized to carbon dioxide and cell biomass, as determined by HPLC. BPA degradation by
25 ACS Paragon Plus Environment
Environmental Science & Technology
444
unidentified strain WH1 followed the main degradation pathway that occurred in Sphingomonas
445
sp. MV1, resulting in transformation to BPA-M and, ultimately, mineralization 48. BPA-M was
446
also detected during BPA degradation by Sphingomonas sp. AO1, though in that study further
447
investigation of the degradation products was not conducted, as it was not the focus of the study
448
42
449
Sphingobium sp. BiD32; other metabolites were not detected suggesting that they were either too
450
small to be detected or they were transformed quickly, resulting in too low concentrations to be
451
detected. While the same metabolite, BPA-M, was identified during BPA degradation by
452
Sphingobium sp. BiD32 and by Sphingomonas sp. AO1, different enzymes were predicted to be
453
involved in the transformation 42 suggesting that the mechanisms responsible for BPA
454
degradation by these organisms likely differ.
. The current study detected BPA-M during degradation of high concentrations of BPA by
455 456
This is the first study using a combination of genomics, proteomics, and metabolite analysis to
457
identify BPA metabolites and enzymes responsible for BPA degradation. Previous studies
458
identified enzymes responsible for the initial transformation of BPA using various inhibitors.
459
This study identified a novel enzyme responsible for the initial transformation of BPA. It also
460
identified a novel BPA transformation pathway with the increased abundance in the
461
protocatechuate transformation genes. The results further suggest how Sphingobium sp. BiD32
462
degrades BPA beyond 4-hydroxybenzoate. Additional research is needed to elucidate the
463
function of multiple enzymes in BPA degradation. Identification of these enzymes may allow for
464
future monitoring of bioaugmented bacteria in WWTPs and identifying other similar BPA-
465
degrading bacteria.
466
26 ACS Paragon Plus Environment
Page 26 of 31
Page 27 of 31
Environmental Science & Technology
467
Acknowledgements
468
This work was funded by the Danish Research Council for Strategic Research via the Research
469
Centre “EcoDesign-MBR” (Grant No. 26-03-0250) and the National Science Foundation (NSF
470
CBET-0829132). Support was provided to N.A. Zhou by the Valle Scholarship and
471
Scandinavian Exchange and the King County Technology Transfer Fellowship Programs at the
472
University of Washington.
473 474
Supporting Information Available
475
Tables and Figures addressing the full protein quantitation lists, protocatechuate transformation
476
gene cluster, the gene sequences of the proteins of interest, and the phylogenetic relationship of
477
identified proteins of interest. This material is available free of charge via the Internet at
478
http://pubs.acs.org.
479
27 ACS Paragon Plus Environment
Environmental Science & Technology
480
References
481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529
(1) (2) (3) (4)
(5)
(6)
(7) (8)
(9)
(10)
(11)
(12)
(13) (14)
(15) (16)
(17)
Bisphenol A (BPA) Action Plan Summary | Existing Chemicals | OPPT | US EPA http://www.epa.gov/oppt/existingchemicals/pubs/actionplans/bpa.html (accessed Apr 28, 2014). Tsai, W.-T. Human health risk on environmental exposure to bisphenol-A: A review. J. Environ. Sci. Health Pt. C-Environ. Carcinog. Ecotoxicol. Rev. 2006, 24 (2), 225–255. Flint, S.; Markle, T.; Thompson, S.; Wallace, E. Bisphenol A exposure, effects, and policy: A wildlife perspective. J. Environ. Manage. 2012, 104, 19–34. Baek, H. J.; Park, M. H.; Lee, Y. D.; Kim, H. B. Effect of in vitro xenoestrogens on steroidogenesis in mature female fish, Chasmichthys dolichognathus. Fish Physiol. Biochem. 2003, 28 (1-4), 413–414. Oehlmann, J.; Schulte-Oehlmann, U.; Bachmann, J.; Oetken, M.; Lutz, I.; Kloas, W.; Ternes, T. A. Bisphenol A induces superfeminization in the ramshorn snail Marisa cornuarietis (Gastropoda : Prosobranchia) at environmentally relevant concentrations. Environ. Health Perspect. 2006, 114, 127–133. Sohoni, P.; Tyler, C. R.; Hurd, K.; Caunter, J.; Hetheridge, M.; Williams, T.; Woods, C.; Evans, M.; Toy, R.; Gargas, M.; et al. Reproductive effects of long-term exposure to bisphenol a in the fathead minnow (Pimephales promelas). Environ. Sci. Technol. 2001, 35 (14), 2917–2925. Lemos, M. F. L.; van Gestel, C. A. M.; Soares, A. M. V. M. Endocrine disruption in a terrestrial isopod under exposure to bisphenol A and vinclozolin. J. Soils Sediments 2009, 9 (5), 492–500. Kolpin, D. W.; Furlong, E. T.; Meyer, M. T.; Thurman, E. M.; Zaugg, S. D.; Barber, L. B.; Buxton, H. T. Pharmaceuticals, Hormones, and Other Organic Wastewater Contaminants in U.S. Streams, 1999−2000: A National Reconnaissance. Environ. Sci. Technol. 2002, 36 (6), 1202– 1211. Rudel, R. A.; Melly, S. J.; Geno, P. W.; Sun, G.; Brody, J. G. Identification of alkylphenols and other estrogenic phenolic compounds in wastewater, septage, and groundwater on Cape Cod, Massachusetts. Environ. Sci. Technol. 1998, 32 (7), 861–869. Zhou, N. A.; Lutovsky, A. C.; Andaker, G. L.; Ferguson, J. F.; Gough, H. L. Kinetics modeling predicts bioaugmentation with Sphingomonad cultures as a viable technology for enhanced pharmaceutical and personal care products removal during wastewater treatment. Bioresour. Technol. 2014, 166, 158–167. Meade, M. J.; Waddell, R. L.; Callahan, T. M. Soil bacteria Pseudomonas putida and Alcaligenes xylosoxidans subsp denitrificans inactivate triclosan in liquid and solid substrates. FEMS Microbiol. Lett. 2001, 204 (1), 45–48. Zhang, C.; Zeng, G.; Yuan, L.; Yu, J.; Li, J.; Huang, G.; Xi, B.; Liu, H. Aerobic degradation of bisphenol A by Achromobacter xylosoxidans strain B-16 isolated from compost leachate of municipal solid waste. Chemosphere 2007, 68 (1), 181–190. Kang, J. H.; Kondo, F. Bisphenol a degradation by bacteria isolated from river water. Arch. Environ. Contam. Toxicol. 2002, 43 (3), 265–269. Pauwels, B.; Wille, K.; Noppe, H.; De Brabander, H.; van de Wiele, T.; Verstraete, W.; Boon, N. 17 alpha-ethinylestradiol cometabolism by bacteria degrading estrone, 17 beta-estradiol and estriol. Biodegradation 2008, 19 (5), 683–693. Murdoch, R. W.; Hay, A. G. Formation of catechols via removal of acid side chains from ibuprofen and related aromatic acids. Appl. Environ. Microbiol. 2005, 71 (10), 6121–6125. Zhou, N. A.; Lutovsky, A. C.; Andaker, G. L.; Gough, H. L.; Ferguson, J. F. Cultivation and characterization of bacterial isolates capable of degrading pharmaceutical and personal care products for improved removal in activated sludge wastewater treatment. Biodegradation 2013, 24 (6), 813–827. Ballesteros-Gomez, A.; Ruiz, F.-J.; Rubio, S.; Perez-Bendito, D. Determination of bisphenols A and F and their diglycidyl ethers in wastewater and river water by coacervative extraction and liquid chromatography-fluorimetry. Anal. Chim. Acta 2007, 603 (1), 51–59. 28 ACS Paragon Plus Environment
Page 28 of 31
Page 29 of 31
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 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579
Environmental Science & Technology
(18)
(19)
(20) (21)
(22)
(23)
(24)
(25) (26)
(27) (28)
(29) (30)
(31)
(32)
(33) (34)
(35)
Duong, C. N.; Ra, J. S.; Cho, J.; Kim, S. D.; Choi, H. K.; Park, J.-H.; Kim, K. W.; Inam, E.; Kim, S. D. Estrogenic chemicals and estrogenicity in river waters of South Korea and seven Asian countries. Chemosphere 2010, 78 (3), 286–293. Almeida, B.; Kjeldal, H.; Lolas, I.; Knudsen, A. D.; Carvalho, G.; Nielsen, K. L.; Crespo, M. T. B.; Stensballe, A.; Nielsen, J. L. Quantitative proteomic analysis of ibuprofen-degrading Patulibacter sp. strain I11. Biodegradation 2013, 24 (5), 615-630. Li, Z.; Nandakumar, R.; Madayiputhiya, N.; Li, X. Proteomic Analysis of 17 beta-Estradiol Degradation by Stenotrophomonas maltophilia. Environ. Sci. Technol. 2012, 46 (11), 5947–5955. Collodoro, M.; Lemaire, P.; Eppe, G.; Bertrand, V.; Dobson, R.; Mazzucchelli, G.; Widart, J.; De Pauw, E.; De Pauw-Gillet, M.-C. Identification and quantification of concentration-dependent biomarkers in MCF-7/BOS cells exposed to 17 beta-estradiol by 2-D DIGE and label-free proteomics. J. Proteomics 2012, 75 (14), 4555–4569. Overbeek, R.; Begley, T.; Butler, R. M.; Choudhuri, J. V.; Chuang, H. Y.; Cohoon, M.; de CrecyLagard, V.; Diaz, N.; Disz, T.; Edwards, R.; et al. The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucleic Acids Res. 2005, 33 (17), 5691–5702. Reasoner, D. J.; Geldreich, E. E. A New Medium for the Enumeration and Sub Culture of Bacteria from Potable Water. Abstracts of the Annual Meeting of the American Society for Microbiology 1979, No. 79, 180. Kjeldal, H.; Pell, L.; Pommerening-Röser, A.; Nielsen, J. L. Influence of p-cresol on the proteome of the autotrophic nitrifying bacterium Nitrosomonas eutropha C91. Arch Microbiol 2014, 196 (7), 497–511. Cox, J.; Mann, M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotech 2008, 26 (12), 1367–1372. Wattam, A. R.; Abraham, D.; Dalay, O.; Disz, T. L.; Driscoll, T.; Gabbard, J. L.; Gillespie, J. J.; Gough, R.; Hix, D.; Kenyon, R.; et al. PATRIC, the bacterial bioinformatics database and analysis resource. Nucl. Acids Res. 2013, gkt1099. Gao, J.; Ellis, L. B. M.; Wackett, L. P. The University of Minnesota Biocatalysis/Biodegradation Database: improving public access. Nucleic Acids Res. 2010, 38, D488–D491. Wissenbach, D. K.; Meyer, M. R.; Remane, D.; Weber, A. A.; Maurer, H. H. Development of the first metabolite-based LC-MS n urine drug screening procedure-exemplified for antidepressants. Anal Bioanal Chem 2011, 400 (1), 79–88. Nojiri, H.; Tsuda, M.; Fukuda, M.; Kamagata, Y. Biodegradative Bacteria: How Bacteria Degrade, Survive, Adapt, and Evolve; Springer Science & Business Media, 2013. Kamimura, N.; Aoyama, T.; Yoshida, R.; Takahashi, K.; Kasai, D.; Abe, T.; Mase, K.; Katayama, Y.; Fukuda, M.; Masai, E. Characterization of the Protocatechuate 4,5-Cleavage Pathway Operon in Comamonas sp. Strain E6 and Discovery of a Novel Pathway Gene. Appl. Environ. Microbiol. 2010, 76 (24), 8093–8101. Masai, E.; Kamimura, N.; Kasai, D.; Oguchi, A.; Ankai, A.; Fukui, S.; Takahashi, M.; Yashiro, I.; Sasaki, H.; Harada, T.; et al. Complete Genome Sequence of Sphingobium sp. Strain SYK-6, a Degrader of Lignin-Derived Biaryls and Monoaryls. J. Bacteriol. 2012, 194 (2), 534–535. Wattiau, P.; Bastiaens, L.; van Herwijnen, R.; Daal, L.; Parsons, J. R.; Renard, M.-E.; Springael, D.; Cornelis, G. R. Fluorene degradation by Sphingomonas sp. LB126 proceeds through protocatechuic acid: a genetic analysis. Research in Microbiology 2001, 152 (10), 861–872. Maddocks, S. E.; Oyston, P. C. F. Structure and function of the LysR-type transcriptional regulator (LTTR) family proteins. Microbiology 2008, 154 (12), 3609–3623. Li, G.; Zu, L.; Wong, P.-K.; Hui, X.; Lu, Y.; Xiong, J.; An, T. Biodegradation and detoxification of bisphenol A with one newly-isolated strain Bacillus sp GZB: Kinetics, mechanism and estrogenic transition. Bioresour. Technol. 2012, 114, 224–230. Kamimura, N.; Aoyama, T.; Yoshida, R.; Takahashi, K.; Kasai, D.; Abe, T.; Mase, K.; Katayama, Y.; Fukuda, M.; Masai, E. Characterization of the Protocatechuate 4,5-Cleavage Pathway Operon
29 ACS Paragon Plus Environment
Environmental Science & Technology
580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619
(36)
(37)
(38)
(39)
(40) (41) (42)
(43) (44)
(45)
(46)
(47)
(48)
in Comamonas sp. Strain E6 and Discovery of a Novel Pathway Gene. Appl. Environ. Microbiol. 2010, 76 (24), 8093–8101. Masai, E.; Kamimura, N.; Kasai, D.; Oguchi, A.; Ankai, A.; Fukui, S.; Takahashi, M.; Yashiro, I.; Sasaki, H.; Harada, T.; et al. Complete genome sequence of Sphingobium sp. strain SYK-6, a degrader of lignin-derived biaryls and monoaryls. J. Bacteriol. 2012, 194 (2), 534–535. Wattiau, P.; Bastiaens, L.; van Herwijnen, R.; Daal, L.; Parsons, J. R.; Renard, M. E.; Springael, D.; Cornelis, G. R. Fluorene degradation by Sphingomonas sp. LB126 proceeds through protocatechuic acid: a genetic analysis. Res. Microbiol. 2001, 152 (10), 861–872. Providenti, M. A.; Mampel, J.; MacSween, S.; Cook, A. M.; Wyndham, R. C. Comamonas testosteroni BR6020 possesses a single genetic locus for extradiol cleavage of protocatechuate. Microbiology (Reading, Engl.) 2001, 147 (Pt 8), 2157–2167. Maruyama, K.; Miwa, M.; Tsujii, N.; Nagai, T.; Tomita, N.; Harada, T.; Sobajima, H.; Sugisaki, H. Cloning, sequencing, and expression of the gene encoding 4-hydroxy-4-methyl-2-oxoglutarate aldolase from Pseudomonas ochraceae NGJ1. Biosci. Biotechnol. Biochem. 2001, 65 (12), 2701– 2709. Eaton, R. W. Plasmid-encoded phthalate catabolic pathway in Arthrobacter keyseri 12B. J. Bacteriol. 2001, 183 (12), 3689–3703. Zhang, W.; Yin, K.; Chen, L. Bacteria-mediated bisphenol A degradation. Appl. Microbiol. Biotechnol. 2013, 97 (13), 5681–5689. Sasaki, M.; Akahira, A.; Oshiman, K. I.; Tsuchido, T.; Matsumura, Y. Purification of cytochrome P450 and ferredoxin, involved in bisphenol A degradation, from Sphingomonas sp strain AO1. Appl. Environ. Microbiol. 2005, 71 (12), 8024–8030. Sasaki, M.; Maki, J.; Oshiman, K.; Matsumura, Y.; Tsuchido, T. Biodegradation of bisphenol A by cells and cell lysate from Sphingomonas sp strain AO1. Biodegradation 2005, 16 (5), 449–459. Kolvenbach, B.; Schlaich, N.; Raoui, Z.; Prell, J.; Zuehlke, S.; Schaeffer, A.; Guengerich, F. P.; Corvini, P. F. X. Degradation pathway of bisphenol A: Does ipso substitution apply to phenols containing a quaternary alpha-carbon structure in the para position? Appl. Environ. Microbiol. 2007, 73 (15), 4776–4784. Roh, H.; Subramanya, N.; Zhao, F.; Yu, C.-P.; Sandt, J.; Chu, K.-H. Biodegradation potential of wastewater micropollutants by ammonia-oxidizing bacteria. Chemosphere 2009, 77 (8), 1084– 1089. Telke, A. A.; Kalyani, D. C.; Jadhav, U. U.; Parshetti, G. K.; Govindwar, S. P. Purification and characterization of an extracellular laccase from a Pseudomonas sp LBC1 and its application for the removal of bisphenol A. J. Mol. Catal. B-Enzym. 2009, 61 (3-4), 252–260. Spivack, J.; Leib, T. K.; Lobos, J. H. Novel pathway for bacterial metabolism of bisphenol A. Rearrangements and stilbene cleavage in bisphenol A metabolism. J. Biol. Chem. 1994, 269 (10), 7323–7329. Ronen, Z.; Abeliovich, A. Anaerobic-aerobic process for microbial degradation of tetrabromobisphenol A. Appl. Environ. Microbiol. 2000, 66 (6), 2372–2377.
30 ACS Paragon Plus Environment
Page 30 of 31
Page 31 of 31
620
Environmental Science & Technology
TOC Art:
621
31 ACS Paragon Plus Environment