Previous Article | Next Article ![]()
Applied and Environmental Microbiology, May 2008, p. 2778-2786, Vol. 74, No. 9
0099-2240/08/$08.00+0 doi:10.1128/AEM.02828-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Department of Bacteriology, Food Research Institute, University of Wisconsin—Madison, Madison, Wisconsin 53706,1 Enteric Diseases Laboratory Branch, Centers for Disease Control and Prevention, Atlanta, Georgia 303332
Received 14 December 2007/ Accepted 21 February 2008
|
|
|---|
|
|
|---|
BoNTs can be immunologically distinguished using homologous antitoxins into seven serotypes, designated A to G. Among these serotypes there is considerable genetic variation, as demonstrated by the recognition of at least 47 subtypes identified thus far (15, 28). Previous studies have determined that the BoNT genes of all strains of C. botulinum and neurotoxigenic strains of Clostridium butyricum and Clostridium baratii have a set of genes located upstream of the BoNT gene that are organized as neurotoxin clusters (4, 7, 10, 11, 13, 14, 17, 18, 19, 23).
Available data indicate that there exist two primary types of neurotoxin clusters: (i) a hemagglutinin (HA) cluster and (ii) an Orfx cluster with genes (orf genes) of unknown function. The HA cluster consists of genes encoding HA17, HA33, HA70, BotR, and NTNH. The Orfx cluster consists of genes encoding Orfx3, Orfx2, Orfx1, P47, P21, and NTNH. The known serotypes associated with HA clusters are BoNT/A1, BoNT/B, BoNT/C, BoNT/D, and BoNT/G, while serotypes BoNT/E, BoNT/F, and BoNT/A2 are associated with Orfx clusters; the E serotype cluster is unique as it lacks p21 (4, 10, 13, 14, 19, 25, 27). However, results from previous studies indicate that different BoNT/A subtypes possess either an HA cluster or an Orfx cluster associated with the expressed bont gene, depending on the subtype and strain (7, 8, 11, 12).
It has been shown that the HA cluster proteins can form toxin complexes of various sizes including LL toxin (
900 kDa), L toxin (
500 kDa), and M toxin (
300 kDa), comprised of various combinations of HA proteins, NTNH, and BoNT (
150 kDa) (17, 18, 23). The biological and structural roles of the complex proteins are not completely characterized, although it is believed that they serve the role of protecting BoNT from digestive enzymes and adverse conditions present in the gastrointestinal tract (23). The function of orfx genes and their role in the neurotoxin cluster are still unknown, but it has been postulated that they perform a function similar to their HA cluster counterparts because the location and orientation of the Orfx and HA cluster genes relative to the BoNT gene in the neurotoxin cluster are analogous (7). A regulatory gene called botR (21, 22) has also been identified in the HA cluster. botR encodes a sigma 70 factor which positively regulates the expression of BoNT and the HA proteins (9). The Orfx cluster has an analogous gene, p21, which appears to be a sigma 70 factor, as determined by amino acid homology BLAST analysis.
BoNT/A is of particular importance and interest since it causes the most severe form of human botulism, is the most significant threat in bioterrorism, and has been increasingly used as a pharmaceutical modality (1, 6). During the past 5 years, a focus of our laboratories has been the study of C. botulinum strains that produce novel BoNT/A subtypes including an A3 strain producing BoNT/A3, an A4 strain that is a bivalent strain producing both BoNT/A4 and BoNT/B (2), and a strain (termed HA– Orfx+ A1) that produces BoNT/A1 with a neurotoxin cluster arrangement normally associated with the BoNT/A2. This paper presents a comprehensive comparison of the individual neurotoxin cluster genes and bont genes among different C. botulinum subtype A strains on both nucleotide and amino acid levels.
In order to study the neurotoxin cluster genes in these subtypes, we have identified and sequenced each of these strains' neurotoxin cluster genes and have compared the results on nucleotide and amino acid levels to the known sequences from ATCC 3502 (a BoNT/A1 producer) (7, 26), Kyoto F (a BoNT/A2 producer) (7), and NCTC 2916 (a BoNT/A1 producer with a unexpressed BoNT/B gene) (10, 16, 24) to determine their genetic relationships and possible genetic lineages. Our laboratory showed that the genes encoding BoNT/A3 and BoNT/A4 are located on large plasmids; this is the first demonstration that genes encoding BoNT/A are present on extrachromosomal elements (20). Prior to our discovery in these strains and the subsequent confirmation of our findings (27), the gene encoding BoNT/A1 was expected to be located on the chromosome by analysis of the genome of C. botulinum ATCC 3502 (26).
Previous studies of BoNT clusters have emphasized the arrangement of the genes and, particularly, the amino acid sequences of the proteins; here, we also address the nucleotide sequence relationships and possible mechanisms of acquisition and rearrangement of the BoNT and neurotoxin cluster genes of the A3 and A4 subtypes. In addition, we describe a novel neurotoxin cluster, HA– Orfx+ A1, in the present study and discuss possible mechanisms of toxin complex acquisition and function.
|
|
|---|
Total genomic DNA isolation.
Total genomic DNA was isolated from C. botulinum by lysozyme and proteinase K treatment as described previously (7). DNA was then diluted to a concentration of 50 ng/µl and used for PCR amplification.
PCR amplification and DNA sequencing.
PCR amplifications were performed using a GeneAmp High Fidelity PCR System (Applied BioSystems). PCR conditions were as follows: 95°C for 2 min, followed by 25 cycles of 95°C for 1 min, annealing for 45 s at 48°C, and extension at 72°C, with a final extension at 72°C for 10 min. Following amplification, PCR products were isolated using a PureLink PCR purification kit (Invitrogen). Sequencing setups were produced under conditions recommended by the University of Wisconsin Biotechnology Center using an ABI Prism BigDye Cycle Sequencing Kit (Applied BioSystems). Primers used for PCR and sequencing are listed in Table 1. PCRs were performed in a staggered manner such that the amplicons produced overlapping products for each of the genes in each neurotoxin cluster. Appropriate primers were then used for sequencing each PCR product. Correct assembly of the contigs was verified by the overlapping sequence data. The toxin genes were sequenced as necessary to confirm the suspected neurotoxin cluster arrangements for each bacterial type. Sequencing analysis was performed at the University of Wisconsin Biotechnology Center, and final sequencing results were analyzed using the Vector NTI Suite program (Invitrogen).
|
View this table: [in a new window] |
TABLE 1. Primers used to identify and characterize the clusters present in the various subtypes of C. botulinum studied
|
|
|
|---|
![]() View larger version (19K): [in a new window] |
FIG. 1. Diagram of the six known cluster arrangements for C. botulinum strains producing BoNT/A neurotoxin. Gene location and direction are indicated by arrows with gene nomenclature given beneath each one. Strain identifications are given to the right of the series of genes in each cluster.
|
|
View this table: [in a new window] |
TABLE 2. Nucleotide identity among the genes of the C. botulinum HA clustera
|
|
View this table: [in a new window] |
TABLE 3. Amino acid similarity and identity for the proteins of the C. botulinum HA clustera
|
Comparison of orfx and p47 genes.
Analysis of the Orfx cluster genes of the A3, A4 and HA– Orfx+ A1 strains gave more diverse results. Comparisons of the nucleotide sequences of the orfx3 genes demonstrated that the highest degree of similarity was between the HA– Orfx+ A1 and the NCTC 2916 strains. orfx3 in strain A3 had the highest degree of divergence with 91.9% identity compared to the gene in HA– Orfx+ A1 (Table 4). On the amino acid level, the relationships among the various strains were consistent with the nucleotide relationship except that the difference between A3 and HA– Orfx+ A1 was less pronounced as many of the nucleotide alterations resulted in conserved amino acid changes (Table 5).
|
View this table: [in a new window] |
TABLE 4. Nucleotide identity (%) between the genes of the C. botulinum Orfx clustera
|
|
View this table: [in a new window] |
TABLE 5. Amino acid similarity and identity for the proteins of the C. botulinum Orfx clustera
|
Comparison of the orfx1 genes gave a different order of relatedness, whereby the A3 strain was most similar to A2, followed by HA– Orfx+ A1, and then A4 (Table 4). At the nucleotide level, the orfx1 gene from A3 and HA– Orfx+ A1 were closely related to that of A2, having 97.4% and 94.4% identity, respectively, while having 92.3% identity to each other (Table 4). However, at the amino acid level, the differences between the A3 and HA– Orfx+ A1 strains compared to A2 were substantial. While A3 was 97.9% and 97.2% similar and identical to A2, HA– Orfx+ A1 showed 91.5% similarity and 90.1% identity to A2, lending slightly more identity than A4 (Table 5).
The relatedness of the p47 genes was interesting as it is the only gene of the Orfx cluster in which HA– Orfx+ A1 strain had 100% nucleotide identity with p47 of A2 (Kyoto F). Based on nucleotide identity, the order of relatedness to A2 was HA– Orfx+ A1, A3, and then A4 with only 80.4% identity (Table 4). The similarities on the amino acid level followed the same pattern as those seen on the nucleotide level (Table 5).
Comparison of the regulatory genes in the HA and Orfx neurotoxin clusters.
Both the HA and Orfx clusters contained a regulatory gene designated botR or p21 based on substantial similarity, but the gene was most homologous to those of the same neurotoxin cluster type. In addition, the two genes were in the opposite orientation to each other within their specific neurotoxin clusters. For botR in the HA cluster, the A4 strain exhibited 98.5% nucleotide identity and 97.8% amino acid similarity and identity to A1 (Tables 2, 3). Compared to the NCTC 2916 HA cluster, A4 demonstrated 99.4% identity on the nucleotide level and was 99.4% identical and similar on the amino acid level.
For the Orfx cluster, p21 in A2 was most similar to p21 in A4 with 95.1% nucleotide identity, similar to HA– Orfx+ A1 with 94.4% identity, and finally similar to A3 with 93.5% identity (Table 4). The amino acid order similarities followed a parallel pattern except that HA– Orfx+ A1 had a slightly higher amino acid identity compared to A2 than A4, indicating that some of the nucleotide substitutions were silent mutations (Table 5).
Comparison of the ntnh genes.
ntnh is the only gene other than bont that is present in both the HA and Orfx clusters. On the nucleotide level, the ntnh genes showed closest homology within similar neurotoxin clusters, and ntnh in the A1 HA cluster had a high level of homology compared to the ntnh genes in the HA clusters of the A4 and NCTC 2916 strains. The ntnh gene in the Orfx cluster of A2 had high homology to the ntnh genes in the Orfx clusters of A3 and HA– Orfx+ A1, as well as the Orfx cluster associated with the BoNT/A in A4 and NCTC 2916 (Table 6). The same trend was observed on the amino acid level (Table 6 and Fig. 2). An insertion of 99 nucleotides in the ntnh genes of the HA cluster was detected, which gave a gene of 3,821 bp in length compared to the 3,722 bp ntnh in the Orfx cluster. As a result, the NTNH in the HA cluster contains an additional 33 amino acids starting at position 113 (Fig. 3).
|
View this table: [in a new window] |
TABLE 6. Comparison of NTNH on the nucleotide and amino acid levels among C. botulinum type A strainsa
|
|
View larger version (4K): [in a new window] |
FIG. 2. Graphical representation of degree of relatedness between the amino acid sequences of the NTNH proteins found in various BoNT/A producing bacteria (drawn to scale). The degree of such relatedness depends on different types of clusters. NTNHs from the HA clusters of different strains are closely related; so are those from Orfx clusters of different strains. However, the degree of relatedness among NTNHs between the two cluster types is lower than that within the same cluster type. NTNH-A1, NTNH in the HA cluster with BoNT/A1; NTNH-A2, NTNH in the Orfx cluster with BoNT/A2; NTNH-A3, NTNH in the Orfx cluster with BoNT/A3; NTNH-A4-A, NTNH in the Orfx cluster with BoNT/A4; NTNH-A4-B, NTNH in the HA cluster with BoNT/B; NTNH-ha-/orfx+ A1, NTNH in the Orfx cluster with BoNT/A1 in the HA– Orfx+ A1 strain; NTNH-NC-A, NTNH in the Orfx cluster with BoNT/A1 in NCTC 2916 strain; NTNH-NC-B, NTNH in the HA cluster with the silent BoNT/B in NCTC 2916.
|
|
View larger version (15K): [in a new window] |
FIG. 3. Alignment of a region of the amino acid sequences of the NTNH genes associated with BoNT in known BoNT/A-producing strains using AlignX software. This figure shows that NTNHs in the Orfx cluster have amino acid deletions compared with those in the HA cluster. The numbers above the alignment indicate the locations of the amino acids of NTNHs in the NTNH genes. All sequences start at the 96th amino acid. Red letters in yellow blocks indicate conserved amino acids, blue letters in light blue blocks indicated conserved amino acids, green letters in white blocks indicate that the amino acids are changed but retain similar properties compared to the amino acids in the other samples, black letters in white blocks indicate distinctness compared to the other samples, and dashes indicate deletions.
|
|
View this table: [in a new window] |
TABLE 7. Comparison of BoNT on the nucleotide and amino acid levels among C. botulinum type A strainsa
|
|
|
|---|
The notion that several unique recombination events led to the neurotoxin cluster organization observed in the strains analyzed is also supported by MLST experiments performed in our laboratory for strain comparison and delineation (M. Jacobson et al., unpublished data). The A4 and HA– Orfx+ A1 strains differed in only one of the seven gene alleles examined. However, the A2 and A3 strains were highly dissimilar as they did not have any allele sequences in common with any of the other C. botulinum strains or with each other. A(B) strains demonstrated a high level of heterogeneity since multiple MLST patterns were observed. Given the lack of evidence for a common clonal relationship among strains containing similar (HA or Orfx) neurotoxin clusters, it is unlikely that the neurotoxin clusters are the result of some shared ancestry among the isolates examined; although there is a high degree of similarity among neurotoxin clusters, the strains themselves demonstrated a high degree of genetic variability. A more plausible explanation for the nucleotide sequence similarity among the genes of similar (HA or Orfx) neurotoxin clusters is that recent recombination events have occurred in the history of such strains, leading to the creation of the BoNT neurotoxin clusters within the snapshot in time observed today.
The ntnh gene is the only gene other than bont present in all neurotoxin clusters from C. botulinum, including both HA and Orfx clusters. NTNHs from HA clusters of different strains are closely related, as are those from Orfx clusters of different strains. However, the ntnh genes were most similar to other ntnh genes of the same neurotoxin cluster. Specifically, the ntnh genes associated with HA clusters were more similar to each other than to the ntnh genes associated with the Orfx clusters. One possible explanation for this difference is that some or all of the ntnh gene was transferred along with the HA or Orfx cluster genes via a recombination event. This hypothesis is supported by the presence of an insertion of 99 nucleotides in the ntnh genes of all the HA clusters examined. This insertion translates to an additional 33 amino acids starting at position 113 in the NTNH of the HA cluster (Fig. 3).
Analysis of the orfx3 gene demonstrated that genes from the HA– Orfx+ A1 and NCTC 2916 strains were most similar to each other and that those from the A2, A4, and A3 were more divergent. A similar trend was observed for the orfx2 and orfx1 genes, with the HA– Orfx+ A1 and NCTC 2916 strains being highly similar while A2 and A3 were similar to each other, and A4 was the most divergent of all the strains tested. This trend was somewhat expected since both the NCTC 2916 and HA– Orfx+ A1 strains exhibited a BoNT/A1 gene associated with an Orfx cluster. However, the two strains are distinct. The NCTC 2916 strain exhibits an unexpressed BoNT/B gene with an HA cluster but the HA– Orfx+ A1 strain completely lacks this neurotoxin cluster and bont/b. Therefore, it is possible that the two strains gained their respective BoNT/A-associated neurotoxin clusters from a similar recombination event separately from the insertion of the BoNT/B neurotoxin cluster in NCTC 2916. This view is further supported when the p21 and ntnh genes of the Orfx clusters are compared as the HA– Orfx+ A1 and NCTC 2916 strains showed the highest degree of similarity for the strains tested for the two genes.
The p47 gene is particularly interesting as it is the only one of the Orfx cluster-specific genes that is in the same orientation and likely to be cotranscribed (8) with ntnh and bont. Comparison of the p47 nucleotide and amino acid sequences indicated that the HA– Orfx+ A1 and A2 strains were most similar, which is not consistent with the results observed with the other genes. Additionally, it was observed that p47 genes in the A4 and NCTC 2916 strains were highly divergent from other strains and each other. The specific role of the p47 gene is unknown, but it may either serve as a regulatory gene or encode a component of the toxin complex. Further studies will be required to elucidate its function and the significance of its genetic lineages within the subtypes.
A clearer relationship was observed when the HA clusters were compared as the A4 and NCTC 2916 BoNT/B-associated neurotoxin clusters demonstrated high levels of similarity for all of the genes tested. This further supports the theory that the neurotoxin clusters arose from genetic events distinct from the evolution of the clostridia since, as mentioned above, the NCTC 2916 and A4 bacteria did not demonstrate the same degree of relatedness when they were compared using an MLST procedure as they did when their neurotoxin cluster genes were compared. Based on these observations and the fact that the BoNTs of each strain are so different, it is highly likely that these neurotoxin clusters also arose from recombination events distinct from those that produced the BoNT genes.
Therefore, the bont gene and its associated neurotoxin cluster genes appear to have been acquired and evolved distinctly in the C. botulinum type A strains examined. The origin and mechanisms of acquisition of the toxin gene neurotoxin clusters, as well as their conservation of gene arrangement and sequence homology, are intriguing. It is likely that their acquisition took place through horizontal gene transfer by plasmid- or phage-mediated gene transfer. Notably, the A3 and A4 gene clusters reside on plasmids; however, the mobility of these plasmids among clostridial strains has not been examined. It appears likely that the acquisition of the neurotoxin clusters occurred as a series of evolutionary events more recent than the divergence of the neurotoxigenic clostridia from a common precursor organism since the neurotoxin cluster genes show relatively high homology. The neurotoxin cluster genes have been shown to be dispensable since the Johnson laboratory has isolated strains that have deletions of the entire A1 neurotoxin cluster and yet have the same physiology and growth properties as the toxigenic parents, to the best of our knowledge based on comparative studies. Lastly, the origin of the bont and neurotoxin cluster genes is a very intriguing mystery. Due to the polyprotein properties of the proteins in the complex and their proteolytic maturation, it has been postulated that the genes may have originated from a mammalian virus, perhaps of a neurotropic source (5, 18).
The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention.
Published ahead of print on 7 March 2008. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»