ABSTRACT

Most cases of severe Staphylococcus aureus disease cannot be explained by the action of a single virulence determinant, and it is likely that a number of factors act in combination during the infective process. This study examined the relationship between disease in humans and a large number of putative virulence determinants, both individually and in combination. S. aureus isolates (n = 334) from healthy blood donors and from patients with invasive disease were compared for variation in the presence of 33 putative virulence determinants. After adjusting for the effect of clonality, seven determinants (fnbA, cna, sdrE, sej, eta, hlg, and ica) were significantly more common in invasive isolates. All seven factors contributed independently to virulence. No single factor predominated as the major predictor of virulence, their effects appearing to be cumulative. No combinations of the seven genes were either more or less likely to cause disease than others with the same number of virulence-associated genes. There was evidence of considerable horizontal transfer of genes on a background of clonality. Our findings also suggested that allelic variants of a polymorphic locus can make different contributions to the disease process, further study of which is likely to expand our understanding of staphylococcal disease pathogenesis.

Staphylococcus aureus is a major cause of severe community-acquired and nosocomial sepsis (12, 36). The relative importance of host factors versus bacterial virulence determinants in disease pathogenesis is unknown, but it is widely held that bacterial factors including toxins, cell wall-associated adhesins, and secreted exoproteins are involved in the process (34). There is clear evidence for this in the case of toxin-mediated S. aureus diseases such as toxic shock syndrome caused by toxic shock syndrome toxin 1, scalded skin syndrome resulting from the action of epidermolytic toxins, and food poisoning following the ingestion of preformed enterotoxins (6). However, the vast majority of cases of severe S. aureus disease cannot be explained by the action of a single virulence determinant and it is likely that a number of factors act in combination during the infective process. Support for this comes from animal models where outcome has been compared following infection with isogenic mutants deficient in a single putative virulence determinant, in which infection was attenuated but not prevented (17, 27). In addition, clinical isolates associated with human infection may be naturally deficient in a range of putative virulence determinants. It seems reasonable to postulate, therefore, that S. aureus strains associated with human infection have variable combinations of pathogenic determinants and that either the presence or the expression of given combinations varies depending on the type of infection and genetic susceptibility of the affected host.

The aim of this study was to examine a large number of putative virulence determinants in a natural population of S. aureus isolates and explore their relationship, both individually and in combination, to invasive S. aureus disease in humans. A panel of isolates was assembled during a prospective case control study in which strains were collected from healthy blood donors and from patients with invasive disease, the quality of clinical information being ensured by the prospective application of defined case definitions. The presence of virulence determinants was compared between the isolate groups associated with carriage and disease, analysis of which included adjustment for the potentially confounding effects of clonality and the hitchhiker effect. This incorporated a framework provided by multilocus sequence typing (MLST), a technique which also gave us insights into the horizontal gene transfer of virulence determinants. Finally, we explored the possibility that allelic variants at a given polymorphic locus (namely, sdrE/bbp) do not make the same contribution to the pathogenic process.

MATERIALS AND METHODS

Bacterial isolates.We examined 33 bacterial determinants in 155 isolates recovered from patients with invasive S. aureus disease (94 hospital acquired and 61 community acquired) and 179 isolates recovered from healthy individuals. These isolates were collected within Oxfordshire, United Kingdom, between 1997 and 1998 by using a prospective case control design. Cases of invasive infection were identified prospectively through the microbiology laboratory serving the main hospitals in the Oxford area. Clinical details were recorded, and the bacterial isolate was stored. Community-acquired disease was defined as admission to a hospital with an illness consistent with invasive S. aureus disease, with isolation of S. aureus from a normally sterile site within 24 h of admission. Hospital-acquired disease was defined as an illness consistent with invasive S. aureus disease and isolation of the organism from a normally sterile site 48 h or more after admission for another condition. Control bacteria were obtained from nasal swabs of blood donors drawn from the same population catchment area as Oxford Radcliffe Hospitals.

Evaluation of bacterial determinants.Twenty-nine determinants were examined for the presence of the gene by PCR, and four determinants were evaluated by phenotypic tests. The bacterial determinants examined incorporated the major putative virulence factors reported in the published literature at the design stages of this study (Table 1).

