Next Article in Journal
DNA Satellites Are Transcribed as Part of the Non-Coding Genome in Eukaryotes and Bacteria
Next Article in Special Issue
The Staphylococcus aureus CC398 Lineage: An Evolution Driven by the Acquisition of Prophages and Other Mobile Genetic Elements
Previous Article in Journal
Precision DNA Mixture Interpretation with Single-Cell Profiling
Previous Article in Special Issue
Antibacterial Peptides Resistance in Staphylococcus aureus: Various Mechanisms and the Association with Pathogenicity
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 12
Issue 11
10.3390/genes12111650
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

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

by
Selvi C. Ersoy
1,†,
Blake M. Hanson
2,3,†,
Richard A. Proctor
4,
Cesar A. Arias
2,3,5,
Truc T. Tran
3,5,
Henry F. Chambers
6 and
Arnold S. Bayer
1,7,*
1
The Lundquist Institute, Torrance, CA 90502, USA
2
Center for Infectious Diseases, School of Public Health, University of Texas Health Science Center, Houston, TX 77030, USA
3
Center for Antimicrobial Resistance and Microbial Genomics, Division of Infectious Diseases, University of Texas Health Science Center at Houston McGovern Medical School, Houston, TX 77030, USA
4
School of Medicine and Public Health, University of Wisconsin, Madison, WI 53706-152, USA
5
Houston Methodist Hospital, Houston, TX 77030, USA
6
UCSF School of Medicine, San Francisco, CA 94143, USA
7
The Geffen School of Medicine at UCLA, Los Angeles, CA 90024, USA
*
Author to whom correspondence should be addressed.
†
Authors contributed equally to this work.
Genes 2021, 12(11), 1650; https://doi.org/10.3390/genes12111650
Submission received: 19 August 2021 / Revised: 6 October 2021 / Accepted: 8 October 2021 / Published: 20 October 2021
(This article belongs to the Special Issue Genetics, Genomics and Pathogenesis of Staphylococcus aureus)
Browse Figures
Versions Notes

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 NaHCO3 (termed ‘NaHCO3-responsiveness’). This increased β-lactam susceptibility phenotype has been verified in both ex vivo and in vivo models. Mechanistic studies to-date have implicated NaHCO3-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 NaHCO3-responsive (MRSA 11/11) vs. non-responsive (COL) strains, in the presence vs. absence of NaHCO3-β-lactam co-exposures. These investigations revealed that NaHCO3 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 NaHCO3-OXA exposure in the NaHCO3-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 NaHCO3-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 [NaHCO3]-responsiveness’, wherein certain MRSA strains become sensitized to β-lactam antibiotics in vitro in the presence of NaHCO3 [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 NaHCO3 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 NaHCO3-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 NaHCO3-responsive vs. non-responsive phenotypes. We concentrated on identifying genes which were selectively expressed during NaHCO3-β-lactam co-exposures, comparing a prototype NaHCO3-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 NaHCO3-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 NaHCO3-responsive and -non-responsive strains, respectively [15]. These two prototype strains have been used in a number of our prior investigations related to NaHCO3-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 NaHCO3 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 NaHCO3 to generate RNA to compare the impact of NaHCO3 alone or the impact of NaHCO3 + 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 NaHCO3) (representing ½ the MIC as determined for this strain by broth microdilution in media without and with NaHCO3, respectively [15]); OXA concentrations for COL were 256 µg/mL for both distinct growth conditions. We have used the above NaHCO3 and OXA concentrations in our prior studies [15,16,17,18,19]. Although this NaHCO3 concentration (44 mM) is above human blood levels, it replicates tissue levels [20], and maximizes the NaHCO3-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 NaHCO3 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 OD600 ~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 NaHCO3 are currently available in the NCBI Sequence Read Archive (SRA) database under the accession numbers SAMN21542115 (MRSA 11/11 CA-MHB Tris + NaHCO3), SAMN21542148 (MRSA 11/11 CA-MHB Tris + OXA + NaHCO3), SAMN21542116 (COL CA-MHB Tris + NaHCO3), and SAMN21542147 (COL CA-MHB Tris + OXA + NaHCO3). Codification of the accession numbers for sequence read data-basing for samples that do not contain NaHCO3 (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 NaHCO3±OXA exposure compared to media without NaHCO3 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 (OD600 = 0.5) grown in CA-MHB 100 mM Tris + 2% NaCl + ½ × MIC OXA ± 44 mM NaHCO3 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−ΔΔCT 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 NaHCO3 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 NaHCO3-responsive vs. Non-responsive Strains

RNA-seq analysis was carried out on two prototypical MRSA strains, MRSA 11/11 (NaHCO3-reponsive) and COL (NaHCO3-non-responsive), for two distinct comparison groups: Group 1—CA-MHB Tris vs. CA-MHB Tris + NaHCO3; and Group 2—CA-MHB Tris + 2% NaCl + OXA vs. CA-MHB Tris + NaHCO3 + 2% NaCl + OXA. RNA-seq analysis of the prototypical NaHCO3-responsive strain MRSA 11/11 revealed that exposure to NaHCO3 alone (Group 1) led to the differential change in abundance (≥ 2-fold change) of 21 transcripts (Figure 1A); in contrast, NaHCO3-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.
NaHCO3 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, NaHCO3 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 NaHCO3 exposure, these tended to be related to cellular metabolism, protein translation, and/or stress response, indicating that NaHCO3 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 NaHCO3 may be stimulating a stress response in NaHCO3-responsive MRSA, NaHCO3 also evoked increased expression of many metabolic genes as compared to cells grown in the absence of NaHCO3. This latter finding (i.e., enhanced metabolic activity) would not support the notion of NaHCO3-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. NaHCO3 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 NaHCO3+OXA (n = 14) in the NaHCO3-responsive strain MRSA 11/11 (Figure 2). Additionally, eight transcriptional units were differentially expressed in both growth conditions (NaHCO3 alone and NaHCO3+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 NaHCO3, 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 NaHCO3 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 NaHCO3 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, NaHCO3 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 NaHCO3 exposure in MRSA 11/11 [40]. Taken together, these data indicate a broad impact of NaHCO3 (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 NaHCO3 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 NaHCO3, 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 NaHCO3-mediated susceptibility to β-lactams was its impact on the expression of genes related to cell wall synthesis. NaHCO3 (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 NaHCO3-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 NaHCO3-OXA, but a substantial decrease in sceD transcription in the presence of NaHCO3 (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 NaHCO3, 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 NaHCO3 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 NaHCO3-responsive phenotype. The qRT-PCR also validated the NaHCO3-mediated repression of pbp2, atl, and fnbA in responsive strain MRSA 11/11 as seen by RNA-seq, while NaHCO3 stimulated increased expression of these genes in COL (Figure 3).

3.3. Limitations of Study

This study was limited by only using a single NaHCO3-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 NaHCO3-responsiveness phenotype. Further, only one β-lactam (OXA) was tested for its impact on gene expression in the presence of NaHCO3. Our previous work has demonstrated that NaHCO3-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 NaHCO3-responsive strains, which could be causing the increased frequency of sensitization to CFZ vs. OXA. Additionally, only one concentration of NaHCO3 (44 mM) was examined for its impact on gene expression. Investigating the impact of NaHCO3 at blood concentrations (~25 mM NaHCO3) may reveal genes that are particularly important to NaHCO3-β-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 NaHCO3-responsive strains may help in determining specific, consensus KEGG pathways associated with the responsive phenotype.

4. Conclusions

NaHCO3 (especially in combination with OXA) impacts a variety of genes specifically and differentially expressed in the NaHCO3-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 NaHCO3 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 NaHCO3-responsive phenotype, using strategic gene knockouts. Finally, future studies will be designed to look at the differential impacts of NaHCO3 upon protein-level expression profiles in NaHCO3-responsive vs. non-responsive MRSA.

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 NaHCO3 exposure (without oxacillin) compared to CA-MHB 100 mM Tris alone (without oxacillin) specifically in NaHCO3-responsive strain MRSA 11/11, (B) Genes altered ≥ 2-fold by CA-MHB 100 mM Tris + 44 mM NaHCO3 with oxacillin exposure compared to CA-MHB 100 mM Tris alone with oxacillin exposure specifically in NaHCO3-responsive strain MRSA 11/11, (C) Genes altered ≥ 2-fold by CA-MHB 100 mM Tris + 44 mM NaHCO3 exposure (± oxacillin) compared to CA-MHB 100 mM Tris alone (± oxacillin) specifically in NaHCO3-responsive strain MRSA 11/11.

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.

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. [Google Scholar]
  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, 778–788. [Google Scholar] [CrossRef]
  4. WHO. Antimicrobial resistance: Global report on surveillance; World Health Organization: Geneva, Switzerland, 2014. [Google Scholar]
  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, 874–885. [Google Scholar] [CrossRef] [Green Version]
  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, 27–34. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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, 1721–1728. [Google Scholar] [CrossRef] [Green Version]
  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, 178–186. [Google Scholar] [CrossRef] [Green Version]
  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, 133–140. [Google Scholar] [CrossRef]
  11. Bamgbola, O. Review of vancomycin-induced renal toxicity: An update. Ther. Adv. Endocrinol. Metab. 2016, 7, 136–147. [Google Scholar] [CrossRef]
  12. Abraham, G.; Finkelberg, D.; Spooner, L.M. Daptomycin-induced acute renal and hepatic toxicity without rhabdomyolysis. Ann. Pharmacother. 2008, 42, 719–721. [Google Scholar] [CrossRef]
  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, 1760–1765. [Google Scholar] [CrossRef] [Green Version]
  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, 192–197. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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 [NaHCO3]-responsiveness’ and β-Lactam sensitization among methicillin-resistant Staphylococcus aureus (MRSA). Antimicrob. Agents Chemother. 2020, 5, e02445-19. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar]
  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. [Google Scholar] [CrossRef]
  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, 173–181. [Google Scholar] [CrossRef] [Green Version]
  20. Fenn, W.O. The carbon dioxide dissociation curve of nerve and muscle. Am. J. Physiol. Leg. Content 1928, 85, 207–223. [Google Scholar] [CrossRef] [Green Version]
  21. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [Green Version]
  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, 562–578. [Google Scholar] [CrossRef] [Green Version]
  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, W317–W325. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed]
  27. Lewis, K. Persister cells, dormancy and infectious disease. Nat. Rev. Genet. 2006, 5, 48–56. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [Green Version]
  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, 5631–5639. [Google Scholar] [CrossRef] [Green Version]
  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, 14036–14041. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [Green Version]
  33. Cheung, A.L.; Zhang, G. Global regulation of virulence determinants in Staphylococcus aureus by the SarA protein family. Front Biosci. 2002, 7, d1825–d1842. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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, 4085–4099. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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, 73–79. [Google Scholar] [CrossRef] [PubMed]
  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, 1832–1836. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef]
  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, 1535–1546. [Google Scholar] [CrossRef]
  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, 1–16. [Google Scholar] [CrossRef]
  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. [Google Scholar]
  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. [Google Scholar] [CrossRef] [Green Version]
  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, 41–52. [Google Scholar] [CrossRef]
  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, 1815–1824. [Google Scholar] [CrossRef] [Green Version]
  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, 6756–6763. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef]
  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, 7316–7325. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef]
