Core-Genome Multilocus Sequence Typing for Epidemiological and Evolutionary Analyses of Phytopathogenic Xanthomonas citri.

Xanthomonas citri subsp. citri is the cause of bacterial citrus canker, responsible for major economic losses to the citrus industry. X. citri subspecies and pathovars are responsible for diseases in soybean, common bean, mango, pomegranate, and cashew. X. citri disease has been tracked using several typing methods, but recent studies using genomic sequencing have been key to understanding the evolutionary relationships within the species, including fundamental differences among X. citri subsp. citri pathotypes. Here, we describe a core-genome multilocus sequence typing (cgMLST) scheme for X. citri based on 250 genomes comprising multiple examples of X. citri subsp. citri pathotypes A, A*, and Aw; X. citri subsp. malvacearum; X. citri pv. aurantifolii, pv. fuscans, pv. glycines, pv. mangiferaeindicae, pv. viticola, and pv. vignicola; and single isolates of X. citri pv. dieffenbachiae and pv. punicae. This data set included genomic sequencing of 100 novel X. citri subsp. citri isolates. cgMLST, based on 1,618 core genes across 250 genomes, is implemented at PubMLST (https://pubmlst.org/organisms/xanthomonas-citri/). GrapeTree minimum-spanning tree and Interactive Tree of Life (iTOL) neighbor-joining phylogenies generated from the cgMLST data resolved almost identical groupings of isolates to a core-genome single nucleotide polymorphism (SNP)-based neighbor-joining phylogeny. These resolved identical groupings of X. citri subsp. citri pathotypes and X. citri subspecies and pathovars. X. citri cgMLST should prove to be an increasingly valuable resource for the study of this key species of plant-pathogenic bacteria. Users can submit genomic data and associated metadata for comparison with previously characterized isolates at PubMLST to allow the rapid characterization of the local, national, and global epidemiology of these pathogens and examine evolutionary relationships. IMPORTANCE Xanthomonas citri is a plant pathogen that causes major economic losses to the citrus industry and sweet orange production in particular. Several subspecies and pathogens are recognized, with host ranges including soybean, common bean, mango, pomegranate, and cashew, among others. Recent genomic studies have shown that host-adapted X. citri subspecies and pathovars and X. citri subsp. citri pathotypes form distinct clades. In this study, we describe a core-genome multilocus sequence typing (cgMLST) scheme for this species that can rapidly and robustly discriminate among these ecologically distinct, host-adapted clades. We have established this scheme and associated databases containing genomic sequences and metadata at PubMLST, which users can interrogate with their own genome sequences to determine X. citri subspecies, pathovars, and pathotypes. X. citri cgMLST should prove to be an invaluable tool for the study of the epidemiology and evolution of this major plant pathogen.

Assessing phage therapy against Pseudomonas aeruginosa using a Galleria mellonella infection model.

The Galleria mellonella infection model was used to assess the in vivo efficacy of phage therapy against laboratory and clinical strains of Pseudomonas aeruginosa. In a first series of experiments, Galleria were infected with the laboratory strain P. aeruginosa PAO1 and were treated with varying multiplicity of infection (MOI) of phages either 2h post-infection (treatment) or 2h pre-infection (prevention) via injection into the haemolymph. To address the kinetics of infection, larvae were bled over a period of 24h for quantification of bacteria and phages. Survival rates at 24h when infected with 10 cells/larvae were greater in the prevention versus treatment model (47% vs. 40%, MOI=10; 47% vs. 20%, MOI=1; and 33% vs. 7%, MOI=0.1). This pattern held true when 100 cells/larvae were used (87% vs. 20%, MOI=10; 53% vs. 13%, MOI=1; 67% vs. 7%, MOI=0.1). By 24h post-infection, phages kept bacterial cell numbers in the haemolymph 1000-fold lower than in the non-treated group. In a second series of experiments using clinical strains to further validate the prevention model, phages protected Galleria when infected with both a bacteraemia (0% vs. 85%) and a cystic fibrosis (80% vs. 100%) isolate. Therefore, this study validates the use of G. mellonella as a simple, robust and cost-effective model for initial in vivo examination of P. aeruginosa-targeted phage therapy, which may be applied to other pathogens with similarly low infective doses.

Combined use of bacteriophage K and a novel bacteriophage to reduce Staphylococcus aureus biofilm formation.

Biofilms are major causes of impairment of wound healing and patient morbidity. One of the most common and aggressive wound pathogens is Staphylococcus aureus, displaying a large repertoire of virulence factors and commonly reduced susceptibility to antibiotics, such as the spread of methicillin-resistant S. aureus (MRSA). Bacteriophages are obligate parasites of bacteria. They multiply intracellularly and lyse their bacterial host, releasing their progeny. We isolated a novel phage, DRA88, which has a broad host range among S. aureus bacteria. Morphologically, the phage belongs to the Myoviridae family and comprises a large double-stranded DNA (dsDNA) genome of 141,907 bp. DRA88 was mixed with phage K to produce a high-titer mixture that showed strong lytic activity against a wide range of S. aureus isolates, including representatives of the major international MRSA clones and coagulase-negative Staphylococcus. Its efficacy was assessed both in planktonic cultures and when treating established biofilms produced by three different biofilm-producing S. aureus isolates. A significant reduction of biofilm biomass over 48 h of treatment was recorded in all cases. The phage mixture may form the basis of an effective treatment for infections caused by S. aureus biofilms.

Bacteriophages for the treatment of Pseudomonas aeruginosa infections.

Bacteriophages were first identified in 1915 and were used as antimicrobial agents from 1919 onwards. Despite apparent successes and widespread application, early users did not understand the nature of these agents and their efficacy remained controversial. As a result, they were replaced in the west by chemical antibiotics once these became available. However, bacteriophages remained a common therapeutic approach in parts of Eastern Europe where they are still in use. Increasing levels of antibiotic-resistant bacterial infections are now driving demand for novel therapeutic approaches. In cases where antibiotic options are limited or nonexistent, the pressure for new agents is greatest. One of the most prominent areas of concern is multidrug-resistant Gram-negative bacteria. Pseudomonas aeruginosa is a prominent member of this class and is the cause of damaging infections that can be resistant to successful treatment with conventional antibiotics. At the same time, it exhibits a number of properties that make it a suitable target for bacteriophage-based approaches, including growth in biofilms that can hydrolyse following phage infection. Pseudomonas aeruginosa provides a striking example of an infection where clinical need and the availability of a practical therapy coincide.

Phage therapy: delivering on the promise.

Bacteriophages are viruses that infect and, in many cases, destroy their bacterial targets. Within a few years of their initial discovery they were being investigated as therapeutic agents for infectious disease, an approach known as phage therapy. However, the nature of these exquisitely specific agents was not understood and much early use was both uninformed and unsuccessful. As a result they were replaced by chemical antibiotics once these became available. Although work on phage therapy continued (and continues) in Eastern Europe, this was not conducted to a standard allowing it to support clinical uses in areas regulated by the European Medicines Agency or the US FDA. To develop phage therapy for these areas requires work carried out in accordance with the requirements of these agencies, and, driven by the current crisis of antibiotic resistance, such clinical trials are now under way. The first Phase I clinical trial of safety was reported in 2005, and the results of the first Phase II clinical trial of efficacy of a bacteriophage therapeutic was published in 2009. While the delivery of these relatively large and complex agents to the site of disease can be more challenging than for conventional, small-molecule antibiotics, bacteriophages are then able to multiply locally even from an extremely low (picogram range) initial dose. This multiplication where and only where they are needed underlies the potential for bacteriophage therapeutics to become a much needed and powerful weapon against bacterial disease.

Comparison of molecular typing methods for characterization of Staphylococcus epidermidis: proposal for clone definition.

In the present study we give some direction on the selection of the most appropriate typing method(s) to be used for the characterization of Staphylococcus epidermidis, in view of the most recent findings on the evolution, population structure, and epidemiology of this species. In order to achieve this aim, quantitative assessment of the correlation of the results of three typing methods--pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), and staphylococcal chromosomal cassette mec (SCCmec) typing, which target different regions of the chromosome that evolve at different rates--was performed. In order to evaluate the discriminatory ability and the strength and direction of the correlation of the different typing methods, Simpson's index of diversity (SID), the adjusted Rand coefficient (AR), and the Wallace coefficient (W) were calculated. PFGE was the most discriminatory method (SID = 99%), followed by MLST (SID = 90%) and SCCmec typing (SID = 75%). The values of AR and W (0.10 < AR < 0.30; 0.50 < W < 0.75) indicated that the partition of the same isolate collection by PFGE, MLST, and SCCmec typing provided results that had only a poor correlation with each other. However, the information provided by the combination of PFGE and SCCmec enabled the prediction of the results obtained by MLST at the level of the clonal complex with a high degree of precision (W > 0.90). We propose that clones of S. epidermidis be defined by the combination of the PFGE type followed by the SCCmec type, which provides reliable information on the short-term epidemiology and the ability to predict with consistency long-term clonal evolution.

