Impact of Bicarbonate-β-Lactam Exposures on Methicillin-Resistant Staphylococcus aureus (MRSA) Gene Expression in Bicarbonate-β-Lactam-Responsive vs. Non-Responsive Strains

Abstract

Methicillin-resistant Staphylococcus aureus (MRSA) infections represent a difficult clinical treatment issue. Recently, a novel phenotype was discovered amongst selected MRSA which exhibited enhanced β-lactam susceptibility in vitro in the presence of NaHCO₃ (termed ‘NaHCO₃-responsiveness’). This increased β-lactam susceptibility phenotype has been verified in both ex vivo and in vivo models. Mechanistic studies to-date have implicated NaHCO₃-mediated repression of genes involved in the production, as well as maturation, of the alternative penicillin-binding protein (PBP) 2a, a necessary component of MRSA β-lactam resistance. Herein, we utilized RNA-sequencing (RNA-seq) to identify genes that were differentially expressed in NaHCO₃-responsive (MRSA 11/11) vs. non-responsive (COL) strains, in the presence vs. absence of NaHCO₃-β-lactam co-exposures. These investigations revealed that NaHCO₃ selectively repressed the expression of a cadre of genes in strain 11/11 known to be a part of the sigB-sarA-agr regulon, as well as a number of genes involved in the anchoring of cell wall proteins in MRSA. Moreover, several genes related to autolysis, cell division, and cell wall biosynthesis/remodeling, were also selectively impacted by NaHCO₃-OXA exposure in the NaHCO₃-responsive strain MRSA 11/11. These outcomes provide an important framework for further studies to mechanistically verify the functional relevance of these genetic perturbations to the NaHCO₃-responsiveness phenotype in MRSA.

1. Introduction

The rise of antibiotic-resistant bacteria has increasingly threatened healthcare initiatives worldwide [1,2,3,4]. In particular, methicillin-resistant Staphylococcus aureus (MRSA) is a major clinical pathogen, causing a wide variety of severe infections, with relatively limited treatment options, poor patient outcomes, and increased healthcare-associated costs [5,6,7,8]. MRSA strains are traditionally considered as “resistant” to nearly all β-lactams (the treatment of choice for methicillin-susceptible S. aureus [MSSA]) [9,10]. Standard treatment options for MRSA, such as vancomycin or daptomycin, tend to be less effective, less well tolerated (in terms of side effects and toxicities) and more expensive than standard-of-care β-lactams as used for MSSA infections [11,12,13,14].

Recently, a novel β-lactam susceptibility phenotype was discovered among MRSA strains, termed ‘bicarbonate [NaHCO₃]-responsiveness’, wherein certain MRSA strains become sensitized to β-lactam antibiotics in vitro in the presence of NaHCO₃ [15,16]. This phenotype was relatively widespread among a small collection (n = 58) of US bloodstream isolates (ranging from 33–75% for oxacillin and cefazolin, respectively [16]). Importantly, the translatability of this phenotype has been verified both ex vivo (in simulated endocarditis vegetations [SEV] models) and in vivo (in experimental infective endocarditis models) [15,17]. To-date, the potential mechanisms of this phenotype have focused on the impact of NaHCO₃ on mecA expression/PBP2a protein production, as well as upon genes associated with PBP2a maturation and functionality [15,18]. However, these mechanistic studies have predominantly centered on targeted investigations of well-characterized genes already known for their involvement in β-lactam resistance in MRSA, rather than a more unbiased and global consideration of genes that are selectively impacted among NaHCO₃-responsive vs. non-responsive MRSA strains.

In the present investigation, we have utilized RNA-sequencing (RNA-seq) methods to identify a broader repertoire of genes likely involved in the NaHCO₃-responsive vs. non-responsive phenotypes. We concentrated on identifying genes which were selectively expressed during NaHCO₃-β-lactam co-exposures, comparing a prototype NaHCO₃-responsive strain (MRSA 11/11) to a prototype non-responsive strain (COL) [15,17,18]. These RNA-seq analyses revealed a number of differentially and selectively expressed genes in the NaHCO₃-responsive strain, with several noteworthy gene cadres, including: (i) genes within the sigB-sarA-agr regulatory axis; (ii) cell wall-anchored/peptidoglycan-associated genes; and (iii) genes involved in either global or divisome–specific autolysis.

2. Materials and Methods

2.1. Bacterial Strains and Growth Conditions

MRSA strains MRSA 11/11 (USA300) and COL (USA100) were utilized as representative NaHCO₃-responsive and -non-responsive strains, respectively [15]. These two prototype strains have been used in a number of our prior investigations related to NaHCO₃-responsiveness in vitro, ex vivo and in vivo [15,17,18]. These strains have the following mecA, agr, spa and CC types: SCCmec IV, agr I, t008, CC8 (MRSA 11/11); and SCCmec I, agr I, t008, CC8 (COL).

Strains were stored at −80 °C until thawed for use, and isolated on tryptic soy agar (TSA). Strains were grown in cation-adjusted Mueller Hinton Broth (CA-MHB; Difco) + 100 mM Tris (pH 7.3 ± 0.1) with or without 44 mM NaHCO₃ or CA-MHB + 100 mM Tris (pH 7.3 ± 0.1) + 2% NaCl + ½ x the minimum inhibitory concentration (MIC) of oxacillin (OXA), with or without 44 mM NaHCO₃ to generate RNA to compare the impact of NaHCO₃ alone or the impact of NaHCO₃ + OXA, respectively, by RNA-seq. OXA concentrations for MRSA 11/11 were 16 µg/mL (CA-MHB Tris) and 0.25 µg/mL (CA-MHB Tris + 44 mM NaHCO₃) (representing ½ the MIC as determined for this strain by broth microdilution in media without and with NaHCO₃, respectively [15]); OXA concentrations for COL were 256 µg/mL for both distinct growth conditions. We have used the above NaHCO₃ and OXA concentrations in our prior studies [15,16,17,18,19]. Although this NaHCO₃ concentration (44 mM) is above human blood levels, it replicates tissue levels [20], and maximizes the NaHCO₃-responsiveness phenotype [15,19].

