Skip to main content

Comparative genomic assessment of Multi-Locus Sequence Typing: rapid accumulation of genomic heterogeneity among clonal isolates of Campylobacter jejuni



Multi-Locus Sequence Typing (MLST) has emerged as a leading molecular typing method owing to its high ability to discriminate among bacterial isolates, the relative ease with which data acquisition and analysis can be standardized, and the high portability of the resulting sequence data. While MLST has been successfully applied to the study of the population structure for a number of different bacterial species, it has also provided compelling evidence for high rates of recombination in some species. We have analyzed a set of Campylobacter jejuni strains using MLST and Comparative Genomic Hybridization (CGH) on a full-genome microarray in order to determine whether recombination and high levels of genomic mosaicism adversely affect the inference of strain relationships based on the analysis of a restricted number of genetic loci.


Our results indicate that, in general, there is significant concordance between strain relationships established by MLST and those based on shared gene content as established by CGH. While MLST has significant predictive power with respect to overall genome similarity of isolates, we also found evidence for significant differences in genomic content among strains that would otherwise appear to be highly related based on their MLST profiles.


The extensive genomic mosaicism between closely related strains has important implications in the context of establishing strain to strain relationships because it suggests that the exact gene content of strains, and by extension their phenotype, is less likely to be "predicted" based on a small number of typing loci. This in turn suggests that a greater emphasis should be placed on analyzing genes of clinical interest as we forge ahead with the next generation of molecular typing methods.


Campylobacter jejuni is the most common cause of acute bacterial enteritis worldwide [1, 2]. Despite significant progress in recent years, critical gaps remain in our understanding of C. jejuni pathogenesis. The lack of a well-defined set of virulence determinants makes it difficult to assess the virulence potential of different strains or to make links between specific genotypes and specific disease manifestations. Similarly, because the majority of infections are sporadic, sources and routes of transmission remain unclear in most cases of campylobacteriosis [3].

Significant effort has been placed on the development of methods for the typing of C. jejuni based on the analysis of polymorphic DNA targets and that have been applied to the study of species diversity and in the context of epidemiology and surveillance [4, 5]. The large number of competing approaches is a reflection on the fact that different methods may be appropriate for investigating short-term outbreak investigations (i.e. local epidemiology) and/or large-scale longitudinal surveillance (i.e. global epidemiology) [6]. Multi-locus sequence typing or MLST [7, 8], which is based on the analysis of DNA sequence polymorphisms in a group of housekeeping genes, has recently emerged as a strong contender for a genotyping "gold standard" for C. jejuni on the strength of several features. These include: a high discriminative power; ease of standardization of data acquisition and analysis across laboratories, and the high portability of the resulting sequence data [4, 9].

MLST benefits from a well-established framework for the phylogenetic analysis of molecular sequences. This has led to the suggestion that the phylogenetic signal contained within the loci analysed by MLST could be successfully used for long-term tracking in population structure studies, global epidemiology and long-term surveillance [10]. However, two outstanding questions need to be addressed in light of emerging data from comparative genomic analyses of C. jejuni. First, C. jejuni is naturally transformable and takes up homologous DNA readily [11], leading to high rates of intraspecies recombination [12, 13] that could distort the genetic relationships inferred from any one genetic locus. Second, a potential weakness which MLST shares with most genotyping approaches is that strain relatedness is inferred based on a very limited sub-sampling of the entire genome [5, 14]. This becomes increasingly relevant given the extensive genomic diversity that has been observed in intraspecies comparisons of C. jejuni through whole genome sequencing [15] and whole-genome microarray-based comparative genomic hybridization (CGH) [1620].

MCGH has recently been successfully applied to the examination of gene conservation dynamics and to the investigation of strain to strain relationships based on whole-genome gene conservation profiles [21]. In this study, we have analyzed a set of strains using both MLST and MCGH in order to assess whether the strain relationships inferred from the seven loci interrogated by MLST are consistent with the phylogenetic signal obtained from the analysis of whole-genome comparative genomic data.


Description of isolate relationships determined by MLST

In order to evaluate relationships among isolates, all 45 strains in this study were analyzed by MLST (Figure 1). The strains were selected from a larger dataset analyzed by MLST and were picked to comprise several levels of genetic similarity, thus enabling us to determine whether relationships assessed by MLST would be supported by CHG data in the short-term vs. long-term epidemiological context.

Figure 1

UPGMA-based clustering of MLST data for the 45 C. jejuni strains included in this study. Clusters representing clonal complexes (CC) are highlighted in red on the dendogram and their corresponding allelic profiles are also boxed in red. Allelic differences with respect to the central sequence type (ST) of the CC are highlighted in blue.