Inferring a population structure for Staphylococcus epidermidis from multilocus sequence typing data.

Despite its importance as a human pathogen, information on population structure and global epidemiology of Staphylococcus epidermidis is scarce and the relative importance of the mechanisms contributing to clonal diversification is unknown. In this study, we addressed these issues by analyzing a representative collection of S. epidermidis isolates from diverse geographic and clinical origins using multilocus sequence typing (MLST). Additionally, we characterized the mobile element (SCCmec) carrying the genetic determinant of methicillin resistance. The 217 S. epidermidis isolates from our collection were split by MLST into 74 types, suggesting a high level of genetic diversity. Analysis of MLST data using the eBURST algorithm revealed the existence of nine epidemic clonal lineages that were disseminated worldwide. One single clonal lineage (clonal complex 2) comprised 74% of the isolates, whereas the remaining isolates were clustered into 8 minor clonal lineages and 13 singletons. According to our evolutionary model, SCCmec was acquired at least 56 times by S. epidermidis. Although geographic dissemination of S. epidermidis strains and the value of the index of association between the alleles, 0.2898 (P < 0.05), support the clonality of S. epidermidis species, examination of the sequence changes at MLST loci during clonal diversification showed that recombination gives rise to new alleles approximately twice as frequently as point mutations. We suggest that S. epidermidis has a population with an epidemic structure, in which nine clones have emerged upon a recombining background and evolved quickly through frequent transfer of genetic mobile elements, including SCCmec.

