Clostridium botulinum group III: a group with dual identity shaped by plasmids, phages and mobile elements
© Skarin et al; licensee BioMed Central Ltd. 2011
Received: 19 January 2011
Accepted: 12 April 2011
Published: 12 April 2011
Clostridium botulinum strains can be divided into four physiological groups that are sufficiently diverged to be considered as separate species. Here we present the first complete genome of a C. botulinum strain from physiological group III, causing animal botulism. We also compare the sequence to three new draft genomes from the same physiological group.
The 2.77 Mb chromosome was highly conserved between the isolates and also closely related to that of C. novyi. However, the sequence was very different from the human C. botulinum group genomes. Replication-directed translocations were rare and conservation of synteny was high. The largest difference between C. botulinum group III isolates occurred within their surprisingly large plasmidomes and in the pattern of mobile elements insertions. Five plasmids, constituting 13.5% of the total genetic material, were present in the completed genome. Interestingly, the set of plasmids differed compared to other isolates. The largest plasmid, the botulinum-neurotoxin carrying prophage, was conserved at a level similar to that of the chromosome while the medium-sized plasmids seemed to be undergoing faster genetic drift. These plasmids also contained more mobile elements than other replicons. Several toxins and resistance genes were identified, many of which were located on the plasmids.
The completion of the genome of C. botulinum group III has revealed it to be a genome with dual identity. It belongs to the pathogenic species C. botulinum, but as a genotypic species it should also include C. novyi and C. haemolyticum. The genotypic species share a conserved chromosomal core that can be transformed into various pathogenic variants by modulation of the highly plastic plasmidome.
Botulism is a paralytic disease caused by some of the most potent toxins known, the botulinum neurotoxins (BoNTs). The toxins are proteins mainly produced by the species Clostridium botulinum but some Clostridium butyricum and Clostridium baratii strains are also capable of producing BoNTs. Due to their extreme potency, the BoNTs are classified as high risk threat agents for bioterrorism . In addition to causing severe intoxications in humans, BoNT-producing bacteria can also produce serious problems in wild and domesticated animals such as birds, cattle, horses, sheep and minks. Outbreaks with high mortal outcome in waterfowl and poultry have become an increasing environmental and economical problem .
C. botulinum is an anaerobic, spore-forming Gram-positive bacterium commonly found in soil and aquatic environments . The species C. botulinum is divided into four physiological groups (I-IV), which produce BoNTs of seven different serotypes (A-G) . The four groups represent distinct phylogenetic clades and are separated by a distance large enough to merit classification as four separate species . Physiological group I (proteolytic) includes C. botulinum strains producing toxins of serotype A, B or F and is mainly associated with human cases. Physiological group II (non-proteolytic) consists of strains producing toxins of serotype B, E or F; these strains also cause human botulism. Group III produces toxins of serotype C or D and is associated with avian and nonhuman mammalian botulism. Genetic studies have shown that strains belonging to physiological group III are more closely related to C. novyi and C. haemolyticum than to C. botulinum serotypes from group I and II . Physiological group IV is a rarer and less characterised group that produces toxin of serotype G.
The neurotoxin genes of C. botulinum type C and D are carried by bacteriophages, which express unstable lysogeny and are frequently lost during cultivation. The BoNT prophage propagates in the bacterium as a large plasmid and strains cured from the phage can be reconverted to toxigenic strains by either type C or D phage particles . The distinction between types C and D is not absolute because chimerical sequences exist . In general, the neurotoxin genes analysed from avian isolates comprise parts from both BoNT/C and BoNT/D genes and are referred to as type C/D or D/C . The chimeric types are more lethal to avian species than either type C or D [9, 10]. C. botulinum types C and D also produce a minor amount of a binary toxin, the C2 toxin, which genes are located on a plasmid . The C2 toxin has a translocation domain as well as an ADP-ribosylating domain that targets actin.
Several completed clostridia genomes have been reported in the last ten years, including C. difficile, C. perfringens, C. novyi-NT and C. botulinum types A, B, E and F [12–16]. The genome sizes of species in the Clostridium genus vary between 2.5 and up to at least 6 Mb and their GC contents are often low. Sequences from three large plasmids from C. botulinum group III strains have been released: the BoNT prophage from strain C-Stockholm ; the C2 toxin plasmid from strains C203U28  and D-1873; and another large plasmid from strain D-1873. Two draft whole genome shotgun (WGS) assemblies from type D strain 1873 and type C strain Eklund are also deposited in the sequence database.
Here we report and analyse the first complete genomic sequence from a C. botulinum group III strain. The sequence comes from a recently isolated C. botulinum type C/D strain originating in 2008 from an outbreak in a Swedish poultry farm. We also release three draft genomes: one from a wildfowl outbreak in 2007, another from a poultry outbreak in 2007, and finally the C-Stockholm strain isolated from a mink outbreak in 1949 . The results have revealed a genome with a conserved core that is shaped by plasmids, phages and mobile elements.
Results and Discussion
General genome features
General genomic features of the strains analysed in this study, the C. novyi-NT genome and the so far released genomic sequences of C. botulinum group III plasmids.
Species and type (strain)
Origin (year of isolation)
GenBank accession number
C. botulinum C/D
2 773 191
3 207 676
C. botulinum C/D
C. botulinum C
C. novyi A
2 547 720
(C. novyi-NT) b
C. botulinum D
C. botulinum C
C. botulinum C
Upon comparing the available genomes of group III strains, we found three different clades (Figure 2B). Our recently isolated strains formed one lineage, the D-1873 and the C-Stockholm strains formed a second one and the C -Eklund (type C/D ) together with the C. novyi-NT strain formed a third lineage. The average genomic variation within group III (C. novyi sensu lato) was large, but not exceptionally so.
Features of the chromosome
The distribution of putative genes coding for virulence proteins, antibiotic resistance-conferring proteins, and CRISPR-proteins found in the C. botulinum group III completed replicons of the strains analysed in this study.
Putative virulence protein
Chemotaxis related protein
Flagella associated protein
Putative antibiotic resistance
Multidrug resistance protein
Microcin self-immunity protein
Tellurium-resistance protein terD
Toxic-anion resistance protein
Features of the Plasmidome
The five plasmids together comprised 13.5% of the total genomic content of strain BKT015925. To date, this is the largest number of plasmids in any completed clostridia genome. In fact, in all the currently completed microbial genomes, less than 5% of them contain five or more plasmids. The copy number, as determined by coverage, was close to one for the BKT015925 plasmids except for the smallest one, where it was approximately two. By mapping the draft genome reads onto the sequences of plasmids from BKT015925, we concluded that p1-p4 also were present and highly similar (99%, 98%, 97% and 99% nucleotide similarity respectively) in V891. However, the p5 plasmid was absent. By PFGE and analysis of de novo contigs unique to V891 and C-Stockholm we could also detect the presence of a sixth plasmid of about 55 kb. This plasmid was completed using the C-Stockholm genome sequence and was therefore called p6CSt. We believe that the plasmids we present here are part of a pan-plasmidome existing in C. botulinum group III. Conjugative transfer of BoNT encoding plasmids has previously been reported between strains of C. botulinum, and it is likely this mechanism that is responsible for exchange of plasmids between strains .
Many of the plasmids contained phage genes and several toxin genes were present in the plasmidome (Table 2). Two of the plasmids also contained clustered, regularly interspaced, short palindromic repeat (CRISPR) regions and genes encoding CRISPR-associated (Cas) proteins, providing acquired immunity against foreign DNA . Typically, CRISPR regions are located on the chromosome and the C. novyi-NT genome is no exception with four chromosomal CRISPR regions. Interestingly, CRISPR loci and cas genes were found exclusively on the plasmidome in BKT015925. The number and identity of the spacer sequences on p2 were identical to matches found in the other poultry strain, BKT028387, confirming that the two poultry strains are very closely related.
The progenitor neurotoxin cluster in Group III strains consists of a regulatory protein (BotR), hemagglutinin (HA) proteins and a non-toxic non-hemagglutinin (NTNH) protein. These proteins are suggested to be involved in protection and transportation of the neurotoxin . The complete gene cluster was found in p1, together with the mosaic neurotoxin gene coding for the C-type light (L)-chain and the C/D chimerical heavy (H)-chain. The L-chain (residues 1-444) showed a 98% sequence identity with the corresponding residues in C-Stockholm, but only 48% identity for the same region in type D-1873. By contrast, the C-terminal part of the H-chain (residues 863-1280) shared 42% sequence identity with phage C-Stockholm, but 95% to D-1873. The N terminal part of the H-chain gene showed high similarity to both C and D types and was probably the place where the recombination event once occurred.
The second-largest plasmid was nearly 100 kb. An interesting feature found in it was an alpha-toxin with conserved domains of the aerolysin superfamily. It showed highest similarity to the alpha-toxin of C. septicum, which has both necrotic and haemolytic properties . Another interesting feature was a homologue to the gene coding for clostripain, which was found downstream of the alpha-toxin. This type of cysteine endopeptidase has been found in many clostridia, for example in C. novyi and C. botulinum group I genomes, but in those cases it was located on the chromosome. The α-clostripain from C. perfringens has been suggested to be involved in pathogenesis . Two genes coding for metallo-beta lactamases and one for penicillin amidase were also found in this plasmid; these may contribute to resistance to beta-lactam antibiotics. Three CRISPR regions were found, consisting of 28, 5 and 7 spacers. We also identified three Cas proteins (Cas3, 5 and 6) in proximity to the CRISPR loci.
We compared the p2 plasmid to other group III plasmids (Figure 4B). Some conserved regions were present, especially in the pCLG2 plasmids from C. botulinum type D strain 1873. Interestingly, two large, conserved blocks contained all putative virulence factors (except for one protease), two of the Cas genes, and the penicillin amidase gene. Small, scattered, and conserved regions were also shared with pC2C203U28 (and pCLG1) and p3, indicating that genetic crossover events have occurred. This is in contrast to the chromosome and BoNT phage where rearrangements have been few.
Among the more interesting features of this 80 kb plasmid were two genes coding for putative lantibiotic ABC transporters, and the two genes coding for C2 toxin. Interestingly, the two C2 toxin genes were also found on plasmid pC2C203U28 (and pCLG1), which otherwise shared little resemblance with the p3 plasmid (Figure 4B). The amino-acid sequence of the C2I component shared 96% sequence identity with that of plasmids C-203U28 and pCLG1 and 100% with the C2I sequence of C. botulinum strain (C) 2300 . In contrast, the C2II component was more diverged compared to strains 2300 and C-203U28 (88% identity to both). The C-terminal extension found in strain 2300 was also missing. Consistent with previous observations , the divergence in the C-terminus of the C2II sequence contributed to most of the differences in sequence identity.
The majority of the coding sequences in this 40 kb plasmid were phage related, suggesting that this is a prophage. However, it also contained genes coding for plasmid segregation and replication. Interestingly, the GC content and GC skew were different for the prophage and the plasmid-replication part. This could indicate that the prophage may have hijacked a plasmid-replication system.
The smallest of the plasmids (12 kb), found only in the poultry isolates, contained two putative Epsilon type B toxins. The sequence identity between the two toxin genes was quite low (~40%). A putative homologue was also found in the pCLG1 plasmid but the conservation was low there as well. Two CDS were phage related and all genes but one were located on the same strand. Interestingly, this was a repressor-coding gene that might regulate the transcription of the opposite, gene-dense, strand.
Among the CDS identified in this 55 kb plasmid found in C-Stockholm (and V891), eleven were phage related, suggesting that it is a prophage as well. Other genes of interest were two for penicillin amidase and one for a bacteriocin. The latter possibly inhibits growth of competing bacterial strains. All CDS, except for four, were located on the same strand.
Mobile elements identified on the chromosome of strains BKT015925, V891 and C. novyi-NT.
Included among the high-copy mobile elements that showed an elevated transposon activity, there were elements related to the IS6, IS21, IS256 and IS1182 families. There was also a trans-located cassette of two, tandem reverse-transcriptase (RT) genes. In total, the mobile elements accounted for 3.4% of the genetic content of the BKT015925 chromosome and tended to be placed in non-coding regions of both replichores (Figure 1). Except for the IS21-related element, the mobile elements were inserted in intergenic regions, thereby avoiding gene disruptions. IS elements were also found on the three largest plasmids. Two types of elements belonging to the IS200/IS605 family were found exclusively in the plasmidome. In fact, the density of mobile elements was highest on p2 and p3, where 8 and 12% respectively of the total genetic content constituted of such elements. This was interesting because these two plasmids seem to be more genetically dynamic than the chromosome and p1 plasmid. The exchange of plasmids between strains further increases the possibility of horizontal gene transfer facilitated by IS elements.
The IS6 element, ISCbo1, was most related to an IS element, ISCpe7, found in C. perfringens. The ISCbo1 elements were found in 25 places (two broken) in the BKT015925 chromosome, of which only eight were on the same place in the V891 draft genome (Figure 5). ISCbo1 elements were also found on p2 and p3. In most cases, the elements were located directly upstream a gene and were surrounded by 8-9-bp non-conserved, direct repeats flanking a 10-bp inverted repeat (Additional file 4, Table S3). Although these elements were not found on the BoNT-carrying prophage in BKT015925, they had previously been identified in a truncated form, between orfX and botR in the bont cluster of C. botulinum group I strains. It has been suggested that these elements, together with other closely located transposons, have been the responsible machinery causing horizontal transfer of the bont cluster between strains of C. botulinum.
The IS21-related element, ISCbo2, consisted of two adjacent genes, integrase and istB. This element was also found in the C. beijerickii NCIMB 8052 genome. IstB has an ATP-binding motif, and is thought to be involved in facilitating the contact between target DNA and a transposase, in this case an integrase . The ISCbo2 element was found on six positions on the chromosome of BKT015925, and was also present on the three largest plasmids. It showed a tendency to be incorporated into another mobile element, a tandem reverse-transcriptase (ISCbo3), thereby disrupting its structure. The incorporations were not at the same positions within the RT gene. Small remnants of the IsCbo3 element were found at several places in the genome which could be a consequence of ISCbo2 insertions.
The IS256 family element, ISCbo4, found in the C. botulinum group III genome showed high amino-acid sequence similarity (71% identity) to ISCpe3 found in C. perfringens. The ISCbo4 element was identified on five locations on the BKT015925 chromosome and on p1-p3. The insertion pattern in V891 was quite different from that in BKT015925 (Figure 5). The element was also found on one location in the C. novyi-NT genome.
The IS1182 element ISCbo5 was previously identified (as ISbma2) in numerous copies in the genome of C. novyi-NT and Burkholderia pseudomallei/B. mallei strains (4 and 48 copies, respectively) [14, 33]. Seventeen copies (including one broken) of this element were located on the BKT015925 chromosome, while only twelve copies were found in V891. Six of the copies in V891 were found on the corresponding loci in BKT015925. However, in some cases the insertion sites were slightly altered, which indicates that they have once left and then been reinserted in the same position. There was only one copy within the plasmidome (on p3). All insertion sites except one were composed of inverted repeats forming a hairpin structure. The hairpin structures were exclusively located at the end of genes and most likely act as transcription stop signals. The ISCbo5 element without a surrounding hairpin seems to indicate that it is transcribed together with the upstream gene. It has previously been suggested that such an element is the ancestral one . This specific ISCbo5 element was also broken, indicating it might not be tolerated by the host. The equivalent element (without a surrounding hairpin) in the C. novyi-NT genome was located next to a resolvase gene, but in strain BKT015925 it was next to a cysteine desulfurase gene. To test the hypothesis that it is a fused transcript, we amplified the mRNA of the cysteine desulfurase gene by RT-PCR. The gene was transcribed but was not fused with ISCbo5 (data not shown).
Seven different elements belonging to the IS200/IS605 family (ISCbo6-12) were found in p1, and one in p2 (ISCbo13), where five of the IS605 elements were located autonomously and two were found adjacent to IS200 elements. ISCbo7 existed in two copies, both autonomously and adjacent to ISCbo6. Interestingly, this class of elements was only found in the plasmidome. One of the IS200/IS605 elements was highly conserved in the phage C-Stockholm (there called ISCbt4) but it was not located at the same position. By comparing flanking sequences, we located all IS200/IS605 elements to the same positions in the V891 genome, suggesting that these elements are less active than the previously discussed ones.
The completion of the genome of C. botulinum group III has revealed a relatively small chromosome indicating a low tolerance for redundant or excessive genetic material. Acceptance for genes on the lagging strand was also low indicating that high transcription efficiency is crucial for this organism. Although there was a low tendency for rearrangements in the chromosome, size variation between strains indicated that there is an on-going insertion and removal of genetic material. It is probable that the abundant mobile elements are contributing to the genome remodelling process. The C. botulinum group III is genetically remote from the types that can cause human botulism, but closely related to C. novyi. There are larger differences between certain pairs of C. botulinum group III strains than between some group III strains and C. novyi. On the basis of these facts, and in analogy with the handling of the Bacillus cereus group, we propose a new genotypic name (genospecies), C. novyi sensu lato, for C. botulinum group III, C. novyi and C. haemolyticum. However, it is unlikely that the pathotypic name (pathospecies) is going to be changed from C. botulinum.
In contrast to the conserved chromosome, we discovered that physiological group III has a remarkably large and variable plasmidome where closely related isolates can have different sets of plasmids. The largest plasmid, the BoNT phage, seemed to obey the same evolutionary rules as the chromosome. In the medium-sized plasmids, and to some extent in the smaller ones, genetic material had been rearranged at a much higher rate and with fewer constrictions. We hypothesise that these plasmids are part of a larger pan-plasmidome and that strains can exchange plasmids.
Mobile elements were abundant also within the plasmidome, where the medium-sized plasmids had the highest transposon density, and they may have been contributing to the mechanism of rearrangements. Interestingly, five out of six toxin types found within the C. botulinum group III genome were located in the plasmidome. Our results suggest a model for genome evolution where a conserved core is located in the chromosome while factors affecting pathogenicity are located on a more plastic plasmidome to allow faster environmental adaption without disturbing core functions.
Bacterial strains and preparation of genomic DNA
Three toxigenic C. botulinum type C strains (07-V891, 08-BKT015925 and reference strain C-Stockholm) and one non-toxigenic strain (07-BKT028387) were used in this study. Strain V891 was isolated by the authors from a wildfowl herring gull that had died from botulism in Sweden in 2007. Strains BKT015925 and BKT028387 were isolated by the authors from poultry broilers that had died during two separate botulism outbreaks in Sweden in 2008 . Bacterial strains were cultured overnight in anaerobic jars (Merck, Darmstadt, Germany) at 37°C in 9 ml pre-reduced TPGY broth. Cells were harvested by centrifugation at 3000 g for 15 min before DNA extraction. Genomic DNA was purified using the DNeasy blood-and-tissue kit (Qiagen, Hilden, Germany). DNA concentration was measured with a NanoDrop 2000 (NanoDrop Technologies, Wilmington, USA) and the DNA was analysed by gel electrophoresis.
Genome sequencing and assembly
Sequencing of the genomic DNA was done by the whole-genome-shotgun sequencing approach on the Roche Genome Sequencer Titanium system (Roche Applied Science, Mannheim, Germany). The C-Stockholm and BKT028387 strains were sequenced with the FLX chemistry. Approximately 5 μg genomic DNA was subjected to standard 454 shotgun sequencing. The BKT015925 strain was sequenced on a full picotiterplate resulting in 853 000 sequences with an average length of 288 nucleotides. The coverage was approximately 80×. The draft genomes (V891, BKT028387 and C-Stockholm strain) were produced on half a picotiterplate and yielded 385 000, 260 000, and 104 000 sequences respectively.
The sequences were assembled de novo using the GS assembler (Newbler, Roche Applied Science). Gaps were closed with help of the Consed package . Sanger sequencing of PCR products and local reassembly were used to resolve gaps, misassemblies and unclear regions. In total, 241 Sanger reads were incorporated into the assembly. Reference guided assemblies were produced with GS mapper or Consed-crossmatch. Genome alignments were made, primarily with Mummer  and with the Artemis comparison tool (ACT) .
Agarose plugs of strains BKT015925, V891, and C-Stockholm were produced  and the gel plugs were incubated with 0.1 units of S1 nuclease (Sigma, St.Louis, MO, USA) for 45 min . The digested DNA was analysed on a 1% agarose gel (Agarose NA; GE Healthcare, Little Chalfont, UK) and electrophoresis was performed in 14°C in HEPES buffer (16 mM HEPES-NaOH, 16 mM sodium acetate, 0.8 mM EDTA, pH 7.5). The settings were 4 V/cm for 26 hours at pulse switch time ramped from 0.5-15 s in a CHEF DRII apparatus (BioRad, Hercules, CA, USA). DNA from Salmonella Branderup strain H9812 digested with Xba I and a lambda PFG marker (New England BioLabs, Ipswich, MA, USA) were used as size markers.
Strain BKT015925 was cultured as described and RNA was extracted using Trizol according to the manufacturer's instructions (Invitrogen, Carlsbad, CA, USA). RNA was converted into cDNA using random oligonucleotides and the SuperScript III First-Strand Synthesis System (Invitrogen) according to the manufacturer's procedure. PCR was carried out with gene-specific primers using AmpliTaq Gold polymerase (Applied Biosystems, Foster City, CA, USA), with denaturation at 95°C for 30 s, annealing at 50°C for 30 s and extension at 72°C for 1 min.
Genes were defined with Glimmer 3 . The annotation process was handled using Artemis software . Circular plots were drawn with DNAplotter . Further comparisons were made to Swissprot and NCBI non-redundant protein database. Well-conserved genes were automatically annotated but less certain annotations were manually assigned. The length, identity, and coverage of the sequences of the subject versus query were inspected. tRNAs were identified with tRNAscan-SE  and rRNA by homology comparison to the C. novyi-NT genome. CRISPR loci were found and analysed using CRISPRFinder .
Mapping positions of mobile elements in draft genomes
A mobile-element-depleted genome sequence was generated by removing all high-copy mobile elements from the chromosome of BKT015925. This depleted genome sequence was then used as a consensus for all three genomes (BKT015925, BKT028387 and V891). The positions of the same multi-copy mobile elements in the draft genomes were identified by making a set of GS reference-guided assemblies for the elements and mapping the reads spanning the elements borders back to the mobile-element-depleted genome sequence.
Phylogenomic distance was measured using the Average Similarity of the conserved Core method (ASC) . When the distance was large and the core was poorly defined on the nucleotide level, total average similarity was used. The dendrogram was created by converting the similarity matrix to a distance matrix and calculating a tree using the neighbor-joining method  and PHYLIP 3.67 through the Mobyle platform http://mobyle.pasteur.fr/. The phylogenetic tree was plotted using PhyloDraw 0.82 http://www.bioinformatics.org/wiki/PhyloDraw.
Nucleotide sequence accession numbers
The sequences and annotations of the 08-BKT015925 chromosome and plasmids have been submitted to GenBank [GenBank:CP002410-CP002415]. The draft genome sequences of C-Stockholm, 07-BKT028387 and V891 were also deposited in GenBank [GenBank:AESA00000000, AESB00000000 and AESC00000000].
We thank Kicki Holmberg at Science for Life Laboratory for help with parallel sequencing and Elin Klarbring for help in gap closure. The work was funded by the Swedish Farmers' Foundation for Agricultural Research (SLF). The work was also supported by the European Union-funded Integrated Project BIOTRACER (contract 036272) and internal funding at National Veterinary Institute (SVA). This research was partly supported by/executed in the framework of the EU-project AniBioThreat (Grant Agreement: Home/2009/ISEC/AG/191) with the financial support from the Prevention of and Fight against Crime Programme of the European Union, European Commission - Directorate General Home Affairs. This publication reflects the views only of the author, and the European Commission cannot be held responsible for any use which may be made of the information contained therein. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
- Arnon SS, Schechter R, Inglesby TV, Henderson DA, Bartlett JG, Ascher MS, Eitzen E, Fine AD, Hauer J, Layton M, et al: Botulinum toxin as a biological weapon: medical and public health management. JAMA. 2001, 285: 1059-1070. 10.1001/jama.285.8.1059.PubMedView ArticleGoogle Scholar
- Skarin H, Lindberg A, Blomqvist G, Aspan A, Baverud V: Molecular characterisation and comparison of Clostridium botulinum type C avian strains. Avian Pathol. 2010, 39: 511-518. 10.1080/03079457.2010.526923.PubMedView ArticleGoogle Scholar
- Hatheway CL: Toxigenic clostridia. Clin Microbiol Rev. 1990, 3: 66-98.PubMed CentralPubMedGoogle Scholar
- Collins MD, East AK: Phylogeny and taxonomy of the food-borne pathogen Clostridium botulinum and its neurotoxins. J Appl Microbiol. 1998, 84: 5-17. 10.1046/j.1365-2672.1997.00313.x.PubMedView ArticleGoogle Scholar
- Hill KK, Smith TJ, Helma CH, Ticknor LO, Foley BT, Svensson RT, Brown JL, Johnson EA, Smith LA, Okinaka RT, et al: Genetic diversity among Botulinum Neurotoxin-producing clostridial strains. J Bacteriol. 2007, 189: 818-832. 10.1128/JB.01180-06.PubMed CentralPubMedView ArticleGoogle Scholar
- Sasaki Y, Takikawa N, Kojima A, Norimatsu M, Suzuki S, Tamura Y: Phylogenetic positions of Clostridium novyi and Clostridium haemolyticum based on 16S rDNA sequences. Int J Syst Evol Microbiol. 2001, 51: 901-904.PubMedView ArticleGoogle Scholar
- Oguma K: The stability of toxigenicity in Clostridium botulinum types C and D. J Gen Microbiol. 1976, 92: 67-75.PubMedView ArticleGoogle Scholar
- Moriishi K, Koura M, Fujii N, Fujinaga Y, Inoue K, Syuto B, Oguma K: Molecular cloning of the gene encoding the mosaic neurotoxin, composed of parts of botulinum neurotoxin types C1 and D, and PCR detection of this gene from Clostridium botulinum type C organisms. Appl Environ Microbiol. 1996, 62: 662-667.PubMed CentralPubMedGoogle Scholar
- Takeda M, Tsukamoto K, Kohda T, Matsui M, Mukamoto M, Kozaki S: Characterization of the neurotoxin produced by isolates associated with avian botulism. Avian Dis. 2005, 49: 376-381. 10.1637/7347-022305R1.1.PubMedView ArticleGoogle Scholar
- Miyazaki S, Sakaguchi G: Experimental botulism in chickens: the cecum as the site of production and absorption of botulinum toxin. Jpn J Med Sci Biol. 1978, 31: 1-15.PubMedView ArticleGoogle Scholar
- Sakaguchi Y, Hayashi T, Yamamoto Y, Nakayama K, Zhang K, Ma S, Arimitsu H, Oguma K: Molecular analysis of an extrachromosomal element containing the C2 toxin gene discovered in Clostridium botulinum type C. J Bacteriol. 2009, 191: 3282-3291. 10.1128/JB.01797-08.PubMed CentralPubMedView ArticleGoogle Scholar
- Sebaihia M, Wren BW, Mullany P, Fairweather NF, Minton N, Stabler R, Thomson NR, Roberts AP, Cerdeno-Tarraga AM, Wang H, et al: The multidrug-resistant human pathogen Clostridium difficile has a highly mobile, mosaic genome. Nat Genet. 2006, 38: 779-786. 10.1038/ng1830.PubMedView ArticleGoogle Scholar
- Shimizu T, Ohtani K, Hirakawa H, Ohshima K, Yamashita A, Shiba T, Ogasawara N, Hattori M, Kuhara S, Hayashi H: Complete genome sequence of Clostridium perfringens, an anaerobic flesh-eater. Proc Natl Acad Sci USA. 2002, 99: 996-1001. 10.1073/pnas.022493799.PubMed CentralPubMedView ArticleGoogle Scholar
- Bettegowda C, Huang X, Lin J, Cheong I, Kohli M, Szabo SA, Zhang X, Diaz LA, Velculescu VE, Parmigiani G, et al: The genome and transcriptomes of the anti-tumor agent Clostridium novyi-NT. Nat Biotechnol. 2006, 24: 1573-1580. 10.1038/nbt1256.PubMedView ArticleGoogle Scholar
- Sebaihia M, Peck MW, Minton NP, Thomson NR, Holden MT, Mitchell WJ, Carter AT, Bentley SD, Mason DR, Crossman L, et al: Genome sequence of a proteolytic (Group I) Clostridium botulinum strain Hall A and comparative analysis of the clostridial genomes. Genome Res. 2007, 17: 1082-1092. 10.1101/gr.6282807.PubMed CentralPubMedView ArticleGoogle Scholar
- Smith TJ, Hill KK, Foley BT, Detter JC, Munk AC, Bruce DC, Doggett NA, Smith LA, Marks JD, Xie G, Brettin TS: Analysis of the neurotoxin complex genes in Clostridium botulinum A1-A4 and B1 strains: BoNT/A3,/Ba4 and/B1 clusters are located within plasmids. PLoS One. 2007, 2: e1271-10.1371/journal.pone.0001271.PubMed CentralPubMedView ArticleGoogle Scholar
- Sakaguchi Y, Hayashi T, Kurokawa K, Nakayama K, Oshima K, Fujinaga Y, Ohnishi M, Ohtsubo E, Hattori M, Oguma K: The genome sequence of Clostridium botulinum type C neurotoxin-converting phage and the molecular mechanisms of unstable lysogeny. Proc Natl Acad Sci USA. 2005, 102: 17472-17477. 10.1073/pnas.0505503102.PubMed CentralPubMedView ArticleGoogle Scholar
- Dinter VZ, Kull KE: Ergänzende Untersuchungen über den Botulismus beim Nerz. Nord Vet Med. 1951, 3: 297-311.Google Scholar
- Carter AT, Paul CJ, Mason DR, Twine SM, Alston MJ, Logan SM, Austin JW, Peck MW: Independent evolution of neurotoxin and flagellar genetic loci in proteolytic Clostridium botulinum. BMC Genomics. 2009, 10: 115-10.1186/1471-2164-10-115.PubMed CentralPubMedView ArticleGoogle Scholar
- Segerman B, De Medici D, Ehling Schulz M, Fach P, Fenicia L, Fricker M, Wielinga P, Van Rotterdam B, Knutsson R: Bioinformatic tools for using whole genome sequencing as a rapid high resolution diagnostic typing tool when tracing bioterror organisms in the food and feed chain. Int J Food Microbiol.
- Eisen JA, Heidelberg JF, White O, Salzberg SL: Evidence for symmetric chromosomal inversions around the replication origin in bacteria. Genome Biol. 2000, 1: RESEARCH0011-10.1186/gb-2000-1-6-research0011.PubMed CentralPubMedView ArticleGoogle Scholar
- Tillier ER, Collins RA: Genome rearrangement by replication-directed translocation. Nat Genet. 2000, 26: 195-197. 10.1038/79918.PubMedView ArticleGoogle Scholar
- Marshall KM, Bradshaw M, Johnson EA: Conjugative botulinum neurotoxin-encoding plasmids in Clostridium botulinum. PLoS One. 2010, 5: e11087-10.1371/journal.pone.0011087.PubMed CentralPubMedView ArticleGoogle Scholar
- Horvath P, Barrangou R: CRISPR/Cas, the immune system of bacteria and archaea. Science. 2010, 327: 167-170. 10.1126/science.1179555.PubMedView ArticleGoogle Scholar
- Maksymowych AB, Reinhard M, Malizio CJ, Goodnough MC, Johnson EA, Simpson LL: Pure botulinum neurotoxin is absorbed from the stomach and small intestine and produces peripheral neuromuscular blockade. Infect Immun. 1999, 67: 4708-4712.PubMed CentralPubMedGoogle Scholar
- Amimoto K, Sasaki Y, Fukuyama S, Tamura Y: Genetic variation and cross-reactivity of Clostridium septicum alpha-toxin. Vet Microbiol. 2006, 114: 51-59. 10.1016/j.vetmic.2005.10.039.PubMedView ArticleGoogle Scholar
- Manabe S, Nariya H, Miyata S, Tanaka H, Minami J, Suzuki M, Taniguchi Y, Okabe A: Purification and characterization of a clostripain-like protease from a recombinant Clostridium perfringens culture. Microbiology. 2010, 156: 561-569. 10.1099/mic.0.031609-0.PubMedView ArticleGoogle Scholar
- Sterthoff C, Lang AE, Schwan C, Tauch A, Aktories K: Functional characterization of an extended binding component of the actin-ADP-ribosylating C2 toxin detected in Clostridium botulinum strain (C) 2300. Infect Immun. 2010, 78: 1468-1474. 10.1128/IAI.01351-09.PubMed CentralPubMedView ArticleGoogle Scholar
- Myers GS, Rasko DA, Cheung JK, Ravel J, Seshadri R, DeBoy RT, Ren Q, Varga J, Awad MM, Brinkac LM, et al: Skewed genomic variability in strains of the toxigenic bacterial pathogen, Clostridium perfringens. Genome Res. 2006, 16: 1031-1040. 10.1101/gr.5238106.PubMed CentralPubMedView ArticleGoogle Scholar
- Liyanage H, Holcroft P, Evans VJ, Keis S, Wilkinson SR, Kashket ER, Young M: A new insertion sequence, ISCb1, from Clostridium beijernickii NCIMB 8052. J Mol Microbiol Biotechnol. 2000, 2: 107-113.PubMedGoogle Scholar
- Hill KK, Xie G, Foley BT, Smith TJ, Munk AC, Bruce D, Smith LA, Brettin TS, Detter JC: Recombination and insertion events involving the botulinum neurotoxin complex genes in Clostridium botulinum types A, B, E and F and Clostridium butyricum type E strains. BMC Biol. 2009, 7: 66-10.1186/1741-7007-7-66.PubMed CentralPubMedView ArticleGoogle Scholar
- Schmid S, Seitz T, Haas D: Cointegrase, a naturally occurring, truncated form of IS21 transposase, catalyzes replicon fusion rather than simple insertion of IS21. J Mol Biol. 1998, 282: 571-583. 10.1006/jmbi.1998.2041.PubMedView ArticleGoogle Scholar
- Song H, Hwang J, Yi H, Ulrich RL, Yu Y, Nierman WC, Kim HS: The early stage of bacterial genome-reductive evolution in the host. PLoS Pathog. 2010, 6: e1000922-10.1371/journal.ppat.1000922.PubMed CentralPubMedView ArticleGoogle Scholar
- Gordon D, Abajian C, Green P: Consed: a graphical tool for sequence finishing. Genome Res. 1998, 8: 195-202.PubMedView ArticleGoogle Scholar
- Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, Salzberg SL: Versatile and open software for comparing large genomes. Genome Biol. 2004, 5: R12-10.1186/gb-2004-5-2-r12.PubMed CentralPubMedView ArticleGoogle Scholar
- Carver TJ, Rutherford KM, Berriman M, Rajandream MA, Barrell BG, Parkhill J: ACT: the Artemis Comparison Tool. Bioinformatics. 2005, 21: 3422-3423. 10.1093/bioinformatics/bti553.PubMedView ArticleGoogle Scholar
- Barton BM, Harding GP, Zuccarelli AJ: A general method for detecting and sizing large plasmids. Anal Biochem. 1995, 226: 235-240. 10.1006/abio.1995.1220.PubMedView ArticleGoogle Scholar
- Delcher AL, Harmon D, Kasif S, White O, Salzberg SL: Improved microbial gene identification with GLIMMER. Nucleic Acids Res. 1999, 27: 4636-4641. 10.1093/nar/27.23.4636.PubMed CentralPubMedView ArticleGoogle Scholar
- Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, Barrell B: Artemis: sequence visualization and annotation. Bioinformatics. 2000, 16: 944-945. 10.1093/bioinformatics/16.10.944.PubMedView ArticleGoogle Scholar
- Carver T, Thomson N, Bleasby A, Berriman M, Parkhill J: DNAPlotter: circular and linear interactive genome visualization. Bioinformatics. 2009, 25: 119-120. 10.1093/bioinformatics/btn578.PubMed CentralPubMedView ArticleGoogle Scholar
- Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964. 10.1093/nar/25.5.955.PubMed CentralPubMedView ArticleGoogle Scholar
- Grissa I, Vergnaud G, Pourcel C: CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2007, 35: W52-57. 10.1093/nar/gkm360.PubMed CentralPubMedView ArticleGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.