PCR analysis.Genomic DNA from S. aureus was extracted by using the Wizard Genomic DNA purification kit (Promega), with the modification that 30 μg of lysostaphin/ml (Ambi) was added at the cell lysis step. sea-i and tst were amplified by multiplex PCR with the primers and conditions described elsewhere (26); sej, eta, and etb were amplified with the primers and conditions provided by G. O'Neill (personal communication). Primers and conditions used to amplify efb, pvl, and hlg were as described in references 5 (efb) and 23 (pvl and hlg). Primer design for the remaining determinants was based on gene sequences available from GenBank (Table 2). This was straightforward with the exception of map/eap, the agr subgroup, and sdrE and bbp. There were three GenBank sequences available for map (AJ243790
, clinical isolate 7; AJ245439
, strain Wood 46; AJ223806
, strain Newman) and one for a gene termed map/eap (AJ290973
, strain Newman). These sequences were aligned, and common forward and reverse primers were designed. To distinguish between agr groups I, II, III, and IV, use was made of sequence differences in or around agrD, the region encoding the autoinducing peptide. The forward primers for agr subgroups I and IV started 9 bases upstream of the region encoding the peptide; primers for agr subgroups II and III started at the ninth and first bases, respectively, of the region encoding the peptide. The reverse primers for all agr subgroups were within agrC.

Primer design for amplification of sdrE and bbp was complicated by the fact that these independently described genes are known to share significant homology (38). We considered whether these genes were alleles of the same locus. The sequences for sdrE (GenBank accession number AJ005647
) and bbp (accession number Y18653
) were blasted against the available genome sequence for two isolates undergoing sequencing at the Sanger Centre, Cambridge, United Kingdom (http://www.sanger.ac.uk/Projects/S_aureus/blast_server.shtml
). These isolates were provided to the Sanger Centre by us and are called isolates 252 (methicillin-resistant isolate) and 476 (methicillin-sensitive isolate) in this paper. Blast analysis demonstrated that sdrE and bbp colocalized to identical positions in the respective isolates. Isolate 252 shared 94% homology at the amino acid level with bbp and 79% homology with sdrE. In contrast, isolate 476 shared 80% homology with bbp and 93% homology with sdrE. The sequences for sdrE and bbp were aligned using Gene Jockey II (Biosoft, Cambridge, United Kingdom), and primers were designed to distinguish between what we assume to be two alleles of a single gene by using a common forward primer but unique reverse primer. These primers were both piloted by using 50 randomly selected isolates. A clear distinction was seen between the primer pairs with the amplification of a single band of the appropriate size for SdrE (766 bp), bbp (1,054 bp), or neither; these primers were used to evaluate the remaining isolates.

PCR amplifications were performed in a PTC-200 DNA engine (MJ Research, Waltham, Mass.) with Taq polymerase (Bioline). The final concentrations of the PCR mixtures were 1× reaction buffer, variable MgCl2 concentrations depending on optimization results (Table 2), 100 pmol of forward and reverse primers, 1 μl of 1:5 dilution template DNA, 200 μM deoxynucleoside triphosphate mix, and 2.5 U of Taq polymerase. The PCR primers and cycling conditions not referred to elsewhere are shown in Table 2. Aliquots of the reaction mixtures were analyzed by 1% agarose gel electrophoresis. A positive control (Table 2) and a negative control (reaction mixture minus DNA) were included in each PCR run.

Phenotypic tests.Production of coagulase was evaluated using standard methodology (8). Production of hemolysins was determined by streaking isolates onto a range of blood agar plates and incubating them in air at 37°C for 18 h. Isolates were considered positive for alpha-toxin if hemolytic on 5% rabbit blood agar, positive for beta-toxin if hemolytic on 5% sheep blood agar following incubation at 37°C for 18 h and then 4°C for 30 min (hot-cold lysis), and positive for delta-toxin if synergy was observed between the test isolate and a beta-toxin producer on 3% washed sheep erythrocytes. Positive controls were NCTC 5655 (alpha-toxin positive), NCTC 7428 (beta-toxin positive), and NCTC 9715 (delta-toxin positive).

Analysis.The isolates used in this study were chosen for two reasons. First, the case control study design helps to avoid the potential confounding effect from comparing isolates from different geographical locations or time periods. Second, having defined the population structure of these isolates, we were able to adjust for the effects of clonality on associations between bacterial determinants and invasive disease. For example, a given determinant could become overrepresented in the invasive group if it is linked to a true virulence gene within a given lineage. We recognize that analysis linked to interpretation of MLST results becomes model dependent. In view of this, in the primary analysis, isolates were evaluated according to origin (nasal carriage isolate or invasive isolate) and further subdivided depending on whether the disease was hospital or community acquired. The second round of analysis used isolates grouped by lineage as defined by MLST results and analysis repeated for associations within each lineage. A lineage (or clonal complex) was defined for the purposes of this analysis as a group of organisms that had sequences at five or more of the seven MLST loci examined that were identical to those from at least one other organism in the group. Further detail of the genetic structure of S. aureus can be found on the MLST website (http://www.mlst.net
).