2.2. RNA Isolation, Library Construction, Sequencing and Analysis

To obtain RNA, strains were grown in CA-MHB 100 mM Tris ± 44 mM NaHCO₃ overnight (O/N), then diluted 1:50 into 25 mL of either the same O/N growth medium or O/N growth medium supplemented with 2% NaCl and ½ x OXA MIC concentrations (as indicated above). Cultures were then grown for 4 h at 37 °C with aeration (to reach an OD₆₀₀ ~0.5), and total RNA was harvested by mechanical disruption by FastPrep disruption (MP Biomedicals, Irvine, CA, USA), followed by column isolation (Qiagen, Germantown, MD, USA), and treatment with Turbo DNAse (Invitrogen, Waltham, MA, USA) to remove DNA from sample [15]. Ribosomal RNA was depleted using the Ribo-Zero Magnetic Bead rRNA Removal Kit (Illumina, San Diego, CA, USA) and prepared for sequencing using the Illumina Stranded Total RNA Prep Kit (Illumina, San Diego, CA, USA). Extracted RNA quality was assessed with an Agilent TapeStation 4200 (Agilent, Santa Clara, CA, USA) using an RNA ScreenTape, assuring an RNA Integrity Number (RIN score) above 6 before moving into library preparation. Following cDNA conversion, fragment sizes were assessed using an Agilent TapeStation 4200 and a High Sensitivity DNA ScreenTape, and DNA concentration was assessed using a Qubit (Thermo Fisher Scientific, Waltham, MA, USA) and dsDNA HS Assay kit. Sequencing libraries were pooled at equimolar concentrations and sequenced on an Illumina HiSeq 4000 (Illumina, San Diego, CA, USA) using 2 × 150 paired end reads. Sequencing adapters and low-complexity reads were removed using Trimmomatic v0.38 [21], and aligned to the ASM1346v1 and ASM1204v1 RefSeq genomes with Bowtie2 v2.3.4.2 for MRSA 11/11 and COL, respectively [22]. Transcripts were assembled using Cufflinks v2.2.1, and CummeRbund v2.32.0 [23]. Two biological replicates were assessed in technical duplicate for each strain and indicated growth condition. Expression data for these samples containing NaHCO₃ are currently available in the NCBI Sequence Read Archive (SRA) database under the accession numbers SAMN21542115 (MRSA 11/11 CA-MHB Tris + NaHCO₃), SAMN21542148 (MRSA 11/11 CA-MHB Tris + OXA + NaHCO₃), SAMN21542116 (COL CA-MHB Tris + NaHCO₃), and SAMN21542147 (COL CA-MHB Tris + OXA + NaHCO₃). Codification of the accession numbers for sequence read data-basing for samples that do not contain NaHCO₃ (MRSA 11/11 CA-MHB Tris, MRSA 11/11 CA-MHB Tris + OXA, COL CA-MHB Tris, COL CA-MHB Tris + OXA), are currently in-progress and will be made available by contacting the lead Author (S.C.E.), at the following email address ([email protected]). “Differentially and selectively expressed genes” in our study represented those genes which exhibited a significant level of differential expression in the presence of NaHCO₃±OXA exposure compared to media without NaHCO₃ for the responsive strain, MRSA 11/11 vs. the non-responsive COL strain (≥ 2-fold). Thus, this definition explicitly eliminated genes which were differentially and concordantly expressed under these conditions in MRSA 11/11 as well as in COL. For selected analyses, we catalogued genes which were differentially expressed by ≥ 5-fold.

Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways and annotations for differentially abundant genes were mapped using KOBAS v3.0 [24].

2.3. qRT-PCR Validation of Select Genes

Total RNA was isolated from MRSA 11/11 and COL log phase cells (OD₆₀₀ = 0.5) grown in CA-MHB 100 mM Tris + 2% NaCl + ½ × MIC OXA ± 44 mM NaHCO₃ as previously described [15]. Briefly, RNA was extracted by column isolation (Qiagen, Germantown, MD, USA) following mechanical disruption by FastPrep (MP Biomedicals, Irvine, CA, USA). Total RNA was treated with Turbo DNAse (Invitrogen, Waltham, MA, USA) and reverse transcribed by random hexamers to generate a cDNA library with SuperScript IV (Invitrogen, Waltham, MA, USA).

Following RNA-seq analyses, we selected five genes from the overall repertoire of differentially expressed genes for qRT-PCR validations, representing broad categories of bacterial factors, including: i) peptidoglycan synthesis (pbp2 and ddh); ii) autolysis (atl and sceD); and iii) virulence (fnbA). The primers listed in Table S1 were used to amplify each gene-of-interest to determine their relative expressions. The gene gyrB was used as a housekeeping gene to normalize transcript abundance. The qRT-PCR was carried out on a StepOne thermocylcer (ThermoFisher, Waltham, MA, USA) and analyzed with StepOne Software. Relative gene expression was calculated using the 2^−ΔΔC_T method from two independent biological replicates performed in triplicate on at least two separate runs for each strain/condition. Relative gene expression in CA-MHB 100 mM Tris+2% NaCl+½ × MIC OXA+44 mM NaHCO₃ was normalized to CA-MHB 100 mM Tris+2% NaCl+½ × MIC OXA for each gene, with expression in the latter condition set equal to 1.0. statistical significances related to relative expression profiles were determined by a Student’s t-test, with a p value of ≤0.05 considered ‘significant’.

