Genetic diversity among Clostridium botulinum strains harboring bont/A2 and bont/A3 genes

Exploring the genetic background of the botulism neurotoxin BoNT/B2 in Spain

To determine whether the neurotoxin BoNT/B2 causing botulism in Spain is clonal, the genetic diversity and phylogenetic relationships of Clostridium botulinum from food-borne episodes and infant cases of the condition were explored. The botulinum toxin gene (bont) subtype, the variable region of the flagellin gene (flaVR), and a seven-gene multi-locus sequence type were examined by sequencing 37 BoNT-positive cultures obtained over the period 2010 to 2022. Out of 37 botulism events, 16 food-borne episodes and 16 infant cases were associated with bont/b2. Eight bont/b2 alleles were detected [nucleotide distance range 0.0259-0.415%, Hunter and Gaston discrimination index (HGDI) 0.71]. The most common bont/b2 allele corresponded to that of strain Prevot 25 NCASE and its single and double locus variations (87.5%). Four known flaVR types were identified (HGDI 0.79), along with one previously unknown (flaVR-15). Sixteen sequence types (STs) (HGDI 0.89) were recorded including seven new STs (ST164-ST170; 10 new alleles) and five new STs (ST171-ST175; with new allele combinations) were also noted. Correlations among some STs and flaVR types were seen. Overall, the present results show that the combined analysis of bont/b2-flaVR-ST at the nucleotide level could be used to track botulism events in Spain. The neurotoxin BoNT/B2 has largely been responsible for human botulism in Spain. The polymorphism analysis of bont/b2, flaVR typing, and sequence type determinations, revealed a wide variety of clones to be responsible for human botulism, ruling out a common source of acquisition. IMPORTANCE Botulism, a potentially fatal disease, is classically characterized by a symmetrical descending flaccid paralysis, which if left untreated can lead to respiratory failure and death. Botulinum neurotoxin (BoNT), produced by certain species of Clostridium, is the most potent biological toxin known, and the direct cause of botulism. This study characterizes the acquisition in Spain of two forms of botulism, i.e., food-borne and infant botulism, which are largely caused by the main neurotoxin BoNT/B2. Polymorphism analysis of the bont/b2 gene, typing of the flagellin variable region sequence (flaVR), and multilocus sequence typing, were used to explore the genetic background of Clostridium botulinum group I. To our knowledge, this is the first phylogenetic and typing study of botulism undertaken in Spain.

The Distinctive Evolution of orfX Clostridium parabotulinum Strains and Their Botulinum Neurotoxin Type A and F Gene Clusters Is Influenced by Environmental Factors and Gene Interactions via Mobile Genetic Elements

Of the seven currently known botulinum neurotoxin-producing species of Clostridium, C. parabotulinum, or C. botulinum Group I, is the species associated with the majority of human botulism cases worldwide. Phylogenetic analysis of these bacteria reveals a diverse species with multiple genomic clades. The neurotoxins they produce are also diverse, with over 20 subtypes currently represented. The existence of different bont genes within very similar genomes and of the same bont genes/gene clusters within different bacterial variants/species indicates that they have evolved independently. The neurotoxin genes are associated with one of two toxin gene cluster types containing either hemagglutinin (ha) genes or orfX genes. These genes may be located within the chromosome or extrachromosomal elements such as large plasmids. Although BoNT-producing C parabotulinum bacteria are distributed globally, they are more ubiquitous in certain specific geographic regions. Notably, northern hemisphere strains primarily contain ha gene clusters while southern hemisphere strains have a preponderance of orfX gene clusters. OrfX C. parabotulinum strains constitute a subset of this species that contain highly conserved bont gene clusters having a diverse range of bont genes. While much has been written about strains with ha gene clusters, less attention has been devoted to those with orfX gene clusters. The recent sequencing of 28 orfX C. parabotulinum strains and the availability of an additional 91 strains for analysis provides an opportunity to compare genomic relationships and identify unique toxin gene cluster characteristics and locations within this species subset in depth. The mechanisms behind the independent processes of bacteria evolution and generation of toxin diversity are explored through the examination of bacterial relationships relating to source locations and evidence of horizontal transfer of genetic material among different bacterial variants, particularly concerning bont gene clusters. Analysis of the content and locations of the bont gene clusters offers insights into common mechanisms of genetic transfer, chromosomal integration, and development of diversity among these genes.

Genomic Characterization of Newly Completed Genomes of Botulinum Neurotoxin-Producing Species from Argentina, Australia, and Africa

