Previous Article | Next Article ![]()
Journal of Clinical Microbiology, December 2005, p. 5888-5898, Vol. 43, No. 12
0095-1137/05/$08.00+0 doi:10.1128/JCM.43.12.5888-5898.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Kate E. Dingle,2,
Martin C. J. Maiden,3
Diane G. Newell,4
Linda van der Graaf-Van Bloois,1,5
Jos P. M. van Putten,5 and
Jaap A. Wagenaar1,5*
Animal Sciences Group, Division of Infectious Diseases (OIE Reference Laboratory for Campylobacteriosis), Lelystad, The Netherlands,1 Nuffield Department of Clinical Sciences, University of Oxford, Department of Microbiology, John Radcliffe Hospital, Oxford, United Kingdom,2 The Peter Medawar Building for Pathogen Research and Department of Zoology, University of Oxford, Oxford, United Kingdom,3 Veterinary Laboratories Agency, Weybridge, United Kingdom,4 Department of Infectious Diseases and Immunology (OIE Reference Laboratory for Campylobacteriosis), Faculty of Veterinary Medicine, Utrecht University, Utrecht, The Netherlands5
Received 15 July 2005/ Returned for modification 15 August 2005/ Accepted 19 September 2005
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The subspecies are conventionally defined with biochemical tests, namely, growth in 1% glycine (32) and H2S production from medium containing 0.02% cysteine (1). Despite these differences and host and niche preferences, studies that have used potentially discriminatory methods, including DNA-DNA hybridization (29, 37), serotyping (27), ribotyping (3), cellular fatty acid analysis (2), and whole-cell protein analysis (37), have been unable to distinguish C. fetus subsp. fetus and C. fetus subsp. venerealis. Therefore, from a taxonomic viewpoint it is questionable whether C. fetus should be divided into subspecies. However, a dichotomy has been detected among C. fetus serotype A and B strains (7, 22, 27). These serotypes reflect differences in both lipopolysaccharide structure and the type of surface layer protein (SLP) and surface array protein (Sap). The two serotypes correlate with the corresponding sap type, sapA or sapB (33). C. fetus subsp. venerealis strains are type A, whereas C. fetus subsp. fetus strains can be either serotype A or serotype B.
Several molecular methods, including pulsed-field gel electrophoresis (PFGE) fingerprinting (25, 31), amplified fragment length polymorphism (AFLP) analysis (39), and PCR amplification (15, 40), have been used for subspecies differentiation; but they give contradictory results. Recently, subspecies-specific random amplification of polymorphic DNA PCRs were described, but they were evaluated with a limited number of strains (34). A new C. fetus subsp. venerealis-specific PCR has been developed based on an AFLP marker (36). Subspecies differentiation is of such statutory importance that additional unambiguous typing methods are still required, preferably those that assist with epidemiological investigations. This may be useful for outbreak control and has been performed by using serotyping and PFGE. The currently available serotyping scheme has a low level of discrimination, as only a few lipopolysaccharide serotypes have been detected. However, PFGE has proven discriminatory and has been applied in epidemiological studies of both human and veterinary infections (11, 20, 23, 24, 28).
Multilocus
sequence typing (MLST) is a nucleotide sequence-based approach to
bacterial typing in which variations in
500-bp fragments of
housekeeping genes (generally seven) are indexed. MLST is useful for
studies of bacterial epidemiology and population genetics, including
studies of species of the genus Campylobacter
(4,
5,
9,
18,
21). MLST has the
advantages that variation in the target gene is accessed directly and
that the technology employed is both readily disseminated and highly
reproducible among laboratories
(19). Furthermore,
comparisons of the population structures and population genetics of
different species can be made by using MLST data, particularly if the
same loci are used for different species.
Here we describe an MLST scheme for C. fetus that uses the same loci as a previous scheme for C. jejuni and C. coli (4, 5). We tested this scheme using a diverse collection of isolates and examined the C. fetus population structure and genetic diversity, together with the suitability of MLST for studies of C. fetus subspecies differentiation and epidemiology.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Species identification and subspecies differentiation. Isolates were identified genotypically as C. fetus by using a previously described PCR protocol (15), with modifications (39). The subspecies were then identified genotypically by Cf C05 PCR (36); the PCR of Hum et al. (15); AFLP analysis (39); and biochemically, according to growth in the presence of 1% glycine (Merck, Darmstadt, Germany) (32) and H2S production from brucella broth (Difco, Becton Dickinson, Franklin Lakes, NJ) containing 0.02% cysteine (Sigma Aldrich, St. Louis, MO) (1, 17).
sap typing. sap typing was performed with PCR primers based on the sapA and sapB 5' conserved regions, as described by Tu et al. (33).
AFLP analysis. AFLP analysis was performed by a previous method, with modifications (6, 39). These enhancements comprised the use of a capillary-based sequencing system (ABI 3100; Applied Biosystems, Foster City, CA) and data analysis with Bionumerics 3.5 software (Applied Maths, St-Martens-Latem, Belgium). Chromosomal DNA was digested with the restriction enzymes HindIII and HhaI, and site-specific adaptors were ligated. A preselective PCR amplification was performed, followed by selective PCR with primers containing a 3' selective nucleotide and one fluorescently labeled (6-carboxyfluorescein) primer. PCR products underwent electrophoresis through an ABI 3100 capillary sequence system by using a performance optimized polymer-4 polymer matrix (Applied Biosystems). AFLP patterns were analyzed by importing the AFLP data (via the conversion software GeneScan 3.7 [Applied Maths]) into the Bionumerics 3.5 software. Patterns were normalized by referring to the molecular mass of the internal DNA marker (GeneScan-500 [carboxy X rhodamine]; Applied Biosystems). The genetic similarity between the patterns was calculated by using the Pearson product-moment correlation coefficient, and clustering analysis was performed by the unweighted pair group method with arithmetic averages for the region corresponding to 84.6% to 87.1% of the band pattern.
PFGE. PFGE was performed as described previously by using the restriction enzymes SmaI and SalI (36a)
Multilocus sequence typing. (i) Primer design. C. fetus sequences corresponding to the seven loci used in a previously described C. jejuni and C. coli MLST scheme were required: aspA, glnA, gltA, glyA, pgm, tkt, and uncA (4, 5). Various C. jejuni and C. coli MLST primers (Sigma-Genosys, Haverhill, United Kingdom) were employed in multiple combinations to identify primers that would amplify these sequences from C. fetus. A reduced annealing temperature of 48°C was used to improve the binding of the C. jejuni primers to heterologous C. fetus DNA. Of the many primer combinations tested, those that provided amplification products were as follows: for C. jejuni aspA, primers A9 and A10; for C. jejuni glnA, primers A1 and A2; for C. jejuni gltA, primers A1 and A2; for C. jejuni glyA, primers A1 and A2; for C. jejuni tkt, primers A5 and A4; and for C. coli uncA, primers A1 and A2 (the sequences of all primers except C. jejuni tkt A5 and C. coli uncA A1 and A2 are available at http://pubmlst.org/campylobacter/; the sequence of tkt A5 is 5'-TTTAAGTGCTGATATGGTGC-3', the sequence of C. coli uncA A1 is 5'-ATGTAGCCATAGGACAAAAGC-3', and the sequence of C. coli uncA A2 is 5'-CATTCTTGTCGCGTTCAGTTG-3'). The C. fetus sequence corresponding to the pgm locus failed to amplify with any C. jejuni- or C. coli-specific primers; however, this sequence was kindly provided by D. Sanchez, Universidad Nacional de General San Martín, Instituto de Investigaciones Biotecnológicas, San Martín, Argentina. The C. fetus amplicons were sequenced directly by using the same primers used for amplification. Nucleotide sequence extension reactions were carried out by using BigDye ready reaction mix (version 3; Applied Biosystems), in accordance with the manufacturer's instructions. The reaction products were separated with an ABI 3730 automated DNA sequencer (Applied Biosystems).
These C. fetus sequences from the seven loci were used to design seven pairs of C. fetus-specific primers for MLST (Table 2). The MLST allele trimming sites for C. fetus were chosen to correspond to the previous C. jejuni and C. coli sites to assist with future phylogenetic comparisons among species. At least 28 nucleotides separated the MLST trimming site and the 3' end of the C. fetus primers. MLST was performed by using the same primer pair used for both amplification and sequencing.
|
(iii) Allele and ST assignment. For each locus, distinct allele sequences were assigned arbitrary allele numbers in the order of identification. Each genotype was therefore designated by seven numbers (e.g., 1-3-2-4-1-1-2) that constituted an allelic profile or sequence type (ST; e.g., ST-3). The STs were assigned arbitrary numbers in the order of description. All data for newly described C. fetus alleles and STs were deposited in a Campylobacter fetus MLST database (http://pubmlst.org/cfetus/). New sequences were assigned allele numbers, and isolates were assigned their STs by interrogating the database. Allele numbers for new sequences and ST numbers for new allelic profiles are available by submission to the database.
(iv) Data analysis. Phylogenetic analysis was performed with MEGA (version 2.1) software, available at http://www.megasoftware.net, by using the concatenated MLST gene sequence fragments for each isolate. Data were also subjected to split decomposition analysis by using SPLITSTREE (version 3.1) software (16).
| RESULTS |
|---|
|
|
|---|
All of the 140 isolates examined were typeable by MLST. The number of alleles per locus ranged from two to four, and the number of variable sites per locus ranged from one to five (Table 3). A total of 14 different STs were identified, some of which shared up to six or seven identical loci. The majority of isolates (91%) were assigned to one of the following STs: ST-2, ST-3, ST-4, ST-5, or ST-6 (Table 1). In this isolate collection, ST-1 was associated with an epidemiologically related group of isolates (n = 5), while the remaining STs (ST-7 to ST-14) each occurred only once.
|
|
|
|
The value of MLST for subspecies determination was assessed. Certain STs did correlate with each of the three AFLP groups. A total of 55 of 57 of the C. fetus subsp. venerealis or C. fetus subsp. venerealis bv. intermedius isolates were ST-4 (Table 1; Fig. 2). The remaining 2 of the 57 isolates were ST-7 and ST-12, respectively; ST-12differed by only one nucleotide from ST-4. The C. fetus subsp. fetus AFLP group showed a greater diversity of STs; but notably, ST-4, ST-7, and ST-12 were absent (Table 1; Fig. 2).
Investigation of C. fetus epidemiology by MLST. MLST data were examined for (i) the reliability of outbreak identification, (ii) any correlation of the ST with the isolation source or geographical region, and (iii) any correlation with and discrimination attained compared to the results of PFGE. Such data were available for a subset of epidemiologically related and unrelated isolates (n = 23).
The collection of 140 isolates studied included isolates from seven known outbreaks (with two strains per outbreak) from different bovine AI stations (Table 1). All but one of the outbreaks were associated with STs frequently identified within this data set (ST-1 to ST-4 and ST-6). The remaining outbreak was associated with a single bovine ST-5 isolate from an additional AI station. Epidemiologically related strains shared identical STs and PFGE types (e.g., ST-1 and ST-2; Table 1). However, some epidemiologically unrelated strains had identical STs, but PFGE typing revealed discriminatory patterns (e.g., ST-3, ST-4, ST-5, and ST-6; Table 1), which indicated that PFGE has a greater discriminatory power.
The MLST data and the isolate data on host and geographical region were evaluated for potential associations. Bovine, ovine, and human isolates were the most numerous; but none of these hosts harbored a single ST. However, ST-1 (n = 5) and ST-4 (n = 48) were isolated almost exclusively from bovine hosts, the exceptions being two human isolates. Seven further ST-4 isolates (indicated by a question mark in Table 1) were considered likely (but unconfirmed) bovine isolates (M. Henton, personal communication). Therefore, with the possible exception of ST-4, which is the C. fetus subsp. venerealis and C. fetus subsp. venerealis bv. intermedius ST and, therefore, by definition, cattle associated, no association between the ST and the host species was detected. The ST and the geographical region of isolation also lacked association, and all the common STs were found in Europe and North America.
| DISCUSSION |
|---|
|
|
|---|
C. fetus was originally divided into two subspecies by using host and niche preferences (38), but subsequent taxonomic studies led investigators to question this approach (14, 29). The difficulty of reliable subspecies differentiation reflects the fact that the two subspecies are extremely closely related genetically. For example, it proved impossible to differentiate subspecies by using the nucleotide sequences of two other housekeeping genes, sapD and recA (33), or the 16S rRNA gene (13, 25, 41). Since subspecies differentiation is of high statutory importance, the suitability of MLST for C. fetus subspecies differentiation was investigated. The C. fetus subsp. venerealis strains were all ST-4, although they were geographically and temporally distinct. The C. fetus subsp. venerealis bv. intermedius strains were also mostly ST-4, with two exceptions (one ST-7 strain and another ST-12 strain). Therefore, the MLST scheme did not allow differentiation between C. fetus subsp. venerealis and C. fetus subsp. venerealis bv. intermedius strains. However, ST-4 was not detected among the C. fetus subsp. fetus isolates. Extension of the MLST scheme to other loci may reveal differences between C. fetus subsp. venerealis and C. fetus subsp. venerealis bv. intermedius strains. The C. fetus subsp. venerealis "ST-4 clone" shows even less diversity within the housekeeping genes than the C. fetus subsp. fetus STs, supporting the theory of Véron and Chatelain, who proposed that C. fetus subsp. venerealis is a mutant of C. fetus subsp. fetus (38).
The data suggest that the MLST approach is reliable for subspecies differentiation since it was performed "blind" without prior knowledge of the subspecies assignments. Although C. fetus subsp. venerealis ST-4 differs from the closest C. fetus subsp. fetus ST (ST-6) by a single nucleotide (Fig. 2), the data generated by this highly reproducible sequencing technique confirm that MLST is an accurate method for the differentiation of the subspecies. Previous studies suggested that differences between C. fetus subsp. venerealis and C. fetus subsp. fetus are plasmid encoded (25, 26), although more recent studies have questioned this proposal (42). This MLST scheme indicates that the differences between the two subspecies are not merely plasmid encoded but are also reflected in the allelic profiles (STs) of the housekeeping genes.
C. fetus isolates from mammals are divided into two sap types, sapA and sapB, on the basis of their SLPs. The MLST data showed a clear correlation among STs and sap types (Fig. 2), supporting the classification of C. fetus into two sap groups. This correlation suggests evolutionary stability and a lack of frequent recombination between the two groups, since in C. jejuni (which recombines frequently) a lack of correlation among STs-clonal complexes and cell surface antigens has been detected.
MLST is very useful for the global epidemiological typing of many bacterial species, and its use is expanding rapidly (18, 35). We demonstrated in the present study that epidemiologically related C. fetus strains had identical STs. However, as a result of the relatively low level of genetic diversity found in C. fetus housekeeping genes, some unrelated isolates had STs identical to those of some of the outbreak strains. Sequencing of additional, more variable loci has proven useful in investigations of outbreaks of Neisseria meningitidis infections and may increase the discriminatory power of the method (10). Whether this approach would also improve C. fetus discrimination requires further investigation. C. fetus has a global distribution; however, no geographically distinct ST cluster could be distinguished by MLST. The major STs were found in multiple host species; however, the C. fetus subsp. venerealis ST-4 was restricted to bovine hosts. It has been suggested that C. fetus subsp. venerealis is a host-restricted mutant clone of C. fetus subsp. fetus that is unable to infect multiple host species (38).
MLST was compared to biochemical and established C. fetus band-based typing methods, AFLP analysis, PFGE, and PCR. Each of these methods measures genotypic characteristics which may change independently; hence, the subtypes that they define may not always be congruent (30). However, MLST subspecies differentiation results showed a high agreement with AFLP, Cf C05 PCR, and H2S results compared to the results of the PCR of Hum et al. (15) and the glycine test. In the present study, a strong correlation was observed between MLST and PFGE data, and both methods gave identical clusters of isolates from known C. fetus outbreaks. PFGE provided a higher level of discrimination than MLST. However, MLST measures variation that accumulates relatively slowly in the absence of selective pressure, and it may be more appropriate for studies of long-term C. fetus epidemiology and phylogenetics. PFGE indexes genome-wide variations and, despite the difficulties of interlaboratory comparisons, may be better suited for short-term investigations of C. fetus epidemiology, such as outbreak investigations, when such investigations are performed by a single laboratory. The level of discrimination of C. fetus isolates was higher by MLST than by the AFLP analysis method described, irrespective of the region of the AFLP pattern analyzed.
In summary, the observations that mammalian C. fetus isolates have a low level of genetic diversity and are genetically homogeneous compared to the homogeneities of other Campylobacter species suggest that C. fetus is young in evolutionary terms and may have arisen from a recent ancestor. In Fig. 1, the long branch length separating C. fetus from the other Campylobacter species suggests that its putative recent ancestor is missing from this tree. This is supported by the lack of any sequences characteristic of these five other Campylobacter species within the C. fetus housekeeping genes. The predominantly stepwise accumulation of point mutations in C. fetus indicates clonal evolution, in which genetic change accumulates slowly only by the vertical transmission of point mutations. The genetic stability of C. fetus is further confirmed by the wide geographic distribution of identical STs and the observation of congruence between the C. fetus subspecies, sap type, and ST. The two C. fetus subspecies were extremely closely genetically related, but ST-4 was associated only with C. fetus subsp. venerealis and was therefore, by definition, a "bovine" subclone. MLST confirms that mammalian C. fetus is genetically stable and clonal, provides a useful additional tool for C. fetus subspecies differentiation and epidemiology, and compares favorably to other bacterial typing methods.
| ACKNOWLEDGMENTS |
|---|
We are very grateful to D. Sanchez (Universidad Nacional de General San Martín, Instituto de Investigaciones Biotecnológicas, San Martín, Argentina) for comparing the C. jejuni pgm sequence with the sequence of C. fetus by use of the BLAST program and providing us with the C. fetus pgm sequence identified for primer design. Furthermore, we thank the following for kindly supplying C. fetus isolates: C. Campero, Instituto Nacional de Tecnoloíga Agropecuaria, Balcarce, Argentina; J. Godfroid, Veterinary and Agrochemical Research Centre, Brussels, Belgium; J. P. Butzler, University of Brussels, Brussels, Belgium; P. Idigoras, Hospital Donostia, San Sebastián, Spain; W. Kalka-Moll, University of Cologne, Germany; J. Varga, Faculty of Veterinary Science, Budapest, Hungary; S. Fujimoto, School of Health Sciences, Kyushu University, Kyushu, Japan; W. Ang, Erasmus University Medical Centre, Rotterdam, The Netherlands; M. Koene and I. Visser, Animal Health Service, Deventer, The Netherlands; M. Hinton, Onderstepoort Veterinary Institute, Onderstepoort, South Africa; L. Guler, Veteriner Kontrol ve Arastirma Enstitusu, Konyn, Turkey; J. Frost, Health Protection Agency, London, United Kingdom; J. Corry, University of Bristol, Bristol, United Kingdom; M. Toszeghy, Veterinary Laboratories Agency, Weybridge, United Kingdom; P. Fields, Centers for Disease Control and Prevention, Atlanta, GA; L. Schroeder-Tucker, National Veterinary Services Laboratories, Ames, IA; I. Wesley, National Animal Disease Centre, Ames, IA; and Z. Tu and M. Blaser, Vanderbilt University School of Medicine, Nashville, TN.
| FOOTNOTES |
|---|
M. van Bergen and K. Dingle contributed equally to the manuscript, and both are considered the first author. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Antimicrob. Agents Chemother. | Clin. Microbiol. Rev. |
|---|---|
| Clin. Vaccine Immunol. | ALL ASM JOURNALS |
|---|