3. Results and Discussion

3.1. Differentially Expressed Genes by RNA-seq in NaHCO₃-responsive vs. Non-responsive Strains

RNA-seq analysis was carried out on two prototypical MRSA strains, MRSA 11/11 (NaHCO₃-reponsive) and COL (NaHCO₃-non-responsive), for two distinct comparison groups: Group 1—CA-MHB Tris vs. CA-MHB Tris + NaHCO₃; and Group 2—CA-MHB Tris + 2% NaCl + OXA vs. CA-MHB Tris + NaHCO₃ + 2% NaCl + OXA. RNA-seq analysis of the prototypical NaHCO₃-responsive strain MRSA 11/11 revealed that exposure to NaHCO₃ alone (Group 1) led to the differential change in abundance (≥ 2-fold change) of 21 transcripts (Figure 1A); in contrast, NaHCO₃-OXA combination (Group 2) yielded a differential change in abundance of 84 transcripts (Figure 1B). Following initial analysis of differentially abundant transcripts for Groups 1 and 2 in MRSA 11/11 and COL, those genes whose expression was altered in the same direction for both MRSA 11/11 and COL in each group were subtracted from subsequent analyses, as presented in Table 1.

NaHCO₃ exposure, in the presence and absence of β-lactam co-exposure, impacted the expression of genes across a range of functional groups, with a large proportion relating to either cell wall synthesis, cell membrane transporters, or cellular metabolism (Table 1). Intriguingly, across most functional categories, NaHCO₃ exposure resulted in significant repression of gene expression, particularly for those genes relating to virulence, cell wall biosynthesis, and membrane transporters (Table 1, Table S2A–C). Of the genes that were upregulated by NaHCO₃ exposure, these tended to be related to cellular metabolism, protein translation, and/or stress response, indicating that NaHCO₃ may mirror the exposure of the responsive MRSA strain to an ‘inhospitable’ host environment.

Oxidative and membrane stress have been implicated in the formation of persister cells [25], which make up approximately 1% of the total cell population; such cells have increased antibiotic tolerance and reduced metabolism [26,27,28]. The RNA-seq data above indicate that, while NaHCO₃ may be stimulating a stress response in NaHCO₃-responsive MRSA, NaHCO₃ also evoked increased expression of many metabolic genes as compared to cells grown in the absence of NaHCO₃. This latter finding (i.e., enhanced metabolic activity) would not support the notion of NaHCO₃-mediated entry of responsive MRSA cells into a persister state (Table 1, Table S2A–C).

We then interrogated only those differentially abundant transcriptional units that were specifically altered in MRSA 11/11, and for which a specific function is known and expression was altered by ≥ 5-fold. NaHCO₃ alone significantly impacted the expression at this magnitude of only a relatively small number of functionally-identifiable genes/transcriptional units (n = 7) as compared to a larger number of genes/transcriptional units differentially impacted by NaHCO₃+OXA (n = 14) in the NaHCO₃-responsive strain MRSA 11/11 (Figure 2). Additionally, eight transcriptional units were differentially expressed in both growth conditions (NaHCO₃ alone and NaHCO₃+OXA) (Figure 2).

Of particular interest, the autolysin sceD and the global virulence regulator agr were both substantially repressed (≥ 5-fold change in expression) by NaHCO₃, with or without β-lactam co-exposure (Figure 2, Table 2). The agr locus is a major virulence regulator in S. aureus [29,30,31], although its role in antibiotic resistance in vivo remains controversial [30,32]. Previously, we observed by qRT-PCR that NaHCO₃ is capable of repressing the global virulence regulons sarA and sigB (both agr-regulated) [15,18]. These RNA-seq data corroborate the downstream effects of sarA and sigB repression, as observed by the repression of agr and several other genes known to be part of the sigB-sarA-agr regulon, including fnbA, fnbB, cap8, clfA, and clpL (Table 2, Table 3) [33,34,35,36]. Of note, one of the few genes upregulated by NaHCO₃ exposure in MRSA 11/11 was clfB, a cell wall-anchored (CWA) protein whose function is linked to nasal colonization and virulence in skin and soft tissue infections [37]. Additionally, NaHCO₃ repressed expression of sdrH, another CWA gene whose repression is associated with decreased biofilm formation [38,39], an impact that corresponds to the observed decrease in biofilm formation during NaHCO₃ exposure in MRSA 11/11 [40]. Taken together, these data indicate a broad impact of NaHCO₃ (especially in combination with OXA) on genes associated with cell surface modifications, as well as cell wall biosynthesis or remodeling.

It should be underscored that NaHCO₃ did increase the expression of several genes in the non-responsive strain, COL, whose expressions were down-regulated in MRSA 11/11, including the cap8 operon, the virulence-associated gene vraX, and the potassium transport operon kdpABCDF (Table 3). Despite these observed differences in such putative virulence gene expressions between COL and MRSA 11/11, we have not observed any differences in the intrinsic virulence of these two strains in vivo [15]. Further, a large cadre of genes was specifically impacted by NaHCO₃, with or without OXA exposure in strain COL, but not in MRSA 11/11. However, the analytics required to determine the role of these specifically and differentially impacted genes in generating the non-responsive MRSA phenotype is beyond the scope of this paper.

