Emerg Infect DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention18826829260988707-110910.3201/eid1410.071109DispatchExtended Sequence Typing of Campylobacter spp., United Kingdom Extended Sequence Typing of Campylobacter spp.DingleKate E.McCarthyNoel D.CodyAlison J.PetoTim E.A.MaidenMartin C. J.John Radcliffe Hospital, Oxford, UK (K.E. Dingle, T.E.A. Peto)University of Oxford, Oxford (K.E. Dingle, N.D. McCarthy, A.J. Cody, T.E.A. Peto, M.C.J. Maiden)Address for correspondence: Martin C.J. Maiden, Department of Zoology, The Peter Medawar Bldg, University of Oxford, South Parks Rd, Oxford OX1 3SY, UK; email: martin.maiden@zoo.ox.ac.uk102008141016201622

Supplementing Campylobacter spp. multilocus sequence typing with nucleotide sequence typing of 3 antigen genes increased the discriminatory index achieved from 0.975 to 0.992 among 620 clinical isolates from Oxfordshire, United Kingdom. This enhanced typing scheme enabled identification of clusters and retained data required for long-range epidemiologic comparisons of isolates.

Keywords: Campylobacter jejuniCampylobacter colimultilocus sequence typingmolecular epidemiologyantigen sequence typingdispatch

Human campylobacteriosis remains a global public health problem. Although many risk factors for this foodborne zoonotic disease are known, the relative contributions of different transmission routes are poorly quantified. Furthermore, the sources of particular infections are frequently obscure and outbreaks are rarely detected (1). The high genetic and antigenic diversity of the 2 major causes of campylobacteriosis, Campylobacter jejuni and C. coli, have proved to be obstacles in routine surveillance, outbreak identification, and source attribution.

Nucleotide sequence–based isolate characterization methods such as multilocus sequence typing (MLST) successfully catalog bacterial pathogens and provide a rational, definitive, and portable typing method with complete reproducibility among laboratories (2). Because many of the sequence types (STs) or their close relatives are observed on multiple occasions with wide geographic distribution, MLST is highly effective for long-range epidemiologic studies (3). However, this characteristic can limit the application of MLST to outbreak identification (4). We combined MLST data with sequence data of the short variable region (SVR) of the flaA and flaB loci, previously used to type Campylobacter spp./isolates (5,6), and a novel typing system based on porA locus encoding the variable outer membrane protein PorA. The resultant high-resolution 10-locus typing scheme was used to characterize 620 isolates obtained from 584 persons with human campylobacteriosis from September 2003 through September 2004 in Oxfordshire (population ≈600,000), United Kingdom; 36 isolates obtained by sampling the same patient more than once were used to confirm reproducibility.

The Study

A comparison of our results with published population-based 7-locus MLST studies showed that the relative abundance of different clonal complexes in northwestern England from April 2003 through March 2004 was similar (7), presumably reflecting widely distributed foods in the United Kingdom. An exception was ST-574 complex, the central genotype, which represented >5% of cases in Oxfordshire but was absent from northwestern England. The clonal complex distribution of 171 isolates collected in New South Wales, Australia (8), was distinct from the 2 English datasets, although many clonal complexes were present in all 3 datasets. Fewer clonal complexes, with different relative abundances, were present in a dataset from Curaçao (9), likely because of different infection sources in the Dutch West Indies (Figure 1).

Relative abundance of clonal complexes of Campylobacter spp. detected in Oxfordshire, United Kingdom, during a 1-year study compared with clonal complexes detected in 3 other studies of human Campylobacter spp. infections in northwestern England (7), New South Wales, Australia (8), and Curaçao (9).

The differences in relative abundance of clonal complexes were mirrored by FST values calculated from concatenated nucleotide sequences of the MLST loci, which indicated that the Australian dataset was 3.6% differentiated and the Curaçao set 9.9% differentiated from the Oxfordshire dataset. One clonal complex detected in Curaçao was absent in the United Kingdom (ST-41 complex) and 13 of the clonal complexes detected in the United Kingdom were absent in Curaçao.

