Potential Associations between Severity of Infection and the Presence of Virulence-Associated Genes in Clinical Strains of Staphylococcus aureus

Abstract

Background

The clinical spectrum of Staphylococcus aureus infection ranges from asymptomatic nasal carriage to osteomyelitis, infective endocarditis (IE) and death. In this study, we evaluate potential association between the presence of specific genes in a collection of prospectively characterized S. aureus clinical isolates and clinical outcome.

Methodology/Principal Findings

Two hundred thirty-nine S. aureus isolates (121 methicillin-resistant S. aureus [MRSA] and 118 methicillin-susceptible S. aureus [MSSA]) were screened by array comparative genomic hybridization (aCGH) to identify genes implicated in complicated infections. After adjustment for multiple tests, 226 genes were significantly associated with severity of infection. Of these 226 genes, 185 were not in the SCCmec element. Within the 185 non-SCCmec genes, 171 were less common and 14 more common in the complicated infection group. Among the 41 genes in the SCCmec element, 37 were more common and 4 were less common in the complicated group. A total of 51 of the 2014 sequences evaluated, 14 non-SCCmec and 37 SCCmec, were identified as genes of interest.

Conclusions/Significance

Of the 171 genes less common in complicated infections, 152 are of unknown function and may contribute to attenuation of virulence. The 14 non-SCCmec genes more common in complicated infections include bacteriophage-encoded genes such as regulatory factors and autolysins with potential roles in tissue adhesion or biofilm formation.

About the Authors:

Steven R. Gill

* E-mail: [email protected]

Current address: Department of Microbiology and Immunology, School of Medicine and Dentistry, University of Rochester, Rochester, New York

Affiliations Department of Oral Biology, University at Buffalo, Buffalo, New York, United States of America, Department of Microbiology and Immunology, University at Buffalo, Buffalo, New York, United States of America, NYS Center of Excellence in Bioinformatics and Life Sciences, University at Buffalo, Buffalo, New York, United States of America

Lauren M. McIntyre

Affiliation: Molecular Genetics and Microbiology, University of Florida, Gainesville, Florida, United States of America

Charlotte L. Nelson

Affiliation: Duke Clinical Research Institute, Duke University Medical Center, Durham, North Carolina, United States of America

Brian Remortel

Current address: Association of Public Health Laboratories, Silver Spring, Maryland

Affiliation: J. Craig Venter Institute, Rockville, Maryland, United States of America

Tom Rude

Affiliation: Division of Infectious Diseases, Department of Medicine, Duke University Medical Center, Durham, North Carolina, United States of America

L. Barth Reller

Affiliations Division of Infectious Diseases, Department of Medicine, Duke University Medical Center, Durham, North Carolina, United States of America, Clinical Microbiology Laboratory, Duke University Medical Center, Durham, North Carolina, United States of America

Vance G. Fowler Jr.

Affiliations Duke Clinical Research Institute, Duke University Medical Center, Durham, North Carolina, United States of America, Division of Infectious Diseases, Department of Medicine, Duke University Medical Center, Durham, North Carolina, United States of America

Introduction

The clinical spectrum of Staphylococcus aureus is diverse, ranging from asymptomatic colonization to osteoarticular infections, endocarditis, and death [1], [2], [3], [4], [5], [6]. Although a growing body of evidence suggests that bacterial genetic characteristics may be associated with distinct clinical manifestations of S. aureus [7], [8], [9], [10], [11], [12], [13], the association between S. aureus genes and severity of illness is incompletely understood.

Using multi-locus sequence typing (MLST), we previously demonstrated a significant association between specific clonal complexes (CC) and severity of infection in S. aureus [9]. Among 371 clinically well-characterized S. aureus from diverse clinical settings, we showed that CC5 and CC30 were significantly associated with hematogenous complications of S. aureus such as endocarditis, septic arthritis, and vertebral osteomyelitis. Isolates within these CCs were also implicated using staphylococcal protein A (spa) [14] and staphylococcal cassette chromosome mec (SCCmec) typing [15]. However, the genetic basis for these associations is unknown.

The current investigation seeks to evaluate potential associations between the presence of specific genes and clinical outcome. To do this, we used a novel array comparative genome hybridization (aCGH) analysis method to identify individual genes in a large collection of clinically well-characterized S. aureus isolates from a wide range of infection severity [16]. We hypothesized that the association between CC5 and CC30 and hematogenous complications was due to the presence of specific bacterial genes contained within these clonal backgrounds, and tested this hypothesis for over 2,000 genes in the genomes of 239 clinical S. aureus isolates.

Methods

Subject characteristics

This study was approved by the Duke University Medical Center (DUMC) Institutional Review Board. Written informed consent was obtained from all participants involved in this study. At Duke Medical Center, hospitalized adult patients with at least one blood culture positive for S. aureus were prospectively identified between September 1994 and October 2003 [17]. Bacterial bloodstream isolates were cataloged, stored and corresponding clinical data were prospectively collected [17]. Twelve weeks after initial positive blood culture, clinical outcomes were prospectively assessed. Patients were excluded from the investigation for the following reasons: outpatient status, age <18 years, polymicrobial infection, neutropenia (absolute neutrophil count <1.0×109/liter), death before evaluation by the investigators, injection drug use, surgery within 30 days, or indwelling prostheses. These exclusion criteria eliminate important non-pathogen characteristics that can influence the severity of staphylococcal disease.

