- Research article
- Open Access
Use of suppression subtractive hybridisation to extend our knowledge of genome diversity in Campylobacter jejuni
BMC Genomicsvolume 8, Article number: 110 (2007)
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).
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.
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.
Infection due to Campylobacter sp. is one of the major causes of diarrheal disease worldwide and is the most common source of bacterial gastroenteritis . Although transmission of Campylobacter occurs mainly through the consumption of livestock, with poultry being the most common source, C. jejuni has been isolated from diverse animal, human and environmental sources. Several recent studies have sought to determine clonal prevalences amongst isolates from these diverse sources by applying multilocus sequence typing (MLST) [2–7]. Whilst some MLST clonal complexes, such as the ST-21 complex, are widespread, others, such as the ST-61 and ST-42 complexes, have a more restricted distribution amongst different host animals, including humans [4, 7].
The complete genome sequences have been published for the C. jejuni strains NCTC11168 , RM1221  and C. jejuni strain 81–176 , including its pTet and pVir plasmids . Genome sequences have also been published for strains of C. lari, C. coli and C. upsaliensis  and further genome sequence projects for several other Campylobacter strains are ongoing . Although unpublished, the genome of C. fetus is also complete. In addition, strain-specific DNA sequences have been identified for C. jejuni ATCC43431 .
There has been considerable interest in characterising genetic variation between isolates of C. jejuni with a view to identifying those genes relevant to severity of disease or host colonisation. Inter-strain variations in loci such as those encoding lipooligosaccharide (LOS) , capsule  or restriction-modification (RM) systems  have already been characterised. Comparative genome analyses using microarrays, based largely upon the strain NCTC11168 genome sequence [16–19], indicate high levels of genome diversity but low levels of genome plasticity in C. jejuni . Recently it has been suggested that this kind of approach can help to identify genetic markers predictive of the source of an infection .
Since these studies were mostly restricted to those genes present in a single strain, NCTC11168, further identification of such markers would be greatly facilitated by the construction of microarrays containing all potentially variable genes. Genes contributing to plasticity amongst C. jejuni populations have been identified from the genomes of strains such as RM1221, 81–176 and ATCC43431 [9, 11, 13, 16]. However, it is not clear whether the widely accessible nucleotide and protein sequence databases are representative of the variable genes that occur in C. jejuni.
Suppression subtractive hybridisation (SSH) is a method designed to identify sequences present in one strain (the TESTER) but absent from a reference strain (the DRIVER) [22, 23] and has been applied previously to identify genetic differences between two C. jejuni human isolates with different colonisation potentials . In this study we used several rounds of SSH between and within C. jejuni MLST clonal complexes, using isolates from various sources, in order to gain a better understanding of the genomic variability that remains uncharacterised in populations of C. jejuni. Furthermore, we tested the hypothesis that the distribution of variable genes involved in the use of alternative terminal electron acceptors may be related to host preference.
Results and discussion
Our initial choice of restriction enzymes for the digestion of DNAs in the SSH was based on a previous study of C. jejuni . A summary of all of the SSH sequences obtained is shown in Table 1. The SSH libraries obtained with Alu I/Dra I-digested DNAs were dominated by very short DNA sequences (92% were <400 bp). In order to obtain longer sequences we carried out two further subtractions using Rsa I-digested DNAs. This reduced the proportion of sequences <400 bp in length. The majority of subtracted sequences from strains 670, 504 and 1967 matched sequences previously found in C. jejuni strains RM1221 or 81–176. Thus, in an attempt to enrich for sequences not previously reported, the final two subtractions were carried out with strain 1967 (ST-42 complex) as an additional driver to strain NCTC11168 (ST-21 complex). This proved successful in that the majority of SSH sequences no longer matched strains RM1221 or 81–176. However, the overall proportion of SSH sequences matching previously reported Campylobacter sp. was not reduced significantly. The combination of enrichment and the choice of an enzyme giving longer SSH fragments resulted in fewer SSH sequences overall (Table 1).
In total, 195 subtracted sequences were obtained. Details of these are given in Additional File 1. These SSH sequences were all confirmed as absent from the driver strain C. jejuni NCTC11168, even though occasionally the best BLASTX match in the database was against this strain. Furthermore, the second driver strain, strain 1967, was PCR-negative for nine of 11 sequences tested from the two dual driver subtractions. Of the 195 SSH sequences obtained only two (1%) had a best BLASTX match with proteins from outside the genus Campylobacter. A further three SSH sequences (2%) had no significant BLASTX match. Thus, the vast majority of SSH sequences matched within the genus Campylobacter. However, it should be noted that eight SSH sequences (4%), although having a best BLASTX match with a Campylobacter sp. protein, shared less than 80% identity with the matching protein. Those SSH sequences matching outside the genus Campylobacter, matching Campylobacter sequences with <80% identity, or used in distribution analysis are shown in Table 2.
Both the mean and the median % GC content for the 195 subtracted sequences was 29%, with a range of 18–42%. The mean % GC contents varied little between the subtractions (28.9–30.8%). Thus, the average % GC contents for the subtracted sequences was only slightly below the values reported for the genomes of the C. jejuni strains NCTC11168 (30.6%)  or RM1221 (30.3%) .
Of the three subtractions conducted using Alu I/Dra I-digested DNAs and a single driver strain (NCTC11168), two were intra-clonal complex and one was between clonal complexes. The numbers of subtracted sequences obtained did not vary significantly between these subtractions (Table 1). However, there were some variations in the putative functions of subtracted sequences. Plasmid and bacteriophage-associated sequences accounted for 73% of those sequences differing between the two ST-21 complex strains 670 and NCTC11168, and 16% of those sequences differing between the two ST-21 complex strains 504 and NCTC11168, but only 2% of those sequences differing between the ST-42 complex strain 1967 and strain NCTC11168. The largest group of subtracted sequences for the inter-clonal complex subtraction were those associated with metabolism/biosynthesis (29%). The SSH data suggest that there are both intra- and inter-clonal complex variations in genes associated with LOS, capsule, flagella/motility, membrane/transport and metabolism (Table 1).
It has been demonstrated that SSH is an effective method for analysing genetic differences between related strains. In a previous study using SSH, Agron et al.  were able to detect most of the 7% of genomic differences between two closely related, fully sequenced strains of Helicobacter pylori. Unlike in our study, the authors used four parallel subtractions with different restriction enzymes. However, they further demonstrated that as tester-specific sequences became limiting the proportion of repeat fragments increased . Because of our use of single restriction digestion conditions for each subtraction, we cannot claim the extent of coverage achieved by Agron et al. , but the increased proportion of repeat sequences obtained with the dual driver subtractions is indicative of a reduction in the overall pool of tester-specific sequences in these subtractions. It is our belief that we have detected a significant proportion of the genetic variation between the strains and that our data indicate that there are very few genes present in C. jejuni strains with best BLASTX matches outside the genus rather than with already sequenced genes/proteins from within the genus.
SSH sequences matching outside the genus Campylobacter
The two sequences matching outside the genus Campylobacter were a putative transport protein (SSH629-D8) and a putative phosphodiesterase (SSH629-D9) (Table 2). These SSH sequences have been submitted to GenBank (EF076761 and EF076762). The predicted protein sequence derived from SSH629-D8 matched a transport protein from Rhodopseudomonas palustris, predicted to be a cytoplasmic component of an ABC-type sugar transport system. The gene encoding the phosphodiesterase enzyme from Delftia acidovorans matching SSH629-D9 has been cloned and expressed in Escherichia coli. It shares sequence similarity to cyclic AMP (cAMP) phosphodiesterase and cyclic nucleotide phosphodiesterases and exhibited activity on cAMP in vivo . Using dot-blot hybridization, we determined the distribution of these two SSH sequences amongst a panel of C. jejuni isolates representative of common MLST clonal complexes and various sources (Table 3). The sequences SSH629-D8 and SSH629-D9 were found only in the ST-45 complex (Tables 3 and 4). It should be noted that SSH629-D8 appears to be part of a pseudogene containing a frame-shift mutation. We used PCR amplification to amplify this region from four of the C. jejuni isolates testing positive for SSH629-D8, including strain 629, and confirmed that this mutation is genuinely present and not an artefact of the SSH procedure.
Genomic islands and mobile elements
The distributions of SSH sequences according to the function of the matching proteins varied considerably between strains, reflecting the divergence amongst those strains chosen (Table 1). The SSH data suggest that the ST-21 wild duck isolate, strain 670, contains at least parts of two of the mobile elements reported previously in C. jejuni RM1221 (CMLP1 and CJIE4) [9, 26]. CMLP1 is a Mu-like bacteriophage and wasn't identified in any of the other tester strains used in this study. CJIE4 carries sequences similar to a putative prophage encoded by C. lari (CLIE1) . Strain 670 alone also carried multiple sequences matching the related pTet or pCC31 plasmids reported previously in C. jejuni 81–176 and C. coli respectively [27, 28]. SSH data also indicated that the ST-21 cattle isolate strain 504 carries the prophage-like element CJIE4 (18 matching SSH sequences). However, none of the strain 504 SSH sequences matched plasmid sequences. In fact, other than those matching CJIE4, most of the SSH sequences from strain 504 matched hypothetical proteins of unknown function (Table 1). The SSH sequences from strains 1967, 629 and 961, each representing a different MLST type, contained no matches to the genomic islands of strain RM1221 and very few matches to either bacteriophage or plasmid-related sequences (Table 1).
Using dot-blot hybridization we determined the distribution of the sequence 670-D7 as an indicator of the prevalence of the pTet/pCC31 plasmid amongst our strain panel (Table 1). The sequence was found in three isolates from the ST-21 complex and four isolates from the ST-61 complex.
Plasticity regions and variable genes
Most of the subtractions identified variations in loci previously reported as variable amongst C. jejuni strains. It has been reported that variable genes (absent or highly divergent) in the Campylobacter genome map to discrete areas, termed variously plasticity regions, hypervariable regions or plasticity zones [16, 17, 21, 29]. Indeed the genotypes of C. jejuni strains have been shown to be inherently unstable, and recombination events occur between unrelated strains both in vitro and in vivo even in the absence of selective pressure . Pearson et al.  identified seven hypervariable plasticity regions (PR1-7) in the genome of C. jejuni NCTC11168. PR6 contains the capsule biosynthesis locus, flanked by capsule transport genes. Twelve SSH sequences in all matched capsule-related proteins. PR7 contains numerous putative outer membrane proteins but also carries a divergent gene encoding a probable flagellar hook protein . Other motility-related genes, including flagellin and flagellin-glycosylation genes, are carried by PR5, along with LOS biosynthesis genes. SSH sequences matching flagellar hook proteins (FlgE), flagellins, flagellin glycosylation proteins and LOS biosynthesis proteins were all identified amongst the SSH sequences in our study. PR4 includes glycosyltransferases and galactosyltransferases of unknown functions as well as genes involved in the sialylation of LOS. SSH sequences matching proteins with similar putative functions were found with both the strain 961 and strain 1967 subtractions. In addition, the 629 subtracted library included one sequence matching both the NCTC11168 gmhA gene (from PR4) and a polysaccharide-related protein (SSH629-21).
RM systems protect bacteria from foreign DNA and may impact on the transfer of genes responsible for virulence or host colonisation. It has also been demonstrated previously that RM genes vary in Campylobacter spp.  and RM genes were amongst those identified in a previous study using SSH . The first three subtractions in our study all identified RM-related SSH sequences. Ahmed et al.  also identified subtracted sequences relating to arsenite-metabolising genes as present in strain 81–116 but absent from strain NCTC11168, noting that phenylarsonic compounds have been used in poultry feed and may have contaminated agricultural lands on which poultry litter has been used as manure. We identified an SSH matching the arsenical resistance protein of RM1221 in the subtracted library of strain 629 (SSH629-23), a rabbit isolate. In addition, Ahmed et al.  identified a putative γ-glutamyl transpeptidase gene in their subtraction. Equivalent genes have been implicated in a role in colonization of the gastric mucosa by Helicobacter pylori [31, 32]. We also identified SSH sequences matching a γ-glutamyl transpeptidase protein in strain 629. It has been reported that C. jejuni 81–176 carries a putative gene encoding a serine protease belonging to the autotransporter family . The SSH indicates that this gene is also present in the cattle isolate 1967 (SSH1967-A5 and SSH1967-G6). C. jejuni 81–176 and isolates 1967 are both members of the ST-42 clonal complex. We determined the distribution of SSH1967-A5 amongst our strain panel. All ST-42 isolates were positive for SSH1967-A5. Of the other isolates, only one ST-45 cattle isolate contained SSH1967-A5 (Tables 3 and 4).
We also determined the distribution of the SSH sequence SSH961-A9, which matches an enzyme from C. jejuni subsp. doylei (Table 2). With the exception of three isolates, the distribution was also restricted to the clonal complex of the tester strain (ST-61 complex; Tables 3 and 4).
When oxygen levels are low C. jejuni has the capacity to utilise a wide range of electron acceptors, including fumarate, nitrate, nitrite, sulfite, trimethylamine-N-oxide (TMAO) and dimethyl sulfoxide (DMSO) [33–35]. In strain NCTC11168 PR1 contains genes encoding the transport apparatus for molybdenum, which has a putative role in the reduction of nitrate as an alternative terminal electron acceptor . Just upstream of PR1 lies the gene (Cj0264c) encoding the reductase responsible for reduction of TMAO and DMSO under oxygen limiting conditions . Whilst lacking the Cj0264c gene or its close homologues, C. jejuni 81–176 carries an alternative DMSO reductase gene cluster (dmsABC-torD), and an additional cytochrome C biogenenesis gene cluster (cytC locus; cju02-09) . Cytochrome C may be an important link between the menaquinine pool and alternative terminal electron acceptors such as DMSO and TMAO . It has been suggested that the capacity of Campylobacter isolates to utilise alternative electron acceptors may contribute to selective advantages in specific ecological niches , and that the presence of additional respiratory capabilities may contribute to the efficiency of colonisation of highly pathogenic strains such as C. jejuni 81–176 . Thus, it is possible that the presence or absence of particular genes or islands contributing to growth in microaerophilic environments may influence host preference. Hofreuter et al. (2006) recently provided some evidence for this when demonstrating that a dmsA mutant of C. jejuni 81–176 colonised mice less well than its wild-type equivalent in a mixed infection model.
We identified several strain-variable SSH sequences matching putative enzymes with roles in electron transport using alternative terminal electron acceptors. The SSH library constructed from strain 629 (rabbit isolate) included SSH sequences with best BLASTX matches against genes in the cytochrome C biogenesis cluster of C. jejuni 81–176 (including SSH629-C10). The SSH library constructed from strain 1967 (cattle isolate) included numerous SSH sequences (including SSH1967-D9) with a best BLASTX match against the dmsABC genes of C. jejuni 81–176. The SSH libraries constructed from strains 1967, 670 and strain 504 each contained SSH sequences (SSH1967-C2, SSH670-B10 and SSH504-C10) matching a pyridine nucleotide disulfide oxidoreductase from C. jejuni 84-25 belonging to a family of enzymes that can play a role in electron transport. SSH sequences SSH670-B10 and SSH504-C10 were identical to each other.
In order to test the hypothesis that the distribution of such metabolic genes may be related to host preference, we chose to screen a panel of isolates varying in MLST type and source for the presence or absence of four representative sequences (Tables 3 and 4). SSH1967-H9 and SSH1967-C10 were chosen to represent the alternative DMSO reductase gene cluster (dmsABC-torD), and the additional cytochrome C biogenenesis gene cluster (cytC locus; cju02-09), of C. jejuni 81–176 respectively. SSH1967-C2 was chosen to represent the pyridine nucleotide disulfide oxidoreductase described in C. jejuni 84-25. In addition, the Cj0264c gene of strain NCTC11168, encoding a reductase, was included.
Interestingly, of the 81 C. jejuni isolates screened, only two (both ST-42) shared the sequences SSH1967-H9 and SSH1967-C2, whilst a further 13 lacked both sequences, including all of the ST-48 isolates. For most of the clonal complexes (ST-21, ST-61, ST-45 and ST-257) the sequences SSH1967-H9 and SSH1967-C2 were mutually exclusive.
There was a strong correlation between the distribution of the two reductase-related sequences (SSH1967-H9 and 1967-C2) or the cytochrome C biogenesis-related sequence (SSH629-C10) and the MLST clonal complex. It has been reported previously that the C. jejuni NCTC11168 ORF Cj0264c, which encodes the sole TMAO and DMSO reductase in this strain , was absent or highly divergent in 10 strains amongst a panel of 18 C. jejuini strains from diverse sources . Our data indicated that this gene was widespread amongst all clonal complexes in our panel of strains, with the exception of the ST-257 complex (Tables 3 and 4).
In our study, there was no evidence for an association between the presence of a particular gene associated with metabolism using alternative electron acceptors and the source of the isolate, suggesting that the presence or absence of these alternative metabolic genes does not play a significant role in niche preferences. However, it should be noted that our strain panel was dominated by isolates from cattle and human sources. Thus we cannot completely rule out such associations with other host species. In a study of gene expression variations by different variants of C. jejuni NCTC11168 it was reported that many of the differences in gene expression were in respiration and metabolism genes . The authors suggested that adaptation to different oxygen tensions may influence colonisation potential. The gene expression profiles compared were those of C. jejuni NCTC 11168-GS, the genome-sequenced isolate, and NCTC 11168-O, the original isolate from which NCTC11168-GS was derived. Of the two, isolate NCTC 11168-O is a much better coloniser of chicks and invades tissue culture cells far more efficiently . Under microaerobic and severely oxygen-limited conditions there were marked difference in the expression of genes associated with metabolism and respiration. Although ORF Cj0264c was not amongst those loci expressed differently between the two variants of NCTC 11168, such observations support the notion that rather than the presence or absence of metabolic and respiratory genes per se, variations in expression may be more relevant to niche preference.
In this study we have demonstrated that by broadening the range of clonal complexes and host sources of C. jejuni isolates submitted to genetic interrogation, we did not greatly increase the pool of identified strain-variable genes. This suggests that the current database already contains most of the diversity within this species. However, it should also be noted that genes currently associated with other closely related species, such as C. coli, can also make a contribution to the diversity within C. jejuni.
It is clear from our subtractions between ST-21 complex isolates that variation in gene content occurs within as well as between clonal complexes. However, our survey of strain-variable sequences, including four associated with genes involved in the use of alternative terminal electron acceptors, indicated a distribution according to clonal complex rather than host source. Thus, we found no evidence that the presence or absence of such genes plays a role in the host preferences of C. jejuni strains.
The bacterial strains used in this study are listed in Table 3 and were isolated in a previous study . Strains for SSH were chosen to enable comparisons within (ST-21 complex) and between MLST clonal complexes. In addition, isolates from diverse animal host sources were chosen (bird, cattle and rabbit; Table 1). These differed from the sources of previously sequenced strains of C. jejuni (human or chicken). A panel of isolates representing common MLST clonal complexes (ST-21, ST-42, ST-45, ST-48, ST-257 and ST-61) and different sources were used to study the distribution of subtracted sequences (Table 3). The bacteria were cultured on blood agar at 37°C under microaerophilic conditions.
MLST alleles, STs and clonal complexes were assigned using the Campylobacter PubMLST database  with sequences submitted for allele designation as appropriate.
Construction and screening of subtraction libraries
Genomic DNA for SSH was isolated from C. jejuni strains 670, 504, 1967, 629, 961 and NCTC11168 as described previously . SSH was carried out using the CLONTECH PCR-Select™ Bacterial Genome Subtraction Kit (Clontech) essentially as recommended by the supplier but with a hybridisation temperature of 58°C. In the first three hybridisations, DNAs from C. jejuni strains 670, 504 and 1967 respectively were used as tester and DNA from C. jejuni NCTC11168 was used as the driver. All DNAs were digested with Alu I and Dra I. A further two hybridisations, using Rsa I-digested DNA from C. jejuni strains 629 and 961 respectively as tester, were carried out with dual Rsa I-digested driver DNAs from strains 1967 and NCTC11168. PCR amplicons obtained following SSH were cloned into pGEM-T (Invitrogen). The subtraction libraries of Alu I/Dra I or Rsa I fragments thus constructed were screened by sequencing of plasmid DNA extracted from individual clones using M13 forward and reverse vector primers (Lark Technologies). BLAST searches at the C. jejuni NCTC11168 genome project web site  were used to determine the presence or absence of sequences in the NCTC11168 genome. Sequences absent from the genome of C. jejuni NCTC11168 were further analysed using BLASTN and BLASTX searches of the general database using the NCBI website  (last accessed 13th October 2006).
PCR amplification and dot-blot screening of strains
Oligonucleotide primers (Sigma-Genosys) for PCR amplifications are listed in Table 5 along with the annealing temperatures used. DNA for PCR amplifications and dot-blots was prepared by using the Wizard Genomic DNA Purification kit (Promega). For PCR amplification, typically, 1 μl of this DNA was used directly in 25 μl volumes containing 1.25 U of Taq DNA polymerase (Promega), 1 × TaqMaster (Helena Biosciences), 300 nM each primer, 1 × Taq buffer, 2.5 mM MgCl2 and 100 μM nucleotides (dATP, dCTP, dGTP, dTTP). Amplifications were carried out in an Eppendorf MasterCycler thermal cycler for 30 cycles consisting of 95°C (1 min), annealing temperature (1 min) and 72°C (2 min) with an additional extension time at 72°C (10 min) following completion of the 30 cycles.
Dot blot hybridisation of genomic DNA was carried out as described previously . Digoxigenin-11-2'-dUTP (DIG) (Roche)-labelled probes were made by carrying out PCR amplification in the presence of 60 μM DIG using vector or internal primers. Hybridisation and subsequent detection of DIG was carried out following the manufacturer's instructions (Roche).
All SSH sequence distributions were determined by dot-blot analysis, with the exception of 629-C10, for which a combination of PCR screening and Southern blots was used. This was because background hybridisation made the 629-C10 dot-blots difficult to interpret. A probe for flaA was used in dot blots to confirm the presence of DNA for each strain tested.
Allos BM: Campylobacter jejuni Infections: update on emerging issues and trends. Clin Infect Dis. 2001, 32: 1201-1206. 10.1086/319760.
Manning G, Dowson CG, Bagnall MC, Ahmed IH, West M, Newell DG: Multilocus sequence typing for comparison of veterinary and human isolates of Campylobacter jejuni. Appl Environ Microbiol. 2003, 69: 6370-6379. 10.1128/AEM.69.11.6370-6379.2003.
Sails AD, Swaminathan B, Fields PI: Clonal complexes of Campylobacter jejuni identified by multilocus sequence typing correlate with strain associations identified by multilocus enzyme electrophoresis. J Clin Microbiol. 2003, 41: 4058-4067. 10.1128/JCM.41.9.4058-4067.2003.
Colles FM, Jones K, Harding RM, Maiden MC: Genetic diversity of Campylobacter jejuni isolates from farm animals and the farm environment. Appl Environ Microbiol. 2003, 69: 7409-7413. 10.1128/AEM.69.12.7409-7413.2003.
Dingle KE, Colles FM, Falush D, Maiden MC: Sequence typing and comparison of population biology of Campylobacter coli and Campylobacter jejuni. J Clin Microbiol. 2005, 43: 340-347. 10.1128/JCM.43.1.340-347.2005.
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. 10.1128/JCM.39.1.14-23.2001.
French N, Barrigas M, Brown P, Ribiero P, Williams N, Leatherbarrow H, Birtles R, Bolton E, Fearnhead P, Fox A: Spatial epidemiology and natural population structure of Campylobacter jejuni colonizing a farmland ecosystem. Environ Microbiol. 2005, 7: 1116-1126. 10.1111/j.1462-2920.2005.00782.x.
Parkhill J, Wren BW, Mungall K, Ketley JM, Churcher C, Basham D, Chillingworth T, Davies RM, Feltwell T, Holroyd S, Jagels K, Karlyshev AV, Moule S, Pallen MJ, Penn CW, Quail MA, Rajandream MA, Rutherford KM, van Vliet AH, Whitehead S, Barrell BG: The genome sequence of the food-borne pathogen Campylobacter jejuni reveals hypervariable sequences. Nature. 2000, 403: 665-668. 10.1038/35001088.
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-10.1371/journal.pbio.0030015.
Hofreuter D, Tsai J, Watson RO, Novik V, Altman B, Benitez M, Clark C, Perbost C, Jarvie T, Du L, Galan JE: Unique features of a highly pathogenic Campylobacter jejuni strain. Infect Immun. 2006, 74: 4694-4707. 10.1128/IAI.00210-06.
Poly F, Threadgill D, Stintzi A: Genomic diversity in Campylobacter jejuni: identification of C. jejuni 81-176-specific genes. J Clin Microbiol. 2005, 43: 2330-2338. 10.1128/JCM.43.5.2330-2338.2005.
Genomes OnLine Database v 2.0. [http://www.genomesonline.org/]
Poly F, Threadgill D, Stintzi A: Identification of Campylobacter jejuni ATCC 43431-specific genes by whole microbial genome comparisons. J Bacteriol. 2004, 186: 4781-4795. 10.1128/JB.186.14.4781-4795.2004.
Karlyshev AV, Ketley JM, Wren BW: The Campylobacter jejuni glycome. FEMS Microbiol Rev. 2005, 29: 377-390. 10.1016/j.femsre.2005.01.003.
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. 10.1099/mic.0.27327-0.
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. 10.1016/S0014-5793(03)01164-5.
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. 10.1101/gr.185801.
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. 10.1086/368268.
SL WO, Dorrell N, Petersen L, Bang DD, Morris S, Forsythe SJ, Wren BW: Numerical analysis of DNA microarray data of Campylobacter jejuni strains correlated with survival, cytolethal distending toxin and haemolysin analyses. Int J Med Microbiol. 2006
Dorrell N, Hinchliffe SJ, Wren BW: Comparative phylogenomics of pathogenic bacteria by microarray analysis. Curr Opin Microbiol. 2005, 8: 620-626. 10.1016/j.mib.2005.08.012.
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. 10.1073/pnas.0503252102.
Agron PG, Macht M, Radnedge L, Skowronski EW, Miller W, Andersen GL: Use of subtractive hybridization for comprehensive surveys of prokaryotic genome differences. FEMS Microbiol Lett. 2002, 211: 175-182. 10.1111/j.1574-6968.2002.tb11221.x.
Winstanley C: Spot the difference: applications of subtractive hybridisation to the study of bacterial pathogens. J Med Microbiol. 2002, 51: 459-467.
Ahmed IH, Manning G, Wassenaar TM, Cawthraw S, Newell DG: Identification of genetic differences between two Campylobacter jejuni strains with different colonization potentials. Microbiology. 2002, 148: 1203-1212.
Tehara SK, Keasling JD: Gene Cloning, purification, and characterization of a phosphodiesterase from Delftia acidovorans. Appl Environ Microbiol. 2003, 69: 504-508. 10.1128/AEM.69.1.504-508.2003.
Parker CT, Quinones B, Miller WG, Horn ST, Mandrell RE: Comparative genomic analysis of Campylobacter jejuni strains reveals diversity due to genomic elements similar to those present in C. jejuni strain RM1221. J Clin Microbiol. 2006
Bacon DJ, Alm RA, Burr DH, Hu L, Kopecko DJ, Ewing CP, Trust TJ, Guerry P: Involvement of a plasmid in virulence of Campylobacter jejuni 81-176. Infect Immun. 2000, 68: 4384-4390. 10.1128/IAI.68.8.4384-4390.2000.
Batchelor RA, Pearson BM, Friis LM, Guerry P, Wells JM: Nucleotide sequences and comparison of two large conjugative plasmids from different Campylobacter species. Microbiology. 2004, 150: 3507-3517. 10.1099/mic.0.27112-0.
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. 10.1128/JCM.42.10.4566-4576.2004.
Boer P, Wagenaar JA, Achterberg RP, Putten JP, Schouls LM, Duim B: Generation of Campylobacter jejuni genetic diversity in vivo. Mol Microbiol. 2002, 44: 351-359. 10.1046/j.1365-2958.2002.02930.x.
Chevalier C, Thiberge JM, Ferrero RL, Labigne A: Essential role of Helicobacter pylori gamma-glutamyltranspeptidase for the colonization of the gastric mucosa of mice. Mol Microbiol. 1999, 31: 1359-1372. 10.1046/j.1365-2958.1999.01271.x.
McGovern KJ, Blanchard TG, Gutierrez JA, Czinn SJ, Krakowka S, Youngman P: gamma-Glutamyltransferase is a Helicobacter pylori virulence factor but is not essential for colonization. Infect Immun. 2001, 69: 4168-4173. 10.1128/IAI.69.6.4168-4173.2001.
Sellars MJ, Hall SJ, Kelly DJ: Growth of Campylobacter jejuni supported by respiration of fumarate, nitrate, nitrite, trimethylamine-N-oxide, or dimethyl sulfoxide requires oxygen. J Bacteriol. 2002, 184: 4187-4196. 10.1128/JB.184.15.4187-4196.2002.
Myers JD, Kelly DJ: A sulphite respiration system in the chemoheterotrophic human pathogen Campylobacter jejuni. Microbiology. 2005, 151: 233-242. 10.1099/mic.0.27573-0.
Pittman MS, Kelly DJ: Electron transport through nitrate and nitrite reductases in Campylobacter jejuni. Biochem Soc Trans. 2005, 33: 190-192. 10.1042/BST0330582.
Gaynor EC, Cawthraw S, Manning G, MacKichan JK, Falkow S, Newell DG: The genome-sequenced variant of Campylobacter jejuni NCTC 11168 and the original clonal clinical isolate differ markedly in colonization, gene expression, and virulence-associated phenotypes. J Bacteriol. 2004, 186: 503-517. 10.1128/JB.186.2.503-517.2004.
PubMLST database. [http://pubmlst.org/campylobacter/]
Winstanley C, Hart CA: Presence of type III secretion genes in Burkholderia pseudomallei correlates with Ara(-) phenotypes. J Clin Microbiol. 2000, 38: 883-885.
Campylobacter jejuni NCTC11168 genome sequence website. [http://www.sanger.ac.uk/Projects/C_jejuni/]
NCBI website. [http://www.ncbi.nlm.nih.gov/]
Multilocus sequence typing website. [http://www.mlst.net/]
This work was funded by DEFRA/HEFC as part of the UK Veterinary Training and Research Initiative. This publication made use of the Multi Locus Sequence Typing website  at Imperial College London developed by David Aanensen and Man-Suen Chan and funded by the Wellcome Trust.
PH carried out most of the laboratory work and participated in data analysis. HL carried out the MLST typing. CAH participated in the design, supervision and coordination of the study. CW participated in the design, supervision and coordination of the study, and drafted the manuscript. All authors read and approved the final manuscript.