Botulinum neurotoxin-producing clostridia are diverse in the types of toxins they produce as well as in their overall genomic composition. They are globally distributed, with prevalent species and toxin types found within distinct geographic regions, but related strains containing the same toxin types may also be located on distinct continents. The mechanisms behind the spread of these bacteria and the independent movements of their bont genes may be understood through examination of their genetic backgrounds. The generation of 15 complete genomic sequences from bacteria isolated in Argentina, Australia, and Africa allows for a thorough examination of genome features, including overall relationships, bont gene cluster locations and arrangements, and plasmid comparisons, in bacteria isolated from various areas in the southern hemisphere. Insights gained from these examinations provide an understanding of the mechanisms behind the independent movements of these elements among distinct species.

Diversity of the Genomes and Neurotoxins of Strains of Clostridium botulinum Group I and Clostridium sporogenes Associated with Foodborne, Infant and Wound Botulism

Clostridium botulinum Group I and Clostridium sporogenes are closely related bacteria responsible for foodborne, infant and wound botulism. A comparative genomic study with 556 highly diverse strains of C. botulinum Group I and C. sporogenes (including 417 newly sequenced strains) has been carried out to characterise the genetic diversity and spread of these bacteria and their neurotoxin genes. Core genome single-nucleotide polymorphism (SNP) analysis revealed two major lineages; C. botulinum Group I (most strains possessed botulinum neurotoxin gene(s) of types A, B and/or F) and C. sporogenes (some strains possessed a type B botulinum neurotoxin gene). Both lineages contained strains responsible for foodborne, infant and wound botulism. A new C. sporogenes cluster was identified that included five strains with a gene encoding botulinum neurotoxin sub-type B1. There was significant evidence of horizontal transfer of botulinum neurotoxin genes between distantly related bacteria. Population structure/diversity have been characterised, and novel associations discovered between whole genome lineage, botulinum neurotoxin sub-type variant, epidemiological links to foodborne, infant and wound botulism, and geographic origin. The impact of genomic and physiological variability on the botulism risk has been assessed. The genome sequences are a valuable resource for future research (e.g., pathogen biology, evolution of C. botulinum and its neurotoxin genes, improved pathogen detection and discrimination), and support enhanced risk assessments and the prevention of botulism.

The Light Chain Defines the Duration of Action of Botulinum Toxin Serotype A Subtypes

Botulinum neurotoxin (BoNT) is the causative agent of botulism and a widely used pharmaceutical to treat a variety of neurological diseases. BoNTs are 150-kDa protein toxins organized into heavy chain (HC) and light chain (LC) domains linked by a disulfide bond. The HC selectively binds to neurons and aids cell entry of the enzymatically active LC. There are seven immunological BoNT serotypes (A to G); each serotype includes genetic variants, termed subtypes. Only two subtypes, BoNT/A1 and BoNT/B1, are currently used as therapeutics. BoNT serotype A (BoNT/A) subtypes A2 to A8 show distinct potency, duration of action, and pathology relative to BoNT/A1. Specifically, BoNT/A3 possesses shorter duration of action and elicits distinct symptoms in mice at high toxin doses. In this report, we analyzed the roles of LC and HC of BoNT/A3 for duration of action, neuronal cell entry, and mouse pathology by using clostridium-derived recombinant hybrid BoNTs consisting of reciprocal LC and HC (BoNTA1/A3 and BoNTA3/A1). Hybrid toxins were processed in their expression host to a dichain BoNT consisting of LC and HC linked via a disulfide bond. The LC and HC defined BoNT potency in mice and BoNT toxicity for cultured neuronal cells, while the LC defined the duration of BoNT action in cell and mouse models. Protein alignment identified a previously unrecognized region within the LC subtype A3 (LC/A3) relative to the other LC serotype A (LC/A) subtypes (low primary acid homology [LPH]) that correlated to intracellular LC localization. This study shows the utility of recombinant hybrid BoNTs with new therapeutic potential, while remaining sensitive to antitoxins and therapies to native BoNT.

Botulinum neurotoxins are the most potent protein toxins for humans and potential bioterrorism threats, but they are also widely used as pharmaceuticals. Within the large family of BoNTs, only two subtypes are currently used as pharmaceuticals, with a large number of BoNT subtypes remaining as untapped potential sources for unique pharmaceuticals. Here, two recombinant hybrid toxins were engineered, consisting of domains from two BoNT subtypes that possess distinct duration of action and activity in human neurons and mice. We define the functional domains responsible for BoNT action and demonstrate creation of functional hybrid BoNTs with new therapeutic potential, while remaining sensitive to antitoxins and therapies to native BoNT.