Next, we used strict definitions to identify three phenotypically distinct clinical groups which represent a progression from healthy individuals to those who are severely affected: 1) nasal carriage, 2) uncomplicated infection, and 3) bacteremia with hematogenous complications. Nasal carriage isolates were obtained from the nares of non-infected subjects in two distinct epidemiologic settings within the DUMC referral area: a) healthy undergraduate students with no current infection or healthcare contact (community carriage isolates) and hospitalized or hemodialysis-dependent subjects with no active infection (healthcare-associated carriage isolates).

Isolates in the uncomplicated infection group were obtained from patients with uncomplicated bacteremia or soft tissue infection. Uncomplicated bacteremia isolates were obtained from the bloodstream of DUMC patients meeting all of the following criteria (in addition to the criteria listed above): a) confirmed intravascular catheter-associated S. aureus bacteremia [17], b) receipt of ≤14 days of antibiotics, c) no evidence of metastatic infection by both clinical evaluation and echocardiogram, and d) alive without complications 12 weeks after the initial blood culture. In addition, to ensure that study findings were not attributable to the single center design, MSSA isolates (n = 48) were obtained from skin and soft tissue infections of anonymous geographically diverse patients participating in a Phase II clinical trial (Vicuron Pharmaceuticals). These patients had a culture-confirmed soft-tissue S. aureus infection with signs of systemic illness, but sterile blood cultures documented at the time of enrollment.

Isolates from the third group, bacteremia with hematogenous complications, were from the bloodstream of patients with definite native aortic or mitral valve IE [18], hematogenous bone and joint infection, or both. Patients with bone and joint infection had S. aureus isolated from a vertebral body, intervertebral disk, or joint space. All patients with bacteremia were followed up 12 weeks after the initial positive blood culture to confirm that infection groupings were accurate.

Microarray design

The S. aureus microarrays used for these experiments were obtained from the Pathogen Functional Genomics Resource Center (PFGRC) at The Institute for Genomic Research (TIGR, now the JCVI) (S. aureus version 3 microarrays). The arrays consist of 10,800 total spots, with 688 empty spots, 1,024 reserved spots containing Arabidopsis control oligonucleotides and 9,088 spots representing duplicate spotted oligos (n = 4544 70 base oligomers) of all known open reading frames in S. aureus at the time of array construction (N315 and Mu50 [18], MSSA476 and MRSA252 [19], COL [20] and MW2 [21]). The 688 empty spots were used as negative controls to measure the signal values that represent background in fluorescence.

Determination of microarray probe specificity

Preliminary examination of hybridization probe specificity suggested that several microarray probes had the potential to hybridize to multiple targets. To evaluate the hybridization specificity of the array probes, we used BLAST to determine the similarity between all individual probe sequences and the open reading frames in the S. aureus genomes included on the array. BLAST was completed using a word size of 30 and the alignment results filtered to require 90% sequence overlap. The e-value and bit score were recorded for each probe. Since some probes matched multiple genes within one genome, the total number of hits with bit scores above 100 was also noted. All genes that had a significant probe hit were aligned using CLUSTAL and the resulting set was then collected and Hidden Markov Models [22] were used to predict function with the Interpro scan program. Probes representing known pathogenicity islands, virulence factors and surface proteins were manually identified and where possible, ordered to identify their relative position in the S. aureus genome.

Identification of the S. aureus core genome and positive microarray controls

BLAST alignment of microarray probes (using criteria as described above) was used to identify which of the 4544 probes are hybridizing to genes shared among sequenced S. aureus genomes (COL, N315, MW2, Mu50, MRSA252 and MSSA476). Sequences with a bit score of greater than 100 in all sequenced genomes, excluding known virulence factors, surface proteins and genes in pathogenicity islands (PI) were designated as core (Table S1). The 1930 core probes represent 1660 genes. These probes were used as positive controls for aCGH analysis of the clinical strains. Of the remaining 2,614 probes, 600 were determined to not be variable in our sample based upon hybridization to all (or none) of the arrays (Table S1). Thus 2,014 variable probes remained, including virulence factors and other factors likely responsible for increased fitness or virulence. We used this definition of core and variable in case some of the known virulence factors were variable in our study population, even though they were not variable in the existing sequenced strains. We found this was the case and seven genes were included in our study in this way. These genes are serine proteases splC (SACOL1867), splD (SA1628) and splF (SA1627) and also cap5M (SACOL0148), cell wall hydrolase (SACOL1264), entertoxin homolog (SA1429) and hypothetical (SA0850).

Selection of strains

A total of 379 isolates (254 MSSA and 125 MRSA) were in our initial study [9] based on narrow clinical groupings described above. The CGH analysis was designed to include all of the MRSA isolates (n = 125) and a 50% random sample of the 244 MSSA isolates with non-missing spa and MLST types, stratifying on clinical severity. Four MRSA isolates and 4 MSSA were excluded during quality control checks. Thus, the final aCGH sample for this study included 239 isolates (121 MRSA; 118 MSSA).

Preparation of S. aureus genomic DNA

S. aureus isolates were grown overnight at 37°C in 1.5 ml Brain Heart Infusion (BHI) media. Genomic DNA was extracted from 500 µl of the cultures using the automated Qiagen M48 extractor and reagents (Qiagen, Valencia, CA), modified by the addition of lysostaphin (Sigma, St. Louis, MO) to the cell lysis buffer (final concentration of lysostaphin at 25 mg/ml).

