eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data

The introduction of multilocus sequence typing (MLST) for the precise characterization of isolates of bacterial pathogens has had a marked impact on both routine epidemiological surveillance and microbial population biology. In both fields, a key prerequisite for exploiting this resource is the ability to discern the relatedness and patterns of evolutionary descent among isolates with similar genotypes. Traditional clustering techniques, such as dendrograms, provide a very poor representation of recent evolutionary events, as they attempt to reconstruct relationships in the absence of a realistic model of the way in which bacterial clones emerge and diversify to form clonal complexes. An increasingly popular approach, called BURST, has been used as an alternative, but present implementations are unable to cope with very large data sets and offer crude graphical outputs. Here we present a new implementation of this algorithm, eBURST, which divides an MLST data set of any size into groups of related isolates and clonal complexes, predicts the founding (ancestral) genotype of each clonal complex, and computes the bootstrap support for the assignment. The most parsimonious patterns of descent of all isolates in each clonal complex from the predicted founder(s) are then displayed. The advantages of eBURST for exploring patterns of evolutionary descent are demonstrated with a number of examples, including the simple Spain(23F)-1 clonal complex of Streptococcus pneumoniae, "population snapshots" of the entire S. pneumoniae and Staphylococcus aureus MLST databases, and the more complicated clonal complexes observed for Campylobacter jejuni and Neisseria meningitidis.     

Global optimal eBURST analysis of multilocus typing data using a graphic matroid approach

Background  

Multilocus Sequence Typing (MLST) is a frequently used typing method for the analysis of the clonal relationships among strains of several clinically relevant microbial species. MLST is based on the sequence of housekeeping genes that result in each strain having a distinct numerical allelic profile, which is abbreviated to a unique identifier: the sequence type (ST). The relatedness between two strains can then be inferred by the differences between allelic profiles. For a more comprehensive analysis of the possible patterns of evolutionary descent, a set of rules were proposed and implemented in the eBURST algorithm. These rules allow the division of a data set into several clusters of related strains, dubbed clonal complexes, by implementing a simple model of clonal expansion and diversification. Within each clonal complex, the rules identify which links between STs correspond to the most probable pattern of descent. However, the eBURST algorithm is not globally optimized, which can result in links, within the clonal complexes, that violate the rules proposed.     

Assessing the reliability of eBURST using simulated populations with known ancestry

Background  

The program eBURST uses multilocus sequence typing data to divide bacterial populations into groups of closely related strains (clonal complexes), predicts the founding genotype of each group, and displays the patterns of recent evolutionary descent of all other strains in the group from the founder. The reliability of eBURST was evaluated using populations simulated with different levels of recombination in which the ancestry of all strains was known.     

Evidence of host-associated populations of Cryptosporidium parvum in Italy

Recent studies have revealed extensive genetic variation among isolates of Cryptosporidium parvum, an Apicomplexan parasite that causes gastroenteritis in both humans and animals worldwide. The parasite's population structure is influenced by the intensity of transmission, the host-parasite interaction, and husbandry practices. As a result, C. parvum populations can be panmictic, clonal, or even epidemic on both a local scale and a larger geographical scale. To extend the study of C. parvum populations to an unexplored region, 173 isolates of C. parvum collected in Italy from humans and livestock (calf, sheep, and goat) over a 10-year period were genotyped using a multilocus scheme based on 7 mini- and microsatellite loci. In agreement with other studies, extensive polymorphism was observed, with 102 distinct multilocus genotypes (MLGs) identified among 173 isolates. The presence of linkage disequilibrium, the confinement of MLGs to individual farms, and the relationship of many MLGs inferred using network analysis (eBURST) suggest a predominantly clonal population structure, but there is also evidence that part of the diversity can be explained by genetic exchange. MLGs from goats were found to differ from bovine and sheep MLGs, supporting the existence of C. parvum subpopulations. Finally, MLGs from isolates collected between 1997 and 1999 were also identified as a distinct subgroup in principal-component analysis and eBURST analysis, suggesting a continuous introduction of novel genotypes in the parasite population.     

Multilocus sequence typing reveals high genetic diversity and epidemic population structure for the fish pathogen Yersinia ruckeri