Staphylococcal chromosome cassette evolution in Staphylococcus aureus inferred from ccr gene complex sequence typing analysis.

Staphylococcal chromosome cassette mec (SCCmec) elements within major lineages of healthcare- and community-associated methicillin-resistant Staphylococcus aureus (MRSA) clones were characterised using intra-SCCmec multilocus sequencing. A strong correlation was observed between sequence- and PCR-based typing methods (p <0.001). However, phylogenetic analysis of the SCCmec locus using concatenated sequences evidenced few recombination events. Sequence type (ST)-SCCmec1 was found in SCCmec elements types I and IV, suggesting the evolution of an SCCmecI element into an SCCmecIV element. This coincided with the spread of the clone harbouring this SCCmec element into the community. No correlation was observed between ST-SCCmec lineage and MRSA lineage, confirming multiple acquisitions of SCCmec by S. aureus. This was exemplified by the SCCmecIV ST-SCCmec10 element, which was detected in all of the clonal complexes examined, including healthcare- and community-associated MRSA. The acquisition of this SCCmec element was five- to ten-fold more common than that of others. Models of MRSA clone evolution suggest that this SCCmec was first found in the paediatric clone.

High interlaboratory reproducibility of DNA sequence-based typing of bacteria in a multicenter study.

Current DNA amplification-based typing methods for bacterial pathogens often lack interlaboratory reproducibility. In this international study, DNA sequence-based typing of the Staphylococcus aureus protein A gene (spa, 110 to 422 bp) showed 100% intra- and interlaboratory reproducibility without extensive harmonization of protocols for 30 blind-coded S. aureus DNA samples sent to 10 laboratories. Specialized software for automated sequence analysis ensured a common typing nomenclature.

Diversity among community isolates of methicillin-resistant Staphylococcus aureus in Australia.