Array Comparative Genome Hybridization (aCGH)

For hybridization, genomic DNA (500 ng) was amplified in Klenow reaction buffer containing random oligonucleotide primers (Invitrogen, Carlsbad, CA), high activity Klenow (New England Biolabs, Ipswich, MA) and dNTP/ddUTP. Amplified DNA was coupled to cy3 or cy5 and purified on Qiagen columns (Qiagen, Valencia, CA) prior to hybridization. The arrays were hybridized, washed and scanned using standard methods [23]. A study on a set of already sequenced strains demonstrated that single replicates where calls of presence/absence were made relative to negative controls showed greater sensitivity and specificity than a reference design (data not shown). Accordingly, we developed a hybridization design where each sample was hybridized once. The groups (control (C), uncomplicated (U), and complicated/sick (S) were each balanced for dye such that half of the samples in a particular group were hybridized with cy3 and half of the samples with cy5 and different clinical groupings represented on each slide. The control group was considered the reference and 36 control strains were repeated. The strains to be repeated were chosen at random from among all control strains. Quality control on these repeated hybridizations found excellent agreement. Images were quantified with Imagene v.7. Local background signals were subtracted from the signal intensity of the spot. Imagene's automated spot flagging feature was used to identify and remove poor quality spots (i.e. those contaminated by small dust particles). Each slide was checked to ensure that less than 5% of the slide was affected by imperfections and poor quality. Spots were determined to be present if they were greater than 90% of the negative controls. Validation of the aCGH method is described in Supporting Information File S1.

Associations

As the primary interest is identifying genes contributing to severe disease, the uncomplicated and control samples were combined into a single group, ‘not-complicated’, and compared to the group, ‘complicated infection’ for significance testing. The presence/absence of a particular probe was compared to this binary clinical outcome (not complicated vs. complicated) using a chi-squared test. P-values for this test were computed exactly as some comparisons had small cells [24]. A false discovery rate (FDR), of 0.20 was used to control for multiple testing [25], [26], [27]. This level was chosen to minimize the chance of failing to identify important loci. Further descriptive comparisons are reported to be significant at a nominal level of 0.05. All unadjusted p-values are reported in Table S2.

PCR verification of microarray data

The presence/absence of the 14 genes associated with severity of illness was verified by PCR amplification of genomic DNA from the corresponding strains. The primers, which flank the 70-mer oligonucleotide probes for these genes, are listed in Table S3. PCRs were said to agree with array results if simple agreement was greater than 0.60.

Results

We used aCGH to survey 239 clinical S. aureus isolates with the goal of identifying genes that may contribute to increased virulence. The genotypic characteristics of the study subjects with asymptomatic S. aureus nasal carriage, soft tissue infection, uncomplicated S. aureus bacteremia and complicated S. aureus bacteremia are listed in Table 1. Six CCs (CC 1, 5, 8, 15, 30, 45) out of 15 CCs total, were observed in 10 or more subjects in either the MRSA or MSSA subset in the current study. The S. aureus microarray used for this aCGH study allows for resolution of variations such as gene presence or absence as well as genomic rearrangements [12], but does not allow detection of individual DNA polymorphisms [11], [28]. Data from this study demonstrate: 1) genetic differences within the 239 isolates, including surface proteins and toxins, 2) genes unique to CC5 and CC30, 3) genes in the mec element associated with disease severity and 4) non-mec genes associated with disease severity (see Table S2).

[Figure omitted. See PDF.]

Table 1. Genotypic characteristics of 239 study subjects.

https://doi.org/10.1371/journal.pone.0018673.t001

Frequency of gene presence or absence in clonal complex

Comparative genomic analysis revealed that the oligonucleotides included on the S. aureus microarray contain 1930 probes to ORFs common in all sequenced genomes at the time, and 2,614 variable probes. The variable component represents strain specific genes including virulence factors and surface proteins frequently present among the sequenced genomes. Analysis of our aCGH data found that 600 of these probes were not variable in our sample of 239 strains, with 569 always present and 31 always absent. The frequency among the 5 major CCs in our study was calculated for each of the remaining 2,014 probes that were variable in our population.

Within CC5 and CC30 previously associated with clinical outcome [9], there were 1232 (61%) sequences that differed in frequency, with 731 (38%) being more frequent in CC30 and 461 (23%) more frequent in CC5 (Table S2). We also compared the presence/absence of each probe sequence in CC5 with other (non CC30) clonal complexes and the presence/absence of each sequence in CC30 with other (non CC5) clonal complexes. We found 626 (31%) genes more frequently present in CC30 and 249 (12%) genes more frequently present in CC5. There were 98 genes more present in both CC5 and CC30 compared to other clonal complexes.

Genes associated with the SCCmec-element (MRSA)

Our initial study demonstrated that MRSA status was associated with clinical outcome [9]. In order to understand the relationship between the 2,014 variable genes and the SCCmec element and MRSA status, we compared the presence/absence of genes in MRSA and MSSA strains. We found 140 genes significantly more frequent in MRSA and 1079 genes more frequent in MSSA (Table S2).

Genes associated with severe infection