Contingency tables were used to compare proportions between groups. To adjust for the hitchhiker effect, analyses were stratified by MLST-defined lineage by either the Mantel-Haenszel method or conditional logistic regression. In this way, the presence of particular determinant(s) was only compared between disease and carriage isolates within lineages, and these separate effects summated to give an overall odds ratio. Strains that were unrelated to any other at five of seven loci were arbitrarily grouped as a lineage for the purpose of this stratification. No corrections were made for multiple comparisons, though in an attempt to recognize this, statistical significance was set at P = 0.01 for the purpose of interpretation and discussion. For the same reason, conservative 99% confidence intervals are used throughout.

RESULTS

Univariate analysis.Eight of the 33 determinants (24%) were significantly more common in invasive isolates than in carriage isolates (Table 3). The eight determinants were as follows: three genes encoding cell wall-associated adhesins (fnbA, encoding fibronectin binding protein A; cna, encoding collagen binding protein; and sdrE, encoding a protein of unknown function), four genes encoding toxins (sea and sej, encoding staphylococcal enterotoxins A and J, respectively; eta, encoding exfoliative toxin A; and hlg, encoding gamma-toxin), and icaA (used here as a marker for the entire ica operon), which is involved in biofilm production.

Presence of putative virulence determinants in carriage and disease isolates as determined by univariate analysis and analysis stratified by clonal complex

Adjusting for the effect of clonality.That such a large proportion of the genes examined were associated with disease (albeit all positively) raises the question of whether the result can be explained by linkage disequilibrium between these genes and one or more true virulence determinants elsewhere in the genome (the hitchhiker effect). In an attempt to adjust for this, we repeated the analysis, this time stratifying for clonality. The clonal structure of this collection of isolates has been defined by MLST (http://www.mlst.net
); using this information and the Mantel-Haenszel method, we compared the prevalence of each gene in the disease and carriage groups within but not between clonal complexes (Table 3 and Fig. 1). Following this stratification, all eight genes remained positively associated with disease. seg was associated with disease on stratified analysis alone. We also repeated the analysis, omitting epidemic methicillin-resistant S. aureus clone 16 (EMRSA-16) (sequence type ST36), by far the single largest clone, containing 22 disease strains from a single group of hospitals but no carriage strains, and found that without this clone sea was not associated with disease. In view of this, sea was not included as a virulence-associated gene in the remaining analysis. The other results were not affected, leaving seven putative virulence determinants positively associated with disease across multiple S. aureus lineages.

Presence or absence of putative virulence determinants. Odds ratios for disease, adjusted for the effects of clonality by Mantel-Haenszel stratification of MLST-defined lineage. Error bars denote 99% confidence intervals (CIs), and an asterisk indicates a significant association with disease (P < 0.01). Genes that were either ubiquitous (such as clfA and clfB) or very rare or absent (such as see and etb) are not shown.

Disease origin: community or hospital.The relative importance of the 33 bacterial factors in hospital-acquired versus community-acquired isolates was assessed by a comparison of carriage isolates and community-acquired disease alone. The presence of genes encoding six of the seven virulence-associated factors from the primary analysis (fnbA, cna, sdrE, sej, eta, and ica) remained associated with disease once hospital-acquired strains were removed. Although hlg was numerically more common in the community-acquired disease group than in the hospital-acquired group (98% versus 95%), comparison with the carriage isolates (88%) was no longer statistically significant.

Multivariate analysis. (i) Linkage disequilibrium between virulence-associated genes.We examined the possibility that identification of one or more of these seven virulence-associated genes was actually the result of linkage disequilibrium within this group. Logistic regression modeling of the relationship of the seven factors to both disease and each other suggested that four of the factors (sej, fnbA, cna, and sdrE) were both independent of each other and independently associated with disease. Two of the remaining three factors, ica and hlg, were in linkage disequilibrium with fnbA, though each contributed significantly to virulence after taking this association into account. eta was in linkage disequilibrium with sej but also contributed independently to virulence. We thus concluded that all seven factors were associated with virulence.