Community methicillin-resistant Staphylococcus aureus (CMRSA) strains are being isolated with increasing frequency around the world. In Western Australia CMRSA are endemic in geographically remote communities and have been found to belong to five different contour-clamped homogeneous electric field (CHEF) electrophoretic patterns. Representatives of each of these CHEF patterns have been compared to CMRSA representative of CHEF patterns from other Australian states and New Zealand. With one exception, all of the isolates were nonmultiresistant and were not resistant to many antimicrobial agents other than the beta-lactams. With one exception, which is not believed to be a CMRSA, all of the isolates harbored a beta-lactamase plasmid. Erythromycin resistance was associated with a 2-kb plasmid. One of the beta-lactamase plasmids was found to be able to acquire additional resistance determinants to become a multiple resistance plasmid. There were 10 multilocus sequence types belonging to eight distantly related clonal complexes of S. aureus. One new sequence type was found. Although most of the CMRSA harbored the type IVa SCCmec, a type IV structural variant was found and two new SCCmec types were identified. Protein A gene (spa) typing revealed two new spa types and, with two exceptions, corresponded to multilocus sequence typing. In contrast to other reports on CMRSA, most of the CMRSA strains studied here did not contain the Panton-Valentine leukocidin genes. The results also demonstrate that nonmultiresistant hospital strains such as UK EMRSA-15 may be able to circulate in the community and could be mistaken for CMRSA based on their resistance profiles.

Multilocus sequence typing and the evolution of methicillin-resistant Staphylococcus aureus.

The prevalence of methicillin-resistant Staphylococcus aureus (MRSA) in many countries is increasing and, in hospitals in some areas, more than half of all S. aureus disease isolates are MRSA. MRSA strains are becoming increasingly multiresistant, and have recently developed resistance to vancomycin, used successfully to treat MRSA for more than 30 years. This review summarises recent studies that have elucidated the evolutionary history of MRSA. The first MRSA isolate evolved from a sensitive, epidemic strain prevalent in Europe, and its progeny-the first MRSA clone-quickly spread to other continents. Analyses of epidemic MRSA isolates from hospitals in different countries by molecular methods, including multilocus sequence typing (MLST) and DNA microarray analysis, reveal that MRSA strains have evolved separately within five distinct epidemic, sensitive lineages. However, resistance has been transferred to S. aureus on many more than five occasions, as some lineages have acquired different structural types of the element carrying the methicillin resistance gene. The emergence of MRSA as a community pathogen has been noted in several countries, and MLST and SCCmec typing have been used to demonstrate that community-acquired MRSA strains are typically related only distantly to hospital MRSA strains, and thus represent novel acquisitions of SCCmec.

Characterization of isolates of methicillin-resistant Staphylococcus aureus from Hong Kong by phage typing, pulsed-field gel electrophoresis, and fluorescent amplified-fragment length polymorphism analysis.

The genetic relatedness of 127 methicillin-resistant Staphylococcus aureus (MRSA) isolates, belonging to five major types as identified by pulsed-field gel electrophoresis (PFGE) and antibiotic resistance profiles, was examined further using phage typing and fluorescent amplified fragment length polymorphism (FAFLP). The MRSA isolates were recovered from patients at the Prince of Wales Hospital (PWH), Hong Kong, over a 13-year period, 1988 to 2000. These strains were also compared with representatives of the well-described MRSA international clones and with epidemic MRSA strains (eMRSA) 1 to 16 from the United Kingdom. Phage typing distinguished two major "clones" at this hospital: all of the phage type 1 (PT1) isolates belonged to PFGE types A, C, D, and E, while most of the PT2 isolates were associated with PFGE type B, which exhibited a unique antibiotic resistance profile. MRSA isolates belonging to PFGE subtype A2 were indistinguishable from the British eMRSA-1, while isolates of PFGE type B were closely related to eMRSA-9 by PFGE. Based on FAFLP, all five predominant PFGE types at the PWH belonged to one group and fell into the same cluster as eMRSA-1, -4, -7, -9, and -11 isolates. Multilocus sequence typing and staphylococcal cassette chromosome mec typing classified representatives of our MRSA isolates as members of the same clone (ST239-MRSA-III). Thus, the predominant MRSA isolates frin the PWH in the last decade are closely related to early United Kingdom eMRSA clones 1, 4, and 11 and are members of a lineage that includes the Brazilian MRSA clone.

Worldwide incidence, molecular epidemiology and mutations implicated in fluoroquinolone-resistant Streptococcus pneumoniae: data from the global PROTEKT surveillance programme.