From the 2014 variable sequences, the presence/absence of 226 sequences were statistically significantly associated with severity of infection (Figure 1 and Table S3). Within these 226 sequences, there were 41 genes known to be part of a SCCmec element. Among these 41 genes, 37 were present more frequently in the complicated group. Of the 185 genes not in the SCCmec element, 14 were present more frequently in the complicated group. For this study we further focused on these 14 non-SCCmec associated genes which were more common in the complicated group (Table 2). The relative proportion of these 14 genes in each of the major clonal complexes in our study is shown in Figure 2.

[Figure omitted. See PDF.]

Figure 1. Disposition of 4544 aCGH probes.

This figure illustrates how we arrived at our candidate genes of interest. Of the 2014 sequences and 239 isolates, a total of 226 sequences were statistically significant after multiple comparisons adjustment. Of these 226 sequences, 185 were not and 41 were in the SCCmec element. Most of the 185 non-SCCmec genes-171-were significantly less common in the complicated infection group. The remaining 14 genes were significantly more common in isolates from the complicated infection group. Among the 41 SCCmec element genes, 37 were significantly more common in the complicated group and 4 were significantly less common in the complicated group. Therefore, a total of 51 of the 2014 sequences, 14 non-mec associated and 37 mec associated, were identified as genes of interest.

https://doi.org/10.1371/journal.pone.0018673.g001

[Figure omitted. See PDF.]

Figure 2. Distribution of 14 genes associated with complicated infections across clonal complexes.

Proportion of isolates that have the 14 genes in each clonal complex shown on a color gradient with 1 being bright red and 0 being bright blue. All 14 genes are more abundant in CC5 and 30, the clonal complexes previously shown to exhibit a significant trend toward increasing levels of hematogenous complications.

https://doi.org/10.1371/journal.pone.0018673.g002

[Figure omitted. See PDF.]

Table 2. Putative virulence genes associated with complicated infections.

https://doi.org/10.1371/journal.pone.0018673.t002

Verification of microarray findings

The results of the aCGH analysis were verified by PCR amplification for 13 of the 14 loci described above. Amplification of PCR primers failed for one gene (SAR1522). For 12 out of the 13 genes (92.3%), the concordance between the array presence/absence data and the PCR validation was excellent. In comparison, the concordance of most large-scale array vs. PCR comparisons [29], [30] is ∼60%, clearly demonstrating the quality of our aCGH data. For genes that were PCR validated, the association with severity of infection is reported based upon the PCR results rather than the array results (Table 2).

Discussion

The ability of S. aureus to cause severe infections is dependent upon a large repertoire of virulence factors, often on mobile genetic elements or MGE (including pathogenicity islands and bacteriophage) that are transferred horizontally through the S. aureus community. The virulence capacity of individual isolates is largely determined by the function and abundance of these variable virulence factors. In our previous work [9], we demonstrated that S. aureus CC5 and CC30 are associated with the presence of hematogenous complications, including endocarditis, septic arthritis and vertebral osteomyelitis. We now undertake aCGH analysis to examine the association of 2014 variable S. aureus genes with severity of infection in a collection of 239 isolates with the goal of identifying genes that contribute to virulence. The 239 S. aureus isolates were selected from three clinical groups in our previous study [9] and included 121 MRSA and 118 MSSA isolates. Of the 226 genes that were identified as significantly associated with severity of infection, 51 were more frequent and 175 less frequent in the complicated infection group (Table S4).

We acknowledge that as a result of the microarray probe design being limited to the six genomes available when the microarray was constructed, our analysis is compromised and will not capture those genes present in many other S. aureus lineages and strains. As examples, genes encoding numerous heavy metal resistance genes on the SCCmec type III element as well as genes encoding the lukSF-PV Panton Valentine Leukotoxin and ACME element from Community Acquired MRSA (CA-MRSA) will not be identified in our analysis. However, a significant strength of our study is identification of the core S. aureus ORFs present in all tested clinical strains representing fifteen clonal complexes.

To minimize SCCmec as a confounding factor, the 226 variable genes were segregated into two groups, based on their presence in the SCCmec element (Figure 1). Of the 51 genes more common in complicated infection, 37 were carried on SCCmec and 14 were not carried on SCCmec. None of these 14 genes encode known virulence factors previously associated with increased virulence and disease severity. Instead, 11 genes include proteins of unknown function, regulatory proteins and autolysins carried on S. aureus bacteriophage. These findings support previous work demonstrating the contribution of S. aureus bacteriophage to pathogenesis of infections [31], [32], [33], [34] through genome variation, mobilization of pathogenicity islands and transmission of phage encoded virulence factors.

Yet to be resolved is the contribution of prophage to potentiate S. aureus virulence by influencing its physiology towards growth as a biofilm, a virulence factor that plays a role in chronic infections, such as native valve endocarditis and osteomyelitis. Evidence of bacteriophage release from S. aureus biofilms has been reported [35], leading to lysis of cells which in turn promotes the persistence and survival of remaining cells and release of extracellular DNA (eDNA) required for extracellular polymeric substance (EPS) development and maturation of biofilms. A second mechanism for cell lysis and release of eDNA required for biofilm development involves autolysis driven by autolysins and regulated by the cid and lrg loci [36]. These autolysins are capable of hydrolyzing the amide bond between N-acetylmuramic acid and L-alanine in bacterial cell-wall peptidoglycan resulting in cell lysis and release of eDNA. Interestingly, two of the bacteriophage genes more frequently present in complicated infections are autolysins carried on bacteriophages Sa2MW (MW1380) and L54a (SACOL0389) (Table 2). Although, no role has yet to be assigned to these two putative autolysin genes, their presence in cases of complicated infection suggests that they may influence biofilm maturation and disease persistence. An alternative function for the S. aureus phage encoded autolysins is suggested by recent work on bacteriophage-encoded autolysins in Streptococcus mitis, which demonstrated their role in binding to human platelets through interaction with fibrinogen [37].