(ii) Role of individual genes in disease.Although upon univariate analysis the odds ratio for disease for fnbA was twice as high as that for any of the other virulence-associated genes, in the full multivariate model, the odds ratios for disease for each gene fell between 2.2 (eta) and 3.3 (fnbA). This suggests that each of the seven genes contributed approximately equally to virulence and that no single bacterial factor predominated as the major predictor of virulence.

(iii) Cumulative effect of genes.There was a linear trend between the number of these seven virulence determinants carried by particular isolates and the proportion of isolates causing disease (P < 0.0001) (Fig. 2A). This suggested that the effects of the factors on virulence were cumulative. The association between the number of determinants and the odds of disease was log linear, consistent with the logistic regression model in which there is a cumulative effect of each subsequent factor on virulence. None of the remaining 26 putative virulence determinants were associated with disease on univariate analysis, either separately or combined, as a count of the number of determinants in each organism.

(A) Association between the number of virulence-associated determinants and the proportion of isolates from cases of disease (rather than carriage). For example, 20% of isolates with one virulence-associated gene were from cases of disease and 80% were carriage strains. The area of each circle is proportional to the number of isolates with that number of determinants (ranging from zero genes to seven genes, with n = 2, 15, 21, 64, 140, 56, 22, and 7, respectively). (B) Association between the number of virulence-associated determinants in a gene combination and the proportion of isolates with that combination from cases of disease (rather than carriage). Each combination is identified by a letter; for the actual gene complement, refer to Table 4. The 12 most common combinations representing 80% of isolates are shown. The area of each circle is proportional to the number of isolates in that gene combination (n for each gene combination is given in Table 4).

Gene combinations and virulence.Given that the effects of the seven genes on the odds of being a disease-causing strain are additive, we examined whether particular gene combinations were either (i) more commonly found together in an organism than by chance or (ii) more virulent.

In the study population, there were 46 (of a possible 128) different combinations of the seven genes, with more than 80% of the isolates accounted for by 12 of these gene combinations (Table 4). We compared the actual frequencies with the predicted frequencies of each gene combination, assuming independent segregation of the genes. One common combination (combination B [Table 4]) was overrepresented in both the carriage and disease groups (19.4% observed versus 11.7% predicted in carriage isolates and 18.6% observed versus 10% predicted in disease isolates), two virulent combinations (combinations F and G, each with six virulence-associated genes) were overrepresented in the disease group (5.2% observed versus 1.3% predicted and 5.2% observed versus 1.1% predicted, respectively), and one less virulent combination (combination J, with one virulence-associated gene) was overrepresented in the carriage group (4.1% observed versus 0.9% predicted). Comparison of observed and predicted frequencies demonstrated that no gene combinations were significantly underrepresented. Are these discrepancies due to the effect of clonality, or do particular combinations of genes have either less or greater virulence than expected? None of the 12 most common combinations were either more or less likely to cause disease than others with the same number of virulence-associated genes, suggesting the effect was due to the clonal population structure (Fig. 2B).

Effect of clonality on genes and gene combinations. (i) Distribution of individual genes between lineages.Using contingency table analysis, we examined whether individual genes were apparently randomly distributed between the major MLST-derived clonal complexes or whether their distribution was significantly influenced by the underlying clonality of the population. clfB, sdrC, and see were excluded from this analysis, as they were present or absent in all isolates. The genes fell into two distinct groups. The distribution of the enterotoxins, tst, the exfoliatins (eta and etb), beta- and delta-toxins, the variable sdr genes (sdrD, sdrE, and bbp), cna, ebpS, and efb within the population were all highly significantly related to clonal complex (P < 0.001 in all cases). In contrast, the distributions of fnbA, clfA, coagulase, spa, map, pvl, hlg, alpha-toxin, ica, and V8 appeared to be unrelated to the underlying clonal structure (P > 0.1 in all cases). Most in this second group are genes which are either very common or very rare, and it is possible that the absence of a significant association with clonal complex is due to lack of power, leading to a type II error. However, the difference in the P value range between the two groups is striking, and even after removing from the second group those genes which are either present or absent in more than 95% of isolates, fnbA, spa, ica, map, and hlg remain apparently unrelated to underlying clonality (and three of these are virulence-associated genes).

Although there is strong statistical evidence for an association with clonality for the first group of genes listed above, there is also evidence to suggest that horizontal transfer of genes plays an important role in determining the distribution of these genes. All the genes with a prevalence of more than 30% were found in isolates from all 11 major lineages. The exceptions to this diversity are agr subgroups I to IV, which were very tightly (but not completely) linked to the underlying clonal structure (Table 5).