Genes 12 01650 g001 550
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 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.
Genes 12 01650 g001
Genes 12 01650 g002 550
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 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.
Genes 12 01650 g002
Genes 12 01650 g003 550
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).
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).
Genes 12 01650 g003
Table
Table 1. Overview of gene functional classes altered ≥ 2-fold in two comparison groups, (Group 1) CA-MHB Tris vs. CA-MHB Tris + NaHCO3, and (Group 2) CA-MHB Tris + 2% NaCl + OXA vs. CA-MHB Tris + NaHCO3 + 2% NaCl + OXA. Genes/transcriptional units indicated in this tables are those altered by NaHCO3 exposure alone (Group 1, n = 11), NaHCO3 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 NaHCO3 compared to media without NaHCO3 in each comparison group.
Table 1. Overview of gene functional classes altered ≥ 2-fold in two comparison groups, (Group 1) CA-MHB Tris vs. CA-MHB Tris + NaHCO3, and (Group 2) CA-MHB Tris + 2% NaCl + OXA vs. CA-MHB Tris + NaHCO3 + 2% NaCl + OXA. Genes/transcriptional units indicated in this tables are those altered by NaHCO3 exposure alone (Group 1, n = 11), NaHCO3 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 NaHCO3 compared to media without NaHCO3 in each comparison group.
Classifications of Gene FunctionsNumber of Genes/Transcriptional Units Per CategoryNumber IncreasedNumber Decreased
Virulence12210
Cell Wall Synthesis1019
Membrane19415
Metabolism221210
Transcription514
Translation651
Transposition220
Stress Response211
Totals:782850
Table
Table 2. Genes of interest selectively altered in NaHCO3-responsive strain MRSA 11/11 (log2 fold change A).
Table 2. Genes of interest selectively altered in NaHCO3-responsive strain MRSA 11/11 (log2 fold change A).
Gene NameLocus Tag+NaHCO3 (no OXA)+NaHCO3 with OXA BFunction−202
clfASAUSA300_0772−2.18−2.97virulence
clfBSAUSA300_25651.291.11
agrSAUSA300_1989−3.21−5.37
sdrHSAUSA300_1985−2.440.00
fnbASAUSA300_24410.00−1.84
fnbBSAUSA300_24400.00−1.47
atlSAUSA300_0955−1.65−2.93autolysins
sceDSAUSA300_2051−6.53−4.41
isaASAUSA300_25060.00−2.12
fmtASAUSA300_09590.00−1.32peptidoglycan/cell wall synthesis
ddhSAUSA300_24630.00−4.00
pbp2SAUSA300_13410.00−1.81
bccTSAUSA300_2549−2.45−2.42osmotic stress
uspSAUSA300_0067−1.30−2.28stress response
A Fold change comparisons are made by comparing the gene expression for the indicated strain and condition (± OXA) exposed to NaHCO3 to the same strain/condition without NaHCO3 exposure. Negative values indicate that exposure to NaHCO3 for the given strain/condition decreases expression of the indicated gene, whereas positive values indicate that NaHCO3 increases expression of the indicated gene. B OXA concentration equal to 1/2 × MIC in CA-MHB Tris ± NaHCO3 for the indicated strain.
Table
Table 3. Genes of interest altered in both NaHCO3-responsive and non-responsive strains (log2 fold change A).
Table 3. Genes of interest altered in both NaHCO3-responsive and non-responsive strains (log2 fold change A).
Gene NameLocus TagMRSA 11/11 +NaHCO3 (no OXA)COL +NaHCO3 (no OXA)MRSA +NaHCO3 with OXA BCOL +NaHCO3 with OXA BFunction
cap8SAUSA300_0152−3.052.36−4.450.00virulence
clpLSAUSA300_24860.000.00−1.44−1.46
sasDSAUSA300_01360.000.902.11−1.82
aaaSAUSA300_04380.00−1.14−2.200.00
vraXSAUSA300_RS030050.000.00−4.451.88
kdpABCDFSAUSA300_20320.000.90−2.071.79osmotic stress
betABSAUSA300_2545−4.44−1.38−5.12−4.28
icaRSAUSA300_25990.001.871.620.00transcriptional regulator
rspSAUSA300_2326−1.60−1.020.000.00
A Fold change comparisons are made by comparing the gene expression for the indicated strain and condition (± OXA) exposed to NaHCO3 to the same strain/condition without NaHCO3 exposure. Negative values indicate that exposure to NaHCO3 for the given strain/condition decreases expression of the indicated gene, whereas positive values indicate that NaHCO3 increases expression of the indicated gene. B OXA concentration equal to 1/2 × MIC in CA-MHB Tris ± NaHCO3 for the indicated strain.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ersoy, S.C.; Hanson, B.M.; Proctor, R.A.; Arias, C.A.; Tran, T.T.; Chambers, H.F.; Bayer, A.S. Impact of Bicarbonate-β-Lactam Exposures on Methicillin-Resistant Staphylococcus aureus (MRSA) Gene Expression in Bicarbonate-β-Lactam-Responsive vs. Non-Responsive Strains. Genes 2021, 12, 1650. https://doi.org/10.3390/genes12111650