Reoccurrence of botulinum neurotoxin subtype A3 inducing food-borne botulism, Slovakia, 2015

A case of food-borne botulism occurred in Slovakia in 2015. Clostridium botulinum type A was isolated from three nearly empty commercial hummus tubes. The product, which was sold in Slovakia and the Czech Republic, was withdrawn from the market and a warning was issued immediately through the European Commission's Rapid Alert System for Food and Feed (RASFF). Further investigation revealed the presence of botulinum neurotoxin (BoNT) subtype BoNT/A3, a very rare subtype implicated in only one previous outbreak (Loch Maree in Scotland, 1922). It is the most divergent subtype of BoNT/A with 15.4% difference at the amino acid level compared with the prototype BoNT/A1. This makes it more prone to evading immunological and PCR-based detection. It is recommended that testing laboratories are advised that this subtype has been associated with food-borne botulism for the second time since the first outbreak almost 100 years ago, and to validate their immunological or PCR-based methods against this divergent subtype.

Pulsotype Diversity of Clostridium botulinum Strains Containing Serotypes A and/or B Genes

Clostridium botulinum strains are prevalent in the environment and produce a potent neurotoxin that causes botulism, a rare but serious paralytic disease. In 2010, a national PulseNet database was established to curate C. botulinum pulsotypes and facilitate epidemiological investigations, particularly for serotypes A and B strains frequently associated with botulism cases in the United States. Between 2010 and 2014 we performed pulsed-field gel electrophoresis (PFGE) using a PulseNet protocol, uploaded the resulting PFGE patterns into a national database, and analyzed data according to PulseNet criteria (UPGMA clustering, Dice coefficient, 1.5% position tolerance, and 1.5% optimization). A retrospective data analysis was undertaken on 349 entries comprised of type A and B strains isolated from foodborne and infant cases to determine epidemiological relevance, resolution of the method, and the diversity of the database. Most studies to date on the pulsotype diversity of C. botulinum have encompassed very small sets of isolates; this study, with over 300 isolates, is more comprehensive than any published to date. Epidemiologically linked isolates had indistinguishable patterns, except in four instances and there were no obvious geographic trends noted. Simpson's Index of Diversity (D) has historically been used to demonstrate species diversity and abundance within a group, and is considered a standard descriptor for PFGE databases. Simpson's Index was calculated for each restriction endonuclease (SmaI, XhoI), the pattern combination SmaI-XhoI, as well as for each toxin serotype. The D values indicate that both enzymes provided better resolution for serotype B isolates than serotype A. XhoI as the secondary enzyme provided little additional discrimination for C. botulinum. SmaI patterns can be used to exclude unrelated isolates during a foodborne outbreak, but pulsotypes should always be considered concurrently with available epidemiological data.

Diversity of Group I and II Clostridium botulinum Strains from France Including Recently Identified Subtypes

In France, human botulism is mainly food-borne intoxication, whereas infant botulism is rare. A total of 99 group I and II Clostridium botulinum strains including 59 type A (12 historical isolates [1947-1961], 43 from France [1986-2013], 3 from other countries, and 1 collection strain), 31 type B (3 historical, 23 recent isolates, 4 from other countries, and 1 collection strain), and 9 type E (5 historical, 3 isolates, and 1 collection strain) were investigated by botulinum locus gene sequencing and multilocus sequence typing analysis. Historical C. botulinum A strains mainly belonged to subtype A1 and sequence type (ST) 1, whereas recent strains exhibited a wide genetic diversity: subtype A1 in orfX or ha locus, A1(B), A1(F), A2, A2b2, A5(B2') A5(B3'), as well as the recently identified A7 and A8 subtypes, and were distributed into 25 STs. Clostridium botulinum A1(B) was the most frequent subtype from food-borne botulism and food. Group I C. botulinum type B in France were mainly subtype B2 (14 out of 20 historical and recent strains) and were divided into 19 STs. Food-borne botulism resulting from ham consumption during the recent period was due to group II C. botulinum B4. Type E botulism is rare in France, 5 historical and 1 recent strains were subtype E3. A subtype E12 was recently identified from an unusual ham contamination. Clostridium botulinum strains from human botulism in France showed a wide genetic diversity and seems to result not from a single evolutionary lineage but from multiple and independent genetic rearrangements.

