Tracking the in vivo evolution of multidrug resistance in Staphylococcus aureus by whole-genome sequencing

Edited by John J. Mekalanos, Harvard Medical School, Boston, MA, and approved April 13, 2007
May 29, 2007
104 (22) 9451-9456


The spread of multidrug-resistant Staphylococcus aureus (MRSA) strains in the clinical environment has begun to pose serious limits to treatment options. Yet virtually nothing is known about how resistance traits are acquired in vivo. Here, we apply the power of whole-genome sequencing to identify steps in the evolution of multidrug resistance in isogenic S. aureus isolates recovered periodically from the bloodstream of a patient undergoing chemotherapy with vancomycin and other antibiotics. After extensive therapy, the bacterium developed resistance, and treatment failed. Sequencing the first vancomycin susceptible isolate and the last vancomycin nonsusceptible isolate identified genome wide only 35 point mutations in 31 loci. These mutations appeared in a sequential order in isolates that were recovered at intermittent times during chemotherapy in parallel with increasing levels of resistance. The vancomycin nonsusceptible isolates also showed a 100-fold decrease in susceptibility to daptomycin, although this antibiotic was not used in the therapy. One of the mutated loci associated with decreasing vancomycin susceptibility (the vraR operon) was found to also carry mutations in six additional vancomycin nonsusceptible S. aureus isolates belonging to different genetic backgrounds and recovered from different geographic sites. As costs drop, whole-genome sequencing will become a useful tool in elucidating complex pathways of in vivo evolution in bacterial pathogens.
Staphylococcus aureus has remained one of the most frequent causes of a wide range of both hospital- and community-acquired infections, from superficial skin and other soft tissue infections to life threatening toxic shock, pneumonia, endocarditis, and septicemia. The spectacular adaptive capacity of this pathogen resulted in the emergence and worldwide spread of lineages that acquired resistance to the majority of available antimicrobial agents. The choice of therapy against such multidrug-resistant S. aureus (MRSA) strains has been narrowed to a few antibacterial agents, among them the glycopeptide antibiotic vancomycin, which has become the mainstay of therapy worldwide. MRSA strains with reduced susceptibility to vancomycin have been reported in clinical specimen since the late 1990s (1). In most of these so-called vancomycin intermediate-resistant S. aureus (VISA) isolates, decrease in drug susceptibility, as expressed by the increase in the minimal inhibitory concentration (MIC) of vancomycin, is sufficient to cause complications in therapy and treatment failure (27). VISA-type resistance has now been identified in each of the globally spread pandemic clones of MRSA (8).
The genetic basis of VISA-type resistance to vancomycin is unknown. Unlike the most recently described and currently still rare VRSA isolates which carry the Tn1546-linked resistance mechanism (9, 10), the VISA-type isolates do not seem to carry acquired genetic elements related to drug resistance: their reduced susceptibility to vancomycin appears to be based on a gradual adaptive process.
Examination of VISA-type isolates recovered from many parts of the world showed a number of different phenotypic alterations, including changes in cell morphology and changes in the composition, thickness, and/or turnover of cell walls (11, 12). Nevertheless, associating these altered properties with the mechanism of resistance has remained problematic because of the lack of availability of an isogenic vancomycin susceptible “parental” isolate that could be used as a valid comparison. For instance, comparing the sequences of the first clinical VISA isolate MU50 to the genetically related vancomycin susceptible strain N315 identified over 174 ORFs that carried nonsynonymous changes (13, 14). However, MRSA strain N315 was isolated 15 years earlier than strain Mu50 and from a different patient. Thus, it is not clear how many of the 174 mutations are related to the mechanism of drug resistance versus the different evolutionary history of the strains.
Recently we obtained a series of MRSA isolates from the blood stream of a patient with congenital heart disease who was treated extensively with vancomycin without success (15). Available clinical data suggests that the primary site of infection was endocarditis.‡‡ In addition to vancomycin, the patient also received a single dose of rifampin and a course of therapy with the β-lactam antibiotic imipenem. After ≈12 weeks of therapy and replacement of a heart valve, the patient died because of complications of the underlying disease.
The first isolate JH1 recovered before the beginning of chemotherapy was fully susceptible to vancomycin (MIC = 1 μg/ml). Vancomycin therapy was begun between the culture isolation of JH1 and JH2. The last isolate JH9 recovered at the end of chemotherapy showed decreased susceptibility to vancomycin (MIC = 8 μg/ml). Comparison of the series of JH isolates by several genetic typing techniques indicated that they were isogenic (15, 16). The JH lineage was also related, although more remotely, to the fully sequenced MRSA strains N315 and MU50 (17).
The availability of these isogenic isolates offered a unique opportunity to identify steps in the in vivo evolution of drug resistance by sequencing the genomes of the initial (drug susceptible) and the terminal (drug resistant) isolates.

Results and Discussion