Yersinia ruckeri is the causative agent of enteric redmouth in fish and one of the major bacterial pathogens causing losses in salmonid aquaculture. Previously typing methods, including restriction enzyme analysis, pulsed-field gel electrophoresis and multilocus enzyme electrophoresis (MLEE) have indicated a clonal population structure. In this work, we describe a multilocus sequence typing (MLST) scheme for Y.ruckeri based on the internal fragment sequence of six housekeeping genes. This MLST scheme was applied to 103 Y.ruckeri strains from diverse geographic areas and hosts as well as environmental sources. Sequences obtained from this work were deposited and are available in a public database (http://publmst.org/yruckeri/). Thirty different sequence types (ST) were identified, 21 of which were represented by a single isolate, evidencing high genetic diversity. ST2 comprised more than one-third of the isolates and was most frequently observed among isolates from trout. Two major clonal complexes (CC) were identified by eBURST analysis showing a common evolutionary origin for 94 isolates forming 21 STs into CC1 and for 6 isolates of 6 STs in the CC2. It was also possible to associate some unique ST with isolates from recent outbreaks in vaccinated salmonid fish.     

The relative contributions of recombination and mutation to the divergence of clones of Neisseria meningitidis.

Multilocus sequence typing (MLST) is a recently developed nucleotide sequence-based method for the definitive assignment of isolates within bacterial populations to specific clones. MLST uses the same principles as multilocus enzyme electrophoresis and provides data that can be used to investigate aspects of the population genetics and evolution of bacterial species. We used an MLST data set consisting of the sequences of approximately 450-bp fragments from seven housekeeping loci from a large strain collection of Neisseria meningitidis to estimate the relative impact of recombination compared with point mutation in the diversification of N. meningitidis clonal complexes. 126 meningococcal isolates were assigned to 10 clonal complexes, 9 of which contained minor clonal variants. The allelic variation within each complex was classified as a recombinational exchange or a putative point mutation through a comparison of the sequences of each variant allele with that of the allele typically found in the clonal complex. The nine clonal complexes contained a total of 23 allelic variants, and analysis of the sequences of these variant alleles revealed that a single nucleotide site in a meningococcal housekeeping gene is at least 80-fold more likely to change as a result of recombination than as a result of mutation. This value is estimated to be 10-50-fold for Escherichia coli and approximately 50-fold for Streptococcus pneumoniae.     

Inferences about the Global Population Structures of Cryptosporidium parvum and Cryptosporidium hominis

Cryptosporidium parvum and Cryptosporidium hominis are two related species of apicomplexan protozoa responsible for the majority of human cases of cryptosporidiosis. In spite of their considerable public health impact, little is known about the population structures of these species. In this study, a battery of C. parvum and C. hominis isolates from seven countries was genotyped using a nine-locus DNA subtyping scheme. To assess the existence of geographical partitions, the multilocus genotype data were mined using a cluster analysis based on the nearest-neighbor principle. Within each country, the population genetic structures were explored by combining diversity statistical tests, linkage disequilibrium, and eBURST analysis. For both parasite species, a quasi-complete phylogenetic segregation was observed among the countries. Cluster analysis accurately identified recently introduced isolates. Rather than conforming to a strict paradigm of either a clonal or a panmictic population structure, data are consistent with a flexible reproductive strategy characterized by the cooccurrence of both propagation patterns. The relative contribution of each pattern appears to vary between the regions, perhaps dependent on the prevailing ecological determinants of transmission.     

Use of amplified-fragment length polymorphism to study the ecology of Campylobacter jejuni in environmental water and to predict multilocus sequence typing clonal complexes

We determined the genetic variability among water isolates of Campylobacter jejuni by using amplified-fragment length polymorphism (AFLP) and multilocus sequence typing (MLST). Across a highly diverse collection of isolates, AFLP clusters did not correlate with MLST clonal complexes, suggesting that AFLP is not reliable for deciphering population genetic relationships and may be problematic for larger epidemiologic analyses.     

Population Structure of Clinical Vibrio parahaemolyticus from 17 Coastal Countries,Determined through Multilocus Sequence Analysis

Vibrio parahaemolyticus is a leading cause of food-borne gastroenteritis worldwide. Although this bacterium has been the subject of much research, the population structure of clinical strains from worldwide collections remains largely undescribed, and the recorded outbreaks of V. parahaemolyticus gastroenteritis highlight the need for the subtyping of this species. We present a broad phylogenetic analysis of 490 clinical V. parahaemolyticus isolates from 17 coastal countries through multilocus sequence analysis (MLST). The 490 tested isolates fell into 161 sequence types (STs). The eBURST algorithm revealed that the 161 clinically relevant STs belonged to 8 clonal complexes, 11 doublets, and 94 singletons, showing a high level of genetic diversity. CC3 was found to be a global epidemic clone of V. parahaemolyticus, and ST-3 was the only ST with an international distribution. recA was observed to be evolving more rapidly, exhibiting the highest degree of nucleotide diversity (0.028) and the largest number of polymorphic nucleotide sites (177). We also found that the high variability of recA was an important cause of differences between the results of the eBURST and ME tree analyses, suggesting that recA has a much greater influence on the apparent evolutionary classification of V. parahaemolyticus based on the current MLST scheme. In conclusion, it is evident that a high degree of genetic diversity within the V. parahaemolyticus population and multiple sequence types are contributing to the burden of disease around the world. MLST, with a fully extractable database, is a powerful system for analysis of the clonal relationships of strains at a global scale. With the addition of more strains, the pubMLST database will provide more detailed and accurate information, which will be conducive to our future research on the population structure of V. parahaemolyticus.     

Evolutionary genetics of the capsular locus of serogroup 6 pneumococci

The evolution of the capsular biosynthetic (cps) locus of serogroup 6 Streptococcus pneumoniae was investigated by analyzing sequence variation within three serotype-specific cps genes from 102 serotype 6A and 6B isolates. Sequence variation within these cps genes was related to the genetic relatedness of the isolates, determined by multilocus sequence typing, and to the inferred patterns of recent evolutionary descent, explored using the eBURST algorithm. The serotype-specific cps genes had a low percent G+C, and there was a low level of sequence diversity in this region among serotype 6A and 6B isolates. There was also little sequence divergence between these serotypes, suggesting a single introduction of an ancestral cps sequence, followed by slight divergence to create serotypes 6A and 6B. A minority of serotype 6B isolates had cps sequences (class 2 sequences) that were approximately 5% divergent from those of other serotype 6B isolates (class 1 sequences) and which may have arisen by a second, more recent introduction from a related but distinct source. Expression of a serotype 6A or 6B capsule correlated perfectly with a single nonsynonymous polymorphism within wciP, the rhamnosyl transferase gene. In addition to ample evidence of the horizontal transfer of the serotype 6A and 6B cps locus into unrelated lineages, there was evidence for relatively frequent changes from serotype 6A to 6B, and vice versa, among very closely related isolates and examples of recent recombinational events between class 1 and 2 cps serogroup 6 sequences.     

Application of Streptococcus uberis multilocus sequence typing: analysis of the population structure detected among environmental and bovine isolates from New Zealand and the United Kingdom

We recently developed a multilocus sequence typing (MLST) scheme to differentiate S. uberis isolates and facilitate an understanding of the population biology of this pathogen. The scheme was initially used to study a collection of 160 bovine milk isolates from the United Kingdom and showed that the majority of isolates were from one clonal complex (designated the ST-5 complex). Here we describe the MLST analysis of a collection of New Zealand isolates. These were obtained from diverse sources, including bovine milk, other bovine anatomical sites, and environmental sources. The complete allelic profiles of 253 isolates were determined. The collection was highly diverse and included 131 different sequence types (STs). The New Zealand and United Kingdom populations were distinct, since none of the 131 STs were represented within the previously studied collection of 160 United Kingdom S. uberis isolates. However, seven of the STs were members of the ST-5 clonal complex, the major complex within the United Kingdom collection. Two new clonal complexes were identified: ST-143 and ST-86. All three major complexes were isolated from milk, other bovine sites, and the environment. Carriage of the hasA gene, which is necessary for capsule formation, correlated with clonal complex and isolation from clinical cases of mastitis.     

Application of Streptococcus uberis Multilocus Sequence Typing: Analysis of the Population Structure Detected among Environmental and Bovine Isolates from New Zealand and the United Kingdom

We recently developed a multilocus sequence typing (MLST) scheme to differentiate S. uberis isolates and facilitate an understanding of the population biology of this pathogen. The scheme was initially used to study a collection of 160 bovine milk isolates from the United Kingdom and showed that the majority of isolates were from one clonal complex (designated the ST-5 complex). Here we describe the MLST analysis of a collection of New Zealand isolates. These were obtained from diverse sources, including bovine milk, other bovine anatomical sites, and environmental sources. The complete allelic profiles of 253 isolates were determined. The collection was highly diverse and included 131 different sequence types (STs). The New Zealand and United Kingdom populations were distinct, since none of the 131 STs were represented within the previously studied collection of 160 United Kingdom S. uberis isolates. However, seven of the STs were members of the ST-5 clonal complex, the major complex within the United Kingdom collection. Two new clonal complexes were identified: ST-143 and ST-86. All three major complexes were isolated from milk, other bovine sites, and the environment. Carriage of the hasA gene, which is necessary for capsule formation, correlated with clonal complex and isolation from clinical cases of mastitis.     

Genetic diversity of Campylobacter jejuni isolates from farm animals and the farm environment

The genetic diversity of Campylobacter jejuni isolates from farm animals and their environment was investigated by multilocus sequence typing (MLST). A total of 30 genotypes, defined by allelic profiles (assigned to sequence types [STs]), were found in 112 C. jejuni isolates originating in poultry, cattle, sheep, starlings, and slurry. All but two of these genotypes belonged to one of nine C. jejuni clonal complexes previously identified in isolates from human disease and retail food samples and one clonal complex previously associated with an environmental source. There was some evidence for the association of certain clonal complexes with particular farm animals: isolates belonging to the ST-45 complex predominated among poultry isolates but were absent among sheep isolates, while isolates belonging to the ST-61 and ST-42 complexes were predominant among sheep isolates but were absent from the poultry isolates. In contrast, ST-21 complex isolates were distributed among the different isolation sources. Comparison with MLST data from 91 human disease isolates showed small but significant genetic differentiation between the farm and human isolates; however, representatives of six clonal complexes were found in both samples. These data demonstrate that MLST and the clonal complex model can be used to identify and compare the genotypes of C. jejuni isolates from farm animals and the environment with those from retail food and human disease.     

Genetic Diversity of Campylobacter jejuni Isolates from Farm Animals and the Farm Environment

The genetic diversity of Campylobacter jejuni isolates from farm animals and their environment was investigated by multilocus sequence typing (MLST). A total of 30 genotypes, defined by allelic profiles (assigned to sequence types [STs]), were found in 112 C. jejuni isolates originating in poultry, cattle, sheep, starlings, and slurry. All but two of these genotypes belonged to one of nine C. jejuni clonal complexes previously identified in isolates from human disease and retail food samples and one clonal complex previously associated with an environmental source. There was some evidence for the association of certain clonal complexes with particular farm animals: isolates belonging to the ST-45 complex predominated among poultry isolates but were absent among sheep isolates, while isolates belonging to the ST-61 and ST-42 complexes were predominant among sheep isolates but were absent from the poultry isolates. In contrast, ST-21 complex isolates were distributed among the different isolation sources. Comparison with MLST data from 91 human disease isolates showed small but significant genetic differentiation between the farm and human isolates; however, representatives of six clonal complexes were found in both samples. These data demonstrate that MLST and the clonal complex model can be used to identify and compare the genotypes of C. jejuni isolates from farm animals and the environment with those from retail food and human disease.     

Molecular phylogenetics of Candida albicans

We analyzed data on multilocus sequence typing (MLST), ABC typing, mating type-like locus (MAT) status, and antifungal susceptibility for a panel of 1,391 Candida albicans isolates. Almost all (96.7%) of the isolates could be assigned by MLST to one of 17 clades. eBURST analysis revealed 53 clonal clusters. Diploid sequence type 69 was the most common MLST strain type and the founder of the largest clonal cluster, and examples were found among isolates from all parts of the world. ABC types and geographical origins showed statistically significant variations among clades by univariate analysis of variance, but anatomical source and antifungal susceptibility data were not significantly associated. A separate analysis limited to European isolates, thereby minimizing geographical effects, showed significant differences in the proportions of isolates from blood, commensal carriage, and superficial infections among the five most populous clades. The proportion of isolates with low antifungal susceptibility was highest for MAT homozygous a/a types and then alpha/alpha types and was lowest for heterozygous a/alpha types. The tree of clades defined by MLST was not congruent with trees generated from the individual gene fragments sequenced, implying a separate evolutionary history for each fragment. Analysis of nucleic acid variation among loci and within loci supported recombination. Computational haplotype analysis showed a high frequency of recombination events, suggesting that isolates had mixed evolutionary histories resembling those of a sexually reproducing species.     

Identification and characterization of insect-specific proteins by genome data analysis

Background

Previous studies have sought to identify a link between the distribution of variable genes amongst isolates of Campylobacter jejuni and particular host preferences. The genomic sequence data available currently was obtained using only isolates from human or chicken hosts. In order to identify variable genes present in isolates from alternative host species, five subtractions between C. jejuni isolates from different sources (rabbit, cattle, wild bird) were carried out, designed to assess genomic variability within and between common multilocus sequence type (MLST) clonal complexes (ST-21, ST-42, ST-45 and ST-61).

Results

The vast majority (97%) of the 195 subtracted sequences identified had a best BLASTX match with a Campylobacter protein. However, there was considerable variation within and between the four clonal complexes included in the subtractions. The distributions of eight variable sequences, including four with putative roles in the use of alternative terminal electron acceptors, amongst a panel of C. jejuni isolates representing diverse sources and STs, were determined.

Conclusion

There was a clear correlation between clonal complex and the distribution of the metabolic genes. In contrast, there was no evidence to support the hypothesis that the distribution of such genes may be related to host preference. The other variable genes studied were also generally distributed according to MLST type. Thus, we found little evidence for widespread horizontal gene transfer between clonal complexes involving these genes.     

Identification of the clonal complexes of Staphylococcus aureus strains by determination of the conservation patterns of small genomic islets

Aims:  To investigate the clonality of Staphylococcus aureus isolates, it is important to identify their clonal complexes (CCs) with multilocus sequence typing (MLST). However, it is expensive to carry out MLST analyses for many isolates. The aim of this study, therefore, was to develop a cost-effective method to identify CCs by determining the conservation pattern of 'small genomic islets' (SGIs). SGIs are nonconserved regions between strains and have single or multiple open-reading frames (ORFs).
Methods and Results:  The whole-genome sequences of nine strains were compared in order to select 16 SGIs. The conservation patterns of the 16 SGIs (islet patterns) were investigated in 136 S. aureus isolates, which were classified into 21 CCs. The islet patterns (IPs) exhibited a one-to-one correspondence with the CCs, except for isolates belonging to CC1, CC5 and CC8. The IPs typical of strains belonging to CC1, CC5 and CC8 differed between those of sequence type 1 (ST1) and ST188 (CC1), ST5 and ST6 (CC5) and ST8 and ST239 (CC8).
Significance and Impact of the Study:  The CCs of many isolates can be identified in an easy and inexpensive manner by detecting these 16 SGIs. Emergent clones, particularly methicillin-resistant ones, can be identified by examining numerous islets by IP analysis.     

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.     

Molecular epidemiology and population structure of the honey bee brood pathogen Melissococcus plutonius

Melissococcus plutonius is the causative agent of European foulbrood (EFB), which is a serious brood disease of the European honey bee (Apis mellifera). EFB remains a threat because of a poor understanding of disease epidemiology. We used a recently published multi-locus sequence typing method to characterise 206 M. plutonius isolates recovered from outbreaks in England and Wales over the course of 2 years. We detected 15 different sequence types (STs), which were resolved by eBURST and phylogenetic analysis into three clonal complexes (CCs) 3, 12 and 13. Single and double locus variants within CC3 were the most abundant and widespread genotypes, accounting for 85% of the cases. In contrast, CCs 12 and 13 were rarer and predominantly found in geographical regions of high sampling intensity, consistent with a more recent introduction and localised spread. K-function analysis and interpoint distance tests revealed significant geographical clustering in five common STs, but pointed to different dispersal patterns between STs. We noted that CCs appeared to vary in pathogenicity and that infection caused by the more pathogenic variants is more likely to lead to honey bee colony destruction, as opposed to treatment. The importance of these findings for improving our understanding of disease aetiology and control are discussed.     

How clonal is Staphylococcus aureus?

Staphylococcus aureus is an important human pathogen and represents a growing public health burden owing to the emergence and spread of antibiotic-resistant clones, particularly within the hospital environment. Despite this, basic questions about the evolution and population biology of the species, particularly with regard to the extent and impact of homologous recombination, remain unanswered. We address these issues through an analysis of sequence data obtained from the characterization by multilocus sequence typing (MLST) of 334 isolates of S. aureus, recovered from a well-defined population, over a limited time span. We find no significant differences in the distribution of multilocus genotypes between strains isolated from carriers and those from patients with invasive disease; there is, therefore, no evidence from MLST data, which index variation within the stable "core" genome, for the existence of hypervirulent clones of this pathogen. Examination of the sequence changes at MLST loci during clonal diversification shows that point mutations give rise to new alleles at least 15-fold more frequently than does recombination. This contrasts with the naturally transformable species Neisseria meningitidis and Streptococcus pneumoniae, in which alleles change between 5- and 10-fold more frequently by recombination than by mutation. However, phylogenetic analysis suggests that homologous recombination does contribute toward the evolution of this species over the long term. Finally, we note a striking excess of nonsynonymous substitutions in comparisons between isolates belonging to the same clonal complex compared to isolates belonging to different clonal complexes, suggesting that the removal of deleterious mutations by purifying selection may be relatively slow.