Genomic Epidemiology of Clostridium botulinum Isolates from Temporally Related Cases of Infant Botulism in New South Wales, Australia

Infant botulism is a potentially life-threatening paralytic disease that can be associated with prolonged morbidity if not rapidly diagnosed and treated. Four infants were diagnosed and treated for infant botulism in NSW, Australia, between May 2011 and August 2013. Despite the temporal relationship between the cases, there was no close geographical clustering or other epidemiological links. Clostridium botulinum isolates, three of which produced botulism neurotoxin serotype A (BoNT/A) and one BoNT serotype B (BoNT/B), were characterized using whole-genome sequencing (WGS). In silico multilocus sequence typing (MLST) found that two of the BoNT/A-producing isolates shared an identical novel sequence type, ST84. The other two isolates were single-locus variants of this sequence type (ST85 and ST86). All BoNT/A-producing isolates contained the same chromosomally integrated BoNT/A2 neurotoxin gene cluster. The BoNT/B-producing isolate carried a single plasmid-borne bont/B gene cluster, encoding BoNT subtype B6. Single nucleotide polymorphism (SNP)-based typing results corresponded well with MLST; however, the extra resolution provided by the whole-genome SNP comparisons showed that the isolates differed from each other by >3,500 SNPs. WGS analyses indicated that the four infant botulism cases were caused by genomically distinct strains of C. botulinum that were unlikely to have originated from a common environmental source. The isolates did, however, cluster together, compared with international isolates, suggesting that C. botulinum from environmental reservoirs throughout NSW have descended from a common ancestor. Analyses showed that the high resolution of WGS provided important phylogenetic information that would not be captured by standard seven-loci MLST.

Historical and current perspectives on Clostridium botulinum diversity

For nearly one hundred years, researchers have attempted to categorize botulinum neurotoxin-producing clostridia and the toxins that they produce according to biochemical characterizations, serological comparisons, and genetic analyses. Throughout this period the bacteria and their toxins have defied such attempts at categorization. Below is a description of both historic and current Clostridium botulinum strain and neurotoxin information that illustrates how each new finding has significantly added to the knowledge of the botulinum neurotoxin-containing clostridia and their diversity.

Molecular subtyping of Clostridium botulinum by pulsed-field gel electrophoresis

Pulsed-field gel electrophoresis (PFGE) has been extensively used to estimate the genetic diversity of Clostridium botulinum. In addition, PFGE is the standard method for investigating foodborne outbreaks associated with various enteric pathogens, including C. botulinum. PFGE can be used to exclude a suspected but not confirmed food source when the patterns of the food and clinical isolates are different. Indistinguishable PFGE patterns may also be useful for linking isolates between patients or to a food source, but results must be interpreted within an epidemiological context to ensure isolates are truly related. Here, we describe a standardized laboratory protocol for molecular subtyping of C. botulinum by PFGE.

First report of an infant botulism case due to Clostridium botulinum type Af

Most infant botulism cases worldwide are due to botulinum toxin types A and B. Rarely, Clostridium botulinum strains that produce two serotypes (Ab, Ba, and Bf) have also been isolated from infant botulism cases. This is the first reported case of infant botulism due to C. botulinum type Af worldwide.

Genomes, neurotoxins and biology of Clostridium botulinum Group I and Group II

Recent developments in whole genome sequencing have made a substantial contribution to understanding the genomes, neurotoxins and biology of Clostridium botulinum Group I (proteolytic C. botulinum) and C. botulinum Group II (non-proteolytic C. botulinum). Two different approaches are used to study genomics in these bacteria; comparative whole genome microarrays and direct comparison of complete genome DNA sequences. The properties of the different types of neurotoxin formed, and different neurotoxin gene clusters found in C. botulinum Groups I and II are explored. Specific examples of botulinum neurotoxin genes are chosen for an in-depth discussion of neurotoxin gene evolution. The most recent cases of foodborne botulism are summarised.

Genetic characterization and comparison of Clostridium botulinum isolates from botulism cases in Japan between 2006 and 2011