The dataset contained representatives from 25 distinct Sequence Types (STs) with 8 STs containing multiple strains. BURST analysis identified two main lineages, clonal complexes ST-21 and ST-45, which figured prominently in the dataset. The strains from the ST-21 complex include strains of the ST-21, Single Locus Variants (SLVs) (ST-50, ST-53, and ST-262), and Double Locus Variants (DLVs) (ST-43 and ST-184). The ST-45 complex includes strains of the ST-45 and SLVs (ST-25 and ST-241).

A subset of strains with the ST-474, which was originally placed in the ST-21 complex based on preliminary BURST analysis on the strength of matches at 5 of the 7 typing loci, was subsequently re-assigned to the ST-48 complex based on lineage assignments obtained from the C. jejuni MLST database, which is based on a much larger dataset including data on over 3000 strains. The ST-48 complex includes strains with the ST-48, SLVs (ST-474) and a Triple Locus Variant (TLV) (ST-475). In addition to the strains from ST-21, ST-45, and ST-48 complexes, groups of strains with identical STs were found from three additional clonal complexes (ST-353, ST-354, and ST-403). The remaining 10 strains did not belong to any of the ST complexes represented in the dataset and shared at most 3 of 7 MLST alleles with their closest matches.

Analysis of isolate relationships determined by MCGH

In order to assess strain-to-strain relationships based on genome similarity, gene conservation profiles derived from CGH data were used to quantify genomic similarity and this data was then used as a measure of the genetic distance between strains. Hierarchical clustering of the strains was performed on the resulting distance matrix of all pair wise distances between strains and bootstrap analysis revealed seven statistically robust clusters of strains (Figure 2) which were highly concordant with those obtained by hierarchical clustering of microarray profiles using the Pearson correlation metric (results not shown).

Figure 2

Concordance between clustering of MLST and whole-genome CGH profiles for the 45 C. jejuni isolates included in this study. Bootstrap support is shown for the statistically robust clusters (shown in red on the dendogram; CGH profiles boxed in blue). Log Ratio data has been colour-coded according to data interpretation thresholds described in Taboada et al. [31]. Strains showing discordant clustering results are boxed in green.

Cluster I can be further divided into 5 sub-types with significant differences in gene content: clusters Ia and Ib include the reference genome strain NCTC 11168 and six other strains with a very small number of genes displaying "significant Log Ratio differences" (SLRDs) with respect to NCTC 11168. These SLRDs correspond to likely gene divergence/gene absence events with respect to the reference strain. The three additional strains in cluster Ia have an average of 12.0 SLRDs whereas strains in cluster Ib have an average of 14.3 SLRDs. Strains in clusters Ic, Id, and Ie have an increasing number of differences with respect to the reference strain (an average of 37.3 SLRDs), although the majority of these are concentrated within four genetic loci: a region spanning Cj0968 to Cj0972 in cluster Ie; the lipo-oligosaccharide biosynthesis locus or LOS (Cj1136-Cj 1146c) in clusters Id and Ie; the capsular polysaccharide biosynthesis locus or CPS (Cj1414c-Cj1449) in clusters Ic, Id, and Ie; and a Type I restriction/modification locus or R-M (Cj1549-Cj1560) in cluster Ie.

Strains from clusters II to VII have an average number of SLRDs with respect to the reference strain that range from 64.5 for cluster IV to 97.9 for cluster V. The distribution and prevalence of various SLRDs across the genome varies substantially for each group, although the bulk can be found within the various hyper-variable loci previously described in C. jejuni [18, 19]. For example: only strains from clusters I and II appear to have a fully conserved (i.e. lacking in SLRDs) region spanning Cj0480 to Cj0490; only strains from cluster V have SLRDs in the region spanning Cj0727 and Cj0741; only strains from clusters I, III and IV lack SLRDs within the Type I restriction/modification locus. It thus appears that the various clusters of strains are characterized by unique patterns of conserved genes at hyper-variable loci.

Comparison of MCGH versus MLST-derived isolate relationships

When results of MCGH-based clustering were compared to the results obtained by MLST-based clustering and BURST analysis (Figure 2) we found that, with a few exceptions, the statistically robust groups obtained from CGH analysis correspond to groups of strains of identical sequence type. For example, strains in cluster Ia are largely of ST-21; strains from cluster 1b are largely of ST-53; strains from cluster 1d are of ST-262; strains from cluster II are of ST-354; strains from cluster IV are of ST-5; strains from cluster V are of ST-45; strains from cluster VII are of ST-403. In general, the congruence observed between CGH and MLST profiles also extends to strains within the same clonal complex (i.e. defined by sharing at least 4 loci) since strains with similar CGH profiles tend to share multiple MLST alleles. For example, strains in clusters Ia, Ib, Ic, and Id form part of clonal complex ST-21 and share 5 or more alleles. Similarly, the eight strains in cluster V share 6 or more alleles.

