PeerJ
PeerJ Inc.
image
Effect of genetic background on the evolution of Vancomycin-Intermediate Staphylococcus aureus (VISA)

DOI: 10.7717/peerj.11764, Volume: 9, Pages: 0-0
Article Type: research-article, Article History
Abstract

Vancomycin-intermediate Staphylococcus aureus (VISA) typically arises through accumulation of chromosomal mutations that alter cell-wall thickness and global regulatory pathways. Genome-based prediction of VISA requires understanding whether strain background influences patterns of mutation that lead to resistance. We used an iterative method to experimentally evolve three important methicillin-resistant S. aureus (MRSA) strain backgrounds—(CC1, CC5 and CC8 (USA300)) to generate a library of 120 laboratory selected VISA isolates. At the endpoint, isolates had vancomycin MICs ranging from 4 to 10 μg/mL. We detected mutations in more than 150 genes, but only six genes (already known to be associated with VISA from prior studies) were mutated in all three background strains (walK, prs, rpoB, rpoC, vraS, yvqF). We found evidence of interactions between loci (e.g., vraS and yvqF mutants were significantly negatively correlated) and rpoB, rpoC, vraS and yvqF were more frequently mutated in one of the backgrounds. Increasing vancomycin resistance was correlated with lower maximal growth rates (a proxy for fitness) regardless of background. However, CC5 VISA isolates had higher MICs with fewer rounds of selection and had lower fitness costs than the CC8 VISA isolates. Using multivariable regression, we found that genes differed in their contribution to overall MIC depending on the background. Overall, these results demonstrated that VISA evolved through mutations in a similar set of loci in all backgrounds, but the effect of mutation in common genes differed with regard to fitness and contribution to resistance in different strains.

Keywords

Introduction

Vancomycin has been a relatively safe and economical drug against MRSA (methicillin-resistant Staphylococcus aureus ) for the past three decades but there have been numerous reports of strains with reduced susceptibility (which we will call “vancomycin-resistant”, using the common terminology) (Gardete & Tomasz, 2014; Zhang et al., 2015; McGuinness, Malachowa & DeLeo, 2017). High-level vancomycin-resistant S. aureus (VRSA; minimum inhibitory concentrations (MIC) ≥ 16 μg/mL) strains have acquired the vanA gene from horizontal transfer from Enterococcus spp (Kobayashi, Musser & DeLeo, 2012). The incidence of VRSA has remained rare, likely due to poor regulation of expression of the vanA gene in S. aureus imposing a fitness burden on the bacterium (Foucault, Courvalin & Grillot-Courvalin, 2009). Vancomycin-intermediate (VISA; MIC 4–8 μg/mL) strains, first reported in 1997 (Hiramatsu et al., 1997), are more commonly encountered in the clinic than VRSA. VISA evolves from vancomycin-susceptible S. aureus (VSSA; MIC ≤ 2 μg/mL) through intrinsic mutagenesis. VISA strains susceptible in vitro to vancomycin (MIC ≤ 2 μg/mL) but containing subpopulations that can grow in the presence of ≥ 4 µg/mL of vancomycin and thus are capable of spontaneous transitions to VISA are termed hVISA (for heterogeneous VISA) (Liu & Chambers, 2003; Sakoulas & Moellering, 2008; Deresinski, 2009; El-Halfawy & Valvano, 2015).

Most intermediate resistance is acquired by within-patient evolution (Klevens et al., 2007; Gardete & Tomasz, 2014). The classic genomics study by Mwangi et al. (2007) showed VISA evolving within a patient on long-term vancomycin therapy through a series of adaptive mutations. VISA has emerged independently in strains from each of the major MRSA lineages (Hiramatsu et al., 1997; Howden et al., 2006, 2008; Klevens et al., 2007; Mwangi et al., 2007; Alam et al., 2014; Gardete & Tomasz, 2014). Some patients with MRSA infections failing therapy showed enhanced vancomycin MIC even when the drug was not used (Horne et al., 2009; Lalueza et al., 2010; Holland & Fowler, 2011), suggesting that VISA may overlap with ‘persister’ phenotypes (Johnson & Levin, 2013) associated with long-term invasive infection.

Standard microbiology laboratory-based phenotypic tests for VISA/hVISA are labor intensive, and the incidence may be higher than currently reported (Marlowe et al., 2001; Charles et al., 2004; Prakash, Lewis & Jorgensen, 2008; Satola et al., 2009; Swenson et al., 2009; Vaudaux et al., 2010). hVISA strains are particularly challenging to detect using current standard clinical microbiology methods (automated broth microdilution assays) as MICs can overlap with VSSA (Swenson et al., 2009). Development of a nucleic acid test using standard PCR based approaches has been confounded to date by the complexity of the genetics of VISA. The genetic basis of VISA has been investigated for several years, primarily through identification of mutations by comparative sequencing of small numbers of isogenic clinical samples (Ohta et al., 2004; Mwangi et al., 2007; Howden et al., 2008) and molecular genetic characterization. VISA cells typically show cell-wall thickening and a reduction of fitness in growth compared to isogenic parents. Mutations in a variety of conserved core S. aureus genes have been reported to be associated with VISA (Howden et al., 2010; Gardete & Tomasz, 2014; Wang et al., 2016), most commonly those involved in regulation of cell wall architecture (e.g., graRS (Neoh et al., 2008; Howden et al., 2008), vraRS (Kato et al., 2010; Baek et al., 2017; Asadpour & Ghazanfari, 2019), yvqF (Kato et al., 2010; Yoo et al., 2013), walK (Shoji et al., 2011; Hafer et al., 2012; McEvoy et al., 2013; Hu et al., 2015), walR (Howden et al., 2011)) or certain global transcriptional regulators (e.g., agr (Sakoulas et al., 2003), stp1 (Cameron et al., 2012; Passalacqua et al., 2012), rpoB (Cui et al., 2010; Watanabe et al., 2011; Hafer et al., 2012; Katayama et al., 2017) and rpoC (Matsuo et al., 2013)). walK controls cell wall autolysis while vraS and yvqF (also known as vraT) regulate cell wall synthesis, thus mutations in these genes contribute to the characteristic thick cell wall of VISA strains.

It is vital for any future gene-based test for VISA to know if there are epistatic interactions between S. aureus strain background and drug resistance mutations. Recent publications have reported important epistatic interactions that drive patterns of antibiotic-resistance (Schubert et al., 2018; Ma et al., 2020). There are also examples of strain-specific effects in S. aureus: clinical CC30 S. aureus strains were found to have elevated (average 100 fold higher) persister formation compared to CC5, CC8, CC30, and CC45 strains (Liu et al., 2020), and under experimental evolution, CC398 strains were found to frequently evolve ciprofloxacin resistance due to a lineage-specific IS element that allowed amplification of the norA efflux pump (Papkou et al., 2020). In this study, we aimed to evaluate the relationship between genetic determinants of intermediate vancomycin resistance and strain background to determine if epistasis can direct the evolutionary trajectory of this phenotype. We evolved VISA strains from three common genetic backgrounds encountered clinically and analyzed the differences in mutation, fitness, and resistance levels.

Materials & methods

Strains

Strains used for the evolution experiments: NRS70 (N315), NRS123 (MW2), and NRS384 (USA300-0114), were obtained from BEI resources (https://www.beiresources.org/).

Experimental evolution of vancomycin resistance

The strategy for the evolution experiment is outlined in Fig. 1. Parent strains were streaked on Brain Heart Infusion (BHI) plates from frozen stocks and single colonies used to establish shaking overnight cultures of BHI broth. We then propagated independent lines in BHI broth shaking at 250 rpm at 37 °C until cultures were turbid (OD600 > 1.5) before diluting cultures at a 1:500 dilution into fresh media. Strains were grown initially on BHI broth containing 1 µg/mL vancomycin then transferred to BHI broth with 2 µg/mL antibiotic. Once turbid, the cultures were similarly diluted to 4 µg/mL, then grown for 3 days before transfer to 8 µg/mL, and grown for a final 3 days. After growth in broth at the final antibiotic concentration, cultures were plated on BHI plates containing 2, 4, and 8 µg/mL of vancomycin. A single colony was picked from the highest concentration at 48 h growth. The colony was passaged twice on selective media and then examined by Gram stain to confirm species.

Experimental evolution of a lineage VISA.
Figure 1
For each of the three S. aureus genetic backgrounds, up to 40 independent lineages were evolved from single colonies isolated from frozen culture streaks. The lineages were cultured in BHI broth containing increasing concentrations of vancomycin. At each stage, single colonies were isolated and frozen and tested for antibiotic MIC using Etest strips (vancomycin and daptomycin) and fitness using a growth curve based assays. Full details in the materials and methods section.Experimental evolution of a lineage VISA.

Genomic analysis

DNA extraction and library prep were performed as manufacturer’s instructions (Wizard Genomic DNA Purification Kit; Promega, Madison, WI, USA; Nextera XT DNA Library Prep Kit; Illumina, San Diego, CA, USA). Genome sequencing was performed on Illumina HiSeq and MiSeq with paired-end reads. Raw read data were deposited in the NCBI Short Read Archive under project accession number PRJNA525705.

Raw read data was processed by BBduk (v37.66) (Bushnell, 2016) to filter out adapters associated with Illumina sequencing and trim reads based on quality. Mutations were called based on comparison to reference sequences using breseq (Deatherage & Barrick, 2014) CONSENSUS mode. Reference genomes used for NRS70, NRS123, and NRS384 VISA strains were NC_002745.2 (BioProject: PRJNA224116, BioSample: SAMD00061099, Assembly: GCF_000009645.1), NC_003923.1 (BioProject: PRJNA224116, BioSample: SAMD00061104, Assembly: GCF_000011265.1), and NZ_CP027476.1 (BioProject: PRJNA224116, BioSample: SAMN07411405, Assembly: GCF_002993865.1) respectively. To compare mutations across the three strains, the protein sequences of NRS123 and NRS384 were aligned using BLAST against that of NRS70 (Camacho et al., 2009). Genes were considered orthologous to NRS70 if the amino acid sequences had at least 95% identity, 95% coverage, or an evalue = 0. To determine if large scale deletions were due to phage loss, PhageWeb was used to identify prophages within the genomes (Sousa et al., 2018).

Measuring fitness based on growth

Strains were grown in duplicate in a 96-well plate beginning at an OD < 0.1 and grown for 24 h at 37 °C with constant shaking in a Biotek Eon Microplate Spectrophotometer, with OD measurements every 10 min. To assess growth curves, OD readings were imported into R (R Core Team, 2016) and maximal growth rate (r) calculated using the growthcurver package (Sprouffske & Wagner, 2016) using the average of two biological replicates. Fitness was calculated as the ratio between the average r of each evolved strain to the average r of the parent strain. Fitness distributions between NRS70 evolved VISA strains and NRS384 evolved VISA strains were compared using a two-sided two-sample Kolmogorov-Smirnov test.

Statistical analysis

To determine if the difference in the prevalence of mutations across different VISA genes (walK, prs, rpoB, rpoC, vraS, yvqF) was significant, a Pearson’s Chi-squared test was performed, with subsequent post hoc pairwise comparisons using Fisher’s exact test and a Bonferroni multiple test correction. To analyze if the prevalence of mutations was also affected by background, binomial generalized linear models (GLMs) were fitted for each gene, and pairwise Wald tests were performed to test whether the difference in the coefficients for each background was non-zero. The statistical significance of gene correlation was assessed using a Spearman correlation with a Holm multiple testing correction. Total SNP distributions were compared using a Kruskal-Wallis test followed by pairwise Mann-Whitney U tests. Synonymous and nonsynonymous SNP distributions were compared with a Pearson’s Chi-squared test followed by Bonferroni corrected Fisher’s exact tests. The vancomycin and daptomycin MIC distributions between NRS70 VISA strains and NRS384 VISA strains were compared using two-sided two-sample Kolmogorov–Smirnov tests. To test if the linear relationship between vancomycin MIC and daptomycin MIC was dependent on background (NRS70 and NRS384), an Analysis of Covariance (ANCOVA) was performed.

To assess the effect of mutations in individual genes on vancomycin MIC, linear regression models were fitted using the lm function in R. walK, prs, rpoB, rpoC, vraS, and yvqF as binary variables (mutated or not) were used as predictors for log transformed vancomycin MICs. All mutations found in these genes were classified as nonsynonymous SNPs. The final model was chosen by backwards selection with the goal of minimizing Akaike information criterion (AIC) as in Eyre et al. (2017).

To investigate differences in the effect sizes of SNPs on vancomycin MIC from background, standard linear models were fitted with interaction terms for SNP presence and background. In these models, strain background was included as a fixed-effects predictor because the mixed-effects version (with background as a vector random effect affecting the main and interaction effects of the SNPs) was not computationally stable. To study potential pairwise interactions between SNPs affecting the MIC of vancomycin, linear mixed-effects models were fit with presence or absence of SNPs as fixed-effects predictors (main and interaction effects) and with background (NRS70, NRS384) as a scalar random effect. Mixed-effects models with vector random effects (affecting intercept and slope of certain interactions) were also tested, but these models were not stable enough for convergence. Given the sparsity of the data matrix, the predictors were restricted to SNPs present in at least 10% of the strains. Predictors included were the six universally mutated genes (walK, prs, rpoB, rpoC, vraS, yvqF) and sdrC. All SNPs with the exception of SNPs in sdrC, which were all synonymous, were nonsynonymous. The response variable (vancomycin MIC) was log transformed. These models were also fitted for another response (log transformed growth rate as a proxy of fitness). Similar to the vancomycin models, pairwise SNP interactions were assessed with mixed-effects models and SNP-background interactions with standard linear models. A Bonferroni corrected p-value (ɑ = 0.00056) was used to assess test significance. R scripts used for these analyses are available at https://github.com/crsl4/staph-visa.

Antimicrobial susceptibility testing

Antimicrobial susceptibility testing for daptomycin, vancomycin, and methicillin were performed according to Clinical Laboratory Standards Institute (CLSI) standards (CLSI, 2013). In brief, cultures were streaked from a frozen stock onto BHI agar and restreaked the next day. From the second day culture, colonies were resuspended in normal saline to a 0.5 McFarland standard. For plate-based testing, Etest for daptomycin and disc diffusion with cefoxitin, cultures were struck to create a lawn on cation-adjusted Mueller-Hinton broth (CAMHB) agar plates before application. Plates were incubated at 35 °C and read at 18 h. For broth microdilution (BMD) to determine strain vancomycin minimum inhibitory concentration (MIC), 0.5 McFarland standards were diluted 1:20 in normal saline, and 10 µL was added to 90 µL of CAMHB with the appropriate concentration of vancomycin. Plates were incubated at 35 °C without shaking and read at 24 h.

Results

Patterns of acquisition of SNPs at VISA endpoints is influenced by genetic background

Three strains, representing major clinical lineages in the US: N315/NRS70 (Kuroda et al., 2001) (ST5/CC5), USA300-0114/NRS384 (Diep et al., 2006) (ST8/CC8) and MW2/NRS123 (Baba et al., 2002) (ST1/CC1) were experimentally evolved under selection pressure to achieve a vancomycin MIC of 4-10 μg/mL from an initial MIC of 1 μg/mL (Fig. 1). Selection of 120 of these mutants, 40 independent lines per genetic background, was achieved by sequentially raising the vancomycin concentration through serial passages of strains in BHI broth. Mutants were sequenced and compared to their isogenic parent using breseq (Deatherage & Barrick, 2014) to ascertain mutations acquired during vancomycin selection. Thirteen evolved strains are excluded from further analyses due to either low sequence quality or culture contamination. The final number of evolved VISA strains per background with high quality DNA sequences was 35 NRS70 strains, 34 NRS123 strains, and 38 NRS384 strains.

The majority of the genetic changes that occurred during the selection for vancomycin-intermediate resistance were single nucleotide polymorphisms (SNPs) and will be the focus of the rest of the analyses. Most deletions, insertions, and substitutions only occurred once (Tables S1S6). Of note, prophage loss (StauST398-4 in NRS70, Pvl108 in NRS123, 23MRA in NRS384) was observed in four strains from each of the genetic backgrounds. Phage loss has been previously observed in vancomycin and daptomycin resistant isolates and has been associated with the induction of the SOS response by antibiotics (Maiques et al., 2006; Boyle-Vavra et al., 2011; Machado et al., 2020).

We found between 5 and 25 SNPs in evolved VISA strains compared to their respective parent strains, including non-synonymous mutations in several genes previously reported to be associated with VISA genes. As the genome sizes for NRS70, NRS123 and NRS384 are 2.81 MB, 2.82 MB, and 2.87 MB, respectively, we did not expect differences in genome size to significantly influence the number of mutations observed during the course of selection.

SNPs fell within 151 coding genes, and no gene was universally mutated in all strains, highlighting the existence of multiple pathways to intermediate vancomycin resistance in S. aureus. However, of the 151 genes, 121 genes had mutations in only one of the three backgrounds; 24 genes were mutated in two backgrounds; and only six genes (walK, prs, rpoB, rpoC, vraS, yvqF) had non-synonymous mutations in all three backgrounds. The six “universally mutated” genes have been previously implicated in VISA from sequencing of clinical isolates, but interestingly, some genes that had been previously reported did not acquire mutations frequently or in all three backgrounds. Notably, we did not find any agrR mutants, and graRS mutants were found in only NRS384 derived strains and walR (Howden, Peleg & Stinear, 2013) only in NRS70 and NRS123 derived strains.

The six universally mutated genes varied in their frequency across strains. On opposite sides of the spectrum, mutations in prs occurred in fewer than 10% of the evolved strains while mutations in walK occurred in 80–95% of the evolved VISA strains (Table 1). A Pearson’s chi-squared test (p < 2.2 × 10−16) and post hoc Fisher’s exact tests (all p < 5.91 × 10−03) determined that the proportion of strains differed significantly between each of these six genes. There was also evidence that some gene mutations occurred more frequently in some backgrounds than others. By fitting binomial generalized linear models (GLMs) to the individual gene mutation distributions with the genetic background as a predictor, proportions of strains with rpoB (NRS70 vs. NRS123, p = 0.03; NRS70 vs. NRS384, p = 0.015), rpoC (NRS70 vs. NRS384, p = 0.019), vraS (NRS70 vs. NRS384, p = 0.04), and yvqF (NRS70 vs. NRS123, p = 0.035; NRS70 vs. NRS384, p = 0.028) were found to be significantly influenced by genetic background. Between NRS70 and NRS123, five other genes were mutated in both backgrounds. Similarly, between NRS123 and NRS384, six other genes were mutated. Finally, between NRS70 and NRS384, 12 other genes were mutated (Table 2). The higher proportion of other shared genes between NRS70 and NRS384 may indicate that these two strains share more similar evolutionary paths toward vancomycin resistance, but mutations in these genes were rare, and differences were not statistically significant.

table-wrap
Table 1
Mutations in these genes were called by breseq. The NRS70 walK prevalence calculation considers the presence of SNPs as well as inframe deletions. All other prevalence calculations are based solely on SNPs.
Prevalence of mutations in six universally mutated genes by genetic background.
Gene/Background NRS70 NRS123 NRS384
walK 80% 85.3% 94.7%
prs 5.7% 5.9% 2.6%
rpoB 45.7% 20.6% 18.4%
rpoC 34.3% 20.6% 10.5%
vraS 17.1% 23.6% 39.5%
yvqF 80% 55.9% 55.3%
table-wrap
Table 2
Genes with SNPs shared between two genetic backgrounds.
NRS70 and NRS123 NRS123 and NRS384 NRS70 and NRS384
Genes shared in at least 1 evolved isolate WP_000101976.1 (walR)
WP_000035320.1 (pdhA)
WP_001032833.1 (gad)
WP_000375864.1
(HTH cro/C1-type domain-containing protein)
WP_000251253.1 (UPF0374 protein)
WP_001081640.1
(Cyclic-di-AMP phosphodiesterase) WP_000260117.1 (pdhD)
WP_001071136.1 (mprF)
WP_000809131.1 (cspA)
WP_000048060.1 (rpsU)
WP_000830380.1 (mgt)
WP_000631969.1 (cysS)
WP_001060462.1 (sdrC)
WP_001273060.1 (tagG) WP_000491755.1 (ykaA)
WP_000120368.1 (pitA)
WP_000431312.1 (greA)
WP_000782121.1 (prsA)
WP_000153535.1 (vraR)
WP_000159960.1 (pyrG)
WP_000347896.1 (dacA)
WP_000004085.1 (rpsQ)
WP_000008673.1 (ureC)

We examined the co-occurrence of mutations in the six universally mutated genes (Figs. 2A2D). Mutations in vraS and yvqF were almost mutually exclusive (coefficient −0.8, p < 0.0001), with minor differences between genetic backgrounds (Figs. 2B2D, NRS70 and NRS384: −0.9; NRS123: −0.6). NRS70 was the only background where we detected significant associations between other genes (Fig. 2B): prs/vraS (coefficient 0.5, p = 0.011) and prs/yvqF (coefficient 0.5, p = 0.035).

Correlation of mutations between VISA genes.
Figure 2
The numbers are correlation coefficients between genes based on the presence of non-synonymous SNPs or frameshift deletions in the same strain (A) Correlation matrix of all evolved VISA strains. (B) Correlation matrix of all NRS70 evolved VISA strains. (C) Correlation matrix of all NRS123 evolved VISA strains. (D) Correlation matrix of all NRS384 evolved VISA strains. Significance codes: ‘***’ <0.0001 ‘**’ <0.001 ‘*’ <0.05.Correlation of mutations between VISA genes.

To examine the distribution of amino acid changes caused by SNPs, we focused on yvqF and walK due to the high frequency of mutation in our VISA strains (Table 1). For yvqF, there were 32 SNPs in 68 strains: 3 SNPs are found in all three backgrounds, 10 SNPs are found in two, and 19 SNPs in only one. Despite walK mutations (41 SNPs across 93 strains) being more common, there are no SNPs in common between the three backgrounds, and 31 SNPs are found in only one. Furthermore, the most common SNPs within each background for these two genes were found in different amino acid residues (yvqF: NRS70/A152V, NRS123/P126S, NRS384/P174A; walK: NRS70/E236D, NRS123/I29M, NRS384/Q216E).

In population genetics, there is a common a priori assumption that synonymous SNPs have little or no functional impact compared to non-synonymous. However, we noted some genes with high rates of synonymous mutations. Strikingly, sdrC was mutated in 7 out of 38 NRS384 VISA strains. The sdr locus encodes sdrC, sdrD, sdrE which are members of the repeat-rich microbial surface components recognizing adhesive matrix molecules (MSCRAMM) family (Josefsson et al., 1998). These genes are not always conserved together (Liu et al., 2015), and their presence has been associated with bone infection (Sabat et al., 2006), resistance to host immunity (Sitkiewicz, Babiak & Hryniewicz, 2011; Askarian et al., 2017), biofilm formation (Barbu et al., 2014), and host-cell adhesion (Corrigan, Miajlovic & Foster, 2009; Cheng et al., 2009; Barbu et al., 2010; Askarian et al., 2016). Belikova et al. (2020) showed that in USA300 sdr genes undergo frequent within-genome recombination during growth in lab conditions and during infection. The high number of synonymous mutations may be the result of repeats expanding or collapsing in evolved strains, causing base-calling artefacts when mapped onto the reference genomes. Ultimately, long-read sequencing approaches are needed to deconvolute these complex structures, but there is a strong possibility that base-calling artefacts inflated synonymous mutation numbers in some strains. We also noted that there were significant differences between backgrounds in the total number of mutations that accrued during the evolution experiments (p = 1.1 × 10−5). NRS123 VISA (median 3 mutations) strains had fewer mutations than NRS70 (median 4, p = 0.004) and NRS384 (median 5, p = 2.8 × 10−5 ) VISA strains (Fig. 3). Furthermore, the number of synonymous SNPs relative to nonsynonymous SNPs was not the same (Table 3, Pearson chi square p = 1.8 × 10−5, NRS70 vs. NRS123 p = 0.046, NRS70 vs. NRS384 p = 0.038, NRS123 vs. NRS384 p = 1.0 × 10−7). However, when excluding sdr mutations, NRS70 and NRS123 are no longer statistically significantly different, but NRS384 remained statistically different from NRS70 and NRS123 (Table 3, Pearson chi square p = 0.00021, NRS70 vs. NRS384 p = 0.014, NRS123 vs. NRS384 p = 0.0011). Furthermore, the number of synonymous SNPs relative to nonsynonymous SNPs was not the same (Table 3, Pearson chi square p = 1.8 × 10−5, NRS70 vs. NRS123 p = 0.046, NRS70 vs. NRS384 p = 0.038, NRS123 vs. NRS384 p = 1.0 × 10−7).

SNP distributions of VISA evolved strains differ by genetic background.
Figure 3
Mutations were called with breseq. Duplicate sdr mutations have been removed. Mode and median of SNPs per background: NRS70 (4), NRS123 (3), NRS384 (5). Significance codes: An asterisk (*) designates p < 0.06.SNP distributions of VISA evolved strains differ by genetic background.
table-wrap
Table 3
Synonymous SNPs and nonsynonymous SNPs by genetic background.
Synonymous Nonsynonymous
NRS70 19 (10) 149
NRS123 4 114
NRS384 45 (32) 166
Note:
Numbers given are the sum of the SNPs of all strains in that background (NRS70: 35 strains, NRS123: 34 strains, NRS384: 38 strains). The number of mutations when excluding sdr operon mutations is given in parenthesis.

Resistance phenotypes of evolved VISA strains

All evolved VISA strains at the endpoint of the iterative selection were able to grow in 8 μg/mL vancomycin BHI broth. However, to accurately measure the level of resistance, we performed BMD on NRS70 and NRS384 strains according to CLSI standards (unfortunately, all NRS123 frozen cultures were accidentally discarded while in −80 °C storage before we could test them). All endpoint NRS70 and NRS384 strains had MICs within the range to be considered VISA (Fig. 4A), but some (20 NRS70 and 28 NRS384 strains) were found to be below 8 μg/mL. The NRS70 VISA strain MICs ranged from 4 to 10 μg/mL, with a median and mode of 6 μg/mL, while the NRS384 VISA strain MICs ranged from 4 to 8 μg/mL, with a median and mode of 6 μg/mL. The difference in the MIC distributions was not statistically significant between NRS70 and NRS384 strains (p = 0.323). Mutated genes and/or mutated gene patterns associated with isolates with the highest level of vancomycin resistance were not appreciably different than those found in isolates with a lower resistance, likely indicating that individual mutations played a large role in the resistance achieved versus genes overall.

Antibiotic MIC distributions of NRS70 and NRS384 evolved strains.
Figure 4
Vancomycin MICs were determined by BMD and daptomycin MICs by Etest. (A) Vancomycin MIC distributions. The comparison of the two distributions by a Kolmogorov–Smirnov test was not statistically significant (p = 0.3236). Parental MICs are indicated by a red triangle at 1 μg/ml. (B) Daptomycin MIC distributions. The comparison of the two distributions by a Kolmogorov–Smirnov test was not statistically significant (p = 0.9996). Parental MICs are indicated by a red triangle at 0.38 ug/ml.Antibiotic MIC distributions of NRS70 and NRS384 evolved strains.

The cross-resistance of vancomycin and daptomycin in VISA strains has been well-studied (Cui et al., 2006; Nam et al., 2018). NRS70 and NRS384 endpoint strains were thus tested for development of daptomycin resistance by Etest. Parental strains had daptomycin MICs of 0.38 μg/mL. Thirty of the 35 NRS70 strains became resistant, and 33 of the 38 NRS384 strains became resistant (Fig. 4B). As observed with the vancomycin resistance genes, there was no clear mutational pattern associated with higher daptomycin resistance. The daptomycin distributions between NRS70 and NRS384 were not statistically different (p = 1), however, only NRS384 strains had a daptomycin MIC of 6 μg/mL versus 4 μg/mL for NRS70 strains. Greater levels of vancomycin resistance correlated to greater levels of daptomycin resistance (Fig. 5, p = 0.006). However, the regressions for each background were not significantly different (p = 0.822).

Linear relationship between vancomycin MIC and daptomycin MIC.
Figure 5
Vancomycin MICs were determined by BMD and daptomycin MICs by Etest. Linear regression of daptomycin MIC with vancomycin MIC was statistically significant (p = 0.006). The comparison of the background specific linear relationships between vancomycin and daptomycin by an Analysis of Covariance (ANCOVA) was not statistically significant (p = 0.822). Overall: Adjusted R-squared: 0.08903 (p = 0.005); N70: Adjusted R-squared: 0.1866 (p = 0.005); N384: Adjusted R-squared: 0.03534 (p = 0.13).Linear relationship between vancomycin MIC and daptomycin MIC.

Collateral sensitivity with beta-lactams is a less frequently reported phenomenon in VISA strains (Sieradzki & Tomasz, 1999; Mwangi et al., 2007; Barber et al., 2014). Within our evolved NRS70 and NRS384 strains, only two NRS70 strains became methicillin sensitive. Known mutations for beta-lactam sensitivity in a MRSA background such as loss of SCCmec or inactivation of mecA were not found, and there the mechanism for collateral sensitivity in these strains was unknown.

Fitness decreases during evolution of higher resistance but is modulated by genetic background

We examined the fitness of intermediate and endpoint strains of each NRS70 and NRS384 experimentally evolved lineage using maximum growth rate in batch culture log(r) as an indicator. Isolates from later stages of selection, able to tolerate greater concentration of vancomycin, generally equated to lower fitness (Figs. 6A6B). Looking at the fitness distributions at each selection stage (4 μg/mL, 6 μg/mL, 8 μg/mL) between NRS70 and NRS384, there was a striking difference at 4 μg/mL (p = 7.08 × 10−6) and 6 μg/mL (p = 1.86 × 10−6), which disappeared at 8 μg/mL (p = 0.273) (Figs. 7A7C). Thus, NRS70 strains appeared to be able to tolerate lower levels of vancomycin resistance with lower fitness cost compared to NRS384 strains. Linear mixed models with SNP presence/absence as fixed effects and background as a scalar random effect were used to determine interaction effects between individual SNPs and log(r). However, no results were statistically significant after Bonferroni correction. The strongest interaction was between sdrC and rpoC whereby mutations in sdrC had a negative effect on growth rate in the presence of rpoC mutation from a previous positive effect (uncorrected p = 0.023, Fig. S1), but as noted above, the caveat is that sdrC mutations may be surrogates for more complex recombination events.

Fitness of NRS70 and NRS384 evolved strains.
Figure 6
Strains were grown in duplicate in a 96-well plate beginning at an OD < 0.1 and grown for 24 h at 37 °C with constant shaking in a Biotek Eon Microplate Spectrophotometer, with OD measurements every 10 min. To assess growth-curves, OD readings were imported into R and maximal growth rate (r) calculated using the growthcurver package. Fitness was calculated as the ratio between the average r of each evolved strain to the average r of the parent strain. (a) Fitness (r) of NRS70 evolved strains plotted against the concentration of vancomycin that the strains were evolved. NRS70-7 at vancomycin selection stage 6 and 8 and NRS70-14 at vancomycin selection stage 6 were excluded due to lack of growth or too large a variance between replicates. (b) Fitness (r) progression of NRS384 evolved strains.Fitness of NRS70 and NRS384 evolved strains.
Fitness distributions at each stage of vancomycin selection for NRS70 and NRS384 strains.
Figure 7
MICs were determined by BMD. Strains were grown in duplicate in a 96-well plate beginning at an OD < 0.1 and grown for 24 h at 37 °C with constant shaking in a Biotek Eon Microplate Spectrophotometer, with OD measurements every 10 min. To assess growth-curves, OD readings were imported into R and maximal growth rate (r) calculated using the growthcurver package. Fitness was calculated as the ratio between the average r of each evolved strain to the average r of the parent strain. Fitness distributions between NRS70 evolved VISA strains and NRS384 evolved VISA strains were compared using a two-sided two-sample Kolmogorov–Smirnov test. (A) Vancomycin selection stage 4 μg/mL. (B) Vancomycin selection stage 6 μg/mL. (C) Vancomycin selection stage 8 μg/mL. NRS70-7 at vancomycin selection stage 6 and 8 and NRS70-14 at vancomycin selection stage 6 were excluded due to lack of growth or too large a variance between replicates.Fitness distributions at each stage of vancomycin selection for NRS70 and NRS384 strains.

Genetic predictors of vancomycin MIC

To investigate the interactions of mutation and S. aureus strain background to vancomycin resistance, we used two modelling approaches. In the first, we used multiple linear regression models to determine which mutated genes were most predictive of vancomycin resistance. To reduce the number of parameters and for ease of comparison, only the six universally mutated genes were used in the initial models (Table 4). To make the models more meaningful, parameters were evaluated for inclusion in the final models using AIC (Eyre et al., 2017). In the final model for NRS70, walK (coefficient estimate 0.298), rpoC (0.228), vraS (0.989), and yvqF (1.058) were statistically significantly associated with vancomycin resistance, defined as log(MIC) (R2 = 0.583, p = 2.15 × 10−6). In comparison, only walK (coefficient estimate 0.633) by itself was statistically significant in the final model of NRS384 (R2 = 0.313, p = 9.95 × 10−4), although vraS (0.340) and yvqF (0.378) were included for improved model fit.

table-wrap
Table 4
An initial model was fitted using walK, prs, rpoB, rpoC, vraS, and yvqF as binary (mutated or not) predictors for log transformed vancomycin MICs. The predictors in the final were chosen by backwards selection.
Linear regression models of vancomycin MIC.
Background Gene Estimate (initial) p (initial) Estimate (final) p (final)
NRS70 walK
prs
rpoB
rpoC
vraS
yvqF
0.288
0.198
0.132
0.309
0.901
0.956
0.015
0.389
0.254
0.012
0.0005
<10−4
0.298


0.228
0.989
1.058
0.012


0.026
<10−4
<10−5
NRS384 walK
prs
rpoB
rpoC
vraS
yvqF
0.629
0.255
0.103
–0.050
0.412
0.450
0.004
0.515
0.460
0.790
0.129
0.073
0.633



0.340
0.378
0.002



0.096
0.053

In the second approach, we fit linear regression models to determine if the effect sizes of SNPs on log(MIC) differed by background. No results were statistically significant after Bonferroni correction, but there was a suggestive result that rpoC may contribute positively to vancomycin MIC in NRS70 strains but not NRS384 strains (uncorrected p = 0.028, Fig. S2), which is consistent with our previous linear regression models. Linear mixed models were used to determine interaction effects between SNPs. No results were statistically significant after Bonferroni correction, but there was a suggestive result that the effect of a sdrC mutation (either increasing or decreasing vancomycin MIC) depended on walK (uncorrected p = 0.0472, Fig. S3). Overall, the models indicated that while the six universally mutated core genes were broadly responsible for vancomycin resistance in NRS70 and NRS384, they were of differential importance within each background. In addition, we did not find strong evidence for pairwise interaction between SNPs influencing MIC.

Discussion

In this project we looked for evidence of strain specificity in patterns of genetic change responsible for transition to vancomycin intermediate resistance in S. aureus . We found evidence that mutations in the same gene had differential effects on strain fitness and MIC depending on the genetic background. There were also statistically significant correlations in the patterns of co-occurrence of mutations linked to VISA. We did not probe the mechanisms of the putative epistatic interactions in this study. It is possible that, since many VISA-linked genes affect global cellular regulatory pathways, there is cross-talk with fixed mutations in the strain background that set global state changes (Priest et al., 2012). The strains also differ in their accessory gene content (Lindsay et al., 2006; McCarthy & Lindsay, 2010; Aanensen et al., 2016), which may also affect gene expression and protein interactions responsible for the VISA phenotype. This work adds to a growing body of literature on the importance of strain background in this pathogen. S. aureus lineages have been shown to be heterogeneous in several significant clinically relevant phenotypes such as toxin production, biofilm formation, and host immune resistance (King et al., 2016; Su et al., 2020). Epistatic interactions between antibiotic resistance and toxicity have also been explored (Yokoyama et al., 2018).

VISA was found to be linked to mutations (primarily SNPs) in a limited number of genes (e.g., walKR, rpoB/C, vraTSR). Six previously implicated VISA-associated genes (walK, prs, rpoB, rpoC, vraS, yvqF ) (Howden et al., 2010) were mutated in our three genetic backgrounds. walK followed by yvqF were the most commonly seen mutations. For rpoB, rpoC, vraS, and yvqF , NRS70 and NRS384 mutation incidence were most different, and NRS123 mutation incidence was intermediate between the others. These results may be partially explained by the closer genetic relationship between CC1 and CC8 than CC5 and CC8 (Petit & Read, 2018). Other genes were found to have been mutated in two of the three genetic backgrounds and may have a role in vancomycin resistance or act as compensatory mutations. Several are already validated as being associated with vancomycin resistance (walR (Howden et al., 2011), vraR (Kato et al., 2010; Baek et al., 2017; Asadpour & Ghazanfari, 2019), mprF (Ruzin et al., 2003; Chen et al., 2018)). NRS70 and NRS384 backgrounds had more of these genes in common than with NRS123, indicating that NRS70 and NRS384 strains may experience more similar selective pressures in the larger mutational landscape outside of the previously discussed six genes.

The overall vancomycin resistance levels achieved after experimental evolution in NRS70 and NRS384 VISA strains were similar, suggesting both backgrounds are equally capable of achieving lower and higher levels of vancomycin intermediate resistance. Linear regression demonstrated that genetic predictors of vancomycin resistance differed between backgrounds. walK, rpoC, vraS, and yvqF were all significantly associated with vancomycin MICs in NRS70, with vraS and yvqF contributing most significantly, but only walK was significant in NRS384. As NRS384 strains carried more mutations than the other two backgrounds, this may indicate that while a subset of the six universally mutated genes (walK, prs, rpoB, rpoC, vraS, yvqF) are necessary for vancomycin resistance, other genes may play a significant role in determining the level of vancomycin resistance achieved in NRS384 VISA strains but not NRS70 strains.

We found that vancomycin-intermediate resistance imposed a fitness cost (measured by maximal growth rate compared to parent strain) to S. aureus that linearly scaled with the level of vancomycin resistance in most cases. Fitness distributions between these two genetic backgrounds were markedly different at stages 4 μg/mL and 6 μg/mL, with vancomycin resistance imposing less of a fitness cost on NRS70 VISA strains than NRS384 VISA strains. This may have clinical implications as it suggests that NRS70 S. aureus strains may evolve vancomycin resistance quickly, not be rapidly selected against due to growth defects, and may reach fixation if vancomycin concentrations are sustained and lead to the extinction of the ancestral vancomycin susceptible strain.

The laboratory evolution approach used in this study gave us the power to evolve VISA multiple times independently and achieve statistical significance using general linear models. With these data, we were able to investigate the effect of different SNPs with the strain background as a fixed effect. However, we cannot draw strong conclusions about clinical relevance as clinical VISA strains are the result of within-host evolution, and immune pressure is not present within our experimental design. Other factors that are present in human infection that were omitted in our experiment include but are not limited to nutrient limiting conditions and microenvironments. These undoubtedly reduce the mutational space S. aureus has available for adaption. In addition, the mutations in our studies were grouped by gene for simplicity, but likely not all mutations in VISA genes contribute equally to vancomycin resistance, and some may have no effect. The sheer number of mutations makes it impossible to recapitulate these mutations individually. We did not run control evolution experiments without vancomycin, so we cannot say whether some of the mutations in the most common six gene were due to random genetic drift, although it is well established that VISA strains do not evolve from VSSA in the laboratory in the absence of antibiotic selection. Another study design limitation was that fitness determinations by measuring maximal growth rate do not always correlate with competitive indices of evolved strains to ancestor strains and have been shown to be sensitive to environmental conditions. Thus, our measured fitness values may not accurately reflect within-host fitness, especially in environments with competing microbes.

The results of this study offer some encouragement that development of VISA strains may be predicted with reasonable accuracy directly from the genome sequence, an important future technique in diagnostic clinical microbiology (Su, Satola & Read, 2019). This is because the most common VISA mutations occurred in all three backgrounds, suggesting a “dictionary” of common mutations could be compiled. However, the finding that individual gene contribution to MIC may be background dependent (Table 4) suggests that prediction of the level of resistance in strains from different clonal complexes may be a much harder problem. In order to obtain the data to make exhaustive catalogs of VISA mutations and their epistatic dependencies or to build an accurate classifier, much more extensive versions of the experiment approach outlined here need to be performed, with a greater number of representative genetic backgrounds.

Conclusions

We found that there was a complex relationship between genetic background and VISA. There was clear evidence for parallel evolution of VISA through mutations in a common set of genes across strains (especially: walK, prs, rpoB, rpoC, vraS, yvqF). Some of these genes (particularly yvqF and vraS) showed negative epistasis in all strains. However, we found evidence the effects of mutations on vancomycin MIC and the relationship between fitness and MIC were dependent on strain background. These results indicate that prediction of the levels of vancomycin resistance based on genome sequence may require extensive databases of mutants from different genetic backgrounds.

Acknowledgements

We would like to thank Jon Moller and Robert Petit for discussion and comments on the manuscript as well as Sujith Swarna for help with the experimental evolution experiments.

Additional Information and Declarations

Competing Interests

Timothy D. Read is an Academic Editor for PeerJ.

Author Contributions

Michelle Su conceived and designed the experiments, performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the paper, and approved the final draft.
Michelle H. Davis performed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.
Jessica Peterson performed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.
Claudia Solis-Lemus analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the paper, and approved the final draft.
Sarah W. Satola analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.
Timothy D. Read conceived and designed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.

DNA Deposition

The following information was supplied regarding the deposition of DNA sequences:

Raw read data are available in the NCBI Short Read Archive: PRJNA525705.

Data Availability

The following information was supplied regarding data availability:

The raw data are available in the Supplemental Files. Antimicrobial susceptibility testing results for VISA strains are given for vancomycin, daptomycin, and cefoxitin. These data were used for analyses regarding MIC distribution and fitness comparisons. The raw data used for fitness/growth rate (r) analyses is given: OD measurements every 10 min for 24 h.

References

Aanensen et al. (2016) 

    Aanensen   DM , Feil   EJ , Holden   MTG , Dordel   J , Yeats   CA , Fedosejev   A , Goater   R , Castillo-Ramírez   S , Corander   J , Colijn   C , Chlebowicz   MA , Schouls   L , Heck   M , Pluister   G , Ruimy   R , Kahlmeter   G , Åhman   J , Matuschek   E , Friedrich   AW , Parkhill   J , Bentley   SD , Spratt   BG , Grundmann   H , European SRL Working Group . Whole-genome sequencing for routine pathogen surveillance in public health: a population snapshot of invasive Staphylococcus aureus in Europe. mBio 2016. 7: 3 e00444-16 doi: 10.1128/mBio.00444-16

Alam et al. (2014) 

    Alam   MT , Petit   RA , Crispell   EK , Thornton   TA , Conneely   KN , Jiang   Y , Satola   SW , Read   TD . Dissecting vancomycin intermediate resistance in Staphylococcus aureus using genome-wide association. Genome Biology and Evolution 2014. 6: 5 1174-1185 doi: 10.1093/gbe/evu092

Asadpour & Ghazanfari (2019) 

    Asadpour   L , Ghazanfari   N . Detection of vancomycin nonsusceptible strains in clinical isolates of Staphylococcus aureus in northern Iran. International Microbiology: the Official Journal of the Spanish Society for Microbiology 2019. 22: 4 411-417 doi: 10.1007/s10123-019-00063-7

Askarian et al. (2016) 

    Askarian   F , Ajayi   C , Hanssen   A-M , van Sorge   NM , Pettersen   I , Diep   DB , Sollid   JUE , Johannessen   M . The interaction between Staphylococcus aureus SdrD and desmoglein 1 is important for adhesion to host cells. Scientific Reports 2016. 6: 1 22134 doi: 10.1038/srep22134

Askarian et al. (2017) 

    Askarian   F , Uchiyama   S , Valderrama   JA , Ajayi   C , Sollid   JUE , van Sorge   NM , Nizet   V , van Strijp   JAG , Johannessen   M . Serine-aspartate repeat protein D increases staphylococcus aureus virulence and survival in blood. Infection and Immunity 2017. 85: 1 e00559-16 doi: 10.1128/IAI.00559-16

Baba et al. (2002) 

    Baba   T , Takeuchi   F , Kuroda   M , Yuzawa   H , Aoki   K-I , Oguchi   A , Nagai   Y , Iwama   N , Asano   K , Naimi   T , Kuroda   H , Cui   L , Yamamoto   K , Hiramatsu   K . Genome and virulence determinants of high virulence community-acquired MRSA. The Lancet 2002. 359: 9320 1819-1827 doi: 10.1016/S0140-6736(02)08713-5

Baek et al. (2017) 

    Baek   JY , Chung   DR , Ko   KS , Kim   SH , Yang   S-J , Kang   C-I , Peck   KR , Song   J-H . Genetic alterations responsible for reduced susceptibility to vancomycin in community-associated MRSA strains of ST72. The Journal of Antimicrobial Chemotherapy 2017. 72: 9 2454-2460 doi: 10.1093/jac/dkx175

Barber et al. (2014) 

    Barber   KE , Ireland   CE , Bukavyn   N , Rybak   MJ . Observation of seesaw effect with vancomycin, teicoplanin, daptomycin and ceftaroline in 150 unique MRSA strains. Infectious Diseases and Therapy 2014. 3: 1 35-43 doi: 10.1007/s40121-014-0023-0

Barbu et al. (2010) 

    Barbu   EM , Ganesh   VK , Gurusiddappa   S , Mackenzie   RC , Foster   TJ , Sudhof   TC , Höök   M . beta-Neurexin is a ligand for the Staphylococcus aureus MSCRAMM SdrC. PLOS Pathogens 2010. 6: 1e1000726 doi: 10.1371/journal.ppat.1000726

Barbu et al. (2014) 

    Barbu   EM , Mackenzie   C , Foster   TJ , Höök   M . SdrC induces staphylococcal biofilm formation through a homophilic interaction. Molecular Microbiology 2014. 94: 1 172-185 doi: 10.1111/mmi.12750

Belikova et al. (2020) 

    Belikova   D , Jochim   A , Power   J , Holden   MTG , Heilbronner   S . Gene accordions cause genotypic and phenotypic heterogeneity in clonal populations of Staphylococcus aureus. Nature Communications 2020. 11: 1 3526 doi: 10.1038/s41467-020-17277-3

Boyle-Vavra et al. (2011) 

    Boyle-Vavra   S , Jones   M , Gourley   BL , Holmes   M , Ruf   R , Balsam   AR , Boulware   DR , Kline   S , Jawahir   S , Devries   A , Peterson   SN , Daum   RS . Comparative genome sequencing of an isogenic pair of USA800 clinical methicillin-resistant Staphylococcus aureus isolates obtained before and after daptomycin treatment failure. Antimicrobial Agents and Chemotherapy 2011. 55: 5 2018-2025 doi: 10.1128/AAC.01593-10

Bushnell (2016) 

Camacho et al. (2009) 

    Camacho   C , Coulouris   G , Avagyan   V , Ma   N , Papadopoulos   J , Bealer   K , Madden   TL . BLAST: architecture and applications. BMC Bioinformatics 2009. 10: 1 421 doi: 10.1186/1471-2105-10-421

Cameron et al. (2012) 

    Cameron   DR , Ward   DV , Kostoulias   X , Howden   BP , Moellering   RC , Eliopoulos   GM , Peleg   AY . Serine/threonine phosphatase Stp1 contributes to reduced susceptibility to vancomycin and virulence in Staphylococcus aureus. The Journal of Infectious Diseases 2012. 205: 11 1677-1687 doi: 10.1093/infdis/jis252

Charles et al. (2004) 

    Charles   PGP , Ward   PB , Johnson   PDR , Howden   BP , Grayson   ML . Clinical features associated with bacteremia due to heterogeneous vancomycin-intermediate Staphylococcus aureus. 2004. Clinical Infectious Diseases: An Official Publication of the Infectious Diseases Society of America 38: 448-451 doi: 10.1086/381093

Cheng et al. (2009) 

    Cheng   AG , Kim   HK , Burts   ML , Krausz   T , Schneewind   O , Missiakas   DM . Genetic requirements for Staphylococcus aureus abscess formation and persistence in host tissues. FASEB Journal: Official Publication of the Federation of American Societies for Experimental Biology 2009. 23: 10 3393-3404 doi: 10.1096/fj.09-135467

Chen et al. (2018) 

    Chen   F-J , Lauderdale   T-L , Lee   C-H , Hsu   Y-C , Huang   I-W , Hsu   P-C , Yang   C-S . Effect of a point mutation in mprF on susceptibility to daptomycin, vancomycin, and oxacillin in an MRSA clinical strain. Frontiers in Microbiology 2018. 9: 1086 doi: 10.3389/fmicb.2018.01086

CLSI (2013) 

    CLSI . Performance standards for antimicrobial susceptibility testing: CLSI approved standard M100-S23. 2013. Wayne CLSI

Corrigan, Miajlovic & Foster (2009) 

    Corrigan   RM , Miajlovic   H , Foster   TJ . Surface proteins that promote adherence of Staphylococcus aureus to human desquamated nasal epithelial cells. BMC Microbiology 2009. 9: 1 22 doi: 10.1186/1471-2180-9-22

Cui et al. (2010) 

    Cui   L , Isii   T , Fukuda   M , Ochiai   T , Neoh   H-M , Camargo   ILBDC , Watanabe   Y , Shoji   M , Hishinuma   T , Hiramatsu   K . An RpoB mutation confers dual heteroresistance to daptomycin and vancomycin in Staphylococcus aureus. Antimicrobial Agents and Chemotherapy 2010. 54: 12 5222-5233 doi: 10.1128/AAC.00437-10

Cui et al. (2006) 

    Cui   L , Tominaga   E , Neoh   H-M , Hiramatsu   K . Correlation between reduced daptomycin susceptibility and vancomycin resistance in vancomycin-intermediate Staphylococcus aureus. Antimicrobial Agents and Chemotherapy 2006. 50: 3 1079-1082 doi: 10.1128/AAC.50.3.1079-1082.2006

Deatherage & Barrick (2014) 

    Deatherage   DE , Barrick   JE . Identification of mutations in laboratory-evolved microbes from next-generation sequencing data using breseq. Methods in Molecular Biology 2014. 1151: 165-188 doi: 10.1007/978-1-4939-0554-6

Deresinski (2009) 

    Deresinski   S . Vancomycin heteroresistance and methicillin-resistant Staphylococcus aureus. The Journal of Infectious Diseases 2009. 199: 5 605-609 doi: 10.1086/596630

Diep et al. (2006) 

    Diep   BA , Gill   SR , Chang   RF , Phan   TH , Chen   JH , Davidson   MG , Lin   F , Lin   J , Carleton   HA , Mongodin   EF , Sensabaugh   GF , Perdreau-Remington   F . Complete genome sequence of USA300, an epidemic clone of community-acquired meticillin-resistant Staphylococcus aureus. The Lancet 2006. 367: 9512 731-739 doi: 10.1016/S0140-6736(06)68231-7

El-Halfawy & Valvano (2015) 

    El-Halfawy   OM , Valvano   MA . Antimicrobial heteroresistance: an emerging field in need of clarity. Clinical Microbiology Reviews 2015. 28: 1 191-207 doi: 10.1128/CMR.00058-14

Eyre et al. (2017) 

    Eyre   DW , De Silva   D , Cole   K , Peters   J , Cole   MJ , Grad   YH , Demczuk   W , Martin   I , Mulvey   MR , Crook   DW , Walker   AS , Peto   TEA , Paul   J . WGS to predict antibiotic MICs for Neisseria gonorrhoeae. The Journal of Antimicrobial Chemotherapy 2017. 72: 7 1937-1947 doi: 10.1093/jac/dkx067

Foucault, Courvalin & Grillot-Courvalin (2009) 

    Foucault   M-L , Courvalin   P , Grillot-Courvalin   C . Fitness cost of VanA-type vancomycin resistance in methicillin-resistant Staphylococcus aureus. Antimicrobial Agents and Chemotherapy 2009. 53: 6 2354-2359 doi: 10.1128/AAC.01702-08

Gardete & Tomasz (2014) 

    Gardete   S , Tomasz   A . Mechanisms of vancomycin resistance in Staphylococcus aureus. The Journal of Clinical Investigation 2014. 124: 7 2836-2840 doi: 10.1172/JCI68834

Hafer et al. (2012) 

    Hafer   C , Lin   Y , Kornblum   J , Lowy   FD , Uhlemann   A-C . Vancomycin-intermediate Staphylococcus aureus: the contribution of selected gene mutations in clinical isolates. Antimicrobial Agents and Chemotherapy 2012. 27: 11 5845-5851 doi: 10.1128/AAC.01139-12

Hiramatsu et al. (1997) 

    Hiramatsu   K , Hanaki   H , Ino   T , Yabuta   K , Oguri   T , Tenover   FC . Methicillin-resistant Staphylococcus aureus clinical strain with reduced vancomycin susceptibility. Journal of Antimicrobial Chemotherapy 1997. 40: 1 135-136 doi: 10.1093/jac/40.1.135

Holland & Fowler (2011) 

    Holland   TL , Fowler   VG . Vancomycin minimum inhibitory concentration and outcome in patients with Staphylococcus aureus bacteremia: pearl or pellet?. Journal of Infectious Diseases 2011. 204: 3 329-331 doi: 10.1093/infdis/jir275

Horne et al. (2009) 

    Horne   KC , Howden   BP , Grabsch   EA , Graham   M , Ward   PB , Xie   S , Mayall   BC , Johnson   PDR , Grayson   ML . Prospective comparison of the clinical impacts of heterogeneous vancomycin-intermediate methicillin-resistant Staphylococcus aureus (MRSA) and vancomycin-susceptible MRSA. Antimicrobial Agents and Chemotherapy 2009. 53: 8 3447-3452 doi: 10.1128/AAC.01365-08

Howden et al. (2010) 

    Howden   BP , Davies   JK , Johnson   PDR , Stinear   TP , Grayson   ML . Reduced vancomycin susceptibility in staphylococcus aureus, including vancomycin-intermediate and heterogeneous vancomycin-intermediate strains: resistance mechanisms, laboratory detection, and clinical implications. Clinical Microbiology Reviews 2010. 23: 1 99-139 doi: 10.1128/CMR.00042-09

Howden et al. (2006) 

    Howden   BP , Johnson   PDR , Ward   PB , Stinear   TP , Davies   JK . Isolates with low-level vancomycin resistance associated with persistent methicillin-resistant Staphylococcus aureus bacteremia. Antimicrobial Agents and Chemotherapy 2006. 50: 9 3039-3047 doi: 10.1128/AAC.00422-06

Howden et al. (2011) 

    Howden   BP , McEvoy   CRE , Allen   DL , Chua   K , Gao   W , Harrison   PF , Bell   J , Coombs   G , Bennett-Wood   V , Porter   JL , Robins-Browne   R , Davies   JK , Seemann   T , Stinear   TP . Evolution of multidrug resistance during staphylococcus aureus infection involves mutation of the essential two component regulator WalKR. PLOS Pathogens 2011. 7: 11e1002359 doi: 10.1371/journal.ppat.1002359

Howden, Peleg & Stinear (2013) 

    Howden   BP , Peleg   AY , Stinear   TP . The evolution of vancomycin intermediate Staphylococcus aureus (VISA) and heterogenous-VISA. Infection, Genetics and Evolution: Journal of Molecular Epidemiology and Evolutionary Genetics in Infectious Diseases 2013. 21: 1-9 doi: 10.1016/j.meegid.2013.03.047

Howden et al. (2008) 

    Howden   BP , Stinear   TP , Allen   DL , Johnson   PDR , Ward   PB , Davies   JK . Genomic analysis reveals a point mutation in the two-component sensor gene graS that leads to intermediate vancomycin resistance in clinical Staphylococcus aureus. Antimicrobial Agents and Chemotherapy 2008. 52: 10 3755-3762 doi: 10.1128/AAC.01613-07

Hu et al. (2015) 

    Hu   J , Zhang   X , Liu   X , Chen   C , Sun   B . Mechanism of reduced vancomycin susceptibility conferred by walK mutation in community-acquired methicillin-resistant Staphylococcus aureus Strain MW2. Antimicrobial Agents and Chemotherapy 2015. 59: 2 1352-1355 doi: 10.1128/AAC.04290-14

Johnson & Levin (2013) 

    Johnson   PJT , Levin   BR . Pharmacodynamics, population dynamics, and the evolution of persistence in Staphylococcus aureus. PLOS Genetics 2013. 9: 1e1003123 doi: 10.1371/journal.pgen.1003123

Josefsson et al. (1998) 

    Josefsson   E , McCrea   KW , Eidhin   DN , O’Connell   D , Cox   J , Höök   M , Foster   TJ . Three new members of the serine-aspartate repeat protein multigene family of Staphylococcus aureus. Microbiology 1998. 144: 12 3387-3395 doi: 10.1099/00221287-144-12-3387

Katayama et al. (2017) 

    Katayama   Y , Azechi   T , Miyazaki   M , Takata   T , Sekine   M , Matsui   H , Hanaki   H , Yahara   K , Sasano   H , Asakura   K , Takaku   T , Ochiai   T , Komatsu   N , Chambers   HF . Prevalence of slow-growth vancomycin nonsusceptibility in methicillin-resistant Staphylococcus aureus. Antimicrobial Agents and Chemotherapy 2017. 61: 11 AAC.00452-17 doi: 10.1128/AAC.00452-17

Kato et al. (2010) 

    Kato   Y , Suzuki   T , Ida   T , Maebashi   K . Genetic changes associated with glycopeptide resistance in Staphylococcus aureus: predominance of amino acid substitutions in YvqF/VraSR. Journal of Antimicrobial Chemotherapy 2010. 65: 1 37-45 doi: 10.1093/jac/dkp394

King et al. (2016) 

    King   JM , Kulhankova   K , Stach   CS , Vu   BG , Salgado-Pabón   W . Phenotypes and virulence among Staphylococcus aureus USA100, USA200, USA300, USA400, and USA600 clonal lineages. mSphere 2016. 1: 3 e00071-16

Klevens et al. (2007) 

    Klevens   RM , Morrison   MA , Nadle   J , Petit   S , Gershman   K , Ray   S , Harrison   LH , Lynfield   R , Dumyati   G , Townes   JM , Craig   AS , Zell   ER , Fosheim   GE , McDougal   LK , Carey   RB , Fridkin   SK . Active Bacterial Core surveillance (ABCs) MRSA investigators, invasive methicillin-resistant Staphylococcus aureus infections in the United States. JAMA: The Journal of the American Medical Association 2007. 298: 15 1763-1771 doi: 10.1001/jama.298.15.1763

Kobayashi, Musser & DeLeo (2012) 

    Kobayashi   SD , Musser   JM , DeLeo   FR . Genomic analysis of the emergence of vancomycin-resistant Staphylococcus aureus. mBio 2012. 3: 4 e00170-12 doi: 10.1128/mBio.00170-12

Kuroda et al. (2001) 

    Kuroda   M , Ohta   T , Uchiyama   I , Baba   T , Yuzawa   H , Kobayashi   I , Cui   L , Oguchi   A , Aoki   K , Nagai   Y , Lian   J , Ito   T , Kanamori   M , Matsumaru   H , Maruyama   A , Murakami   H , Hosoyama   A , Mizutani-Ui   Y , Takahashi   NK , Sawano   T , Inoue   R , Kaito   C , Sekimizu   K , Hirakawa   H , Kuhara   S , Goto   S , Yabuzaki   J , Kanehisa   M , Yamashita   A , Oshima   K , Furuya   K , Yoshino   C , Shiba   T , Hattori   M , Ogasawara   N , Hayashi   H , Hiramatsu   K . Whole genome sequencing of meticillin-resistant Staphylococcus aureus. The Lancet 2001. 357: 9264 1225-1240 doi: 10.1016/S0140-6736(00)04403-2

Lalueza et al. (2010) 

    Lalueza   A , Chaves   F , San Juan   R , Daskalaki   M , Otero   JR , Aguado   JM . Is high vancomycin minimum inhibitory concentration a good marker to predict the outcome of methicillin-resistant Staphylococcus aureus bacteremia?. The Journal of Infectious Diseases 2010. 201: 311-312 doi: 10.1086/649572

Lindsay et al. (2006) 

    Lindsay   JA , Moore   CE , Day   NP , Peacock   SJ , Witney   AA , Stabler   RA , Husain   SE , Butcher   PD , Hinds   J . Microarrays reveal that each of the ten dominant lineages of Staphylococcus aureus has a unique combination of surface-associated and regulatory genes. Journal of Bacteriology 2006. 188: 2 669-676 doi: 10.1128/JB.188.2.669-676.2006

Liu & Chambers (2003) 

    Liu   C , Chambers   HF . Staphylococcus aureus with heterogeneous resistance to vancomycin: epidemiology, clinical significance, and critical assessment of diagnostic methods. Antimicrobial Agents and Chemotherapy 2003. 47: 3040-3045

Liu et al. (2015) 

    Liu   H , Lv   J , Qi   X , Ding   Y , Li   D , Hu   L , Wang   L , Yu   F . The carriage of the serine-aspartate repeat protein-encoding sdr genes among Staphylococcus aureus lineages. The Brazilian Journal of Infectious Diseases: An Official Publication of the Brazilian Society of Infectious Diseases 2015. 19: 5 498-502 doi: 10.1016/j.bjid.2015.07.003

Liu et al. (2020) 

    Liu   L , Wang   Y , Bojer   MS , Andersen   PS , Ingmer   H . High persister cell formation by clinical Staphylococcus aureus strains belonging to clonal complex 30. Microbiology 2020. 166: 7 654-658 doi: 10.1099/mic.0.000926

Machado et al. (2020) 

    Machado   H , Seif   Y , Sakoulas   G , Olson   CA , Szubin   R , Palsson   BO , Nizet   V , Feist   AM . Environmental conditions dictate differential evolution of vancomycin resistance in Staphylococcus aureus. BioRxiv 2020. doi: 10.1101/2020.06.07.138933

Maiques et al. (2006) 

    Maiques   E , Ubeda   C , Campoy   S , Salvador   N , Lasa   I , Novick   RP , Barbé   J , Penadés   JR . beta-lactam antibiotics induce the SOS response and horizontal transfer of virulence factors in Staphylococcus aureus. Journal of Bacteriology 2006. 188: 7 2726-2729 doi: 10.1128/JB.188.7.2726-2729.2006

Ma et al. (2020) 

    Ma   KC , Mortimer   TD , Hicks   AL , Wheeler   NE , Sánchez-Busó   L , Golparian   D , Taiaroa   G , Rubin   DHF , Wang   Y , Williamson   DA , Unemo   M , Harris   SR , Grad   YH . Increased antibiotic susceptibility in Neisseria gonorrhoeae through adaptation to the cervical environment. BioRxiv 2020. doi: 10.1038/s41467-020-17980-1

Marlowe et al. (2001) 

    Marlowe   EM , Cohen   MD , Hindler   JF , Ward   KW , Bruckner   DA . Practical strategies for detecting and confirming vancomycin-intermediate Staphylococcus aureus: a tertiary-care hospital laboratory’s experience. Journal of Clinical Microbiology 2001. 39: 7 2637-2639 doi: 10.1128/JCM.39.7.2637-2639.2001

Matsuo et al. (2013) 

    Matsuo   M , Cui   L , Kim   J , Hiramatsu   K . Comprehensive identification of mutations responsible for heterogeneous vancomycin-intermediate Staphylococcus aureus (hVISA)-to-VISA conversion in laboratory-generated VISA strains derived from hVISA clinical strain Mu3. Antimicrobial Agents and Chemotherapy 2013. 57: 12 5843-5853 doi: 10.1128/AAC.00425-13

McCarthy & Lindsay (2010) 

    McCarthy   AJ , Lindsay   JA . Genetic variation in Staphylococcus aureus surface and immune evasion genes is lineage associated: implications for vaccine design and host-pathogen interactions. BMC Microbiology 2010. 10: 1 173 doi: 10.1186/1471-2180-10-173

McEvoy et al. (2013) 

    McEvoy   CRE , Tsuji   B , Gao   W , Seemann   T , Porter   JL , Doig   K , Ngo   D , Howden   BP , Stinear   TP . Decreased vancomycin susceptibility in Staphylococcus aureus caused by IS256 tempering of WalKR expression. Antimicrobial Agents and Chemotherapy 2013. 57: 7 3240-3249 doi: 10.1128/AAC.00279-13

McGuinness, Malachowa & DeLeo (2017) 

    McGuinness   WA , Malachowa   N , DeLeo   FR . Vancomycin resistance in Staphylococcus aureus. The Yale Journal of Biology and Medicine 2017. 90: 269-281

Mwangi et al. (2007) 

    Mwangi   MM , Wu   SW , Zhou   Y , Sieradzki   K , de Lencastre   H , Richardson   P , Bruce   D , Rubin   E , Myers   E , Siggia   ED , Tomasz   A . Tracking the in vivo evolution of multidrug resistance in Staphylococcus aureus by whole-genome sequencing. Proceedings of the National Academy of Sciences of the United States of America 2007. 104: 22 9451-9456 doi: 10.1073/pnas.0609839104

Nam et al. (2018) 

    Nam   EY , Yang   S-J , Kim   ES , Cho   JE , Park   K-H , Jung   S-I , Yoon   N , Kim   D-M , Lee   C-S , Jang   H-C , Park   Y , Lee   KS , Kwak   YG , Lee   JH , Park   SY , Hwang   J-H , Kim   M , Song   K-H , Kim   HB . Emergence of Daptomycin-nonsusceptible methicillin-resistant staphylococcus aureus clinical isolates among daptomycin-naive patients in Korea. Microbial Drug Resistance 2018. 24: 5 534-541 doi: 10.1089/mdr.2017.0212

Neoh et al. (2008) 

    Neoh   H-M , Cui   L , Yuzawa   H , Takeuchi   F , Matsuo   M , Hiramatsu   K . Mutated response regulator graR is responsible for phenotypic conversion of Staphylococcus aureus from heterogeneous vancomycin-intermediate resistance to vancomycin-intermediate resistance. Antimicrobial Agents and Chemotherapy 2008. 52: 1 45-53 doi: 10.1128/AAC.00534-07

Ohta et al. (2004) 

    Ohta   T , Hirakawa   H , Morikawa   K , Maruyama   A , Inose   Y , Yamashita   A , Oshima   K , Kuroda   M , Hattori   M , Hiramatsu   K , Kuhara   S , Hayashi   H . Nucleotide substitutions in Staphylococcus aureus strains, Mu50, Mu3, and N315. DNA Research: An International Journal for Rapid Publication of Reports on Genes and Genomes 2004. 11: 1 51-56 doi: 10.1093/dnares/11.1.51

Papkou et al. (2020) 

    Papkou   A , Hedge   J , Kapel   N , Young   B , Craig MacLean   R . Efflux pump activity potentiates the evolution of antibiotic resistance across S. aureus isolates. Nature Communications 2020. 11: 1 1-15 doi: 10.1038/s41467-020-17735-y

Passalacqua et al. (2012) 

    Passalacqua   KD , Satola   SW , Crispell   EK , Read   TD . A mutation in the PP2C phosphatase gene in a Staphylococcus aureus USA300 clinical isolate with reduced susceptibility to vancomycin and daptomycin. Antimicrobial Agents and Chemotherapy 2012. 56: 10 5212-5223 doi: 10.1128/AAC.05770-11

Petit & Read (2018) 

    Petit   RA , Read   TD . Staphylococcus aureus viewed from the perspective of 40,000+ genomes. PeerJ 2018. 6: 5e5261 doi: 10.7717/peerj.5261

Prakash, Lewis & Jorgensen (2008) 

    Prakash   V , Lewis   JS , Jorgensen   JH . Vancomycin MICs for methicillin-resistant Staphylococcus aureus isolates differ based upon the susceptibility test method used. Antimicrobial Agents and Chemotherapy 2008. 52: 12 4528 doi: 10.1128/AAC.00904-08

Priest et al. (2012) 

    Priest   NK , Rudkin   JK , Feil   EJ , van den Elsen   JMH , Cheung   A , Peacock   SJ , Laabei   M , Lucks   DA , Recker   M , Massey   RC . From genotype to phenotype: can systems biology be used to predict Staphylococcus aureus virulence?. Nature Reviews Microbiology 2012. 10: 11 791-797 doi: 10.1038/nrmicro2880

R Core Team (2016) 

    R Core Team R: A language and environment for statistical computing 2016. Vienna, Austria The R Foundation for Statistical Computing

Ruzin et al. (2003) 

    Ruzin   A , Severin   A , Moghazeh   SL , Etienne   J , Bradford   PA , Projan   SJ , Shlaes   DM . Inactivation of mprF affects vancomycin susceptibility in Staphylococcus aureus. Biochimica et Biophysica Acta 2003. 1621: 2 117-121 doi: 10.1016/S0304-4165(03)00028-X

Sabat et al. (2006) 

    Sabat   A , Melles   DC , Martirosian   G , Grundmann   H , van Belkum   A , Hryniewicz   W . Distribution of the serine-aspartate repeat protein-encoding sdr genes among nasal-carriage and invasive Staphylococcus aureus strains. Journal of Clinical Microbiology 2006. 44: 3 1135-1138 doi: 10.1128/JCM.44.3.1135-1138.2006

Sakoulas et al. (2003) 

    Sakoulas   G , Eliopoulos   GM , Moellering   RC , Novick   RP , Venkataraman   L , Wennersten   C , DeGirolami   PC , Schwaber   MJ , Gold   HS . Staphylococcus aureus accessory gene regulator (agr) group II: is there a relationship to the development of intermediate-level glycopeptide resistance?. The Journal of Infectious Diseases 2003. 187: 6 929-938 doi: 10.1086/368128

Sakoulas & Moellering (2008) 

    Sakoulas   G , Moellering   RC . Increasing antibiotic resistance among methicillin-resistant Staphylococcus aureus strains. Clinical Infectious Diseases: An Official Publication of the Infectious Diseases Society of America 2008. 5: S360-S367 doi: 10.1086/533592

Satola et al. (2009) 

    Satola   SW , Caliendo   AM , Farley   MM , Patel   JB , Burd   EM . Lack of heteroresistance among Staphylococcus aureus isolates with vancomycin MICs of 2 micrograms per milliliter by automated testing. Journal of Clinical Microbiology 2009. 47: 8 2680-2681 doi: 10.1128/JCM.01184-09

Schubert et al. (2018) 

    Schubert   B , Maddamsetti   R , Nyman   J , Farhat   MR , Marks   DS . Genome-wide discovery of epistatic loci affecting antibiotic resistance in Neisseria gonorrhoeae using evolutionary couplings. Nature Microbiology 2018. 4: 2 328-338 doi: 10.1038/s41564-018-0309-1

Shoji et al. (2011) 

    Shoji   M , Cui   L , Iizuka   R , Komoto   A , Neoh   H-M , Watanabe   Y , Hishinuma   T , Hiramatsu   K . walK and clpP mutations confer reduced vancomycin susceptibility in Staphylococcus aureus. Antimicrobial Agents and Chemotherapy 2011. 55: 8 3870-3881 doi: 10.1128/AAC.01563-10

Sieradzki & Tomasz (1999) 

    Sieradzki   K , Tomasz   A . Gradual alterations in cell wall structure and metabolism in vancomycin-resistant mutants of Staphylococcus aureus. Journal of Bacteriology 1999. 181: 24 7566-7570 doi: 10.1128/JB.181.24.7566-7570.1999

Sitkiewicz, Babiak & Hryniewicz (2011) 

    Sitkiewicz   I , Babiak   I , Hryniewicz   W . Characterization of transcription within sdr region of Staphylococcus aureus. Antonie van Leeuwenhoek 2011. 99: 2 409-416 doi: 10.1007/s10482-010-9476-7

Sousa et al. (2018) 

    Sousa   ALD , Maués   D , Lobato   A , Franco   EF , Pinheiro   K , Araújo   F , Pantoja   Y , Costa da Silva   ALD , Morais   J , Ramos   RTJ . PhageWeb—web interface for rapid identification and characterization of prophages in bacterial genomes. Frontiers in Genetics 2018. 9: W16 doi: 10.3389/fgene.2018.00644

Sprouffske & Wagner (2016) 

    Sprouffske   K , Wagner   A . Growthcurver: an R package for obtaining interpretable metrics from microbial growth curves. BMC Bioinformatics 2016. 17: 1 172 doi: 10.1186/s12859-016-1016-7

Su et al. (2020) 

    Su   M , Lyles   JT , Petit III   RA , Peterson   J , Hargita   M , Tang   H , Solis-Lemus   C , Quave   CL , Read   TD . Genomic analysis of variability in Delta-toxin levels between Staphylococcus aureus strains. PeerJ 2020. 8: 22e8717 doi: 10.7717/peerj.8717

Su, Satola & Read (2019) 

    Su   M , Satola   SW , Read   TD . Genome-based prediction of bacterial antibiotic resistance. Journal of Clinical Microbiology 2019. 57: 3 e01405-18 doi: 10.1128/JCM.01405-18

Swenson et al. (2009) 

    Swenson   JM , Anderson   KF , Lonsway   DR , Thompson   A , McAllister   SK , Limbago   BM , Carey   RB , Tenover   FC , Patel   JB . Accuracy of commercial and reference susceptibility testing methods for detecting Vancomycin-Intermediate Staphylococcus aureus. Journal of Clinical Microbiology 2009. 47: 7 2013-2017 doi: 10.1128/JCM.00221-09

Vaudaux et al. (2010) 

    Vaudaux   P , Huggler   E , Bernard   L , Ferry   T , Renzoni   A , Lew   DP . Underestimation of Vancomycin and Teicoplanin MICs by broth microdilution leads to under-detection of glycopeptide-intermediate isolates of Staphylococcus aureus. Antimicrobial Agents and Chemotherapy 2010. 54: 9 3861-3870 doi: 10.1128/AAC.00269-10

Wang et al. (2016) 

    Wang   Y , Li   X , Jiang   L , Han   W , Xie   X , Jin   Y , He   X , Wu   R . Novel mutation sites in the development of vancomycin—intermediate resistance in Staphylococcus aureus. Frontiers in Microbiology 2016. 7: 2163 doi: 10.3389/fmicb.2016.02163

Watanabe et al. (2011) 

    Watanabe   Y , Cui   L , Katayama   Y , Kozue   K , Hiramatsu   K . Impact of rpoB mutations on reduced vancomycin susceptibility in Staphylococcus aureus. Journal of Clinical Microbiology 2011. 49: 7 2680-2684 doi: 10.1128/JCM.02144-10

Yokoyama et al. (2018) 

    Yokoyama   M , Stevens   E , Laabei   M , Bacon   L , Heesom   K , Bayliss   S , Ooi   N , O’Neill   AJ , Murray   E , Williams   P , Lubben   A , Reeksting   S , Meric   G , Pascoe   B , Sheppard   SK , Recker   M , Hurst   LD , Massey   RC . Epistasis analysis uncovers hidden antibiotic resistance-associated fitness costs hampering the evolution of MRSA. Genome Biology 2018. 19: 1 94 doi: 10.1186/s13059-018-1469-2

Yoo et al. (2013) 

    Yoo   JI , Kim   JW , Kang   GS , Kim   HS , Yoo   JS , Lee   YS . Prevalence of amino acid changes in the yvqF, vraSR, graSR, and tcaRAB genes from vancomycin intermediate resistant Staphylococcus aureus. The Journal of Microbiology 2013. 51: 2 160-165 doi: 10.1007/s12275-013-3088-7

Zhang et al. (2015) 

    Zhang   S , Sun   X , Chang   W , Dai   Y , Ma   X . Systematic review and meta-analysis of the epidemiology of vancomycin-intermediate and heterogeneous vancomycin-intermediate Staphylococcus aureus isolates. PLOS ONE 2015. 10: 8e0136082 doi: 10.1371/journal.pone.0136082
https://creativecommons.org/licenses/by/4.0/This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, reproduction and adaptation in any medium and for any purpose provided that it is properly attributed. For attribution, the original author(s), title, publication source (PeerJ) and either DOI or URL of the article must be cited.
https://www.researchpad.co/tools/openurl?pubtype=article&doi=10.7717/peerj.11764&title=Effect of genetic background on the evolution of Vancomycin-Intermediate <i>Staphylococcus aureus</i> (VISA)&author=&keyword=Epistasis,Vancomycin resistance,Experimental evolution,Staphylococcus aureus,&subject=Genomics,Microbiology,