Genetic characterization was performed for 10 group I Clostridium botulinum strains isolated from botulism cases in Japan between 2006 and 2011. Of these, 1 was type A, 2 were type B, and 7 were type A(B) {carrying a silent bont/B [bont/(B)] gene} serotype strains, based on botulinum neurotoxin (BoNT) production. The type A strain harbored the subtype A1 BoNT gene (bont/A1), which is associated with the ha gene cluster. The type B strains carried bont/B5 or bont/B6 subtype genes. The type A(B) strains carried bont/A1 identical to that of type A(B) strain NCTC2916. However, bont/(B) genes in these strains showed single-nucleotide polymorphisms (SNPs) among strains. SNPs at 2 nucleotide positions of bont/(B) enabled classification of the type A(B) strains into 3 groups. Pulsed-field gel electrophoresis (PFGE) and multiple-locus variable-number tandem-repeat analysis (MLVA) also provided consistent separation results. In addition, the type A(B) strains were separated into 2 lineages based on their plasmid profiles. One lineage carried a small plasmid (5.9 kb), and another harbored 21-kb plasmids. To obtain more detailed genetic information about the 10 strains, we sequenced their genomes and compared them with 13 group I C. botulinum genomes in a database using whole-genome SNP analysis. This analysis provided high-resolution strain discrimination and enabled us to generate a refined phylogenetic tree that provides effective traceability of botulism cases, as well as bioterrorism materials. In the phylogenetic tree, the subtype B6 strains, Okayama2011 and Osaka05, were distantly separated from the other strains, indicating genomic divergence of subtype B6 strains among group I strains.

Comparison of assembled Clostridium botulinum A1 genomes revealed their evolutionary relationship

Clostridium botulinum encompasses bacteria that produce at least one of the seven serotypes of botulinum neurotoxin (BoNT/A-G). The availability of genome sequences of four closely related Type A1 or A1(B) strains, as well as the A1-specific microarray, allowed the analysis of their genomic organizations and evolutionary relationship. The four genomes share >90% core genes and >96% functional groups. Phylogenetic analysis based on COG shows closer relations of the A1(B) strain, NCTC 2916, to B1 and F1 than A1 strains. Alignment of the genomes of the three A1 strains revealed a highly similar chromosomal structure with three small gaps in the genome of ATCC 19397 and one additional gap in the genome of Hall A, suggesting ATCC 19379 as an evolutionary intermediate between Hall A and ATCC 3502. Analyses of the four gap regions indicated potential horizontal gene transfer and recombination events important for the evolution of A1 strains.

Analysis of the genetic distribution among members of Clostridium botulinum group I using a novel multilocus sequence typing (MLST) assay

Clostridium botulinum is the etiological agent of botulism. Due to food-borne poisoning and the potential use of the extremely toxic botulinum neurotoxin (BoNT) from C. botulinum in bioterror or biocrime related actions, reliable high resolution typing methods for discriminating C. botulinum strains are needed. Partial sequencing of the adk, atpH, gyrB, proC, rpoD and spo0A genes from 51 various C. botulinum/sporogenes isolates was performed, resulting in 37 different sequence types (STs). Analysis of the sequence data revealed a genetic distribution in five larger clusters with a loose correlation to the BoNT serotypes. The developed MLST assay had a slightly lower resolution ability when compared to the MLVA (multilocus variable number of tandem repeat analysis), but the two methods resulted in similar subclusters of the strains possessing the BoNT serotypes A, B and F. The current work presents the development of a novel MLST assay useful for genotyping C. botulinum related to basic phylogenetic research and trace-back analysis in microbial forensic studies.

Multi-locus variable number tandem repeat analysis for Clostridium botulinum type B isolates in Japan: comparison with other isolates and genotyping methods

Clostridium botulinum produces botulinum neurotoxin (BoNT) and causes botulism in humans and animals. Recently, 15-loci multi-locus variable number tandem repeat analysis (MLVA) for C. botulinum was developed for high-resolution and inter-lab comparative genotyping. This study examines the relation between MLVA and other genotyping methods such as pulsed-field gel electrophoresis (PFGE), multi-locus sequence typing (MLST), BoNT/B subtyping and bont/b gene location to evaluate MLVA as a method applicable to the genetic markers for C. botulinum type B. Japanese isolates were genotyped using MLVA and were compared with strains from other sources reported previously. Results show that the discriminatory power of MLVA was comparable to that of PFGE and higher than that of MLST. The topology of the minimum spanning tree (MST) constructed using MLVA data was very consistent with the phylogenetic classifications of PFGE and MLST. The MST topology also represented genetic diversity between the strains possessing bont/b gene on chromosomes and plasmids. Some Japanese isolates including those associated with infant botulism were inferred to be related to isolates of Europe origin from MLVA genotyping results. The MLVA scheme used for this study is apparently useful not only for high-resolution molecular typing, but also for phylogenetic characterization of C. botulinum type B.