Next Article in Journal
Phenolic Contents and Antioxidant Activity of Citrullus colocynthis Fruits, Growing in the Hot Arid Desert of the UAE, Influenced by the Fruit Parts, Accessions, and Seasons of Fruit Collection
Next Article in Special Issue
Reciprocal Regulation of Shh Trafficking and H2O2 Levels via a Noncanonical BOC-Rac1 Pathway
Previous Article in Journal
Exploring the Diversity of the Thioredoxin Systems in Cyanobacteria
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Antioxidants
Volume 11
Issue 4
10.3390/antiox11040655
antioxidants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptomic Analysis of E. coli after Exposure to a Sublethal Concentration of Hydrogen Peroxide Revealed a Coordinated Up-Regulation of the Cysteine Biosynthesis Pathway

by
Myriam Roth
1,*,
Vincent Jaquet
1,2,
Sylvain Lemeille
1,
Eve-Julie Bonetti
3,
Yves Cambet
2,
Patrice François
3,† and
Karl-Heinz Krause
1,†
1
Department of Pathology and Immunology, Medical School, University of Geneva, 1211 Geneva, Switzerland
2
REaders, Assay Development & Screening Unit (READS Unit), Faculty of Medecine, University of Geneva, 1211 Geneva, Switzerland
3
Genomic Research Laboratory, Infectious Diseases Service, University Hospitals Geneva Medical Center, Michel-Servet 1, 1211 Geneva, Switzerland
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Antioxidants 2022, 11(4), 655; https://doi.org/10.3390/antiox11040655
Submission received: 8 March 2022 / Revised: 24 March 2022 / Accepted: 25 March 2022 / Published: 28 March 2022
(This article belongs to the Special Issue Hydrogen Peroxide Signaling in Physiology and Pathology)
Browse Figures
Review Reports Versions Notes

Abstract

:
Hydrogen peroxide (H2O2) is a key defense component of host-microbe interaction. However, H2O2 concentrations generated by immune cells or epithelia are usually insufficient for bacterial killing and rather modulate bacterial responses. Here, we investigated the impact of sublethal H2O2 concentration on gene expression of E. coli BW25113 after 10 and 60 min of exposure. RNA-seq analysis revealed that approximately 12% of bacterial genes were strongly dysregulated 10 min following exposure to 2.5 mM H2O2. H2O2 exposure led to the activation of a specific antioxidant response and a general stress response. The latter was characterized by a transient down-regulation of genes involved in general metabolism, such as nucleic acid biosynthesis and translation, with a striking and coordinated down-regulation of genes involved in ribosome formation, and a sustained up-regulation of the SOS response. We confirmed the rapid transient and specific response mediated by the transcription factor OxyR leading to up-regulation of antioxidant systems, including the catalase-encoding gene (katG), that rapidly degrade extracellular H2O2 and promote bacterial survival. We documented a strong and transient up-regulation of genes involved in sulfur metabolism and cysteine biosynthesis, which are under the control of the transcription factor CysB. This strong specific transcriptional response to H2O2 exposure had no apparent impact on bacterial survival, but possibly replenishes the stores of oxidized cysteine and glutathione. In summary, our results demonstrate that different stress response mechanisms are activated by H2O2 exposure and highlight the cysteine synthesis as an antioxidant response in E. coli.

1. Introduction

Escherichia coli, a Gram-negative rod in the family Enterobacteriaceae, is a commensal member of the normal gut flora. However, E. coli also represents a major source of infection. For example, E. coli causes 75–85% of urinary tract infections [1], and is also the most frequent pathogen responsible for bacteremia (30% of the bacteremia in Switzerland) [2]. The emergence of antibiotic resistance, especially against fluoroquinolones and 3rd and 4th generations of cephalosporins [3], as well as against beta-lactam antibiotics through the production of extended spectrum beta-lactamases (ESBLs) and carbapenemases, imply that E. coli is a priority target for the development of new classes of antibiotics [4].
Reactive oxygen species (ROS) are oxygen derivatives that play important roles in chemistry and biology [5]. ROS can be generated by physico-chemical interactions (e.g., radiation), by eukaryotic cells [6], where they play a role in signaling, biosynthesis and host defense, as well as by bacteria, where they might participate in microbial warfare. In many instances, the initial step in the generation of ROS is the superoxide radical anion O2•− (generated through single electron transfer to molecular oxygen), but the latter is rapidly dismutated to the more stable and hence quantitatively predominant hydrogen peroxide (H2O2). ROS play an important role in the regulation of bacterial colonization and pathogenicity, as exemplified by the innate immune system. Phagocytes, especially neutrophils and macrophages, produce large amounts of ROS following the activation of the NADPH oxidase 2 (NOX2), a process necessary for efficient bacterial killing [7]. Loss of function of NOX2 leads to a primary immunodeficiency called chronic granulomatous disease (CGD), where patients suffer from severe and chronic infections [8]. Other ROS-producing NADPH oxidases (NOX1, DUOX1/2) are found in digestive and respiratory epithelia where they control bacterial colonization [9] and bacterial gene expression [10]. Interestingly, the regulation of bacterial colonization by ROS is not limited to enzymes of the host defense, as the production of ROS by the natural human flora limits the colonization by pathogenic bacteria. Indeed, H2O2-producing Lactobacillus inhibit the growth of uropathogenic E. coli strains [11]. Such H2O2-producing Lactobacillus are found in the vaginal microbiota [12,13] and in the urobiome [14,15,16]. Colonization with H2O2-producing bacteria limits the expansion of potentially pathogenic E. coli and therefore represents a protective factor against urinary tract infections [13,15].
While it is evident that ROS plays a role in controlling bacterial colonization and pathogenicity, the underlying mechanisms are only partially understood. For example, H2O2 is able to kill bacteria in vitro, but only at high (millimolar) concentrations. The resistance of bacteria to cytotoxic concentrations of H2O2 is explained by the fact that bacteria synthesize H2O2-scavenging enzymes such as catalases (KatG, KatE) and alkyl peroxidase (AhpC) [17]. Therefore, bactericidal H2O2 concentrations are probably not easily reached in the human body [18].
Previous studies in E. coli have characterized ROS-specific regulators of gene expression such as the “hydrogen peroxide-inducible genes activator” OxyR and the “superoxide response regulon” SoxRS [19]. OxyR is a transcription factor that senses H2O2 at a concentration as small as 100 nM [20]. Oxidation-sensitive cysteine residues are oxidized and form an intramolecular disulfide bridge, modifying the conformation of the protein and allowing the binding to functional DNA-binding sites and induction of the transcription of specific antioxidant genes. The last detailed study of the OxyR regulon in E. coli showed that OxyR regulates 38 genes in 28 transcription units while RNA-seq analysis showed that 110 dysregulated genes were differentially expressed under oxidative stress in a oxyR deleted strain [21]. However the sole regulation of this transcription factor cannot explain the massive gene dysregulation under H2O2 exposure [22,23]. The improvements in next-generation sequencing techniques, bioinformatics analytical tools, and genome annotations allow detailed exploration of these questions.
In this study, we performed RNA-seq analysis of E. coli following exposure to sublethal concentration of H2O2. We observed major gene dysregulation 10 min after the addition of H2O2, in H2O2-specific regulons as well as in genes involved in more general stress responses. We identified sulfur metabolism, including the sulfate assimilation and the cysteine biosynthesis pathways, as the most significantly up-regulated metabolic pathways at 10 min. The H2O2-enhanced expression of genes from these pathways was dependent on the transcription factor CysB. Our results revealed the importance of the sulfur metabolism under oxidative stress and the role of CysB in the immediate H2O2 response in E. coli.

2. Materials and Methods

2.1. Bacterial Strains and Growth Conditions

All bacterial strains used in this study appear in Table 1. E. coli strains were cultured at 37 °C in Luria-Bertani (LB) (Becton & Dickinson, Sparks, MD, USA) broth or on Luria-Bertani Agar (Becton & Dickinson, Sparks, MD, USA). Wild-type (WT) strain refers to BW25113.

2.2. Growth Curves—Exposure to H2O2