An inspection of local gene conservation patterns shows that strains from the same ST tend to share attributes that are nearly exclusive to the group (Figure 3). For example, whereas the ten strains with ST-45 have SLRDs at Cj0057 and Cj0058, this specific pattern is observed in only 5 of the remaining 35 strains in the dataset. Similarly, the pattern found among ST-45 strains in the region from Cj0177 to Cj0181 is found in only 7 of the remaining strains in the dataset. Other genomic regions where SLRDs are found almost exclusively among ST-45 strains are the multi-gene loci spanning Cj0296c to Cj0299, Cj0617-Cj0618, Cj0727 to Cj0733, and eight other single-gene loci (Cj0030, Cj0380c, Cj0690c, Cj0753c, Cj0794, Cj1305, Cj1585c, Cj1668c). A small number of genes (Cj0246c, Cj0859c, Cj0860, and Cj0970) appear to be fully conserved in all ST-45 strains but have SLRDs in a number of strains in the remainder of the dataset. More broadly, such differences in gene content can be used to differentiate the various CGH clusters and ST complexes.

Figure 3

Shared genomic attributes in strains from the same MLST clonal complex. The strains of ST-45 show significant differences in gene conservation rates at the loci shown in the first column with respect to all other strains in the dataset and differentiate this group of genetically related strains from other groups of strains.

Although our data suggests that strains with similar MLST profiles share similar CGH profiles, we have also found evidence for strains with significant levels of genomic similarity despite sharing few or none of the alleles used for MLST. For example, strain GBS_PG836 shows high overall genome similarity to the strains in the ST-48 complex (Cluster Ie) despite having different alleles at 4 of 7 MLST loci. In an extreme example, strains NCTC_13260 and NCTC_13261 show significant congruence in overall CGH profiles despite sharing no mutual alleles at any of the seven MLST loci. Thus, despite the fact that similarity in MLST profiles is generally a good predictor for genomic similarity, it is not always indicative of overall genome similarity between strains.

Genomic heterogeneity within groups of strains with the same MLST sequence type

Although genetic relatedness is reflected by shared gene content and high similarity in overall CGH profiles, when global and local gene conservation profiles are examined it is also possible to observe significant differences in gene content between strains of the same CGH cluster/clonal complex, particularly within hyper-variable genomic loci. For example, an examination of the strains in the ST-45 complex reveals that each strain appears to contain a heterogeneous mixture of conserved and divergent/absent genes and thus a "mosaic" pattern of gene conservation is apparent even among strains with high overall genomic similarity that are members of the same MLST clonal complex (Figure 4a, 4b). Although, on average, strains from the ST-45 complex have more SLRDs with respect to other strains in the dataset (μ = 93 ± 13), they bear a significant number of SLRDs with respect to one another (μ = 65 ± 14). Similar observations can be made for the strains in the ST-354 complex (μ = 83 ± 15 vs. μ = 35 ± 15). [For a complete set of all pair-wise SLRDs consult Additional file 1]

Figure 4

An examination of genomic mosaicism within clonal complex ST-45. Although strains within ST-45 have similar overall CGH profiles (A), significant genomic heterogeneity can be observed across various hyper-variable loci in the C. jejuni genome (B). Mosaicism observed in the CGH data is consistent with that observed in newly sequenced C. jejuni genomes (C). (note: Log Ratio data in (A) and (B) and sequence identity data in (C) were colour coded using a common scale reflecting the likelihood of gene presence/absence).

To examine whether these apparent mosaic patterns of gene content do not merely represent an artifact of the MCGH technique, we examined and visualized gene conservation patterns of various hyper-variable loci among newly sequenced C. jejuni genomes (Figure 4c). This examination revealed that the hypervariable loci in these strains have highly heterogeneous gene content with conserved and absent/divergent genes interspersed, a pattern that is consistent with our microarray-CGH data.

Disruption of genetic linkage in a genomic region flanked by two MLST loci

We have exploited the relatively close proximity between two MLST loci in combination with microarray-derived comparative genomic data to examine genetic linkage in C. jejuni. The tkt (Cj1645) and gltA (Cj1682c) loci are located approximately 36 Kb apart in the C. jejuni NCTC 11168 genome and variability has been observed in several genes contained within the region flanked by these two genes [19]. We thus set out to examine whether we could find evidence for an association between specific allele combinations of tkt and gltA and gene conservation profiles in the intervening region among members of the same clonal complex. Although it is possible to observe similar gene conservation profiles for members of the same clonal complex sharing the same tkt-gltA allele combinations (data not shown), some genomic heterogeneity is also apparent (Figure 5).

Figure 5

An examination of genetic linkage in groups of genetically related strains. Although strains with identical tkt-gltA alleles can also share similar gene conservation profiles within the intervening genomic region, disruption of linkage is apparent among members the same CC that share the same tkt-gltA loci.


