Skip to main content

The OxyR and SoxR transcriptional regulators are involved in a broad oxidative stress response in Paraburkholderia xenovorans LB400



Aerobic metabolism generates reactive oxygen species that may cause critical harm to the cell. The aim of this study is the characterization of the stress responses in the model aromatic-degrading bacterium Paraburkholderia xenovorans LB400 to the oxidizing agents paraquat and H2O2.


Antioxidant genes were identified by bioinformatic methods in the genome of P. xenovorans LB400, and the phylogeny of its OxyR and SoxR transcriptional regulators were studied. Functionality of the transcriptional regulators from strain LB400 was assessed by complementation with LB400 SoxR of null mutant P. aeruginosa ΔsoxR, and the construction of P. xenovorans pIZoxyR that overexpresses OxyR. The effects of oxidizing agents on P. xenovorans were studied measuring bacterial susceptibility, survival and ROS formation after exposure to paraquat and H2O2. The effects of these oxidants on gene expression (qRT-PCR) and the proteome (LC–MS/MS) were quantified.


P. xenovorans LB400 possesses a wide repertoire of genes for the antioxidant defense including the oxyR, ahpC, ahpF, kat, trxB, dpsA and gorA genes, whose orthologous genes are regulated by the transcriptional regulator OxyR in E. coli. The LB400 genome also harbors the soxR, fumC, acnA, sodB, fpr and fldX genes, whose orthologous genes are regulated by the transcriptional regulator SoxR in E. coli. The functionality of the LB400 soxR gene was confirmed by complementation of null mutant P. aeruginosa ΔsoxR. Growth, susceptibility, and ROS formation assays revealed that LB400 cells were more susceptible to paraquat than H2O2. Transcriptional analyses indicated the upregulation of the oxyR, ahpC1, katE and ohrB genes in LB400 cells after exposure to H2O2, whereas the oxyR, fumC, ahpC1, sodB1 and ohrB genes were induced in presence of paraquat. Proteome analysis revealed that paraquat induced the oxidative stress response proteins AhpCF and DpsA, the universal stress protein UspA and the RNA chaperone CspA. Both oxidizing agents induced the Ohr protein, which is involved in organic peroxide resistance. Notably, the overexpression of the LB400 oxyR gene in P. xenovorans significantly decreased the ROS formation and the susceptibility to paraquat, suggesting a broad OxyR-regulated antioxidant response.


This study showed that P. xenovorans LB400 possess a broad range oxidative stress response, which explain the high resistance of this strain to the oxidizing compounds paraquat and H2O2.


Aerobic metabolism induces partial oxygen reduction, generating reactive oxygen species (ROS: \({\text{O}}_{2}^{ - }\), OH· and H2O2). Formation of ROS may cause critical harm to the cell. However, antioxidant mechanisms of the cell counteract the toxicity of ROS. The transcriptional regulators OxyR and SoxRS from Escherichia coli respond to oxidative stress caused by hydrogen peroxide (H2O2) and superoxide (\({\text{O}}_{2}^{ - }\)), respectively, activating the expression of diverse antioxidant enzymes (e.g., superoxide dismutase, alkyl hydroperoxide reductase, catalase). In E. coli, the OxyR transcriptional regulator is commonly activated by H2O2, while SoxR is activated by superoxide producing redox-cycling compounds, such as paraquat (PQ; methyl viologen) and phenazine methosulfate (PMS) [29]. In non-enteric Proteobacteria and Actinobacteria such as Pseudomonas and Streptomyces strains, SoxR is activated by endogenous redox-active compounds (e.g., phenazine, actinorhodin) or their precursors, and exogenous redox compounds such as paraquat [22, 30]. Therefore, SoxR mediates primarily the response to redox-cycling molecules but not to superoxide. However, the DNA-binding property of SoxR is conserved in enteric bacteria and non-enterics [80]. In non-enteric bacteria, the SoxR may regulate the antibiotic production and export, the oxidative stress response and the morphological development [22, 30, 80].

Oxidative stress and ROS accumulation in the cell may be enhanced by environmental factors, such as the presence of aromatic compounds, heavy metals and quaternary ammonium compounds [2, 3, 17, 43, 51, 61, 67]. Paraburkholderia xenovorans LB400T (previously classified as Burkholderia xenovorans LB400) is a model and versatile aromatic-degrading bacterium that has been widely studied [14]. The metabolism of aromatic compounds in strain LB400 induces general and oxidative stress and may produce toxic compounds [2, 3, 9, 16, 48, 61]. Oxidative stress is detrimental for aromatic biodegradation [15, 61]. Notably, the addition of antioxidant compounds such as α-tocoferol enhances the degradation of chlorobiphenyls by P. xenovorans LB400 [61]. Mechanisms involved in P. xenovorans LB400 response to oxidative stress caused by oxidizing agents are barely known. The aim of this study is the characterization of the stress responses of the model aromatic-degrading bacterium P. xenovorans LB400 to the oxidizing agents paraquat and H2O2.

Materials and methods


Paraquat dichloride hydrate (PQ; > 99% purity) and phenazine methosulfate (PMS) were obtained from Sigma-Aldrich (Saint Louis, MO, USA) and MP Biomedicals (Irvine, CA, USA), respectively. Hydrogen peroxide (H2O2) solution (9% v/v) was obtained from Diphem Pharma S.A. (Santiago, Chile).

Bacterial strains and culture conditions

P. xenovorans strains were cultivated in modified Luria–Bertani at 30 °C or in M9 mineral medium with trace solutions and glucose (5 mM) or biphenyl as the sole carbon and energy source [50]. For recombinant P. xenovorans strains, gentamicin (10 µg µl−1) was added. The effect of oxidizing agents on P. xenovorans growth was assessed by adding different paraquat and H2O2 concentrations to exponential-growing cells (turbidity at 600 nm of 0.5) in M9 mineral medium with glucose (5 mM). E. coli S17λpir was grown in LB medium at 37 °C. Pseudomonas aeruginosa cells were cultivated in LB or M63 minimal media at 37 °C. Gentamicin (10 µg µl−1) was added for recombinant P. aeruginosa strains. Bacterial growth in liquid media was determined by measuring turbidity at 600 nm.

Bioinformatic analyses

