In this study, we sequenced bont/A, bont/B, and the associated neurotoxin cluster genes of C. botulinum type Ab strain CDC41370 and the bont/A gene of C. botulinum type Ab strain CDC588. These nucleotide sequences were compared with bont/A and bont/B sequences from other strains, including types Ab, A(B), and Ba. In addition, we analyzed the organization of the neurotoxin gene clusters in C. botulinum type Ab strains CDC588 and CDC1436.
The
bont/A nucleotide sequence of strain CDC588 was 99.9% identical to that of subtype A1. In contrast, the
bont/A gene from strain CDC1436 is identical to that of BoNT/A2 strains (
12). The
bont/A nucleotide sequence from strain CDC41370 was most similar to
bont/A1, although it differed from this subtype by 2.2%, and the predicted amino acid sequence differed by 4.3%. Among strains of subtype A1, little nucleotide sequence variation has been shown; i.e.,
bont/A nucleotide sequences in A(B) strains harbor only two nucleotide differences from those of other A1 strains (
12), and
bont/A nucleotide sequences in HA
− Orfx
+ A1 strains differ from those of other A1 strains by 4 to 5 nucleotides (
18). In comparison, the difference of 87 nucleotides observed in the
bont/A nucleotide sequence of strain CDC41370 is remarkable; however, none of the changes in the predicted amino acid sequence observed in strain CDC41370 localized to regions encoding the known conserved functional motifs (
1). These variations could cause modifications in the structure of the toxin, with consequent implications such as differences in the binding affinity for neuronal cell receptors or in BoNT antitoxin binding, thereby having an effect on the potency of a given antitoxin. Traditionally, botulinum toxin types (serotypes A through G) have been defined by specific neutralization with antitoxins derived from hyperimmunized animals, and some subtypes have been shown to have differential binding of monoclonal antibodies (
10,
16,
17,
20). For example, subtypes A1 and A2 present differences in monoclonal antibody binding affinity, despite the fact that these subtypes share 90% of amino acid sequence identity.
Toxin subtypes were historically identified through monoclonal antibody analyses, but this approach is labor-intensive. Currently, subtypes are defined based on the identification of separate clusters by using phylogenetic analysis of the
bont gene sequence (
12). For example,
bont/E1 and
bont/E2 are classified as different subtypes, although they differ by only 1% at the amino acid level. Carter et al. (
2) recently described a
bont/A gene nucleotide sequence in strains from heroin-associated wound botulism cases in the United Kingdom that differs from that from subtype A1 by 1.4%, and they proposed that this unique sequence be classified as subtype A5. According to these criteria, the
bont/A gene sequence of strain CDC41370 could represent a different BoNT/A subtype since the
bont gene nucleotide sequence differs by 2.2% from the reported A1 sequence. However, the impact of the variation observed in the
bont/A sequence of this strain on functional properties of the toxin is unknown. Moreover, the probability of describing additional unique sequences will likely increase as more diverse strains are examined, questioning the value of clustering neurotoxin gene sequences into subtypes.
Comparative analysis of
bont/A from strain CDC41370 with
bont/A1 to
bont/A4 suggests that this gene is derived from recombinant events involving
bont/A1,
bont/A2, and probably other
bont/A genes. Similarly, subtype A2 could represent a recombination event between
bont/A1 and
bont/A3 (
12). In addition, recombination in the
ntnh gene of certain type A and B strains has been reported (
3). These findings indicate that recombination events have likely contributed to the diversity among the BoNT serotypes. Further study of the mechanisms producing such diversity in the neurotoxin gene cluster may shed light on the role of such genes, if any, on toxin production and potency.
Previously reported
bont/B nucleotide sequences from bivalent strains harboring two
bont genes, i.e., A(B), Ba, Bf, and Ab, are highly similar and represent a separate
bont/B subtype (bivalent B) (
12). Interestingly, the
bont/B nucleotide sequence from strain CDC41370 showed higher identity to the sequence from subtype B1 than to other bivalent strains. The
bont/B in strain CDC41370 differed from the bivalent B subtype by up to 2.3%, suggesting a different origin of this gene in strain CDC41370 than that in other bivalent strains. In addition, two
C. botulinum type Ab strains from Italy exhibited B1 and B3 subtypes using a PCR-restriction fragment length polymorphism method (
9). These findings support the hypothesis of mobilization of the
bont/B gene among bivalent strains.
Among A(B) strains, the nucleotide sequence analysis of
bont/A and
bont/B genes showed a relatively low variability. The A(B) strains sequenced during this study harbored identical
bont/A nucleotide sequences to that of previously published sequences in other A(B) strains (
12). All of the
bont/B nucleotide sequences from A(B) strains examined in this study and others (
13) harbored a substitution of a T at nucleotide 384, resulting in a premature stop codon. Although accumulation of nucleotide changes after nonsense mutations may be expected, the nucleotide changes in the
bont/B sequences were conserved among the A(B) strains examined. Moreover, most of the nucleotide changes observed in these silent
bont/B genes were also identified in the previously reported
bont/B gene of type Ab and Ba strains. These data suggest that A(B) strains have highly conserved
bont/A and
bont/B gene sequences. In contrast, strains producing two toxin types (e.g., types Ab and Ba) showed comparatively more variability because at least four different
bont/A sequences and three different
bont/B sequences have been described in these bivalent strains.
In addition to
bont/A and
bont/B genes, we sequenced the neurotoxin cluster genes of the
C. botulinum strain CDC41370 and compared these sequences to known sequences from types A and B. To our knowledge, this is the first report of the complete nucleotide sequences of both neurotoxin gene clusters in a type Ab strain. In addition, we analyzed the content of the neurotoxin gene cluster in two other type Ab strains (strains CDC588 and CDC1436) and seven A(B) strains (CDC28184, CDC1727, CDC5178, CDC4893, and CDC5277) by PCR assays. As expected (
4,
15, 21), all 10 strains contained
bont/A,
bont/B,
ntnh,
ha70,
ha17,
ha33,
botR,
p21,
p47,
orfX1,
orfX2, and
orfX3. However, a difference in intergenic spacing of
p21-orfX1 and
orfX2-orfX3 was observed among these strains. Strain CDC588 and all six A(B) strains contained a 0.6-kb insertion between
orfX2 and
orfX3. This has been also described in other A(B) and HA
− OrfX
+ A1 strains (
5,
14,
18). These intergenic spacing sequences contain partial insertion sequence (IS) elements, suggesting a role in gene mobility (
4). Various partial IS elements have been described within the
orfX cluster in A2, A3, Ba, and A(B) strains (
5,
19). However, two type Ab strains analyzed in this study (strains CDC1436 and CDC41370) do not have the IS elements residing within these intergenic sequences. The lack of these potential IS elements has been also reported for the subtype A2 strain Mascarpone (
8). The significance of these partial IS elements in the evolution of the neurotoxin gene cluster remains unclear.
Analysis of the nucleotide sequence of the
orfX cluster genes associated with
bont/A in strain CDC41370 showed variable results compared to the
orfX cluster genes of various subtypes. The nucleotide sequence of the
orfX3 gene of strain CDC41370 was more similar to that of subtype A4, whereas the
orfX2 and
orfX1 nucleotide sequences had the highest similarity with those of A(B) strains. The sequence of the
p21 gene was the most divergent among the sequences within the
orfX cluster. The
p47 gene of strain CDC41370 was more similar to the sequence of subtype A3, and the
ntnh gene sequence showed high identity with that of subtype A2. Genetic variability in the
orfX cluster sequences has been also demonstrated among all the
bont/A subtypes (
14). Differences in toxin gene cluster arrangement among various subtypes have been described only for type A subtypes. In contrast, the neurotoxin cluster genes associated with the
bont/B gene seemed to be less diverse than the
orfX cluster. The nucleotide and amino acid sequences of the neurotoxin cluster genes associated with the
bont/B gene in strain CDC41370 were highly similar to those of the
ha cluster in A(B), Ba, and B1 strains.
BoNT-producing
Clostridium strains are genetically and phenotypically diverse (
3).
C. botulinum strains can be grouped into four phylogenetically distinct lineages, and rare strains of
Clostridium butyricum and
Clostridium baratii are able to produce BoNT/E and BoNT/F, respectively. Based on these observations, a revision of the nomenclature and classification within this species has been proposed (
3). More recently, various molecular analyses have demonstrated a high degree of variability among the BoNT-producing strains, the
bont genes, and the neurotoxin-associated protein-encoding genes (
12,
14). The description of novel arrangements of the neurotoxin cluster genes and variations in the
bont gene sequence, as described in this study and elsewhere, suggests that the nomenclature for BoNT subtypes requires revision. The value of clustering neurotoxin gene sequences into subtypes may diminish as additional variants of the
bont sequences are identified and differences among clusters become less evident. In addition, the description of novel toxin gene cluster arrangements increases the complexity of subtype nomenclature, e.g., HA
− Orfx
+ A1 strains. Thus, the use of only neurotoxin gene sequences as the basis for subtype nomenclature may not be accurate. Additional molecular characterization of BoNT-producing strains and analysis of the functional impact of differences in toxin gene sequences may help to improve the classification of such diverse strains and their associated toxins.