Of primary interest in regard to NaHCO₃-mediated susceptibility to β-lactams was its impact on the expression of genes related to cell wall synthesis. NaHCO₃ (with or without OXA co-exposure) repressed the expression of two key autolysins, atl and sceD; in addition, the expressions of several other key genes involved in cell wall synthesis were repressed by NaHCO₃-OXA co-exposures (isaA, fmtA, ddh, pbp2) (Table 2). These latter six genes encompass those involved in the formation of nascent (new) peptidoglycan and cell wall restructuring (pbp2, ddh, fmtA) [41,42,43,44], as well as peptidoglycan hydrolases, involved in cell wall turnover and division (atl, sceD, isaA) [45,46,47]. Previous studies have demonstrated that the deletion or inactivation of isaA and fmtA results in increased susceptibility to β-lactams, an event independent of impacts on the expression of PBP2a [45,48]. Furthermore, the regulation of sceD and isaA transcription is linked, since deletion in isaA results in enhanced expression of sceD [45,46]. We observed a modest decrease in isaA transcription in the presence of NaHCO₃-OXA, but a substantial decrease in sceD transcription in the presence of NaHCO₃ (with or without OXA). Although previous reports did not find that the deletion of sceD alone resulted in increased β-lactam susceptibility [45], it is possible that within the context of the gene expression profile stimulated by NaHCO₃, reduced sceD expression may have a larger impact on β-lactam susceptibility than previously understood.

3.2. qRT-PCR Validation

The qRT-PCR validation of specific genes identified by RNA-seq revealed that sceD expression was strongly repressed by NaHCO₃ in MRSA 11/11, while being substantially upregulated in the non-responsive strain COL (Figure 3); this finding further supports the notion that altered sceD expression may be of key importance to the NaHCO₃-responsive phenotype. The qRT-PCR also validated the NaHCO₃-mediated repression of pbp2, atl, and fnbA in responsive strain MRSA 11/11 as seen by RNA-seq, while NaHCO₃ stimulated increased expression of these genes in COL (Figure 3).

3.3. Limitations of Study

This study was limited by only using a single NaHCO₃-responsive/non-responsive strain pair for RNA-seq analyses. A larger number of strains exhibiting each phenotype would increase the likelihood of identifying a more consensus cadre of differentially expressed genes involved in the NaHCO₃-responsiveness phenotype. Further, only one β-lactam (OXA) was tested for its impact on gene expression in the presence of NaHCO₃. Our previous work has demonstrated that NaHCO₃-induced sensitization to cefazolin (CFZ) is much more frequent among MRSA strains than sensitization to OXA [16]. This implies that CFZ and OXA may be impacting somewhat different genetic targets in NaHCO₃-responsive strains, which could be causing the increased frequency of sensitization to CFZ vs. OXA. Additionally, only one concentration of NaHCO₃ (44 mM) was examined for its impact on gene expression. Investigating the impact of NaHCO₃ at blood concentrations (~25 mM NaHCO₃) may reveal genes that are particularly important to NaHCO₃-β-lactam-induced changes in expression within that specific host microenvironment.

Although KEGG analysis was carried out, the relatively small number of differentially expressed genes in MRSA 11/11 limited the number that could be functionally categorized to a selected, definable metabolic or biosynthetic pathway. The relative lack of descriptive KEGG pathways to ascribe to the differentially transcribed genes in these analyses underscores the need for separate functional characterizations of those genes we have identified by RNA-seq. Furthermore, additional investigations with a larger collection of NaHCO₃-responsive strains may help in determining specific, consensus KEGG pathways associated with the responsive phenotype.

4. Conclusions

NaHCO₃ (especially in combination with OXA) impacts a variety of genes specifically and differentially expressed in the NaHCO₃-responsive strain MRSA 11/11 relating to: cell wall synthesis, turn-over and division; and cell surface-associated genes regulated by the sigB-sarA-agr regulatory axis. These analyses point to a mechanism of β-lactam sensitization by NaHCO₃ that involves global impacts on the cell surface (where β-lactams initially bind), as well as altered expression of genes required to synthesize the cell wall and maintain its homeostasis. Moreover, further work is needed to validate and quantify the specific impact of those genes identified herein on the NaHCO₃-responsive phenotype, using strategic gene knockouts. Finally, future studies will be designed to look at the differential impacts of NaHCO₃ upon protein-level expression profiles in NaHCO₃-responsive vs. non-responsive MRSA.

Author Contributions

Conceptualization, S.C.E., B.M.H., C.A.A., R.A.P., and A.S.B.; methodology, S.C.E., B.M.H., and T.T.T.; software, B.M.H.; formal analysis, S.C.E., B.M.H., and R.A.P.; writing—original draft preparation, S.C.E., B.M.H., and A.S.B.; writing—review and editing, S.C.E., B.M.H., R.A.P., H.F.C., and A.S.B.; supervision, A.S.B.; project administration, C.A.A., A.S.B.; funding acquisition, B.M.H. and A.S.B. All authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by the following grants from the National Institutes of Health: 1R01-AI146078 (to ASB), K01-AI148593 (to BMH), 1R01-AI134637 (to CAA), 1R01-AI148342 (to CAA), K24-AI121296 (to CAA), and P01-AI152999 (to CAA).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in the NCBI Sequence Read Archive (SRA) as indicated in the manuscript materials and methods.

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Tables

Figure 1. Volcano plots of the differential expression of genes in the presence of: (A) NaHCO3 exposure alone; or (B) NaHCO3 +½xMIC OXA co-exposures in NaHCO3-responsive strain MRSA 11/11. NaHCO3 alone led to an observed 21 differentially expressed transcripts. NaHCO3 and OXA co-exposures together yielded 84 differentially expressed transcripts. These numbers represent the total number of transcripts identified, featuring genes of both known and unknown functions.

Figure 2. Venn diagram of genes altered ≥ 5-fold by NaHCO3 exposure without (left) or with (right) OXA exposure specifically in NaHCO3-responsive strain MRSA 11/11. Arrows indicate whether expression was decreased (downward arrow) or increased (upward arrow) by NaHCO3 exposure. OXA exposure is equivalent to ½ × MIC of OXA in CA-MHB Tris + 2% NaCl ± NaHCO3.