Of the remaining 9 genes carried on prophage, three are DNA binding proteins with potential regulatory functions (SAS0928, SAR1502 and SAR1555) and six (SACOL0329, SAR1546, SAR1522, MW1385, SACOL0338 and SA0397) are genes of unknown function. The function of three remaining genes (VRA0002, SA1483 and VRA0005) relative to virulence is also not clear. Association of these genes with increased disease severity may be a result of their linkage to unidentified genes not on the microarray and associated with virulence. This study did not start with candidate genes, and as such uncovered a list of genes for which little is known about the function. Detailed functional analysis of these and other S. aureus bacteriophage genes of unknown function will likely lead to a greater understanding of their role in virulence and enabling mutation and adaptation of S. aureus to new environments.

The 37 genes that were more frequently present in complicated infection and carried by SCCmec span the entire length of the type II SCCmec element, which is found in the majority (97/120; 81%) of MRSA isolates included in this study. While there is no direct evidence that SCCmec contributes to S. aureus virulence, recent work by Queck et. al. [38] has shown the type II SCCmec elements encode a new phenol soluble modulin (PSM) that alters the capacity of S. aureus to cause disease. Although the arrays used in this study do not contain a probe for this PSM, two genes flanking the PSM encoding region in SCCmec type II, mecI (N315-SA0039) and xylR (N315-SA0041) are associated with severity of infection, suggesting that the new PSM may also be a significant component in our study. However, sequence analysis of this region showed that the PSM is present in all members of a subset of 69 isolates from the study group representing all 3 clinical groups and student controls. Outside of the SCCmec element, there were 140 genes associated with SCCmec and severity of infection (Table S4). This is intriguing as it suggests the possibility of co-evolution of genes in the S. aureus genome with SCCmec and the possible preferential insertion of SCCmec into specific strains.

An unexpected finding in our study was the identification of 171 genes more frequent in strains isolated in subjects with uncomplicated infection. These genes (included in Table S4) include 152 genes of unknown function, some of which may be involved in attenuation of virulence. Approximately 38% of S. aureus genes have no assigned function and an additional 10% have a putative function based on functional domain matches [39]. Their roles as genes that contribute to S. aureus specific functions such as pathogenicity or host niche adaptation, remains an unexplored area of this significant bacterial pathogen.

In summary, we have demonstrated a statistical association between the presence/absence of individual genes and clinical outcome. Our data strongly supports the concept that differences in disease potential are rooted in the S. aureus genotype. However, the use of aCGH in our current study limits the results to identification of large-scale genomic changes, such as gain or loss of genes, responsible for phenotypic changes. Additional genotypic factors, including genomic polymorphisms, novel genes and differential gene expression, associated with complicated infections, will be determined by resequencing of representative members of the SABG collection and global transcriptional analysis from in vitro and in vivo environments. Correlation of these genomic and transcriptional differences with clinical outcome will contribute significantly to our understanding of S. aureus virulence.

Supporting Information

[Figure omitted. See PDF.]

File S1.

Supplementary Text and Description of Tables S1, S2 and S3.

https://doi.org/10.1371/journal.pone.0018673.s001

(DOC)

Table S1.

Ubiquitous genes described in Figure 1 (1930 core family genes and 600 not variable in our cohort).

https://doi.org/10.1371/journal.pone.0018673.s002

(XLS)

Table S2.

Overall association of S. aureus genes with clinical outcome.

https://doi.org/10.1371/journal.pone.0018673.s003

(XLS)

Table S3.

Primers used in this study.

https://doi.org/10.1371/journal.pone.0018673.s004

(XLS)

Table S4.

Genes statistically significantly associated with severity of infection.

https://doi.org/10.1371/journal.pone.0018673.s005

(XLS)

Author Contributions

Conceived and designed the experiments: SRG LMM CLN BR LBR VGF. Performed the experiments: SRG BR TR. Analyzed the data: LMM CLN SRG VGF. Contributed reagents/materials/analysis tools: SRG LMM CLN VGF LBR. Wrote the paper: SRG LMM CLN VGF.

Citation: Gill SR, McIntyre LM, Nelson CL, Remortel B, Rude T, Reller LB, et al. (2011) Potential Associations between Severity of Infection and the Presence of Virulence-Associated Genes in Clinical Strains of Staphylococcus aureus. PLoS ONE6(4): e18673. https://doi.org/10.1371/journal.pone.0018673

References

1. Cabell CH, Jollis JG, Peterson GE, Corey GR, Anderson DJ, et al. (2002) Changing patient characteristics and the effect on mortality in endocarditis. Arch Intern Med 162: 90–94.CH CabellJG JollisGE PetersonGR CoreyDJ Anderson2002Changing patient characteristics and the effect on mortality in endocarditis.Arch Intern Med1629094

