Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Research and News from ASM Journals
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSphere
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Research and News from ASM Journals
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Clinical Science and Epidemiology

Population Genomics of Reduced Vancomycin Susceptibility in Staphylococcus aureus

Lavanya Rishishwar, Colleen S. Kraft, I. King Jordan
Brandi M. Limbago, Editor
Lavanya Rishishwar
aSchool of Biology, Georgia Institute of Technology, Atlanta, Georgia, USA
bApplied Bioinformatics Laboratory, Atlanta, Georgia, USA
cPanAmerican Bioinformatics Institute, Cali, Valle del Cauca, Colombia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Colleen S. Kraft
dPathology and Laboratory Medicine, Emory University School of Medicine, Atlanta, Georgia, USA
eAntibiotic Resistance Center, Emory University School of Medicine, Atlanta, Georgia, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
I. King Jordan
aSchool of Biology, Georgia Institute of Technology, Atlanta, Georgia, USA
bApplied Bioinformatics Laboratory, Atlanta, Georgia, USA
cPanAmerican Bioinformatics Institute, Cali, Valle del Cauca, Colombia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Brandi M. Limbago
Centers for Disease Control and Prevention
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSphere.00094-16
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

The increased prevalence of vancomycin-intermediate Staphylococcus aureus (VISA) is an emerging health care threat. Genome-based comparative methods hold great promise to uncover the genetic basis of the VISA phenotype, which remains obscure. S. aureus isolates were collected from a single individual that presented with recurrent staphylococcal bacteremia at three time points, and the isolates showed successively reduced levels of vancomycin susceptibility. A population genomic approach was taken to compare patient S. aureus isolates with decreasing vancomycin susceptibility across the three time points. To do this, patient isolates were sequenced to high coverage (~500×), and sequence reads were used to model site-specific allelic variation within and between isolate populations. Population genetic methods were then applied to evaluate the overall levels of variation across the three time points and to identify individual variants that show anomalous levels of allelic change between populations. A successive reduction in the overall levels of population genomic variation was observed across the three time points, consistent with a population bottleneck resulting from antibiotic treatment. Despite this overall reduction in variation, a number of individual mutations were swept to high frequency in the VISA population. These mutations were implicated as potentially involved in the VISA phenotype and interrogated with respect to their functional roles. This approach allowed us to identify a number of mutations previously implicated in VISA along with allelic changes within a novel class of genes, encoding LPXTG motif-containing cell-wall-anchoring proteins, which shed light on a novel mechanistic aspect of vancomycin resistance.

IMPORTANCE The emergence and spread of antibiotic resistance among bacterial pathogens are two of the gravest threats to public health facing the world today. We report the development and application of a novel population genomic technique aimed at uncovering the evolutionary dynamics and genetic determinants of antibiotic resistance in Staphylococcus aureus. This method was applied to S. aureus cultures isolated from a single patient who showed decreased susceptibility to the vancomycin antibiotic over time. Our approach relies on the increased resolution afforded by next-generation genome-sequencing technology, and it allowed us to discover a number of S. aureus mutations, in both known and novel gene targets, which appear to have evolved under adaptive pressure to evade vancomycin mechanisms of action. The approach we lay out in this work can be applied to resistance to any number of antibiotics across numerous species of bacterial pathogens.

INTRODUCTION

Staphylococcus aureus infections are a major cause of mortality and morbidity worldwide (1). Incidence rates range from 20 to 50 cases per 100,000 individuals with 10 to 30% mortality (2, 3). Elimination of S. aureus infections typically requires a prolonged course of antibiotics (4), and evolved antibiotic resistance is a major challenge to the effective treatment of S. aureus infections (5). S. aureus isolates resistant to the methicillin antibiotic were first identified in the early 1960s (6, 7). Since that time, methicillin-resistant S. aureus (MRSA) has become increasingly widespread and is now a leading cause of hospital-acquired infections (8). The mechanism of methicillin resistance is well understood and involves the acquisition of a single-mobile-element-borne gene mecA (9, 10).

MRSA is so common that hospital-acquired infections are typically assumed to be methicillin resistant, and patient treatment is accordingly initiated with vancomycin, which has emerged as the mainstay of S. aureus infection therapy (11). Despite the overall efficacy of vancomycin as an antibiotic, S. aureus resistance to vancomycin is becoming increasingly prevalent (12, 13). The first cases of intermediate resistance to vancomycin were identified in 1996 (14), and fully resistant strains were later found in 2002 (15). As is the case for methicillin resistance, full vancomycin resistance is based on the acquisition of a single gene, vanA (16, 17). Vancomycin-resistant S. aureus (VRSA) is extremely rare, with only 12 cases reported in the United States since 2002 (15), and therefore does not represent an urgent public health threat.

Vancomycin-intermediate resistance is defined on the basis of the range of MICs of the antibiotic needed to inhibit growth. S. aureus strains that show MICs from 3 to 8 µg/ml are characterized as vancomycin-intermediate S. aureus (VISA) (13). Unlike VRSA, the incidence of VISA infections is steadily rising (12, 13), a phenomenon referred to as vancomycin MIC creep (18), and VISA therefore does pose a serious potential threat to the effective antibiotic treatment of S. aureus.

At this time, the precise genetic basis of the VISA phenotype is unknown. VISA does not appear to be a single gene phenomenon, as seen for MRSA and VRSA, and different VISA cases may be caused by different gene sets (19). The lack of a clear-cut gene-to-VISA relationship, combined with the public health relevance of increased vancomycin resistance, has stimulated numerous attempts to uncover the genetic basis of VISA (20–26). One specific approach that can be taken to address this question entails the comparison of genome sequences between S. aureus isolates with different vancomycin susceptibility profiles (22, 27–33). The goal of such studies is to identify mutated genes that are (i) exclusively found in isolates with reduced vancomycin susceptibility and (ii) encode proteins with plausible roles in the VISA phenotype (e.g., cell-wall-related functions).

The genomic approach to studying VISA is typically performed by comparing genome consensus sequences that are assembled from overlapping sequence reads that cover the entire S. aureus genome multiple times. The consensus sequence approach implicitly assumes that patient isolates are monoclonal, since it collapses all observed site-specific variation among sequence reads to a single consensus sequence, i.e., it considers sequence variation among colocated reads as noise that should be removed. Nevertheless, as the phenomenon of vancomycin creep is based on evolution, it is highly likely that many S. aureus patient isolates are in fact polyclonal and characterized by subpopulations that bear VISA-related mutations at low frequencies. In fact, it is known that there are many cases of heterogeneous VISA (hVISA) that correspond precisely to this description (34). hVISA isolates have overall MICs that fall within the vancomycin-susceptible range, but they also contain low-frequency subpopulations that are demonstrably less susceptible (35). For such hVISA populations, treatment with vancomycin could cause less-susceptible subpopulations to increase in frequency along with their relevant VISA-related mutations (36).