Figure 3. qRT-PCR confirmation of the expression of key genes-of-interest from RNA-seq analysis in: (A) MRSA 11/11, NaHCO3-responsive strain; and (B) COL, non-responsive strain. Expression was determined in log phase MRSA 11/11 and COL strains, each grown in the presence of ½ x MIC OXA. For each gene, expression is normalized to that observed for the indicated strain grown in CA-MHB Tris+2% NaCl + ½ × MIC OXA (CA-MHB Tris + 2% NaCl + ½ × MIC OXA expression set at 1.0). Statistics were calculated by Student’s t-test (*** P < 0.001, **** P < 0.0001).

Table 1

Overview of gene functional classes altered ≥ 2-fold in two comparison groups, (Group 1) CA-MHB Tris vs. CA-MHB Tris + NaHCO₃, and (Group 2) CA-MHB Tris + 2% NaCl + OXA vs. CA-MHB Tris + NaHCO₃ + 2% NaCl + OXA. Genes/transcriptional units indicated in this tables are those altered by NaHCO₃ exposure alone (Group 1, n = 11), NaHCO₃ and OXA co-exposure (Group 2, n = 57), or common to both groups (Group 1 and 2, n = 10) for a total number of functionally identifiable genes/transcriptional units of 78. "Increase" and "decrease" refer to expression being increased or decreased by exposure to NaHCO₃ compared to media without NaHCO₃ in each comparison group.

Classifications of Gene Functions	Number of Genes/Transcriptional Units Per Category	Number Increased	Number Decreased
Virulence	12	2	10
Cell Wall Synthesis	10	1	9
Membrane	19	4	15
Metabolism	22	12	10
Transcription	5	1	4
Translation	6	5	1
Transposition	2	2	0
Stress Response	2	1	1
Totals:	78	28	50

Table 2

Genes of interest selectively altered in NaHCO₃-responsive strain MRSA 11/11 (log₂ fold change ^A).

Gene Name	Locus Tag	+NaHCO3 (no OXA)	+NaHCO3 with OXA B	Function	−2	0	2
clfA	SAUSA300_0772	−2.18	−2.97	virulence
clfB	SAUSA300_2565	1.29	1.11
agr	SAUSA300_1989	−3.21	−5.37
sdrH	SAUSA300_1985	−2.44	0.00
fnbA	SAUSA300_2441	0.00	−1.84
fnbB	SAUSA300_2440	0.00	−1.47
atl	SAUSA300_0955	−1.65	−2.93	autolysins
sceD	SAUSA300_2051	−6.53	−4.41
isaA	SAUSA300_2506	0.00	−2.12
fmtA	SAUSA300_0959	0.00	−1.32	peptidoglycan/cell wall synthesis
ddh	SAUSA300_2463	0.00	−4.00
pbp2	SAUSA300_1341	0.00	−1.81
bccT	SAUSA300_2549	−2.45	−2.42	osmotic stress
usp	SAUSA300_0067	−1.30	−2.28	stress response

^A Fold change comparisons are made by comparing the gene expression for the indicated strain and condition (± OXA) exposed to NaHCO₃ to the same strain/condition without NaHCO₃ exposure. Negative values indicate that exposure to NaHCO₃ for the given strain/condition decreases expression of the indicated gene, whereas positive values indicate that NaHCO₃ increases expression of the indicated gene. ^B OXA concentration equal to 1/2 × MIC in CA-MHB Tris ± NaHCO3 for the indicated strain.

Table 3

Genes of interest altered in both NaHCO₃-responsive and non-responsive strains (log₂ fold change ^A).

Gene Name	Locus Tag	MRSA 11/11 +NaHCO3 (no OXA)	COL +NaHCO3 (no OXA)	MRSA +NaHCO3 with OXA B	COL +NaHCO3 with OXA B	Function
cap8	SAUSA300_0152	−3.05	2.36	−4.45	0.00	virulence
clpL	SAUSA300_2486	0.00	0.00	−1.44	−1.46
sasD	SAUSA300_0136	0.00	0.90	2.11	−1.82
aaa	SAUSA300_0438	0.00	−1.14	−2.20	0.00
vraX	SAUSA300_RS03005	0.00	0.00	−4.45	1.88
kdpABCDF	SAUSA300_2032	0.00	0.90	−2.07	1.79	osmotic stress
betAB	SAUSA300_2545	−4.44	−1.38	−5.12	−4.28
icaR	SAUSA300_2599	0.00	1.87	1.62	0.00	transcriptional regulator
rsp	SAUSA300_2326	−1.60	−1.02	0.00	0.00

^A Fold change comparisons are made by comparing the gene expression for the indicated strain and condition (± OXA) exposed to NaHCO₃ to the same strain/condition without NaHCO₃ exposure. Negative values indicate that exposure to NaHCO₃ for the given strain/condition decreases expression of the indicated gene, whereas positive values indicate that NaHCO₃ increases expression of the indicated gene. ^B OXA concentration equal to 1/2 × MIC in CA-MHB Tris ± NaHCO3 for the indicated strain.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12111650/s1, Table S1: qRT-PCR primers, Table S2: (A) Genes altered ≥ 2-fold by CA-MHB 100 mM Tris + 44 mM NaHCO₃ exposure (without oxacillin) compared to CA-MHB 100 mM Tris alone (without oxacillin) specifically in NaHCO₃-responsive strain MRSA 11/11, (B) Genes altered ≥ 2-fold by CA-MHB 100 mM Tris + 44 mM NaHCO₃ with oxacillin exposure compared to CA-MHB 100 mM Tris alone with oxacillin exposure specifically in NaHCO₃-responsive strain MRSA 11/11, (C) Genes altered ≥ 2-fold by CA-MHB 100 mM Tris + 44 mM NaHCO₃ exposure (± oxacillin) compared to CA-MHB 100 mM Tris alone (± oxacillin) specifically in NaHCO₃-responsive strain MRSA 11/11.