2. Fernandez-Guerrero ML, Verdejo C, Azofra J, de Gorgolas M (1995) Hospital-acquired infectious endocarditis not associated with cardiac surgery: an emerging problem. Clin Infect Dis 20: 16–23.ML Fernandez-GuerreroC. VerdejoJ. AzofraM. de Gorgolas1995Hospital-acquired infectious endocarditis not associated with cardiac surgery: an emerging problem.Clin Infect Dis201623

3. Ing M, Baddour L, Bayer AS (1997) Bacteremia and infective endocarditis: pathogenesis, diagnosis, and complications. . In: Crossley K, Archer GL, editors. The Staphylococci in human disease. New York: Churchill Livingstone. pp. 331–354.M. IngL. BaddourAS Bayer1997Bacteremia and infective endocarditis: pathogenesis, diagnosis, and complications. .K. CrossleyGL ArcherThe Staphylococci in human diseaseNew YorkChurchill Livingstone331354

4. Jensen AG, Espersen F, Skinhoj P, Rosdahl VT, Frimodt-Moller N (1997) Increasing frequency of vertebral osteomyelitis following Staphylococcus aureus bacteraemia in Denmark 1980–1990. J Infect 34: 113–118.AG JensenF. EspersenP. SkinhojVT RosdahlN. Frimodt-Moller1997Increasing frequency of vertebral osteomyelitis following Staphylococcus aureus bacteraemia in Denmark 1980–1990.J Infect34113118

5. Jensen AG, Wachmann CH, Poulsen KB, Espersen F, Scheibel J, et al. (1999) Risk factors for hospital-acquired Staphylococcus aureus bacteremia. Arch Intern Med 159: 1437–1444.AG JensenCH WachmannKB PoulsenF. EspersenJ. Scheibel1999Risk factors for hospital-acquired Staphylococcus aureus bacteremia.Arch Intern Med15914371444

6. Sanabria TJ, Alpert JS, Goldberg R, Pape LA, Cheeseman SH (1990) Increasing frequency of staphylococcal infective endocarditis. Experience at a university hospital, 1981 through 1988. Arch Intern Med 150: 1305–1309.TJ SanabriaJS AlpertR. GoldbergLA PapeSH Cheeseman1990Increasing frequency of staphylococcal infective endocarditis. Experience at a university hospital, 1981 through 1988.Arch Intern Med15013051309

7. Bhattacharya D, Carleton H, Tsai CJ, Baron EJ, Perdreau-Remington F (2007) Differences in clinical and molecular characteristics of skin and soft tissue methicillin-resistant Staphylococcus aureus isolates between two hospitals in Northern California. J Clin Microbiol 45: 1798–1803.D. BhattacharyaH. CarletonCJ TsaiEJ BaronF. Perdreau-Remington2007Differences in clinical and molecular characteristics of skin and soft tissue methicillin-resistant Staphylococcus aureus isolates between two hospitals in Northern California.J Clin Microbiol4517981803

8. Campbell SJ, Deshmukh HS, Nelson CL, Bae IG, Stryjewski ME, et al. (2008) Genotypic characteristics of Staphylococcus aureus isolates from a multinational trial of complicated skin and skin structure infections. J Clin Microbiol 46: 678–684.SJ CampbellHS DeshmukhCL NelsonIG BaeME Stryjewski2008Genotypic characteristics of Staphylococcus aureus isolates from a multinational trial of complicated skin and skin structure infections.J Clin Microbiol46678684

9. Fowler VG Jr, Nelson CL, McIntyre LM, Kreiswirth BN, Monk A, et al. (2007) Potential associations between hematogenous complications and bacterial genotype in Staphylococcus aureus infection. J Infect Dis 196: 738–747.VG Fowler JrCL NelsonLM McIntyreBN KreiswirthA. Monk2007Potential associations between hematogenous complications and bacterial genotype in Staphylococcus aureus infection.J Infect Dis196738747

10. Herron-Olson L, Fitzgerald JR, Musser JM, Kapur V (2007) Molecular Correlates of Host Specialization in Staphylococcus aureus. PLoS ONE 2: e1120.L. Herron-OlsonJR FitzgeraldJM MusserV. Kapur2007Molecular Correlates of Host Specialization in Staphylococcus aureus.PLoS ONE2e1120

11. Highlander SK, Hulten KG, Qin X, Jiang H, Yerrapragada S, et al. (2007) Subtle genetic changes enhance virulence of methicillin resistant and sensitive Staphylococcus aureus. BMC Microbiol 7: 99.SK HighlanderKG HultenX. QinH. JiangS. Yerrapragada2007Subtle genetic changes enhance virulence of methicillin resistant and sensitive Staphylococcus aureus.BMC Microbiol799

12. Lindsay JA, Moore CE, Day NP, Peacock SJ, Witney AA, et al. (2006) Microarrays reveal that each of the ten dominant lineages of Staphylococcus aureus has a unique combination of surface-associated and regulatory genes. J Bacteriol 188: 669–676.JA LindsayCE MooreNP DaySJ PeacockAA Witney2006Microarrays reveal that each of the ten dominant lineages of Staphylococcus aureus has a unique combination of surface-associated and regulatory genes.J Bacteriol188669676

