- Research article
Multilocus sequence typing (and phylogenetic analysis) of Campylobacter jejuni and Campylobacter coli strains isolated from clinical cases in Greece
BMC Research Notesvolume 6, Article number: 359 (2013)
The molecular epidemiology of C. jejuni and C. coli clinical strains isolated from children with gastroenteritis, was investigated using the multilocus sequence typing method (MLST). This analysis establishes for the first time in Greece and constitutes an important tool for the epidemiological surveillance and control of Campylobacter infection in our country.
The MLST genotypes were compared with those gained by other typing methods (HS-typing, PFGE and FlaA typing) and were also phylogenetically analyzed, in order to uncover genetic relationships.
Among 68 C. jejuni strains, 41 different MLST-Sequence Types (MLST-STs) were found. Fifty six strains or 34 MLST-STs could be sorted into 15 different MLST-Sequence Type Complexes (MLST-STCs), while twelve strains or seven MLST-STs did not match any of the MLST-STCs of the database. Twenty C. coli strains belonged to 14 different MLST-STs. Eleven MLST-STs were classified in the same MLST-STC (828), and three were unclassifiable. There was no significant association between the MLST-STs and the results of the other typing methods.
Phylogenetic analysis revealed that some strains, classified to the species of C. jejuni, formed a separate, phylogenetically distinct group. In eight strains some alleles belonging to the taxonomic cluster of C. jejuni, were also detected in C. coli and vice versa, a phenomenon caused by the genetic mosaic encountered inside the genus Campylobacter.
The MLST-ST determination proved to be a very useful tool for the typing as well as the identification of Campylobacter on the species level.
Campylobacter is one of the leading causes of human bacterial gastroenteritis worldwide, with a very wide distribution in food animals. The genus Campylobacter comprises many species, among which the most commonly isolated in humans are Campylobacter jejuni and C. coli. Sporadic cases have been associated with consumption of undercooked poultry meat, while small or large outbreaks have been associated with raw milk and contaminated water . Most cases of human campylobacteriosis are self limiting, requiring no medical treatment. Therapeutic intervention is often warranted in children or in severe or long-lasting infections, the latter occurring mainly in immunocompromised individuals .
Due to the impact of Campylobacter on public health, epidemiological investigations analyzing the clonality of the isolated strains are very important, in order to trace the sources and routes of transmission, to follow up the temporal and geographic distribution of important phenotypic characteristics and to develop effective strategies for the control and prevention of the pathogen spread, especially inside the food chain .
Therefore, numerous phenotyping and genotyping methods are available . Although for the choice of the appropriate method, the rapidity, the cost and the ease of implementation are essential, the main point should be above all the ability in providing accurate, reliable and highly reproducible results . One method meeting the latter condition very well is the multilocus sequence typing (MLST), which has already been developed and used for the molecular typing of Campylobacter spp., based on the sequence identification of multiple genetic loci, located in seven housekeeping genes: aspA, glnA, gltA, glyA, tkt, pgm, uncA[7–9]. It has proved to be a useful tool for discriminating isolates, for the definition of the structures of the various Campylobacter populations and for disclosing associations among sequence types or lineages with the isolation sources or specific virulence factors, thus contributing to the identification of outbreaks and to the appropriate public health intervention . Results from various reports using MLST are imported in international databases, which are open to the enrichment by newly discovered MLST-types, providing an on-going updated global standpoint of Campylobacter epidemiology.
The present study aimed to investigate for the first time in Greece the molecular epidemiology of C. jejuni and C. coli clinical strains isolated from children with gastroenteritis, using the MLST. The bacterial strains included in the study were heat- stable serotyped (HS), genotyped by pulsed field gel electrophoresis (PFGE) and by flagellin A typing (flaA typing). The MLST genotypes were compared with those gained by the other three typing methods (HS-typing, PFGE and FlaA typing) and statistically analyzed. In addition the MLST types were phylogenetically identified on the species level and analyzed, in order to uncover genetic relationships [7, 11, 12].
A total of 68 C. jejuni and 20 C. coli clinical isolates were collected from infected children up to 14 years from 5 general hospitals of the area of Attica. All the strains were isolated from sporadic cases with gastroenteritis. Bacteria were identified by means of conventional bacteriological standard procedures. Stool samples were cultured for three days at 42°C in Campy enrichment Thioglycolate Medium, supplemented with trimethoprim, vancomycin, polymyxin B, cephalothin & ampotericin B [(Thioglycolate Medium w/0.16% Agar), (REMEL Inc. Kansas U.S.A.)]. Liquid cultures were centrifuged, the supernatants were discarded and the sediments were inoculated in Skirrow agar (Brucella Agar of w/5% Sheep Blood, Trimethoprim, Vancomycin, Polymyxin B, REMEL Inc. Kansas U.S.A.) and incubated for 48 h at 42°C, in a microaerophilic atmosphere for Campylobacter. Identification of typical colony was performed by hippurate hydrolysis and a commercially available identification kit (Api Campy, bioMerieux). Each strain was numerically designated with an internal number, corresponding to the respective clinical specimen sent to the laboratory. Ethical approval was not required under Greek research regulations and all samples were taken as part of standard patient care.
Serotyping was performed as reported previously  with a commercially available set of antisera (Campylobacter Antisera, “Seiken Set”, Denka Seiken, Japan) based on Penner’s heat-stable (HS) serogroups, and containing 25 absorbed single or group antisera against the most common HS-serotypes. HS-antigen extracts were prepared by the original Penner boiling-method.
Pulsed-field gel electrophoresis
PFGE was performed as reported previously .
Flagellin gene typing
Bacterial DNA isolation for MLST
Template DNA for the MLST was extracted with the commercially available extraction reagent “instagene matrix” (Bio Rad Laboratories, California, U.S.A.), according to the manufacturer’s instructions.
Multilocus sequence typing
For the PCR, the commercial ready mix “Jumpstart Red Taq” (SIGMA Laboratories, St. Louis, U.S.A.) was used. The final mix contained 25 μl Jumpstart Red Taq, 1 μM of each primer, 18 μl distilled water and 5 μl template. Primers were provided by Sigma Genosys St. Louis USA in lyophilized form. For each PCR reaction for the detection of each of the seven housekeeping genes (aspA, glnA, gltA, glyA, tkt, pgm, uncA) the primer concentration was 10 pmol/μl. The composition of the primer pairs used for the initial PCR for the detection of the housekeeping genes, as well as for the subsequent PCR for the sequencing in order to determine the alleles, has been described previously elsewhere [7, 9].
An aliquot of initial PCR sample was analyzed by conventional 1.2% agarose gel electrophoresis. DNA was stained with ethidium bromide and visualized and photographed with a computerized transluminator system (Chemidoc, Bio Rad Laboratories, California, U.S.A.).
The sequencing of the PCR products was performed by Macrogen Incorporation (908 World Merdian Center #60-24 Gasan-dong, Geumchum-gu Seoul, Republic of Korea).
Due to the existence of many independent variables (serotypes and genotypes) and a relative small number of observations, the penalized logistic regression with the LASSO-Least Absolute Shrinkage and Selection Estimator was used.
The electrophoretic patterns gained by the PFGE were expressed with a binary system and symbolized as follows: “A” for the presence of a band and “T” for the absence. The distances between the nucleotide clusters were measured after the model of Jukes and Cantor and the dendrograms (phylogenetic trees) were constructed with the UPGMA (Unweighted Pair Group Method with Arithmetic Mean) method.
For the MLST results, the alignment of all the gene sequences and loci was performed with the algorithm Clustal W using the MEGA (Molecular Evolutionary Genetics Analysis) 5.05 software. The genetic distances between the sequences were calculated using the Kimura 2-parameter model . The phylogenetic trees were constructed using the neighbor – joining method and their reliability was checked by bootstrapping analysis (1000 replicates). One cluster has considered as significant, if it was present in over than 75% of the permuted trees.
For C. jejuni 22 out of the 68 strains were not typed, while 20 different serotypes or serotype complexes were found in the 46 serotyped strains. From the 20 C. coli strains 12 were not typed and 7 different serotypes were detected in the 8 serotyped strains. Some strains gave positive agglutination reactions with more than one antisera (Tables 1 and 2). There was no significant association of a certain serotype with a certain MLST type for both Campylobacter species. All the serotypes of C. coli were shared also by C. jejuni, while the serotype complex HS 37 + HS 57 was found only in C. coli.
Flagellin gene typing
Among the 68 C. jejuni and the 20 C. coli strains, 28 and 10 different flaA types, were detected respectively (Tables 1 and 2). The two species bore five flaA types in common: 7, 10, 23, 27 and 39. There was no significant association between a certain flaA type and a certain MLST type.
Pulsed-field gel electrophoresis
Multilocus sequence typing
All the strains were designated to a certain MLST sequence type (MLST-ST) by the combination of the seven allelic housekeeping genes. Moreover, some of them could be classified to a certain MLST sequence type complex (MLST-STC), which represented a clonal complex of a specific combination of four or more allelic genes . Each MLST-ST and MLST-STC is marked with a unique number and the corresponding gene combinations are registered and can be accessed in an international electronic database PubMLST (http://pubmlst.org/campylobacter/).
Among the 68 C. jejuni strains, 41 different MLST-STs were found. Fifty six strains or 34 MLST-STs could be sorted into 15 different MLST-STCs, while the remaining 12 strains or seven MLST-STs did not match any of the MLST-STCs of the MLST database.
The 20 C. coli strains belonged to 14 different MLST-STs. Eleven MLST-STs were classified in the same MLST-STC-828 were is common as described previously , while the remaining three could not match any MLST-STC.
The numbers of different alleles for each housekeeping gene found in C. jejuni were as follows: 13 for aspA, 14 for glnA, 10 for gltA, 10 for glyA, 16 for pgm, 14 for tkt and 9 for uncA. For C. coli the corresponding numbers were: 1 for aspA, 3 for glnA, 4 for gltA, 5 for glyA, 6 for pgm, 8 for tkt and 5 for uncA. The distribution of C. jejuni and C. coli MLST- STs and STCs from Greece, among strains isolated from various countries (data available on PubMLST) is presented in Tables 3 and 4. There was no significant correlation between the PFGE and the MLST profiles of both Campylobacter species.
It must be pointed out, that there was a complete agreement of the MLST-STs of the database, with the conventionally performed identification of the strains on the species level (C. jejuni or C. coli).
Phylogenetic analysis aimed to classify the various genotypes by constructing phylogenetic trees based on the patterns gained by PFGE and MLST (Figures 1 and 2). PFGE results are representative for the entire genome, but the method is biased by selective forces, while the MLST patterns are suitable for phylogenesis and global epidemiology as they derive from the combination of the allelic forms of the seven housekeeping genes. For that purpose the electrophoretic patterns were converted to binary numbers. The results arising from both analytical methods were in complete agreement with each other and the strains were grouped together by species (C.jejuni and C. coli) in two distinct clusters. Phylogenetic analysis of the MLST patterns revealed that some strains, although classified to the species of C. jejuni, they formed a separate, phylogenetically distinct group, a point of epidemiological interest (Figure 1). The paradox arising from the MLST phylogenetic patterns, was that in eight strains some alleles belonging to the taxonomic cluster of C. jejuni, were also detected in C. coli and vice versa, a phenomenon already described in previous reports [7, 11, 12], caused by the genetic mosaic encountered inside the genus Campylobacter. The genetic mosaic of these eight particular strains is shown in Figure 3. Regarding the number of mosaic alleles, for C. jejuni there was only one allele recorded in the uncA locus and was specific for ST-61, and for C. coli, there were three alleles in the tkt locus and one in each of the uncA, pgm and gltA loci. The distribution of mosaic genes was asymmetrical between C. jejuni and C. coli, and only one of the mosaic genes showed horizontal transfer from C. coli to C. jejuni STs, while four mosaic genes from C. jejuni into C. coli.
The results confirm the high discriminatory power of the MLST analysis for Campylobacter genotyping, which is applicable worldwide [6, 19–21]. The advantage of MLST over the other described genotyping techniques arises from its standardization by the use of seven housekeeping genes and the subsequent sequence analysis of the products, warranting very reproducible and reliable results. Moreover the MLST-ST determination proved to be a very useful tool for the identification of Campylobacter on the species level. This is not applicable regarding the MLST-STC determination. C. jejuni and C. coli share approximately 86.5% identity at the nucleotide sequence level within the MLST loci . Therefore strains of these two species can belong to the same MLST-STC but not to the same MLST-ST. The comparison of results obtained by other genotyping techniques is often difficult, due to the diversity of the used performance protocols among the various laboratories. In the present study the MLST analysis yielded more genotypes than the other two methods (PFGE and flaA typing).
Among the 88 Campylobacter strains included in the study, 55 different MLST genotypes were found and some of them are reported for the first time (Tables 3 and 4). The absence of prevalence of a particular genotype is consistent with the aspect, that there is considerable diversity among the Campylobacter strains isolated in the area of Attica . As far as we know, there is no report from Greece about campylobacteriosis outbreaks during the last thirty years, all the clinical cases being sporadic [2, 13, 14, 22]. This observation is enforced in our investigation by the phenotypic and genotypic diversity found by the other typing methods, as well by the absence of association of any specific phenotype or genotype (serotype, PFGE or flaA) with a certain MLST type. In addition, the phylogenetic analysis provided a very fine discrimination with respect to the identification and the clonal distribution of the strains, revealing also mosaic genotypes derived through interspecies recombination of the housekeeping genes between C. jejuni and C. coli strains. These findings represent an asymmetric gene flow, (4.69% in C. jejuni vs 25% in C. coli), also described elsewhere [11, 12].
The MLST typing as well as the phylogenetic analysis of Campylobacter clinical strains isolated from Greece is presented for the first time in the present study and significantly completes the epidemiological database kept in our laboratory. This database includes, besides the already described phenotypic and genotypic features, all the demographic and clinical data of the patients and constitutes an important tool for the epidemiological surveillance and control of Campylobacter infection in our country.
- FlaA typing:
Flagellin A typing
Molecular evolutionary genetics analysis
Multilocus sequence typing
Multilocus sequence typing-sequence type
Multilocus sequence typing-sequence type complex
Unweighted pair group method with arithmetic mean.
Man SM: The clinical importance of emerging Campylobacter species. Nat Rev Gastroenterol Hepatol. 2011, 8: 669-685. 10.1038/nrgastro.2011.191.
Ioannidis A, Nicolaou C, Chatzipanagiotou S: Correlation between flagellin A (flaA) genotypes and antimicrobial susceptibility patterns of Campylobacter jejuni strains isolated from children with gastroenteritis in Athens, Greece. Mol Diagn Ther. 2009, 13: 389-395. 10.1007/BF03256345.
Barton Behravesh C, Jones TF, Vugia DJ, Long C, Marcus R, Smith K, Thomas S, Zansky S, Fullerton KE, Henao OL, Scallan E: Deaths associated with bacterial pathogens transmitted commonly through food: foodborne diseases active surveillance network (FoodNet), 1996–2005. J Infect Dis. 2011, 204: 263-267. 10.1093/infdis/jir263.
Sheppard SK, Colles F, Richardson J, Cody AJ, Elson R, Lawson A, Brick G, Meldrum R, Little CL, Owen RJ, et al: Host association of Campylobacter genotypes transcends geographic variation. Appl Environ Microbiol. 2010, 76: 5269-5277. 10.1128/AEM.00124-10.
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.
Maiden MC, Bygraves JA, Feil E, Morelli G, Russell JE, Urwin R, Zhang Q, Zhou J, Zurth K, Caugant DA, et al: Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci USA. 1998, 95: 3140-3145. 10.1073/pnas.95.6.3140.
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.
Miller WG, On SL, Wang G, Fontanoz S, Lastovica AJ, Mandrell RE: Extended multilocus sequence typing system for Campylobacter coli, C. lari, C. upsaliensis, and C. helveticus. J Clin Microbiol. 2005, 43: 2315-2329. 10.1128/JCM.43.5.2315-2329.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.
Quintana-Hayashi MP, Thakur S: Phylogenetic analysis reveals common antimicrobial resistant campylobacter coli population in antimicrobial-free (ABF) and commercial swine systems. PLoS One. 2012, 7: e44662-10.1371/journal.pone.0044662.
Sheppard SK, McCarthy ND, Jolley KA, Maiden MC: Introgression in the genus Campylobacter: generation and spread of mosaic alleles. Microbiology. 2011, 157: 1066-1074. 10.1099/mic.0.045153-0.
Sheppard SK, McCarthy ND, Falush D, Maiden MC: Convergence of Campylobacter species: implications for bacterial evolution. Science. 2008, 320: 237-239. 10.1126/science.1155532.
Chatzipanagiotou S, Papavasileiou E, Lakumenta A, Makri A, Nicolaou C, Chantzis K, Manganas S, Legakis N: Heat-stable antigen serotyping of Campylobacter jejuni strains isolated from hospitalized children in Athens, Greece. Eur J Epidemiol. 2003, 18: 1097-1100. 10.1023/A:1026108702971.
Ioannidis A, Nicolaou C, Legakis NJ, Ioannidou V, Papavasileiou E, Voyatzi A, Chatzipanagiotou S: Genotyping of human Campylobacter jejuni isolates in Greece by pulsed-field gel electrophoresis. Mol Diagn Ther. 2006, 10: 391-396. 10.1007/BF03256216.
Nachamkin I, Bohachick K, Patton CM: Flagellin gene typing of Campylobacter jejuni by restriction fragment length polymorphism analysis. J Clin Microbiol. 1993, 31: 1531-1536.
Wassenaar TM, Newell DG: Genotyping of Campylobacter spp. Appl Environ Microbiol. 2000, 66: 1-9. 10.1128/AEM.66.1.1-9.2000.
Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.
Sheppard SK, Didelot X, Jolley KA, Darling AE, Pascoe B, Meric G, Kelly DJ, Cody A, Colles FM, Strachan NJ, et al: Progressive genome-wide introgression in agricultural Campylobacter coli. Mol Ecol. 2013, 22: 1051-1064. 10.1111/mec.12162.
Urwin R, Maiden MC: Multi-locus sequence typing: a tool for global epidemiology. Trends Microbiol. 2003, 11: 479-487. 10.1016/j.tim.2003.08.006.
Kinana AD, Cardinale E, Tall F, Bahsoun I, Sire JM, Garin B, Breurec S, Boye CS, Perrier-Gros-Claude JD: Genetic diversity and quinolone resistance in Campylobacter jejuni isolates from poultry in Senegal. Appl Environ Microbiol. 2006, 72: 3309-3313. 10.1128/AEM.72.5.3309-3313.2006.
Meinersmann RJ, Berrang ME, Little E: Campylobacter spp. recovered from the upper Oconee River Watershed, Georgia in a 4-year study. Microb Ecol. 2012, 65 (1): 22-27.
Ioannidis A, Nicolaou C, Legakis NJ, Ioannidou V, Chatzipanagiotou S: The first database comprised of flagellin gene (flaA) types of Campylobacter jejuni human clinical isolates from Greece. Eur J Epidemiol. 2006, 21: 823-829. 10.1007/s10654-006-9073-2.
The authors declared that they have no competing interests.
VI performed the experimental part, AI and EM worked out the methodology, PB performed the statistical analysis, CN, NL and SC designed the study, and SC wrote the manuscript. All authors read and approved the final manuscript.
Vasiliki Ioannidou, Anastasios Ioannidis, Emmanouil Magiorkinis, Pantelis Bagos, Chryssoula Nicolaou, Nicolaos Legakis contributed equally to this work.