The most clonal of all the genes characterized as either present or absent was tst, with 92% of positive isolates occurring in clonal complexes CC30 or CC39. These two clonal complexes are closely related (the only two which are) and are probably derived from the same progenitor strain. One interpretation is that CC30 and CC39 represent a relatively recent but highly successful clonal expansion (perhaps fueled by tst) which has had less time to diversify through horizontal transfer than other lineages. That other genes (sea and bbp) were also heavily polarized with respect to this complex lends support to this view. The distributions of the seven virulence-related genes plus sea and tst are shown in Fig. 3; this selection includes genes which are chromosomally encoded (ica, fnbA, and sdrE) and associated with plasmids (sej), bacteriophage (sea), and the SaPI 1 pathogenicity island (tst).

Distribution of genes within the seven largest MLST-defined clonal complexes (denoted by numbers along the x axes in the bottom row). The overall height of each bar denotes the total number of isolates in the complex. The height of the shaded area represents the number of isolates positive for the determinant. The population genetic structure of S. aureus can be found on the MLST website (http://www.mlst.net
).

(ii) Distribution of gene combinations.Underlying clonality had an important effect on the distribution of certain combinations of genes (linkage disequilibrium). The overrepresentation in both disease and carriage populations of the virulence-associated gene combination B (containing ica, fnbA, hlg, and cna but lacking sdrE, sej, and eta) (Table 4) is entirely because it was present in 51% of the isolates from clonal complexes CC30 and CC39; its frequency in other clonal complexes was only 6%, not significantly different from its predicted frequency. However, this clonality was not the explanation for the association of these genes with virulence; combination B was not itself associated with disease, and three of the four genes (ica, fnbA, and hlg) were significantly associated with disease within CC30/C39, even after excluding those isolates with combination B. Despite the underlying clonality, there was still considerable variation within clonal complexes and even within MLST sequence types. The 22 isolates of nosocomial disease-causing EMRSA-16 (ST36, a single MLST) were represented by 10 different gene combinations, despite all being isolated from within the same city over a 2-year period (data not shown). ST36 was not unique in showing this diversity; the five most common combinations were all present in over half of the MLST-defined clonal complexes (data not shown), suggesting that horizontal transfer of genes between lineages is a relatively frequent occurrence.

DISCUSSION

These data represent a comprehensive analysis of virulence determinants in a large, natural population of S. aureus isolates of clearly defined clinical provenance. Our finding that several bacterial factors were associated with invasive disease may directly reflect their involvement in the pathogenesis of severe disease. It is unlikely that one or more of the determinants have become overrepresented due to the hitchhiker effect, as the positive associations remained despite stratification of the analysis by genetic lineage. This was demonstrated against a background in which many determinants were highly clonal, as shown here for several determinants.

The study was assumption free with respect to the presumed function and biological effect during infection of the bacterial factors examined. However, having identified seven virulence-associated genes, we now consider the biological plausibility of our findings. Fibronectin, the host ligand for the adhesin encoded by fnbA, is a ubiquitous host protein present in soluble form in the blood and in fibrillar form in cellular matrices, bacterial adhesion to which may be important during several steps of the disease process (31). Likewise, toxins that modulate the immune response through superantigen activity probably make a significant contribution to disease manifestations. However, it is also important to be somewhat circumspect when ascribing a given function to the pathogenic process for the following reasons. First, many bacterial determinants are multifunctional (for example, fibronectin binding protein has recently been found to be an adhesin for fibrinogen) (39) and the critical interaction of a given factor with the host may not yet have been defined. Second, the functions ascribed to a determinant by in vitro assays may not accurately reflect their behavior in vivo.

Our study design enriched the chances of finding virulence-associated determinants by examining those previously thought to be involved in disease and excluding those with purely housekeeping functions. However, it is unlikely that we have identified all the genes that are responsible for disease manifestations, an enterprise well suited to microarray analysis. In addition, some of the non-virulence-associated determinants as defined here may indeed play a role in pathogenesis but, because of their essential nature, were not flagged as virulence associated by our study design. For example, genes encoding the fibrinogen binding proteins ClfA and ClfB were ubiquitous, regardless of the origin of the strain.