Thirty-Five Point Mutations Were Identified Between the First (JH1) and the Last (JH9) Blood Isolates from the Patient.

When the sequence data were analyzed, both JH1 and JH9 were found to contain one 2.9-Mb chromosome and a 30-kb plasmid, and the two isolates were predicted to differ by only 35 point mutations, 33 of which were on the chromosome, and two of which were on the plasmid. With the exception of two of the predicted mutations [see mutations numbered 34 and 35 in Table 1 and supporting information (SI) Appendix], all of the predicted mutations were confirmed by PCR sequencing. (In the two exceptional cases, the loci harboring the putative mutations could not be PCR sequenced because of technical reasons; see SI Appendix.) The list of 35 mutations is predicted to be nearly exhaustive, with <2 point mutations on the chromosome and zero mutations on the plasmid that are expected to have gone undetected (see SI Appendix).
Table 1.
Sequential appearance of 35 point mutations in the blood isolates
Date of isolation, month/day/yearIsolateMIC,* μg/mlNumeric identifier of mutation
VancomycinRifampinOxacillinDaptomycin12§3–678**9††10‡‡11§§1213141516171819–3334 and 35
7/20/2000JH11.00.0120.750.01               ND
9/20/2000JH24.016250.05           ND
10/1/2000JH56.0160.750.05          ND
10/6/2000JH68.0161.51.0 ND
*, E-test.
†, Numbers 1–35 are assigned to mutations in the time order of their appearance. The presence of a mutation is indicated by ●, and the absence of a mutation due to a possible reversion is indicated by ○. ND, not determined. Mutations 34 and 35 between JH1 and JH9 were predicted to be in JH9 but could not be PCR sequenced (see SI Appendix for a discussion).
‡, in blaR1 on the plasmid;
§, in SA1702 in the vraR operon;
¶, in rpoB;
‖, in rpoC;
**, in gene SA1129 with unknown function;
††, in SA1249;
‡‡, in agrC in the agr locus;
§§, in yycH in the yyc gene cluster with yycF and yycG. SAXXXX: N315 identifiers (17).
These data confirm and refine the isogenic nature of the JH isolates. The 35 mutations identified in the terminal vancomycin nonsusceptible isolate were presumably generated by the selective pressures of antibiotic treatment with the three antibiotics in the complex in vivo environment of the invaded host.

Number of Loci Effected and Categories of Mutations.

The 35 point mutations that appeared in JH9 by the end of antibiotic treatment occurred in 31 distinct loci and included 25 substitutions, 1 insertion, and 9 deletions. The first 18 of these mutations are listed in Table 2 in the time order of their appearance. All 35 mutations are described in more detail in SI Appendix. Only 6 of the 35 mutations were synonymous substitutions.
Table 2.
Descriptions of the point mutations in Table 1 found in JH1-JH6
No.*Isolate in which mutation appears firstMutated locusLocus previously implicated in resistance toDescription of mutated locusMutation
1JH1SAP011 (blaR1) (on the plasmid)β-lactams (21)Involved in regulation of β-lactamase gene blaZ and broad spectrum β-lactam resistance gene mecADeletion in A× 8 frameshifted last 70% of gene.
2JH2SA1702Vancomycin and β-lactams (23, 39, 46, 47)In operon with gene vraR, which is possibly involved in the regulation of cell wall synthesisAmino acid change H164R
3JH2SA0500 (rpoB)Rifampin (18)Codes for β-subunit of RNA polymeraseAmino acid change D471Y
4    Amino acid change A473S
5    Amino acid change A477S
6    Amino acid change E478D
7JH2SA0501 (rpoC)Daptomycin (38)Codes for β'-subunit of RNA polymeraseAmino acid change E854K
8JH2SA1129 Contains match to RNA binding motifAmino acid change D296Y
9JH5SA1249 Putative gene possibly in an operon with the gene murG, which is involved in cell wall synthesisDeletion in G×9 frameshifted last 80% of gene.
10JH6SA1843 (agrC)Vancomycin (24, 25)Part of agr locus involved in qorum sensing and regulation of the expression of virulence and cell surface proteinsDeletion in T×9 frameshifted last 70% of gene.
11JH6SA0019 (yycH)Daptomycin (38)In a gene cluster with the gene yycF and possibly involved in the regulation of the autolysin gene lytMSubstitution introduced stop codon, truncating protein to 10% of its length.
12JH6Between divergently transcribed genes SAS014 and SA0411 (ndhF) NdhF is the F-subunit of NADH dehydrogenase.Deletion in T×7 upstream of both SAS014 and SA0411
13JH6SA0582 Similar to Na+/H+ antiporter subunit MrpE in B. subtilisSynonymous substitution
14JH6SA0980 (isdE) Involved in passage of heme-iron to cytoplasm during pathogenesisAmino acid change A84V
15JH6SA1659 (prsA) Chaperone that assists post-translocational folding of proteins at the cytoplasmic/cell wall interfaceDeletion in A×7 frameshifted last 15% of gene.
16JH6SA2094 Similar to malic/Na+-lactate antiporter MleN in B. subtilisAmino acid change A94T
17JH6Between divergently transcribed genes SA2125 and SA2126 SA2125 matches protein family consisting of arginases, agmatinases, and formiminoglutamases.Substitution upstream of both SA2125 and SA2126
18JH6SA2320 (pfoR) Contains a match to a domain of a sugar specific permeaseSynonymous substitution
Because JH6 and JH9 have nearly identical MICs to the four antibiotics that were tested, the mutations 19–33 found only in JH9 are not described here. For descriptions of all the mutations, see SI Appendix, which gives the exact positions of the mutations in nucleotide coordinates and many references to relevant research articles.
*, Lists the numeric identifiers of the mutations used in Table 1.
†, The mutated locus is on the chromosome unless otherwise indicated. The SAXXXX and SAPXXX identifiers are the N315 identifiers for the genes (17). If a more descriptive identifier has been assigned to a gene, it is listed in parentheses.
‡, The nomenclature N × α is used to describe a run of identical bases: α is the number of bases in the run, and N is the repeated base.