The antigen loci added further resolution to the Oxfordshire dataset: 575 (98.6% of isolates) flaA SVR sequences contained 130 distinct SVR alleles, and 567 (97.1%) flaB sequences contained 111 SVR alleles. The allele fragment ≈630bp used for porA typing was amplified and sequenced with primers MOMP-1 (5′-GAT GGT TTA ACT CTA GCT GC-3′) or MOMP-3 (5′-GAT GGT TTA GTW GGM ACA GG-3′) and MOMP-2 (5′-TGA GAA GTT AAG TTT TGG AG AG-3′). PorA allele and variant numbers were assigned and sequences were deposited in a database (http://hercules.medawar.ox.ac.uk/momp). A description of amplification conditions is also available at this website. Of the 558 (95.5% of the isolates) porA sequences assigned, 1 occurred in 65 isolates whereas 135 alleles occurred only once.

The discriminatory index (DI) (10) was calculated for various subsets of these data. Each additional antigen gene increased discrimination relative to MLST data alone. The combination of flaA SVR and flaB SVR (DI = 0.976) provide a similar level of discrimination as porA (DI = 0.972). The DI obtained with the porA gene fragment alone was similar to that obtained previously with a larger fragment of the same gene (DI = 0.973) (11). The 10-locus combination provided a degree of discrimination (DI = 0.992) higher than those published for pulsed-field gel electrophoresis fingerprinting, antigen typing, or MLST (4,6,9,12,13). These studies calculated DI, ignoring the fact that some isolates probably shared a common source and may have underestimated the true DI, the capacity to discriminate between epidemiologically unrelated isolates. For the same reason, our study is also likely to have underestimated the true DI.

There were 68 groups of >2 isolates with identical 10-locus types, ranging in size from 2 to 34 isolates and accounting for 283 (48.5%) of the independent isolates. Of the remaining isolates, 290 (49.7%) had unique types, typing data were incomplete for the remaining 10 isolates (1.7%). A permutation test with 283 isolates belonging to a cluster showed highly significant temporal clustering of Campylobacter spp. isolates of identical genotype (p<0.0001). The extent of clustering was independent of group size (data not shown). Of 16 groups of >5 identical isolates, 5 exhibited significant temporal clustering (Table, Figure 2). Isolates belonging to the largest of these groups, comprising 34 isolates of ST-257, flaA 16, flaB 301, and porA 1, were submitted mainly in the last part of the study year with a peak of 5 isolates in week 39 (Figure 2). The second largest group, comprising 13 isolates, shared ST-51 and an identical porA type with a genotypic group comprising 7 isolates, but the 2 groups were distinct at both flaA and flaB loci. All members of these 2 groups were isolated from week 14 through week 40. Of the 12 isolates comprising the third largest group, all but 1 were isolated over a 17-week period (week 2 through week 18) at the beginning of the year (Figure 2). The smallest group of isolates to show evidence for temporal clustering comprised 5 identical isolates obtained during weeks 18–29. Some of the other genotypic groups were seen throughout the year, with no evidence for temporal clustering; for others, weak evidence of clustering was found (Table).

Temporal association of genotypically identical isolates of <italic>Campylobacter</italic> spp., United Kingdom*
ST10-locus genotype
No. isolatesp value of temporal association
flaA SVRflab SVRporA
4911115350.0005
2061496750.36
5832391774350.85
45884460.92
48321031460.44
35418185760.42
5121211070.0007
4751051056770.16
503636680.35
8272552363380.49
193636790.15
658552590.19
104363614110.14
574105105112<0.0001
513162951013<0.0001
25716301134<0.0001

*ST, sequence type; SVR, short variable region.

Clusters of related 10-locus types of Campylobacter spp. detected in Oxfordshire, United Kingdom, during a 1-year study. Five groups of isolates with identical genotypes show statistically significant clustering in time (p values are shown in the Table). Each group is indicated by 1 color. White bars indicate other isolates that share the same sequence type (ST) but that are differentiated by their different antigen type. Numbers of isolates of each genotype are shown on a weekly basis; week 1 corresponds to the start of the study on September 15, 2003.

Conclusions

Clusters were detected by this method but represented a small part of the overall disease incidence (14). Temporal association within these clusters suggests that they may represent undetected outbreaks. Further epidemiologic information unavailable to this study, which was based solely on laboratory isolates, would be needed to confirm or refute this possibility.

The clusters occurred over periods longer than the typical duration of outbreaks of gastroenteritis, which was consistent with episodes of contamination entering the food chain rather than single proximate-point source events. These findings indicated that many of these clusters may be associated with widely distributed foods. The observed temporal association of groups of identical isolates could also be caused by certain genotypes having different seasonality or environmental sources from sporadic cases. In either scenario, the data supported the interpretation that cases sharing an identical genotype were more likely to be epidemiologically linked than those infected with different genotypes.

In conclusion, the 10-locus typing scheme is highly discriminatory for C. jejuni and C. coli isolates and provides information that can be used flexibly to support long-range comparisons and short-term epidemiology. The scheme can be applied in real time or near real time, enabling the data to be used to identify outbreaks and inform public health interventions. When combined with improved genetic methods of attributing source of bacterial isolates (15), this approach will contribute to refining the epidemiology of these enigmatic pathogens.

Suggested citation for this article: Dingle KE, McCarthy ND, Cody AJ, Peto TEA, Maiden MCJ. Extended sequence typing of Campylobacter spp., United Kingdom. Emerg Infect Dis [serial on the Internet]. 2008 Oct [date cited]. Available from http://www.cdc.gov/EID/content/14/10/1620.htm

Acknowledgments

We thank the staff of the Clinical Microbiology Laboratory, John Radcliffe Hospital, Oxford, UK, for assistance in collecting Campylobacter spp. isolates; Keith Jolley for establishing the major outer membrane protein (MOMP) database; and Lynne Richardson for assistance with sequencing.

This study was supported in part by the Oxford University Medical Research Fund. This study used the Internet-accessible databases (http://pubmlst.org/campylobacter), the maintenance, development, and curation of which is supported by the Wellcome Trust and the Department for Environment, Food and Rural Affairs Research Contract no. OZ0611.

Dr Dingle is a senior research scientist at the Nuffield Department of Clinical Laboratory Sciences, University of Oxford, United Kingdom. Her research interests include the molecular epidemiology and population structure of Campylobacter species and norovirus epidemiology.

ReferencesHumphrey T, O’Brien S, Madsen M Campylobacters as zoonotic pathogens: a food production perspective.Int J Food Microbiol 2007;117:23757 10.1016/j.ijfoodmicro.2007.01.00617368847Dingle KE, Colles FM, Wareing DRA, Ure R, Fox AJ, Bolton FJ, Multilocus sequence typing system for Campylobacter jejuni.J Clin Microbiol 2001;39:1423 10.1128/JCM.39.1.14-23.200111136741Dingle KE, Colles FM, Ure R, Wagenaar J, Duim B, Bolton FJ, Molecular characterization of Campylobacter jejuni clones: a rational basis for epidemiological investigations.Emerg Infect Dis 2002;8:9495512194772Sails 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:47339 10.1128/JCM.41.10.4733-4739.200314532212Meinersmann RJ, Helsel LO, Fields PI, Hiett KL Discrimination of Campylobacter jejuni isolates by fla gene sequencing.J Clin Microbiol 1997;35:281049350739Mellmann A, Mosters J, Bartelt E, Roggentin P, Ammon A, Friedrich AW, Sequence-based typing of flaB is a more stable screening tool than typing of flaA for monitoring of Campylobacter populations.J Clin Microbiol 2004;42:48402 10.1128/JCM.42.10.4840-4842.200415472357Sopwith W, Birtles A, Matthews M, Fox A, Gee S, Painter M, Campylobacter jejuni multilocus sequence types in humans, northwest England, 2003–2004.Emerg Infect Dis 2006;12:1500717176563Mickan L, Doyle R, Valcanis M, Dingle KE, Unicomb L, Lanser J Multilocus sequence typing of Campylobacter jejuni isolates from New South Wales, Australia.J Appl Microbiol 2007;102:14452 10.1111/j.1365-2672.2006.03049.x17184329Duim B, Godschalk PC, van den Braak N, Dingle KE, Dijkstra JR, Leyde E, Molecular evidence for dissemination of unique Campylobacter jejuni clones in Curacao, Netherlands Antilles.J Clin Microbiol 2003;41:55937 10.1128/JCM.41.12.5593-5597.200314662946Hunter PR, Gaston MA Numerical index of discriminatory ability of typing systems: an application of Simpson’s index of diversity.J Clin Microbiol 1988;26:246563069867Huang S, Luangtongkum T, Morishita TY, Zhang Q Molecular typing of Campylobacter strains using the cmp gene encoding the major outer membrane protein.Foodborne Pathog Dis 2005;2:1223 10.1089/fpd.2005.2.1215992295Schouls LM, Reulen S, Duim B, Wagenaar JA, Willems RJ, Dingle KE, 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:1526 10.1128/JCM.41.1.15-26.200312517820Thakur S, Morrow WE, Funk JA, Bahnson PB, Gebreyes WA Molecular epidemiologic investigation of Campylobacter coli in swine production systems, using multilocus sequence typing.Appl Environ Microbiol 2006;72:56669 10.1128/AEM.00658-0616885327Pebody RG, Ryan MJ, Wall PG Outbreaks of campylobacter infection: rare events for a common pathogen.Commun Dis Rep CDR Rev 1997;7:R3379080726McCarthy ND, Colles FM, Dingle KE, Bagnall MC, Manning G, Maiden MC, Host-associated genetic import in Campylobacter jejuni.Emerg Infect Dis 2007;13:2677217479890