Multi-locus Sequence Typing (MLST) has emerged as a leading molecular typing method for examining strain relationships which has been effectively applied to a number of different bacterial species [10], including C. jejuni [7, 22]. Although MLST has been successfully applied to the study of the population structure of the C. jejuni it has also, paradoxically, provided compelling evidence that C. jejuni populations are subject to high rates of horizontal genetic exchange, with recombinational events contributing to a significant proportion of the allelic diversity observed [12, 13].

The effects of high rates of intraspecies recombination observed include: a) conflicting phylogenetic signal obtained from different genes due to their different evolutionary trajectories, and b) a panmictic population structure for which clonal evolution is not the predominant trend. Both of these effects could pose limitations on the reliability of genetic relationships inferred from one or a small number of molecular markers as this small number of loci could themselves be subject to recombination. We have used microarray-based comparative genomic hybridization (MCGH) to analyze a collection of strains for which MLST analysis suggests varying levels of genetic relatedness. The dataset also includes clonal clusters comprised of strains with and without apparent epidemiological links in order to provide a "whole-genome" context for examining strain relationships inferred from MLST data.

An interesting finding from this study is the fact that global patterns of gene conservation obtained by MCGH are well correlated to MLST data. We have found that the robust clusters predicted by MCGH analysis (> 75% bootstrap support) and by MLST analysis (clonal complexes sharing 4 or more alleles) display excellent agreement (Figure 3). Of the 35 strains in our dataset that fall within six MLST clonal complexes only 4 fail to cluster robustly within a corresponding CGH cluster. Similarly, of the 37 strains that form statistically robust CGH clusters only 3 lack support from MLST data. Although the relationship between MLST genotypes and global-gene conservation profiles might be expected for strains with shared epidemiology, surprisingly this relationship appears to be evident for strains that do not share an obvious epidemiological connection. It is worth noting that the dendograms obtained for this dataset from the phylogenetic analysis of individual MLST loci are largely incongruent with one another (data not shown), likely due to recombination at the various alleles. Overall congruence between gene content and MLST has recently been observed for Streptococcus pneumoniae, another highly recombinogenic species [23]. In this light, the high degree of congruence between MLST and CGH data is not surprising and suggests that the multilocus approach appears to significantly mitigate the effects of lateral exchange in the examination of strain relationships.

It is important to stress that congruence between global gene conservation patterns and MLST genotypes does not preclude significant differences in gene content between related strains. Although, as expected, gene content differences tend to be highest between unrelated strains with greatest genetic distance (i.e. based on analysis of the MLST loci), our data provides evidence for significant genomic mosaicism between closely related strains through the accumulation of gene content differences. Thus, while related strains may share an increased number of genomic features, including a similar profile of significant SLRDs at any of a number of hyper-variable loci spread throughout the genome, their specific gene content in terms of absent and divergent genes may differ considerably. Although the widespread extent of this mosaic pattern of gene content might appear to be an artefact of the CGH technique we have employed in our analysis, evidence from comparative genomic sequencing would suggest otherwise. For example, mosaic patterns of gene conservation have been previously observed at a number of hyper-variable loci (LOS: [24, 25]; CPS: [26]; and RM: [27]). Similarly, our preliminary examination of additional hyper-variable loci among newly sequenced C. jejuni genomes (Figure 5c) has yielded similar observations. An examination of the genomic region bracketed by the MLST loci tkt and gltA also demonstrates that genomic events altering local gene conservation profiles can occur among members of the same clonal complex that share the same alleles at these loci (Figure 5), implying that great care must be made in extrapolating the gene content of strains based on indirect observations made at a different set of genetic loci.

A common theme among comparative sequencing studies is the suggestion that recombination is a potent driving force shaping the gene conservation patterns of hyper-variable loci of C. jejuni through events such as allelic replacements, gene fusions, gene duplications and gene deletions. Analysis of MLST allele patterns further suggests that housekeeping loci are also targeted by recombination [12, 13]. It thus appears that recombinational exchange in C. jejuni is not only widespread but that it must occur at significant frequencies consistent with the rapid accumulation of gene content differences we have observed among closely related strains in this study.


Although our data suggest that reliable strain relationships can be inferred despite the rapid pace of genetic change due to recombination, our ability to couple molecular typing data to phenotypes of interest (e.g. virulence, drug resistance) may be restricted by the shifting gene content among related strains. An advantage of molecular characterization methods based on the comparison of gene content is illustrated by our recent analysis of C. jejuni strains implicated in Guillain-Barré and Miller Fisher syndromes [28]. Neuropathogenic C. jejuni have been highly refractory to analysis by conventional molecular typing because of their diverse lineage and due to the lack of association between conventional molecular typing markers and their clinical phenotype [29]. Although our whole-genome CGH analysis merely confirmed earlier observations regarding the population structure of neuropathogenic strains, it correctly identified a small number of genes whose presence among strains of diverse lineage is though to be highly correlated with a neuropathic clinical outcome [28].