Sequential Appearance of Antibiotic Resistance and Mutations in the JH Isolates During Antibiotic Treatment.

The time course of antibiotic treatment is shown in Fig. 1, and Table 1 summarizes the parallel changes observed in the antibiotic susceptibility profiles of the isolates. To explore parallels between the changes observed in antibiotic susceptibility and genetic changes, loci found mutated between JH1 and JH9 were PCR sequenced in the intermittent blood isolates JH2, JH5, and JH6. The same loci were also PCR sequenced in an additional fully vancomycin susceptible isolate JH15 that was recovered from the nares of a healthy contact of the patient (15).
Fig. 1.
Dates of antibiotic exposure and recoveries of MRSA isolates (month/day in year 2000).
The results summarized in Table 1 allowed us to make several conclusions.
The first fully antibiotic susceptible blood isolate JH1 and the also fully susceptible contact isolate JH15 (data not shown) had identically low MICs to the antibiotics tested and carried only mutation 1.
Once a mutation appeared in an early blood isolate, it was retained in all subsequent blood isolates (with the single exception of mutation 9). Thus, the genetic changes appeared in a sequential order in parallel with the increasing vancomycin MIC values.
Following the brief rifampin therapy, the rifampin MIC of isolate JH1 (0.012 μg/ml) increased to 16 μg/ml in JH2 and remained at that level in all subsequent isolates.
In parallel with the treatment with the beta-lactam antibiotic imipenem, the beta-lactam (oxacillin) MIC of strain JH1 (0.75 μg/ml) increased to 25 μg/ml in JH2 but then declined to between 0.75 and 1.5 μg/ml in all subsequent isolates.
The vancomycin susceptibility of the isolates declined gradually in several discrete steps: the vancomycin MIC of strain JH1 (1 μg/ml) increased to 4 μg/ml in JH2, 6 μg/ml in JH5, and its maximum value of 8 μg/ml in JH6 and JH9.
The MIC value for daptomycin also increased in the JH isolates, although this newly introduced antimicrobial agent had not been used in the therapy. The daptomycin MIC of JH1 (0.01 μg/ml) increased to 0.05 μg/ml in JH2 and JH5 and then underwent a further jump to 1.0 μg/ml in JH6 and JH9.

Appearance of Mutations Associated with Rifampin Resistance.

The >1,000-fold increase in the rifampin MIC from 0.012 to 16 μg/ml between JH1 and JH2 is likely related to mutations 3–6 (Table 1), involving four amino acid changes (D471Y, A473S, A477S, and E478D) in the β-subunit RpoB of RNA polymerase (Table 2). The change D471Y has alone been shown to confer rifampin resistance, and all four changes occur in the region of amino acids 463–550 found to harbor the majority of mutations responsible for rifampin resistance in S. aureus (see ref. 18 for review). Most rifampin resistant S. aureus mutants are reported to have only one or two amino acid changes in RpoB. However, these changes frequently reduce fitness. It is possible that some of the four changes we report here in RpoB are compensatory mutations that helped to offset a loss in fitness. In prior work, rifampin resistant Escherichia coli mutants were passaged through hundreds of generations in vitro and were observed to evolve increased fitness by compensatory mutations in RpoB rather than by reversion to drug sensitivity (19).

Appearance of a Mutation Associated with β-Lactam Resistance.