An overnight culture of E. coli was normalized to 1.0 Mc Farland using a Densimat (bioMérieux, Marcy-l’Étoile, France) and further diluted 1:10 in fresh LB. This bacterial culture was grown in a volume of 1 mL in a 24-well plate (ThermoFisher, Waltham, MA, USA) and incubated at 37 °C with 5 mm orbital shaking in an Infinite 200PRO plate reader (Tecan, Männedorf, Switzerland). Absorbance was measured every 10 min at an optical density of 595 nm (OD595nm). At early exponential phase (OD595nm = 0.2), H2O2 was added from a 35% w/w (Acros Organics, Geel, Belgium) solution freshly diluted at 11-fold concentration in LB and serially diluted 1:2 in LB. A volume of 100 μL was added to the 1 mL bacterial culture (total volume per wells: 1100 μL) to reach final concentrations ranging from 80 mM to 0.078 mM H2O2.
The concentration of H2O2 stock was assessed before the experiments by using a WPA Biowave DNA spectrophotometer (Biochrom, Cambridge, UK) at 240 nm with quartz cuvette with a molar extinction coefficient of 43.6.
The Doubling Time Software v3.1.0 (http://www.doublingtime.com, accessed on 18 February 2022) was used to calculate the growth rate after H2O2 addition.

2.3. H2O2 Measurments by Amplex Red

Bacterial cultures of E. coli (see conditions above) were grown to OD595nm = 0.2 and H2O2 was added at a final concentration of 2.5 mM. At indicated time points, 400 µL were taken from each well and centrifuged for 4 min at 6000× g to remove bacteria. The supernatant was diluted 1:100 in LB. 100 µL of each sample were transferred into a 96-well black plate with clear bottom (Corning, Corning, NY, USA) in technical triplicates. Amplex Red (Invitrogen, Carlsbad, CA, USA) was used to detect H2O2 according to manufacturer’s instructions, by the addition of 100 µL of Amplex Red mix to each well for a final concentration of 55 µM Amplex Red and 10−2 UI/mL Horseradish peroxidase. The plate was incubated 1 h at 37 °C and the absorbance at 575 nm was read in a Spectramax Paradigm (Molecular Devices, San Jose, CA, USA).
A H2O2 calibration curve was generated by 1:2 serial dilutions of H2O2 in sterile LB (from 0.125 mM to 2.4 × 10−4 mM) and used to calculate the H2O2 concentration of the samples using linear regression.

2.4. RNA Isolation, RNA-Seq

2 mL of bacterial cultures (see conditions above) were harvested 10 and 60 min after H2O2 addition. RNA was extracted using the TRIzol® Max™ Bacterial RNA Isolation Kit (ThermoFisher, Waltham, MA, USA) according to the manufacturer’s instructions. After elution in RNA-free water, DNA digestion was performed in liquid by adding 2.5 µL DNase I and 10 µL Buffer RDD (Qiagen) in a total volume of 100 µL for 10 min at room temperature. Potential contaminants were removed by using the RNeasy mini kit plus (Qiagen, Germantown, MD, USA) according to the manufacturer’s instructions. Analysis of total RNA quality and quantity was performed using Nanodrop and Agilent 2100 BioAnalyzer (Agilent Technologies, Santa Clara, CA, USA).
For RNA sequencing, 1.2 µg of total RNA was ribo-depleted with the RiboMinus™ Bacteria 2.0 Transcriptome Isolation kit (ThermoFisher Scientific, Waltham, MA, USA). The truseq total RNA stranded kit (Illumina, San Diego, CA, USA) was then used for the library preparation. Library quantity was measured by Qubit and quality was assessed with a Tapestation on a DNA High sensitivity chip (Agilent Technologies, Santa Clara, CA, USA). Libraries were pooled at equimolarity for clustering. Single-read sequencing (100 bases) was performed using the SBS chemistry on an Illumina HiSeq 4000 sequencer. Raw data have been deposited at ENA (https://www.ebi.ac.uk/ena, accessed on 4 March 2022) under the following accession number: PRJEB51098.
The RNA-seq data were analyzed as follows. Fastq reads were mapped to the ENSEMBL reference genome (Ecolibw25113.48) using STAR (version 2.4.0j) [25] with standard settings, except that any reads mapping to more than one location in the genome (ambiguous reads) were discarded (m = 1). All reads were reported using featureCounts version 1.4.6-p1 [26]. Levels of expression were reported as raw counts and, in parallel, normalized in RPKM in order to filter out genes with low expression values (<1 RPKM) before calling genes as differentially expressed. Library size normalizations and differential gene expression calculations were performed using the package edgeR [27] designed for the R software [22]. Only genes having a significant fold-change (Benjamini-Hochberg corrected p-value < 0.01) were considered for the rest of the RNA-seq analysis.
Genes under the regulation (regulon) of selected transcription factors were investigated in the Ecocyc database [28] and used to generate heatmaps.
Genes were categorized in different expression patterns depending on if they were significantly up-regulated, down-regulated, or unchanged at the 2 time points, leading to 8 different groups. STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) web software [29] was used to analyze physical or functional interactions in each of the 8 groups of genes.
GO term enrichment was performed using homemade scripts for the R software [30]. Gene set enrichment analysis (GSEA) was performed, without threshold of p-value or fold-change, using the KEGG metabolic pathways (KEGG http://www.genome.jp/kegg/, accessed on 15 December 2021) relative to Ecolibw25113.48 to generate gene sets. Genes were ranked by their calculated fold-changes (decreasing ranking). A gene set analysis using the GSEA package Version 2.2 [31,32] from the Broad Institute (MIT, Cambridge, MA, USA) was used to analyze the pattern of differential gene expression between the two groups. Gene set permutations were performed 1000 times for each analysis. The Normalized Enrichment Score (NES) was calculated for each gene set. GSEA results with a nominal FDR < 0.05 and abs(NES) > 1 were considered significant.
ISMARA analysis [33] was performed to retrieve transcription factors implicated in causing changes in gene expression. Raw fastq data were uploaded to the ISMARA web server (Swiss Institute of Bioinformatics, Basel, Switzerland). Four replicates were averaged for each condition.

2.5. qRT-PCR

qRT-PCR was performed on RNA samples prepared as described above. The cDNA was produced by reverse-transcribing 400 ng of total RNA using a mix of random hexamers and oligo d(T) primers and Primescript reverse transcriptase enzyme (Takara Bio, San Jose, CA, USA). The efficiency of each pair of primers was tested with serial dilutions of cDNA. Oligonucleotides are indicated in Table 2. PCR reactions (10 µL volume) contained 1:20 diluted cDNA, 2 × Power SYBR Green Master Mix (Applied Biosystems by ThermoFisher, Waltham, MA, USA), and 300 nM of forward and reverse primers. PCRs were performed on a SDS 7900 HT instrument (Applied Biosystems by ThermoFisher, Waltham, MA, USA) with the following parameters: 50 °C for two minutes, 95 °C for ten minutes, and 45 cycles of 95 °C 15 s, 60 °C one minute. Each reaction was performed in three replicates on 384-well plate. Raw Ct values obtained with SDS 2.2 (Applied Biosystems by Thermo Fisher, Waltham, MA, USA) were imported into Excel and normalization factors were calculated using the GeNorm method as described by Vandesompele et al. [34].

2.6. P1 Transduction

cysB and katG genes were deleted in E. coli BW25113 using P1 transduction from the Keio library strain as previously described [36]. The mutants were validated using PCR with appropriate gene-specific primers (Table 2).

2.7. H2O2 Sensivity Testing by Disk Diffusion Assay

To assess the sensitivity to H2O2, we adapted the EUCAST methods of disc diffusion assay to H2O2 testing [37]. Briefly, an overnight culture of bacteria was diluted in LB to McFarland 0.5 using a Densimat (bioMérieux Marcy-l’Étoile, France) and this was used to inoculate LB agar plates. A sterile disk of cellulose of 6 mm diameter from blotting paper was placed on the plate and 10 mL of 1 M H2O2 diluted in sterile water was added in the center of the disk. Plates were incubated at 37 °C for 18 h and the diameter of inhibition was measured.

2.8. Softwares

Graphpad Prism v.9.3.1 for Windows (GraphPad Software, San Diego, CA, USA) was used for data processing, graph plotting and statistical analysis. Inkscape v.1.1.1 for Windows was used for image editing (https://inkscape.org).

3. Results

3.1. Sublethal H2O2 Concentration Determination

To determine the sublethal H2O2 concentration, the WT E. coli BW25113 strain was grown in LB medium until an OD of 0.2 and was then exposed to a range of H2O2 concentrations (Figure 1A). H2O2 was added during the exponential phase of growth where the bacterial population is increasing and where quorum-sensing signaling that could mask gene dysregulation is minimal. We defined the sublethal concentration as the highest concentration with a similar profile to the untreated control (Figure 1B). This corresponded to 2.5 mM H2O2 and the induction of a growth lag time of approximatively 30 min.
We followed the degradation of 2.5 mM H2O2 by E. coli BW25113 using the H2O2-sensitive probe Amplex Red/HRP. We observed rapid and complete degradation of the exogenous H2O2 after 25 min (Figure 1C). This rapid degradation is relevant for our experiments, as results obtained 10 min after the addition of H2O2 reflected direct responses to high H2O2 concentrations, while results obtained after 60 min reflected the long-term impact of transient H2O2 exposure.
We next investigated the impact of increasing H2O2 concentrations on a katG deficient mutant. The bifunctional hydroperoxidase I (HPI) encoded by katG possesses both catalase and peroxidase. The simplified term “catalase” is further used in the manuscript. This enzyme is E. coli’s principal H2O2 scavenger at high H2O2 concentrations. As expected [38], the katG deficient mutant was at least 10 times more sensitive to H2O2 compared to wild-type (Figure 1D).

3.2. RNA-seq under Sublethal H2O2 Exposure

Figure 2 shows the changes in gene expression caused by H2O2 exposure after 10 min and 60 min compared to a H2O2-free condition at identical time points (Figure 2A,B). Each condition was compared to an untreated control at identical time points. A massive gene dysregulation was observed 10 min after the addition of H2O2. In order to focus on the most relevant phenomena, we used stringent criteria to determine significance (adjusted p-value < 0.01, Fold-change ≥ 5). With such criteria, 12.2% of all transcripts of E. coli were significantly dysregulated (5.3% up- and 6.8% down-regulated) in the H2O2-treated condition compared to untreated cells after 10 min, and 2.4% of transcripts were significantly dysregulated in the respective conditions at 60 min (1.5% up- and 0.9% down-regulated).
Multi-dimensional Scaling (MDS) was performed on the top 500 genes that distinguished the conditions (Figure 2C). In this analysis, distances on the plot represent the variation of expression between samples. This analysis shows that replicates of our experiments (n = 4) were near each other, ensuring a high reproducibility.
There were marked differences in gene expression after 10 and 60 min under control conditions. These differences could be explained by bacterial growth from exponential phase to mid-exponential or early stationary phase. This is supported by the fact that quorum-sensing is one of the most up-regulated pathways while flagellar assembly and ribosome formation were the most down-regulated pathways in the 60 min control condition compared to the 10 min control condition (Figure S3).
Gene expression 10 min after H2O2 exposure was highly distinct from other conditions. In contrast, at 60 min, gene expression was more similar to that of the control condition. This is most likely due to the fact that H2O2 was rapidly degraded by E. coli (Figure 1C). This is corroborated by results depicted in Figure 1A, showing bacterial regrowth after a lag period of approximatively 30 min after addition of 2.5 mM H2O2. In addition, time-dependent changes of media composition due to bacterial metabolism may also affect gene expression.
OxyR is a well described regulator of gene expression known to be activated by H2O2. To validate our experimental approach, we therefore investigated H2O2–induced changes in OxyR-regulated genes [28]. A total of 19 (58%) of the OxyR-up-regulated genes were found significantly up-regulated upon H2O2 exposure compared to the untreated control, while there was no marked down-regulation of OxyR-down-regulated genes. Interestingly, most of the OxyR-regulated genes returned to baseline at 60 min (Figure 2D). Thus, the activity of the OxyR regulon appears to be correlated with H2O2 concentration (Figure 1C). These results demonstrate that—as expected—the OxyR pathway was activated under our experimental conditions. However, the absence of a dysregulation of some of the OxyR genes (e.g., dsbG, rscC, uxuA etc.) suggests that some of these genes depend on other transcription factors or are not activated at the H2O2 concentration used in this study.
In terms of oxidative stress response, the SoxRS system, which is implicated in the defense against superoxide, redox cycling compounds and H2O2, was also activated [19]. SoxR is activated by the oxidation of its iron-sulfur cluster and allows the transcription of the small RNA SoxS which activates the expression of approximatively 40 genes [39] (Figure S1). However, HypT, a recently discovered transcription factor activated by hypochlorite (HOCl) [40], a key ROS of the neutrophil host defense, did not induce its target genes, confirming the specificity of the OxyR and SoxRS activation by H2O2.
One of the best documented genes of the OxyR regulon is katG, which encodes E. coli’s primary H2O2 scavenger at high H2O2 concentrations. We validated katG induction by H2O2 by measuring its mRNA expression by qRT-PCR following treatment with increasing concentrations of H2O2 (Figure 2E). Sub-millimolar concentrations of H2O2 were sufficient for activation of katG gene expression and katG expression increased in a concentration-dependent manner until it was massively up-regulated at the sublethal concentration of 2.5 mM H2O2. These results validated the activation of OxyR under our experimental conditions, and, together with Figure 1D, confirmed that the relative H2O2 resistance of the BW25113 strain of E. coli is—at least in part—due to induction of the expression of katG.

3.3. Comparison between Gene Expression at 10 and 60 min after Sublethal H2O2 Treatment

In order to better define the global response to H2O2, we analyzed the time dependence of the impact of H2O2 on gene expression. Figure 3 shows the impact of H2O2 on differential gene expression after 10 min (H2O2 vs. control) compared to 60 min (H2O2 vs. control). A vast majority of genes (87%) showed absence of dysregulation (see white square of Figure 3A). Each colored square represents H2O2-dependent gene dysregulated at different time points. For example, the blue square on the left top represents genes that were down-regulated at 10 min but up-regulated after 60 min. Figure 3B provides the number of dysregulated genes at the two different time points. The most represented patterns were genes that were either up- or down-regulated at 10 min and displayed no difference at 60 min. Only 13 genes displayed a sustained up-regulation and only one gene (carB, implicated in the first reaction of pyrimidine biosynthesis) was down-regulated at 10 and 60 min. Selected genes from each pattern are represented in Figure S2.

3.4. Determination of the Functionnal Interaction of Genes Following the Same Kinetic Pattern

We determined whether each group within the Venn diagram depicted in Figure 3B included genes that shared physical or functional interactions using the STRING database [29].
Genes that were down-regulated at 10 min and whose expression then normalized after 60 min were massively enriched for genes coding for ribosomal proteins. The group of genes displaying an up-regulation at 10 min followed by no differences at 60 min represented the largest category with 180 genes. There was a significant enrichment of genes from the gene ontology response to oxidative stress, which is consistent with the OxyR regulon analysis (Figure 1D), and iron homeostasis. Genes that showed a down-regulation only at 60 min were implicated in nitrate assimilation. Interestingly, the genes that had sustained up-regulation were almost all implicated in the SOS-response pathway.
The fact that most changes occurred at 10 min and that genes already related to oxidative stress defense, such as the OxyR regulon depicted in Figure 1D, were up-regulated at this time point suggests that other genes following the same expression pattern can also be implicated in the response to oxidative stress.
This analysis highlights the fact that most changes occurred within minutes following H2O2 exposure and that, although it was relevant to analyze gene expression at the 60 min time point, several of the observed changes are in fact not directly due to H2O2 treatment. Thus, we mostly focused the following analyses on the direct response to H2O2 (10 min time point).

3.5. Identification of Transcription Factors Involved in ROS-Dependent Gene Regulation

The ISMARA (Integrated System for Motif Activity Response Analysis) software (Swiss Institute of Bioinformatics, Basel, Switzerland) was used on RNA-seq data to detect motifs in the promoters of dysregulated genes and deduce the involvement of given regulators (such as transcription factors, small RNA, RNA-polymerase subunit). Figure 4 shows 10 transcription factors that were identified as the most significant hits by the ISMARA analysis at the 10 min time point. The significance of the expression of the respective target genes at 10 min (dark purple) and 60 min (light purple) after the addition of H2O2 are illustrated. After 10 min, the target genes of 8 transcription factors were up-regulated, while the target genes of two regulators were down-regulated. At the 60 min time point, target genes of 8 of the regulators had returned to baseline levels of expression. Interestingly, the LexA target genes (63 genes distributed on 42 transcription units [28]) remained significantly up-regulated even after 60 min, while the RpoS target genes (up to 10% of the genes in E. coli [41]) were up-regulated after 10 min, but down-regulated after 60 min.
Note that the ISMARA results reflect the changes of targeted genes under the regulation of a given transcription factor and not necessarily its activity. For example, Fur target genes are up-regulated. However, as Fur is a transcription repressor, our results suggest that Fur is inhibited by H2O2.

3.6. Pathway Enrichment Analysis of H2O2—Dysregulated Genes after 10 min

As metabolism seemed highly impacted, we also performed a Gene Set Enrichment Analysis (GSEA) using the KEGG pathway database (Figure 5). Ribosome, fatty acid metabolism and oxidative phosphorylation pathways were the most down-regulated pathways. On the other side, the most up-regulated metabolic pathways were the sulfur metabolism and the siderophore biosynthesis pathways, which concurs with the up-regulation of the Fur and CysB regulated genes documented in Figure 4.

3.7. Investigation of the Sulfur Metabolism of E. coli under Sublethal H2O2 Concentration

The regulation of the sulfate assimilation and cysteine biosynthesis pathways have so far received little attention in the context of oxidative stress. We observed a strong up-regulation of genes regulated by CysB, the regulator of the cysteine regulon, 10 min after H2O2 exposure, which returned to baseline after 60 min (Figure 6A). Genes involved in sulfate assimilation (cysA, cysW, cysU, cysP, spb, cysN, cysD, cysI, cysJ, cysH) and cysteine biosynthesis (cysM, cysK), that are distributed in 4 distinct transcription units, showed the highest up-regulation, suggesting coordinated transcriptional regulation.
In order to test the significance of the CysB regulon in the presence of H2O2, a single deletion mutant of cysB was tested in a H2O2 disk diffusion assay (Figure 6B,C). The cysB mutant did not exhibit an increased sensitivity toward H2O2 compared to the katG deleted mutant used as control. However, the cysB deletion led to an impaired fitness, indicated by the formation of small and translucent colonies on agar plates and a slower growth rate in liquid medium, confirming the importance of CysB for optimal growth (Figure S3).
We used qRT-PCR to investigate the expression of several CysB-regulated genes in response to increasing concentration of H2O2 (Figure 7). We observed a strong and dose-dependent up-regulation of cysI, cysJ, cysH, cysN and tcyP, which was completely mitigated in the cysB mutant at 2.5 mM H2O2. Thus, the induction of these genes is specific for H2O2 and is totally dependent on cysB. Note that the cysE gene, which is known to be independent of CysB [42], was indeed not influenced by cysB deletion. The cysB expression level did not show a H2O2 dose-dependent regulation and was similar to the control at 2.5 mM, suggesting that its activity does not depend on itself or OxyR.
Interestingly, the expression level of katG appeared to be up-regulated in the cysB deletion mutant compared to the WT after exposition to 2.5 mM H2O2.
Altogether, our study provides, for the first time, a quantitative snapshot of genes affected by H2O2 in E. coli and identified that the CysB regulon represents a specific response to H2O2.

4. Discussion

In this study, we analyzed the impact of sublethal H2O2 concentrations on gene expression in E. coli. Consistent with previous studies, the most important changes were the activation of the general stress responses and redox-specific stress responses. In particular, we highlighted that several transcription factors were regulating the H2O2 response. Most of these transcription factor activities were transient; however, the response to DNA damage (LexA) was maintained over a prolonged period. The main finding of this study was that sulfur metabolism was the most significantly up-regulated pathway and that H2O2 strongly and specifically induced the activation of the transcription factor CysB.
A role of ROS in the regulation of bacterial growth and survival is beyond doubt. H2O2 is the most abundant product of the phagocyte NADPH oxidase, as well as of certain H2O2 -generating Lactobacilli. Thus, the medical relevance of the impact of H2O2 on bacterial pathogenesis is clear. However other reactive species are also involved in host response, in particular halogenated oxidants generated by myeloperoxidase and lactoperoxidase, such as hypocyanite, hypobromite and hypochlorite [43]. The lethal toxicity of H2O2 to E. coli occurred at relatively high concentrations in our experimental system (concentrations ≥ 5 mM). However, alterations of E. coli gene expression were already detectable at 280 μM as observed by the up-regulation of katG, and as little as 100 nM of H2O2 has been documented to activate the redox-dependent transcription factor OxyR [20]. It has been calculated that H2O2 concentrations generated by neutrophils are in the micromolar order [44]. Thus, it is likely that H2O2 concentrations generated in vivo do not induce direct toxicity, but rather regulate gene expression.
In a previous study, using comparable experimental protocols, we found that Staphylococcus aureus resisted up to 10-fold higher H2O2 concentrations [45] as compared to E. coli. While E. coli has powerful defense mechanisms against H2O2, such as the KatG catalase, the antioxidant defense is even more developed or efficient in some other bacteria. The example of S. aureus is particularly interesting, as despite its high ROS resistance, it is the organism that is the most predominant pathogen in CGD patients, devoid of the oxidant-generating phagocyte NADPH oxidase [8].
The analysis of the time-dependent impact of H2O2 exposure was complicated by the fact that the RNA-seq data showed massive differences in gene expression between 10 min and 60 min even under control conditions (i.e., without the addition of H2O2). The difference in gene expression between 10 min and 60 min suggested a transition toward late exponential phase (Figure S3). Thus, we principally concentrated our analysis on the changes observed at 10 min.
We used pathway enrichment analysis and the ISMARA algorithm for detection of gene regulators. ISMARA analysis fully depends on the quality of genome annotation to establish links between gene expression data and transcription factors. As E. coli is one of the most studied organisms and possesses the most complete genomic annotations, the results of this analysis are highly credible. As validation of the significance of this analysis, the widely characterized H2O2-sensitive transcription factors OxyR and SoxRS were indeed among the most activated transcription factors 10 min after the addition of H2O2.
Our results demonstrate that H2O2 treatment not only induced specific ROS-sensitive mechanisms but also induced responses that have previously been described for other stressors. However, the stress responses may be more activated by the H2O2-induced damage to DNA and other macromolecules than by the direct sensing of H2O2 itself. The stringent response is a well-described mechanism commonly activated by a large variety of environmental stresses, such as heat shock [46], amino acid, fatty acid or iron starvation [47,48]. This response profoundly modifies bacterial metabolism essentially through the inhibition of translation, DNA replication and nucleotide biosynthesis. The stringent response is highly conserved among bacteria and has previously been described to be induced by H2O2 [49]. It leads to a slow-growing, stress-resistant phenotype also implicated in bacterial pathogenicity (reviewed in [47]).
The group of Imlay et al. has extensively studied the response to H2O2 in E. coli and has described the activation of the SOS-response in response to 1 to 3 mM H2O2 [50]. The SOS-response is triggered by DNA damage such as UV-induced DNA lesions [51]. It regulates DNA damage repair but also promotes an elevated mutation rate, generating genetic diversity and adaptation, and is a key factor of antibiotic resistance emergence (reviewed in [52]). As H2O2 exposure can also cause DNA damage, it triggers the SOS-response in E. coli [53] and in other Gram-negative [54,55] as well as Gram-positive bacteria [56]. In our analysis, LexA, the transcription factor that regulates the SOS-response, was among the top 10 identified regulators at 10 min. It is particularly interesting that, as opposed to other stress responses, only genes from the SOS-response were up-regulated at both 10 min and 60 min. This suggests that the concentration of H2O2 used in this study led to DNA damage, which required a sustained response for efficient repair, lasting many minutes after H2O2 was fully degraded.
The alternative subunits of the RNA polymerase RpoS, the general stress response regulator [57], and RpoH that controls the heat shock response [58], have also been described to be activated by oxidative stress in E. coli [59,60] and were indeed among the top 10 regulators implicated after 10 min of H2O2 exposure. Deletion of these genes was previously shown to increase sensitivity to H2O2 in diverse Gram-negative bacteria, suggesting a conserved importance of their transcriptional regulation following H2O2 exposure [61,62,63]. In Gram-positive bacteria, another alternative polymerase subunit, SigB, is involved in H2O2 resistance [64,65,66]. RpoS is also known to be active during the transition between the exponential and stationary phases, or even in the mid-exponential phase [67,68], which can explain why its targeted genes appeared more highly expressed in the untreated condition at 60 min (Figure 4).
Compatible with the activation of general stress responses, we observed a massive and concerted down-regulation of genes encoding ribosomal proteins at 10 min and a transient arrest of bacterial growth for 30 min. Thus, exposure to 2.5 mM H2O2 most likely slowed bacterial protein synthesis. This is at least in part compensated by increased expression of genes relevant for the antioxidant defense. Data from the literature confirm that these genes are also up-regulated at the protein levels [38]. Similar stress responses have been shown to be activated by H2O2 in Salmonella enterica serovar Enteritidis [55].
Other transcriptional factors have previously been associated with the response to oxidative stress; however, the evidence is less obvious. Fur, the transcription factor regulating iron homeostasis, is up-regulated by H2O2 via OxyR activation [69,70]. As Fur acts as a repressor, it is thought to be up-regulated to inhibit excessive iron uptake that would otherwise exacerbate H2O2 damage. However, at 10 min, we observed a massive up-regulation of genes normally repressed by Fur. Other transcriptomic data of H2O2 treated E. coli accessible on the Gene Expression Omnibus database (GSE20305, GSE56133) and a RNA-seq experiment in H2O2 treated Salmonella confirmed this observation [55]. Similarly, in our analysis, the Fur-regulated siderophore biosynthesis pathway was the second most up-regulated pathway 10 min after H2O2 addition. Other studies suggest that a siderophore called enterobactin has a role in protection against H2O2 [71,72,73,74].
FNR, the transcription factor which regulates transition from aerobic to anaerobic growth appeared as one of the most altered transcription factors in our analysis. Although FNR has been suggested to sense and respond to oxidative stress [75], its mode of action has not been extensively studied in this context. In summary, the transcriptomic changes observed in our study align with the available literature and confirmed the implication of transcription factors that are less studied in the context of H2O2 exposure.
In our study, we observed a striking up-regulation of CysB regulated genes 10 min after H2O2 exposure and CysB was among the most significantly implicated regulator at 10 min. Similarly, sulfur metabolism was the most up-regulated metabolic pathway, especially genes from the sulfate assimilation and the cysteine biosynthesis pathways that were strongly up-regulated in a strikingly coordinated manner. The up-regulation of these pathways after H2O2 addition were mentioned in the literature [22,23,76], but not specifically addressed. Navigating through E. coli transcriptomic datasets available on the Gene Expression Omnibus database (GSE135556, GSE56133) indeed indicated an up-regulation of these pathways following H2O2 exposure. Further, the sulfur metabolism pathway is induced by H2O2 in both Gram-negative and Gram-positive bacteria, suggesting a highly conserved function [77,78].
CysB belongs to the family of LysR-type transcriptional regulators, a family of transcription factors characterized by an N-terminal helix-turn-helix (HTH) DNA binding motif and a similar amino-acid sequence [79]. This family, which contains 45 members, including OxyR, are thought to have evolved from a common ancestor due to their high homology and the conservation of their DNA-binding domain. Common features include a similar size, the formation of either homodimers or homotetramers (which is the case for CysB), the presence of the helix-turn-helix motif, and the requirement for a small molecule that acts as co-inducer. CysB is activated by sulfur limitation and the presence of N-acetyl-l-serine and leads to up-regulation of the assimilatory sulfate pathway and promotes the synthesis of the cysteine amino acid [79].
Cysteine is among the least abundant residues in proteins in both eukaryotes and prokaryotes; however, it is often highly conserved in proteins due to its reversible regulatory oxidative characteristics [80]. The biological significance of the up-regulation of cysteine biosynthesis following H2O2 is broad. Cysteine oxidation is a key regulatory element for many cellular processes such as metal binding and formation of di-sulfur bridges. However excessive oxidative stress can lead to uncontrolled and irreversible oxidation of the cysteine residues [81]. In our experimental conditions of excessive oxidative stress, irreversibly oxidized cysteine most likely needed to be renewed. Cysteine is also needed to replace the sulfur atom of H2O2-destroyed iron-sulfur clusters (Fe-S), an essential cofactor of various enzymes [82]. Diverse defense mechanisms against oxidative stress depend on cysteine; these mechanisms include the synthesis of thioredoxin and glutathione that maintain the cell redox-potential [83], and the production of hydrogen sulfide [76] and the L-cystine/L-cysteine shuttle system that protect the periplasm against oxidative damage [84].
The deletion mutant of cysB displayed a reduced fitness; however, its sensitivity to H2O2 was surprisingly not different from that of the WT. As the cysB deleted mutant is auxotrophic for cysteine, the sensitivity was assessed in media that contains amino acids (including cysteine), which may have attenuated their sensitivity to H2O2.
The mechanism leading to the up-regulation of sulfate assimilation and cysteine biosynthesis under oxidative stress is still unknown. Cysteine biosynthesis has not been described to be regulated by the prototypical H2O2 responses induced by OxyR and SoxRS [21]. H2O2-dependent up-regulation of CysB target genes was dose-dependent and the deletion of cysB completely abolished their up-regulation following H2O2 exposure. This confirms the specificity of CysB as a yet undescribed signaling pathway regulating sulfur metabolism under oxidative stress. In terms of mechanism, a cysteine regulator (CymR) of Staphylococcus aureus is known to sense oxidative stress through the oxidation of its unique cysteine residue and to up-regulate the cysteine biosynthesis pathway under H2O2 exposure [78]. Although CysB protein structure is completely different from CymR, CysB also contains a unique cysteine residue (Cys163) in its sequence. Our findings shall stimulate further molecular studies to understand the H2O2-sensing mechanism of CysB and the implication of Cys163 in this response.

5. Conclusions

This study aimed at improving the knowledge of how H2O2 exposure affects gene expression in E. coli. The typical H2O2-sensors (OxyR, SoxRS) were confirmed to regulate genes in response to H2O2 exposure and could be considered as a positive control for our experimental conditions. We observed the implication of other stress regulators, confirming the complex response triggered by exposure to H2O2 on the transcriptome of E. coli. The striking up-regulation of the cysteine biosynthesis pathway through the transcription factor CysB was, for the first time, highlighted as an antioxidant response in E. coli.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/antiox11040655/s1, Figure S1. Heatmaps of genes regulated by the regulators SoxRS and HypT. Figure S2. Expression levels of genes reresentative of the expression pattern represented in Figure 3. Figure S3. Dysregulated pathway comparing the untreated conditions at 60 min vs. 10 min. Figure S4 Impaired fitness of cysB deleted mutant. Table S1: differential gene expression table of control compared to H2O2 10 min after H2O2. Table S2: differential gene expression table of control compared to H2O2 60 min after H2O2. Table S3: differential gene expression table of control 60 min compared to control 10 min.

Author Contributions

Conceptualization, K.-H.K., P.F. and M.R.; methodology, M.R. and V.J.; software, S.L.; formal analysis S.L. and Y.C.; investigation, M.R. and E.-J.B.; data curation, S.L.; writing—original draft preparation, M.R. and V.J.; writing—review and editing, Y.C.; visualization, M.R. and V.J.; project administration, P.F.; funding acquisition, K.-H.K. All authors have read and agreed to the published version of the manuscript.

Funding

This study was granted by Swiss National Science Foundation to Karl-Heinz Krause (Funding number 31003A-179478).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Raw data have been deposited at ENA (https://www.ebi.ac.uk/ena) under the following accession number: PRJEB51098.

Acknowledgments

The authors gratefully acknowledge the technical support from Nadia Gaïa, Melanie Roch, Roberto Sierra, Gaël Panis, and the collaborators of the genomic platform iGE3 of the faculty of Medicine of the University of Geneva, Didier Chollet and Mylène Docquier.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Barber, A.E.; Norton, J.P.; Spivak, A.M.; Mulvey, M.A. Urinary tract infections: Current and emerging management strategies. Clin. Infect. Dis. 2013, 57, 719–724. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Buetti, N.; Atkinson, A.; Marschall, J.; Kronenberg, A. Incidence of bloodstream infections: A nationwide surveillance of acute care hospitals in Switzerland 2008–2014. BMJ Open 2017, 7, e013665. [Google Scholar] [CrossRef] [Green Version]
  3. Federal office of public health FOPH, Communicable Disease Division, Bern Switzerland: How Is Antibiotic Resistance Evolving? Available online: https://www.bag.admin.ch/bag/en/home/krankheiten/infektionskrankheiten-bekaempfen/antibiotikaresistenzen/wie-entwickelt-sich-die-antibiotikaresistenzlage---.html (accessed on 5 November 2021).
  4. Tacconelli, E. Global Priority List of Antibiotic-Resistant Bacteria to Guide Research, Discovery, and Development of New Antibiotics; World Health Organisation: Geneva, Switzerland, 2017; p. 7.
  5. Sies, H.; Belousov, V.V.; Chandel, N.S.; Davies, M.J.; Jones, D.P.; Mann, G.E.; Murphy, M.P.; Yamamoto, M.; Winterbourn, C. Defining roles of specific reactive oxygen species (ROS) in cell biology and physiology. Nat. Rev. Mol. Cell Biol. 2022, 1–17. [Google Scholar] [CrossRef] [PubMed]
  6. Nauseef, W.M. The Phagocyte NOX2 NADPH oxidase in microbial killing and cell signaling. Curr. Opin. Immunol. 2019, 60, 130–140. [Google Scholar] [CrossRef]
  7. Nauseef, W.M. Assembly of the phagocyte NADPH oxidase. Histochem. Cell Biol. 2004, 122, 277–291. [Google Scholar] [CrossRef]
  8. Segal, B.H.; Veys, P.; Malech, H.; Cowan, M.J. Chronic granulomatous disease: Lessons from a rare disorder. Biol. Blood Marrow Transpl. 2011, 17, S123–S131. [Google Scholar] [CrossRef] [Green Version]
  9. Miller, B.M.; Liou, M.J.; Zhang, L.F.; Nguyen, H.; Litvak, Y.; Schorr, E.-M.; Jang, K.K.; Tiffany, C.R.; Butler, B.P.; Bäumler, A.J. Anaerobic respiration of NOX1-derived hydrogen peroxide licenses bacterial growth at the colonic surface. Cell Host Microbe 2020, 28, 789–797.e5. [Google Scholar] [CrossRef]
  10. Alvarez, L.A.; Kovačič, L.; Rodríguez, J.; Gosemann, J.-H.; Kubica, M.; Pircalabioru, G.G.; Friedmacher, F.; Cean, A.; Ghişe, A.; Sărăndan, M.B.; et al. NADPH oxidase-derived H2O2 subverts pathogen signaling by oxidative phosphotyrosine conversion to PB-DOPA. Proc. Natl. Acad. Sci. USA 2016, 113, 10406–10411. [Google Scholar] [CrossRef] [Green Version]
  11. Hudson, P.L.; Hung, K.J.; Bergerat, A.; Mitchell, C. Effect of vaginal Lactobacillus species on Escherichia coli growth. Female Pelvic Med. Reconstr. Surg. 2020, 26, 146–151. [Google Scholar] [CrossRef]
  12. Eschenbach, D.A.; Davick, P.R.; Williams, B.L.; Klebanoff, S.J.; Young-Smith, K.; Critchlow, C.M.; Holmes, K.K. Prevalence of hydrogen peroxide-producing Lactobacillus species in normal women and women with bacterial vaginosis. J. Clin. Microbiol. 1989, 27, 251–256. [Google Scholar] [CrossRef] [Green Version]
  13. Gupta, K.; Stapleton, A.E.; Hooton, T.M.; Roberts, P.L.; Fennell, C.L.; Stamm, W.E. Inverse Association of H2O2-producing lactobacilli and vaginal Escherichia coli colonization in women with recurrent urinary tract infections. J. Infect. Dis. 1998, 178, 446–450. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Hilt, E.E.; McKinley, K.; Pearce, M.M.; Rosenfeld, A.B.; Zilliox, M.J.; Mueller, E.R.; Brubaker, L.; Gai, X.; Wolfe, A.J.; Schreckenberger, P.C. Urine is not sterile: Use of enhanced urine culture techniques to detect resident bacterial flora in the adult female bladder. J. Clin. Microbiol. 2014, 52, 871–876. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Thomas-White, K.J.; Gao, X.; Lin, H.; Fok, C.S.; Ghanayem, K.; Mueller, E.R.; Dong, Q.; Brubaker, L.; Wolfe, A.J. Urinary microbes and post-operative urinary tract infection risk in urogynecologic surgical patients. Int. Urogynecol. J. 2018, 29, 1797–1805. [Google Scholar] [CrossRef]
  16. Price, T.K.; Hilt, E.E.; Thomas-White, K.; Mueller, E.R.; Wolfe, A.J.; Brubaker, L. The urobiome of continent adult women: A cross-sectional study. BJOG 2020, 127, 193–201. [Google Scholar] [CrossRef] [PubMed]
  17. Imlay, J.A. Cellular defenses against superoxide and hydrogen peroxide. Annu. Rev. Biochem. 2008, 77, 755–776. [Google Scholar] [CrossRef] [Green Version]
  18. Winterbourn, C.C.; Kettle, A.J. Redox reactions and microbial killing in the neutrophil phagosome. Antioxid. Redox Signal. 2013, 18, 642–660. [Google Scholar] [CrossRef]
  19. Imlay, J.A. Transcription factors that defend bacteria against reactive oxygen species. Annu. Rev. Microbiol. 2015, 69, 93–108. [Google Scholar] [CrossRef] [Green Version]
  20. Seaver, L.C.; Imlay, J.A. Hydrogen peroxide fluxes and compartmentalization inside growing Escherichia coli. J. Bacteriol. 2001, 183, 7182–7189. [Google Scholar] [CrossRef] [Green Version]
  21. Seo, S.W.; Kim, D.; Szubin, R.; Palsson, B.O. Genome-wide reconstruction of OxyR and SoxRS transcriptional regulatory networks under oxidative stress in Escherichia coli K-12 MG1655. Cell Rep. 2015, 12, 1289–1299. [Google Scholar] [CrossRef] [Green Version]
  22. Zheng, M.; Wang, X.; Templeton, L.J.; Smulski, D.R.; LaRossa, R.A.; Storz, G. DNA microarray-mediated transcriptional profiling of the Escherichia coli response to hydrogen peroxide. J. Bacteriol. 2001, 183, 4562–4570. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Wang, S.; Deng, K.; Zaremba, S.; Deng, X.; Lin, C.; Wang, Q.; Tortorello, M.L.; Zhang, W. Transcriptomic response of Escherichia coli O157:H7 to oxidative stress. Appl. Environ. Microbiol. 2009, 75, 6110–6123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Baba, T.; Ara, T.; Hasegawa, M.; Takai, Y.; Okumura, Y.; Baba, M.; Datsenko, K.A.; Tomita, M.; Wanner, B.L.; Mori, H. Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: The Keio collection. Mol. Syst. Biol. 2006, 2, 50. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-Seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  26. Liao, Y.; Smyth, G.K.; Shi, W. Feature counts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014, 30, 923–930. [Google Scholar] [CrossRef] [Green Version]
  27. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. EdgeR: A bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [Green Version]
  28. Keseler, I.M.; Gama-Castro, S.; Mackie, A.; Billington, R.; Bonavides-Martínez, C.; Caspi, R.; Kothari, A.; Krummenacker, M.; Midford, P.E.; Muñiz-Rascado, L.; et al. The EcoCyc database in 2021. Front. Microbiol. 2021, 12, 711077. [Google Scholar] [CrossRef]
  29. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING V11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef] [Green Version]
  30. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2011. [Google Scholar]
  31. Mootha, V.K.; Lindgren, C.M.; Eriksson, K.-F.; Subramanian, A.; Sihag, S.; Lehar, J.; Puigserver, P.; Carlsson, E.; Ridderstråle, M.; Laurila, E.; et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat. Genet. 2003, 34, 267–273. [Google Scholar] [CrossRef]
  32. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [Green Version]
  33. Balwierz, P.J.; Pachkov, M.; Arnold, P.; Gruber, A.J.; Zavolan, M.; van Nimwegen, E. ISMARA: Automated modeling of genomic signals as a democracy of regulatory motifs. Genome Res. 2014, 24, 869–884. [Google Scholar] [CrossRef] [Green Version]
  34. Vandesompele, J.; De Preter, K.; Pattyn, F.; Poppe, B.; Van Roy, N.; De Paepe, A.; Speleman, F. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3, research0034.1. [Google Scholar] [CrossRef] [Green Version]
  35. Rocha, D.J.P.G.; Castro, T.L.P.; Aguiar, E.R.G.R.; Pacheco, L.G.C. Gene expression analysis in bacteria by RT-QPCR. In Quantitative Real-Time PCR: Methods and Protocols; Methods in Molecular Biology; Biassoni, R., Raso, A., Eds.; Springer: New York, NY, USA, 2020; pp. 119–137. ISBN 978-1-4939-9833-3. [Google Scholar]
  36. Thomason, L.C.; Costantino, N.; Court, D.L. E. coli genome manipulation by P1 transduction. Curr. Protoc. Mol. Biol. 2007, 79, 1.17.1–1.17.8. [Google Scholar] [CrossRef]
  37. EUCAST: Disk Diffusion Methodology. Available online: https://www.eucast.org/ast_of_bacteria/disk_diffusion_methodology/ (accessed on 5 November 2021).
  38. Rodríguez-Rojas, A.; Kim, J.J.; Johnston, P.R.; Makarova, O.; Eravci, M.; Weise, C.; Hengge, R.; Rolff, J. Non-lethal exposure to H2O2 boosts bacterial survival and evolvability against oxidative stress. PLoS Genet. 2020, 16, e1008649. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Hidalgo, E.; Ding, H.; Demple, B. Redox signal transduction: Mutations shifting [2Fe-2S] centers of the SoxR sensor-regulator to the oxidized form. Cell 1997, 88, 121–129. [Google Scholar] [CrossRef] [Green Version]
  40. Gebendorfer, K.M.; Drazic, A.; Le, Y.; Gundlach, J.; Bepperling, A.; Kastenmüller, A.; Ganzinger, K.A.; Braun, N.; Franzmann, T.M.; Winter, J. Identification of a hypochlorite-specific transcription factor from Escherichia coli. J. Biol. Chem. 2012, 287, 6892–6903. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Weber, H.; Polen, T.; Heuveling, J.; Wendisch, V.F.; Hengge, R. Genome-wide analysis of the general stress response network in Escherichia coli: SigmaS-dependent genes, promoters, and sigma factor selectivity. J. Bacteriol. 2005, 187, 1591–1603. [Google Scholar] [CrossRef] [Green Version]
  42. Kredich, N.M. Biosynthesis of cysteine. EcoSal Plus 2008, 3. [Google Scholar] [CrossRef] [PubMed]
  43. Nauseef, W.M. Myeloperoxidase in human neutrophil host defense. Cell Microbiol. 2014, 16, 1146–1155. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Li, X.; Imlay, J.A. Improved measurements of scant hydrogen peroxide enable experiments that define its threshold of toxicity for Escherichia coli. Free Radic. Biol. Med. 2018, 120, 217–227. [Google Scholar] [CrossRef]
  45. Buvelot, H.; Roth, M.; Jaquet, V.; Lozkhin, A.; Renzoni, A.; Bonetti, E.-J.; Gaia, N.; Laumay, F.; Mollin, M.; Stasia, M.-J.; et al. Hydrogen peroxide affects growth of S. Aureus through downregulation of genes involved in pyrimidine biosynthesis. Front. Immunol. 2021, 12, 673985. [Google Scholar] [CrossRef]
  46. Gallant, J.; Palmer, L.; Pao, C.C. Anomalous synthesis of PpGpp in growing cells. Cell 1977, 11, 181–185. [Google Scholar] [CrossRef]
  47. Irving, S.E.; Choudhury, N.R.; Corrigan, R.M. The stringent response and physiological roles of (Pp)PGpp in bacteria. Nat. Rev. Microbiol. 2021, 19, 256–271. [Google Scholar] [CrossRef]
  48. Vinella, D.; Albrecht, C.; Cashel, M.; D’Ari, R. Iron limitation induces SpoT-dependent accumulation of PpGpp in Escherichia coli. Mol. Microbiol. 2005, 56, 958–970. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. VanBogelen, R.A.; Kelley, P.M.; Neidhardt, F.C. differential induction of heat shock, SOS, and oxidation stress regulons and accumulation of nucleotides in Escherichia coli. J. Bacteriol. 1987, 169, 26–32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Imlay, J.A.; Linn, S. Mutagenesis and stress responses induced in Escherichia coli by hydrogen peroxide. J. Bacteriol. 1987, 169, 2967–2976. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Krishna, S.; Maslov, S.; Sneppen, K. UV-induced mutagenesis in Escherichia coli SOS response: A quantitative model. PLoS Comput. Biol. 2007, 3, e41. [Google Scholar] [CrossRef]
  52. Podlesek, Z.; Žgur Bertok, D. The DNA damage inducible SOS Response is a key player in the generation of bacterial persister cells and population wide tolerance. Front. Microbiol. 2020, 11, 1785. [Google Scholar] [CrossRef] [PubMed]
  53. Goerlich, O.; Quillardet, P.; Hofnung, M. Induction of the SOS response by hydrogen peroxide in various Escherichia coli mutants with altered protection against oxidative DNA damage. J. Bacteriol. 1989, 171, 6141–6147. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Chang, W.; Small, D.A.; Toghrol, F.; Bentley, W.E. Microarray analysis of Pseudomonas aeruginosa reveals induction of pyocin genes in response to hydrogen peroxide. BMC Genom. 2005, 6, 115. [Google Scholar] [CrossRef] [Green Version]
  55. Liu, X.; Omar, M.; Abrahante, J.E.; Nagaraja, K.V.; Vidovic, S. Insights into the oxidative stress response of Salmonella Enterica serovar enteritidis revealed by the next generation sequencing approach. Antioxidants 2020, 9, 849. [Google Scholar] [CrossRef] [PubMed]
  56. Chang, W.; Small, D.A.; Toghrol, F.; Bentley, W.E. Global transcriptome analysis of Staphylococcus Aureus response to hydrogen peroxide. J. Bacteriol. 2006, 188, 1648–1659. [Google Scholar] [CrossRef] [Green Version]
  57. Hengge, R. The general stress response in gram-negative bacteria. In Bacterial Stress Responses, 2nd ed.; Hengge, R., Storz, G., Eds.; American Society for Microbiology: Washington, DC, USA, 2010; pp. 251–289. ISBN 978-1-68367-121-3. [Google Scholar]
  58. Arsène, F.; Tomoyasu, T.; Bukau, B. The heat shock response of Escherichia coli. Int. J. Food Microbiol. 2000, 55, 3–9. [Google Scholar] [CrossRef]
  59. Lange, R.; Hengge-Aronis, R. Identification of a central regulator of stationary-phase gene expression in Escherichia coli. Mol. Microbiol. 1991, 5, 49–59. [Google Scholar] [CrossRef]
  60. Müller, A.; Hoffmann, J.H.; Meyer, H.E.; Narberhaus, F.; Jakob, U.; Leichert, L.I. Nonnative disulfide bond formation activates the Σ32-dependent heat shock response in Escherichia coli. J. Bacteriol. 2013, 195, 2807–2816. [Google Scholar] [CrossRef] [Green Version]
  61. Jørgensen, F.; Bally, M.; Chapon-Herve, V.; Michel, G.; Lazdunski, A.; Williams, P.; Stewart, G.S.A.B. RpoS-dependent stress tolerance in Pseudomonas aeruginosa. Microbiology 1999, 145, 835–844. [Google Scholar] [CrossRef] [Green Version]
  62. Kojic, M.; Degrassi, G.; Venturi, V. Cloning and characterisation of the RpoS gene from plant growth-promoting Pseudomonas Putida WCS358: RpoS is not involved in siderophore and homoserine lactone production. Biochim. Biophys. Acta (BBA) Gene Struct. Exp. 1999, 1489, 413–420. [Google Scholar] [CrossRef]
  63. Kogoma, T.; Yura, T. Sensitization of Escherichia coli cells to oxidative stress by deletion of the RpoH gene, which encodes the heat shock sigma factor. J. Bacteriol. 1992, 174, 630–632. [Google Scholar] [CrossRef] [Green Version]
  64. Antelmann, H.; Engelmann, S.; Schmid, R.; Hecker, M. General and oxidative stress responses in Bacillus Subtilis: Cloning, expression, and mutation of the Alkyl hydroperoxide reductase operon. J. Bacteriol. 1996, 178, 6571–6578. [Google Scholar] [CrossRef] [Green Version]
  65. Chan, P.F.; Foster, S.J.; Ingham, E.; Clements, M.O. The Staphylococcus Aureus alternative sigma factor ΣB controls the environmental stress response but not starvation survival or pathogenicity in a mouse abscess model. J. Bacteriol. 1998, 180, 6082–6089. [Google Scholar] [CrossRef]
  66. Kazmierczak, M.J.; Mithoe, S.C.; Boor, K.J.; Wiedmann, M. Listeria monocytogenes ΣB regulates stress response and virulence functions. J. Bacteriol. 2003, 185, 5722–5734. [Google Scholar] [CrossRef] [Green Version]
  67. Hengge-Aronis, R. Survival of hunger and stress: The role of RpoS in early stationary phase gene regulation in E. coli. Cell 1993, 72, 165–168. [Google Scholar] [CrossRef]
  68. Maciąg, A.; Peano, C.; Pietrelli, A.; Egli, T.; De Bellis, G.; Landini, P. In vitro transcription profiling of the ΣS subunit of bacterial RNA polymerase: Re-definition of the ΣS regulon and identification of ΣS-specific promoter sequence elements. Nucleic Acids Res. 2011, 39, 5338–5355. [Google Scholar] [CrossRef] [Green Version]
  69. Zheng, M.; Doan, B.; Schneider, T.D.; Storz, G. OxyR and SoxRS regulation of fur. J. Bacteriol. 1999, 181, 4639–4643. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Varghese, S.; Wu, A.; Park, S.; Imlay, K.R.C.; Imlay, J.A. Submicromolar hydrogen peroxide disrupts the ability of fur protein to control free-iron levels in Escherichia coli. Mol. Microbiol. 2007, 64, 822–830. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Adler, C.; Corbalan, N.S.; Peralta, D.R.; Pomares, M.F.; de Cristóbal, R.E.; Vincent, P.A. The alternative role of enterobactin as an oxidative stress protector allows Escherichia coli colony development. PLoS ONE 2014, 9, e84734. [Google Scholar] [CrossRef] [PubMed]
  72. Peralta, D.R.; Adler, C.; Corbalán, N.S.; Paz García, E.C.; Pomares, M.F.; Vincent, P.A. Enterobactin as part of the oxidative stress response repertoire. PLoS ONE 2016, 11, e0157799. [Google Scholar] [CrossRef] [Green Version]
  73. Yeoh, B.S.; Saha, P.; Xiao, X.; Singh, V.; Vijay-Kumar, M. Enterobactin, a metallophore, mitigates the immune responses of macrophages. J. Immunol. 2017, 198, 121.13. [Google Scholar]
  74. Saha, P.; Yeoh, B.S.; Olvera, R.A.; Xiao, X.; Singh, V.; Awasthi, D.; Subramanian, B.C.; Chen, Q.; Dikshit, M.; Wang, Y.; et al. Bacterial Siderophores hijack neutrophil functions. J. Immunol. 2017, 198, 4293–4303. [Google Scholar] [CrossRef]
  75. Green, J.; Scott, C.; Guest, J.R. Functional versatility in the CRP-FNR superfamily of transcription factors: FNR and FLP. In Advances in Microbial Physiology; Academic Press: Cambridge, MA, USA, 2001; Volume 44, pp. 1–34. [Google Scholar]
  76. Mironov, A.; Seregina, T.; Nagornykh, M.; Luhachack, L.G.; Korolkova, N.; Lopes, L.E.; Kotova, V.; Zavilgelsky, G.; Shakulov, R.; Shatalin, K.; et al. Mechanism of H2S-mediated protection against oxidative stress in Escherichia coli. Proc. Natl. Acad. Sci. USA 2017, 114, 6022–6027. [Google Scholar] [CrossRef] [Green Version]
  77. Turnbull, A.L.; Surette, M.G. Cysteine biosynthesis, oxidative stress and antibiotic resistance in Salmonella typhimurium. Res. Microbiol. 2010, 161, 643–650. [Google Scholar] [CrossRef] [PubMed]
  78. Ji, Q.; Zhang, L.; Sun, F.; Deng, X.; Liang, H.; Bae, T.; He, C. Staphylococcus aureus CymR is a new thiol-based oxidation-sensing regulator of stress resistance and oxidative response. J. Biol. Chem. 2012, 287, 21102–21109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Lochowska, A.; Iwanicka-Nowicka, R.; Plochocka, D.; Hryniewicz, M.M. Functional dissection of the LysR-type CysB transcriptional regulator: Regions important for DNA binding, inducer response, oligomerization, and positive control. J. Biol. Chem. 2001, 276, 2098–2107. [Google Scholar] [CrossRef] [Green Version]
  80. Miseta, A.; Csutora, P. Relationship between the occurrence of cysteine in proteins and the complexity of organisms. Mol. Biol. Evol. 2000, 17, 1232–1239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Chung, H.S.; Wang, S.-B.; Venkatraman, V.; Murray, C.I.; Van Eyk, J.E. Cysteine oxidative post-translational modifications: Emerging regulation in the cardiovascular system. Circ. Res. 2013, 112, 382–392. [Google Scholar] [CrossRef] [Green Version]
  82. Jang, S.; Imlay, J.A. Hydrogen peroxide inactivates the Escherichia coli Isc iron-sulphur assembly system, and OxyR induces the suf system to compensate. Mol. Microbiol. 2010, 78, 1448–1467. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Toledano, M.B.; Kumar, C.; Le Moan, N.; Spector, D.; Tacnet, F. The system biology of thiol redox system in Escherichia coli and yeast: Differential functions in oxidative stress, iron metabolism and DNA synthesis. FEBS Lett. 2007, 581, 3598–3607. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Ohtsu, I.; Wiriyathanawudhiwong, N.; Morigasaki, S.; Nakatani, T.; Kadokura, H.; Takagi, H. The L-Cysteine/l-cystine shuttle system provides reducing equivalents to the periplasm in Escherichia coli. J. Biol. Chem. 2010, 285, 17479–17487. [Google Scholar] [CrossRef] [Green Version]
Antioxidants 11 00655 g001 550
Figure 1. Effect of H2O2 on E. coli BW25113. (A) Bacterial growth of E. coli BW25113 exposed to increasing concentrations of H2O2, added during the exponential phase (OD595nm = 0.2). The color code corresponds to the different H2O2 concentrations depicted in the graph in (B); (means of technical duplicate +/− SD; one representative experiment); (B) Representation of the growth rate (% of the control condition) as a function of H2O2 concentrations. The sublethal concentration was defined as the highest concentration where the growth rate was not significantly different from the control condition, namely 2.5 mM; (mean +/− SD; N = 3; statistical test was ANOVA with Turkey multiple comparison; *** p ≤ 0.001); (C) Quantification of H2O2 degradation using Amplex Red/HRP after addition of 2.5 mM H2O2 to bacterial culture at OD595nm = 0.2 (mean +/− SD; N = 3); (D) Growth rate (% of the control condition) of the katG deleted mutant as a function of the H2O2 concentration, WT from panel B shown as dotted line (mean +/− SD, N = 3).
Figure 1. Effect of H2O2 on E. coli BW25113. (A) Bacterial growth of E. coli BW25113 exposed to increasing concentrations of H2O2, added during the exponential phase (OD595nm = 0.2). The color code corresponds to the different H2O2 concentrations depicted in the graph in (B); (means of technical duplicate +/− SD; one representative experiment); (B) Representation of the growth rate (% of the control condition) as a function of H2O2 concentrations. The sublethal concentration was defined as the highest concentration where the growth rate was not significantly different from the control condition, namely 2.5 mM; (mean +/− SD; N = 3; statistical test was ANOVA with Turkey multiple comparison; *** p ≤ 0.001); (C) Quantification of H2O2 degradation using Amplex Red/HRP after addition of 2.5 mM H2O2 to bacterial culture at OD595nm = 0.2 (mean +/− SD; N = 3); (D) Growth rate (% of the control condition) of the katG deleted mutant as a function of the H2O2 concentration, WT from panel B shown as dotted line (mean +/− SD, N = 3).
Antioxidants 11 00655 g001
Antioxidants 11 00655 g002 550
Figure 2. Gene expression at 10 and 60 min after exposure of E. coli BW25113 to 2.5 mM H2O2. Volcano plots representing the dysregulated genes at 10 (A) and 60 min (B) after H2O2 addition compared to respective untreated cells. Blue dots represent significantly up-regulated genes and red dots represent significantly down-regulated genes; thresholds were set as p-value = 0.01 and fold-change ≥ 5; (C) Multi-Dimensional Scaling plot of all samples, each color represent a condition, each dot a sample. Distances between dots represent the coefficient of variation of expression between samples for the top 500 dysregulated genes that best distinguish the samples; (D) Heatmap of the expression of genes of the OxyR regulon with up-regulated genes in blue and down-regulated genes in red [28] (fold change); (E) Catalase katG mRNA expression levels in response to increasing H2O2 concentrations, as measured by qRT-PCR (mean +/− SEM, N = 3).
Figure 2. Gene expression at 10 and 60 min after exposure of E. coli BW25113 to 2.5 mM H2O2. Volcano plots representing the dysregulated genes at 10 (A) and 60 min (B) after H2O2 addition compared to respective untreated cells. Blue dots represent significantly up-regulated genes and red dots represent significantly down-regulated genes; thresholds were set as p-value = 0.01 and fold-change ≥ 5; (C) Multi-Dimensional Scaling plot of all samples, each color represent a condition, each dot a sample. Distances between dots represent the coefficient of variation of expression between samples for the top 500 dysregulated genes that best distinguish the samples; (D) Heatmap of the expression of genes of the OxyR regulon with up-regulated genes in blue and down-regulated genes in red [28] (fold change); (E) Catalase katG mRNA expression levels in response to increasing H2O2 concentrations, as measured by qRT-PCR (mean +/− SEM, N = 3).
Antioxidants 11 00655 g002
Antioxidants 11 00655 g003 550
Figure 3. Comparison of transcriptomic changes at 10 and 60 min; (A) Each gene is represented by a black dot, x-axis: Fold Change (log2) at 10 min (H2O2 vs. untreated control), y-axis: Fold Change (log2) at 60 min (H2O2 vs. untreated control), pink: up-regulated at 10 min, blue: down-regulated at 10 min, green: up-regulated at 60 min, yellow: down-regulated at 60 min; (B) Venn diagram of conditions presented in A documenting the number of genes in each condition.
Figure 3. Comparison of transcriptomic changes at 10 and 60 min; (A) Each gene is represented by a black dot, x-axis: Fold Change (log2) at 10 min (H2O2 vs. untreated control), y-axis: Fold Change (log2) at 60 min (H2O2 vs. untreated control), pink: up-regulated at 10 min, blue: down-regulated at 10 min, green: up-regulated at 60 min, yellow: down-regulated at 60 min; (B) Venn diagram of conditions presented in A documenting the number of genes in each condition.
Antioxidants 11 00655 g003
Antioxidants 11 00655 g004 550
Figure 4. The top 10 most significant regulators implicated in gene expression changes 10 min after H2O2 exposure, using ISMARA. The significance (z-score) is illustrated for the two time points, 10 min (dark purple) and 60 min (light purple); bars extend to the left if target genes are down-regulated or to the right if they are up-regulated in the H2O2-exposed condition.
Figure 4. The top 10 most significant regulators implicated in gene expression changes 10 min after H2O2 exposure, using ISMARA. The significance (z-score) is illustrated for the two time points, 10 min (dark purple) and 60 min (light purple); bars extend to the left if target genes are down-regulated or to the right if they are up-regulated in the H2O2-exposed condition.
Antioxidants 11 00655 g004
Antioxidants 11 00655 g005 550
Figure 5. Dysregulated metabolic pathways 10 min after H2O2 addition, determined using gene set enrichment analysis. Each dot represents a metabolic pathway from the KEGG database; thresholds: FDR = 0.05, NES (normalized enrichment score) >1 and <−1.
Figure 5. Dysregulated metabolic pathways 10 min after H2O2 addition, determined using gene set enrichment analysis. Each dot represents a metabolic pathway from the KEGG database; thresholds: FDR = 0.05, NES (normalized enrichment score) >1 and <−1.
Antioxidants 11 00655 g005
Antioxidants 11 00655 g006 550
Figure 6. Dysregulation of the sulfur assimilation and cysteine biosynthesis pathways; (A) Heatmap of the differential expression of genes of the CysB regulon (fold-change (log2)); (B) Disk diffusion assay performed with 10 µL of 1 M H2O2, showing one representative experiment of the WT and the cysB and katG deletion mutants; (C) Zone of bacterial growth inhibition assessed by the disk diffusion assay on the WT and the cysB and katG deletion mutants (mean+/− SD; N = 6; statistical test was Kruskal-Wallis with Dunn multiple comparison, ns: p > 0.05, *: p ≤ 0.05, **: p ≤ 0.01).
Figure 6. Dysregulation of the sulfur assimilation and cysteine biosynthesis pathways; (A) Heatmap of the differential expression of genes of the CysB regulon (fold-change (log2)); (B) Disk diffusion assay performed with 10 µL of 1 M H2O2, showing one representative experiment of the WT and the cysB and katG deletion mutants; (C) Zone of bacterial growth inhibition assessed by the disk diffusion assay on the WT and the cysB and katG deletion mutants (mean+/− SD; N = 6; statistical test was Kruskal-Wallis with Dunn multiple comparison, ns: p > 0.05, *: p ≤ 0.05, **: p ≤ 0.01).
Antioxidants 11 00655 g006
Antioxidants 11 00655 g007 550
Figure 7. Expression of katG and sulfur metabolism-related genes (cysB, cysE, cysH, cysI, cysJ, cysN. tcyP) by qRT-PCR; the WT and the cysB deleted strain were exposed to the indicated concentrations of H2O2 (N = 3; mean +/− SEM).
Figure 7. Expression of katG and sulfur metabolism-related genes (cysB, cysE, cysH, cysI, cysJ, cysN. tcyP) by qRT-PCR; the WT and the cysB deleted strain were exposed to the indicated concentrations of H2O2 (N = 3; mean +/− SEM).
Antioxidants 11 00655 g007
Table
Table 1. Strains used in this study.
Table 1. Strains used in this study.
NameGenotypeReference
BW25113F-Δ(araD-araB)567ΔlacZ4787(::rrnB-3)rph-1Δ(rhaDrhaB) 568hsdR514[24]
JW3914BW25113, Δkat::kan[24]
JW1267BW25113, ΔcysB::kan[24]
Table
Table 2. Primers used in this study.
Table 2. Primers used in this study.
NameSequenceEfficiency
(for qPCR Primers)		Reference
secA_qPCR_FGGTAGTCGTAACGATCGCA1.95 (95.06%)[35]
secA_qPCR_RTTTCCATCTCCGGTTCCAT[35]
gyrB_N_qPCR_FGTCCTGAAAGGGCTGGATG1.89 (89.37%)[35]
gyrB_N_qPCR_RCGAATACCATGTGGTGCAGA[35]
gyrB_V_qPCR_FGAAATTCTCCTCCCAGACCA1.83 (82.56%)[35]
gyrB_V_qPCR_RGCAGTTCGTTCATCTGCTGT[35]
katG_qPCR_FGGGCCGACCTGTTTATCCTC1.92 (92.09%)This study
katG_qPCR_RATCCAGATCCGGTTCCCAGAThis study
cysB_qPCR_FACTGTATATCGCCACCACGC2.15 (115.12%)This study
cysB_qPCR_RCAGCAATTTGTGTCGGCGAGThis study
cysE_qPCR_FCCTGCAAGGCATGAAAACCC1.92 (91.55%)This study
cysE_qPCR_RCCTCTGCGGCCTGTGThis study
cysH_qPCR_FCAGCACCGGTAAATTGGCAC1.93 (92.60%)This study
cysH_qPCR_RGCACGCTACGGAAAACTGTGThis study
cysI_qPCR_FCTTGTTGCCGTGTTCGATGG1.93 (93.35%)This study
cysI_qPCR_RCGCCAACGACATGAACTTCGThis study
cysJ_qPCR_FCTCAACGTCTTTCGCCATGC2.03 (103.25%)This study
cysJ_qPCR_RGCTACGTCAAAGATGGCGTGThis study
cysN_qPCR_FGCGACGTTTGATTCCACACC1.95 (94.84%)This study
cysN_qPCR_RCTCGATTTTCGTGGTTACGCCGThis study
tcyP_qPCR_FTCGCCGCACTGATTGTACTG2.05 (104.82%)This study
tcyP_qPCR_RGCCACTAACGTTTAACGCCGThis study
katG_seq_FACCCTTTTTTATAAAGCATTTGTCCG This study
katG_seq_RGGGTTGCTCTTTCCTGCG This study
cysB_seq_3_FATGTTGATGGCAAATGGGTTGAAGG This study
cysB_seq_3_RGCCATCACTTATCAGCAAGACG This study
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Roth, M.; Jaquet, V.; Lemeille, S.; Bonetti, E.-J.; Cambet, Y.; François, P.; Krause, K.-H. Transcriptomic Analysis of E. coli after Exposure to a Sublethal Concentration of Hydrogen Peroxide Revealed a Coordinated Up-Regulation of the Cysteine Biosynthesis Pathway. Antioxidants 2022, 11, 655. https://doi.org/10.3390/antiox11040655

AMA Style

Roth M, Jaquet V, Lemeille S, Bonetti E-J, Cambet Y, François P, Krause K-H. Transcriptomic Analysis of E. coli after Exposure to a Sublethal Concentration of Hydrogen Peroxide Revealed a Coordinated Up-Regulation of the Cysteine Biosynthesis Pathway. Antioxidants. 2022; 11(4):655. https://doi.org/10.3390/antiox11040655

Chicago/Turabian Style

Roth, Myriam, Vincent Jaquet, Sylvain Lemeille, Eve-Julie Bonetti, Yves Cambet, Patrice François, and Karl-Heinz Krause. 2022. "Transcriptomic Analysis of E. coli after Exposure to a Sublethal Concentration of Hydrogen Peroxide Revealed a Coordinated Up-Regulation of the Cysteine Biosynthesis Pathway" Antioxidants 11, no. 4: 655. https://doi.org/10.3390/antiox11040655

APA Style

Roth, M., Jaquet, V., Lemeille, S., Bonetti, E. -J., Cambet, Y., François, P., & Krause, K. -H. (2022). Transcriptomic Analysis of E. coli after Exposure to a Sublethal Concentration of Hydrogen Peroxide Revealed a Coordinated Up-Regulation of the Cysteine Biosynthesis Pathway. Antioxidants, 11(4), 655. https://doi.org/10.3390/antiox11040655

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop