Human infection is usually acquired by the consumption of contaminated food (especially poultry) or water¹. Motile campylobacters colonize the intestines of a wide range of animals, but in immunologically naive humans infection frequently results in an inflammatory enterocolitis. The number of cases ofCampylobacter infection reported in England and Wales in 1998 increased by 17% from the previous year, with the number of reported cases now more than double that due toSalmonella³. Despite its importance, effective control ofCampylobacter in the food chain and the design of disease prevention strategies are hindered by a poor understanding of the genetics, physiology and virulence of this organism.

The genome ofC. jejuni NCTC11168 is 1,641,481 base pairs (bp) in length. Of the 1,654 predicted coding sequences (CDS), at least 20 probably represent pseudogenes; the average gene length is 948 bp, and 94.3% of the genome codes for proteins, making it the densest bacterial genome sequenced to date. The bias towards G on the leading strand of the chromosome⁴ indicates that the origin of replication is near to the start of thednaA gene. Strand bias is also evident in the CDSs; overall, 61.1% are transcribed in the same direction as replication (Fig. 1). We discovered two large regions of lower G+C content that encompass CDSs Cj1135–Cj1148 (25.4%) and Cj1421–Cj1442 (26.5%); these correspond to genes within the lipooligosaccharide (LOS) and extracellular polysaccharide (EP) biosynthesis clusters, respectively. Functional information (matches to genes of known function, or informative hydrophobicity profiles) could be deduced for 77.8% of the 1,654 CDSs, wheras 13.5% matched genes of unknown function in the database and 8.7% had no database match, or other functional information. The unusually low number of unknowns reflects the preponderance of predicted membrane, periplasmic and lipoproteins; these make up 10.3%, 7.8% and 2.3 % of the CDSs, respectively, and many of these have no database matches.

From the outside to the inside, the first circle shows coding sequences transcribed in the clockwise direction in dark green; the second shows coding sequences in the anticlockwise direction in pale green. The putative origin of replication is marked. The third shows the positions of hypervariable sequences in black, and the fourth and fifth show genes involved in the production of surface structures: clockwise in dark red and anticlockwise in pale red. The innermost histogram shows the similarity of each gene to itsH. pylori orthologue, where present; the height of the bar, and the intensity of the colour, are proportional to the degree of similarity. The clusters of genes responsible for LOS biosynthesis, EP biosynthesis and flagellar modification are marked.

Full size image

One surprising feature of theC. jejuni genome is the almost complete lack of repetitive DNA sequences. In fact, there are only four repeated sequences within the entire genome; three copies of the ribosomal RNA operon (6 kilobases (kb)) and three duplicated or triplicated CDSs. Apart from Cj0752, which is similar to part of IS605tnpB fromH. pylori, there is no evidence of any functional inserted sequence (IS) elements, transposons, retrons or prophages in the genome. Another intriguing feature of the genome is that, apart from salient exceptions such as the two ribosomal protein operons and gene clusters involved in LOS biosynthesis, EP biosynthesis and flagellar modification, there appears to be little organization of genes into operons or clusters. Although genes do fall into long, apparently linked sets, generally the genes within these sets appear to be functionally unrelated. The distribution of the genes involved in amino-acid biosynthesis, for example, reflects this organization: some of thehis, leu andtrp genes are apparently organized into operons, but thearo, asp, dap, gln, gly, ilv, met, phe, pro, ser, thr and tyr genes are scattered randomly throughout the genome.

The shotgun assembly revealed regions in which the sequences of otherwise identical clones varied at a single point. These were mainly, but not exclusively, length variations in polyG:C tracts (Table 1). The degree of variation differs between sites, and there are homopolymeric tracts that do not display variability within the observed shotgun sequences (Table 2). Variation in the length of polyG:C tracts is frequently associated with contingency genes in other pathogenic bacteria, and may be produced by slipped-strand mispairing during replication⁵; it can affect translation and has been shown to be responsible for phase variation of surface properties or antigenicity⁶. The appearance of variants due to slipped-strand mispairing occurs inNeisseria menigitidis with a frequency of 10^-3 per cell per generation⁷; these variants are therefore likely to be rare in the sampled clones of a 10-fold shotgun library derived from a clonal population. TheC. jejuni sequence, by contrast, shows some regions where three or more variants are present in almost equal proportions, in addition to many regions where two variants are present. The absolute rate of variation is difficult to estimate from these data, although the number of changes seen suggests that the frequency is much higher than in other organisms.

The variation in homopolymeric tracts was not an artefact of the sequencing process, as fourfold resequencing of several variants did not show any difference from the original sequence. A shotgun sequence of an appropriate lambda clone demonstrated that the variation did not occur during subcloning inEscherichia coli. Support for the existence of rapid phase variation inC. jejuni is also provided byH. pylori J99, for which similar variation was seen⁸. This rapid sequence variation suggests thatC. jejuni may be lacking in DNA repair functions, indeed, many DNA repair genes studied inE. coli cannot be found inC. jejuni, including the direct repair genesada andphr; the glycosylasestag, alkA, mutM andnfo; the mismatch repair genesvsr, mutH, mutL and, sbcB; and the SOS response geneslexA, umuC and umuD. Significantly, transcription-coupled repair may influence the rate of phase variation inN. meningitidis⁷. The high levels of variation seen in the shotgun sequences mean that it is not possible to produce a single definitive sequence for theC. jejuni genome. As such, it possesses some of the properties of a quasispecies; a phenomenon that is well described in RNA viruses⁹.

Most of the hypervariable sequences cluster on the genome and are coincident with the clusters of genes responsible for LOS biosynthesis, EP biosynthesis and flagellar modification (Fig. 1). Some of the variable genes can be ascribed putative functions, such as glycosyl-transferases; several others belong to two families (designated as 617 and 1318; Fig. 2,Table 3). The 617 family of genes have no homologues outsideC. jejuni, while the 1318 family has two homologues inH. pylori (HP0114 and HP0465). Sixteen additional members of the 1318 family occur in various bacterial and archaeal species, none with a well-defined function; however, family members within the enterobacteriaceae are found within lipopolysaccharide gene clusters, supporting our hypothesis that these proteins are involved in the synthesis of surface structures. The rapid variation in surface properties implied by these results may have more relevance to the colonization of a dynamic intestinal environment than to immune avoidance.

Amino acids encoded by the hypervariable homopolymeric tracts are in boxes. Cj0617 and Cj0618, and Cj1335 and Cj1336 represent coding sequences frameshifted at the homopolymeric tract.

Full size image

C. jejuni has been reported to produce a variety of toxins whose activity and/or role in pathogenesis remain controversial¹⁰. The genome of NCTC11168 does not contain a cholera-like toxin gene, although genes encoding the cytolethal distending toxin (cdtA-C) are present. A member of the family of contact-dependent haemolysins found in pathogenicSerpulina andMycobacterium species¹¹ (Cj0588), a putative integral membrane protein with a haemolysin domain (Cj0183) and a phospholipase (pldA) were also identified.

In contrast to most lipid A and some inner core biosynthesis genes, many LOS biosynthesis proteins are encoded by a large gene cluster (Cj1119–Cj1152) which has a role in core biosynthesis and also protein glycosylation^12,13. There is increasing evidence thatC. jejuni synthesizes an EP that is not attached to lipid A¹³, and thus the presence of genes similar to those involved in bacterial capsule biogenesis is significant. Two groups containingkps orthologues (involved in transport; A.V.K., manuscript submitted) were found and the region between the groups contains many different polysaccharide biosynthetic genes, some required for the biosynthesis of bacterial capsules. As might be expected with the presence ofkps genes, no orthologues towzx,wzy andwzz genes, which are essential for heteropolymeric O-chain biosynthesis, are present. Unusually, C. jejuni has three sets ofneu genes involved in sialic acid biosynthesis. Sialic acid is an uncommon constituent of bacterial surface structures and, through molecular mimicry, may be important for evasion of host immunity and in post-infection autoimmune diseases such as Guillain–Barré syndrome. A complete cluster (neuB1, C1, A1) has a role in LOS sialylation (D. Lintonet al., personal communication). Another set (neuB2 andC2) is in close proximity toptmAB and may be involved with these genes in post-translational modification of flagellin¹⁴.

There are no obvious orthologues of the extensive Hop porin family of H. pylori¹⁵ and, contrary to a recent report¹⁶, no type III secretion systems were identified other than the flagellin export apparatus. In the absence of such a system, it is possible that the CiaB protein¹⁶ is secreted by the flagellin export apparatus¹⁷. In contrast toCampylobacter fetus andCampylobacter rectus, genes encoding S-layer proteins were apparently absent inC. jejuni. Under certain conditions,C. jejuni produce 4–7 nm wide filaments resembling pili¹⁸; although structural pilin orthologues appear to be absent, there are several type 4 pilus related genes such as a pre-pilin peptidase (Cj0825) and several putative type II export genes (Cj1470c–Cj1474c). Flagella are one of the best-defined virulence factors, asflaA andflaB mutants are markedly reduced in virulence¹⁹. Whether this is due to a direct involvement of the flagellin, or perhaps an inability to secrete proteins through the flagella export apparatus, is unknown¹⁷. Two further structural flagellin paralogues were identified (flaC Cj0720c andflaD Cj0887c), and orthologues of the majority of flagellar-associated genes of clearly understood function in bacteria such asS. typhimurium are represented; however, several genes involved in regulation are absent, includingflhCD,flgM, fliT andfliK. These differences from the enterobacterial paradigm mirror those found inH. pylori. However, the gene encoding the H. pylori flagellar sheath proteinhpaA (ref. 20) is absent inC. jejuni, as is a flagellar sheath.

Chemotaxis is important for intestinal colonization byC. jejuni. Like H. pylori, C. jejuni produces three proteins containing the response regulatory domain of CheY. One is CheY (Cj1118)²¹, and the others are found fused to a histidine protein kinase domain (CheA) or a CheW-like domain (CheV). In contrast toH. pylori, only one CheV orthologue is present. The regulation of CheY activity inC. jejuni is different from that ofH. pylori as, unlikeH. pylori, orthologues of both CheB and CheR are present. No orthologues of CheZ or other chemotaxis genes were found. Ten genes contain methyl-accepting chemotaxis protein domains; these are candidate chemoreceptor genes, some of which may transduce signals to non-taxis-associated pathways; only three (Cj1506, Cj0448, Cj0019) have orthologues inH. pylori.

The genome encodes five major iron-acquisition systems, which are mostly organized in operons under the control of the Fur protein²². Two of these iron-acquisition systems, the enterochelin uptake operon ceuBCDE and the siderophore receptor orthologuecfrA, have been described previously. In addition, there is a predicted haemin uptake operon chu (Cj1614–1617)²², a periplasmic binding protein dependent system (Cj0173c–0175c) and a siderophore receptor (Cj0178) with accessory genes. Three copies of the accessorytonB, exbB and exbD genes, which form the energy transduction machinery for transport of iron compounds over the outer membrane, were found. One of theexbB-exbD-tonB triplets follows the Cj0178 siderophore receptor, whereas another tonB gene is divergently orientated tocfrA.C. jejuni might therefore use separate energy transduction mechanisms for transport of the different iron substrates, as has been shown forVibrio cholerae.

C. jejuni appears to have a broader repertoire of regulatory systems thanH. pylori, which has a similar sized genome. Given thatC. jejuni is found in a more diverse range of ecological niches than H. pylori, this might be expected. The apparent lack of operon organization raises fundamental issues concerning transcriptional regulation inCampylobacter. Although it is possible that each gene is independently transcribed, strand-specific gene grouping suggests co-transcription. LikeH. pylori, theCampylobacter genome contains only three predicted sigma factors (rpoD,rpoN andfliA). The largest proportion of regulatory genes consists of members of the two-component regulator family. As inBacillus subtilis,C. jejuni has an additional member of the Fur²² family, PerR, which regulates the peroxide stress regulon²³. UnlikeH. pylori,C. jejuni contains a possiblecrp/fnr family member (Cj0466) with both helix–turn–helix and cNMP-binding motifs.

As might be expected from the inability ofC. jejuni to use carbohydrates as carbon or energy sources, very few genes for degradation of carbohydrates or amino acids were detected. The glycolytic pathway is also apparently incomplete; orthologues of the glucokinase and 6-phosphofructokinase genes were not found, although these functions may well be supplied by non-orthologous genes. Despite this,C. jejuni does appear to have all the genes necessary for gluconeogenesis and, unlikeH. pylori, C. jejuni appears to encode an intact tricarboxylic acid cycle (likeH. pylori,C. jejuni uses 2-oxoglutarate ferredoxin oxidoreductase (OorDABC)²⁴, rather than 2-oxoglutarate dehydrogenase to interconvert 2-oxoglutarate and succinyl CoA).

C. jejuni andH. pylori are closely related by 16S rRNA phylogeny, share many biological properties and were previously classified within the Campylobacter genus.Figure 1 (and Supplementary Information) indicates whichC. jejuni genes are present or absent inH. pylori 26695 (ref.15). The three polysaccharide biosynthetic loci relating to surface structure stand out as being unique toC. jejuni and are correlated with a high density of polymorphic sequences.C. jejuni also contains a large number of biosynthetic genes not present in H. pylori, such as those directing the synthesis of purines, thiamine and many amino acids. Genes present inH. pylori but absent inC. jejuni include the urease operon, nickel transport system, vacuolating toxin and the Cag pathogenicity island¹⁵. These features are consistent with the unique niche ofH. pylori in the stomach, and its pathophysiology and propensity for chronic infection.