The sharp rise in resistance to oxacillin (a β-lactam) between JH1 and JH2 is accompanied by the reversal of mutation 1. Found in both JH1 and JH15, mutation 1 involves a deletion of an adenine in a run of eight adenines that frameshifted 70% of blaR1 on the plasmid. The transmembrane sensor BlaR1 and its cognate cytosolic repressor BlaI are part of a signaling pathway that senses β-lactams and induces the expression of the β-lactamase gene blaZ and the broad spectrum β-lactam resistance gene mecA. It is known that β-lactams bind to and thereby activate BlaR1, which then in turn deactivates BlaI to alleviate the BlaI repression of blaZ and mecA (20). Therefore, the inactivation of BlaR1 would be expected to lead to the BlaI mediated constitutive repression of both blaZ and mecA. Hence, it is no surprise that a frameshift in blaR1 has been previously shown to abolish resistance to β-lactams (including to β-lactams, like imipenem and oxacillin, which are not susceptible to cleavage by β-lactamase) (21). The reversal of the frameshift in blaR1 may therefore be responsible for the increase in oxacillin resistance in JH2. Even though blaR1 remains in frame in JH2 onwards, the oxacillin resistance drops in subsequent isolates (beginning with JH5) as the vancomycin resistance increases. Such an inverse relationship between oxacillin and vancomycin resistance has been demonstrated in vitro (22). The mechanism of this effect is not known.

Stages in Vancomycin Resistance.

Between JH1 and JH2, the vancomycin MIC increases from 1.0 to 4.0 μg/ml. Along with this increase in resistance appears mutation 2, involving the amino acid change H164R in the protein SA1702. Although the gene SA1702 has an unknown function, it is in an operon with the gene vraR, which codes for a response regulator of a two-component system. Up-regulation of vraR has been shown to increase the vancomycin MIC by 4-fold in one S. aureus strain (23), and increased transcription of vraR in JH9 compared with JH1 was observed in a study with DNA microarrays (16). A further increase (from 4 to 6 μg/ml) in vancomycin MIC occurs between JH2 and JH5 and is accompanied by a single additional mutation in SA1249, a gene of unknown function.

Between JH5 and JH6.

The vancomycin MIC increases further from 6 to 8 μg/ml, associated with the appearance of mutations 10–18. Mutation 10 appears to frameshift 70% of agrC in the agr locus. Involved in quorum sensing, the agr locus regulates many virulence and cell surface genes. Loss of the locus has already been proposed to be associated with the VISA phenotype (24, 25). The increase in the vancomycin MIC between JH5 and JH6 is also linked to mutation 11, involving a substitution that introduces a premature stop codon in yycH, truncating the gene to 10% of its length. The gene yycH is in the yyc gene cluster with the gene yycF, which codes for the response regulator of a two-component system. In Bacillus subtilis, deletion of yycH results in the up-regulation of genes controlled by yycF (26). In S. aureus, the regulator YycF was shown to bind upstream of the gene lytM, coding for an autolysin (27). Hence, loss of yycH might be expected to lead to up-regulation of lytM, and indeed, transcriptome profiling has shown that the expression of lytM is 5-fold greater in JH9 compared with JH1 (16). Vancomycin blocks access of LytM to its cell-wall substrate (28), and the up-regulation of lytM may be a response to this inhibition. Although lytM is overexpressed in JH9, JH9 was shown to have decreased susceptibility to autolysis, which may be related to an abnormality and/or over-production of teichoic acids (11).

Reversions of Mutation 9 and the Potential Role of Homopolymeric Tracts in the Evolution of Resistance in the JH Isolates.

There are two possible explanations why mutation 9 (involving a frameshift in SA1249) reverts between the blood isolates JH5, JH6, and JH9 (Tables 1 and 2). It is possible that the population of S. aureus in the blood of the patient was heterogeneous, with some individuals carrying the wild-type allele of SA1249 and other individuals carrying the mutant allele. Alternatively, mutation 9 may have appeared and fixed in all of the bacteria, then disappeared entirely, and finally reappeared and fixed again. Both scenarios are plausible, because mutation 9 is a deletion of a cystine in a string of nine cystines. Because of slippage of the DNA polymerase, a homopolymeric tract can rapidly increase or decrease in length (29). In our list of 33 confirmed point mutations, changes in runs of identical nucleotides are statistically overrepresented (P = 10−7) (SI Appendix). Our 33 confirmed point mutations include 1-bp expansions or contractions of eight distinct homopolymeric segments, all initially ≥6 bp (mutations 1, 9, 10, 12, 15, 21, 22, and 33). Because rapid changes in phenotype have been associated with variable length simple repeats (e.g., phase variation in Neisseria meningitidis) (30), it is tempting to speculate that runs of identical nucleotides in S. aureus can promote quick switching between resistant phenotypes. In particular, the VISA phenotype has been noted to be unstable (31).

Sequential Appearance of Mutations.