OBJECTIVES: To analyse the mutations and epidemiology associated with fluoroquinolone-resistant pneumococci collected as part of the PROTEKT global surveillance programme during 1999-2000. METHODS: Sixty-nine centres in 25 countries submitted a total of 3362 Streptococcus pneumoniae isolates, for which the MICs of antimicrobial agents were determined using NCCLS methodology. RESULTS: Levofloxacin resistance was low overall (1% worldwide), with higher rates in: Hong Kong (14.3%), South Korea (2.9%), USA (1.8%), Mexico (1.5%), Canada (1.4%) and Japan (1.3%). Levofloxacin resistance was very low or absent in European countries, and absent in Australia. Worldwide, there was a total of 35 levofloxacin-resistant isolates, of which 22 (63%) were resistant and 10 (29%) were intermediate to moxifloxacin. All levofloxacin-resistant isolates were susceptible to telithromycin (< or =0.5 mg/L), linezolid (< or =2 mg/L) and quinupristin/dalfopristin (< or =1 mg/L). One or more mutations in the topoisomerase genes were identified in all levofloxacin-resistant isolates; most of these isolates (33/35) had a mutation in one of the DNA gyrase encoding genes (gyrA, gyrB) and one of the topoisomerase IV encoding genes (parC, parE). Eighteen (51%) isolates carried the same combination of amino acid substitutions: Ser-81-->Phe in GyrA and Ser-79-->Phe in ParC. Isolates displaying a levofloxacin MIC of 2-4 mg/L generally had no mutation or one mutation in either a DNA gyrase or a topoisomerase IV gene, although most mutations were in parC. CONCLUSIONS: Most levofloxacin-resistant isolates possess two mutations (one in DNA gyrase and one in topoisomerase IV). Although multilocus sequence typing demonstrated that most of these isolates were unrelated, 12 (34%) were the Spain23F-1 clone: 10 from Hong Kong and one each from Saskatchewan, Canada and Sao Paulo, Brazil.

Directional gene movement from human-pathogenic to commensal-like streptococci.

Group A streptococci (GAS) are highly pathogenic for humans, and their closest genetic relatives, group C and G streptococci (GCS and GGS, respectively), are generally regarded as commensals, although they can be found in association with human disease. As part of an effort to better understand the evolution of virulence, the phylogenetic relationships between GAS, GCS, and GGS were examined. The nucleotide sequence was determined for an internal portion of seven housekeeping (neutral) loci among >200 isolates of GAS and 34 isolates of GCS or GGS obtained from human subjects. Genotypic analysis failed to show support for the separation of GCS and GGS into two distinct populations. Unlike GAS, there was poor concordance between emm type and genetic relatedness among GCS and GGS. All housekeeping genes within GAS displayed relatively low levels of sequence diversity. In contrast, individual GCS and GGS strains had mosaic genomes, containing alleles at some loci that were similar or identical to GAS alleles, whereas the alleles at other loci were about 10 to 30% diverged. The data provide evidence for a history of recent interspecies transfer of neutral genes that exhibits a strong net directionality from GAS donors to GCS and GGS recipients. A model for the evolution of GAS and of GCS and GGS is described.

A link between virulence and ecological abundance in natural populations of Staphylococcus aureus.

Staphylococcus aureus is a major cause of severe infection in humans and yet is carried without symptoms by a large proportion of the population. We used multilocus sequence typing to characterize isolates of S. aureus recovered from asymptomatic nasal carriage and from episodes of severe disease within a defined population. We identified a number of frequently carried genotypes that were disproportionately common as causes of disease, even taking into account their relative abundance among carriage isolates. The existence of these ecologically abundant hypervirulent clones suggests that factors promoting the ecological fitness of this important pathogen also increase its virulence.

Multilocus sequence typing of Streptococcus pyogenes and the relationships between emm type and clone.