The rapid assessment of neuropathogenic potential of strains has since been achieved by the directly targeting polymorphisms within the genes of interest [30]. Ultimately, the development of clinically relevant molecular typing approaches may be better served by comparative genomic methods that directly survey the genetic differences responsible for the phenotypes of interest rather than through indirect evidence from comparison of molecular typing targets unrelated to phenotype


Bacterial strains

Background on the 45 strains we analyzed by microarray CGH and MLST is presented in Table 1. The strains were selected from a larger dataset analyzed by MLST and picked to comprise several levels of genetic similarity, which would thus enable us to determine whether relationships assessed by MLST would be supported by CHG data in the short-term vs. long-term epidemiological context. Strains were picked to comprise several levels of genetic similarity, which would thus enable us to determine whether relationships assessed by MLST would be supported by CHG data.

Table 1 List of strains used for this study.

DNA isolation

Cells were grown on Mueller-Hinton agar plates (BACTO, Oakville, ON) for 36 hours at 42°C under microaerophilic conditions prior to genomic DNA isolation. Genomic DNA was isolated by phenol:chloroform extraction as previously described [19]. For MLST, genomic DNA was prepared using the Qiagen Tissue Kit (Qiagen, Mississauga, ON) according to the manufacturer's instructions.

Microarray Hybridizations

Details of the microarray, including primer selection, the parameters for primer synthesis, selection of amplicons, as well as the purification and printing of DNA onto slides were previously described elsewhere [19]. Hybridizations were performed using protocols described previously [31]. Briefly, for each tester strain equivalent amounts of Cy-3 labelled tester and Cy-5 labelled control genomic DNAs (i.e. strain NCTC 11168) with similar dye incorporation efficiencies were pooled and co-hybridized to our microarray.

Microarray data acquisition and analysis

Microarrays were scanned using a Chipreader laser scanner (BioRad, Mississauga, ON) according to the manufacturer's recommendations. Spot quantification, visual inspection of potential outliers, and flagging of anomalous spots was performed using the program ArrayPro Analyzer (version 4.5; Media Cybernetics). The microarray data exported from ArrayPro was imported into the BioArray Software Environment (BASE version 1.2) [32] and is available at NCBI's Gene Expression Omnibus [33] under accession number GSE9919. Spots flagged due to poor spot morphology or low signal intensity (less than 3 X local background) were filtered out. After print-tip Loess normalization, data was used to calculate the average Log Ratio or LR (i.e. log2 [Signal Tester/Signal Control]) from the two replicates for each gene represented on the microarray. The filtered data exported from BASE contains the average LR data for 1606 genes.

MCGH data analysis and visualization

LR data was visualized and analyzed in TIGR's MultiExperiment Viewer (MEV version 3.0) [34] with high-resolution heat maps of LR data generated using a custom-script written in VBA for MS-Excel; all CGH data was organized assuming synteny with C. jejuni NCTC 11168 in order to examine mapping of variable genes to genomic regions. Clustering of strains based on LR profile similarities was performed by the average linkage hierarchical clustering method [35], as implemented in TMEV, using Pearson correlation coefficient as a distance metric with the Support Tree method of bootstrapping implemented in TMEV used to test the reliability of the clustering patterns (500 bootstrap re-samplings). A second method for clustering strains was developed based on calculating pair-wise genetic similarities in gene conservation profiles by using trinary thresholding of LR data [31], with a score of 1 given to all gene conservation matches (i.e. conserved, divergent of absent in both strains), a score of 0.5 given to absent/divergent pairs, and a score of 0 given to all other mismatches. The genetic similarity was then calculated by dividing the total score of all genes in the array by the total possible score. A custom VBA script for MS-Excel was written to calculate all pair-wise genetic distances (i.e. genetic distance = 1 – genetic similarity) and to calculate bootstrapped distance matrices which were used to create a consensus tree using the programs Neighbor and Consense from the phylogenetic inference package Phylip v.3.6 [36]. Phylogenetic trees were visualized using Treeview v.1.6.6 [37].

MLST analysis

MLST was performed using the methods of Dingle et al. [7, 22]. PCR amplification of the seven target genes was performed using primers described in the references above. Amplicons were purified with the Qiaquick column purification kit (Qiagen, Mississauga, ON) followed by Autoseq-96 (Molecular Dynamics). Sequencing was performed using the MegaBACE Long Read Matrix (Amersham Biosciences) according to the manufacturer's instructions and reaction products were separated on a MegaBACE 500 sequencer (GE Health Care, Piscataway, NJ). Sequence traces were analyzed using the MegaBACE software. MLST alleles and sequence types (ST) were determined for each strain by querying the C. jejuni MLST database at the University of Oxford [38] with the edited sequence data. Dendrograms based on the MLST sequence types were obtained using the method of unweighted pair group with arithmetic mean (UPGMA) implemented in the program START2 version 0.5.10 [39]. BURST analysis [40], as implemented in START2, was used to identify potential clonal complexes composed of strains sharing 5 or more identical alleles. Additional clonal complex assignments were determined by querying the C. jejuni MLST database with the allelic profiles of the strains.