The data summarized in Table 1 strongly suggest that we are observing genetic alterations in a single S. aureus lineage under the selective pressure of the antimicrobial agents used for the therapy. The selection for these mutations and the emergence of genetically altered bacterial populations most likely occurred at the primary site of infection: in the endocarditic vegetation of the infected heart valve, which is known to be capable of carrying large (109–1011 bacteria per g of tissue) and physiologically heterogeneous populations of bacteria (32), which then “seed” the blood stream followed by rapid clearance (33).
We also had access to a single additional isolate (JH14), which was recovered from the heart-valve after replacement surgery performed on 17 November 2000 (Fig. 1). The heart valve isolate JH14 had exactly the same antibiotic resistance and mutational profiles as JH6: resistance to rifampin, decreased susceptibility to daptomycin (MIC = 1 μg/ml) and to vancomycin (MIC = 8 μg/ml), and the same set of confirmed mutations (28 and 1018) identified in isolate JH6 and also carried in JH9. These observations suggest that the “extra” 16 confirmed mutations (9 and 1933) identified only in JH9 reflect some unknown selective pressure operative in the complex and heterogeneous environment of the endocarditic lesion but are not essential for the decreased susceptibility to vancomycin. Isolate JH9 may have originated from an area of the endocarditic vegetation where growth and replication of bacteria can be slow (32), leading to the appearance of traits, such as antibiotic tolerance and/or persistence, which are often associated with slow growth. Decreased growth rate and slow loss of viability during antibiotic treatment has been documented in JH9 (11).

Decreased Susceptibility to Daptomycin in the JH Isolates.

We were surprised to find that the MIC value for daptomycin also increased in the JH isolates, although this antibiotic was not used in the therapy.
The antibiotic daptomycin was introduced in 2003 as a therapeutic agent against MRSA infections, including strains with decreased susceptibility to vancomycin (34). The current clinical breakpoint for resistance to daptomycin is defined as an MIC ≥ 2 μg/ml (7). Although the effectiveness of daptomycin against a number of clinical VISA strains was demonstrated (35), more recent studies found a strong positive correlation between reduced susceptibility to vancomycin and daptomycin among VISA isolates (36, 37). In one of these studies (36), 70 independent clinical and laboratory VISA isolates with vancomycin MICs between 4–16 μg/ml were tested, and >80% were observed to have daptomycin MICs ≥ 2 μg/ml.
The genetic mechanisms underlying decreased susceptibility to both vancomycin and daptomycin are unknown. It was therefore of interest that the late JH isolates harbor mutation 7 in rpoC coding for the β′-subunit of RNA polymerase and mutation 11 in the yyc gene cluster (Table 1). These loci were two of only several loci found genome-wide to be mutated in laboratory generated S. aureus mutants with reduced susceptibility to daptomycin (38). In our study, the daptomycin MIC of JHI (0.01 μg/ml) increased to 0.05 μg/ml in JH2-JH5 followed by a sharp increase to 1 μg/ml in JH6 onwards. These increases in resistance coincided with the mutations in rpoC and the yyc gene cluster.

A Limited Number of Mutations Can Cause Extensive Changes in the Gene Expression Profile.

Although only 35 point mutations in 31 loci are observed between JH1 and JH9, prior transcriptome profiling work found 224 genes that were differentially expressed by 2-fold or more between the isolates (16). We assessed the overlap between this list of 224 genes and three other lists of transcriptional changes described in the literature: the changes in the expression of 48 genes controlled by vraR (39); 244 genes controlled by positive regulators of the agr locus (4042); and 32 genes controlled by yycF (27). Each of the loci vraR, agr, and yycF codes for a transcriptional regulator that appears to be involved with decreased susceptibility to vancomycin. In each case, the overlap was statistically significant (P = 10−11, 10−2, and 10−2, respectively) (see SI Appendix). Thus, many of the mRNA expression changes between JH1 and JH9 are likely due to the mutations 2, 10, and 11 (Table 1). Additionally, mutations 3–6 (Table 1) may effect the expression of many genes, because nonsynonymous substitutions in rpoB conferring rifampin resistance in B. subtilis were shown to have pleitropic effects (43).

The Relationship of the JH Isolates to Previously Sequenced S. aureus.

In prior work (15), multilocus sequence typing (MLST) (44) showed that the JH isolates (each with an MLST 1-4-1-4-12-1-28) were closely related to two sequenced clinical S. aureus isolates: the MRSA strain N315 and the rifampin resistant VISA MU50 (both with a MLST 1-4-1-4-12-1-10). Similar to JH1 and JH9, N315 and MU50 were each found when sequenced to contain one 2.8- to 2.9-Mb chromosome and a 20- to 30-kb plasmid (17). In this work, multialignments of orthologous JH1, N315, and MU50 chromosomal sequence were constructed, and a randomly chosen subset of 50 cases were examined where the N315 and MU50 sequences differed from each other by an isolated nucleotide difference. Among these cases, JH1 was found to agree with N315 three-fourths of the time, indicating that JH1 is more closely related to N315 than MU50. Despite the similar MLSTs of JH1 and N315, a full genome-wide comparison indicated that the two isolates differ on the chromosome by >500 mutations ranging in size from 1 to 50,000 bp (see SI Appendix). Almost all of the >500 reported mutations are expected to be real, with at most only a very small fraction representing sequencing errors (see SI Appendix). As for the JH1 plasmid, it was found to be a chimera of the N315 and MU50 plasmids (SI Appendix). The sheer number of mutations between JH1 and N315 again highlights the need for the careful selection of isolates for an intergenomic comparison to be informative.