Multilocus sequence typing (MLST) is a tool that can be used to study the molecular epidemiology and population genetic structure of microorganisms. A MLST scheme was developed for Streptococcus pyogenes and the nucleotide sequences of internal fragments of seven selected housekeeping loci were obtained for 212 isolates. A total of 100 unique combinations of housekeeping alleles (allelic profiles) were identified. The MLST scheme was highly concordant with several other typing methods. The emm type, corresponding to a locus that is subject to host immune selection, was determined for each isolate; of the >150 distinct emm types identified to date, 78 are represented in this report. For a given emm type, the majority of isolates shared five or more of the seven housekeeping alleles. Stable associations between emm type and MLST were documented by comparing isolates obtained decades apart and/or from different continents. For the 33 emm types for which more than one isolate was examined, only five emm types were present on widely divergent backgrounds, differing at four or more of the housekeeping loci. The findings indicate that the majority of emm types examined define clones or clonal complexes. In addition, an MLST database is made accessible to investigators who seek to characterize other isolates of this species via the internet (http://www.mlst.net).

Recombination within natural populations of pathogenic bacteria: short-term empirical estimates and long-term phylogenetic consequences.

The identification of clones within bacterial populations is often taken as evidence for a low rate of recombination, but the validity of this inference is rarely examined. We have used statistical tests of congruence between gene trees to examine the extent and significance of recombination in six bacterial pathogens. For Neisseria meningitidis, Streptococcus pneumoniae, Streptococcus pyogenes, and Staphylococcus aureus, the congruence between the maximum likelihood trees reconstructed using seven house-keeping genes was in most cases no better than that between each tree and trees of random topology. The lack of congruence between gene trees in these four species, which include both naturally transformable and nontransformable species, is in three cases supported by high ratios of recombination to point mutation during clonal diversification (estimates of this parameter were not possible for Strep. pyogenes). In contrast, gene trees constructed for Hemophilus influenzae and pathogenic isolates of Escherichia coli showed a higher degree of congruence, suggesting lower rates of recombination. The impact of recombination therefore varies between bacterial species but in many species is sufficient to obliterate the phylogenetic signal in gene trees.

Molecular typing of bacteria directly from cerebrospinal fluid.

Using Streptococcus pneumoniae as an example, the ability of multilocus sequence typing (MLST) to characterise isolates directly from cerebrospinal fluid (CSF) was investigated. A nested multiplex polymerase chain reaction method that amplifies the seven housekeeping gene fragments used for pneumococcal MLST was applied to 30 CSF samples from suspected cases of bacterial meningitis. The fragments were amplified from all 14 samples from which Streptococcus pneumoniae was cultured, and, after direct sequencing, the allelic profiles obtained from ten of the samples corresponded to those of clones previously associated with invasive pneumococcal disease. MLST could also predict the penicillin susceptibility and serotype of the CSF isolates.

Estimating the relative contributions of mutation and recombination to clonal diversification: a comparison between Neisseria meningitidis and Streptococcus pneumoniae.

Both Neisseria meningitidis and Streptococcus pneumoniae are naturally transformable species and are known to be freely recombining in the wild. Large multilocus sequence typing (MLST) datasets have been generated for these species. Here we outline an approach which exploits these data sets in order to quantify the extent of recombination, thus enabling meaningful comparisons between the two species. Two parameters are estimated; the rate at which recombination changes alleles, compared to point mutation, and the rate at which recombination changes individual nucleotide sites, compared to point mutation. Estimates for the former parameter are 4:1 in the meningococcus (i.e. alleles are changed four-fold more frequently by recombination than by mutation), and 10:1 in the pneumococcus. However, estimates for the latter parameter are at least 80:1 in the meningococcus (i.e. an individual nucleotide site is at least 80-fold more likely to change by recombination than by mutation) and 50:1 in the pneumococcus. These data imply that recombination events, compared to mutational events, may be more common in the pneumococcus than in the meningococcus. However, because it is a more diverse species, each recombinational exchange in the meningococcus results in more nucleotide changes on average.

Estimating recombinational parameters in Streptococcus pneumoniae from multilocus sequence typing data.

Multilocus sequence typing (MLST) is a highly discriminatory molecular typing method that defines isolates of bacterial pathogens using the sequences of approximately 450-bp internal fragments of seven housekeeping genes. This technique has been applied to 575 isolates of Streptococcus pneumoniae and identifies a number of discrete clonal complexes. These clonal complexes are typically represented by a single group of isolates sharing identical alleles at all seven loci, plus single-locus variants that differ from this group at only one out of the seven loci. As MLST is highly discriminatory, the members of each clonal complex can be assumed to have a recent common ancestor, and the molecular events that give rise to the single-locus variants can be used to estimate the relative contributions of recombination and mutation to clonal divergence. By comparing the sequences of the variant alleles within each clonal complex with the allele typically found within that clonal complex, we estimate that recombination has generated new alleles at a frequency approximately 10-fold higher than mutation, and that a single nucleotide site is approximately 50 times more likely to change through recombination than mutation. We also demonstrate how to estimate the average length of recombinational replacements from MLST data.