Studies of staphylococcal pathogenesis often focus on the presence or absence of a given determinant. Our study results were analyzed on the same basis. However, the nature of the methodology used here means that we cannot differentiate between the presence or absence of the entire gene and the presence or absence of an allelic variant at a polymorphic locus. Southern hybridization may have been more sensitive than PCR in determining the presence or absence of a given gene. However, we believe that allelic variation may be important in defining virulence and that as such, PCR provides interesting clues. For example, bbp and sdrE appear to be allelic variants of each other yet only sdrE was associated with virulence. This observation is currently undergoing further study.

The variation in virulence determinants in a given clone (as demonstrated for EMRSA-16), together with the widespread nature of combinations of virulence determinants, suggests that horizontal transfer of genes is a common event. The preservation of the associations between certain genes and virulence after stratification for bacterial lineage is in itself strong evidence for either loss or gain through horizontal transfer of genetic elements. The mechanism of transfer can be readily explained for mobile determinants such as sej, which is plasmid mediated (40). However, the majority of the virulence-associated genes are chromosomal and it is not clear how such genes are being transmitted. In addition, not all mobile elements appeared to undergo frequent horizontal transmission. An example of this is tst, which resides on a pathogenicity island termed SaPI 1 (35). This gene was very common (>60%) in two closely related clonal complexes in our study (CC30 and CC39) but appeared either very infrequently or not at all in the remainder.

S. aureus usually behaves as a harmless human commensal, as reflected by the large proportion of the healthy population who sometimes or usually carry this organism, the transition to disease representing the exception rather than the rule. So what is the relationship between virulence-associated genes and bacterial fitness? The answer to this question is unknown but is likely to be complex and could vary depending on the factor in question. For example, a determinant that confers the ability to invade the host may be associated with a fitness disadvantage since invasion could lead to bacterial death following antibiotic treatment or death of the host. Conversely, the presence of genes that enhance the ability to cause conditions such as superficial skin infection could lead to enhanced host-to-host transmission. A further consideration is that accessory genes that are important to disease pathogenesis but which are not critical to survival may impose a cost to fitness through additional gene replication and protein secretion if expressed during periods of carriage. We postulate, therefore, that the proportion of strains circulating in the community which carry virulence determinants and the number of virulence-associated genes carried per strain are a product of the interplay between rates of gene acquisition, the cost to biological fitness, and the rate of decay of strains causing human disease.

Host factors for S. aureus disease are likely to include a genetic predisposition via one or more susceptibility genes and acquired factors such as the presence of intravenous devices, surgical wounds, and other events that perturb normal host defenses. It seems plausible to speculate that invasive disease occurs in two overlapping host populations. Individuals affected in the community (where the rate of disease is low) may represent those most genetically predisposed to S. aureus infection, while individuals infected in the hospital (where the burden of disease is higher) may represent a much larger at-risk group who may or may not carry susceptibility genes. Given this scenario, it is possible that the pattern or number of bacterial determinants associated with disease in the hospitalized host would differ from that seen in the community. However, we found little difference between strains associated with disease in the two settings, suggesting that bacterial factors play a role in causing disease, even in the compromised host. Whether this is the case for truly immunocompromised individuals such as neutropenics and those with AIDS requires further study. A larger study would also be required to address whether a given determinant(s) is associated with a particular clinical syndrome(s).

In conclusion, this study has demonstrated the variable presence of virulence genes in natural populations of S. aureus, providing evidence that bacterial factors play a role in determining invasive disease in both community and hospital settings. We have demonstrated that the effect of these genes was cumulative, each independently multiplying the odds of disease. There was also evidence of considerable horizontal transfer of genes on a background of clonality. This study also indicates that it may be an oversimplification to consider virulence in relation to the presence or absence of a given bacterial factor. Our findings suggest that allelic variants of a polymorphic locus can make different contributions to the disease process, further study of which is likely to expand our understanding of staphylococcal disease pathogenesis.

ACKNOWLEDGMENTS

This project was supported by Royal Society grant 20955 and Wellcome Trust grant 059064 to S. Peacock and Wellcome Trust Career Development award 049310 to N. Day.

We are grateful to B. Cameron for technical assistance during pilot studies.

Hemker, H. C., B. M. Bas, and A. D. Muller.1975. Activation of a pro-enzyme by a stoichiometric reaction with another protein. The reaction between prothrombin and staphylocoagulase. Biochim. Biophys. Acta379:180-188.

Supersac, G., G. Prevost, and Y. Piemont.1993. Sequencing of leucocidin R from Staphylococcus aureus P83 suggests that staphylococcal leucocidins and gamma-hemolysin are members of a single, two-component family of toxins. Infect. Immun.61:580-587.