Analysis of gene conservation patterns in newly sequenced C. jejunigenomes

Available sequence data from completed (strains NCTC 11168, RM1221, 81–176) and ongoing C. jejuni sequencing projects (strains 84–25, HB93-13, 260.94, 269.97, CF93-6) was obtained from NCBI's prokaryotic genome sequencing resource [41] and homology searching of genes in selected loci was performed using the program BLASTP [42]. Visualization of sequence identity levels was performed via heat maps generated using a custom-script written in VBA for MS-Excel.



comparative genomic hybridization


highly variable


highly divergent


moderately divergent


multi-locus sequence typing


sequence type


clonal complex


single-locus variant


double-locus variant


triple-locus variant


Log Ratio


significant Log Ratio difference


lipo-oligosaccharide biosynthesis locus


flagellar biosynthesis locus


capsular polysaccharide biosynthesis locus


type I restriction-modification locus.


  1. 1.

    Allos BM: Campylobacter jejuni Infections: update on emerging issues and trends. Clin Infect Dis. 2001, 32: 1201-1206.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Wassenaar TM, Blaser MJ: Pathophysiology of Campylobacter jejuni infections of humans. Microbes Infect. 1999, 1: 1023-1033.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Blaser MJ: Epidemiologic and clinical features of Campylobacter jejuni infections. J Infect Dis. 1997, 176 Suppl 2: S103-S105.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Maiden MC, Bygraves JA, Feil E, Morelli G, Russell JE, Urwin R, Zhang Q, Zhou J, Zurth K, Caugant DA, Feavers IM, Achtman M, Spratt BG: Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci U S A. 1998, 95: 3140-3145.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  5. 5.

    Wassenaar TM, Newell DG: Genotyping of Campylobacter spp. Appl Environ Microbiol. 2000, 66: 1-9.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  6. 6.

    van Belkum A, Struelens M, de Visser A, Verbrugh H, Tibayrenc M: Role of genomic typing in taxonomy, evolutionary genetics, and microbial epidemiology. Clin Microbiol Rev. 2001, 14: 547-560.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  7. 7.

    Dingle KE, Colles FM, Wareing DR, Ure R, Fox AJ, Bolton FE, Bootsma HJ, Willems RJ, Urwin R, Maiden MC: Multilocus sequence typing system for Campylobacter jejuni. J Clin Microbiol. 2001, 39: 14-23.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  8. 8.

    Sails AD, Swaminathan B, Fields PI: Utility of multilocus sequence typing as an epidemiological tool for investigation of outbreaks of gastroenteritis caused by Campylobacter jejuni. J Clin Microbiol. 2003, 41: 4733-4739.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  9. 9.

    Maiden MC: Multilocus sequence typing of bacteria. Annu Rev Microbiol. 2006, 60: 561-588.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Urwin R, Maiden MC: Multi-locus sequence typing: a tool for global epidemiology. Trends Microbiol. 2003, 11: 479-487.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Wang Y, Taylor DE: Natural transformation in Campylobacter species. J Bacteriol. 1990, 172: 949-955.

    PubMed Central  CAS  PubMed  Google Scholar 

  12. 12.

    Suerbaum S, Lohrengel M, Sonnevend A, Ruberg F, Kist M: Allelic diversity and recombination in Campylobacter jejuni. J Bacteriol. 2001, 183: 2553-2559.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  13. 13.

    Schouls LM, Reulen S, Duim B, Wagenaar JA, Willems RJ, Dingle KE, Colles FM, Van Embden JD: Comparative genotyping of Campylobacter jejuni by amplified fragment length polymorphism, multilocus sequence typing, and short repeat sequencing: strain diversity, host range, and recombination. J Clin Microbiol. 2003, 41: 15-26.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  14. 14.

    van Belkum A, Troesch A: High-density DNA arrays for comparative genomics and epidemiological studies in clinical microbiology. Expert Rev Mol Diagn. 2003, 3: 1-4.

    Article  PubMed  Google Scholar 

  15. 15.

    Fouts DE, Mongodin EF, Mandrell RE, Miller WG, Rasko DA, Ravel J, Brinkac LM, DeBoy RT, Parker CT, Daugherty SC, Dodson RJ, Durkin AS, Madupu R, Sullivan SA, Shetty JU, Ayodeji MA, Shvartsbeyn A, Schatz MC, Badger JH, Fraser CM, Nelson KE: Major structural differences and novel potential virulence mechanisms from the genomes of multiple campylobacter species. PLoS Biol. 2005, 3: e15-

    PubMed Central  Article  PubMed  Google Scholar 

  16. 16.

    Dorrell N, Mangan JA, Laing KG, Hinds J, Linton D, Al-Ghusein H, Barrell BG, Parkhill J, Stoker NG, Karlyshev AV, Butcher PD, Wren BW: Whole genome comparison of Campylobacter jejuni human isolates using a low-cost microarray reveals extensive genetic diversity. Genome Res. 2001, 11: 1706-1715.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  17. 17.

    Leonard EE, Takata T, Blaser MJ, Falkow S, Tompkins LS, Gaynor EC: Use of an open-reading frame-specific Campylobacter jejuni DNA microarray as a new genotyping tool for studying epidemiologically related isolates. J Infect Dis. 2003, 187: 691-694.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Pearson BM, Pin C, Wright J, I'Anson K, Humphrey T, Wells JM: Comparative genome analysis of Campylobacter jejuni using whole genome DNA microarrays. FEBS Lett. 2003, 554: 224-230.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Taboada EN, Acedillo RR, Carrillo CD, Findlay WA, Medeiros DT, Mykytczuk OL, Roberts MJ, Valencia CA, Farber JM, Nash JH: Large-scale comparative genomics meta-analysis of Campylobacter jejuni isolates reveals low level of genome plasticity. J Clin Microbiol. 2004, 42: 4566-4576.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  20. 20.

    Champion OL, Gaunt MW, Gundogdu O, Elmi A, Witney AA, Hinds J, Dorrell N, Wren BW: Comparative phylogenomics of the food-borne pathogen Campylobacter jejuni reveals genetic markers predictive of infection source. Proc Natl Acad Sci U S A. 2005, 102: 16043-16048.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  21. 21.

    Taboada EN, Luebbert CC, Nash JHE: Studying Bacterial Genome Dynamics Using Microarray-Based Comparative Genomic Hybridization. In Comparative Genomics (Vol.2) Volume 15. 396th edition. Edited by: Bergman NH. Humana Press; 2007:223-254.

    Google Scholar 

  22. 22.

    Dingle KE, Colles FM, Ure R, Wagenaar JA, Duim B, Bolton FJ, Fox AJ, Wareing DR, Maiden MC: Molecular characterization of Campylobacter jejuni clones: a basis for epidemiologic investigation. Emerg Infect Dis. 2002, 8: 949-955.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  23. 23.

    Dagerhamn J, Blomberg C, Browall S, Sjostrom K, Morfeldt E, Henriques-Normark B: Pattern of accessory genes predicts the same relatedness among strains of Streptococcus pneumoniae as sequencing house keeping genes: a novel approach in molecular epidemiology. J Clin Microbiol. 2007

    Google Scholar 

  24. 24.

    Gilbert M, Karwaski MF, Bernatchez S, Young NM, Taboada E, Michniewicz J, Cunningham AM, Wakarchuk WW: The genetic bases for the variation in the lipo-oligosaccharide of the mucosal pathogen, Campylobacter jejuni. Biosynthesis of sialylated ganglioside mimics in the core oligosaccharide. J Biol Chem. 2002, 277: 327-337.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Parker CT, Horn ST, Gilbert M, Miller WG, Woodward DL, Mandrell RE: Comparison of Campylobacter jejuni lipooligosaccharide biosynthesis loci from a variety of sources. J Clin Microbiol. 2005, 43: 2771-2781.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  26. 26.

    Karlyshev AV, Champion OL, Churcher C, Brisson JR, Jarrell HC, Gilbert M, Brochu D, St MF, Li J, Wakarchuk WW, Goodhead I, Sanders M, Stevens K, White B, Parkhill J, Wren BW, Szymanski CM: Analysis of Campylobacter jejuni capsular loci reveals multiple mechanisms for the generation of structural diversity and the ability to form complex heptoses. Mol Microbiol. 2005, 55: 90-103.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Miller WG, Pearson BM, Wells JM, Parker CT, Kapitonov VV, Mandrell RE: Diversity within the Campylobacter jejuni type I restriction-modification loci. Microbiology. 2005, 151: 337-351.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Taboada EN, van Belkum AF, Yuki N, Acedillo RR, Godschalk PCR, Koga M, Endtz HP, Gilbert M, Nash JH: Comparative genomic analysis of Campylobacter jejuni associated with Guillain-Barre and Miller Fisher syndromes: neuropathogenic and enteritis-associated isolates can share high levels of genomic similarity. BMC Genomics. 2007, 8: 359-

    PubMed Central  Article  PubMed  Google Scholar 

  29. 29.

    Dingle KE, van den Braak N, Colles FM, Price LJ, Woodward DL, Rodgers FG, Endtz HP, van Belkum A, Maiden MC: Sequence typing confirms that Campylobacter jejuni strains associated with Guillain-Barre and Miller-Fisher syndromes are of diverse genetic lineage, serotype, and flagella type. J Clin Microbiol. 2001, 39: 3346-3349.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  30. 30.

    Godschalk PC, van Belkum A, van den Braak N, van Netten D, Ang CW, Jacobs BC, Gilbert M, Endtz HP: PCR-Restriction Fragment Length Polymorphism Analysis of Campylobacter jejuni Genes Involved in Lipooligosaccharide Biosynthesis Identifies Putative Molecular Markers for Guillain-Barre Syndrome. J Clin Microbiol. 2007, 45: 2316-2320.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  31. 31.

    Taboada EN, Acedillo RR, Luebbert CC, Findlay WA, Nash JH: A new approach for the analysis of bacterial microarray-based Comparative Genomic Hybridization: insights from an empirical study. BMC Genomics. 2005, 6: 78-

    PubMed Central  Article  PubMed  Google Scholar 

  32. 32.

    Saal LH, Troein C, Vallon-Christersson J, Gruvberger S, Borg A, Peterson C: BioArray Software Environment (BASE): a platform for comprehensive management and analysis of microarray data. Genome Biol. 2002, 3: SOFTWARE0003-

    PubMed Central  Article  PubMed  Google Scholar 

  33. 33.

    NCBI's Gene Expression Omnibus . 2008, []

  34. 34.

    Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, Sturn A, Snuffin M, Rezantsev A, Popov D, Ryltsov A, Kostukovich E, Borisovsky I, Liu Z, Vinsavich A, Trush V, Quackenbush J: TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003, 34: 374-378.

    CAS  PubMed  Google Scholar 

  35. 35.

    Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998, 95: 14863-14868.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  36. 36.

    Felsenstein J: PHYLIP - Phylogeny Inference Package (Version 3.2). Cladistics. 1989, 5: 164-166.

    Google Scholar 

  37. 37.

    Page RD: TreeView: an application to display phylogenetic trees on personal computers. Comput Appl Biosci. 1996, 12: 357-358.

    CAS  PubMed  Google Scholar 

  38. 38.

    C. jejuni MLST database at the University of Oxford . 2008, []

  39. 39.

    Jolley KA, Feil EJ, Chan MS, Maiden MC: Sequence type analysis and recombinational tests (START). Bioinformatics. 2001, 17: 1230-1231.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Feil EJ, Li BC, Aanensen DM, Hanage WP, Spratt BG: eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data. J Bacteriol. 2004, 186: 1518-1530.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  41. 41.

    NCBI's prokaryotic genome sequencing resource . 2008, []

  42. 42.

    Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  43. 43.

    Pope C, Wilson J, Taboada EN, Mackinnon J, Felipe Alves CA, Nash JH, Rahn K, Tannock GW: Epidemiology, relative invasive ability, molecular characterization, and competitive performance of Campylobacter jejuni strains in the chicken gut. Appl Environ Microbiol. 2007, 73: 7959-7966.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