Common Mutations in Geographically Diverse VISA Isolates.

It was of interest to check whether any of the loci mutated in JH9 were also mutated in the fully sequenced rifampin resistant VISA strain Mu50. MU50 was found to harbor mutations in the genes vraS in the vraR operon, agrA in the agr locus, and rpoB. In MU50, the mutation in rpoB (a single amino acid change H481Y in RpoB) was shown to confer rifampin resistance (18). It has already been shown that many geographically diverse VISA isolates carry loss of function mutations in the agr locus (24, 25).

The Status of the vraR Operon in Geographically Diverse VISA Isolates.

Because the vraR operon was mutated in both JH9 and MU50, we proceeded to PCR sequence the vraR operon in other VISA isolates. In total, we knew or determined the sequence of the vraR operon in six geographically diverse clinical VISA isolates (vancomycin MICs = 4–8 μg/ml) from different genetic backgrounds: JH9 isolated in Maryland,‡‡ MU50 from Japan (17), PC3 from New York (3), VNJ from New Jersey (45), VMI from Michigan (45), and HSMB1 from Portugal (H.d.L., unpublished data). We also determined the sequence of the vraR operon in a laboratory mutant VM3 (vancomycin MIC = 6 μg/ml) (22). In all seven VISA isolates, a mutation was found in the vraR operon, and this was shown to be highly statistically significant [P < 0.001 (possibly much less)] (see SI Appendix).
In Table 1, mutation 2 in SA1702 in the vraR operon is one of only eight genetic changes that occur between JH1 and JH2, during which time the vancomycin MIC value increases from 1 to 4 μg/ml. Of the other mutations, mutation 1 in blaR1 and mutations 3–6 in rpoB are likely to be related to β-lactam and rifampin resistance respectively. Thus, the preponderance of evidence strongly suggests that genetic alterations in the vraR operon, a key regulatory system involved with monitoring cell-wall biosynthesis in S. aureus (39, 46, 47), plays an important role in the VISA type vancomycin resistance.

From Genetic Change to Altered Phenotype.

Our discovery of an ordered series of mutations leading to a level of VISA type vancomycin resistance that can compromise therapy provides a short list of mutations that should enable us to design experiments to reconstruct the vancomycin resistant phenotype in susceptible strains of S. aureus. Polygenic traits are very hard to trace by association studies, and even expression array data may be hard to interpret when regulators with pleiotropic functions are involved. Our short list of mutant loci can easily be screened in other VISA strains to determine their prevalence.
With many new inexpensive differential sequencing technologies (48) on the horizon, subsequent clinical studies could collect more samples to temporally resolve the appearance of single mutations and the population structure of the bacteria as has been done for HIV (49). The spread and emergence of virulent strains among patients could be followed.


For full details, see SI Appendix.
For JH1 (JH9), the whole-genome shotgun sequencing was carried out to a mean coverage of 8.5× (9.5×), with 98.5% (97.7%) of the chromosome having a coverage ≥1×. The JH1 and JH9 assembled sequences and reads were multialigned, and a Bayesian probabilistic model was used to identify mutations in the multialignment. In 94% of the genome, the read coverage and quality was sufficient to call mutations between JH1 and JH9 with an expected error rate (both false positive and negative) of zero. The number of mutations estimated to have gone undetected in the remaining 6% of the genome is <2 point mutations. PCR sequencing was done to confirm the model's predictions and to rule out marginal cases.


minimal inhibitory concentration
multidrug-resistant S. aureus
vancomycin intermediate-resistant S. aureus.


We thank Seth Darst, Lu Bai, and Alex Morozov for valuable discussions on the genes rpoB and rpoC and Neal Steigbigel for discussion on medical/microbiological aspects of the paper. This work was supported by United States Public Health Service Grant RO1 AI37275NIH (to A.T.) and National Science Foundation Grant DMR 129848 (to E.D.S.).

Supporting Information

Adobe PDF - 09839Appendix.pdf
Adobe PDF - 09839Appendix.pdf