To search for stress response genes in P. xenovorans LB400 genome, sequences of proteins from several bacteria were retrieved from protein databases (UniProt, NCBI). For the LB400 genome analyses, two databases were employed ( and BLAST (BLASTP and TBLASTX) search tools [4] were used to compare query sequences with the LB400 genome. Swiss-Prot protein database was employed for protein searches [5]. An amino acid sequence identity ≥ 30% was used for protein identification.

Phylogenetic analysis of OxyR and SoxR transcriptional regulators

Several OxyR and SoxR amino acid sequences with experimental evidence were retrieved from the UniProtKB Swiss-Prot/TrEMBL database [6]. All the amino acid sequences were aligned, using the M-Coffee server (using the following methods: Mpcma_msa Mmafft_msa, Mdialigntx_msa, Mpoa_msa, Mmuscle_msa, Mprobcons_msa, Mt_coffee_msa) [52]. Ambiguous positions were trimmed using the gappyout strategy of the TrimAl software [11]. The MSA were evaluated using the editor AliView version 1.24 [40]. The best partitioning scheme was identified, using the program PartitionFinder version 2.1.1 [39]. A distribution of probable trees was obtained, by Bayesian Inference as implemented in MrBayes 3.2.6 [64]. Two separate runs of 250.000 generations were executed (two chains each run,sampling every 1000 generations). Visualization and editing of phylogenetic trees were performed using the FigTree v. 1.4.2 software ( The bootstrap values (percentage) > 50% were shown for each branch point.

Molecular biology techniques

Recombinant DNA techniques were performed according to standard methods [66]. Plasmid DNA was isolated with the E.Z.N.A. Plasmid Mini Kit I (Omega Bio-tek, Norcross, USA) according to manufacturer recommendations. DNA was sequenced using an ABI Prism 377 automated DNA sequencer (Applied Biosystems Inc., Foster, CA, USA).

Construction of P. aeruginosa strain ∆soxR::BxeC1217

To study the functionality of the soxR gene (BxeC1217) of strain LB400, a complementation assay was performed with P. aeruginosa ∆soxR [21] using the LB400 BxeC1217 gene. As positive control, P. aeruginosa ∆soxR was complemented with the PA2273 gene (soxR) of P. aeruginosa PA14. For the plasmid construction, the P. xenovorans LB400 BxeC1217 and P. aeruginosa PA14 PA2273 genes were cloned into the suicide vector pMQ30, along with a region upstream and a region downstream of the PA2273 gene of strain PA14. An identity of 28.6% between SoxR regulators of P. xenovorans LB400 (BxeC1217) and P. aeruginosa PA14 (PA2273 gene) was determined by a pairwise sequence alignment ( The recombinant plasmids were transferred to P. aeruginosa ∆soxR strain by conjugation using E. coli S17ʎpir.

Construction of P. xenovorans strain pIZoxyR

A recombinant P. xenovorans strain overexpressing oxyR gene was constructed, which allowed to study the protective effect of this transcriptional regulator under oxidative stress. The oxyR gene (BxeB3987) was amplified by PCR from LB400 genomic DNA using the oligonucleotides oxyR-5′ (CCTCTAGAGCGCGGCCGTCAGTTGAC) and oxyR-3′ (CCAAGCTTTGCCCTAAGGAGGTAAAACATGACCCTCACCGAACTGAAATACATC). Primer oxyR-5′ contains XbaI restriction site, while primer oxyR-3′ contains HindIII restriction site. The DNA fragment was cloned in the broad-host range plasmid pIZ1016, a derivative of vector pBBR1MCS-5 [35], which harbors a gentamicin resistance marker. For conjugal transfers, E. coli S-17λ-pir was used as donor strain. After conjugation, clones were selected in M9 mineral medium with biphenyl crystals and gentamicin (10 µg µl−1). P. xenovorans harboring the plasmid vector was constructed as a control strain.

Morphological characterization of macrocolonies of P. aeruginosa strains

P. aeruginosa strains were grown for 7 h; thereafter, the turbidity at 600 nm was standarized. The cultures (10-μL) were plated on tryptone agar (1% w/v) in absence and presence of the antibiotic PMS (600 µM), which generates superoxide. Plates were incubated at room temperature and the growth was monitored using a high-resolution scanner (EPSON, 600 dpi) [24].

Susceptibility to oxidizing agents

The effect of oxidizing agents on the growth of P. xenovorans strains was determined by using a disk diffusion assay [21]. P. xenovorans was grown in modified LB medium in absence or presence of antibiotics for 16 h at 30 °C. 100 μl-culture were added to 4 ml of melted soft agar (1% tryptone, 0.5% agar), then plated on 1% tryptone agar plates. For P. aeruginosa strains, cells were grown in LB medium in presence of gentamicin for 14 h at 37 °C. Aliquots of oxidizing agents (15 μl) were deposited in 6 mm-diffusion disks (Whatman) at the required concentration onto the agar. P. xenovorans plates were incubated for 24 h at 30 °C. The growth inhibition zones were recorded and determined using the ImageJ software ( The diameter of the Petri dish (85 mm) was established as a length reference. The oxidizing agents employed were hydrogen peroxide and paraquat (10 and 20 mM). Paraquat is a redox-cycling compound that constitutes a continuous source of superoxide radical in the cytoplasm, which is able to activate the SoxR transcriptional regulator [23].

Reactive oxygen species quantification

To measure hydroxyl radical, the fluorescent probe 3′-(p-hydroxyphenyl) fluorescein (HPF) (Life Technologies, Carlsbad, California, USA) was employed following manufacturer instructions. P. xenovorans strains grown on glucose until early exponential phase were incubated under agitation with the HPF probe in a ratio of 1:1000. After 1 h of incubation, glucose-grown cells were exposed to paraquat or H2O2 (1 mM). The formation of ROS was monitored for 3 h with agitation at 30 °C. Fluorescence was measured at excitation/emission maxima of 490/515 nm with a fluorescence reader (Tecan Trading AG, Männedorf, Switzerland).

Isolation of total RNA and quantitative RT-PCR

Total RNA was isolated from LB400 cells grown on glucose (5 mM) until stationary growth phase and incubated for 1 h in absence and presence of paraquat or H2O2 (1 mM) using a RNeasy mini kit (Qiagen, Hilden, Germany) according to the manufacturer recommendations. DNase treatment was carried out using the TURBO DNA-free Kit (Ambion, LifeTechologies, Carlsbad, CA, USA) to degrade residual DNA. Amplification of the 16S rRNA gene was used as control for DNA contamination using the primers 27f (5′-AGAGTTTGATCMTGGCTCAG-3′) and 1492r (5′-TACGGYTACCTTGTTACGACTT-3′). RNA concentration was quantified using a QubitTM fluorometer (Invitrogen, Carlsbad, California, USA). For RT-qPCR, total RNA (300 ng) was transcribed with First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, Massachusetts, USA) according to the manufacturer recommendations. Quantitative RT-PCR reactions were performed using the Kapa Sybr Fast qPCR Master Mix Universal kit (Hoffmann-La Roche, Basel, Switzerland). Primers employed in the analysis are listed in Additional file 1: Table S1. The gyrB gene was amplified as a reference gene. Quantitative RT-PCR analysis was performed on Mx3000P qPCR system (Stratagene, Agilent Technologies, Santa Clara, California, USA). To quantify gene expression the method of ∆∆CT was employed [44].

Proteome analysis

P. xenovorans LB400 cells were grown until early exponential phase in mineral M9 medium with glucose (5 mM) as sole carbon source. Cells were washed with a NaCl solution (0.9%) and exposed to paraquat and H2O2 (1 mM) for 1 h. Cultures without oxidizing agents were used as control. Three biological replicates were performed under these conditions. For the extraction of total proteins, the cells were suspended in one volume of buffer (10 mM Tris-HCI pH 7.4, 5 mM MgCl2 and 50 pg mL−1 pancreatic RNAase) [70], and sonicated on ice in pulses of 30 s. The cells were centrifuged at 15,022×g at 4 °C, followed by protein precipitation with cold acetone. Samples were incubated at − 20 °C for 1 h, centrifuged at 15,022×g at 4 °C and the supernatant was carefully discarded. The precipitated proteins were stored at − 20 °C. Protein reduction was performed by adding 5 μl of 200 mM dithiothreitol (DTT) to each sample and incubation at room temperature for 30 min. For the carboamidomethylation, 20 μl of 200 mM iodoacetamide were added to each sample and incubated in darkness for 30 min at room temperature. To remove the remnant of iodoacetamide without reacting, 10 μl of 200 mM DTT were added and incubated in darkness for 30 min at room temperature. Protein digestion was performed using lysil-endopeptidase (LEP) and trypsin. For this purpose, the digestion was performed with LEP for 4 h at 30 °C. Digestion with trypsin was performed for 16 h at 37 °C. Finally, the samples were dried by vacuum centrifugation at Speed-Vac concentrator (Thermo Fisher Scientific, Waltham, Massachusetts, USA).

Desalted peptides were loaded onto a reversed phase column packed in house with ReprosilPur C18 Acqua (1.9 µm diameter; Dr. Maisch, Ammerbuch, Germany) with a length of 30 cm and inner diameter of 75 μm using a nano-liquid chromatography system Easy nLC-1000 (Proxeon Biosystems, Odense, Denmark). Peptides were eluted with a gradient of 2–60% of phase B (0.1% v/v formic acid in acetonitrile) and phase A (0.1% formic acid) for 160 min at a flow rate of 250 nL/min. Injections were made by duplicate.

Spectra were acquired in an LTQ Orbitrap XL mass spectrometer (LC–MS/MS) (Thermo Fisher Scientific Waltham, Massachusetts, USA) by data dependent acquisition (DDA), automatically switching between full scan MS (m/z 300–2000) at resolution of 60,000 (m/z 100) and MS/MS with dynamic exclusion of 90 s. The five most intense ions with + 2 and + 3 charges were isolated and fragmented by collision induced dissociation (CID). The mass spectrometer and the gradients in the nLC were controlled by Xcalibur 2.0 software (Thermo Fisher Scientific).

Changes in protein levels were established when changes of ≥ twofold were observed on treated cells compared with untreated cells. The signals corresponding to each peptide obtained from the analysis by mass spectrometry were quantified using the software Mascot Distiller (MatrixScience, Boston, MA, USA) [59]. Subsequently, proteins were identified using P. xenovorans LB400 protein databases. Heatmap visualization of proteins classified according to their function was performed using the pheatmap V 1.0.8 R package [37].

Statistical analyses

Statistical analyses were performed using one-way ANOVA and LSD Fisher test to assess differences in mean values from each experiment. Differences are significant at p < 0.05.


P. xenovorans LB400 possess genes for a broad antioxidant response

In order to study the oxidative stress response in P. xenovorans strain LB400, a search of ortholog oxidative stress genes belonging to the OxyR regulon from E. coli that is activated by H2O2 was performed. Several oxidative stress genes were identified in the genome of strain LB400, including genes encoding 11 alkyl hydroperoxide reductase system subunits (AhpC, AhpD and AhpF), 4 catalases, 3 superoxide dismutases, and a high number of peroxidases and peroxiredoxins (Additional file 1: Table S2). The highest identities (> 80%) were observed with proteins from bacteria of the order Burkholderiales. Two subunits of the enzyme alkyl hydroperoxide reductase system, AhpC1 and AhpD1, are encoded by the gene cluster ahpC1D1 (BxeA2309 and BxeA2310) on the LB400 major chromosome (C1). The minor chromosome (C2) possesses 7 copies of the gene that encodes the enzyme alkyl hydroperoxide reductase subunit D (AhpD), and the gene cluster ahpC2F (BxeB1205 and BxeB1206). On the C1 and C2 chromosomes of the strain LB400, the katA, katG, katN and katE genes encoding for catalases (one Mn-catalase and three heme-catalases) were identified. The redundancy in the antioxidant systems in strain LB400 is reflected by the high number of copies of oxidative stress response genes. An analysis of the genomic context of the oxyR transcriptional regulator gene indicated that it is clustered with other oxidative stress genes, such as the katA gene. Upstream of the oxyR gene, the recG gene (BxeA3988) encoding an ATP-dependent DNA helicase RecG was identified (Fig. 1). Other genes of the LB400 OxyR regulon are the oxyR gene encoding the OxyR transcriptional regulator (BxeA3987) and the gstA1 gene that encodes a glutathione S-transferase (BxeA0624).

Fig. 1
figure 1

Genomic context of genes involved in oxidative stress response in P. xenovorans LB400. C1 and C2 indicate the major and the minor chromosome, respectively

A search of genes involved in oxidative stress response in P. xenovorans LB400, whose orthologous genes are associated to the SoxRS regulon in E. coli was performed. The SoxRS regulon in E. coli activates the expression of oxidative stress genes, such as sodA (superoxide dismutase), fumC (oxidative stress resistant fumarase), acnA (ROS resistant aconitate hydratase), fpr (ferredoxin NADP reductase) and fldA/fldB (flavodoxins A and B) [28]. LB400 genome possesses genes encoding the enzyme aconitate hydratase (acnA1, acnA2), Mn-Fe superoxide dismutase and Cu–Zn superoxide dismutase (sodB1, sodB2, sodC), and a flavodoxin/ferredoxin NADP reductase (fpr). The oxidative stress genes identified in the LB400 genome possess high identity with the corresponding proteins of bacteria belonging to the class Betaproteobacteria. LB400 strain possesses the fumC gene encoding for fumarate hydratase and the trxB1 and trxB2 genes (BxeA3442 and BxeA3962) encoding thioredoxin reductases. Two LB400 genes encode long-chain and short-chain flavodoxins FldX1 and FldX2, which protect cells during oxidative stress [63]. In addition, two homologs of the organic hydroperoxide resistance protein Ohr (BxeB2195 and BxeB2843), regulated in E. coli by the organic hydroperoxide resistance transcriptional regulator OhrR, was observed in the LB400 genome. The regulator OhrR was also identified in LB400 strain (BxeB2842). The MerR-type transcriptional regulator SoxR was searched in the LB400 genome. In the LB400 genome, 15 genes encoding MerR-type family regulators were identified. The BxeC1217 (soxR1) gene possesses the highest identity with the soxR genes from P. aeruginosa (31%) and E. coli (33%). The soxS gene is not present in the LB400 genome, which is expected for non-enteric Proteobacteria. This is in accordance with the observation that SoxR directly activates antioxidant genes in non-enteric bacteria [23, 30]. Other interesting genes involved in stress response are described in Additional file 1: Table S2. The high number of genes identified in the LB400 genome encoding antioxidant enzymes, in some cases with several copies, suggests a broad enzyme repertoire for the cellular response to oxidative stress.

Phylogenetic analysis of transcriptional regulators

To investigate the function of the OxyR and SoxR transcriptional regulators in P. xenovorans LB400, their phylogenetic relationships were studied. For OxyR analysis, 17 amino acid sequences with experimental evidence in bacteria from Alphaproteobacteria, Betaproteobacteria, Gammaproteobacteria and Actinobacteria classes were used. A Bayesian Inference tree evidenced four main clades within OxyR orthologs, mainly clustered by taxonomic relatedness of the sequences (Fig. 2A). The first clade (group I, Fig. 2A) consists of Betaproteobacteria and Gammaproteobacteria. P. xenovorans LB400 OxyR (BxeB3987) clustered with the OxyR regulators of other Burkholderiales strains (Ralstonia solanacearum K60 and Alcaligenes aquatilis QD168) (group IA, Fig. 2A). The second cluster in this clade (group IB, Fig. 2A) is formed by bacteria from Pseudomonadales order, Gammaproteobacteria class (P. aeruginosa PA01 and Pseudomonas chlororaphis GP72). In the same clade, OxyR of Neisseria gonorrhoeae 1291 and Neisseria meningitidis MC58 (Neisseriales order, Betaproteobacteria) grouped together (group IC, Fig. 2A). The second clade (group II, Fig. 2A) is composed by Gammaproteobacteria Haemophilus influenzae KW20 (Pasteurellales order) and Enterobacterales order strains (E. coli K12, Pectobacterium carotovorum subsp. carotovorum AE1202, Dickeya chrysanthemi EC16 and Pantoea stewartii subsp. stewartii DC283). The third clade consists only of OxyR from Actinobacteria, clustered by Mycobacteriales order (Mycobacterium avium ATCC 19075, Mycobacterium marinum ATCC 15069 and Corynebacterium diphtheriae C7) and Streptomyces viridosporus T7A (Streptomycetales order). OxyR of Brucella abortus from Alphaproteobacteria class, and OxyR of Acinetobacter oleivorans DR1 (Gammaproteobacteria) formed the fourth clade (group IV, Fig. 2A).

Fig. 2
figure 2

Evolutionary relationships of OxyR and SoxR of P. xenovorans LB400 and other bacteria. Mid-rooted Bayesian Inference trees calculated by MrBayes of the OxyR and SoxR regulators with experimental evidence. P. xenovorans LB400 OxyR (PxOxyR) and SoxR (PxSoxR) are shown in bold letters and marked with a dark circle. Node values represent the bootstrapping value (%) of the analysis. Each group determined was identified at the class rank level. If the amino acid sequences within each group correspond to the same family rank, the taxon name is included in parentheses. A Phylogenetic tree of the H2O2-sensing transcriptional regulator OxyR homologs. Four major groups were identified, predominantly clustered by taxonomic relatedness. PxOxyR is encompassed within the Group IA, represented by Burkholderiales. B Phylogenetic tree of the redox-sensitive transcriptional regulator SoxR homologs. Four major groups are identified, including three singleton sequences, not observing a clear taxonomic distribution within the groups. PxSoxR is a singleton closer to the SoxR regulator of Chromobacterium violaceum DSM 30191 than to other groups

For SoxR analysis,16 amino acid sequences with experimental evidence in bacteria from Betaproteobacteria, Gammaproteobacteria and Actinobacteria classes were used. In contrast to OxyR transcriptional regulators that a similar role in diverse bacteria, the SoxR regulators from diverse bacteria have been involved in different functions such as the production and export of redox cycling compounds, the oxidative stress protection, and the development [22, 30, 80]. Therefore, a higher diversity of SoxR regulators is expected due to their association to different functions. The Bayesian inference phylogeny revealed four clades and three singletons, observing an overall robust taxonomic profiling (Fig. 2B). Unlike LB400 OxyR, which clustered with OxyR from other Burkholderiales order, SoxR of strain LB400 (BxeC1217) formed a singleton suggesting a new clade of SoxR regulators. The SoxR of Pseudomonas putida KT2440 clustered within the first clade with SoxR of A. oleivorans DR1 (Gammaproteobacteria,Pseudomonadales order), and SoxR1 and SoxR2 of A. aquatilis QD168 (Betaproteobacteria; Burkholderiales order) (group I, Fig. 2B). In contrast, P. aeruginosa PAO1 SoxR (Pseudomonadales order; Betaproteobacteria class) formed a singleton. The second clade is composed of SoxR of the Enterobacterales order (E. coli K12, Salmonella typhi CT18, Salmonella typhimurium LT2, Shigella flexneri 301, P. stewartii subsp. stewartii DC283 and Klebsiella pneumoniae KPBj5 E) (group II, Fig. 2B). SoxR of the Gammaproteobacteria Stenotrophomonas maltophilia K279a (Lysobacterales order) and Xanthomonas campestris pv. campestris P25 (Xanthomonodales order) form the third clade (group III, Fig. 2B). The fourth clade consisted of SoxR of Actinobacteria from the Streptomycetales order, Streptomyces coelicolor M145 and Streptomyces avermitilis MA-4680 (group IV, Fig. 2B). SoxR of Chromobacterium violaceum DSM_30191 (Neisseriales order, Betaproteobacteria) formed a singleton.

P. xenovorans LB400 SoxR has a protective role in P. aeruginosa upon exposure to oxidizing agents

In order to determine the functionality of the P. xenovorans LB400 soxR gene (BxeC1217), BxeC1217 was complemented in the null mutant strain P. aeruginosasoxR. Susceptibility assays were performed using the disk diffusion method in the presence of different concentrations of paraquat (10 and 20 mM). Paraquat is a redox-cycling compound that is widely used as herbicide, which is transported inside the bacterial cell and constitutes a continuous source of superoxide radical. Inside the cell, paraquat is reduced by electron donors such as NAD(P)H and is oxidized via the transfer of electrons to the electron acceptor dioxygen, producing superoxide. Redox-cycling agents such as paraquat are able to activate the SoxR transcriptional regulator [23].

The complemented strains P. aeruginosasoxR::BxeC1217 and P. aeruginosasoxR::PA2273 were less susceptible to paraquat than the mutant strain P. aeruginosasoxR, with a smaller zone of inhibition (Fig. 3A). These results suggest that LB400 soxR gene (BxeC1217) encodes a functional transcriptional regulator SoxR that has a protective effect during oxidative stress induced by paraquat.

Fig. 3
figure 3

The soxR gene of P. xenovorans LB400 complements faulty oxidative response in P. aeruginosa ∆soxR. A Susceptibility assays of P. aeruginosa strains PA14 wt, and mutant strains (∆soxR, ∆soxR::BxeC1217 and ∆soxR::PA2273) exposed to paraquat (PQ). Paper disks soaked with paraquat (10 or 20 mM) were placed on P. aeruginosa bacterial lawns. Growth inhibition zones around the disks were recorded after 24 h at 37 °C. The dotted line at 6 mm corresponds to the disk diameter and indicates no detectable inhibition. Inhibition zone values were calculated as the mean ± SD of three independent experiments. Letters under the error bars indicate significant differences between strains in each treatment. B Growth of P. aeruginosa wt and mutants strains exposed to phenazine methosulfate (PMS) (600 µM). Control cells were grown in absence of PMS. Images are representative of three independent experiments

To further characterize the role of the P. xenovorans LB400 soxR gene, the formation of complemented P. aeruginosa macrocolonies in presence of the toxic antibiotic PMS was evaluated (Fig. 3B). After 46 h of incubation growth of the null mutant strain ∆soxR on plates supplemented with PMS (600 µM) was scarcely observed, while P. aeruginosa strains ∆soxR::BxeC1217, ∆soxR::PA2273 and the wild type strain PA14 grew forming macrocolonies, indicating that BxeC1217 can restore the growth-associated phenotype when exposed to PMS. A yellow-green pigmentation started to appear around 118 h in ∆soxR::BxeC1217, which was not observed in PA14 or ∆soxR::PA2273 strains (Fig. 3B). The pigmentation suggests that SoxR of LB400 strain regulates also the synthesis in P. aeruginosa of the yellow phenazine-1-carboxylic acid and the blue phenazine pyocyanin [21].

Effect of oxidizing compounds on P. xenovorans LB400 growth

To study the physiological response of the strain LB400 to oxidative stress, growth studies were performed using cells exposed to different concentrations of paraquat and H2O2. Strain LB400 cells grown on glucose as sole carbon and energy source until early exponential phase were exposed to oxidizing agents in the concentration range from 1 to 16 mM. Growth of strain LB400 was not affected after exposure to paraquat and H2O2 (1 mM). These results revealed that the oxidizing agents at the concentration 1 mM are not toxic for strain LB400. In contrast, a significant decrease in growth was observed at concentrations ≥ 4 mM of paraquat and H2O2 (Fig. 4A). Interestingly, an inhibitory effect on growth was observed 30 min upon H2O2 (4–16 mM) exposure, whereas a negative effect on growth with cell reduction was observed after 2 h exposure to paraquat (4–16 mM) (Fig. 4). Hydrogen peroxide-treated cells (4–16 mM H2O2) recovered their growth after 120 min. In contrast, paraquat-exposed cells (4–16 mM paraquat) after 120 min reduced progressively the cell turbidity, indicating cell death.

Fig. 4
figure 4

Effects of paraquat and H2O2 on P. xenovorans LB400 growth and ROS formation. Strain LB400 grown until exponential growth phase on glucose was exposed to paraquat (PQ) (A) or H2O2 (B). Growth was monitored by measuring turbidity at 600 nm. C Paper disks soaked with paraquat or H2O2 solutions (1–20 mM) were placed on P. xenovorans LB400 bacterial lawns. Growth inhibition zones around the disks were recorded after 24 h incubation at 30 °C. The dotted line at 6 mm corresponds to the disk diameter and indicates no inhibition. Letters under the error bars indicate significant differences between each concentration of each oxidizing agent. D Exponential glucose-grown cells were exposed to paraquat or H2O2. After incubation with HPF, fluorescence was monitored at 515 nm (emission) and 490 nm (excitation). Cells incubated in absence of oxidizing agents were used as control. Letters under the error bars indicate significant differences between treatments in each time. All values were calculated as the mean ± SD of three independent experiments

LB400 susceptibility to paraquat and H2O2 (1, 5, 10 and 20 mM) was assessed during growth in solid medium, using diffusion disks (Fig. 4C). In the case of H2O2, raising concentrations of the oxidizing agent showed no differences in the inhibition zone (~ 8 mm). An increase in paraquat concentration showed an increase in the inhibition zone within 5 and 20 mM, indicating that strain LB400 is more susceptible to paraquat than to H2O2.

P. xenovorans LB400 is susceptible to paraquat, which increases intracellular ROS levels

The formation of endogenous ROS in strain LB400 during exposure to paraquat and H2O2 were studied using a specific probe for the hydroxyl radical. LB400 cells grown to exponential phase on glucose were incubated with a ROS indicator for 1 h. Subsequently, the cells were exposed to paraquat and H2O2. As shown in Fig. 4D, ROS accumulation was higher in cells exposed to H2O2 compared to non-exposed cells at short incubation times (0.5 h). An increase in ROS accumulation in paraquat-exposed cells was observed after 1 h incubation, which is later than ROS accumulation in hydrogen peroxide-exposed cells (Fig. 4D).

Paraquat and H 2 O 2 activate oxidative stress response genes

Expression of selected genes by RT-qPCR under mild oxidative stress conditions was studied. In order to study the role of selected antioxidant genes, their expression levels were measured in LB400 cells incubated in absence and presence of sublethal concentrations of paraquat or H2O2 (1 mM). Transcriptional analysis by RT-qPCR of the oxyR (BxeA3987), katE (BxeB1215), fumC (BxeA1038), acnA2 (BxeB2903), ahpC1 (BxeA2309), sodB1 (BxeA0769) genes was performed. Transcriptional expression of additional genes encoding two thioredoxin reductases (trxB1 and trxB2; BxeA3442 and BxeA3962, respectively), a high potential Fe-S protein (hpf; BxeA4333), the OhrB protein (ohrB; BxeB2843), a flavodoxin/ferredoxin NADP reductase (fpr; BxeA4345) and a glutathione S-transferase enzyme (gstA1; BxeA0624) was also analyzed. In presence of paraquat, the expression of the genes that encode OxyR, AhpC1 and OhrB increased > sixfold compared to non-treated cells (Fig. 5). The fumC, trxB1 and trxB2 genes were upregulated in presence of paraquat ≥ fourfold, whereas the fumC, hpf, and gstA1 genes increased their expression ≥ twofold. Upon exposure to H2O2, genes that encode the OxyR transcriptional regulator and the OhrB protein increased > fourfold, while AhpC1-coding gene increased > sixfold. Interestingly, expression of the gene encoding a high potential Fe-S protein (Hpf) increased ≥ 24-fold in presence of H2O2 compared to the control condition (Fig. 5). In addition, in presence of hydrogen peroxide the transcription levels of the acnA2, katE genes, and genes encoding a ferredoxin (fpr) and thioredoxin reductases (trxB1 and trxB2) increased ≥ twofold. The increase in the expression of antioxidant genes encoding AhpCF and the OxyR transcriptional regulator in the presence of both oxidizing agents, indicate under these oxidizing conditions, a general type oxidative stress response against paraquat and H2O2.

Fig. 5
figure 5

Expression of oxidative stress genes in P. xenovorans LB400 upon exposure to paraquat or H2O2. Expression levels of oxidative stress genes were measured by RT-qPCR after strain LB400 exposure to paraquat (PQ) or H2O2 (1 mM) for 1 h. The gyrB gene was used as the reference gene. Fold-change levels were calculated as the mean ± SD of three independent experiments. The dotted lines indicate the level of expression from which a significant increase (2) or decrease (-2) of expression is observed

Paraquat and H 2 O 2 induce a stress-associated response

In order to further study the physiological response of P. xenovorans LB400 to paraquat and H2O2, a proteome analysis was performed. LB400 cells grown on glucose until exponential phase were exposed to paraquat and H2O2 at a sublethal concentration (1 mM) for 1 h. Tables S3 and S4 indicate the significantly induced and repressed (twofold) proteins during exposure to these oxidizing agents.

During exposure to paraquat (Fig. 6A; Additional file 1: Table S3), the levels of stress-related proteins were increased. Interestingly, the alkyl hydroperoxide reductase AhpC2 and AhpF subunits, and a ferritin protein of the DPS family (DpsA) were induced. The BxeA3984 gene encoding the DpsA protein is present in the oxyR gene neighborhood of strain LB400. The organic hydroperoxide resistance protein (OhrB, BxeB2843) was also induced (Fig. 6A; Additional file 1: Table S3). The response of strain LB400 to paraquat included the increase of detoxifying and repair enzymes, along with general stress proteins (Fig. 6A; Additional file 1: Table S3). Induction of the universal stress protein UspA (BxeB0607) was observed during paraquat exposure, as well as three RNA chaperones from the cold-shock DNA-binding protein family (BxeA0430, BxeA0798, BxeB2951). The induction of these chaperones suggests a condition of cellular damage at protein and DNA levels. Interestingly, phasin protein (BxeA1544), a surface protein related to polyhydroxyalkanoates (PHAs) granules, was induced in the presence of paraquat. Paraquat also induced two oxidoreductases (BxeA2478 and BxeA4053), a methylenetetrahydrofolate reductase (MetF) and a transcriptional regulator of the LysR family (BxeA2466), in whose genomic context the genes encoding the CysP, CysT, CysW and CysA subunits of the sulfate/thiosulfate ABC transporter (BxeA2467-BxeA2470) are present. The repression of the stress chaperones GroEL, GroES, HtpG and Hsp20 was observed during exposure to paraquat (Fig. 6B; Additional file 1: Table S4). Moreover, the expression level of the protein UvrA of the UvrABC system (BxeA4108) decreased. An RNA helicase ATP-dependent (BxeA1652) was also downregulated during paraquat exposure (Fig. 6B; Additional file 1: Table S4).

Fig. 6
figure 6

Differentially expressed proteins of P. xenovorans LB400 after exposure to paraquat or H2O2. Heatmap showing upregulated (A) or downregulated (B) protein levels of P. xenovorans LB400 after exposure to paraquat (PQ) or H2O2 (1 mM) for 1 h compared to untreated cells. Changes in protein levels were established when changes of ≥ twofold were observed on treated cells versus control cells. Red represents higher expression, whereas blue represents lower expression. Three biological replicates were performed. AhpC, alkyl hydroperoxide reductase subunit C2; organic hydroperoxide resistance protein OhrB; Dps, ferritin DPS family DNA-binding protein DpsA. OxRd oxidoreductase, Tas Tas oxidoreductase, TR transcriptional regulator of LysR family, OEP outer membrane efflux protein related to copper resistance, TP transport protein, HDH homoserine dehydrogenase, GK glycerate kinase, SRd sulphite reductase (NADPH) beta-subunit, PBP penicillin binding-protein. Additional file 1: Tables S3 and S4 provide more details of the proteins

Upon H2O2 exposure, the organic hydroperoxide resistance protein Ohr (BxeB2843) was strongly induced (Fig. 6A; Additional file 1: Table S3). Two transporters were induced, an outer membrane efflux pump related to copper resistance (BxeB2295) and a transport protein (BxeB1679). Several proteins of the general metabolism were downregulated, including the sulfite reductase beta subunit (sulfur metabolism), a NADH quinone oxidoreductase (oxidative phosphorylation), and a glycerate kinase that participates in serine/glycine/threonine metabolism, glycolipids metabolism and glyoxylate carboxylate metabolism. In addition, the induction and repression of proteins involved in the RNA metabolism, as well as protein synthesis was observed (Fig. 6; Additional file 1: Tables S3 and S4). The induction of the antitermination factor NusG and homoserine dehydrogenase was observed in H2O2-exposed LB400 cells. The expression levels of the ribosomal proteins S6 and L20 were upregulated, whereas the L7/L12 protein of the 30S and 50S ribosomal subunits were downregulated (Fig. 6; Additional file 1: Tables S3 and S4). The general stress chaperones GroEL, GroES and HtpG as well as the stress chaperones ClpB and HSP20 were repressed.

OxyR has a protective role in strain LB400 against oxidative stress

In order to study the role of OxyR in strain LB400 during exposure to oxidizing agents, the LB400 oxyR gene was overexpressed in P. xenovorans (strain pIZoxyR). Susceptibility assays to oxidizing compounds and ROS formation was determined in cells that overexpressed the oxyR gene.

Susceptibility assays were performed using different concentrations of paraquat and H2O2 (10 and 20 mM). As shown in Fig. 7A, strain pIZoxyR was less sensitive to paraquat, showing a slighter zone of inhibition than the control P. xenovorans strain (pIZ1016). However, no significant differences upon exposure to H2O2 were observed between inhibition zones of the pIZoxyR and the reference strain (data not shown).

Fig. 7
figure 7

Susceptibility and ROS formation in P. xenovorans pIZoxyR exposed to paraquat. A Susceptibility assays of P. xenovorans strains exposed to paraquat (PQ). Paper disks soaked with paraquat (1–20 mM) were placed on P. xenovorans strain bacterial lawns. Zones of growth inhibition around the disks were recorded after 24 h incubation at 30 °C. The diameter of 6 mm corresponds to the disk diameter and indicates no inhibition. Inhibition zone values were calculated as the mean ± SD of three independent experiments. Letters under the error bars indicate significant differences between strains in each treatment. B ROS quantification assay. P. xenovorans pIZoxyR was grown on glucose (5 mM) until exponential growth phase and exposed to paraquat (20 mM). ROS formation of non-treated cells was also monitored. Fluorescence was monitored after incubation with HPF at 490 nm (excitation) and 515 nm (emission). Fluorescence values were calculated as the mean ± SD of three independent experiments. Letters under the error bars indicate significant differences between strains in each time

Subsequently, ROS formation was measured in oxyR-overexpressing cells exposed to paraquat and H2O2 (20 mM). A decrease in ROS accumulation was observed in P. xenovorans strain pIZoxyR during exposure to paraquat compared to reference strain (Fig. 7B). However, no significant differences in ROS formation were observed between pIZoxyR and pIZ1016 cells exposed to H2O2 (data not shown). These results, as well as those observed on susceptibility assays, suggest that OxyR contributes significantly to the oxidative stress response of P. xenovorans upon exposure to paraquat.


Oxidative stress in P. xenovorans LB400. In this study, we characterized the oxidative stress response of the model aromatic-degrader P. xenovorans LB400 to the redox-cycling aromatic compound paraquat and the ROS hydrogen peroxide. The results revealed by genomic, transcriptomic, and proteomic analyses indicate that strain LB400 possess a robust antioxidant enzymatic repertoire. Oxidative stress response in P. xenovorans LB400 and other bacteria in presence of aromatic compounds has been previously reported. The aerobic degradation of aromatic compounds is catalyzed by mono-and dioxygenases that uses dioxygen as substrate may produce ROS as side products [33, 61]. During p-cymene degradation, LB400 proteome showed an increase of the alkyl hydroperoxide reductase AhpC, the organic hydroperoxide resistance protein Ohr, and the molecular chaperones DnaK, GroEL and ClpB [3]. p-Cymene also induced diverse proteins of the energy metabolism, such as aconitase AcnA, succinyl-CoA synthetase, ATP synthase, enolase and pyruvate kinase, which are required for the high-energy demand during the adaptation of glucose-grown LB400 cells to aromatic compounds. (Chloro)biphenyls degradation by P. xenovorans LB400 involves oxidative stress, inducing the synthesis of an alkyl hydroperoxide reductase AhpC and molecular chaperones DnaK, GroEL, HtpG and elongation factor G [2, 61]. Chlorobenzoates also generate a general stress response in LB400 cells, inducing the expression of the molecular chaperones DnaK and HtpG [48]. Burkholderia sp. BNS showed during phenol catabolism an increase in superoxide dismutase, DNA repair enzymes RecN and RadA, and chaperones Hsp15 and Hsp12 [47]. Fluoranthene catabolism induces expression of catalase and superoxide dismutase in Mycobacterium sp. JS14 [42]. Bacillus subtilis 168 showed a oxidative stress response to catechol, inducing the alkyl hydroperoxide reductase genes [74]. The addition of antioxidants and the reconfiguration of metabolic fluxes are useful strategies to counteracting oxidative stress of the model aromatic-degraders P. xenovorans, P. putida and other bacterial strains, which may increase the degradation and bioremediation of persistent organic pollutants [26, 54, 61].

High redundancy of antioxidant genes. The broad genetic repertoire of antioxidant genes identified in P. xenovorans LB400 (Fig. 1; Additional file 1: Table S2) correlates with its high resistance to hydrogen peroxide and paraquat. LB400 has two clusters encoding two complete alkyl hydroperoxide reductase systems (ahpC1D1 and ahpC2F), and 7 additional ahpD gene copies in its genome, which are unusual high ahp gene copies in bacteria. In aerobic bacteria, the alkyl hydroperoxide reductase detoxifies organic hydroperoxides and H2O2 during oxidative stress. In E. coli, B. subtilis, Salmonella enterica sv. typhimurium and Xanthomonas spp., this enzyme consists of two subunits, a minor subunit with peroxidase activity (AhpC, 22 kDa), and a major subunit with disulfide reductase activity (AhpF, 57 kDa) [62]. AhpC has been reported as the main antioxidant mechanism against endogenous H2O2 production in E. coli [68]. In S. coelicolor and Mycobacterium species, the ahpF gene is absent, whereas downstream of the ahpC gene is located the ahpD gene, which encodes a thioredoxin-like protein (19 kDa) involved as electron donor in the disulfide reduction of AhpC. Notably, P. xenovorans LB400 encodes both AhpCF and AhpCD systems, which has been rarely observed in bacteria. Bioinformatic analyses indicate that Paraburkholderia hospital BS001 has genes encoding five AhpC, one AhpD and one AhpF, whereas Paraburkholderia terrae possess four AhpC, one AhpD and one AhpF. In Burkholderia thailandensis, the AphCD system has been recently characterized [82]. E. coli K12, P. putida KT2440 and B. subtilis 168 possess only one or two copies of AhpC and one AhpF. S. coelicolor ATCC BAA-471, Bradyrhizobium diazoefficiens JCM 10833 and Mycobacterium tuberculosis ATCC 25618 possess one or two copies of AhpC and one or two copies of AhpD.

Four genes encoding catalases (three heme-catalases and one Mn-catalase) are present in the LB400 genome. LB400 KatE possesses a sequence identity of 27.6, 15.3 and 8.8% with LB400 KatA, KatG and KatN, respectively. Uniprot database indicate that the toluene degrader P. putida KT2440 possesses three catalases (KatA, KatE and KatG), E. coli strain K12 has two catalases (KatE and KatG), and B. subtilis 168 possess four catalases (KatA, KatE, KatG and KatX). Three superoxide dismutase-encoding genes were identified in strain LB400. Superoxide dismutase detoxifies superoxide radical, transforming it in H2O2 and O2 [29]. Strain LB400 showed high gene redundancy for the oxidative stress response, redundancy that has been previously described in this bacterium for genes involved in aromatic degradation [14]. A lower redundancy of antioxidant genes (4 ahp, 3 kat, 2 sod) has been reported in P. aeruginosa [55]. The high number of genes encoding antioxidant enzymes in P. xenovorans LB400 may contribute to its fitness and tolerance to oxidative stress in environments where bacteria are subjected to a range of adverse environmental conditions.

The transcriptional regulator oxyR gene of P. xenovorans LB400 grouped together with other antioxidant genes in the recG-oxyR-katA-dpsA gene cluster on the major chromosome. A similar gene organization has been observed in Betaproteobacteria and Gammaproteobacteria. This is in accordance with the phylogenetic analysis of OxyR, where LB400 OxyR clustered with OxyR of members from the Burkholderiales order (Betaproteobacteria class). The gene context of the LB400 oxyR gene suggests that the ahpC1D1, katA and dpsA genes are activated by OxyR in response to oxidizing compounds, which has been reported in E. coli [28]. Interestingly, our study indicates that the LB400 DpsA protein and the oxyR gene were induced upon exposure to superoxide-producing compound paraquat. In B. pseudomallei, the OxyR transcriptional regulator activates the katG-dpsA operon under oxidative stress [32]. The oxyR-recG locus is essential in P. aeruginosa for the defense to oxidative stress generated by H2O2 and paraquat [55].

Functionality of OxyR and SoxR regulators. In this study, the functionality of OxyR and SoxR regulators in P. xenovorans LB400 was demostrated, which showed a protective role against oxidizing agents. The role of the soxR gene of P. xenovorans LB400 was determined by a complementation assay. In our study, the mutant strain P. aeruginosa ΔsoxR complemented with the LB400 soxR gene (BxeC1217) was less sensitive to paraquat and phenazine methosulphate, two superoxide-producing compounds. This is particularly interesting considering that SoxR phylogenetic analysis revealed that on one side SoxR of LB400 strain is not closely related to other well characterized SoxR, forming a distinct clade, and on the other side, SoxR showed only 33% identity with its closed bacterial SoxR. The complemented strain ∆soxR::BxeC1217 also changed P. aeruginosa pigmentation. These results indicate that the regulator SoxR of P. xenovorans LB400 is functional, playing an important role to protect bacteria from superoxide-producing agents (paraquat and PMS), and may activate a specific regulon of not enteric bacteria involved in phenazine production in P. aeruginosa. The SoxR regulator is a redox-sensitive transcriptional activator during oxidative stress produced by endogenous or environmental factors. SoxR of P. aeruginosa is activated by compounds that produce superoxide, such as paraquat, or also by its own synthesized phenazines [22, 30, 36]. P. aeruginosa SoxR activates the PumA monooxygenase and the MexGHI-OpmD efflux pump, which are involved in the modulation of cellular intracellular redox state, and the resistance to exogenous and endogenous phenazines [21, 22, 65, 72]. PumA converts the exogenous phenazine PMS into a highly oxidized compound, a mechanism also proposed for intracellular phenazines produced by P. aeruginosa PA14 to control the internal phenazine pool and the redox intracellular state [72]. SoxR from P. putida, A. oleivorans, S. coelicolor, and E. coli are differentially activated by redox-active compounds [34, 42], which could also be the case of LB400 SoxR, inducing an divergent activation of the SoxR regulon in P. aeruginosa that may explain ∆soxR::BxeC1217 yellow-green pigmentation, probably associated with the synthesis of the yellow phenazine-1-carboxylic acid and the blue phenazine pyocyanin [21].

Notably, we also showed that oxyR (BxeB3987) overexpression in P. xenovorans conferred higher resistance to the oxidizing agent paraquat (Fig. 7). Paraquat is a continuous source of superoxide radical in the cytoplasm. Paraquat and H2O2 may activate the OxyR transcriptional regulator in non-enteric bacteria. Interestingly, we observed that paraquat induces stronger the expression of LB400 OxyR than hydrogen peroxide. OxyR activates the H2O2 scavenging enzymes catalase and alkyl hydroperoxide reductase, which may explain the protective effects observed during paraquat exposure in the P. xenovorans strain overexpressing the transcriptional regulator OxyR.

Response to paraquat and H2O2. The physiological and molecular response of strain LB400 upon exposure to paraquat and H2O2 is summarized in Fig. 8. H2O2 showed after shorter incubation growth inhibition and increased cytoplasmic ROS accumulation probably due to its oxidizing nature, whereas paraquat showed negative effects on growth and increased cytoplasmic ROS formation only after a longer exposure likely because it requires first to generate radicals. Paraquat showed a higher reduction of cell density than H2O2, probably due to its continuous production of superoxide in the cytoplasm and the membrane diffusion limitation of H2O2. In E. coli the membrane permeability coefficient for H2O2 is 1.6 × 10−3 cm/s, indicating that H2O2 diffusion inside the cell is limited [69]. This could explain that different concentration of H2O2 showed no significant differences on growth in liquid medium (Fig. 4B) and in solid medium (Fig. 4C). During the period between 120 and 300 min, the growth rates of cells exposed to different H2O2-concentration and control condition were similar (~ 10–4 h−1). In presence of paraquat, higher concentration of the oxidizing agent increased the negative effect on bacterial density, probably to the higher production of superoxide radical inside the cell. Proteomic and transcriptional analyses also indicated that paraquat induced a stronger oxidative stress response than hydrogen peroxide. During exposure to paraquat, the alkyl hydroperoxide reductase system AhpC2F was induced. In addition, the upregulation of the oxyR and ahpC1 genes of strain LB400 was observed in the presence of paraquat and H2O2, whereas an increase in the katE gene expression was observed in cells incubated with H2O2. The increase of the ahpC1 gene expression during exposure to paraquat suggests that the superoxide radical produced by paraquat is reduced by superoxide dismutase into H2O2. These results suggest that P. xenovorans displays a broad antioxidant response that includes AhpCF, a catalase and the DpsA protein to protect from the oxidative stress generated by both paraquat and H2O2, which probably is regulated by OxyR. The hydrocarbon-degrading bacterium A. aquatilis QD168 also showed an upregulation of its two ahpC gene copies upon paraquat exposure [25]. In bacteria, detoxification of H2O2 is carried out mainly by the action of alkyl hydroperoxide reductases and catalases [29, 58]. In E. coli, the AhpCF enzyme is activated by the regulator OxyR [28]. This indicates that the presence of H2O2 and paraquat may be triggering the expression of LB400 oxidative stress genes that are probably regulated by the transcriptional regulator OxyR. The oxyR gene also increased its expression in strain LB400 in presence of paraquat (> 14-fold) and H2O2 (fourfold). In E. coli, the OxyR regulator activates a specific antioxidant response against H2O2 with transcription of several oxidative stress genes [28]. The OxyR regulator in E. coli is sensitive to H2O2, whereas in P. aeruginosa other non-enteric models the OxyR regulator is activated by different oxidizing agents such as paraquat and H2O2 [23]. The results of our study suggest that the non-enteric bacterial strain LB400 unfolds a broad stress response against superoxide-producing compound paraquat and H2O2 mediated by the OxyR regulator.

Fig. 8
figure 8

Molecular oxidative stress response of P. xenovorans LB400 during exposure to paraquat and H2O2. Green arrows indicate upregulation of genes and/or proteins, while red arrows indicate downregulation. Dashed arrows indicate the transport of O2, H2O2 and paraquat through the cell membranes. Dotted arrow indicates Fenton reaction, a process during which metal ions (mainly Fe2+) react with H2O2 forming hydroxyl radical. Glucose is converted into gluconate 6-phosphate (gluconate-6P) by glucokinase and the Pentose Phosphate (PP) pathway. Blue arrow indicates the transformation of the gluconate-6P to glyceraldehyde-3-phosphate (G3P) and pyruvate via the Entner–Doudoroff pathway. Purple arrow indicates the conversion of G3P into pyruvate via the lower Embden–Meyerhof–Parnas (EMP) pathway. OM outer membrane, PS periplasmic space, IM inner membrane, FumC/fumC fumarate hydratase C, AcnB aconitate hydratase, hpf high potential Fe-S protein, trxB1 thioredoxin 1, trxB2 thioredoxin 2, fpr flavodoxin/ferredoxin NADP oxidoreductase, sodB1 superoxide dismutase, DnaK GroEL and GroES, molecular chaperones, Hsp20 heat shock protein, UspA universal stress protein, ahpC1D1 and ahpC2F alkyl hydroperoxide reductases, katE catalase, DpsA ferritin protein of the DPS family, Gst/gstA1 glutathione S-transferase, Ohr organic hydroperoxide resistance protein, OhrR organic hydroperoxide resistance transcriptional regulator, MetF 5,10-methylenetetrahydrofolate reductase

Other antioxidant genes that may be part of the OxyR regulon in LB400 were upregulated upon paraquat and H2O2 exposure. The hpf gene that encodes for a high potential Fe-S protein, increased > 24-fold in presence of these compounds, whereas the trxB1 and trxB2 genes increased twofold. The upregulation of thioredoxin reductases and the Fe-S protein of high potential Hpf in presence of H2O2 and paraquat, are probably contributing to the LB400 redox state of cells providing reducing power in these oxidizing conditions. Fe-S proteins of high potential are electron-donors that contain a cytochrome tetraheme in their active center [78]. Thioredoxin reductases are widely distributed homodimeric flavoproteins. Bacterial thioredoxin reductases are low molecular weight proteins and their function may vary between organisms. In E. coli, thioredoxin reductases transfer reducing equivalents to a wide range of enzymes, which is critical for DNA synthesis and defense against oxidative stress [45]. In addition, they can catalyze the reduction of disulfide bonds of Fe-S cluster proteins [27]. The thioredoxin reductase Trx2, along with the glutathione reductase Grx1, can reduce the disulfide bonds of OxyR protein [28]. At the same time, the OxyR regulator activates the transcription of these and other genes of oxidative stress, which triggers its own self-regulation [45].

Strain LB400 induced the sodB1 and fumC gene expression (fourfold) during paraquat exposure. Paraquat is a continuous source of superoxide radical, which explains the upregulation of the sodB1 gene in strain LB400. In E. coli, genes encoding for the scavenging enzyme SodB1 and ROS-resistant isoform of fumarate hydratase FumC are regulated by SoxRS in response to superoxide [29]. P. putida KT2440 exposed to this oxidizing agent also showed an expression increase of the sodA and fumC1 genes [57]. In strain LB400, the ROS-sensitive aconitase AcnB was downregulated in presence of paraquat (Fig. 6). Aconitate hydratase catalyzes the reversible isomerization of citrate and isocitrate via cis-aconitate in the Krebs cycle. In E. coli, the aconitase AcnB is susceptible to oxidative stress, whereas the aconitase AcnA, which is not dependent on a [4Fe–4S] cluster, plays an important role in survival during oxidative stress [18]. Superoxide radical produced by paraquat could damage [4Fe–4S] dehydratases such as AcnB [79]. This ROS could remove an Fe2+ ion from the [4Fe–4S] cluster, inactivating this enzyme [29]. The continuous generation of superoxide radical by paraquat leads to a recurring inactivation of the ROS-sensitive aconitase AcnB. To synthetize an enzyme that could be constantly inactivated represent an unnecessary metabolic burden for the cell, which may explain the downregulation of AcnB when strain LB400 was exposed to paraquat. The results of our study suggest a highly regulated network of key enzymatic components of the oxidative stress response in strain LB400.

Induction of the organic hydroperoxide resistance protein Ohr (BxeB2843) was observed upon paraquat and H2O2 exposure in strain LB400. A gene encoding the OhrR transcriptional regulator of the MarR family was identified upstream of the LB400 ohr gene (BxeB2843). Interestingly, the Ohr protein in P. aeruginosa is essential for resistance to organic hydroperoxides but not to H2O2 or paraquat, whereas its induction is independent of OxyR [56]. Exposure of Cupriavidus pinatubonensis AEO106 to H2O2 or sublethal concentrations of copper increased ROS accumulation, and upregulated the Ohr protein [73]. These results suggest that in these bacteria, the Ohr protein plays a key role on stress response during ROS accumulation.

ROS damage biomolecules including proteins [43, 63]. Threonine, arginine, proline, histidine, lysine and tryptophan are the main amino acids targets for oxidation by ROS [41]. LB400 cells exposed to paraquat showed an increase of 5,10-methylenetetrahydrofolate reductase (MetF) gene expression. MetF catalyzes the reduction of 5,10-methylenetetrahydrofolate to 5-methyltetrahydrofolate, a reaction involved in the synthesis of methionine from homocysteine [71]. During H2O2 exposure, strain LB400 upregulated homoserine reductase, which is involved in the synthesis of threonine from aspartate [81].

Other metabolic processes were triggered upon paraquat exposure. Proteomic analyses showed the upregulation of the phasin PhaP1 (BxeA1544), which is associated to the synthesis of PHA granules. PHAs are linear polyesters produced from sugars and fatty acids, which store carbon and energy in cytoplasmic granules under unbalanced conditions and increase the stress response [1, 19, 76]. The induction of the phaP1 gene expression during the synthesis of PHA granules by P. xenovorans LB400 has been reported [77].

High resistance to paraquat and H2O2. Strain LB400 showed a high resistance to the oxidizing agents paraquat and H2O2 compared to other bacteria. Strain LB400 was capable to resist until 16 mM H2O2 (Fig. 4B). In contrast, H2O2 1.25 mM and 2.5 mM inhibit the growth of E. coli and B. subtilis, respectively [8]. P. pseudomonas MPAO1 growth is inhibited by a H2O2 (0.6 mM) [20]. P. xenovorans LB400 showed much higher survival to paraquat compared to other bacteria. Strain LB400 exposed during 1 h to paraquat (16 mM) showed ~ 96% survival. In comparison, E. coli K-12 after 45 min-exposure to paraquat (0.1 mM), showed ~ 75% survival [49]. P. aeruginosa PAO1 was also less resistant to 0.2 mM paraquat [46]. P. xenovorans strain LB400 exposed to paraquat (20 mM) revealed an inhibition zone of ~ 38 mm2, while Synechocystis sp. PCC 6803 treated with paraquat (7.8 mM) showed an inhibition halo > 900 mm2 [53]. Nevertheless, B. subtilis CU1065 is more resistant than strain LB400, showing in presence of paraquat (0.5 M) an inhibition zone of 27 mm2 [10]. The reducing power (e.g., NADH, NADPH) generated by the metabolism in the cells is crucial to generate a response to oxidative stress. Interestingly, P. xenovorans LB400 and most of the Burkholderia sensu lato strains metabolize glucose and other sugars through the Entner-Douderoff (ED) pathway and the pentose phosphate (PP) pathway but not through the classical upper Embden-Meyerhof-Parnas (EMP) pathway because they lack the 6-phosphofructokinase pfk gene [1]. The ED pathway and the PP pathway are linked to the lower EMP pathway, producing NAD(P)H reducing power. The activation of the ED route increased the levels of NAD(P)H, which allowed to generate an effective response against the oxidative stress generated by the aerobic metabolism in P. putida KT2440 [15, 54]. The induction of NAD reducing enzymes, such as glyceraldehyde-3-phosphate dehydrogenase in E. coli, caused an accumulation of NADH via ferredoxin NADP reductase, resulting in the activation of the SoxRS response and an increased tolerance to paraquat [38]. Interestingly, completing the EMP pathway of P. putida KT2440 by expressing the E. coli pfkA, leads to a decrease in the resistance of this strain to the oxidizing compound diamide [15].

Overall, a broad antioxidant response was observed in P. xenovorans strain LB400 against paraquat and H2O2, in which the OxyR and the SoxR transcriptional regulators and other associated genes play crucial roles.


Diverse studies have shown that the model aromatic-degrading bacterium P. xenovorans LB400 suffers oxidative stress during the catabolism of aromatic compounds. In this study, we have shown that P. xenovorans strain LB400 possesses a broad repertoire of genes involved in oxidative and general stress response. In response to the oxidizing aromatic agent paraquat and hydrogen peroxide, strain LB400 displays differential broad oxidative stress responses. Notably, the functionality of OxyR and the SoxR transcriptional regulators were determined through their expression in P. xenovorans and P. aeruginosa strains. The broad range antioxidant response of P. xenovorans LB400, which involves the OxyR and the SoxR transcriptional regulators, may explain in part its metabolic versatility to degrade a wide range of aromatic compounds that causes oxidative stress.

Availability of data and materials

All data generated or analyzed during this study are included in this published article.



OxyR transcriptional regulator


SoxR transcriptional regulator


Reactive oxygen species




Phenazine methosulfate


3′-(p-Hydroxyphenyl) fluorescein






Collision induced dissociation


  1. Álvarez-Santullano N, Villegas P, Mardones MS, Durán RE, Donoso R, González A, et al. Genome-wide metabolic reconstruction of the synthesis of polyhydroxyalkanoates from sugars and fatty acids by Burkholderia sensu lato species. Microorganisms. 2021;9:1290.

    PubMed  PubMed Central  Google Scholar 

  2. Agulló L, Cámara B, Martínez P, Latorre-Reyes V, Seeger M. Response to (chloro)biphenyls of the PCB-degrader Burkholderia xenovorans LB400 involves stress proteins induced also by heat shock and oxidative stress. FEMS Microbiol Lett. 2007;267:167–75.

    PubMed  Google Scholar 

  3. Agulló L, Romero-Silva MJ, Domenech M, Seeger M. p-Cymene promotes its catabolism through the p-cymene and the p-cumate pathways, activates a stress response and reduces the biofilm formation in Burkholderia xenovorans LB400. PLoS ONE. 2017;12:e0169544.

    PubMed  PubMed Central  Google Scholar 

  4. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.

    CAS  Google Scholar 

  5. Apweiler R, Bairoch A, Wu CH. Protein sequence databases. Curr Opin Che Biol. 2004;8:76–80.

    CAS  Google Scholar 

  6. Bateman A, Martin MJ, O’Donovan C, Magrane M, Apweiler R, Alpi E, et al. UniProt: a hub for protein information. Nucleic Acids Res. 2015;43:D204–12.

    Google Scholar 

  7. Brown NL, Stoyanov JV, Kidd SP, Hobman JL. The MerR family of transcriptional regulators. FEMS Microbiol Rev. 2003;27:145–63.

    CAS  PubMed  Google Scholar 

  8. Brudzynski K, Abubaker K, Castle A. Re-examining the role of hydrogen peroxide in bacteriostatic and bactericidal activities of honey. Front Microbiol. 2011;2:213.

    PubMed  PubMed Central  Google Scholar 

  9. Cámara B, Herrera C, Gonzales M, Couve E, Hofer B, Seeger M. From PCBs to highly toxic metabolites by the biphenyl pathway. Env Microbiol. 2004;6:842–50.

    Google Scholar 

  10. Cao M, Moore CM, Helmann JD. Bacillus subtilis paraquat resistance is directed by σM, an extracytoplasmic function sigma factor, and is conferred by YqjL and BcrC. J Bacteriol. 2005;187:2948–56.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009;25:1972–3.

    PubMed  PubMed Central  Google Scholar 

  12. Caswell CC, Baumgartner JE, Martin DW, Roop RM. Characterization of the organic hydroperoxide resistance system of Brucella abortus 2308. J Bacteriol. 2012;194:5065–72.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. Cárdenas F. Biotransformation and effect of flavonoids on oxidative stress during the degradation of aromatic compounds in Burkholderia xenovorans LB400. PhD thesis in Biotechnology. Universidad Técnica Federico Santa María & Pontificia Universidad Católica de Valparaíso, Valparaíso, Chile. 2015.

  14. Chain PS, Denef VJ, Konstantinidis KT, Vergez LM, Agulló L, Latorre-Reyes V, et al. Burkholderia xenovorans LB400 harbors a multi-replicon, 9.73-Mbp genome shaped for versatility. Proc Natl Acad Sci USA. 2006;103:15280–7.

    PubMed  PubMed Central  Google Scholar 

  15. Chavarría M, Nikel PI, Pérez-Pantoja D, de Lorenzo V. The Entner-Doudoroff pathway empowers Pseudomonas putida KT2440 with a high tolerance to oxidative stress. Environ Microbiol. 2013;15:1772–85.

    PubMed  Google Scholar 

  16. Chirino B, Strahsburger E, Agulló L, González M, Seeger M. Genomic and functional analyses of the 2-aminophenol catabolic pathway and partial conversion of its substrate into picolinic acid in Burkholderia xenovorans LB400. PLoS ONE. 2013;8:e75746.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Coba de la Peña T, Redondo FJ, Fillat MF, Lucas MM, Pueyo JJ. Flavodoxin overexpression confers tolerance to oxidative stress in beneficial soil bacteria and improves survival in the presence of the herbicides paraquat and atrazine. J Appl Microbiol. 2013;115:236–46.

    PubMed  Google Scholar 

  18. Cunningham L, Gruer MJ, Guest JR. Transcriptional regulation of the aconitase genes (acnA and acnB) of Escherichia coli. Microbiology. 1997;143:3795–805.

    CAS  PubMed  Google Scholar 

  19. de Almeida A, Nikel PI, Giordano AM, Pettinari MJ. Effects of granule-associated protein PhaP on glycerol-dependent growth and polymer production in poly (3-hydroxybutyrate)-producing Escherichia coli. Appl Environ Microbiol. 2007;73:7912–6.

    PubMed  PubMed Central  Google Scholar 

  20. Deng X, Liang H, Ulanovskaya OA, Ji Q, Zhou T, Sun F, et al. Steady-state hydrogen peroxide induces glycolysis in Staphylococcus aureus and Pseudomonas aeruginosa. J Bacteriol. 2014;196:2499–513.

    PubMed  PubMed Central  Google Scholar 

  21. Dietrich LE, Price-Whelan A, Petersen A, Whiteley M, Newman DK. The phenazine pyocyanin is a terminal signalling factor in the quorum sensing network of Pseudomonas aeruginosa. Mol Microbiol. 2006;61:1308–21.

    CAS  PubMed  Google Scholar 

  22. Dietrich LE, Teal TK, Price-Whelan A, Newman DK. Redox-active antibiotics control gene expression and community behavior in divergent bacteria. Science. 2008;321:1203–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Dietrich LE, Kiley PJ. A shared mechanism of SoxR activation by redox-cycling compounds. Mol Microbiol. 2011;79:1119–22.

    CAS  PubMed  Google Scholar 

  24. Dietrich LE, Okegbe C, Price-Whelan A, Sakhtah H, Hunter C, Newman DK. Bacterial community morphogenesis is intimately linked to the intracellular redox state. J Bacteriol. 2013;195:1371–80.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Durán RE, Méndez V, Rodríguez-Castro L, Barra-Sanhueza B, Salvà Serra F, Moore ER, et al. Genomic and physiological traits of the marine bacterium Alcaligenes aquatilis QD168 isolated from Quintero Bay, Central Chile, reveal a robust adaptive response to environmental stressors. Front Microbiol. 2019;10:528.

    PubMed  PubMed Central  Google Scholar 

  26. Fuentes S, Méndez V, Aguila P, Seeger M. Bioremediation of petroleum hydrocarbons: catabolic genes, microbial communities and applications. Appl Microbiol Biotechnol. 2014;98:4781–94.

    CAS  PubMed  Google Scholar 

  27. Gleason FK, Holmgren A. Thioredoxin and related proteins in procaryotes. FEMS Microbiol Rev. 1988;4:271–97.

    CAS  PubMed  Google Scholar 

  28. Imlay JA. Cellular defenses against superoxide and hydrogen peroxide. Annu Rev Biochem. 2008;77:755–76.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Imlay JA. The molecular mechanisms and physiological consequences of oxidative stress: lessons from a model bacterium. Nature. 2013;11:443–54.

    CAS  Google Scholar 

  30. Imlay JA. Transcription factors that defend bacteria against reactive oxygen species. Ann Rev Microbiol. 2015;69:93–108.

    CAS  Google Scholar 

  31. Imlay JA. Where in the world do bacteria experience oxidative stress? Environ Microbiol. 2019;21:521–30.

    CAS  PubMed  Google Scholar 

  32. Jangiam W, Loprasert S, Smith DR, Tungpradabkul S. Burkholderia pseudomallei RpoS regulates OxyR and the katG-dpsA operon under conditions of oxidative stress. Microbiol Immunol. 2010;54:389–97.

    CAS  PubMed  Google Scholar 

  33. Kim J, Park W. Oxidative stress response in Pseudomonas putida. Appl Microbiol Biotechnol. 2014;98:6933–46.

    CAS  PubMed  Google Scholar 

  34. Kim J, Park C, Imlay JA, Park W. Lineage-specific SoxR-mediated regulation of an endoribonuclease protects non-enteric bacteria from redox-active compounds. J Biol Chem. 2017;292:121–33.

    CAS  PubMed  Google Scholar 

  35. Kovach ME, Elzer PH, Hill DS, Robertson GT, Farris MA, Martin RI, Peterson KM. Four new derivatives of the broad-host-range cloning vector pBBR1MCS, carrying different antibiotic-resistance cassettes. Gene. 1995;166:175–6.

    CAS  PubMed  Google Scholar 

  36. Kobayashi K, Tagawa S. Activation of SoxR-dependent transcription in Pseudomonas aeruginosa. J Biochem. 2004;136:607–15.

    CAS  PubMed  Google Scholar 

  37. Kolde R. pheatmap: Pretty Heatmaps in R package. R Package. version 1.0.8. 2015. Available at:

  38. Krapp AR, Humbert MV, Carrillo N. The soxRS response of Escherichia coli can be induced in the absence of oxidative stress and oxygen by modulation of NADPH content. Microbiology. 2011;157:957–65.

    CAS  PubMed  Google Scholar 

  39. Lanfear R, Frandsen PB, Wright AM, Senfeld T, Calcott B. Partitionfinder 2: new methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol Biol Evol. 2017;34:772–3.

    CAS  PubMed  Google Scholar 

  40. Larsson A. AliView: a fast and lightweight alignment viewer and editor for large datasets. Bioinformatics. 2014;30:3276–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  41. Lascano R, Muñoz N, Robert G, Rodriguez M, Melchiorre M, Trippi V, et al. Paraquat: an oxidative stress inducer. In: Mohammed NH, editor., et al., Herbicides-properties, synthesis and control of weeds. Shanghai: InTechChina; 2012. p. 135–48.

    Google Scholar 

  42. Lee SE, Seo JS, Keum YS, Lee KJ, Li QX. Fluoranthene metabolism and associated proteins in Mycobacterium sp. JS14. Proteomics. 2007;7:2059–69.

    CAS  PubMed  Google Scholar 

  43. Lemire J, Alhasawi A, Appanna VP, Tharmalingam S, Appanna VD. Metabolic defence against oxidative stress: the road less travelled so far. J Appl Microbiol. 2017;123:798–809.

    CAS  PubMed  Google Scholar 

  44. Livak JK, Schmittgen TD. Analysis of relative gene expression data using real time quantitative PCR and the 2-∆∆CT method. Methods. 2001;25:402–8.

    CAS  Google Scholar 

  45. Lu J, Holmgren A. The thioredoxin antioxidant system. Free Radic Biol Med. 2014;66:75–87.

    CAS  PubMed  Google Scholar 

  46. Ma JF, Hager PW, Howell ML, Phibbs PV, Hassett DJ. Cloning and characterization of the Pseudomonas aeruginosa zwf gene encoding glucose-6-phosphate dehydrogenase, an enzyme important in resistance to methyl viologen (paraquat). J Bacteriol. 1998;180:1741–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. Ma Y, Li L, Awasthi MK, Tian H, Lu M, Megharaj M, et al. Time-course transcriptome analysis reveals the mechanisms of Burkholderia sp. adaptation to high phenol concentrations. Appl Microbiol Biotechnol. 2020;104:5873–87.

    CAS  PubMed  Google Scholar 

  48. Martínez P, Agulló L, Hernández M, Seeger M. Chlorobenzoate inhibits growth and induces stress proteins in the PCB-degrading bacterium Burkholderia xenovorans LB400. Arch Microbiol. 2007;188:289–97.

    PubMed  Google Scholar 

  49. Membrillo-Hernández J, Kim SO, Cook GM, Poole RK. Paraquat regulation of hmp (flavohemoglobin) gene expression in Escherichia coli K-12 is SoxRS independent but modulated by sigma S. J Bacteriol. 1997;179:3164–70.

    PubMed  PubMed Central  Google Scholar 

  50. Méndez V, Agulló L, González M, Seeger M. The homogentisate and homoprotocatechuate central pathways are involved in 3- and 4-hydroxyphenylacetate degradation by Burkholderia xenovorans LB400. PLoS ONE. 2011;6:e17583.

    PubMed  PubMed Central  Google Scholar 

  51. Montero-Silva F, Durán N, Seeger M. Synthesis of extracellular gold nanoparticles using Cupriavidus metallidurans CH34 cells. IET Nanobiotech. 2018;12:40–6.

    Google Scholar 

  52. Moretti S, Armougom F, Wallace IM, Higgins DG, Jongeneel CV, Notredame C. The M-Coffee web server: a meta-method for computing multiple sequence alignments by combining alternative alignment methods. Nucleic Acids Res. 2007;35:W645–8.

    PubMed  PubMed Central  Google Scholar 

  53. Nefedova LN, Fantin YS, Zinchenko VV, Babykin MM. The prqA and mvrA genes encoding drug efflux proteins control resistance to methyl viologen in the cyanobacterium Synechocystis sp. PCC 6803. Russ J Genet. 2003;39:264–8.

    CAS  Google Scholar 

  54. Nikel PI, Fuhrer T, Chavarría M, Sánchez-Pascuala A, Sauer U, de Lorenzo V. Reconfiguration of metabolic fluxes in Pseudomonas putida as a response to sub-lethal oxidative stress. ISME J. 2021;15:1751–66.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. Ochsner UA, Vasil ML, Alsabbagh E, Parvatiyar K, Hassett DJ. Role of the Pseudomonas aeruginosa oxyR-recG operon in oxidative stress defense and DNA repair: OxyR-dependent regulation of katB-ankB, ahpB, and ahpC-ahpF. J Bacteriol. 2000;182:4533–44.

    CAS  PubMed  PubMed Central  Google Scholar 

  56. Ochsner UA, Hassett DJ, Vasil ML. Genetic and physiological characterization of ohr, encoding a protein involved in organic hydroperoxide resistance in Pseudomonas aeruginosa. J Bacteriol. 2001;183:773–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  57. Park W, Pena-Llopis S, Lee Y, Demple B. Regulation of superoxide stress in Pseudomonas putida KT2440 is different from the SoxR paradigm in Escherichia coli. Biochem Biophys Res Commun. 2006;341:51–6.

    CAS  PubMed  Google Scholar 

  58. Parsonage D, Youngblood DS, Sarma GN, Wood ZA, Karplus PA, Poole LB. Analysis of the link between enzymatic activity and oligomeric state in AhpC, a bacterial peroxiredoxin. Biochemistry. 2005;9:10583–92.

    Google Scholar 

  59. Perkins DN, Pappin DJ, Creasy DM, Cottrell JS. Probability-based protein identification by searching sequence databases using mass spectrometry data. Electrophoresis. 1999;20:3551–67.

    CAS  PubMed  Google Scholar 

  60. Pomposiello PJ, Demple B. Redox-operated genetic switches: the SoxR and OxyR transcription factors. Trends Biotechnol. 2001;19:109–14.

    CAS  PubMed  Google Scholar 

  61. Ponce BL, Latorre VK, González M, Seeger M. Antioxidant compounds improved PCB-degradation by Burkholderia xenovorans strain LB400. Enzyme Microb Technol. 2011;49:509–16.

    CAS  PubMed  Google Scholar 

  62. Rocha ER, Smith CJ. Role of the alkyl hydroperoxide reductase (ahpCF) gene in oxidative stress defense of the obligate anaerobe Bacteroides fragilis. J Bacteriol. 1999;181:5701–10.

    CAS  PubMed  PubMed Central  Google Scholar 

  63. Rodríguez-Castro L, Méndez V, Durán RE, Seeger M. Long-chain flavodoxin FldX1 improves Paraburkholderia xenovorans LB400 tolerance to oxidative stress caused by paraquat and H2O2. PLoS ONE. 2019;14:e0221881.

    PubMed  PubMed Central  Google Scholar 

  64. Ronquist F, Teslenko M, Van Der Mark P, Ayres DL, Darling A, Höhna S, et al. Mrbayes 3.2: efficient bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42.

    PubMed  PubMed Central  Google Scholar 

  65. Sakhtah H, Koyama L, Zhang Y, Morales DK, Fields BL, Price-Whelan A, Hogan DA, Shepard K, Dietrich LE. The Pseudomonas aeruginosa efflux pump MexGHI-OpmD transports a natural phenazine that controls gene expression and biofilm development. Proc Natl Acad Sci USA. 2016;113:E3538–47.

    CAS  PubMed  PubMed Central  Google Scholar 

  66. Sambrook J, Fritch EF, Maniatis T. Molecular cloning: a laboratory manual. 2nd ed. USA: Cold Spring Harbor Laboratory Press; 1989.

    Google Scholar 

  67. Sandoval JM, Arenas FA, Vásquez CC. Glucose-6-phosphate dehydrogenase protects Escherichia coli from tellurite-mediated oxidative stress. PLoS ONE. 2011;6:e25573.

    CAS  PubMed  PubMed Central  Google Scholar 

  68. Seaver LC, Imlay JA. Alkyl hydroperoxide reductase is the primary scavenger of endogenous hydrogen peroxide in Escherichia coli. J Bacteriol. 2001a;183:7173–81.

    CAS  PubMed  PubMed Central  Google Scholar 

  69. Seaver LC, Imlay JA. Hydrogen peroxide fluxes and compartmentalization inside growing Escherichia coli. J Bacteriol. 2001b;183:7182–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  70. Seeger M, Osorio G, Jerez CA. Phosphorylation of GroEL, DnaK and other proteins from Thiobacillus ferrooxidans grown under different conditions. FEMS Microbiol Lett. 1996;138:129–34.

    CAS  PubMed  Google Scholar 

  71. Sheppard CA, Trimmer EE, Matthews RG. Purification and properties of NADH-dependent 5,10-methylenetetrahydrofolate reductase (MetF) from Escherichia coli. J Bacteriol. 1999;181:718–25.

    CAS  PubMed  PubMed Central  Google Scholar 

  72. Sporer AJ, Beierschmitt C, Bendebury A, Zink KE, Price-Whelan A, Buzzeo MC, Sanchez LM, Dietrich L. Pseudomonas aeruginosa PumA acts on an endogenous phenazine to promote self-resistance. Microbiology. 2018;164:790–800.

    CAS  PubMed  PubMed Central  Google Scholar 

  73. Svenningsen NB, Damgaard M, Rasmussen M, Pérez-Pantoja D, Nybroe O, Nicolaisen MH. Cupriavidus pinatubonensis AEO106 deals with copper-induced oxidative stress before engaging in biodegradation of the herbicide 4-chloro-2-methylphenoxyacetic acid. BMC Microbiol. 2017;17:211.

    PubMed  PubMed Central  Google Scholar 

  74. Tam LT, Eymann C, Albrecht D, Sietmann R, Schauer F, Hecker M, Antelmann H. Differential gene expression in response to phenol and catechol reveals different metabolic activities for the degradation of aromatic compounds in Bacillus subtilis. Environ Microbiol. 2006;8:1408–27.

    CAS  Google Scholar 

  75. Tocheva EI, Fortin PD, Eltis LD, Murphy ME. Structures of ternary complexes of BphK, a bacterial glutathione S-transferase that reductively dechlorinates polychlorinated biphenyl metabolites. J Biol Chem. 2006;281:30933–40.

    CAS  PubMed  Google Scholar 

  76. Urtuvia V, Villegas P, González M, Seeger M. Bacterial production of the biodegradable plastics polyhydroxyalkanoates. Int J Biol Macromol. 2014;70:208–13.

    CAS  PubMed  Google Scholar 

  77. Urtuvia V, Villegas P, Fuentes S, González M, Seeger M. Burkholderia xenovorans LB400 possesses a functional polyhydroxyalkanoate anabolic pathway encoded by the pha genes and synthesizes poly (3-hydroxybutyrate) under nitrogen-limiting conditions. Int Microbiol. 2018;21:47–57.

    CAS  PubMed  Google Scholar 

  78. Van Driessche G, Vandenberghe I, Devreese B, Samyn B, Meyer TE, Leigh R, et al. Amino acid sequences and distribution of high-potential iron-sulfur proteins that donate electrons to the photosynthetic reaction center in phototropic proteobacteria. J Mol Evol. 2003;57:181–99.

    PubMed  Google Scholar 

  79. Varghese S, Tang Y, Imlay JA. Contrasting sensitivities of Escherichia coli aconitases A and B to oxidation and iron depletion. J Bacteriol. 2003;185:221–30.

    CAS  PubMed  PubMed Central  Google Scholar 

  80. Wang Q, Lu X, Yang H, Yan H, Wen Y. Redox-sensitive transcriptional regulator SoxR directly controls antibiotic production, development and thiol-oxidative stress response in Streptomyces avermitilis. Microb Biotechnol. 2021.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Zhang Y, Meng Q, Ma H, Liu Y, Cao G, Zhang X, et al. Determination of key enzymes for threonine synthesis through in vitro metabolic pathway analysis. Microb Cell Fact. 2015;14:86.

    PubMed  PubMed Central  Google Scholar 

  82. Zhang B, Gu H, Yang Y, Bai H, Zhao C, Si M, Su T, Shen X. Molecular mechanisms of AhpC in resistance to oxidative stress in Burkholderia thailandensis. Front Microbiol. 2019;10:1483.

    PubMed  PubMed Central  Google Scholar 

Download references


Authors acknowledge Dr. Jonas Perales and Dr. Patricia Cuervo (Fundación Oswaldo Cruz-Fiocruz, Rio de Janeiro, Brazil) for the proteomic analyses, Dr. Lars Dietrich (Department of Biological Sciences, Columbia University, USA) for providing bacterial strains and the support during the construction of P. aeruginosa recombinant strains in his laboratory, and Dr. Franco Cárdenas for the kind donation of some oligonucleotides employed in the gene expression analyses.


M.S., V.M. and L.R-C. dedicated this article in memory of our inspiring professor, generous colleague and friend Professor Dr. Claudio C. Vásquez Guzmán (1952–2020), who inspired and advised our research in microbiology and oxidative stress until his final departure. Professor Dr. Claudio C. Vásquez was member of the PhD thesis committees of Valentina Méndez and Laura Rodríguez-Castro, whose PhD theses originated part of the results reported in this article.


This study was supported by RIABIN and ANID PhD fellowships (V.M., L.R-C.), ANID Gastos Operacionales scholarships (V.M., L.R-C.), FONDECYT 1200756 & 1151174 (M.S., V.M., L.R-C., R.E.D.) (, ANID Programa de Investigación Asociativa (PIA) Anillo GAMBIO ACT172128 and USM PI_M_2020_43 (R.E.D., M.S.) ( grants.

Author information




VM, MS and GP designed the experiments. V.M. performed the experiments. All authors analyzed and interpreted the data. VM, LR-C, RED and MS were main contributors in writing the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Michael Seeger.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent to publish

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1

. Oligonucleotides used as primers for oxidative stress genes of P. xenovorans LB400. Table S2. Identification of oxidative stress resistance and iron metabolism proteins in P. xenovorans LB400 genome. Table S3. Proteins upregulated in P. xenovorans LB400 during exposure to oxidizing agents. Table S4. Proteins downregulated in P. xenovorans LB400 during exposure to oxidizing agents.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Méndez, V., Rodríguez-Castro, L., Durán, R.E. et al. The OxyR and SoxR transcriptional regulators are involved in a broad oxidative stress response in Paraburkholderia xenovorans LB400. Biol Res 55, 7 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Oxidative stress
  • Paraburkholderia xenovorans
  • Paraquat
  • Hydrogen peroxide
  • Superoxide
  • OxyR
  • SoxR