In this study, a novel population genomic approach was developed and applied in an effort to characterize the genomic determinants of the emergence of reduced vancomycin susceptibility in a single patient (Fig. 1). This patient was admitted to Emory Healthcare with recurrent S. aureus infections several times over the course of 1 year, and at each successive time point, the patient’s isolates showed increased vancomycin MIC levels. The population genomic approach employed here does not assume that the patient’s S. aureus isolates are monoclonal. Instead, it uses site-specific S. aureus sequence variation uncovered by deep sequencing to model the allele frequency dynamics between patient isolates that were taken from different time points and showed different vancomycin MICs. This approach provides additional resolution for the study of the evolution of vancomycin resistance and for the detection of VISA-implicated mutations compared to the consensus sequence approach. A number of mutations in candidate genes previously implicated in the VISA phenotype were found using this approach along with mutations distributed among members of a novel class of VISA-related genes that encode cell-wall-related functions.

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Overview of the S. aureus population genomic approach used in this study.

RESULTS AND DISCUSSION

Presenting patient.In January 2012, a 50-year-old African-American male with a history of uncontrolled diabetes and end-stage renal disease requiring hemodialysis presented at Emory Healthcare and was found to have an S. aureus bloodstream infection. The patient was treated with intravenous vancomycin for 22 days for this bacteremia. The patient’s S. aureus isolate from this time point was determined to be methicillin resistant and vancomycin susceptible with an MIC of 1.5 µg/ml by Etest on Mueller-Hinton agar. The patient was discharged and readmitted in April with recurrent S. aureus bloodstream infection, complicated by septic arthritis of the shoulder. Vancomycin treatment was reinitiated after the infection was cultured and continued for 44 days. The patient’s second time point isolate was found to have a vancomycin MIC of 2.0 µg/ml. The patient returned in August with recurrent S. aureus infection and received treatment with 40 additional days of intravenous vancomycin treatment. In November, the patient was again readmitted with recurrent S. aureus bloodstream infection, complicated by discitis and osteomyelitis of his T7 and T8 vertebral bodies. The patient’s isolate from this time point was found to have a vancomycin MIC of 3.0 µg/ml, initially by automated susceptibility testing, and then determined by Etest on Mueller-Hinton agar, which is considered to be at the lower end of the VISA MIC range (13). The patient was therefore treated with daptomycin after this bacteremia. Patient isolates from three time points (January, April, and November; Table 1) were preserved and subsequently treated for DNA extraction (see Materials and Methods).

View this table:
  • View inline
  • View popup
TABLE 1

Patient’s S. aureus isolate sites, phenotypes, and genome sequences

Isolate genome sequences and phylogeny.Patient S. aureus isolate genome sequences from the three time points were sequenced to ~500× coverage as described in Materials and Methods (Table 1). De novo assembly was performed in order to yield sets of consensus sequence contigs for each of the three sets of sequence reads. Sequences of the resulting three unfinished assemblies were compared against all complete S. aureus genome sequences available in the NCBI RefSeq database, along with a set of previously published S. aureus genome sequences from clinical isolates taken at the same Emory Healthcare laboratory (33), to yield all-against-all pairwise nucleotide distances. Phylogenetic analysis of these distance data was used to choose the closest complete S. aureus genome sequence for subsequent read-to-genome mapping and single nucleotide polymorphism (SNP) calling.

S. aureus genome sequences from the patient isolates (VSSA-T1, VSSA-T2, and VISA-T3 where T1, T2, and T3 stand for time points 1, 2, and 3) form a single monophyletic cluster and are most closely related to two of the previously characterized clinical isolates (Fig. 2). The most closely related complete S. aureus reference sequences correspond to a pair of genomes from a previous study (NCBI RefSeq accession no. NC_009487 and NC_009632 ), which also documented and traced the emergence of VISA over time in a single patient (32). The tight phylogenetic clustering seen for the three time point isolates characterized here is consistent with a single infection followed by in situ evolution of the patient’s strains leading to reduced vancomycin susceptibility, as opposed to multiple infections with distinct S. aureus strains having different vancomycin susceptibilities. The phylogenetic analysis does however show substantial variation among the patient’s three time point sequences, far more than seen for the two strains from the previous study, which were collected 3 months apart in time. Furthermore, the VISA isolate sequence from time point three diverges markedly compared to the sequences from the first two time points. This suggests the possibility of a particularly rapid period of evolution, after an initial slower phase, leading to the VISA phenotype at the third time point. However, it likely also reflects, to some extent, the fact that there are only 3 months between time points one and two compared to 7 months between time points two and three.

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Phylogeny of S. aureus isolates based on genome-wide average nucleotide distances (scale bar). Patient isolates from three time points are shown as VSSA-T1, VSSA-T2, and VISA-T3.

Overall reduction in variation.Having defined the closest complete S. aureus reference sequences to the patient’s three isolate sequences, individual reads for each time point isolate were mapped to the reference sequences. As the two reference sequences are highly similar (>99.99% identity), all results based on the read-to-genome mapping were qualitatively identical for each reference sequence. Here, results for the NC_009487 reference sequence are reported.

The high sequencing depth achieved for each time point isolate (Table 1) provided substantial resolution to search for genome sequence variation within individual isolate samples. In other words, individual patient isolates could be evaluated to assess their clonality or lack thereof. To do this, a technique developed to evaluate the clonality of cancer samples was applied to the isolate genome sequence reads (37). This approach treats individual reads and their site-specific base calls as alleles in a population. The variant allele frequencies, i.e., the frequency of alleles that differ from the reference sequence, are compared to the sequencing depth in order to identify colocated clusters that are likely to correspond to individual clones within a sample of cells. Applying this technique to the patient’s time point genome sequence reads revealed that the patient isolates are polyclonal and that there was a monotonic reduction in the total number of clones from time point one to time point three (Fig. 3). This is consistent with an overall reduction of variation over time and can be attributed to the bottleneck that bacterial populations experience when exposed to antibiotic treatment.

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Clonality of S. aureus patient isolates from three time points. Clusters of variants with similar variant allele frequencies (VAFs) characterize individual clones. The number of clones, along with their average VAFs, are shown for each time point isolate.

The overall levels of variation within each time point’s isolate population were further quantified using distributions of both the site-specific variant allele frequencies and site-specific heterozygosity values (see Materials and Methods). Both of these metrics reveal a monotonic reduction in the overall sequence variation across the three time points, consistent with the results from the clonality analysis (Fig. 4). There is a particularly marked reduction in variation seen for the VISA isolate from time point three. Differences in the levels of variation between the time points are all highly significant (P ≈ 0 by Mann-Whitney U test).

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Levels of S. aureus genome variation seen for patient isolates from three time points. Box plots are shown for site-specific distributions of variant allele frequency and heterozygosity. For a given mapped position in a time point isolate, the variant allele frequency is defined as the number of nonreference bases divided by the total number of mapped bases. Heterozygosity is a measure of allelic diversity at each of the mapped positions (see Materials and Methods).

Allele frequency spectrum shift.The allele frequency spectra for the patient S. aureus isolate sequences were characterized in order to further interrogate the population evolutionary dynamics that accompanied the decrease in vancomycin susceptibility over time. To do this, variant allele counts for the three time points were calculated across an increasingly stringent series of variant allele frequency cutoff levels (Fig. 5). At low cutoff levels, all or most of the variant sites are considered, whereas at higher cutoff levels, only variant alleles that are found at high frequency are counted. At the 0 cutoff level, which is equivalent to Fig. 4B, all variant alleles are considered. At this level, and for other low cutoff levels (0 to 0.2), time point one shows the most variant alleles, followed closely by time point two; then, there is a precipitous drop in the number of variant alleles at time point three. This is consistent with the overall reduction in variation across time points documented in the previous section. As the variant allele frequency cutoff increases, this pattern begins to shift. At higher cutoff levels, i.e., when only high-frequency variants are counted (>0.25), time point three shows the highest number of variants, followed by time points two and one, respectively (Fig. 5). This result indicates that, despite the overall reduction in variation over time that was caused by antibiotic treatment, a number of low-frequency variants were swept to high frequency levels (or fixed) in the VISA time point population.

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

S. aureus variant allele count distribution for patient isolates from three time points. Variant allele counts are shown across a range of variant allele frequency cutoffs (top panel). Insets are shown for the two most extreme cutoff values (bottom panel).

These kinds of VISA-characteristic mutations are likely to point to genomic loci (genes) that are implicated in the emergence of the VISA phenotype. Furthermore, the overall reduction in variation over time, i.e., the population bottleneck, provides a background against which VISA site-specific allele frequency shifts can be parameterized to look for VISA-characteristic mutations that are statistically significant outliers. Two statistical tests based on this logic were applied to the time point-specific allelic variant data to look for sites that are most likely to be related to the emergence of VISA in this patient. The first test is based on the distribution of site-specific variant allele frequency differences, and the second test is based on the distribution of site-specific population differentiation levels (polarized FST). All ~2 million genomic positions are considered in these tests, yielding high statistical power. These tests are further described in Materials and Methods.

The variant allele frequency difference distribution and the FST were computed for the vancomycin-susceptible S. aureus isolate at time point one (VSSA-T1) versus the VISA time point and for the VSSA-T2 versus VISA time points. The results are qualitatively similar, and results of the comparisons between the VSSA-T1 and VISA time points are reported here (Fig. 6). The overall variant allele frequency difference distribution is shifted toward VSSA, consistent with the overall reduction in variation over time. However, the VISA-specific allele frequency differences have a broader distribution with a long tail, consistent with the allele frequency distribution shift previously described. The statistical significance values of the site-specific allele frequency differences were used to rank VISA-implicated mutations for subsequent interrogation (Fig. 6).

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

Site-specific values of S. aureus variant allele frequency differences and variant allele population differentiation between VSSA and VISA patient isolates. (Left) Distributions showing site-specific values of variant allele frequency change and population differentiation (FST). Site-specific values are colored according to whether they are greater in VSSA (blue) or VISA (red). (Right) Manhattan plots showing the corresponding statistical significance levels for site-specific changes in variant allele frequencies and variant allele population differentiation in the VISA isolate.

The VSSA-to-VISA FST distribution shows a different pattern with the majority of site-specific polarized FST values being VISA related. Nevertheless, the VISA side of the distribution is broader and shows a longer tail as seen for the variant allele frequency difference distribution. The statistical significance values of the site-specific polarized FST values were also used to rank VISA-implicated mutations for subsequent interrogation (Fig. 6).

Candidate gene VISA mutations.The utility of the population genomic approach employed here for the detection of VISA-related mutations was evaluated by searching for statistically significant mutations in candidate genes that had previously been implicated in the VISA phenotype. To do this, a database of 24 VISA-implicated candidate genes based on the results of 17 previous studies (reviewed in reference 19) was created. Of these 24 genes, 12 bear statistically significant nonsynonymous mutations (Table 2; see Table S1 in the supplemental material); furthermore, the total number of significant mutations found in the candidate genes is higher than the expected number based on the distribution of significant mutations across the entire genome (P = 8.1 × 10−9 by Fisher exact test). A similar approach was used to interrogate a larger set of 237 genes, which had previously been implicated in the VISA phenotype by virtue of their differential expression in VISA compared to related VSSA strains (26). Of the 237 genes, 145 bear statistically significant nonsynonymous mutations, and there is also a significant increase in the total number of significant mutations that map to these genes (P = 2.8 × 10−3 by Fisher exact test). These results serve as a proof of principle for the population genomic approach employed here and underscore the power of the statistical tests employed.

Table S1

List of previously reported vancomycin resistance mutations detected in the isolates analyzed here. For each mutation listed, the genomic position in the reference genome NC_009487 , the nucleotide change, the NCBI identifier (ID), the gene name, the function of the gene that the mutation is present in, the amino acid change, and the Grantham score of the effect of mutation are provided. Download Table S1, XLSX file, 0.01 MB.
Copyright © 2016 Rishishwar et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

View this table:
  • View inline
  • View popup
TABLE 2

Genes reported in the literature to be implicated in the VISA phenotype

Novel gene VISA mutations.Having established the relevance of the statistically significant VISA-related mutations by virtue of their overabundance in previously identified VISA candidate genes, the presence of such mutations in novel genes not previously implicated in the VISA phenotype was interrogated for clues as to potentially novel mechanisms that underlie the acquisition of VISA in this patient. To do this, mutations were ranked by their statistical significance, and functional enrichment analysis was performed on genes corresponding to the top 500 SNPs. This was done independently for each test, variant allele frequency difference and FST, and for both tests combined.

One particularly interesting set of genes, which turned up as enriched in all three instances of this approach, encodes the LPXTG cell wall anchor domain protein family. Genes encoding 7 out of 20 members of this family contain statistically significant mutations, far more than is expected by chance (Fig. 7A and B). Members of this family are surface proteins that perform a variety of functions related to pathogen-host interactions, including antigen presentation, induction of platelet aggregation, and maintenance of cell wall integrity (38–40). These proteins all contain an LPXTG sequence motif at their C terminus, which is a cleavage signal that leads to the covalent binding of the proteins to the cell wall (39, 40) (Fig. 7C). The presence of covalently bound LPXTG proteins has been shown to contribute to the regular cross-linked structure of the cell wall peptidoglycan layer (41).

FIG 7
  • Open in new tab
  • Download powerpoint
FIG 7

The LPXTG cell wall anchor domain protein family. (A) S. aureus genomic distribution of LPXTG motif-containing protein-coding genes. (B) Observed versus expected numbers of LPXTG motif-containing protein-coding genes that overlap with the top 500 statistically significant mutations (Fig. 6). (C) Schematic representation of cell-wall-related LPXTG motif-containing protein functions.

We evaluated the domain architectures of the seven LPXTG proteins implicated by the three statistical tests in order to better understand the functional characteristics that are shared among them (Fig. 8). Each protein can be seen to contain the YSIRK type signal peptide at the N-terminal domain, followed by low-complexity regions (black lines) and, in most cases, an SdrG glycoprotein adhesion domain. Some of the LPXTG proteins contain additional domains such as carboxypeptidase regulatory-like and/or Cna type B domains. The proteins typically end with another low-complexity region followed by a Gram-positive anchor (shown at the right of the diagram) at the C terminus, which contains the characteristic LPXTG domain. Transmembrane segments are found to be colocated with N- and C-terminal domains. All of these domains point to a functional role for these proteins related to the peptidoglycan cell wall layer.

FIG 8
  • Open in new tab
  • Download powerpoint
FIG 8

Domain architectures of LPXTG cell wall anchor domain protein family members encoded by genes with VISA-implicated mutations. The NCBI accession number, length (in amino acids [aa]), gene symbol, and domain architecture of each protein are shown. The domains are color coded and shown along with their names. The transmembrane domains (TM), carboxypeptidase regulatory-like (CarboxypepD reg), Gram-positive anchor (Gram_pos_anchor), and the other domain names are described in the text. Black lines above each protein diagram correspond to low-complexity regions.

Mutations to members of this family, such as those seen in the time point three isolate population for this patient, could function to reduce the cross-linking of the peptidoglycan cell wall layer. This is relevant to the VISA phenotype, since the extent to which vancomycin can permeate the S. aureus cell wall has previously been shown to be related to the integrity of peptidoglycan cross-linking (42). Specifically, aberrations in the cross-linked peptidoglycan structure could provide false binding sites for vancomycin, thereby reducing its permeation across the cell wall layer. In other words, mutations to members of the LPXTG cell wall anchor domain protein family may result in a lowering of the effective concentration of vancomycin and accordingly reduced susceptibility.

Previous results showing a relationship between the VISA phenotype and an overall reduction in the expression of surface proteins (43) points to another potential mechanism by which mutations to LPXTG cell wall anchor domain proteins could yield reduced vancomycin susceptibility. Since LPXTG motifs are responsible for mediating the proper translocation and anchoring of surface proteins to the peptidoglycan layer, mutations to these proteins could lead to lower overall levels of surface proteins and thus be functionally analogous to the previously observed lower levels of surface protein expression in VISA strains.

A case of heterogeneous VISA.The results of the population genomic comparison between VSSA and VISA isolates from the same patient reveal this to be a likely case of heterogeneous VISA (hVISA). hVISA isolates are polyclonal populations wherein the entire isolate population shows a vancomycin MIC within the susceptible range (<2 µg/ml) but a small proportion of clones from the same isolate show an MIC in the intermediate range (3 to 8 µg/ml) (35). In the case of the patient case studied here, the entire isolate population from time point three shows an MIC at the low end of the intermediate range (3 µg/ml). Nevertheless, there are likely to be a number of clones for this same time point that have higher MIC values. This is supported by the fact that most of the statistically significant mutations detected by the population genomic tests employed here have been swept to relatively high frequency but have not been fixed. Furthermore, the patterns of allelic change observed here indicate that there are likely to be clones from the time point one and time point two VSSA isolates, characterized by relatively low-frequency mutations, which also have MICs in the intermediate range.

hVISA is typically detected using a culture-based test referred to as a vancomycin population analysis profile (PAP) (44). This approach is both labor-intensive and time intensive, since it requires a series of continuous MIC readings at different concentrations over 2 to 3 days. The results here point to the possibility of using a population genomic approach, which is made possible by the increasing throughput and decreasing costs of next-generation sequencing techniques, for the characterization of hVISA. Adaptation and standardization of such a genome sequence-based approach to the detection of hVISA could be particularly relevant, given the fact that the retrospective studies have revealed 50% of VSSA isolates to represent cases of hVISA (34). The utility of the population genomic approach is discussed further below.

Utility of the population genomic approach.A novel population genomic approach was taken for this study on the genomic basis of the emergence of the VISA phenotype in a single patient over time. Typically, studies of this kind have relied on the comparison of pairs of single genome sequences, each corresponding to a particular level of vancomycin susceptibility, which are in reality consensus sequences generated from the merging of numerous sequence reads covering each position of the genome. This approach assumes that patient isolates are monoclonal, and as such, any observed site-specific variation is considered to be noise that is removed in the process of generating a single genomic consensus sequence.

The population genomic approach employed here relies on deep genomic sequencing that yields high sequence read coverage across the genome (Fig. 1 and Table 1). This high coverage can be used to detect bona fide sequence variation at specific sites that may indicate the presence of more than one clone in a patient isolate. Site-specific variations uncovered by this approach can be considered allelic differences between multiple clones within time points and can also be used to model the population genome dynamics of acquired antibiotic resistance over time. Much of the power of this approach rests on the fact that it considers site-specific variation levels genome-wide, yielding ~2 million data points, which can be used to create distributions of variation that can parameterize statistical tests as shown here (Fig. 6).

Comparison of results based on the more-traditional consensus sequence-based approach with the novel tests employed here underscores the power of the population genomic approach and the extent to which distinct statistical tests are complementary. The genome consensus sequence approach reveals only 12 differences between the VSSA isolate from time point one and the VISA isolate from time point three, whereas the variant allele frequency difference and FST tests yield numerous statistically significant results. The large number of significant results returned by these tests could be taken to suggest that it is overly sensitive, but there is a simple solution to this potential problem that was used here, which is known as the outlier approach. In the outlier approach, mutations are ranked accordingly to their statistical significance, and only the most significant mutations are considered. In addition to increasing the stringency of the statistical tests employed, this approach also allows for changing of the threshold for consideration in order to assess the robustness of any functional enrichment analysis performed on the results.

MATERIALS AND METHODS

Patient isolates and susceptibility testing.Patient bloodstream samples at three time points (Table 1) were collected and cultured in the Emory Healthcare microbiology laboratory (under Institutional Review Board approval 50685). For clinical care of the patient, the S. aureus isolates were subjected to automated susceptibility testing using the MicroScan WalkAway 96 Plus system (Siemens Healthcare Diagnostics Inc., Tarrytown, NY), and vancomycin susceptibility was confirmed by Etest (bioMérieux, Inc., Durham, NC). The initial subculture plate was entirely swept using a plastic loop (after at least 48 h of incubation) and frozen at −80°C for DNA isolation in order to ensure the representation of isolate subpopulations for subsequent deep sequencing and analysis.