Download references


The authors would like to thank M. Karmali, G. Tannock, J. Farber, and the Campy-On-Ice Consortium for bacterial strains. The authors would also like to thank Andrew Kropinski for helpful discussions and critical reading of this manuscript. Funding for this work has been provided through the National Research Council's Genomics and Health Initiative and through grants from Health Canada's Office of Biotechnology and Science.

Author information



Corresponding author

Correspondence to Eduardo N Taboada.

Additional information

Authors' contributions

ENT contributed to study design, designed MCGH experiments, wrote custom Excel VBA scripts, carried out downstream data analysis and visualization, and drafted the manuscript; CCL performed all MCGH experiments and performed preliminary data analysis. JMM performed MLST analysis and assisted with preliminary MCGH analysis. VPJG assisted in data interpretation and drafting of the manuscript. JHEN and KR conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors submitted comments on drafts and read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Pairwise matrix of SLRDs for the 45 strains in the dataset. The (partial) gene content of strains was assessed by analyzing the CGH data using two sets of empirically determined thresholds [31]: Figure S1: < 1% error rate on "likely conserved" and "likely absent" calls. ; Figure S2: < 1% error rate on "likely conserved" and "likely divergent/absent". Strains have been arranged based on the UPGMA dendogram from analysis of MLST data. Boxes in bold represent pairwise distances between members of the same clonal complex (red branches on dendograms). (XLS 79 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Taboada, E.N., MacKinnon, J.M., Luebbert, C.C. et al. Comparative genomic assessment of Multi-Locus Sequence Typing: rapid accumulation of genomic heterogeneity among clonal isolates of Campylobacter jejuni. BMC Evol Biol 8, 229 (2008).

Download citation


  • Comparative Genomic Hybridization
  • Clonal Complex
  • Gene Conservation
  • Miller Fisher Syndrome
  • Comparative Genomic Hybridization Data