13. van Belkum A, Melles DC, Snijders SV, van Leeuwen WB, Wertheim HF, et al. (2006) Clonal distribution and differential occurrence of the enterotoxin gene cluster, egc, in carriage- versus bacteremia-associated isolates of Staphylococcus aureus. J Clin Microbiol 44: 1555–1557.A. van BelkumDC MellesSV SnijdersWB van LeeuwenHF Wertheim2006Clonal distribution and differential occurrence of the enterotoxin gene cluster, egc, in carriage- versus bacteremia-associated isolates of Staphylococcus aureus.J Clin Microbiol4415551557

14. Mathema B, Mediavilla J, Kreiswirth BN (2008) Sequence analysis of the variable number tandem repeat in Staphylococcus aureus protein A gene: spa typing. Methods Mol Biol 431: 285–305.B. MathemaJ. MediavillaBN Kreiswirth2008Sequence analysis of the variable number tandem repeat in Staphylococcus aureus protein A gene: spa typing.Methods Mol Biol431285305

15. Chen L, Mediavilla JR, Oliveira DC, Willey BM, de Lencastre H, et al. (2009) Multiplex real-time PCR for rapid Staphylococcal cassette chromosome mec typing. J Clin Microbiol 47: 3692–3706.L. ChenJR MediavillaDC OliveiraBM WilleyH. de Lencastre2009Multiplex real-time PCR for rapid Staphylococcal cassette chromosome mec typing.J Clin Microbiol4736923706

16. Fowler VG Jr, Olsen MK, Corey GR, Woods CW, Cabell CH, et al. (2003) Clinical identifiers of complicated Staphylococcus aureus bacteremia. Arch Intern Med 163: 2066–2072.VG Fowler JrMK OlsenGR CoreyCW WoodsCH Cabell2003Clinical identifiers of complicated Staphylococcus aureus bacteremia.Arch Intern Med16320662072

17. Benjamin DK Jr, Miro JM, Hoen B, Steinbach WJ, Fowler VG Jr, et al. (2004) Candida endocarditis: contemporary cases from the International Collaboration of Infectious Endocarditis Merged Database (ICE-mD). Scand J Infect Dis 36: 453–455.DK Benjamin JrJM MiroB. HoenWJ SteinbachVG Fowler Jr2004Candida endocarditis: contemporary cases from the International Collaboration of Infectious Endocarditis Merged Database (ICE-mD).Scand J Infect Dis36453455

18. Kuroda M, Ohta T, Uchiyama I, Baba T, Yuzawa H, et al. (2001) Whole genome sequencing of meticillin-resistant Staphylococcus aureus. Lancet 357: 1225–1240.M. KurodaT. OhtaI. UchiyamaT. BabaH. Yuzawa2001Whole genome sequencing of meticillin-resistant Staphylococcus aureus.Lancet35712251240

19. Holden MT, Feil EJ, Lindsay JA, Peacock SJ, Day NP, et al. (2004) Complete genomes of two clinical Staphylococcus aureus strains: evidence for the rapid evolution of virulence and drug resistance. Proc Natl Acad Sci U S A 101: 9786–9791.MT HoldenEJ FeilJA LindsaySJ PeacockNP Day2004Complete genomes of two clinical Staphylococcus aureus strains: evidence for the rapid evolution of virulence and drug resistance.Proc Natl Acad Sci U S A10197869791

20. Gill SR, Fouts DE, Archer GL, Mongodin EF, Deboy RT, et al. (2005) Insights on evolution of virulence and resistance from the complete genome analysis of an early methicillin-resistant Staphylococcus aureus strain and a biofilm-producing methicillin-resistant Staphylococcus epidermidis strain. J Bacteriol 187: 2426–2438.SR GillDE FoutsGL ArcherEF MongodinRT Deboy2005Insights on evolution of virulence and resistance from the complete genome analysis of an early methicillin-resistant Staphylococcus aureus strain and a biofilm-producing methicillin-resistant Staphylococcus epidermidis strain.J Bacteriol18724262438

21. Baba T, Takeuchi F, Kuroda M, Yuzawa H, Aoki K, et al. (2002) Genome and virulence determinants of high virulence community-acquired MRSA. Lancet 359: 1819–1827.T. BabaF. TakeuchiM. KurodaH. YuzawaK. Aoki2002Genome and virulence determinants of high virulence community-acquired MRSA.Lancet35918191827

22. Eddy SR (2004) What is a hidden Markov model? Nat Biotechnol 22: 1315–1316.SR Eddy2004What is a hidden Markov model?Nat Biotechnol2213151316

23. Vickerman MM, Iobst S, Jesionowski AM, Gill SR (2007) Genome-wide transcriptional changes in Streptococcus gordonii in response to competence signaling peptide. J Bacteriol 189: 7799–7807.MM VickermanS. IobstAM JesionowskiSR Gill2007Genome-wide transcriptional changes in Streptococcus gordonii in response to competence signaling peptide.J Bacteriol18977997807

24. Edgington ES (1995) Randomization Tests. New York: Marcel Dekker. ES Edgington1995Randomization TestsNew YorkMarcel Dekker

25. Benjamini Y, Yekutieli D (2005) Quantitative trait Loci analysis using the false discovery rate. Genetics 171: 783–790.Y. BenjaminiD. Yekutieli2005Quantitative trait Loci analysis using the false discovery rate.Genetics171783790

26. Storey JD, Tibshirani R (2003) Statistical significance for genomewide studies. Proc Natl Acad Sci U S A 100: 9440–9445.JD StoreyR. Tibshirani2003Statistical significance for genomewide studies.Proc Natl Acad Sci U S A10094409445