PC Appelbaum Clin Microbiol Infect 12, 16–23 (2006).
K Hiramatsu, H Hanaki, T Ino, K Yabuta, T Oguri, FC Tenover J Antimicrob Chemother 40, 135–136 (1997).
K Sieradzki, RB Roberts, SW Haber, A Tomasz N Engl J Med 340, 517–523 (1999).
PG Charles, PB Ward, PD Johnson, BP Howden, ML Grayson Clin Infect Dis 38, 448–451 (2004).
BP Howden, PB Ward, PG Charles, TM Korman, A Fuller, P du Cros, EA Grabsch, SA Roberts, J Robson, K Read, et al. Clin Infect Dis 38, 521–528 (2004).
BP Howden Intern Med J 35, S136–S140 (2005).
Performance Standards for Antimicrobial Susceptibility Testing: 16th informational supplement M100–S16 (Clinical and Laboratory Standards Institute, Wayne, PA, 2006).
RA Howe, A Monk, M Wootton, TR Walsh, MC Enright Emerg Infect Dis 10, 855–857 (2004).
LM Weigel, DB Clewell, SR Gill, NC Clark, LK McDougal, SE Flannagan, JF Kolonay, J Shetty, GE Killgore, FC Tenover Science 302, 1569–1571 (2003).
A Severin, K Tabei, F Tenover, M Chung, N Clarke, A Tomasz J Biol Chem 279, 3398–3407 (2004).
K Sieradzki, A Tomasz J Bacteriol 185, 7103–7110 (2003).
RF Pfeltz, BJ Wilkinson Curr Drug Targets Infect Disord 4, 273–294 (2004).
MB Avison, PM Bennett, RA Howe, TR Walsh J Antimicrob Chemother 49, 255–260 (2002).
T Ohta, H Hirakawa, K Morikawa, A Maruyama, Y Inose, A Yamashita, K Oshima, M Kuroda, M Hattori, K Hiramatsu, et al. DNA Res 11, 51–56 (2004).
K Sieradzki, T Leski, J Dick, L Borio, A Tomasz J Clin Microbiol 41, 1687–1693 (2003).
F McAleese, SW Wu, K Sieradzki, P Dunman, E Murphy, S Projan, A Tomasz J Bacteriol 188, 1120–1133 (2006).
M Kuroda, T Ohta, I Uchiyama, T Baba, H Yuzawa, I Kobayashi, L Cui, A Oguchi, K Aoki, Y Nagai, et al. Lancet 357, 1225–1240 (2001).
AJ O'Neill, T Huovinen, CW Fishwick, I Chopra Antimicrob Agents Chemother 50, 298–309 (2006).
MG Reynolds Genetics 156, 1471–1481 (2000).
HZ Zhang, CJ Hackbarth, KM Chansky, HF Chambers Science 291, 1962–1965 (2001).
CJ Hackbarth, C Miick, HF Chambers Antimicrob Agents Chemother 38, 2568–2571 (1994).
K Sieradzki, A Tomasz J Bacteriol 181, 7566–7570 (1999).
M Kuroda, K Kuwahara-Arai, K Hiramatsu Biochem Biophys Res Commun 269, 485–490 (2000).
G Sakoulas, GM Eliopoulos, RC Moellering, C Wennersten, L Venkataraman, RP Novick, HS Gold Antimicrob Agents Chemother 46, 1492–1502 (2002).
G Sakoulas, GM Eliopoulos, VG Fowler, RC Moellering, RP Novick, N Lucindo, MR Yeaman, AS Bayer Antimicrob Agents Chemother 49, 2687–2692 (2005).
H Szurmant, K Nelson, EJ Kim, M Perego, JA Hoch J Bacteriol 187, 5419–5426 (2005).
S Dubrac, T Msadek J Bacteriol 186, 1175–1181 (2004).
K Sieradzki, A Tomasz Antimicrob Agents Chemother 50, 527–533 (2006).
RM Schaaper, RL Dunn Proc Natl Acad Sci USA 84, 6220–6224 (1987).
ER Moxon, PB Rainey, MA Nowak, RE Lenski Curr Biol 4, 24–33 (1994).
S Boyle-Vavra, SK Berke, JC Lee, RS Daum Antimicrob Agents Chemother 44, 272–277 (2000).
DT Durack, PB Beeson Br J Exp Pathol 53, 50–53 (1972).
PB Beeson, ES Brannon, JV Warren J Exp Med 81, 9–23 (1945).
JD Alder Drugs Today (Barc) 41, 81–90 (2005).
PJ Petersen, PA Bradford, WJ Weiss, TM Murphy, PE Sum, SJ Projan Antimicrob Agents Chemother 46, 2595–2601 (2002).
JB Patel, LA Jevitt, J Hageman, LC McDonald, FC Tenover Clin Infect Dis 42, 1652–1653 (2006).
L Cui, E Tominaga, HM Neoh, K Hiramatsu Antimicrob Agents Chemother 50, 1079–1082 (2006).
L Friedman, JD Alder, JA Silverman Antimicrob Agents Chemother 50, 2137–2145 (2006).
M Kuroda, H Kuroda, T Oshima, F Takeuchi, H Mori, K Hiramatsu Mol Microbiol 49, 807–821 (2003).
PM Dunman, E Murphy, S Haney, D Palacios, G Tucker-Kellogg, S Wu, EL Brown, RJ Zagursky, D Shlaes, SJ Projan J Bacteriol 183, 7341–7353 (2001).
M Korem, Y Gov, MD Kiran, N Balaban Infect Immun 73, 6220–6228 (2005).
X Liang, L Zheng, C Landwehr, D Lunsford, D Holmes, Y Ji J Bacteriol 187, 5486–5492 (2005).
H Maughan, B Galeano, WL Nicholson J Bacteriol 186, 2481–2486 (2004).
MC Enright, NP Day, CE Davies, SJ Peacock, BG Spratt J Clin Microbiol 38, 1008–1015 (2000).
Morbid Mortal Wkly Rep 46, 813–815 (1997).
S Yin, RS Daum, S Boyle-Vavra Antimicrob Agents Chemother 50, 336–343 (2006).
S Gardete, SW Wu, S Gill, A Tomasz Antimicrob Agents Chemother 50, 3424–3434 (2006).
RF Service Science 311, 1544–1546 (2006).
A Rambaut, D Posada, KA Crandall, EC Holmes Nat Rev Genet 5, 52–61 (2004).

