To identify genes with similar transcription patterns during the time course, we classified the 343 differentially regulated genes into seven groups on the basis of their transcription directions. Table shows the genes in each group and their n-fold changes and P values in response to 10- and 20-min exposures. Note that genes belonging to the functional classes of “hypothetical proteins,” “hypothetical proteins-conserved,” and “unknown function” are not included in Table . The complete description of all 343 genes is displayed in Table S1 in the supplemental material. Briefly, group I contained 20 genes upregulated upon both exposure times, while group II had 92 genes with increased expression levels at 10 min and no significant changes upon 20 min of exposure. Further, group III possessed seven genes that were induced and repressed in response to 10- and 20-min exposures, respectively. Group IV contained 132 genes downregulated after 10 min, whereas 12 genes of group V exhibited decreased mRNA levels upon both exposure times. Finally, groups VI and VII had 68 and 12 genes that were induced and repressed, respectively, upon 20 min of exposure.
| TABLE 1.S. aureus genes that showed statistically significant mRNA level changes upon either 10 or 20 min of exposure to hydrogen peroxideaa |
Group I was composed of 20 genes that were induced upon both 10- and 20-min exposures. One of the characteristics of this group was the induction of DNA repair-related genes COL-SA0823 (
uvrB), -SA0824 (
uvrA), -SA1374 (
lexA), -SA1400, and -SA2131 (Table ). Specifically,
uvrA and
uvrB encode proteins required during the early steps of nucleotide excision repair to form a DNA-protein complex at the damaged site that allows incision to occur (
24). The
lexA gene encodes a repressor protein regulating the SOS response genes, which include DNA repair and recombination genes (
2,
3). In relation to this gene, we also observed that COL-SA1304, which codes for RecA, was upregulated upon 20 min of exposure (group VI). This gene was also induced by 2.3-fold in response to 10 min of exposure, but the
P value (0.079) was higher than the cutoff. The RecA protein stimulates the autocatalytic cleavage of LexA and thus increases the expression of the genes of the SOS regulon (
3). Congruent with this finding, increased expression of
recA and
lexA because of hydrogen peroxide treatment was previously reported in
E. coli and
P. aeruginosa (
46,
63). Besides, COL-SA1400 codes for an ImpB/MucB/SamB family protein which is involved in UV protection, whereas COL-SA2131 encodes a Dps family protein, the DNA-binding ferritin-like protein, which plays a central role in protecting DNA from oxidative damage by directly binding to DNA (
61). Consequently, our data corroborate previous studies that have associated oxidative stress response genes with hydrogen peroxide and other reactive oxidants and reinforce the conclusion that DNA repair proteins may be among the most central mechanisms that
S. aureus uses to counteract lethal effects of reactive oxygen intermediates. More importantly, our data suggest that the DNA repair system was continuously activated even after the growth of
S. aureus, which had been initially inhibited by hydrogen peroxide, resumed at the same rate as that of untreated cells.
Group II consisted of 92 genes that were upregulated at 10 min; however, upon 20 min of exposure, the expression level changes of these genes became statistically insignificant. As shown in Table , this group also had a number of genes belonging to the functional class of “DNA metabolism.” Specifically, COL-SA1241 (
recG) is involved in the repair of DNA damage resulting from quinolone treatment in
S. aureus (
43). The Nth protein (endonuclease III) encoded by COL-SA1492 is a DNA glycosylase involved in the first step of base excision repair of DNA damage in
E. coli (
16,
31). Moreover, the DnaD protein, putatively encoded by COL-SA1493, is essential for the initiation step in DNA replication and is also involved in DNA repair (
34). As discussed above, group I also had DNA repair-related genes which exhibited expression level increases upon both 10- and 20-min exposures, whereas the genes here showed increases only at 10 min of exposure. Therefore, this result suggests that DNA repair mechanisms are selectively induced to maintain DNA integrity for the synthesis of proteins vital for cell survival. For example, since the
nth gene is related to oxidative pyrimidine damage (
31), no significant change in the expression level of this gene at 20 min of exposure might indicate that this lesion was already repaired by the base excision pathway at that time. Moreover, this versatile repair capability might be one of the schemes that allow
S. aureus to resume growing even while part of the damage was apparently still being restored.
Of particular interest was the finding that several genes directly associated with the virulence of
S. aureus were included in group II. Specifically, the Cap5C protein encoded by COL-SA0138 is involved in capsular polysaccharides synthesis, which enhances staphylococcal virulence by impeding phagocytosis (
45). COL-SA1522 encodes an elastin binding protein which promotes bacterial adherence to extracellular matrix and thus the colonization of host tissues during infection (
49). Furthermore, COL-SA1179 and -SA1180 reportedly encode exotoxins which are involved in food poisoning and toxic shock syndrome (
22,
59). This finding is congruent with the previous outcome that hydrogen peroxide induces several virulence factor-related genes in
P. aeruginosa (
46). Indeed, virulence-related enzymes are involved in microbial defense systems against oxidants by damaging phagocytes and/or impairing oxidants (
40). Further, it was previously revealed that virulence factors scavenge reactive oxygen species (
7,
8,
32). Consequently, the possibility that the virulence-related genes mentioned above, such as
cap5C and COL-SA1178 to -SA1180, contributed to cellular oxidative defense against hydrogen peroxide in
S. aureus should not be excluded.
Group IV comprised 132 gene that exhibited mRNA level decreases at 10 min of exposure and no significant changes at 20 min of exposure (Table ). This group represents the largest portion of the statistically significant 343 genes in our study. One of the most distinctive features of group IV was that 40 out of the 132 genes belonged to the functional class of “transport and binding proteins.” These genes are primarily related to permeases and ATP-binding cassette (ABC) proteins. As presented in Fig. , this functional class consists of a total of 288 genes in S. aureus; thus, 14% of its genes were repressed at10 min in response to hydrogen peroxide exposure. This finding possibly implies that membrane components of S. aureuswere altered and that active and/or facilitated transport through the cell membrane was initially attenuated upon exposure to hydrogen peroxide. Furthermore, the finding that most of these genes exhibited no expression level changes at 20 min of exposure indicates that the transport system of S. aureus was restored, which might be linked to the resumption of growth.
Particularly important was the finding that many of the genes in the class of “transport and binding proteins” were also members of the “signal transduction” class (Table ). Intriguingly, the genes were all involved in the bacterial phosphoenolpyruvate:sugar phosphotransferase system, which mediates the uptake and phosphorylation of carbohydrates and controls metabolism in response to their availability (
25). The system is composed of several types of proteins; however, the genes here exclusively encode the carbohydrate-specific enzymes IIA and IIB and/or the membrane permease IIC, which recognizes and transports the sugar molecules (
25). Besides, we found that COL-SA0403 and -SA2147, genes adjacent to COL-SA0402 and -SA2146, code for enzymes belonging to the BglG family of transcriptional antiterminators that regulate the expression of bacterial genes and operons, whose products are required for utilization of phosphoenolpyruvate:sugar phosphotransferase system carbohydrates (
21). Hence, the repression of these genes can directly deteriorate carbohydrate uptake and subsequent metabolism in
S. aureus, which might be associated with the growth arrest effect of hydrogen peroxide.
To our surprise, the repression of the intercellular adhesion locus (
icaADBC), which is associated with the virulence activity of
S. aureus, was found in group IV. The
icaADBC locus mediates polysaccharide intercellular adhesion (PIA) production in
S. aureus and
Staphylococcus epidermidis, which leads to cell-cell adhesion and is required for biofilm formation (
13,
41,
56). PIA is synthesized by the expression of the
icaADBC genes, which encode three membrane proteins (IcaA, IcaD, and IcaC) with enzymatic activity and one extracellular protein (IcaB) (
17). Prior studies demonstrated that PIA production is involved in the pathogenesis of
S. epidermidis (
50,
51) and is also induced by subinhibitory concentrations of certain antibiotics (
48). Considering these previous conclusions, it was striking that hydrogen peroxide-driven oxidative stress repressed the transcription of the
icaADBC locus in our study. However, our result may propose the possibility that hydrogen peroxide insult attenuated biofilm formation, which depends on the activity of the
icaADBC locus.
Table also shows that many genes related to primary metabolic pathways, including the classes of “energy metabolism” and “fatty acid and phospholipid metabolism,” were repressed in response to hydrogen peroxide. This phenomenon might reflect general changes in cellular physiology and a metabolic repression as a result of oxidative damage. Related to this finding, we also discovered that a number of genes in the class of “biosynthesis of cofactors, prosthetic groups, and carriers” were repressed. All these genes are responsible for the synthesis of various cofactors such as nicotinamide, pantothenate, riboflavin, and thiamine, which are essential for many enzymatic reactions in respiration. Therefore, it can be postulated that the repression of these genes interfered with part of respiratory metabolic pathways, which may be associated with the growth inhibition seen upon 10 min of exposure.
Group VI consisted of 68 genes whose expression levels increased only in response to 20 min of exposure (Table ). The most dominant class was “transport and binding proteins,” which possessed one-fourth of the genes in group VI. Moreover, many of the genes in this class encode proteins that convey cations and iron-carrying compounds. In particular, COL-SA0099 (
sirA), -SA0104, and -SA0105 aroused our attention because they are involved in iron uptake system in
S. aureus. First, COL-SA0104 and -SA0105 code for siderophore (iron-chelating compound) biosynthesis proteins. Siderophore-mediated iron uptake is one of the most important mechanisms that bacteria use to acquire iron from the environment (
57). Second, the SirA protein encoded by COL-SA0099 is involved in iron-siderophore import in
S. aureus (
14). Iron metabolism is coordinately regulated with oxidative stress defenses because iron promotes the formation of hydroxyl radicals, which indiscriminately damage all cellular components (
62). Further, superoxide, generated during the process of oxygen reduction, releases free iron from iron-sulfur proteins, thus increasing the levels of intracellular free iron (
30). Supporting this hypothesis, COL-SA0665 and -SA0666 of group IV, iron compound transport proteins, were repressed upon 10 min of exposure to hydrogen peroxide. Consequently, our results may indicate that the iron uptake system was attenuated to prevent further oxidative damage and/or was initially inactivated by the increased concentration of intracellular free iron resulting from the oxidative damage, but the uptake resumed subsequently for cellular growth after the normal iron level was restored.
In relation to iron metabolism, we also found genes responsible for heme synthesis and iron storage in group VI. COL-SA1887 (
hemG) and -SA1888 (
hemH) code for proteins involved in heme synthesis; in particular, HemH catalyzes the final step of heme biosynthesis, which involves the insertion of ferrous iron into protoporphyrin IX (
19). Heme is essential for respiration and defense against oxygen intermediates because heme compounds are cofactors for cytochromes and catalases (
19). Consistent with our finding, it was previously demonstrated that control of heme biosynthesis is attuned more to oxidative stress than to iron levels and HemH is induced by hydrogen peroxide-driven oxidative stress in
Salmonella enterica (
19). Next, the protein encoded by COL-SA0799 and -SA1952 was homologous to a ferritin family protein. This outcome is intriguing because ferritins are the major iron storage proteins that contribute to scavenging intracellular iron, which lessens cellular oxidative damage, as discussed above (
42). Hence, this result suggests that the induction of these genes possibly helped protect against hydrogen peroxide-caused oxidative stress by controlling intracellular iron levels. However, the reason that these genes were upregulated only upon 20 min of exposure awaits further investigation.
In Table , another notable finding was the presence of COL-SA0245 (
lytS) and -SA0246 (
lytR) in group VI. The
lytS and
lytR genes, whose products are members of the two-component regulatory family of proteins, are involved in the control of autolysis by affecting murein hydrolase activity, which is important in the biological processes of antibiotic resistance, cell division, cell-to-surface adhesion, and biofilm formation (
4,
20,
27,
39). Specifically, a
lytS mutant strain exhibits an increased propensity for spontaneous lysis and an increased rate of penicillin- and Triton X-100-induced lysis in
S. aureus (
4,
27), whereas a
lytR mutant strain shows defective cell division and attenuated autolytic activity (
11). In conjunction with this finding, we also observed that the
scdA gene of group I, which is immediately upstream of the
lytSR genes and important for staphylococcal cell division (
5), showed much stronger expression (10-fold) upon 20 min of exposure. Consequently, our data suggest that the induction of
lytSR might be involved in the regulation of cell division, which apparently occurred more vigorously upon 20 min of exposure in our study.
Table also shows that several genes related to pathogenesis of
S. aureus were present in group VI. For example, COL-SA1168 (
efb) encodes a virulence factor that binds to both the complement C3b and fibrinogen, inhibits complement activation, and blocks opsonophagocytosis (
33). Further, COL-SA1173 (
hly) and -SA2006 are likely associated with alpha-toxins, which cause membrane damage to many types of mammalian cells (
6). It should be emphasized that many of the virulence-related genes of
S. aureus were differently regulated in response to oxidative stress; that is, such virulence-related genes as
cap5C, COL-SA1179, and -SA1180 of group II were induced only at 10 min, whereas others including
efb,
hly, and COL-SA2006 of group VI were induced only at 20 min. On the other hand, the
icaADBC genes were downregulated upon exposure to hydrogen peroxide.
The last and perhaps most striking result is that a number of genes of group VI encode proteins involved in anaerobic metabolism, most of which belonged to the functional class of “energy metabolism.” Table shows that COL-SA0204 (pflB), -SA0205 (pflA), -SA1094 (cydA), -SA1095 (cydB), -SA1181 (arcB), -SA1182 (arbC), -SA2618 (ldh), -SA2634 (nrdG), and -SA2635 (nrdD) were classified into that category. Note that the expression level changes of all these genes were also validated by using quantitative real-time PCR analysis. Table shows that our microarray results were corroborated by real-time PCR analysis, which provides independent verification of transcript level changes of the genes discussed here.
First, the
pflBA genes, which exhibited 25- and 33-fold increases in transcription levels upon 20 min, respectively, code for enzymes homologous to pyruvate formate-lyases that catalyze the nonoxidative dissimilation of pyruvate to acetyl coenzyme A and formate when
E. coli grows under oxygen-limiting conditions (
54). The proteins encoded by
arcBC are responsible for the arginine deiminase pathway, which enables arginine-dependent anaerobic growth (
36). Further, the
ldh gene codes for a protein that shares considerable homology to a lactate dehydrogenase that converts pyruvate to lactate in
E. coli under anaerobic conditions (
29). The
nrdDG-encoded enzymes are the class III ribonucleotide reductases that are responsible for the synthesis of deoxyribonucleotides needed for DNA synthesis under oxygen-limiting conditions (
38).
Moreover, the
cydAB genes, which together encode cytochrome
d oxidase, were strongly induced by eight- and ninefold, respectively, upon 20 min of exposure. Cytochrome
d oxidase catalyzes the last step of oxygen respiration and prevails under oxygen-limiting conditions (
26). Interestingly, it was speculated that cytochrome
d oxidase is required under conditions of environmental stress and may have crucial roles in cellular physiology other than acting as an oxidase (
12). Further, prior studies revealed that cytochrome
d oxidase plays an imperative part in cellular protection against oxidative stress, at least under microaerobic growth conditions, by showing that mutation or deletion of the genes encoding the enzyme increases sensitivity to oxidative stress in
E. coli and
Azotobacter vinelandii (
18,
23,
35). It was also suggested that the ability of cytochrome
d oxidase to reduce dioxygen to water might minimize the generation of reactive oxygen species (
15). Therefore, the result that the
cydAB genes were strongly induced upon exposure to hydrogen peroxide strengthens the confidence of the prior assignments about the role of cytochrome
d oxidase in oxidative protection processes.
As mentioned above, in addition to this oxidative protective role of cytochrome
d oxidase, it is also known to be associated with microaerobic dioxygen respiration (
15). That is, the transcription of the
cydAB genes is activated when oxygen becomes limiting (
26). Further, Alexeeva et al. proposed that the rapid consumption of oxygen by cytochrome
d oxidase may contribute to the activity of pyruvate formate-lyase under microaerobic conditions by demonstrating that increased expression of genes coding for cytochrome
d oxidase and pyruvate formate-lyase is coordinated in
E. coli (
1). This hypothesis might account for the phenomenon in our study that the pyruvate formate-lyase genes (
pflBA) were induced upon 20 min of exposure, in chorus with the strong expression of the cytochrome
d oxidase genes (
cydAB).
Consequently, our result described here possibly implies that
S. aureus experienced an oxygen-limiting state in response to hydrogen peroxide-driven oxidative stress. Supporting this possibility is the finding that genes responsible for fermentative metabolism (
pflBA,
arcBC,
ldh, and
nrdGD), as well as genes encoding cytochrome
d oxidase (
cydAB), were upregulated upon 20 min of exposure. Further, our observation with respect to the transcription level changes of the
pflBA and
ldh genes is congruent with the previous result that in
E. coli pfl and
ldhA are induced by more than 10-fold and 2-fold, respectively, by shifting the culture condition from an aerobic to a microaerobic state (
47). Therefore, it seems that
S. aureus underwent similar conditions upon exposure to hydrogen peroxide in our study. Notably, the finding that these genes were significantly induced only upon 20 min of exposure suggests that fermentative or microaerobic respiration, which had not been initially activated, was stimulated afterward in response to oxidative stress. Moreover, despite the activation of fermentative metabolism, which provides significantly less energy,
S. aureus was able to resume growing at the same rate as untreated controls (Fig. ). Considering the fact that our cultures were provided with sufficient aeration for growth, the reason for this phenomenon is obscure. However, the possibility that the cells might strive to avoid further cytotoxicity arising from reactive oxidants produced during normal oxygen respiration should not be excluded. Indeed, this speculation is in line with the outcome of a prior study by Sabra et al. which demonstrated that
Pseudomonas aeruginosa prefers microaerobic conditions for growth and for the formation of some of its virulence factors under oxidative stress (
52). Most strikingly, a similar phenomenon was also observed in mammalian cells. That is, several species of parasites, such as
Schistosoma mansoni,
Angiostrongylus cantonensis, and
Dirofilaria immitis, show a reduction of their aerobic respiration along their developmental cycles on vertebrate blood, relying on fermentation to achieve their energy requirements (
44). Further, the study proposed that the arrest of respiration constitutes an adaptation to avoid the toxic effects of reactive oxygen species (
44).