Previous Article | Next Article ![]()
Journal of Clinical Microbiology, December 2008, p. 3957-3964, Vol. 46, No. 12
0095-1137/08/$08.00+0 doi:10.1128/JCM.01548-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

D. B. Hunter,1
V. R. Parreira,1
R. R. Kulkarni,1
Y.-F. Jiang,1
J. F. Prescott,1 and
P. Boerlin1,3*
Department of Pathobiology, Ontario Veterinary College, University of Guelph, Guelph, Ontario N1G 2W1, Canada,1 Nutreco Canada Agresearch, 473-Sixth Concession Road, R.R. #3, Burford, Ontario N0E 1A0, Canada,2 Laboratory for Foodborne Zoonoses, Public Health Agency of Canada, Guelph, Ontario N1G 3W4, Canada3
Received 11 August 2008/ Returned for modification 13 September 2008/ Accepted 10 October 2008
|
|
|---|
|
|
|---|
There are various molecular typing techniques used for assessing bacterial population diversity and for outbreak investigation of C. perfringens. These methods include, among others, pulsed-field gel electrophoresis (4, 7, 23), multilocus variable number of tandem repeats analysis (5, 25), and multilocus sequence typing (MLST) (16). MLST is a procedure for characterizing bacteria by using the sequences of several "house-keeping" genes (20). Generally, 400 to 500 base pairs of each gene are amplified by PCR and sequenced, and differences in the DNA sequences define the allelic profile for each isolate, known as the sequence type (ST). Genetic events such as point mutations or recombinations result in different alleles and subsequently different allele combinations (8, 32). Most bacterial species have sufficient variation within house-keeping genes to provide many alleles per locus, allowing multiple STs to be identified using a limited number of loci. An MLST protocol has been developed for C. perfringens by Jost and collaborators (16); however, it is not currently available as a standard protocol on MLST.net (http://www.mlst.net) or PubMLST (http://pubmlst.org). Although MLST is generally expensive, it can be set up as a high-throughput technique which generates data that can be shared between laboratories and allows for the creation of a global database of sequences. This global approach would be of great value for the study of NE epidemiology and of other diseases caused by C. perfringens.
In this study, our aim was to characterize the genetic diversity of C. perfringens isolates from outbreaks of NE in broiler chicken populations by using the MLST protocol described by Jost and collaborators (16). Our focus was on investigating the level of clustering of MLST types in NE outbreaks to assess the diversity of isolates within and between outbreaks. Others have shown that only a few pulsed-field gel electrophoresis types of C. perfringens can be isolated from field cases of NE, whereas a higher diversity can be found in the intestinal flora of healthy birds (7, 12, 23). With its population genetics implications, MLST data will help provide a clearer picture of the epidemiology of C. perfringens and of NE.
|
|
|---|
From each bird, cecal, mid-small-intestinal, and duodenal samples were removed aseptically, and the contents were suspended in 5 ml Difco brain heart infusion base (Becton Dickinson Microbiology Systems, Sparks, MD) plus 20% glycerol (Fisher Scientific, Fair Lawn, NJ) and frozen at –70°C. Duodenal samples were used in NE cases, and mid-small-intestinal and cecal samples were used only when C. perfringens isolates could not be recovered from the duodenum. A 100-µl aliquot of each sample was plated onto Shahidi-Ferguson perfringens selective medium plates (Becton Dickinson Microbiology Systems). Plates were incubated anaerobically for 24 h at 37°C. Two single colonies with dark centers were subcultured on Shahidi-Ferguson perfringens plates and subcultured again for purity. At least three isolates collected from sick or dead birds from each outbreak were used for MLST analysis. Whenever possible, isolates from the matched healthy birds were obtained from the duodenum, but cecal isolates were used if no duodenal isolates could be successfully cultured.
A total of 61 C. perfringens isolates were cultured from 45 broiler chickens, and these isolates were used for MLST analysis. These included 41 isolates from birds with NE, 9 isolates from matched healthy birds, and 11 isolates from nonoutbreak birds (Table 1). An additional C. perfringens isolate, CP4, which is a known virulent strain isolated from an NE case in Ontario (30), was also included for MLST analysis and toxin typing.
|
View this table: [in a new window] |
TABLE 1. Complete list of Clostridium perfringens isolates typed by MLST, including outbreak attributes, ST, bacitracin MIC determined by Etest, and presence or absence of netB, tpeL, and cpb2 genes detected by PCR
|
|
View this table: [in a new window] |
TABLE 2. MLST locus properties
|
eBURST analysis of isolates was performed, using eBURSTv3, with STs assigned based on unique allele profiles (http://eburst.mlst.net) (10). Groups were defined by isolates with at least seven of the nine alleles being identical. All ST numbers were based on this analysis (data not shown).
Presence of toxin genes. Detection of various toxin genes was carried out by PCR amplification of each C. perfringens isolate (n = 61). This included toxin typing (alpha, beta, epsilon, and iota toxin genes), which was performed as a multiplex reaction according to Yoo and collaborators (35). The beta-2 toxin gene (cpb2) and the enterotoxin gene (cpe) were detected using protocols described by Meer and Songer (21) and Herholz and collaborators (14), respectively. Additionally, a recently described NE beta-like toxin, NetB (netB) (17), was detected by the primers AKP78 and AKP79 described in the publication.
The presence of the tpeL (toxin C. perfringens large cytotoxin) gene was detected by PCR amplification of an internal fragment, using primers designed from the C. perfringens toxin type C sequence of tpeL (AB262081) supplied by Amimoto and collaborators (2). Amplifications were performed in a 25-µl total volume containing the following: 5 µl of template DNA, 1x PCR buffer with Mg2+ (New England BioLabs, Ipswich, MA), 0.2 mM deoxynucleoside triphosphate mixture, 2.5 units of Taq DNA polymerase (New England BioLabs), and 800 nM of each primer (T2875-F [ATAAGTGCCAATCAATATGA] and T3348-R [TTCTCCTGAATTAATATTAA]). After an initial denaturing step of 3 min at 94°C, the samples were subjected to 30 cycles of denaturation at 94°C for 1 min, annealing at 55°C for 1 min, extension at 72°C for 1 min, and final extension at 72°C for 10 min. All PCR products were visualized by horizontal agarose gel electrophoresis.
One netB and three tpeL PCR products were sequenced to confirm the specificity of the PCR. The entire tpeL gene from C. perfringens strain CP4 was amplified using primers specific to sequences 298 bp upstream of the first potential tpeL start site to 370 bp downstream of the stop codon for tpeL in C. perfringens type C (AB262081). Both strands of this entire PCR product were sequenced by primer walking, and the DNA sequence was deposited in the GenBank database (EU848493). A Southern blot to determine if tpeL is plasmid-borne was performed on plasmid preparations (plasmid midi kit; Qiagen Inc.) from four isolates (three tpeL positive and one negative), using a PCR digoxigenin-labeled probe (Roche Diagnostics, Mannheim, Germany) synthesized using the primers and PCR conditions described for tpeL detection by Amimoto and collaborators (2).
Susceptibility testing. Susceptibility to bacitracin was determined for each isolate by using Etest strips (AB Biodisk, Solna, Sweden) following the manufacturer's instructions. Isolates were grown overnight on blood agar plates, and a 1.0 McFarland standard was suspended in Mueller-Hinton broth (Sensititre, Westlake, OH). Brucella agar (150 mm) with vitamin K and Hemin agar plates (Medox Diagnostics, Ottawa, Ontario, Canada) were swabbed with the suspension, and single strips were applied to the surface. Readings were taken after 24 h of incubation under anaerobic conditions, and end points were determined using Etest guidelines for anaerobes. C. perfringens strain ATCC 13124 was used as a control. Isolates with MICs of >16 µg/ml were considered resistant, according to Johansson and collaborators (15).
Statistical analysis. Linkage disequilibrium analysis was performed using Arlequin version 3.11 (9) population genetics software. Polymorphic sites (n = 116) from each ST were analyzed (only one isolate per ST was used for this analysis), and linkage disequilibrium testing was estimated for all pairs of polymorphic sites. Significance was set at a P value of 0.05, with 10,000 steps in the Markov chain and 1,000 dememorization steps, resulting in a matrix of P values. A second measure of linkage disequilibrium was performed using the tool provided by MLST.net (http://linux.mlst.net/link_dis/index.htm). This method calculates an index of association (IA), which measures the degree of association between loci (28).
Nucleotide sequence accession numbers. Each unique sequence for each locus found in this study (n = 53) has been uploaded to the GenBank database (http://www.ncbi.nlm.nih.gov) under the following accession numbers: EU605734-EU605739 for plc, the alpha toxin gene; EU605740-EU605748 for ddlA, the D-alanine-D-alanine ligase gene; EU605749-EU605753 for dut, the deoxyuridine triphosphatase gene; EU605754-EU605759 for glpK, the glycerol kinase gene; EU605760-EU605764 for gmk, the deoxyguanylate kinase gene; EU605765-EU605768 for recA, the recombinase gene; EU605769-EU605773 for sod, the superoxide dismutase gene; EU605774-EU605782 for tpiA, the triose phosphate isomerase gene; and EU605783-EU605786 for pfoS, the putative membrane protein gene.
|
|
|---|
Amplification and sequence analysis of the 61 isolates were successful using the protocol outlined by Jost and collaborators (16). The addition of the pfoS locus did not create more STs and did not increase the discrimination of the method even though four different alleles were found (Table 2). The pfoS amplicon also showed some sequence length polymorphism of up to three base pairs, which was not observed in the other loci (Table 2). One pfoS allele (EU605786) showed a deletion within a long poly(A) sequence, introducing an internal stop codon, and resulted in a truncated PfoS protein.
The average number of alleles per locus was 5.9, with a minimum of 4 and a maximum of 9. Twenty-two STs were identified. With a few exceptions (outbreaks 4, 5, and 11), clustering of isolates within each outbreak was observed (Fig. 1). Isolates obtained from healthy birds in the same outbreak flock (outbreaks 5, 6, 7, and 8) also clustered with isolates from dead birds affected by NE. Multiple isolates obtained from a single bird were either identical (n = 8) or very similar (n = 2) (Fig. 1). The isolates from healthy birds originating from the seven farms of processor 2 all clustered tightly together. Although a major cluster contained isolates from eight different outbreaks from processor 1 (outbreaks 1, 2, 4, 6, 8, 9, 10, and 11), the isolates from the three remaining outbreaks were part of unrelated clusters. Neighbor-joining cluster analysis produced identical major clusters of isolates compared to the UPGMA clustering.
![]() View larger version (21K): [in a new window] |
FIG. 1. Dendrogram of concatenated MLST sequences of 61 avian isolates from Southern Ontario generated with pairwise similarities and UPGMA clustering. It includes eight previously published avian C. perfringens MLST sequences (16). Isolates from the study by Jost and collaborators are labeled JGS and shown in gray. CP4 is a virulent NE-associated isolate from another previous study in Ontario (30). The three C. perfringens strains with published genome sequences (strains 13, ATCC 13124, and SM101) were also included, and their MLST types were generated in silico (22, 27). Selected bootstrap values are indicated at nodes. Susceptibility to bacitracin was determined using a 16-µg/ml breakpoint. STs were determined through eBURSTv3 analysis (10), and ST01 and 18 are labeled in boldface type (founders).
|
The isolate JGS4135 from Jost and collaborators (16), which was the only one studied by these authors that clustered tightly with our own isolates, was traced back to a poultry sample from Ontario originally collected in Guelph (Fig. 1). Strain CP4, which was isolated about 10 years ago from a chicken with NE in Ontario, is also part of the major ST01 group (Fig. 1).
Results from the Arlequin linkage disequilibrium analysis (Fig. 2) showed a random distribution of significant values. The extent of significant disequilibrium (P
0.05) between polymorphic sites ranged between 7.2% (tpiA and ddlA) and 43.9% (glpK and ddlA), with an overall average of 23.6%. Linkage disequilibrium within each locus was higher (38.0% of polymorphic sites), as was expected. The IA value calculated using the MLST.net tool was 3.140, which implies significant linkage disequilibrium as well (28).
![]() View larger version (46K): [in a new window] |
FIG. 2. Linkage disequilibrium matrix of polymorphic sites in 22 STs, generated by Arlequin version 3.11 (9). Darker cells refer to more-significant P values (black cells represent values of 0.05 and the lightest-gray cells up to 0.20). Locus names are shown on each axis, in sequential order on the chromosome.
|
DNA sequencing of the gene amplified from CP4 showed that the tpeL gene consisted of 4,955 bp (EU848493). The translated amino acid sequence of TpeL in strain CP4 has 97% identity with the sequence of TpeL from a C. perfringens toxin type B strain (NZ_ABDV01000024), and both are 128 amino acids longer than that in the toxin type C strain MC18 (AB262081). The complete tpeL genes from two NE isolates, 23 and 28, were sequenced and found to have 100% identity with tpeL from strain CP4. Among the 61 isolates investigated here, tpeL was detected in 11 outbreak-associated isolates (Table 1) which included eight STs and four out of eleven outbreaks. All isolates from processor 2 were negative for tpeL (n = 13). Southern blot results showed that the tpeL gene was present on a large plasmid of identical size (
100 kb) in all three tpeL-positive isolates (13, 30, and CP4) examined (data not shown).
Etest values for susceptibility to bacitracin ranged between 1.5 and 256 µg/ml (Table 1). A very distinct bimodal distribution was present with MICs being either between 1.5 and 3 µg/ml (susceptible) or >256 µg/ml (resistant). Seventeen resistant and 33 susceptible isolates were obtained from processor 1, whereas only bacitracin-resistant isolates (n = 11) were detected in birds from processor 2. All the bacitracin-resistant isolates were part of the two major clusters containing ST01 and ST18, but not all the isolates in these two clusters were resistant (Fig. 1).
|
|
|---|
The MLST method developed by Jost and collaborators was originally developed for a broad array of C. perfringens isolates of very diverse origins (16) and may not have been sufficient for this study, where only chicken isolates were used. Therefore, we initially attempted to increase the discriminatory power of this typing scheme by investigating potential additional loci. The loci for gyrA and rplL (24) and rpoD (primers derived from the published C. perfringens SM101 sequence [22]) showed little or no polymorphism at all in a subset of our isolates (data not shown). However, the pfoS locus did show some diversity and was added to the established protocol. Ultimately, this locus did not contribute any additional discrimination despite the presence of four different alleles at this locus in our isolates. It may also not be a suitable MLST locus because of its sequence length variability; one allele of the pfoS locus appeared to have resulted from the deletion of an "ATT" tandem repeat in the gene. Another pfoS allele encodes a truncated and probably nonfunctional protein. Therefore, it appears that the protocol described by Jost and collaborators (16) already provides a good estimate of C. perfringens diversity and was not improved with the addition of pfoS. When isolates from the three sequenced genomes of C. perfringens were included (Fig. 1), as well as the eight avian isolates from the Jost et al. study (16), these isolates generally belonged to different sequence types than to those from this study. Interestingly though, the JGS4135 isolate typed by Jost and collaborators was identical to isolate 10 in our study, and it was later revealed that the former isolate had originated from the University of Guelph about 10 years earlier. This is a convincing indication that this MLST system is epidemiologically viable and also showed that the same C. perfringens clone may have persisted for many years in Southern Ontario.
In previous studies (7, 12, 23), researchers have proposed that there are distinct clones within C. perfringens populations which cause NE outbreaks and that birds within an outbreak are usually affected by only one of these clones. Our results, which showed strong genetic clustering in isolates from NE outbreaks, both within a single bird and between birds of the same outbreak, are in agreement with this hypothesis. However, we did not observe any major differences between isolates from healthy birds and those from NE cases within any outbreak. This stresses once more the fact that NE is a multifactorial disease and that the mere presence of a specific C. perfringens strain in a bird is not a sufficient cause for the development of the disease.
The very strong clustering of isolates originating from birds of seven different farms from geographically separate areas affiliated with processor 2 is striking and was unexpected. Similarly, despite the higher diversity observed among isolates from farms with NE outbreaks and those associated with processor 1, the majority of these isolates belonged to a single major cluster (cluster containing ST01, -02, -10, -15, and -16) (Fig. 1). In the case of this processor, the clustering may be related to a higher propensity of a specific clone for causing NE. These clusters also suggest a limited degree of genetic diversity of C. perfringens throughout each processor's poultry operations. This could be the result of the dissemination of clones throughout integrated producer networks, starting at hatcheries and chick distribution systems. Selection of a limited number of resistant clones through antimicrobial agent use, either for prevention (cluster containing ST18) or therapy (cluster that contains ST01 and is associated with NE outbreaks), may also be part of the reason for the observed clustering. This latter hypothesis is supported by the fact that bacitracin resistance was clearly associated with these two major clusters and lineages.
Locus linkage was found to be high in the isolates characterized in this study, with significant linkage disequilibrium found between all loci (IA = 3.140); panmictic populations have an IA value approaching zero. This also shows that recombination rates for the core genome of C. perfringens are low or possibly absent (28). Additionally, only 22 STs were found in the sample population. When considering the number of possible allele combinations (over 5 million), the 22 STs found in 61 isolates may be considered low. These data suggest that the C. perfringens population under investigation had a highly clonal structure (28).
Since MLST examines core genes encoding house-keeping functions not usually involved in pathogenesis, there is not necessarily a correlation between virulence potential and the overall genomic relationships evidenced by MLST (20, 31). It is therefore possible that more virulent clones do exist within these outbreaks, which carry accessory genes and virulence factors not detected using MLST.
The complete absence of cpb2 in the NE outbreak isolates was surprising, since this toxin had originally been associated with virulence (11, 13, 34). However, there is a growing body of evidence that suggests that the beta-2 toxin is not a critical factor of C. perfringens infection in this particular case of poultry (4, 6, 30). The recently discovered netB toxin gene, however, was found exclusively in isolates which were all associated with NE outbreaks and was never detected in healthy control samples from processor 2. Although netB was detected in every outbreak, it was still found in healthy matched birds from outbreak flocks as well. Additionally, netB was not always detected in isolates from birds clearly affected with NE (Table 1, isolates 03 and 33). Overall, this study provided a good opportunity to observe the prevalence of the netB gene in a field situation where NE- and non-NE-associated isolates were used. The data suggest that this toxin is an important factor in C. perfringens infection, facilitating the development of NE. However, our results also show that it is neither a necessary nor a sufficient cause for the disease to occur. As the authors of the netB discovery stated (17), there may be many other unidentified factors which are also important for virulence.
The C. perfringens toxin TpeL (2) was identified in the culture supernatant of C. perfringens strain CP4, when the secreted proteins of CP4 were electrophoretically separated and reacted with the serum from broiler chickens immune to NE (18, 19), and was further identified by mass spectrometry (unpublished data). This is the first report describing the presence of tpeL genes in C. perfringens type A strains and the first report of the presence of this gene in isolates from chickens with NE. The tpeL gene identified in these toxin type A isolates is very similar to the tpeL gene in the toxin type B strain (NZ_ABDV01000024). The toxin type C strain (2), however, has a cytosine insert at position 4941, which creates a premature stop codon and thus truncates the protein. The variable presence of this toxin gene within clones and even within clonally related isolates from a single farm suggested that it is located on a mobile genetic element. Our Southern blot results confirmed that it is indeed located on large plasmids. The identification of TpeL among the secreted proteins of CP4 may be related to the greater virulence of this strain in an experimental model of NE in comparison to some other strains in which NE originated (30). Since tpeL was detected in isolates and STs associated with NE outbreaks only, it is reasonable to speculate that it contributes to the pathogenesis of NE, but further work is required to assess its contribution to the virulence of NE isolates.
In conclusion, the MLST protocol described by Jost and collaborators (16) was successfully applied in this study for typing avian C. perfringens isolates. The resulting data suggest that clustering of MLST types occurs in association with the source of the isolates and with bacitracin resistance. A major clonal lineage containing isolates from eight different outbreaks was identified, which may be either more prone to cause outbreaks than others or selected by the use of antimicrobial agents. Toxin typing suggests that NetB is a far more significant toxin than the beta-2 toxin, although there may be many more unexplored factors that are important for NE development, such as TpeL, which was identified here for the first time in NE isolates.
We thank the Animal Health Laboratory at the University of Guelph, the poultry processors, and the Canadian Food Inspection Agency for their assistance in this project. We also thank Teresa Crease for her help with the linkage disequilibrium analysis.
Published ahead of print on 22 October 2008. ![]()
Present address: Department of Agriculture, Food and Nutritional Science, University of Alberta, Edmonton, Alberta T6G 2P5, Canada. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»