Information & Authors


Published in

Go to Proceedings of the National Academy of Sciences
Go to Proceedings of the National Academy of Sciences
Proceedings of the National Academy of Sciences
Vol. 104 | No. 22
May 29, 2007
PubMed: 17517606


Submission history

Received: November 6, 2006
Published online: May 29, 2007
Published in issue: May 29, 2007


We thank Seth Darst, Lu Bai, and Alex Morozov for valuable discussions on the genes rpoB and rpoC and Neal Steigbigel for discussion on medical/microbiological aspects of the paper. This work was supported by United States Public Health Service Grant RO1 AI37275NIH (to A.T.) and National Science Foundation Grant DMR 129848 (to E.D.S.).


This article is a PNAS Direct Submission.
Flayhart, D., Hanlon, A., Wakefield, T., Ross, T., Borio, L., Dick, J. (2001) in Abstracts of the 101st General Meeting of the American Society of Microbiology, May 20–24, 2001, Orlando, FL, Abstr. A-39.
This article contains supporting information online at



Michael M. Mwangi
Physics Department, Cornell University, Ithaca, NY 14850;
Center for Studies in Physics and Biology and
Shang Wei Wu
Laboratory of Microbiology, The Rockefeller University, New York, NY 10021;
Department of Microbiology, Tianjin Medical University, Tianjin 300070, People's Republic of China;
Yanjiao Zhou
Laboratory of Microbiology, The Rockefeller University, New York, NY 10021;
Department of Microbiology, Tianjin Medical University, Tianjin 300070, People's Republic of China;
Krzysztof Sieradzki
Laboratory of Microbiology, The Rockefeller University, New York, NY 10021;
Herminia de Lencastre
Laboratory of Microbiology, The Rockefeller University, New York, NY 10021;
Laboratory of Molecular Genetics, Instituto de Tecnologia Química e Biológica, Universidade Nova de Lisboa, Oeiras, Portugal;
Paul Richardson
United States Department of Energy Joint Genomic Institute, Walnut Creek, CA 94598; and
David Bruce
United States Department of Energy Joint Genomic Institute, Walnut Creek, CA 94598; and
Edward Rubin
United States Department of Energy Joint Genomic Institute, Walnut Creek, CA 94598; and
Eugene Myers
Howard Hughes Medical Institute, Janelia Farm Research Campus, Ashburn, VA 20146
Eric D. Siggia
Physics Department, Cornell University, Ithaca, NY 14850;
Center for Studies in Physics and Biology and
Alexander Tomasz†† [email protected]
Laboratory of Microbiology, The Rockefeller University, New York, NY 10021;


To whom correspondence should be addressed. E-mail: [email protected]
Author contributions: A.T. designed research; M.M.M., S.W.W., and Y.Z. performed research; K.S., P.R., D.B., and E.R. contributed new reagents/analytic tools; M.M.M. H.d.L., E.M., and E.D.S. analyzed data; and M.M.M., E.D.S., and A.T. wrote the paper.

Competing Interests

The authors declare no conflict of interest.

Metrics & Citations


Note: The article usage is presented with a three- to four-day delay and will update daily once available. Due to ths delay, usage data will not appear immediately following publication. Citation information is sourced from Crossref Cited-by service.

Citation statements



If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

Cited by


    View Options

    View options

    PDF format

    Download this article as a PDF file


    Get Access

    Login options

    Check if you have access through your login credentials or your institution to get full access on this article.

    Personal login Institutional Login

    Recommend to a librarian

    Recommend PNAS to a Librarian

    Purchase options

    Purchase this article to get full access to it.

    Single Article Purchase

    Tracking the in vivo evolution of multidrug resistance in Staphylococcus aureus by whole-genome sequencing
    Proceedings of the National Academy of Sciences
    • Vol. 104
    • No. 22
    • pp. 9101-9546







    Share article link

    Share on social media