AMA Style

Ersoy SC, Hanson BM, Proctor RA, Arias CA, Tran TT, Chambers HF, Bayer AS. Impact of Bicarbonate-β-Lactam Exposures on Methicillin-Resistant Staphylococcus aureus (MRSA) Gene Expression in Bicarbonate-β-Lactam-Responsive vs. Non-Responsive Strains. Genes. 2021; 12(11):1650. https://doi.org/10.3390/genes12111650

Chicago/Turabian Style

Ersoy, Selvi C., Blake M. Hanson, Richard A. Proctor, Cesar A. Arias, Truc T. Tran, Henry F. Chambers, and Arnold S. Bayer. 2021. "Impact of Bicarbonate-β-Lactam Exposures on Methicillin-Resistant Staphylococcus aureus (MRSA) Gene Expression in Bicarbonate-β-Lactam-Responsive vs. Non-Responsive Strains" Genes 12, no. 11: 1650. https://doi.org/10.3390/genes12111650

APA Style

Ersoy, S. C., Hanson, B. M., Proctor, R. A., Arias, C. A., Tran, T. T., Chambers, H. F., & Bayer, A. S. (2021). Impact of Bicarbonate-β-Lactam Exposures on Methicillin-Resistant Staphylococcus aureus (MRSA) Gene Expression in Bicarbonate-β-Lactam-Responsive vs. Non-Responsive Strains. Genes, 12(11), 1650. https://doi.org/10.3390/genes12111650

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

Article Metrics

Back to TopTop