References

1. CDC. Antibiotic Resistance Threats in the United States, 2019; U.S. Deppartment of Health and Human Services, CDC: Atlanta, GA, USA, 2019.

2. Ventola, C.L. The antibiotic resistance crisis: Part 1: Causes and threats. Pharm. Ther.; 2015; 40, 277.

3. Abadi, A.T.B.; Rizvanov, A.A.; Haertlé, T.; Blatt, N.L. World Health Organization report: Current crisis of antibiotic resistance. BioNanoScience; 2019; 9, pp. 778-788. [DOI: https://dx.doi.org/10.1007/s12668-019-00658-4]

4. WHO. Antimicrobial resistance: Global report on surveillance; World Health Organization: Geneva, Switzerland, 2014.

5. Grundmann, H.; Aires-de-Sousa, M.; Boyce, J.; Tiemersma, E. Emergence and resurgence of meticillin-resistant Staphylococcus aureus as a public-health threat. Lancet; 2006; 368, pp. 874-885. [DOI: https://dx.doi.org/10.1016/S0140-6736(06)68853-3]

6. Goetghebeur, M.; Landry, P.-A.; Han, D.; Vicente, C. Methicillin-resistant Staphylococcus aureus: A public health issue with economic consequences. Can. J. Infect. Dis. Med. Microbiol.; 2007; 18, pp. 27-34. [DOI: https://dx.doi.org/10.1155/2007/253947]

7. Hassoun, A.; Linden, P.K.; Friedman, B. Incidence, prevalence, and management of MRSA bacteremia across patient populations—a review of recent developments in MRSA management and treatment. Crit. Care; 2017; 21, 211. [DOI: https://dx.doi.org/10.1186/s13054-017-1801-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28807042]

8. Gould, I.; Reilly, J.; Bunyan, D.; Walker, A. Costs of healthcare-associated methicillin-resistant Staphylococcus aureus and its control. Clin. Microbiol. Infect.; 2010; 16, pp. 1721-1728. [DOI: https://dx.doi.org/10.1111/j.1469-0691.2010.03365.x]

9. Purrello, S.; Garau, J.; Giamarellos, E.; Mazzei, T.; Pea, F.; Soriano, A.; Stefani, S. Methicillin-resistant Staphylococcus aureus infections: A review of the currently available treatment options. J. Glob. Antimicrob. Resist.; 2016; 7, pp. 178-186. [DOI: https://dx.doi.org/10.1016/j.jgar.2016.07.010]

10. Edwards, B.; Andini, R.; Esposito, S.; Grossi, P.; Lew, D.; Mazzei, T.; Novelli, A.; Soriano, A.; Gould, I. Treatment options for methicillin-resistant Staphylococcus aureus (MRSA) infection: Where are we now?. J. Glob. Antimicrob. Resist.; 2014; 2, pp. 133-140. [DOI: https://dx.doi.org/10.1016/j.jgar.2014.03.009]

11. Bamgbola, O. Review of vancomycin-induced renal toxicity: An update. Ther. Adv. Endocrinol. Metab.; 2016; 7, pp. 136-147. [DOI: https://dx.doi.org/10.1177/2042018816638223]

12. Abraham, G.; Finkelberg, D.; Spooner, L.M. Daptomycin-induced acute renal and hepatic toxicity without rhabdomyolysis. Ann. Pharmacother.; 2008; 42, pp. 719-721. [DOI: https://dx.doi.org/10.1345/aph.1K579]

13. McConeghy, K.W.; Bleasdale, S.C.; Rodvold, K.A. The empirical combination of vancomycin and a β-lactam for staphylococcal bacteremia. Clin. Infect. Dis.; 2013; 57, pp. 1760-1765. [DOI: https://dx.doi.org/10.1093/cid/cit560]

14. Kim, S.-H.; Kim, K.-H.; Kim, H.-B.; Kim, N.-J.; Kim, E.-C.; Oh, M.-D.; Choe, K.-W. Outcome of vancomycin treatment in patients with methicillin-susceptible Staphylococcus aureus bacteremia. Antimicrob. Agents Chemother.; 2008; 52, pp. 192-197. [DOI: https://dx.doi.org/10.1128/AAC.00700-07]

15. Ersoy, S.C.; Abdelhady, W.; Li, L.; Chambers, H.F.; Xiong, Y.Q.; Bayer, A.S. Bicarbonate resensitization of methicillin-resistant Staphylococcus aureus to β-lactam antibiotics. Antimicrob. Agents Chemother.; 2019; 63, e00496-19. [DOI: https://dx.doi.org/10.1128/AAC.00496-19] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31010857]

16. Ersoy, S.C.; Otmishi, M.; Milan, V.T.; Li, L.; Pak, Y.; Mediavilla, J.; Chen, L.; Kreiswirth, B.; Chambers, H.F.; Proctor, R.A. et al. Scope and predictive genetic/phenotypic signatures of ‘bicarbonate [NaHCO₃]-responsiveness’ and β-Lactam sensitization among methicillin-resistant Staphylococcus aureus (MRSA). Antimicrob. Agents Chemother.; 2020; 5, e02445-19. [DOI: https://dx.doi.org/10.1128/AAC.02445-19]

17. Rose, W.E.; Bienvenida, A.M.; Xiong, Y.Q.; Chambers, H.F.; Bayer, A.S.; Ersoy, S.C. Ability of bicarbonate supplementation to sensitize selected methicillin-resistant Staphylococcus aureus (MRSA) strains to β-lactam antibiotics in an Ex Vivo simulated endocardial vegetation model. Antimicrob. Agents Chemother.; 2019; 64, e02072-19.

18. Ersoy, S.C.; Chambers, H.F.; Proctor, R.A.; Rosato, A.E.; Mishra, N.N.; Xiong, Y.Q.; Bayer, A.S. Impact of bicarbonate on PBP2a production, maturation, and functionality in methicillin-resistant Staphylococcus aureus. Antimicrob. Agents Chemother.; 2021; 65, e02621-20. [DOI: https://dx.doi.org/10.1128/AAC.02621-20]

19. Ersoy, S.C.; Heithoff, D.M.; Barnes, L.; Tripp, G.K.; House, J.K.; Marth, J.D.; Smith, J.W.; Mahan, M.J. Faculty opinions recommendation of correcting a fundamental flaw in the paradigm for antimicrobial susceptibility testing. EBioMedicine; 2019; 20, pp. 173-181. [DOI: https://dx.doi.org/10.1016/j.ebiom.2017.05.026]

20. Fenn, W.O. The carbon dioxide dissociation curve of nerve and muscle. Am. J. Physiol. Leg. Content; 1928; 85, pp. 207-223. [DOI: https://dx.doi.org/10.1152/ajplegacy.1928.85.2.207]

21. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics; 2014; 30, pp. 2114-2120. [DOI: https://dx.doi.org/10.1093/bioinformatics/btu170] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24695404]

22. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods; 2012; 9, pp. 357-359. [DOI: https://dx.doi.org/10.1038/nmeth.1923]

23. Trapnell, C.; Roberts, A.; Goff, L.; Pertea, G.; Kim, D.; Kelley, D.R.; Pimentel, H.; Salzberg, S.L.; Rinn, J.L.; Pachter, L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc.; 2012; 7, pp. 562-578. [DOI: https://dx.doi.org/10.1038/nprot.2012.016]

24. Bu, D.; Luo, H.; Huo, P.; Wang, Z.; Zhang, S.; He, Z.; Wu, Y.; Zhao, L.; Liu, J.; Guo, J. et al. KOBAS-i: Intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis. Nucleic Acids Res.; 2021; 49, pp. W317-W325. [DOI: https://dx.doi.org/10.1093/nar/gkab447]

25. Sahukhal, G.S.; Pandey, S.; Elasri, M.O. MsaABCR operon is involved in persister cell formation in Staphylococcus aureus. BMC Microbiol.; 2017; 17, 218. [DOI: https://dx.doi.org/10.1186/s12866-017-1129-9]

26. Grassi, L.; Di Luca, M.; Maisetta, G.; Rinaldi, A.C.; Esin, S.; Trampuz, A.; Batoni, G. Generation of persister cells of Pseudomonas aeruginosa and Staphylococcus aureus by chemical treatment and evaluation of their susceptibility to membrane-targeting agents. Front. Microbiol.; 2017; 8, 1917. [DOI: https://dx.doi.org/10.3389/fmicb.2017.01917] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29046671]

27. Lewis, K. Persister cells, dormancy and infectious disease. Nat. Rev. Genet.; 2006; 5, pp. 48-56. [DOI: https://dx.doi.org/10.1038/nrmicro1557]

28. Tuchscherr, L.; Löffler, B.; Proctor, R.A. Persistence of Staphylococcus aureus: Multiple metabolic pathways impact the expression of virulence factors in small-colony variants (SCVs). Front. Microbiol.; 2020; 11, 1028. [DOI: https://dx.doi.org/10.3389/fmicb.2020.01028]

29. Traber, K.E.; Lee, E.; Benson, S.; Corrigan, R.; Cantera, M.; Shopsin, B.; Novick, R.P. Agr function in clinical Staphylococcus aureus isolates. Microbiology; 2008; 154, 2265. [DOI: https://dx.doi.org/10.1099/mic.0.2007/011874-0]

30. Seidl, K.; Chen, L.; Bayer, A.S.; Hady, W.A.; Kreiswirth, B.N.; Xiong, Y.Q. Relationship of agr expression and function with virulence and vancomycin treatment outcomes in experimental endocarditis due to methicillin-resistant Staphylococcus aureus. Antimicrob. Agents Chemother.; 2011; 55, pp. 5631-5639. [DOI: https://dx.doi.org/10.1128/AAC.05251-11]

31. Gupta, R.K.; Luong, T.T.; Lee, C.Y. RNAIII of the Staphylococcus aureus agr system activates global regulator MgrA by stabilizing mRNA. Proc. Natl. Acad. Sci. USA; 2015; 112, pp. 14036-14041. [DOI: https://dx.doi.org/10.1073/pnas.1509251112]

32. Montgomery, C.P.; Boyle-Vavra, S.; Daum, R.S. Importance of the global regulators Agr and SaeRS in the pathogenesis of CA-MRSA USA300 infection. PLoS ONE; 2010; 5, e15177. [DOI: https://dx.doi.org/10.1371/journal.pone.0015177]

33. Cheung, A.L.; Zhang, G. Global regulation of virulence determinants in Staphylococcus aureus by the SarA protein family. Front Biosci.; 2002; 7, pp. d1825-d1842. [DOI: https://dx.doi.org/10.2741/cheung] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12133812]

34. Bischoff, M.; Dunman, P.; Kormanec, J.; Macapagal, D.; Murphy, E.; Mounts, W.; Berger-Bächi, B.; Projan, S. Microarray-based analysis of the Staphylococcus aureus σB regulon. J. Bacteriol.; 2004; 186, pp. 4085-4099. [DOI: https://dx.doi.org/10.1128/JB.186.13.4085-4099.2004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15205410]

35. Van Wamel, W.; Xiong, Y.-Q.; Bayer, A.; Yeaman, M.; Nast, C.; Cheung, A. Regulation of Staphylococcus aureus type 5 capsular polysaccharides by agr and sarA In Vitro and in an experimental endocarditis model. Microb. Pathog.; 2002; 33, pp. 73-79. [DOI: https://dx.doi.org/10.1006/mpat.2002.0513] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12202106]

36. Xiong, Y.-Q.; Bayer, A.S.; Yeaman, M.R.; Van Wamel, W.; Manna, A.C.; Cheung, A.L. Impacts of sarA and agr in Staphylococcus aureus strain newman on fibronectin-binding protein A gene expression and fibronectin adherence capacity In Vitro and in experimental infective endocarditis. Infect. Immun.; 2004; 72, pp. 1832-1836. [DOI: https://dx.doi.org/10.1128/IAI.72.3.1832-1836.2004]

37. Lacey, K.A.; Mulcahy, M.E.; Towell, A.M.; Geoghegan, J.A.; McLoughlin, R.M. Clumping factor B is an important virulence factor during Staphylococcus aureus skin infection and a promising vaccine target. PLoS Pathog.; 2019; 15, e1007713. [DOI: https://dx.doi.org/10.1371/journal.ppat.1007713]

38. McCrea, K.W.; Hartford, O.; Davis, S.; Eidhin, D.N.; Lina, G.; Speziale, P.; Foster, T.J.; Höök, M. The serine-aspartate repeat (Sdr) protein family in Staphylococcus epidermidis. Microbiology; 2000; 146, pp. 1535-1546. [DOI: https://dx.doi.org/10.1099/00221287-146-7-1535]

39. Kong, C.; Chee, C.-F.; Richter, K.; Thomas, N.; Rahman, N.A.; Nathan, S. Suppression of Staphylococcus aureus biofilm formation and virulence by a benzimidazole derivative, UM-C162. Sci. Rep.; 2018; 8, pp. 1-16. [DOI: https://dx.doi.org/10.1038/s41598-018-21141-2]

40. Ersoy, S.C.; Abdelhady, W.; Li, L.; Xiong, Y.Q.; Bayer, A.S. Bicarbonate represses the sarA/mecA regulatory axis, re-sensitizing selected methicillin-resistant Staphylococcus aureus (MRSA) to β-lactam antibiotics. Gordon Research Conference on Staphylococcal Diseases; poster 36 Castelldefels, Spain, 4–9 August 2019.

41. Rahman, M.M.; Hunter, H.N.; Prova, S.; Verma, V.; Qamar, A.; Golemi-Kotra, D. Faculty opinions recommendation of the Staphylococcus aureus methicillin resistance factor FmtA is a d-amino esterase that acts on teichoic acids. mBio; 2016; 7, 02070-15. [DOI: https://dx.doi.org/10.1128/mBio.02070-15]

42. Straume, D.; Piechowiak, K.W.; Kjos, M.; Håvarstein, L.S. Class A PBPs: It is time to rethink traditional paradigms. Mol. Microbiol.; 2021; 116, pp. 41-52. [DOI: https://dx.doi.org/10.1111/mmi.14714]

43. Łęski, T.A.; Tomasz, A. Role of penicillin-binding protein 2 (PBP2) in the antibiotic susceptibility and cell wall cross-linking of Staphylococcus aureus: Evidence for the cooperative functioning of PBP2, PBP4, and PBP2A. J. Bacteriol.; 2005; 187, pp. 1815-1824. [DOI: https://dx.doi.org/10.1128/JB.187.5.1815-1824.2005]

44. Boyle-Vavra, S.; De Jonge, B.; Ebert, C.C.; Daum, R.S. Cloning of the Staphylococcus aureus ddh gene encoding NAD+-dependent d-lactate dehydrogenase and insertional inactivation in a glycopeptide-resistant isolate. J. Bacteriol.; 1997; 179, pp. 6756-6763. [DOI: https://dx.doi.org/10.1128/jb.179.21.6756-6763.1997]

45. Lopes, A.-A.; Yoshii, Y.; Yamada, S.; Nagakura, M.; Kinjo, Y.; Mizunoe, Y.; Okuda, K.-i. Roles of lytic transglycosylases in biofilm formation and β-lactam resistance in methicillin-resistant Staphylococcus aureus. Antimicrob. Agents Chemother.; 2019; 63, e01277-19. [DOI: https://dx.doi.org/10.1128/AAC.01277-19]

46. Stapleton, M.R.; Horsburgh, M.J.; Hayhurst, E.J.; Wright, L.; Jonsson, I.-M.; Tarkowski, A.; Kokai-Kun, J.F.; Mond, J.J.; Foster, S.J. Characterization of IsaA and SceD, two putative lytic transglycosylases of Staphylococcus aureus. J. Bacteriol.; 2007; 189, pp. 7316-7325. [DOI: https://dx.doi.org/10.1128/JB.00734-07]

47. Biswas, R.; Martinez, R.E.; Göhring, N.; Schlag, M.; Josten, M.; Xia, G.; Hegler, F.; Gekeler, C.; Gleske, A.-K.; Götz, F. Proton-binding capacity of Staphylococcus aureus wall teichoic acid and its role in controlling autolysin activity. PLoS ONE; 2012; 7, e41415. [DOI: https://dx.doi.org/10.1371/journal.pone.0041415]

48. Zhao, Y.; Verma, V.; Belcheva, A.; Singh, A.; Fridman, M.; Golemi-Kotra, D. Staphylococcus aureus methicillin-resistance factor fmtA is regulated by the global regulator SarA. PLoS ONE; 2012; 7, e43998. [DOI: https://dx.doi.org/10.1371/journal.pone.0043998]