Despite the close phylogenetic relationship ofC. jejuni and H. pylori, strong similarities between them are mainly confined to housekeeping functions; only 55.4% ofC. jejuni genes have orthologues inH. pylori. In most functions related to survival, transmission and pathogenesis, the organisms have remarkably little in common. This indicates that selective pressures have driven profound evolutionary changes to create two very different and specific pathogens appropriate to their niches, from a relatively close common ancestor. Overall, 28.0% of genes show closest similarity to genes fromE. coli, 27.0% to genes fromB. subtilis, 4.6% to genes fromArcheoglobus fulgidus and 2.1% to genes fromSaccharomyces cerevisiae. Several genes were found which have homologues only in the eukaryotic domain; these include the dUTPase (dut) gene (Cj1451), which is similar to those ofLeishmania major andTrypanosoma cruzi only; there is no orthologue of theE. coli dut gene. Taken together, these statistics suggest that the evolution of theC. jejuni genome does not neatly mirror that of its small subunit ribosomal RNA, and that the placement ofC. jejuni within the Gram-negative proteobacteria may rely on a simplified view of the evolutionary origin of its genome.

This genomic sequence provides the resources for a complete and detailed analysis of the pathogenic potential of this enigmatic pathogen. New insights into the biology ofC. jejuni include the identification of hypervariable sequences, lack of classical operon structure and repetitive DNA, and an unexpected capacity for polysaccharide production.

A single colony ofC. jejuni (NCTC11168, human origin, serotype O2, minimally passaged) was spread on one petri dish of CCDH agar (Campylobacter selective agar, Oxoid), and re-streaked on 10 CCDH agar plates. Cells were harvested and total DNA (10 mg) was isolated using proteinase K treatment followed by a phenol extraction procedure. The DNA was fragmented by sonication, size-fractionated on an agarose gel, and seven libraries were generated in pUC18 using size fractions ranging from 1.0 kb to 2.2 kb. Roughly 19,400 pUC clones were sequenced from both ends, using Dye-terminator chemistry on ABI 373 and 377 sequencing machines. The final assembly was generated from 33,900 reads, giving an∼10-fold coverage of the genome. Sequence assembly was accomplished using Phrap (P. Green, unpublished), and the sequencing was finished using GAP4²⁵. The assembly was verified by genomic PCR reactions across all repeats, in addition to 130 forward and reverse reads from a random library of∼10–16-kb fragments of genomic DNA cloned in lambdaFixII (Stratagene). In the final assembly, 0.13% of the genome was covered by a single clone only, and 0.11% was not sequenced on both strands, or with complementary sequencing chemistries

The DNA was compared with sequences in the EMBL database using BLASTN and BLASTX²⁶. Transfer RNAs were predicted by tRNAscan-SE²⁷. Potential CDSs were predicted using ORPHEUS²⁸ and GLIMMER²⁹ (both trained on an initial open reading frame set generated by ORPHEUS), and also stop-to-stop prediction; the results were combined. The predicted protein sequences were searched against a non-redundant protein database using WUBLASTP and FASTA. The complete six-frame translation was used to search PROSITE, and the predicted proteins compared against the PFAM³⁰ database of protein domain hidden Markov models. The results of all these analyses were assembled together using the Artemis sequence viewer (K.M.R., unpublished) and used to inform a manual annotation of the sequence and predicted proteins. Annotation was based, wherever possible, on characterized proteins or genes.

We would like to thank N. Loman for theCampylobacter Genome Browser; A. Brás, C. Constantinidou, D. Linton, J. Marchant, N. Oldfield and S. Park for assistance with interpretation of sequence data; A. Bateman for CDS clustering; and S. Bowman, D. Kelly, and D. Maskell for a critical reading of the manuscript. Research in the laboratories of B.W.W. and J.M.K. was funded by the BBSRC and sequencing at the Sanger Centre by the Wellcome Trust through its Beowulf Genomics initiative.

1. J. Parkhill: The complete sequence and annotation can be obtained from the EMBL database with the ID CJ11168 (accession number AL111168).

Authors and Affiliations

1. The Sanger Centre, The Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK

   J. Parkhill, K. Mungall, C. Churcher, D. Basham, T. Chillingworth, R. M. Davies, T. Feltwell, S. Holroyd, K. Jagels, S. Moule, M. A. Quail, M-A. Rajandream, K. M. Rutherford, S. Whitehead & B. G. Barrell

2. Department of Infectious and Tropical Diseases, The London School of Hygiene and Tropical Medicine, Keppel St, London, WCIE 7HT, UK

   B. W. Wren & A. V. Karlyshev

3. Department of Genetics, University of Leicester, University Road, Leicester, LE1 7RH, UK

   J. M. Ketley

4. Department of Microbiology and Immunobiology, Queen's University Belfast, Grosvenor Road, Belfast , BT12 6BN, UK

   M. J. Pallen

5. School of Biosciences, University of Birmingham, Birmingham, B15 2TT, UK

   C. W. Penn

6. Departments of Medical Microbiology and Gastroenterology, Faculty of Medicine, Vrije Universiteit, Van der Boechorststraat 7, Amsterdam, 1081 BT, The Netherlands

   A. H. M. van Vliet

Corresponding author

Correspondence toJ. Parkhill.

Reprints and permissions