DNA isolation and genome sequencing.Each time point isolate colony set was incubated with 50 mg/ml lysozyme and 5 mg/ml lysostaphin at 37°C for 1 h prior to DNA extraction. DNA extraction was performed with the EZ1 Advanced XL system using the DNA Tissue kit and the Bacteria Card (Qiagen, Valencia, CA). Time point isolate DNA sequencing was performed on the Illumina MiSeq instrument using the paired-end read protocol. This resulted in high-coverage (~500×) draft genome sequences used for subsequent analyses.

Genome sequence analysis.Genome sequence reads were evaluated for quality using the FastQC program (45). The mean Phred score for all reads was 35, and a Phred cutoff score of 20 was used for removal of low-quality reads and read trimming. Removal of low-quality reads and read trimming were done with the FASTX-toolkit (46). De novo genome assemblies for the three time point isolates were performed using the ABySS program (47). Isolate assemblies were compared to 14 complete S. aureus genome sequences taken from the NCBI RefSeq database (48), along with 5 S. aureus clinical isolate genome sequences previously characterized by Emory Healthcare (33), using the MUMmer program (49). The all-against-all pairwise average nucleotide identities computed with MUMmer were converted to p-distances and used to reconstruct an S. aureus whole-genome phylogeny with the neighbor-joining algorithm (50) implemented in the MEGA program (51). This phylogeny was used to select the most closely related complete S. aureus genome sequence(s) for subsequent population genomic analyses.

Population genomic analyses.Individual sequence reads from each of the three time point isolates were mapped to S. aureus reference sequence(s) using the BWA program (52). Individual nucleotide variants (SNPs) were called on the basis of the resulting read-to-genome alignments using SAMtools (53) and GATK (54). Results from these two programs were virtually identical, and the results from SAMtools are presented here. Site-specific variant calls were used to model allelic variation within each time point isolate population. Quality control measures based on the mapping quality and read depth were used to ensure the reliability of the site-specific variant calls used to represent allelic differences. The mean BWA (Phred-like) mapping quality score was 56.4, and the mean read depth was 481. A lower bound depth cutoff of 98 reads per position (corresponding to 1.2 standard deviations from the mean) was used for calculating site-specific allelic diversity. A number of population genomic parameters were calculated on the basis of the site-specific allelic diversity values as described below.

Site-specific variant allele frequencies (VAFi) are computed as VAFi = vi/ni where vi is the number of variant nucleotides at site i (i.e., nucleotides that differ from the consensus sequence) and ni is the total number of aligned nucleotides at that site. VAFi values were used together with site-specific read depth values in order to compute the number of clones present in each of the three time point isolate populations using a kernel density method implemented in the SciClone program (37). Site-specific heterozygosity (Hi) values were computed as Hi=1−∑ATCGp2 where p is the frequency of each nucleotide at site i.

Population genomic differentiation between time point isolates was computed using site-specific variant allele frequency differences (ΔVAFi) and polarized site-specific fixation indices (FSTi). The site-specific variant allele frequency differences were calculated as follows: ΔVAFi = t3VAFi − txVAFi where t3 is time point three (VISA isolate) and tx can be time point one or time point two (VSSA isolate). The polarized site-specific fixation indices were calculated as follows: FSTi = (HTi − HSi)/HTi where HTi is the heterozygosity computed with both time points considered as a single metapopulation, and HSi is the heterozygosity computed individually for both time point isolate subpopulations. The FSTi values were polarized by assigning a negative value to sites with variant allele frequencies higher at time point one or two and a positive value to sites with variant allele frequencies higher at time point three. ΔVAF and polarized FSTi distributions were parameterized in order to identify statistically significant outliers using a Z test. A Bonferroni’s correction, based on the approximate number of sites analyzed genome-wide (2 million), was used to compute a P value significance threshold of 10−8.

Functional enrichment analysis.The genomic distributions of variant nucleotide sites with statistically significant time point isolate population differentiation values (for ΔVAFi, polarized FSTi, or both) were evaluated with respect to their presence in a set of 24 candidate genes previously implicated in the VISA phenotype (19). To do this, a Fisher exact test was used to compare the observed number of such sites that overlap candidate genes compared to the expected number of sites computed on the basis of their genomic background density. An outlier approach was used to search for functional enrichment of novel genes not previously implicated in the VISA phenotype. To do this, genes that overlap with the top 500 most statistically significant sites for each test were analyzed. These genes were expected to identify functionally coherent gene sets (families), and enrichment values for these families were manually computed using Fisher exact test in the same way as described above.

Protein domain architecture characterization.Protein domain architectures for members of the LPXTG cell wall anchor domain protein family were characterized using the SMART tool with both the default SMART and Pfam domain databases (55, 56). The protein domain architecture diagrams were drawn using the MyDomains tool from ExPASy PROSITE (57).

Ethics statement.This study and patient confidentiality were covered under the Emory Institutional Review Board (IRB) approval (approval no. 50685). Patient isolates were collected under this study as discarded clinical samples based on the IRB approval to perform retrospective chart review and analysis.

Accession numbers.The raw sequence reads of the three time point isolates can be accessed from NCBI SRA (Table 1). The SRA accession numbers of the three time point isolates VSSA-T1, VSSA-T2, and VISA-T3 are SRX689719 , SRX689725 , and SRX689726 , respectively.

FOOTNOTES

    • Received April 8, 2016.
    • Accepted June 23, 2016.
  • Copyright © 2016 Rishishwar et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license .

