Quorum-sensing agr system of Staphylococcus aureus primes gene expression for protection from lethal oxidative stress
Abstract
The agr quorum-sensing system links Staphylococcus aureus metabolism to virulence, in part by increasing bacterial survival during exposure to lethal concentrations of H2O2, a crucial host defense against S. aureus. We now report that protection by agr surprisingly extends beyond post-exponential growth to the exit from stationary phase when the agr system is no longer turned on. Thus, agr can be considered a constitutive protective factor. Deletion of agr resulted in decreased ATP levels and growth, despite increased rates of respiration or fermentation at appropriate oxygen tensions, suggesting that Δagr cells undergo a shift towards a hyperactive metabolic state in response to diminished metabolic efficiency. As expected from increased respiratory gene expression, reactive oxygen species (ROS) accumulated more in the agr mutant than in wild-type cells, thereby explaining elevated susceptibility of Δagr strains to lethal H2O2 doses. Increased survival of wild-type agr cells during H2O2 exposure required sodA, which detoxifies superoxide. Additionally, pretreatment of S. aureus with respiration-reducing menadione protected Δagr cells from killing by H2O2. Thus, genetic deletion and pharmacologic experiments indicate that agr helps control endogenous ROS, thereby providing resilience against exogenous ROS. The long-lived ‘memory’ of agr-mediated protection, which is uncoupled from agr activation kinetics, increased hematogenous dissemination to certain tissues during sepsis in ROS-producing, wild-type mice but not ROS-deficient (Cybb−/−) mice. These results demonstrate the importance of protection that anticipates impending ROS-mediated immune attack. The ubiquity of quorum sensing suggests that it protects many bacterial species from oxidative damage.
eLife assessment
This important study outlines how the agr quorum sensing system in Staphylococcus aureus confers long-lived protection against oxidative stress, thereby linking bacterial metabolism to virulence in this pathogen. While the findings, which are supported by solid data, seem at first glance to contradict earlier findings that show increased fitness of agr mutants under oxidative stress, the core conclusions of the study are well-substantiated. The topic of the paper holds broad relevance to microbiologists, especially those focusing on host-pathogen interactions and bacterial responses to ROS.
https://doi.org/10.7554/eLife.89098.4.sa0Important: Findings that have theoretical or practical implications beyond a single subfield
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Solid: Methods, data and analyses broadly support the claims with only minor weaknesses
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Introduction
Innate, bactericidal immune defenses and antimicrobials act, at least in part, by stimulating the accumulation of ROS in bacteria (Spaan et al., 2013; Drlica and Zhao, 2021). Thus, understanding how Staphylococcus aureus and other bacterial pathogens manage ROS-mediated stress has important implications for controlling infections.
Knowledge of factors that govern the biology of ROS has advanced considerably in recent years. For example, studies have centered on how specific metabolic features, such as aerobic respiration, affect killing by ROS (Kohanski et al., 2007; Lobritz et al., 2015), and small-molecule enhancers of ROS-mediated lethality are emerging (Shatalin et al., 2021; Shee et al., 2022). Less well characterized is how defense against ROS and metabolism changes integrate with the virulence regulatory network that promotes S. aureus pathogenesis. The agr quorum-sensing system provides a way to study this dynamic: agr is a major virulence regulator that responds to oxidative stress (H2O2). The response occurs through a redox sensor in AgrA that attenuates agr activity, thereby increasing the expression of glutathione peroxidase (BsaA), an enzyme that detoxifies ROS (Sun et al., 2012). Whether protection from ROS also occurs from positive agr action is unknown and likely to be an important issue in the development of Agr-targeted therapies (Khan et al., 2015).
In cultured S. aureus, agr governs the expression of ~200 genes. Its two-part regulatory role is characterized by (1) increased post-exponential-phase production of toxins and exoenzymes that facilitate dissemination of bacteria via tissue invasion, and (2) decreased production of cell surface and other proteins that facilitate adherence, attachment, biofilm production, and evasion of host defenses (Novick, 2003; Novick and Geisinger, 2008). Thus, agr coordinates a switch from an adherent state to an invasive state at elevated bacterial population density. The invasive state would be facilitated by protection from host defense.
The agr locus consists of two divergent transcription units driven by promoters P2 and P3 (Novick et al., 1995). The P2 operon encodes the quorum-signaling module, which contains four genes, agrB, agrD, agrC, and agrA. AgrC is a receptor histidine kinase, and AgrA is a DNA-binding response regulator. AgrD is an autoinducing, secreted peptide derived from a pro-peptide processed by AgrB. The autoinducing peptide binds to and causes autophosphorylation of the AgrC histidine kinase, which phosphorylates and activates the DNA-binding AgrA response regulator. AgrA then stimulates transcription from the P2 (RNAII) and P3 (RNAIII) promoters. RNAIII is a regulatory RNA that additionally contains the gene for delta-hemolysin (hld). The DNA-binding domain of AgrA contains an intramolecular disulfide switch (Sun et al., 2012). Oxidation leads to dissociation of AgrA from DNA, thereby preventing an AgrA-mediated down-regulation of the BsaA peroxidase.
When we used antimicrobials to study bacterial responses to lethal stress involving the accumulation of ROS, we found that inactivation (deletion) of agr reduces lethality arising from treatment with antimicrobials, such as fluoroquinolones, in a largely bsaA-dependent manner (Kumar et al., 2017). Thus, oxidation sensing appears to be an intrinsic checkpoint that ameliorates the endogenous oxidative burden generated by certain antimicrobials. Surprisingly, deletion of agr increases the lethal effects of exogenous H2O2 (Kumar et al., 2017), in contrast to the expected expression of the protective bsaA system (Sun et al., 2012). Thus, agr must help protect S. aureus from exogenous ROS, a principal host defense, through mechanisms other than bsaA.
In the present work we found that protection by wild-type agr against lethal concentrations of H2O2 was unexpectedly long-lived and (1) associated with decreased expression of respiration genes, and (2) potentially aided by defense systems that suppress the oxidative surge triggered by subsequent, high-level H2O2 exposure. The redox switch in AgrA, plus these additional protective properties, indicate that agr increases resilience to oxidative stress in S. aureus both when it is present and when it is absent. Thus, agr integrates protection from host defense into the regulation of staphylococcal virulence.
Results
agr protects S. aureus from lethal concentrations of H2O2 throughout the growth cycle
Because agr is a quorum-sensing regulon, maximal agr activity occurs during exponential growth (Figure 1—figure supplement 1) and is followed by a sharp drop during stationary phase (Kumar et al., 2017; Geisinger et al., 2012). Surprisingly, protection from H2O2 toxicity by wild-type agr, assessed by comparison with an agr deletion mutant, was observed throughout the growth cycle (Figure 1A). Indeed, maximal protection occurred shortly after overnight growth, long after induction and expression of agr transcripts. Comparison of survival rates of Δagr mutant and wild-type cells, following dilution of overnight cultures and regrowth for 1 hr prior to challenge with 20 mM H2O2, revealed an initial rate of killing that was ~1000 fold faster for the Δagr mutant (Figure 1B). Peroxide concentration dependence was observed up to 10 mM during a 60 min treatment; at that point, mutant survival was about 100-fold lower (Figure 1C). Complementation tests confirmed that the agr deletion elevated killing by H2O2 (Figure 1D).
We also monitored the time required for the wild-type agr survival advantage against H2O2 to manifest itself (Figure 1—figure supplement 2). Overnight cultures were not readily killed by H2O2, as expected from previous results with other lethal stressors (Conlon et al., 2016). Following dilution to fresh medium, wild-type survival dropped gradually, while mutant survival, although lower, was constant for 20 min. By 40 min, mutant survival exhibited a precipitous 10-fold drop not seen with wild-type cells (Figure 1—figure supplement 2). This drop in mutant survival correlated temporally with changes in cell density (Figure 1—figure supplement 2); i.e., the first cell division following dilution to fresh medium. Overall, the agr-mediated survival advantage during H2O2 exposure was absent in stationary-phase cells and small during lag phase (before exponential growth resumes), but it increased markedly during early growth.
Lag-time differences between strains were more obvious in experiments using less complex, chemically defined medium (CDM) with highly diluted starting cultures and automated growth analysis (Figure 1—figure supplement 3). In CDM, wild-type cells divided within ~150 min, while the lag times with the Δagr mutant were more than 205 min (in Tryptic Soy Broth the lag time is 30 min for both). These observations suggest a novel agr-mediated decrease in time to enter exponential growth following dilution of stationary phase cultures. The poor killing of agr mutant cells by H2O2 early in lag phase is consistent with other work in which cells experiencing long lag times are less readily killed (Fridman et al., 2014), presumably due to remaining longer in a dormant, protected state. To focus on effects during growth, subsequent experiments were performed after incubation of overnight cultures for 1 hr in fresh Tryptic Soy Broth unless otherwise specified.
The elevating effect of agr inactivation on H2O2-mediated lethality was observed across a variety of S. aureus strains, although differences in wild-type survival were observed (Figure 1—figure supplement 4). Thus, agr-mediated protection from H2O2 appears to be common among S. aureus lineages.
Expression of RNAIII and repression of Rot is required for protection from H2O2-mediated lethality
ΔrnaIII and Δagr mutants showed identical loss of protection from H2O2-mediated killing (Figure 2A), indicating that protection is RNAIII-dependent. Since RNAIII represses translation of the downstream regulator Rot (Geisinger et al., 2006), a transcription factor having a key role in agr regulation of staphylococcal virulence, we also examined the effects of rot on the protective action of agr against H2O2. When the wild-type strain, a Δagr mutant, a Δrot mutant, and a Δagr Δrot double mutant were compared for survival following treatment with 20 mM H2O2, survival of the Δagr Δrot double mutant phenocopied that of the wild-type strain (Figure 2B): the rot deletion reversed the effect of an agr deficiency. These data are consistent with agr activity allowing induction of rot-repressed genes important for protection from peroxide (RNAIII repression of the Rot repressor).
When a low-copy-number plasmid expressing rot was introduced into a wild-type strain, the transformant was more readily killed by H2O2, indicating that the expression of rot is sufficient for increased lethality (Figure 2C–D). These data suggest that wild-type Rot down-regulates expression of protective genes. The observed epistatic effect of agr and rot did not apply to other downstream, potentially epistatic regulators, such as saeRS, mgrA, and sigB (Figure 2—figure supplement 1; Bronesky et al., 2016). Thus, the epistatic relationship between agr and protection from H2O2 appears to be rot-specific.
Agr-mediated protection from H2O2 stress is kinetically uncoupled from agr activation
Since agr-mediated protection from H2O2 occurs throughout the growth cycle, it was possible that protection arises from constitutive, low-level agr expression rather than from autoinduction and thereby quorum sensing. To test for a requirement of quorum in the agr-mediated oxidative-stress phenotypes, we characterized the role of agr activation using a mixed culture strategy in which one strain, an in-frame deletion mutant of agrBD, is activated in trans by AIP produced by a second, ΔrnaIII mutant strain (Figure 3A). The AIP-responsive ΔagrBD strain carried an intact RNAIII, while the ΔrnaIII mutant was wild-type for agrBD. As shown in Figure 3B, hemolytic activity (a marker for RNAIII) of the ΔagrBD mutant was restored by mixing it with the ΔrnaIII mutant strain that secreted AIP into the surrounding medium. This result confirmed that agrCA-directed trans-activation of RNAIII by AIP remained intact in the ΔagrBD mutant.
Mixed culture tests using these mutants, scored by differential plating for the presence of an erythromycin resistance marker in the ΔagrBD mutant, showed no protection from lethality of H2O2 when the two strains were mixed 1:1 immediately prior to growth from stationary phase (Figure 3C). Autoinducer accumulated during subsequent growth, activating agr expression and commencing protection from exogenous H2O2 (Figure 3D–E). During H2O2 treatment, the percentage of the ΔagrBD mutant (rnaIII+) increased while the percentage of the ΔrnaIII mutant decreased; this cis-acting result is consistent with the idea that pathways downstream from RNAIII, such as those regulated by rot, are the primary drivers of agr-mediated protection from H2O2. These results confirm an intimate link between agr-mediated protection and the quorum-controlled agr gene expression program of late exponential phase. However, after an overnight co-culture of the ΔrnaIII and ΔagrBD mutant strains, the ΔagrBD mutant demonstrated the same degree of protection expected for wild-type cells during exposure to H2O2 (Figure 3F–H). Thus, protection by agr after overnight co-culture extends to growth resumption from stationary phase, prior to reaching quorum, and therefore protection is uncoupled from the constraint of strict cell-density dependence. These results indicate that protection lasts long after maximal transcription of agr, when agr expression has largely halted (Kumar et al., 2017; Geisinger et al., 2012). This phenomenon is a critical feature of the agr system not appreciated in previous analyses of agr activation kinetics.
agr deficiency increases transcription of genes involved in respiration and overflow metabolism in the absence of stress
To explore mechanisms underlying protection from H2O2, we performed RNA-seq with the Δagr and wild-type strains after growth to late exponential growth phase, a point when agr expression is maximal. As expected, agr up-regulated the transcription of many known virulence genes (Supplementary file 1). The Δagr strain showed elevated expression of genes involved in respiration (cydA, qoxA-D) and fermentation (Fuchs et al., 2007; Pagels et al., 2010), including nrdGD, alcohol dehydrogenases (adhE and adh1), and lactate dehydrogenases (ldh, ddh) (Figure 4A and Supplementary file 1). Increased respiration and fermentation are expected to increase energy generation. However, metabolic modeling of transcriptomic data showed a ~30% reduction in tricarboxylic acid (TCA) cycle and lactate flux per unit of glucose taken up by the Δagr mutant (Figure 4B, Supplementary file 1). Additionally, intracellular ATP levels were ~50% lower in the Δagr mutant compared to the wild-type control, suggesting reduced metabolic efficiency during exponential growth (Figure 5A). Moreover, although the agr deletion has little effect on growth in the rich medium in which RNA-seq was performed (Somerville et al., 2002a), analysis in nutrient-constrained medium (CDM) revealed decreased growth rate and yield of the Δagr mutant relative to wild-type S. aureus (Figure 1—figure supplement 3). Collectively, these data suggest that Δagr increases respiration and fermentation to compensate for low metabolic efficiency. Consistent with this idea, agr deficiency also increases ATP-yielding carbon ‘overflow’ pathways, as evidenced by increased acetate production (Figure 5B; Sadykov et al., 2013; Somerville et al., 2002b). The increase in accumulated acetate in the culture medium during exponential growth was largely consumed after 24 hr of growth (Figure 5B). Thus, Δagr mutants exhibit TCA cycle proficiency (Somerville et al., 2002a) and, despite some expense of efficiency, an increased catabolism of acetate.
Differential transcription of selected genes was confirmed by RT-qPCR measurements (Figure 4—figure supplement 1) We also confirmed that respiration levels were lower (15%) in wild-type compared to Δagr (Figure 5E, F). Although the stimulatory effect of the agr deletion on production of the fermentation product lactate was not observed in optimally aerated broth cultures after growth to late exponential growth phase, it was confirmed for organisms grown in broth under more metabolically demanding suboptimal aeration conditions (limitations in the rate of respiration when oxygen is limiting are expected to increase overall levels of fermentation) (Figure 5C). Overall, these results are consistent with transcription-level up-regulation of respiratory and fermentative pathways in agr-deficient strains.
Since respiration and fermentation generally increase NAD+/NADH ratios and since these activities are increased in Δagr strains (Figure 5C and E–F), we expected a higher NAD+/NADH ratio relative to wild-type cells. However, we observed a decrease in the NAD+/NADH ratio due to an increase in NADH accompanied by relative stability in NAD+ compared to wild-type. Collectively, these observations suggest that a surge in NADH accumulation and reductive stress in the Δagr strain induces a burst in respiration, but levels of NADH are saturating, thereby driving fermentation under microaerobic conditions.
To help determine the metabolic fate of glucose, we measured glucose consumption and intracellular levels of pyruvate and TCA-cycle metabolites fumarate and citrate in the wild-type and Δagr mutant strains. At 4 hr of growth to late-exponential phase, intracellular pyruvate, and acetyl-CoA levels were increased in the Δagr mutant compared to wild-type strain, but levels of fumarate and citrate were similar (Figure 5—figure supplement 1D–E). Glucose was depleted after 4 hr of growth, but glucose consumption after 3 hr of growth (exponential phase) was increased in the Δagr mutant compared to the wild-type strain (Figure 5—figure supplement 1A). These observations, together with the decrease in the NAD+/NADH ratio and increase in acetate and lactate production described above, are consistent with a model in which respiration in Δagr mutants is inadequate for (1) energy production, resulting in an increase in acetogenesis, and (2) maintenance of redox balance, resulting in an increase in fermentative metabolism, lactate production, and conversion of NADH to NAD+. Increased levels of acetate compared to lactate under optimal aeration conditions suggests that demand for ATP is in excess of demand for NAD+.
Elevated respiratory activity of Δagr is expected to increase endogenous ROS (Lobritz et al., 2015). To test this idea, we assessed ROS accumulation in bulk culture by flow cytometry of Δagr and wild-type stains using carboxy-H2DCFDA, a dye that becomes fluorescent in the presence of several forms of ROS. As shown in Figure 6, ROS levels increased with agr deficiency, indicating correlation between agr activity, lower ROS levels, and increased bacterial survival in response to exogenous H2O2. These data help explain the elevated lethality of peroxide in the absence of agr. Since lower ROS accumulation in wild-type cells correlates with decreased respiration and protection from killing by H2O2, the data also support the idea that suppression of endogenous ROS is key to agr-mediated protection from exogenous H2O2-mediated lethality.
Transcriptional changes due to Δagr mutation are long-lived and result in down-regulation of H2O2-stimulated genes relative to those in an agr wild-type
We reasoned that the transcriptional changes due to the Δagr mutation likely persist, as does this strain’s susceptibility to killing by H2O2, after growth from overnight culture. With this in mind, and to determine whether agr-mediated changes act through rot, we performed RNA-seq experiments after 1 hr growth from overnight cultures of a Δagr Δrot double mutant that phenocopies wild-type with respect to H2O2-mediated death and with respect to its parental Δagr strain (Supplementary file 3). Fold-changes and number of genes differentially expressed were lower in the Δagr mutant relative to the wild-type culture, potentially because a significant portion of the population, even after an hour of growth (early exponential phase), still consisted of cells experiencing stationary phase at the time of sampling. Nevertheless, we did observe a shift in the expression of fermentation-associated genes (ilvA, pflAB, aldh1, ddh, lctp2) in the Δagr strain (Figure 4C and Supplementary file 3). Thus, up-regulation of metabolic genes in the Δagr mutant extends beyond post-exponential growth to the exit from stationary phase and into subsequent cell proliferation, as does the long-lived protection from H2O2-mediated killing seen with the wild-type strain.
To examine the induction of genes by lethal levels of H2O2, our gene expression analysis included a comparison between untreated and H2O2-treated cells after growth from overnight culture (Supplementary file 3). The Δagr Δrot double mutant that phenocopies wild-type had elevated expression of many genes involved in lowering oxidative stress compared to the Δagr mutant. Those genes are involved in the regulation of misfolded proteins (mcsA, mcsB, clpC, clpB), Fe-S cluster repair (iscS), DNA protection and repair (dps), and genes regulated by the protein-damage repair gene bshA (fhuB/G, queC-E) (Posada et al., 2014; Figure 7, Figure 7—figure supplement 1, and Supplementary file 3). Elevated expression of protective genes suggests that the double mutant survives damage from H2O2 better because protective genes are rendered inducible (loss of Rot-mediated repression). Overall, the data show that agr wild-type cells assume a long-lived stage after activation at high cell density in which they are primed to express genes (e.g. clpB/C, dps) that protect against high levels of exogenous oxidative stress.
Endogenous ROS is involved in agr-mediated protection from lethal, exogenous H2O2 stress
We next monitored the effect of reducing respiration and ATP levels by adding subinhibitory doses of the redox cycling agent menadione (Rowe et al., 2020) to cultures of Δagr and wild-type cells prior to lethal levels of H2O2. Addition of menadione for 30 min, which induces a burst of ROS that inactivates the TCA cycle and thereby respiration (Rowe et al., 2020), protected the Δagr mutant but had little effect on the wild-type strain (Figure 8A). Menadione’s effect on respiration and ATP can be reversed by N-acetyl cysteine (Rowe et al., 2020). Addition of N-acetyl cysteine in the presence of menadione restored H2O2 susceptibility to the agr mutant (Figure 8A). Thus, blocking endogenous ROS production/accumulation reverses the lethal effect of an agr deficiency with respect to a subsequent exogenous challenge with H2O2.
Rowe et al., 2020 showed that menadione exerts its effects on endogenous ROS by inactivating the TCA cycle in S. aureus. To determine whether this mechanism can induce protection in the Δagr mutant, we inactivated the TCA cycle gene acnA in agr wild-type and Δagr strains (Figure 8—figure supplement 2). We found that ΔacnA mutation completely protected the Δagr mutant from peroxide killing after growth to late exponential growth phase but had little effect on the wild-type agr strain. This finding supports the idea that TCA cycle activity contributes to an imbalance in endogenous ROS homeostasis in the Δagr mutant, and that this shift is a critical factor for Δagr hyperlethality. When we evaluated long-lived protection by comparing survival rates of Δagr ΔacnA mutant and Δagr cells following dilution of overnight cultures and regrowth prior to challenge with H2O2, ΔacnA remained protective, but less so (Figure 8—figure supplement 2). These partial effects of an ΔacnA deficiency suggest that Δagr stimulates long-lived lethality for peroxide through both TCA-dependent and TCA-independent pathways.
S. aureus has multiple enzymes that control the endogenous production and detoxification of ROS. SodA and SodM dismutate superoxide (O2•-) to H2O2, and catalase and AhpC then convert H2O2 to water, limiting the formation of toxic hydroxyl radical (OH•). Accordingly, we asked whether mutations in these pathways affect agr-dependent phenotypes with respect to lethal H2O2 exposure. A deficiency in the sodA superoxide dismutase (Clements et al., 1999) resulted in lower survival of the wild-type strain, similar to that observed with the Δagr mutant (Figure 8B–C). The effect was reversed by complementation with sodA on a low-copy-number plasmid (Figure 8D). The ΔsodA mutation had no effect on killing with the Δagr strain. Moreover, sodA expression (Supplementary file 1) and activity levels (Figure 8E) were similar for wild-type and the Δagr mutant. Together, these observations suggest that the contribution of sodA toward protective priming by wild-type involves dismutation of low levels of endogenous superoxide generated by respiration. In contrast, endogenous levels of ROS are saturating for sodA in Δagr cells. Inactivation of sodM, which is thought to be primarily induced by exogenous oxidative stress (Gaupp et al., 2012), had no noticeable effect on the H2O2 susceptibility of the wild-type or the Δagr mutant. We conclude that scavenging enzymes, such as SodA, are better able to control the threat posed by endogenous ROS in wild-type than in Δagr cells. They render the former better able to survive a subsequent lethal dose of H2O2, a compound that freely enters cells (Imlay, 2008) and would add to endogenous ROS levels.
Other oxidative-stress-response mutations in genes encoding catalase, thiol-dependent peroxidases, and bacillithiol showed little effect on the relative lethality of H2O2 between wild-type and Δagr mutant strains (Figure 8—figure supplement 1). Thus, protection against H2O2 lethality by these genes is not agr-specific. Paradoxically, a deficiency of ahpC (ahpC::bursa), which encodes a peroxidase (Cosgrove et al., 2007), almost completely reversed the elevated killing associated with the Δagr mutation (Figure 8F). An ahpC deficiency had no effect on the response of the otherwise wild-type strain. A deficiency in other downstream genes in the ahpC operon (ahpF, SAUSA300_0377–0378) showed no effect, indicating that the protective behavior of mutant ahpC was not caused by polar effects (Figure 8—figure supplement 3).
Results with ahpC deficiencies were initially surprising, because reduced ROS detoxification should increase rather than decrease killing. Compensatory expression of other protective genes, such as katA in the ΔahpC Δagr double mutant (Cosgrove et al., 2007), might enable cells to better survive damage from subsequent stress-stimulated ROS increases. Indeed, katA expression increased >10 fold in the ΔahpC and ΔahpC Δagr double mutants (Figure 8G). Thus, ΔkatA overcomes ΔahpC-mediated protection, consistent with the idea expressed previously that katA is more protective than ahpC against high levels of exogenous oxidative stress (Cosgrove et al., 2007; Seaver and Imlay, 2001). We conclude that the protective action of an ahpC-deficient mutant is due to a pre-induced, compensatory increase in the expression of another protective catalase.
Importance of the long-lived ‘memory’ of agr-mediated protection in a murine intraperitoneal infection model
To determine whether long-lived agr-mediated protection is important for S. aureus pathogenesis, we used the mixed infection strategy (outlined in Figure 3) in which a ΔagrBD mutant is ‘primed’ in response to AIP produced by a ΔrnaIII mutant after overnight co-culture containing an equal ratio of the two mutant strains (Figure 9A and Figure 9—figure supplement 1). Then mice were infected via intraperitoneal inoculation; 2 hr later, we lavaged the peritoneal cavity and harvested organs for determination of colony forming units (CFU). By 2 hr after bacterial administration, the number of S. aureus cells injected as inoculum had declined by 1000-fold (Figure 9B and Figure 9—figure supplement 1). Mutant proportions, identified by differential plating, demonstrated that ΔagrBD cells were enriched by ~30% in both peritoneum and organs compared to the ΔrnaIII mutant. The fraction of ΔagrBD (rnaIII+) mutants in sites of bacterial dissemination (heart, kidney, liver, lungs, and spleen) was similar to their elevated fraction in the peritoneum, thereby suggesting that agr enhances intraperitoneal infection and access to, rather than entry into extraperitoneal organs. In a control infection in which ΔagrBD was ‘unprimed’ by mixing ΔagrBD and ΔrnaIII mutants immediately before growth from stationary phase, the proportion of ΔagrBD bacterial burden was lower at all tissue sites (Figure 9A and Figure 9—figure supplement 1). This drop represented a decline in long-lived agr induction of virulence.
To study agr-ROS effects during infection, we repeated in vivo studies using Cybb−/− mice deficient in enzymes associated with host phagocyte production of ROS (the gp91 [phox] component of the phagocyte NADPH oxidase)(Pollock et al., 1995). We found that agr-mediated priming (mixing ΔagrBD and ΔrnaIII before overnight co-culture) failed to increase hematogenous dissemination to lung and spleen tissues following infection of Cybb−/− mice (Figure 9). Thus, when the host makes little ROS, long-lived agr-mediated protection has little effect in these tissues. The data also indicate that agr-mediated protection against ROS enhances fitness in lung and spleen, but it is dispensable for full virulence in other organs. Collectively, the murine experiments indicate that the long-lived ‘memory’ of agr induction enhances overall pathogenicity of S. aureus during sepsis. They also support data previously published indicating that clearance of disseminating bloodstream pathogens (Yipp et al., 2017) and protection from ROS buildup (Beavers et al., 2021) are tuned to particular sites in the host organism.
Discussion
We report that agr, a quorum-sensing regulator of virulence in S. aureus, provides surprisingly long-lived protection from the lethal action of exogenous H2O2. The protection, which is uncoupled from agr activation kinetics, arises in part from limiting the accumulation of endogenous ROS. This apparent tolerance to lethal stress derives from an RNAIII-rot regulatory connection that couples virulence-factor production to metabolism and thereby to levels of ROS. Collectively, the results suggest that agr anticipates and protects the bacterium from increases in ROS expected from the host during S. aureus infection.
Details of agr-mediated protection are sketched in Figure 10. At low levels of ROS, agr is activated by a redox sensor in AgrA, RNAIII is expressed and represses the Rot repressor, thereby rendering protective genes (e.g. clpB/C, dps) inducible via an unknown mechanism (induction, candidate protective gene(s), and their connection to endogenous ROS levels are being pursued, independent of the current report). Superoxide dismutase and scavenging catalases/peroxidases detoxify superoxide and peroxide, respectively (scavenging deficiencies reduce the protective effect of wild-type). Deletion of agr eliminates expression of RNAIII and repression of Rot, resulting in a metabolic instability associated with a 100-fold increase in H2O2-mediated death.
The agr system directly reduces H2O2-mediated killing by reducing levels of endogenous ROS, much like intrinsic tolerance to lethal antimicrobial stress (Zeng et al., 2022). However, the protective system we describe is distinct in that it primes cells for induction of genes (e.g. clpB/C, dps) that mitigate damage upon subsequent exposure to high levels of ROS. Still unidentified protective genes exist; thus, agr-mediated protection may be further shaped by both known (ahpC) and unidentified pathways and factors that modify the redox state. Another distinctive feature of agr-mediated protection is its manifestation even in early log-phase cultures, long after the maximal transcription of agr at high cell density, i.e., quorum. In a sense, S. aureus has a ‘memory’ of the agr-activated state.
Transcriptional profiling during growth from diluted, overnight cultures revealed that the Δagr mutation elevated the expression of several respiration and fermentation genes. Acceleration of cellular respiration is likely a source of ROS, as it appears to be for bactericidal antibiotics (Kohanski et al., 2007). Our work supports this idea by showing that increased respiration caused by deletion of agr is associated with increased ROS-mediated lethality. How agr deficiency is connected to the corruption of downstream processes that result in metabolic inefficiency and increased endogenous ROS levels is unknown. Given that Δagr mutants are unable to downregulate surface proteins during stationary phase (Morfeldt et al., 1995; Novick et al., 1993), it is possible that deletion of agr perturbs the cytoplasmic membrane or the machinery that sorts proteins across the cell wall. In support of this notion, jamming SecY translocation machinery of E. coli results in downstream events shared with antibiotic lethality, including accelerated respiration and accumulation of ROS (Takahashi et al., 2017). In this scenario, the formation of a futile macromolecular cycle may accelerate cellular respiration to meet the metabolic demand of unresolvable problems caused by elevated surface sorting.
As noted above, agr is inactivated by oxidation, which elevates levels of the antioxidant BsaA during exposure to H2O2 (Sun et al., 2012). That would make our finding that H2O2-mediated killing is increased in the Δagr mutant paradoxical. This apparent inconsistency can be explained by a focus of prior work on growth-related phenotypes (Sun et al., 2012) rather than on lethality (the underlying mechanisms are distinct Drlica and Zhao, 2021). Additionally, we note that bsaA was not upregulated in either our RNA-seq experiments (Supplementary file 1) or in previous transcriptional profiling data (George et al., 2019). Thus, an alternative, but not mutually exclusive, hypothesis is that the growth-related effect of bsaA on agr-mediated responses to stress is strain-dependent. Another complexity involves test conditions, as indicated by consideration of previous work in which wild-type cells exhibited greater oxidative stress than the agr-deficient mutant due to agrA-mediated production of ROS-inducing phenol-soluble modulins (George et al., 2019). The present experiments were performed in highly diluted cultures in which levels of these modulins are likely low (Queck et al., 2008; Wang et al., 2007). The complex relationship between agr, ROS-mediated lethality, and physiological state illustrates the importance of understanding agr biology before applying therapies that inactivate agr (Khan et al., 2015).
We also note that although the absence of agr increases killing by high levels of H2O2, it has the opposite effect on lethal concentrations of ciprofloxacin (Kumar et al., 2017). In the latter case, the absence of agr upregulates the expression of bsaA in the strain examined; bsaA counters endogenous ROS induced by ciprofloxacin (Kumar et al., 2017). The present work shows that excess endogenous ROS is generated during agr deficiency. Thus, protection from endogenous lethal stress via agr inactivation may not be only through the redox-dependent bsaA but also by a second pathway involving increased respiratory metabolism. The present work also supports the idea that exposing bacteria to exogenous H2O2 does not fully recapitulate the intracellular environment created by antibiotics and other stresses that act via ROS-mediated cell death (Takahashi et al., 2017), emphasizing that inactivation of agr can be either destructive or protective, depending on the type of lethal stress. Similar results have been reported with other bacteria: mazF, lepA, and cpx are destructive or protective based on the level of lethal stress (Dorsey-Oresto et al., 2013; Wu et al., 2011).
The protective activity of agr was carried over to in vivo studies using mice, as it was largely absent if the mice were deficient in host phagocyte production of ROS (Cybb−/− mice with a null allele of NADPH oxidase). The benefits of agr to S. aureus fitness seen with NADPH oxidase-proficient mice were observed largely in lungs, a key host defense niche for neutrophil-mediated clearance of disseminating pathogens (Yipp et al., 2017). The redox switch in AgrA, plus the protective properties associated with agr activation, lead to a clinical model in which agr links virulence-factor expression to an intrinsic protection against a lethal, H2O2-mediated immune response during infection (Figure 11). In this model, agr quorum-sensing renders cells better prepared to respond to lethal, exogenous oxidative stress. We note, however, that agr-mediated fitness benefits were present in certain tissues even in NADPH oxidase-deficient mice, indicating the existence of long-lived factors other than those that suppress oxidative stress. Thus, such a pre-emptive defense system may apply to many challenges experienced by S. aureus during infection, especially during bloodstream dissemination and conditions within inflamed tissues (Richardson et al., 2008; Vitko et al., 2015).
In conclusion, uncoupling of agr-mediated tolerance from bacterial population density anticipates increases in exogenous ROS expected during S. aureus-host interactions, thereby contributing to virulence. The ubiquity of quorum sensing suggests that it protects many bacterial species from oxidative damage. The next step is to find RNA, protein, and/or epigenetic markers underlying the agr-mediated ‘memory’ that improves protection against subsequent H2O2 exposure, since that will provide insights into the role of agr in cellular survival and adaptation during infection. Discovering ways to manipulate the lethal stress response, as seen with supplementation of antimicrobials with N-acetyl-cysteine during treatment of Mycobacterium tuberculosis (Vilchèze and Jacobs, 2023) and development of inhibitors of enzymes that produce protective H2S (Shatalin et al., 2021), could reveal novel approaches for enhancing antimicrobial therapy and host defense systems (Cao et al., 2017; Gusarov et al., 2009; Shatalin et al., 2011).
Materials and methods
Bacterial strains, plasmids, primers, and growth conditions
Request a detailed protocolS. aureus strains, plasmids, and primers used in the study are described in Tables 1 and 2. Bacterial cells were grown in Tryptic Soy Broth (TSB, glucose concentration at 2.5 g/L) at 37 °C with rotary shaking at 180 rpm. For suboptimal aeration, broth cultures were grown in a closed-capped 15 mL conical tube with 10 mL of TSB. Colony formation was on Tryptic Soy Agar (TSA) with or without defibrinated sheep blood, incubated at 37 °C or 30 °C. Phages 80α and Φ11-mediated transduction was used for strain construction (Novick, 1991); transductants were selected on TSA plates containing the appropriate antimicrobial.
For analysis of in vitro growth curves, overnight cultures grown in TSB were diluted 1:1000 in CDM (Hussain et al., 1991), and growth was monitored at 37 °C in 96-well plates (100 µL/well) using an Agilent LogPhase 600 Microbiology Reader (Santa Clara, CA) with 1 mm orbital shaking, measuring OD600 at 40 min intervals. The curves represent averaged values from five biological replicates. The exponential phase was used to determine growth rate (μ) from two datapoints, OD1 and OD2 flanking the linear portion of the growth curve, following the equation lnOD2-lnOD1/t2-t1, as described (Grosser et al., 2016).
Measurement of bioluminescence
Request a detailed protocolOvernight cultures were diluted to OD600 ~0.05 and grown in TSB at 37 °C with rotary shaking at 180 rpm. Aliquots (100 μL) were inoculated into flat bottom 96-well microtiter plates (Corning, Corning, NY), and bioluminescence was detected using a BioTek Synergy Neo2 plate reader (Agilent, Santa Clara, CA).
Antimicrobials and chemicals
Request a detailed protocolAntimicrobials, chemicals, and reagents were obtained from MilliporeSigma (Burlington, MA) or Thermo Fisher Scientific (Waltham, MA).
Construction of mutants
Request a detailed protocolTransposon mutants were generated by transducing Bursa aurealis insertions, obtained from the University of Nebraska transposon mutant (ΦNE) library (Fey et al., 2013), into LAC or LAC agr::tetM using phages 80α and Φ11.
Construction of the ΔagrBD mutant: S. aureus ΔagrBD mutant GAW128 was generated by a chromosomal integration strategy outlined in Chen et al., 2014 in an agr-null background strain, BS687. Plasmid pJC1111, a suicide plasmid containing a cadmium resistance (CdR) cassette and the SaPI1 attS site that enables single-copy insertion into the corresponding chromosomal SaPI1 attC site, was used as the backbone vector for the S. aureus agrBD construct. pJC1000 contains the RN6734 agr locus cloned into pUC18. Inverse PCR of pJC1000 was performed using agrBD primers GWO#27 and GWO#28, re-ligated following treatment with polynucleotide kinase, and designated pGAW98. The SphI-EcoRI fragment of pGAW98 was ligated into the SaPI1 integration vector pJC1111 and designated pGAW119. Strain RN9011 (RN4220 with pRN7023 [vector (CmR) containing SaPI1 integrase]) was electroporated with plasmid pGAW119 and plated on GL agar containing 0.1 mM CdCl2. Phage 80α lysates of CdR colonies were used to transduce BS687 (RN6734 Δagr::ermC, Erm), generating GAW128 (ΔagrBD).
To construct agr mutant BS687, agr flanking regions were amplified with primer pairs JCO#339, JCO#340, and JCO#342, JCO#343 and cloned into the HincII site of pUC18 to generate pJC1527 and pJC1528, respectively. pJC1530 was generated by four-way ligation of the KpnI-BamHI fragment of pJC1527 (agr left flank), XhoI-SphI fragment of pJC1528 (agr right flank), and BamHI-XhoI fragment from pJC1073 (Erm cassette) to KpnI-SphI digested pJC1202 (replacement vector). Plasmid pJC1530 was electroporated into strain RN4220 with selection on GL agar containing 10 µg/mL of chloramphenicol at 30 °C. Phage 80α lysates of CmR colonies were used to transduce strain JCSA18 (rpsL* mutant of RN6734 that results in streptomycin resistance) and then allelic exchange of the EmR SmS CmR colonies was performed as previously described. Phage 80 a lysates of EmR SmR CmS colonies were then used to transduce RN6734 with selection for EmR, generating BS687. sodA complementation: Plasmid PsarA-sodA-pJC1111, expressing sodA under the control of the constitutive promoter PsarA, was integrated into the S. aureus chromosome at the SaPI1 attC site of strain LAC (Geisinger et al., 2008), LAC sodA::tetM, and LAC agr::ermC. Complementation plasmid PsarA-sodRBS-sodA was generated by Gibson assembly and inserted into the SaPI1 integration vector pJC1111. Wild-type sodA and the sarA promoter were amplified from S. aureus gDNA using primers MPsodA#1–2 (sodA gene and RBS) and MPsodA#3–4 (PsarA). Primers MPsodA#5–6 were used to linearize pJC1111. Primers introduced relevant oligonucleotide overlaps that enabled Gibson assembly (Shee et al., 2022), generating PsarA-sodRBS-sodA. PsarA-sodRBS-sodA was transformed into E. coli DH5α for amplification, purification, and sequence validation via primers pJC1111 FOR and pJC1111 REV. Purified PsarA-sodRBS-sodA was electroporated into RN9011 and positive chromosomal integrants at the SaPI1 chromosomal attachment (attC) site were selected with 0.1 mM CdCl2. Phage 80 a lysates of positive integrants were used to transduce BS1422 (LAC sod::tetM) and BS1348 (LAC agr::tetM), generating BS1707 and BS1708, respectively.
Survival measurements
Request a detailed protocolTo measure lethal action, overnight cultures were diluted (OD600∼0.05) in fresh medium and grown with shaking to early exponential (OD600∼0.15) or late log (OD600∼4) phase, conditions when agr expression is largely absent (Kumar et al., 2017) or maximally activated, respectively. Early (undiluted) and late exponential phase cultures (diluted into fresh TSB medium to OD600∼0.15) were incubated with H2O2 under aerobic conditions either at a fixed concentration for one or more time points or at various concentrations for a fixed time. At the end of treatment, aliquots were removed, concentrated by centrifugation and serially diluted in phosphate-buffered saline to remove H2O2, and plated for determination of viable counts at 24 hr. The percentage of survival was calculated relative to a sample taken at the time of H2O2 addition. When menadione and N-acetylcysteine were used to inhibit or potentiate killing by H2O2, they were added prior to lethal treatments as described previously (Conlon et al., 2016). For experiments involving menadione pretreatment, cultures were grown for 3.5 hr, and menadione (40 mM solution in 96% EtOH, final concentration 80 μM) was added for the last 0.5 hr of culture, preceding the H2O2 treatment at 4 hr. N-acetylcysteine was used to counter the action of menadione; it was added simultaneously with menadione, at a final concentration of 30 mM (640 mM stock in sterile ddH2O was used). All experiments were repeated at least three times; similar results were obtained from the biological replicates.
Measurement of glucose consumption
Request a detailed protocolOvernight cultures were diluted into fresh TSB (OD6000.05) and grown for 4 hr with shaking at 180 rpm (OD6004) at 37 °C. Glucose was assayed in the supernatant fluids of bacterial cultures following centrifugation at 12,000 × g, using Centricon-10 concentrators (MilliporeSigma, Burlington, MA), and pH adjustment to 6.5–7.0 using NaOH. Cells were assayed and plated hourly for determination of viable counts as indicated in figures. Glucose content was measured from serial dilutions of supernatants using the UV method (cat. no. 10-716-251-035) following manufacturer’s instructions (R-Biopharm, Darmstadt, Germany). Glucose consumption was expressed asμg of glucose consumed over 3 hr of culture per 108 bacterial cells. There was no detectable glucose in culture supernatants at 4 hr of culture (data not shown).
Measurement of excreted metabolites
Request a detailed protocolExcreted metabolites were assayed in the supernatant fluids of bacterial cultures following centrifugation at 12,000 × g for 10 min for late exponential (4 hr, OD600~4) or multiple time points (acetate), as indicated in figures. Aliquoted supernatants were stored at −80 °C and thawed on ice prior to analysis. Cells were plated for determination of viable counts; L(+)-lactate and acetate concentrations were measured using commercially available colorimetric and fluorometric kits (cat. no. MAK065, ab204719), according to manufacturer’s recommendations (MilliporeSigma, Burlington, MA and Abcam, Cambridge, UK, respectively).
Measurement of intracellular metabolites
Request a detailed protocolOvernight cultures were diluted into fresh TSB (OD600∼0.05), grown for 4 hr at 37 °C with shaking at 180 rpm (OD600∼4), and plated for determination of viable counts at 4 hr. The remaining cells were concentrated by centrifugation at 12,000 x× g for 10 min, and resuspended in lysis buffer provided by the assay kit. Cells were lysed by repeated homogenization (two cycles of 45 s homogenization time at 6 M/s followed by a 5 min pause on ice) using Lysing Matrix B tubes in a FastPrep-24 homogenizer (MP Biomedicals, Irvine, CA). After lysis, cell debris was removed by centrifugation (12,000 × g, 10 min) and the supernatant was used for determination of pyruvate, fumarate, citrate, and acetyl-CoA levels using colorimetric (pyruvate, fumarate), or flurometric (citrate, acetyl-CoA) assays (cat. no. KA1674, ab102516, KA3791, and MAK039, respectively) and a microplate reader (BioTek Synergy Neo2, Agilent, Santa Clara, CA) according to the manufacturer’s instructions (Abnova, Taipei City, Taiwan and Abcam, Cambridge, UK and MilliporeSigma, Burlington, MA, respectively). Assayed metabolites were measured in μg and normalized to cell count.
Measurement of oxygen consumption
Request a detailed protocolOvernight cultures were diluted into fresh TSB (OD600∼0.05), grown for 5 hr at 37 °C with shaking at 180 rpm (OD600∼4), diluted (OD600∼0.025) in fresh TSB, and added to a microtiter plate (200 μL/well). Oxygen consumption rate (OCR) was measured using a Seahorse XF HS Mini Analyzer (Agilent, Santa Clara, CA) according to the manufacturer’s instructions. The Seahorse XF sensor cartridge was hydrated in a non-CO2 37 °C incubator with sterile water (overnight) and pre-warmed XF calibrant for 1 hr prior to measurement. OCR measurements were recorded in 15 measurement cycles with 3 min of measurement and 3 min of mixing per cycle. CFU were enumerated to confirm equal concentrations of agr-deficient mutant and wild-type cells.
Measurement of ATP, NAD+, and NADH
Request a detailed protocolFor ATP, overnight cultures were diluted into fresh TSB (OD600∼0.05), grown for 4 hr at 37 °C with shaking at 180 rpm (OD600∼4), diluted (OD600∼1.0) in fresh TSB, and incubated at room temperature with reagent for determination of ATP using BacTiter-Glo Microbial Cell Viability Assay (cat. no. G8232; Promega, Madison, WI), according to the manufacturer’s instructions. Luminescence was detected in a BioTek Synergy Neo2 plate reader (Agilent, Santa Clara, CA). The amount of ATP was calculated and normalized to cell count.
For NAD+ and NADH, overnight cultures were diluted into fresh TSB (OD600∼0.05), grown for 4 hr at 37 °C with shaking at 180 rpm (OD600∼4), and plated for viable counts at 4 hr or concentrated by centrifugation at 12,000 × g for 10 min and resuspended in lysis buffer provided by the assay kit. Cells were lysed by repeated homogenization (two cycles of 45 s homogenization time at 6 M/s followed by a 5 min pause on ice) using Lysing Matrix B tubes in a FastPrep-24 homogenizer (MP Biomedicals, Irvine, CA). After lysis, cell debris was removed by centrifugation (12,000 × g, 10 min), and the supernatant was used for determination of NAD+ and NADH levels using a colorimetric assay kit (cat. no. KA1657; Abnova, Taipei City, Taiwan) and a microplate reader (BioTek Synergy Neo2, Agilent, Santa Clara, CA) according to the manufacturer’s instructions.
Measurement of baseline ROS levels
Request a detailed protocolOvernight cultures were diluted into fresh TSB (OD600∼0.05), and grown with shaking to early exponential phase (OD600∼0.2). 200 µL of culture was removed and cell density was normalized before staining with carboxy-H2DCFDA fluorescent dye (final concentration 10 µM) (Invitrogen, Waltham, MA). Samples were incubated at room temperature for 5 min, then 800 µL of PBS + EDTA buffer (100 mM) was added to each sample, and ROS levels were measured by fluorescence-based flow cytometry (BD Fortessa, BD Biosciences, San Jose, CA). All tubes with cultures were wrapped with aluminum foil to avoid light. A sample containing LAC wild-type cells lacking carboxy-H2DCFDA was included as a control for auto-fluorescence. Forward and side scatter parameters were acquired with logarithmic amplification. ROS was detected using the 488 laser and a 530/30 nm bandpass filter. Data were analyzed using FlowJo software version 10.8.1 (BD Biosciences, San Jose, CA).
Measurement of superoxide dismutase (SOD) activity
Request a detailed protocolOvernight cultures were diluted (OD600∼0.05) into fresh TSB, grown to late exponential phase (OD600~4), diluted to OD600=1, centrifuged at 12,000 × g for 5 min, and the cell pellet was homogenized in 300 μL of ice-cold lysis buffer (100 mM Tris-HCl pH 7.4 + 0.5% Triton + 5 mM 2-mercaptoethanol + 0.2 mM PMSF). SOD activity was measured using a commercially available kit (cat. no CS0009-1KT), according to the manufacturers’ instructions (MilliporeSigma, Burlington, MA). The experiment was repeated three times with similar results.
RNA sequencing and data analysis
Request a detailed protocolOvernight cultures were diluted (OD600∼0.05) into fresh TSB medium and grown at 37 °C to early exponential phase (OD600~0.5) (Δagr single mutant and Δagr Δrot double mutant) or late exponential phase (OD600~4) (wild-type and Δagr strains). Samples of Δagr and Δagr Δrot were divided into two 3 mL aliquots, and the aliquots were incubated at 37 °C for another 30 min, with or without treatment with H2O2. Peroxide concentrations for Δagr and Δagr Δrot were normalized to expected killing at the time of harvest (Figure 7—figure supplement 1).
Three independent cultures for each sample were used for determination of transcriptional profiles. Briefly, cultures were concentrated by centrifugation (12,000 × g for 5 min), and resuspended cells were disrupted using Lysing Matrix B tubes in a FastPrep-24 homogenizer (MP Biomedicals, Irvine, CA) at 6 M/s, for 30 s, three times (samples were resting on ice between homogenizer runs), and RNAs were extracted from the collected bacterial cells using TRIzol reagent (Thermo Fisher Scientific, Waltham, MA). RNA was isolated using RNeasy (Qiagen, Germantown, MD) mini spin columns. Sequence libraries were generated using the TruSeq Stranded Total RNA Library Prep kit (Illumina, San Diego, CA) following the manufacturer’s recommendations. The rRNAs were removed by the Ribo-zero Kit (Illumina) to enrich mRNA, using 13 cycles of PCR amplification of the final library. Amplified libraries were purified using AMPure beads (Beckman Coulter, Brea, CA), quantified by Qubit (Thermo Fisher Scientific, Waltham, MA) and qPCR, and visualized in an Agilent Bioanalyzer (Santa Clara, CA). Pooled libraries were sequenced as paired-end 50 bp reads using an Illumina NovaSeq instrument.
Reads were initially trimmed using Trimmomatic version 0.39 (Bolger et al., 2014) to remove adaptors as well as leading or trailing bases with a quality score less than 3, filtering reads with a minimum length of 36. Reads were mapped to reference strain USA300 FPR3757 (RefSeq identifier GCF_000013465.1) using Bowtie2 version 2.2.5 (Langmead and Salzberg, 2012). Using gene annotations from the same assembly, reads mapped to each gene were counted with featureCounts version 2.0.1 (Liao et al., 2014), producing a counts matrix. Additional analysis was performed in R (R Core Team 2021) using the package DESeq2 version 1.32 (Love et al., 2014).
Normalization to account for inter-sample library size variation was performed using the built-in normalization function of DESeq2. All RNA-seq heatmaps were colored according to row (gene) z-scores of DESeq2 normalized counts. For differential expression testing, the Wald test was used with a log2 fold-change threshold of 0.5 and an FDR of 0.1. For simple pairwise comparisons (e.g. the effect of strain under control conditions), datasets were split so that analysis was performed independently for strains used in the comparison. To determine the interaction between strain and condition variables, all samples were included with the experimental design (formula ‘expression ~condition + strain + condition:strain,’ where condition:strain is the interaction between variables).
Metabolic flux prediction
Request a detailed protocolThe SPOT (Simplified Pearson cOrrelation with Transcriptomic data) computational method (Kim et al., 2016) was used to analyze the difference in intracellular metabolic fluxes between wild-type LAC and agr::tetM mutant grown in TSB to late exponential phase (OD600~4). SPOT is similar to the E-Flux2 method described previously (Balasubramanian et al., 2016), but a recent validation study (Bhadra-Lobo et al., 2020) shows that SPOT generally outperforms E-Flux2. SPOT infers metabolic flux distribution by integrating transcriptomic data in a genome-scale metabolic model of S. aureus (Becker and Palsson, 2005) that was adapted for use with strain LAC. For a list of the metabolic reactions ranked by unit of flux per 100 units of glucose uptake flux, see Supplementary file 3.
Real-time qRT-PCR assays
Request a detailed protocolBriefly, RNA was purified as described above from late exponential (OD600~4.0) cells, cDNAs were synthesized using Maxima First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA), and real-time reverse transcription quantitative PCR (qRT-PCR) was performed using QuantiNova SYBR Green PCR Kit (Qiagen, Hilden, Germany). Primers were synthesized by IDT Inc (Coralville, IA). Three independent biological samples were run in triplicate and rpoB was used to normalize gene expression. 2–∆∆Ct method was used to calculate the relative fold gene expression (Livak and Schmittgen, 2001).
Peritoneal infection of mice
Request a detailed protocolC57BL/6 mice and C57BL/6 Cybb-/- (also known as gp91phox/nox2) were purchased from the Jackson Laboratory and bred onsite to generate animals for experimentation. Age and gender-matched, 8–10 week-old mice were used. S. aureus strains harboring RNAIII or agrBD deletion in the NCTC 8325 background were grown overnight in TSB (37 °C, 180 rpm) separately or mixed at a 1:1 ratio. Overnight cultures were diluted (OD600∼0.05) into fresh TSB medium (subcultured separately for the cultures mixed overnight (‘primed’) or mixed 1:1 for RNAIII or agrBD mutant single cultures (‘unprimed’) and grown at 37 °C to early exponential phase (OD600~0.5)). Bacteria were washed one time by centrifugation with PBS and adjusted to 109 CFU/mL. Twenty C57BL/6 WT mice and 17 Cybb-/- mice were injected intraperitoneally with 100 μL of either ‘primed’ or ‘unprimed’ inoculum. After 2 hr, internal organs, peritoneal lavage, and blood were collected. The organs were homogenized in sterile PBS and serial dilutions were plated for viable counts on TS agar. Collected blood was lysed with saponin and plated for viable counts on TSA plates. Peritoneal lavage fluid was serially diluted and plated for viable counts. All animal studies were performed as per an NYU Grossman School of Medicine Institutional Animal Care and Use Committee (IACUC) approved protocol for the Shopsin Lab.
Statistical analysis
Request a detailed protocolPrism software (GraphPad, Inc) was used to perform statistical analyses.
Statistical significance was determined using the Student’s t-test, Mann–Whitney U test, one-way analysis of variance (ANOVA), or the Kruskal-Wallis test, depending on the data type. Statistical significance was considered to be represented by p values of <0.05.
Data availability
Sequencing data have been deposited in GEO under accession code GSE207045.
-
NCBI Gene Expression OmnibusID GSE207045. The quorum-sensing agr system protects Staphylococcus aureus from oxidative stress.
References
-
Microbial stress priming: a meta-analysisEnvironmental Microbiology 18:1277–1288.https://doi.org/10.1111/1462-2920.13223
-
Staphylococcus aureus peptide methionine sulfoxide reductases protect from human whole-blood killingInfection and Immunity 89:e0014621.https://doi.org/10.1128/IAI.00146-21
-
Rot and SaeRS cooperate to activate expression of the staphylococcal superantigen-like exoproteinsJournal of Bacteriology 194:4355–4365.https://doi.org/10.1128/JB.00706-12
-
Trimmomatic: a flexible trimmer for Illumina sequence dataBioinformatics 30:2114–2120.https://doi.org/10.1093/bioinformatics/btu170
-
Bacterial death from treatment with fluoroquinolones and other lethal stressorsExpert Review of Anti-Infective Therapy 19:601–618.https://doi.org/10.1080/14787210.2021.1840353
-
Staphylococcal coagulase; mode of action and antigenicityJournal of General Microbiology 6:95–107.https://doi.org/10.1099/00221287-6-1-2-95
-
Anaerobic gene expression in Staphylococcus aureusJournal of Bacteriology 189:4275–4289.https://doi.org/10.1128/JB.00081-07
-
Staphylococcal response to oxidative stressFrontiers in Cellular and Infection Microbiology 2:33.https://doi.org/10.3389/fcimb.2012.00033
-
Inhibition of rot translation by RNAIII, a key feature of agr functionMolecular Microbiology 61:1038–1048.https://doi.org/10.1111/j.1365-2958.2006.05292.x
-
Identification of ligand specificity determinants in AgrC, the Staphylococcus aureus quorum-sensing receptorThe Journal of Biological Chemistry 283:8930–8938.https://doi.org/10.1074/jbc.M710227200
-
Regulatory requirements for Staphylococcus aureus nitric oxide resistanceJournal of Bacteriology 198:2043–2055.https://doi.org/10.1128/JB.00229-16
-
A chemically defined medium for slime production by coagulase-negative staphylococciJournal of Medical Microbiology 34:143–147.https://doi.org/10.1099/00222615-34-3-143
-
Cellular defenses against superoxide and hydrogen peroxideAnnual Review of Biochemistry 77:755–776.https://doi.org/10.1146/annurev.biochem.77.061606.161055
-
Investigational therapies targeting quorum-sensing for the treatment of Staphylococcus aureus infectionsExpert Opinion on Investigational Drugs 24:689–704.https://doi.org/10.1517/13543784.2015.1019062
-
Fast gapped-read alignment with Bowtie 2Nature Methods 9:357–359.https://doi.org/10.1038/nmeth.1923
-
Transcription profiling of the mgrA Regulon in Staphylococcus aureusJournal of Bacteriology 188:1899–1910.https://doi.org/10.1128/JB.188.5.1899-1910.2006
-
Genetic systems in staphylococciMethods in Enzymology 204:587–636.https://doi.org/10.1016/0076-6879(91)04029-n
-
The agr P2 operon: an autocatalytic sensory transduction system in Staphylococcus aureusMolecular & General Genetics 248:446–458.https://doi.org/10.1007/BF02191645
-
Autoinduction and signal transduction in the regulation of staphylococcal virulenceMolecular Microbiology 48:1429–1449.https://doi.org/10.1046/j.1365-2958.2003.03526.x
-
Quorum sensing in staphylococciAnnual Review of Genetics 42:541–564.https://doi.org/10.1146/annurev.genet.42.110807.091640
-
Importance of bacillithiol in the oxidative stress response of Staphylococcus aureusInfection and Immunity 82:316–332.https://doi.org/10.1128/IAI.01074-13
-
Inactivation of the Pta-AckA pathway causes cell death in Staphylococcus aureusJournal of Bacteriology 195:3035–3044.https://doi.org/10.1128/JB.00042-13
-
Moxifloxacin-mediated killing of mycobacterium tuberculosis involves respiratory downshift, reductive stress, and accumulation of reactive oxygen speciesAntimicrobial Agents and Chemotherapy 66:e0059222.https://doi.org/10.1128/aac.00592-22
-
Neutrophils versus Staphylococcus aureus: a biological tug of warAnnual Review of Microbiology 67:629–650.https://doi.org/10.1146/annurev-micro-092412-155746
-
The promises and limitations of N-acetylcysteine as a potentiator of first-line and second-line tuberculosis drugsAntimicrobial Agents and Chemotherapy 65:e01703-20.https://doi.org/10.1128/AAC.01703-20
Article and author information
Author details
Funding
National Institute of Allergy and Infectious Diseases (R01AI137336)
- Itai Yanai
- Victor J Torres
- Bo Shopsin
National Institute of Allergy and Infectious Diseases (R01AI140754)
- Victor J Torres
- Bo Shopsin
National Institute of Allergy and Infectious Diseases (R01AI150893)
- Jeffrey N Weiser
National Institute of Allergy and Infectious Diseases (R01AI038446)
- Jeffrey N Weiser
National Institute of Allergy and Infectious Diseases (R01AI149350)
- Victor J Torres
National Institute of Allergy and Infectious Diseases (K08AI163457)
- Robert J Ulrich
National Institute of Allergy and Infectious Diseases (T32AR064184)
- Theodora K Karagounis
National Institute of Allergy and Infectious Diseases (R21AI153646)
- Dane Parker
New Jersey Health Foundation (PC 142-22)
- Dane Parker
New Jersey Commission on Cancer Research (COCR22RBG005)
- Dane Parker
NYU Langone Health Antimicrobial-Resistant Pathogens Program
- Alejandro Pironti
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Andrew Darwin for his critical comments on the manuscript. This work was supported by NIH National Institute of Allergy and Infectious Diseases grants R01AI137336 (BS, IY, and VJT); R01AI140754 (BS and VJT); R01AI150893 and R01AI038446 (JNW); R01AI149350 (VJT); K08AI163457 (RJU), T32AR064184 (TKK), and R21AI153646 (DP); New Jersey Health Foundation PC 142–22 and New Jersey Commission on Cancer Research COCR22RBG005 grants (DP); and funds from the NYU Langone Health Antimicrobial-Resistant Pathogens Program (BS, AP, and VJT). The NYU Langone Health Genome Technology Center, and the Cytometry and Cell Sorting Laboratory are shared resources that are partially supported by the Cancer Center Support Grant P30CA016087 at the Laura and Isaac Perlmutter Cancer Center.
Ethics
This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled according to approved institutional animal care and use committee (IACUC) protocol (IA16-01941) of NYU Langone Health. The protocol was approved by the The Animal Care and Use Program at the NYU Grossman School of Medicine (Assurance number: D16-00274). Every effort was made to minimize suffering.
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Reviewed Preprint version 3:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.89098. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2023, Podkowik et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,755
- views
-
- 174
- downloads
-
- 3
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Microbiology and Infectious Disease
Most malaria rapid diagnostic tests (RDTs) detect Plasmodium falciparum histidine-rich protein 2 (PfHRP2) and PfHRP3, but deletions of pfhrp2 and phfrp3 genes make parasites undetectable by RDTs. We analyzed 19,313 public whole-genome-sequenced P. falciparum field samples to understand these deletions better. Pfhrp2 deletion only occurred by chromosomal breakage with subsequent telomere healing. Pfhrp3 deletions involved loss from pfhrp3 to the telomere and showed three patterns: no other associated rearrangement with evidence of telomere healing at breakpoint (Asia; Pattern 13-TARE1); associated with duplication of a chromosome 5 segment containing multidrug-resistant-1 gene (Asia; Pattern 13-5++); and most commonly, associated with duplication of a chromosome 11 segment (Americas/Africa; Pattern 13-11++). We confirmed a 13–11 hybrid chromosome with long-read sequencing, consistent with a translocation product arising from recombination between large interchromosomal ribosome-containing segmental duplications. Within most 13-11++ parasites, the duplicated chromosome 11 segments were identical. Across parasites, multiple distinct haplotype groupings were consistent with emergence due to clonal expansion of progeny from intrastrain meiotic recombination. Together, these observations suggest negative selection normally removes 13-11++pfhrp3 deletions, and specific conditions are needed for their emergence and spread including low transmission, findings that can help refine surveillance strategies.
-
- Immunology and Inflammation
- Microbiology and Infectious Disease
The balanced gut microbiota in intestinal mucus layer plays an instrumental role in the health of the host. However, the mechanisms by which the host regulates microbial communities in the mucus layer remain largely unknown. Here, we discovered that the host regulates bacterial colonization in the gut mucus layer by producing a protein called Chitinase 3-like protein 1 (Chi3l1). Intestinal epithelial cells are stimulated by the gut microbiota to express Chi3l1. Once expressed, Chi3l1 is secreted into the mucus layer where it interacts with the gut microbiota, specifically through a component of bacterial cell walls called peptidoglycan. This interaction between Chi3l1 and bacteria is beneficial for the colonization of bacteria in the mucus, particularly for Gram-positive bacteria like Lactobacillus. Moreover, a deficiency of Chi3l1 leads to an imbalance in the gut microbiota, which exacerbates colitis induced by dextran sodium sulfate. By performing fecal microbiota transplantation from Villin-cre mice or replenishing Lactobacillus in IEC∆Chil1 mice, we were able to restore their colitis to the same level as that of Villin-cre mice. In summary, this study shows a ‘scaffold model’ for microbiota homeostasis by interaction between intestinal Chi3l1 and bacteria cell wall interaction, and it also highlights that an unbalanced gut microbiota in the intestinal mucus contributes to the development of colitis.