27. Verhoeven KJ, Simonsen KL (2005) Genomic haplotype blocks may not accurately reflect spatial variation in historic recombination intensity. Mol Biol Evol 22: 735–740.KJ VerhoevenKL Simonsen2005Genomic haplotype blocks may not accurately reflect spatial variation in historic recombination intensity.Mol Biol Evol22735740

28. Kennedy AD, Otto M, Braughton KR, Whitney AR, Chen L, et al. (2008) Epidemic community-associated methicillin-resistant Staphylococcus aureus: recent clonal expansion and diversification. Proc Natl Acad Sci U S A 105: 1327–1332.AD KennedyM. OttoKR BraughtonAR WhitneyL. Chen2008Epidemic community-associated methicillin-resistant Staphylococcus aureus: recent clonal expansion and diversification.Proc Natl Acad Sci U S A10513271332

29. Dallas PB, Gottardo NG, Firth MJ, Beesley AH, Hoffmann K, et al. (2005) Gene expression levels assessed by oligonucleotide microarray analysis and quantitative real-time RT-PCR - how well do they correlate? BMC Genomics 6: 59.PB DallasNG GottardoMJ FirthAH BeesleyK. Hoffmann2005Gene expression levels assessed by oligonucleotide microarray analysis and quantitative real-time RT-PCR - how well do they correlate?BMC Genomics659

30. Morey JS, Ryan JC, Van Dolah FM (2006) Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR. Biol Proced Online 8: 175–193.JS MoreyJC RyanFM Van Dolah2006Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR.Biol Proced Online8175193

31. Bae T, Schneewind O (2006) Allelic replacement in Staphylococcus aureus with inducible counter-selection. Plasmid 55: 58–63.T. BaeO. Schneewind2006Allelic replacement in Staphylococcus aureus with inducible counter-selection.Plasmid555863

32. Goerke C, Wirtz C, Fluckiger U, Wolz C (2006) Extensive phage dynamics in Staphylococcus aureus contributes to adaptation to the human host during infection. Mol Microbiol 61: 1673–1685.C. GoerkeC. WirtzU. FluckigerC. Wolz2006Extensive phage dynamics in Staphylococcus aureus contributes to adaptation to the human host during infection.Mol Microbiol6116731685

33. Rolain JM, Francois P, Hernandez D, Bittar F, Richet H, et al. (2009) Genomic analysis of an emerging multiresistant Staphylococcus aureus strain rapidly spreading in cystic fibrosis patients revealed the presence of an antibiotic inducible bacteriophage. Biol Direct 4: 1.JM RolainP. FrancoisD. HernandezF. BittarH. Richet2009Genomic analysis of an emerging multiresistant Staphylococcus aureus strain rapidly spreading in cystic fibrosis patients revealed the presence of an antibiotic inducible bacteriophage.Biol Direct41

34. Wirtz C, Witte W, Wolz C, Goerke C (2009) Transcription of the phage-encoded Panton-Valentine leukocidin of Staphylococcus aureus is dependent on the phage life-cycle and on the host background. Microbiology 155: 3491–3499.C. WirtzW. WitteC. WolzC. Goerke2009Transcription of the phage-encoded Panton-Valentine leukocidin of Staphylococcus aureus is dependent on the phage life-cycle and on the host background.Microbiology15534913499

35. Resch A, Fehrenbacher B, Eisele K, Schaller M, Gotz F (2005) Phage release from biofilm and planktonic Staphylococcus aureus cells. FEMS Microbiol Lett 252: 89–96.A. ReschB. FehrenbacherK. EiseleM. SchallerF. Gotz2005Phage release from biofilm and planktonic Staphylococcus aureus cells.FEMS Microbiol Lett2528996

36. Rice KC, Bayles KW (2008) Molecular control of bacterial death and lysis. Microbiol Mol Biol Rev 72: 85–109.KC RiceKW Bayles2008Molecular control of bacterial death and lysis.Microbiol Mol Biol Rev7285109

37. Seo HS, Xiong YQ, Mitchell J, Seepersaud R, Bayer AS, et al. (2010) Bacteriophage lysin mediates the binding of streptococcus mitis to human platelets through interaction with fibrinogen. PLoS Pathog 6: e1001047.HS SeoYQ XiongJ. MitchellR. SeepersaudAS Bayer2010Bacteriophage lysin mediates the binding of streptococcus mitis to human platelets through interaction with fibrinogen.PLoS Pathog6e1001047

38. Queck SY, Khan BA, Wang R, Bach TH, Kretschmer D, et al. (2009) Mobile genetic element-encoded cytolysin connects virulence to methicillin resistance in MRSA. PLoS Pathog 5: e1000533.SY QueckBA KhanR. WangTH BachD. Kretschmer2009Mobile genetic element-encoded cytolysin connects virulence to methicillin resistance in MRSA.PLoS Pathog5e1000533

39. Holden M, Lindsay J (2008) Whole genomes: sequence, microarray and systems biology. In: Lindsay J, editor. Staphylococcus molecular genetics. Norfolk, UK: Caister Academic Press. pp. 1–28.M. HoldenJ. Lindsay2008Whole genomes: sequence, microarray and systems biology.J. LindsayStaphylococcus molecular geneticsNorfolk, UKCaister Academic Press128