REFERENCES

  1. 1.↵
    1. Lowy FD
    . 1998. Staphylococcus aureus infections. N Engl J Med339:520–532. doi:10.1056/NEJM199808203390806.
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.↵
    1. Laupland KB,
    2. Ross T,
    3. Gregson DB
    . 2008. Staphylococcus aureus bloodstream infections: risk factors, outcomes, and the influence of methicillin resistance in Calgary, Canada, 2000–2006. J Infect Dis198:336–343. doi:10.1086/589717.
    OpenUrlCrossRefPubMedWeb of Science
  3. 3.↵
    1. van Hal SJ,
    2. Jensen SO,
    3. Vaska VL,
    4. Espedido BA,
    5. Paterson DL,
    6. Gosbell IB
    . 2012. Predictors of mortality in Staphylococcus aureus bacteremia. Clin Microbiol Rev25:362–386. doi:10.1128/CMR.05022-11.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Liu C,
    2. Bayer A,
    3. Cosgrove SE,
    4. Daum RS,
    5. Fridkin SK,
    6. Gorwitz RJ,
    7. Kaplan SL,
    8. Karchmer AW,
    9. Levine DP,
    10. Murray BE,
    11. Rybak MJ,
    12. Talan DA,
    13. Chambers HF
    . 2011. Clinical practice guidelines by the Infectious Diseases Society of America for the treatment of methicillin-resistant Staphylococcus aureus infections in adults and children: executive summary. Clin Infect Dis52:285–292. doi:10.1093/cid/cir034.
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.↵
    1. Chambers HF,
    2. Deleo FR
    . 2009. Waves of resistance: Staphylococcus aureus in the antibiotic era. Nat Rev Microbiol7:629–641. doi:10.1038/nrmicro2200.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    1. Barber M
    . 1961. Methicillin-resistant staphylococci. J Clin Pathol14:385–393. doi:10.1136/jcp.14.4.385.
    OpenUrlAbstract/FREE Full Text
  7. 7.↵
    1. Jevons MP
    . 1961. “Celbemom”-resistant staphylococci. Br Med J1:113–114.
    OpenUrlFREE Full Text
  8. 8.↵
    1. Otto M
    . 2012. MRSA virulence and spread. Cell Microbiol14:1513–1521. doi:10.1111/j.1462-5822.2012.01832.x.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Beck WD,
    2. Berger-Bächi B,
    3. Kayser FH
    . 1986. Additional DNA in methicillin-resistant Staphylococcus aureus and molecular cloning of mec-specific DNA. J Bacteriol165:373–378.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    1. Katayama Y,
    2. Ito T,
    3. Hiramatsu K
    . 2000. A new class of genetic element, staphylococcus cassette chromosome mec, encodes methicillin resistance in Staphylococcus aureus. Antimicrob Agents Chemother44:1549–1555. doi:10.1128/AAC.44.6.1549-1555.2000.
    OpenUrlAbstract/FREE Full Text
  11. 11.↵
    1. Pope SD,
    2. Roecker AM
    . 2007. Vancomycin for treatment of invasive, multi-drug resistant Staphylococcus aureus infections. Expert Opin Pharmacother8:1245–1261. doi:10.1517/14656566.8.9.1245.
    OpenUrlCrossRefPubMedWeb of Science
  12. 12.↵
    1. Tenover FC,
    2. Biddle JW,
    3. Lancaster MV
    . 2001. Increasing resistance to vancomycin and other glycopeptides in Staphylococcus aureus. Emerg Infect Dis7:327–332. doi:10.3201/eid0702.700327.
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    1. Tenover FC,
    2. Moellering RC Jr.
    2007. The rationale for revising the Clinical and Laboratory Standards Institute vancomycin minimal inhibitory concentration interpretive criteria for Staphylococcus aureus. Clin Infect Dis44:1208–1215. doi:10.1086/513203.
    OpenUrlCrossRefPubMedWeb of Science
  14. 14.↵
    1. Hiramatsu K,
    2. Hanaki H,
    3. Ino T,
    4. Yabuta K,
    5. Oguri T,
    6. Tenover FC
    . 1997. Methicillin-resistant Staphylococcus aureus clinical strain with reduced vancomycin susceptibility. J Antimicrob Chemother40:135–136.
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.↵
    1. Centers for Disease Control and Prevention
    . 2002. Staphylococcus aureus resistant to vancomycin--United States, 2002. MMWR Morb Mortal Wkly Rep51:565–567.
    OpenUrlPubMed
  16. 16.↵
    1. Arthur M,
    2. Molinas C,
    3. Depardieu F,
    4. Courvalin P
    . 1993. Characterization of Tn1546, a Tn3-related transposon conferring glycopeptide resistance by synthesis of depsipeptide peptidoglycan precursors in Enterococcus faecium BM4147. J Bacteriol175:117–127.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    1. Noble WC,
    2. Virani Z,
    3. Cree RG
    . 1992. Co-transfer of vancomycin and other resistance genes from Enterococcus faecalis NCTC 12201 to Staphylococcus aureus. FEMS Microbiol Lett72:195–198.
    OpenUrlPubMed
  18. 18.↵
    1. Dhand A,
    2. Sakoulas G
    . 2012. Reduced vancomycin susceptibility among clinical Staphylococcus aureus isolates (‘the MIC Creep’): implications for therapy. F1000 Med Rep4:4. doi:10.3410/M4-4.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Howden BP,
    2. Davies JK,
    3. Johnson PD,
    4. Stinear TP,
    5. Grayson ML
    . 2010. Reduced vancomycin susceptibility in Staphylococcus aureus, including vancomycin-intermediate and heterogeneous vancomycin-intermediate strains: resistance mechanisms, laboratory detection, and clinical implications. Clin Microbiol Rev23:99–139. doi:10.1128/CMR.00042-09.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Trotonda MP,
    2. Xiong YQ,
    3. Memmi G,
    4. Bayer AS,
    5. Cheung AL
    . 2009. Role of mgrA and sarA in methicillin-resistant Staphylococcus aureus autolysis and resistance to cell wall-active antibiotics. J Infect Dis199:209–218. doi:10.1086/595740.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Neoh HM,
    2. Cui L,
    3. Yuzawa H,
    4. Takeuchi F,
    5. Matsuo M,
    6. Hiramatsu K
    . 2008. Mutated response regulator graR is responsible for phenotypic conversion of Staphylococcus aureus from heterogeneous vancomycin-intermediate resistance to vancomycin-intermediate resistance. Antimicrob Agents Chemother52:45–53. doi:10.1128/AAC.00534-07.
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    1. Cui L,
    2. Neoh HM,
    3. Shoji M,
    4. Hiramatsu K
    . 2009. Contribution of vraSR and graSR point mutations to vancomycin resistance in vancomycin-intermediate Staphylococcus aureus. Antimicrob Agents Chemother53:1231–1234. doi:10.1128/AAC.01173-08.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    1. Holmes NE,
    2. Turnidge JD,
    3. Munckhof WJ,
    4. Robinson JO,
    5. Korman TM,
    6. O’Sullivan MV,
    7. Anderson TL,
    8. Roberts SA,
    9. Warren SJ,
    10. Coombs GW,
    11. Tan HL,
    12. Gao W,
    13. Johnson PD,
    14. Howden BP
    . 2014. Genetic and molecular predictors of high vancomycin MIC in Staphylococcus aureus bacteremia isolates. J Clin Microbiol52:3384–3393. doi:10.1128/JCM.01320-14.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. Schulthess B,
    2. Meier S,
    3. Homerova D,
    4. Goerke C,
    5. Wolz C,
    6. Kormanec J,
    7. Berger-Bächi B,
    8. Bischoff M
    . 2009. Functional characterization of the sigmaB-dependent yabJ-spoVG operon in Staphylococcus aureus: role in methicillin and glycopeptide resistance. Antimicrob Agents Chemother53:1832–1839. doi:10.1128/AAC.01255-08.
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    1. Howden BP,
    2. Stinear TP,
    3. Allen DL,
    4. Johnson PD,
    5. Ward PB,
    6. Davies JK
    . 2008. Genomic analysis reveals a point mutation in the two-component sensor gene graS that leads to intermediate vancomycin resistance in clinical Staphylococcus aureus. Antimicrob Agents Chemother52:3755–3762. doi:10.1128/AAC.01613-07.
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    1. Howden BP,
    2. Smith DJ,
    3. Mansell A,
    4. Johnson PD,
    5. Ward PB,
    6. Stinear TP,
    7. Davies JK
    . 2008. Different bacterial gene expression patterns and attenuated host immune responses are associated with the evolution of low-level vancomycin resistance during persistent methicillin-resistant Staphylococcus aureus bacteraemia. BMC Microbiol8:39. doi:10.1186/1471-2180-8-39.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. Katayama Y,
    2. Murakami-Kuroda H,
    3. Cui L,
    4. Hiramatsu K
    . 2009. Selection of heterogeneous vancomycin-intermediate Staphylococcus aureus by imipenem. Antimicrob Agents Chemother53:3190–3196. doi:10.1128/AAC.00834-08.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. Matsuo M,
    2. Hishinuma T,
    3. Katayama Y,
    4. Cui L,
    5. Kapi M,
    6. Hiramatsu K
    . 2011. Mutation of RNA polymerase beta subunit (rpoB) promotes hVISA-to-VISA phenotypic conversion of strain Mu3. Antimicrob Agents Chemother55:4188–4195. doi:10.1128/AAC.00398-11.
    OpenUrlAbstract/FREE Full Text
  29. 29.↵
    1. Howden BP,
    2. McEvoy CR,
    3. Allen DL,
    4. Chua K,
    5. Gao W,
    6. Harrison PF,
    7. Bell J,
    8. Coombs G,
    9. Bennett-Wood V,
    10. Porter JL,
    11. Robins-Browne R,
    12. Davies JK,
    13. Seemann T,
    14. Stinear TP
    . 2011. Evolution of multidrug resistance during Staphylococcus aureus infection involves mutation of the essential two component regulator WalKR. PLoS Pathog7:e1002359. doi:10.1371/journal.ppat.1002359.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Gardete S,
    2. Kim C,
    3. Hartmann BM,
    4. Mwangi M,
    5. Roux CM,
    6. Dunman PM,
    7. Chambers HF,
    8. Tomasz A
    . 2012. Genetic pathway in acquisition and loss of vancomycin resistance in a methicillin resistant Staphylococcus aureus (MRSA) strain of clonal type USA300. PLoS Pathog8:e1002505. doi:10.1371/journal.ppat.1002505.
    OpenUrlCrossRefPubMed
  31. 31.↵
    1. Alam MT,
    2. Petit RA III,
    3. Crispell EK,
    4. Thornton TA,
    5. Conneely KN,
    6. Jiang Y,
    7. Satola SW,
    8. Read TD
    . 2014. Dissecting vancomycin-intermediate resistance in Staphylococcus aureus using genome-wide association. Genome Biol Evol6:1174–1185. doi:10.1093/gbe/evu092.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Mwangi MM,
    2. Wu SW,
    3. Zhou Y,
    4. Sieradzki K,
    5. de Lencastre H,
    6. Richardson P,
    7. Bruce D,
    8. Rubin E,
    9. Myers E,
    10. Siggia ED,
    11. Tomasz A
    . 2007. Tracking the in vivo evolution of multidrug resistance in Staphylococcus aureus by whole-genome sequencing. Proc Natl Acad Sci U S A104:9451–9456. doi:10.1073/pnas.0609839104.
    OpenUrlAbstract/FREE Full Text
  33. 33.↵
    1. Rishishwar L,
    2. Petit RA III,
    3. Kraft CS,
    4. Jordan IK
    . 2014. Genome sequence-based discriminator for vancomycin-intermediate Staphylococcus aureus. J Bacteriol196:940–948. doi:10.1128/JB.01410-13.
    OpenUrlAbstract/FREE Full Text
  34. 34.↵
    1. Horne KC,
    2. Howden BP,
    3. Grabsch EA,
    4. Graham M,
    5. Ward PB,
    6. Xie S,
    7. Mayall BC,
    8. Johnson PD,
    9. Grayson ML
    . 2009. Prospective comparison of the clinical impacts of heterogeneous vancomycin-intermediate methicillin-resistant Staphylococcus aureus (MRSA) and vancomycin-susceptible MRSA. Antimicrob Agents Chemother53:3447–3452. doi:10.1128/AAC.01365-08.
    OpenUrlAbstract/FREE Full Text
  35. 35.↵
    1. Liu C,
    2. Chambers HF
    . 2003. Staphylococcus aureus with heterogeneous resistance to vancomycin: epidemiology, clinical significance, and critical assessment of diagnostic methods. Antimicrob Agents Chemother47:3040–3045. doi:10.1128/AAC.47.10.3040-3045.2003.
    OpenUrlFREE Full Text
  36. 36.↵
    1. Howden BP,
    2. Peleg AY,
    3. Stinear TP
    . 2014. The evolution of vancomycin intermediate Staphylococcus aureus (VISA) and heterogenous-VISA. Infect Genet Evol21:575–582. doi:10.1016/j.meegid.2013.03.047.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Miller CA,
    2. White BS,
    3. Dees ND,
    4. Griffith M,
    5. Welch JS,
    6. Griffith OL,
    7. Vij R,
    8. Tomasson MH,
    9. Graubert TA,
    10. Walter MJ,
    11. Ellis MJ,
    12. Schierding W,
    13. DiPersio JF,
    14. Ley TJ,
    15. Mardis ER,
    16. Wilson RK,
    17. Ding L
    . 2014. SciClone: inferring clonal architecture and tracking the spatial and temporal patterns of tumor evolution. PLoS Comput Biol10:e1003665. doi:10.1371/journal.pcbi.1003665.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Hendrickx AP,
    2. Budzik JM,
    3. Oh SY,
    4. Schneewind O
    . 2011. Architects at the bacterial surface—sortases and the assembly of pili with isopeptide bonds. Nat Rev Microbiol9:166–176. doi:10.1038/nrmicro2520.
    OpenUrlCrossRefPubMedWeb of Science
  39. 39.↵
    1. Marraffini LA,
    2. Dedent AC,
    3. Schneewind O
    . 2006. Sortases and the art of anchoring proteins to the envelopes of Gram-positive bacteria. Microbiol Mol Biol Rev70:192–221. doi:10.1128/MMBR.70.1.192-221.2006.
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    1. Navarre WW,
    2. Schneewind O
    . 1999. Surface proteins of Gram-positive bacteria and mechanisms of their targeting to the cell wall envelope. Microbiol Mol Biol Rev63:174–229.
    OpenUrlAbstract/FREE Full Text
  41. 41.↵
    1. Mazmanian SK,
    2. Liu G,
    3. Ton-That H,
    4. Schneewind O
    . 1999. Staphylococcus aureus sortase, an enzyme that anchors surface proteins to the cell wall. Science285:760–763. doi:10.1126/science.285.5428.760.
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    1. Ton-That H,
    2. Schneewind O
    . 1999. Anchor structure of staphylococcal surface proteins. IV. Inhibitors of the cell wall sorting reaction. J Biol Chem274:24316–24320. doi:10.1074/jbc.274.34.24316.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    1. McAleese F,
    2. Wu SW,
    3. Sieradzki K,
    4. Dunman P,
    5. Murphy E,
    6. Projan S,
    7. Tomasz A
    . 2006. Overexpression of genes of the cell wall stimulon in clinical isolates of Staphylococcus aureus exhibiting vancomycin-intermediate-S. aureus-type resistance to vancomycin. J Bacteriol188:1120–1133. doi:10.1128/JB.188.3.1120-1133.2006.
    OpenUrlAbstract/FREE Full Text
  44. 44.↵
    1. Wootton M,
    2. Howe RA,
    3. Hillman R,
    4. Walsh TR,
    5. Bennett PM,
    6. MacGowan AP
    . 2001. A modified population analysis profile (PAP) method to detect hetero-resistance to vancomycin in Staphylococcus aureus in a UK hospital. J Antimicrob Chemother47:399–403. doi:10.1093/jac/47.4.399.
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    1. Andrews S
    . 2012. FastQC. A quality control tool for high throughput sequence data. Bioinformatics Group, Babraham Institute, Cambridge, United Kingdom. http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc .
  46. 46.↵
    1. Gordon A,
    2. Hannon G
    . 2010. FASTX-toolkit. FASTQ/A short-reads pre-processing tools. Cold Spring Harbor Laboratory, Cold Spring Harbor, NY. http://hannonlab.cshl.edu/fastx_toolkit .
  47. 47.↵
    1. Simpson JT,
    2. Wong K,
    3. Jackman SD,
    4. Schein JE,
    5. Jones SJ,
    6. Birol I
    . 2009. ABySS: a parallel assembler for short read sequence data. Genome Res19:1117–1123. doi:10.1101/gr.089532.108.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    1. Pruitt KD,
    2. Tatusova T,
    3. Maglott DR
    . 2007. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res35:D61–D65. doi:10.1093/nar/gkl842.
    OpenUrlCrossRefPubMedWeb of Science
  49. 49.↵
    1. Kurtz S,
    2. Phillippy A,
    3. Delcher AL,
    4. Smoot M,
    5. Shumway M,
    6. Antonescu C,
    7. Salzberg SL
    . 2004. Versatile and open software for comparing large genomes. Genome Biol5:R12. doi:10.1186/gb-2004-5-2-r12.
    OpenUrlCrossRefPubMed
  50. 50.↵
    1. Saitou N,
    2. Nei M
    . 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol4:406–425.
    OpenUrlCrossRefPubMedWeb of Science
  51. 51.↵
    1. Tamura K,
    2. Stecher G,
    3. Peterson D,
    4. Filipski A,
    5. Kumar S
    . 2013. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol30:2725–2729. doi:10.1093/molbev/mst197.
    OpenUrlCrossRefPubMedWeb of Science
  52. 52.↵
    1. Li H,
    2. Durbin R
    . 2010. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics26:589–595. doi:10.1093/bioinformatics/btp698.
    OpenUrlCrossRefPubMedWeb of Science
  53. 53.↵
    1. Li H,
    2. Handsaker B,
    3. Wysoker A,
    4. Fennell T,
    5. Ruan J,
    6. Homer N,
    7. Marth G,
    8. Abecasis G,
    9. Durbin R,
    10. Genome Project Data Processing Subgroup
    . 2009. The Sequence Alignment/Map format and SAMtools. Bioinformatics25:2078–2079. doi:10.1093/bioinformatics/btp352.
    OpenUrlCrossRefPubMedWeb of Science
  54. 54.↵
    1. DePristo MA,
    2. Banks E,
    3. Poplin R,
    4. Garimella KV,
    5. Maguire JR,
    6. Hartl C,
    7. Philippakis AA,
    8. del Angel G,
    9. Rivas MA,
    10. Hanna M,
    11. McKenna A,
    12. Fennell TJ,
    13. Kernytsky AM,
    14. Sivachenko AY,
    15. Cibulskis K,
    16. Gabriel SB,
    17. Altshuler D,
    18. Daly MJ
    . 2011. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet43:491–498. doi:10.1038/ng.806.
    OpenUrlCrossRefPubMedWeb of Science
  55. 55.↵
    1. Finn RD,
    2. Coggill P,
    3. Eberhardt RY,
    4. Eddy SR,
    5. Mistry J,
    6. Mitchell AL,
    7. Potter SC,
    8. Punta M,
    9. Qureshi M,
    10. Sangrador-Vegas A,
    11. Salazar GA,
    12. Tate J,
    13. Bateman A
    . 2016. The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res44:D279–D285. doi:10.1093/nar/gkv1344.
    OpenUrlCrossRefPubMed
  56. 56.↵
    1. Schultz J,
    2. Milpetz F,
    3. Bork P,
    4. Ponting CP
    . 1998. SMART, a simple modular architecture research tool: identification of signaling domains. Proc Natl Acad Sci U S A95:5857–5864. doi:10.1073/pnas.95.11.5857.
    OpenUrlAbstract/FREE Full Text
  57. 57.↵
    1. Hulo N,
    2. Bairoch A,
    3. Bulliard V,
    4. Cerutti L,
    5. Cuche BA,
    6. de Castro E,
    7. Lachaize C,
    8. Langendijk-Genevaux PS,
    9. Sigrist CJ
    . 2008. The 20 years of PROSITE. Nucleic Acids Res36:D245–D249. doi:10.1093/nar/gkm977.
    OpenUrlCrossRefPubMedWeb of Science
PreviousNext
Back to top
Download PDF
Citation Tools
Population Genomics of Reduced Vancomycin Susceptibility in Staphylococcus aureus
Lavanya Rishishwar, Colleen S. Kraft, I. King Jordan
mSphere Jul 2016, 1 (4) e00094-16; DOI: 10.1128/mSphere.00094-16

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSphere article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Population Genomics of Reduced Vancomycin Susceptibility in Staphylococcus aureus
(Your Name) has forwarded a page to you from mSphere
(Your Name) thought you would be interested in this article in mSphere.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Population Genomics of Reduced Vancomycin Susceptibility in Staphylococcus aureus
Lavanya Rishishwar, Colleen S. Kraft, I. King Jordan
mSphere Jul 2016, 1 (4) e00094-16; DOI: 10.1128/mSphere.00094-16
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS AND DISCUSSION
    • MATERIALS AND METHODS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

Staphylococcus aureus
antibiotic resistance
genomics
population genetics
vancomycin

Related Articles

Cited By...

About

  • About mSphere
  • Board of Editors
  • Policies
  • For Reviewers
  • For the Media
  • Embargo Policy
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSphereJ

@ASMicrobiology

       

 

Website feedback

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5042