An Integrative Approach to Computational Modelling of the Gene Regulatory Network Controlling Clostridium botulinum Type A1 Toxin Production.

Clostridium botulinum produces botulinum neurotoxins (BoNTs), highly potent substances responsible for botulism. Currently, mathematical models of C. botulinum growth and toxigenesis are largely aimed at risk assessment and do not include explicit genetic information beyond group level but integrate many component processes, such as signalling, membrane permeability and metabolic activity. In this paper we present a scheme for modelling neurotoxin production in C. botulinum Group I type A1, based on the integration of diverse information coming from experimental results available in the literature. Experiments show that production of BoNTs depends on the growth-phase and is under the control of positive and negative regulatory elements at the intracellular level. Toxins are released as large protein complexes and are associated with non-toxic components. Here, we systematically review and integrate those regulatory elements previously described in the literature for C. botulinum Group I type A1 into a population dynamics model, to build the very first computational model of toxin production at the molecular level. We conduct a validation of our model against several items of published experimental data for different wild type and mutant strains of C. botulinum Group I type A1. The result of this process underscores the potential of mathematical modelling at the cellular level, as a means of creating opportunities in developing new strategies that could be used to prevent botulism; and potentially contribute to improved methods for the production of toxin that is used for therapeutics.

Distinguishing highly-related outbreak-associated Clostridium botulinum type A(B) strains.

BACKGROUND: In the United States, most Clostridium botulinum type A strains isolated during laboratory investigations of human botulism demonstrate the presence of an expressed type A botulinum neurotoxin (BoNT/A) gene and an unexpressed BoNT/B gene. These strains are designated type A(B). The most common pulsed-field gel electrophoresis (PFGE) pattern in the C. botulinum PulseNet database is composed of A(B) strains. The purpose of this study was to evaluate the ability of genome sequencing and multi-loci variable number of tandem repeat analysis (MLVA) to differentiate such strains. RESULTS: The genome sequences of type A(B) strains evaluated in this study are closely related and cluster together compared to other available C. botulinum Group I genomes. In silico multilocus sequence typing (MLST) analysis (7-loci) was unable to differentiate any of the type A(B) strains isolated from seven different outbreak investigations evaluated in this study. A 15-locus MLVA scheme demonstrated an improved ability to differentiate these strains, however, repeat unit variation among the strains was restricted to only two loci. Reference-free single nucleotide polymorphism (SNP) analysis demonstrated the ability to differentiate strains from all of the outbreaks examined and a non-outbreak associated strain. CONCLUSIONS: This study confirms that type A(B) strains that share the same PFGE pattern also share closely-related genome sequences. The lack of a complete type A(B) strain representative genome sequence hinders the ability to assemble genomes by reference mapping and analysis of SNPs at pre-identified sites. However, compared to other methods evaluated in this study, a reference-free SNP analysis demonstrated optimal subtyping utility for type A(B) strains using de novo assembled genome sequences.

Towards an international standard for detection and typing botulinum neurotoxin-producing Clostridia types A, B, E and F in food, feed and environmental samples: a European ring trial study to evaluate a real-time PCR assay.

A real-time PCR method for detection and typing of BoNT-producing Clostridia types A, B, E, and F was developed on the framework of the European Research Project "Biotracer". A primary evaluation was carried out using 104 strains and 17 clinical and food samples linked to botulism cases. Results showed 100% relative accuracy, 100% relative sensitivity, 100% relative specificity, and 100% selectivity (inclusivity on 73 strains and exclusivity on 31 strains) of the real-time PCR against the reference cultural method combined with the standard mouse bioassay. Furthermore, a ring trial study performed at four different European laboratories in Italy, France, the Netherlands, and Sweden was carried out using 47 strains, and 30 clinical and food samples linked to botulism cases. Results showed a concordance of 95.7% among the four laboratories. The reproducibility generated a relative standard deviation in the range of 2.18% to 13.61%. Considering the high level of agreement achieved between the laboratories, this real-time PCR is a suitable method for rapid detection and typing of BoNT-producing Clostridia in clinical, food and environmental samples and thus support the use of it as an international standard method.

Detection and differentiation of Clostridium botulinum type A strains using a focused DNA microarray.

A focused oligonucleotide microarray featuring 62 probes targeting strain variable regions of the Clostridium botulinum strain ATCC 3502 genome sequence was developed to differentiate C. botulinum type A strains. The strain variable regions were selected from deletions identified among a panel of 10 type A strains compared to the strain ATCC 3502 genome sequence using high density comparative genomic hybridization microarrays. The focused microarray also featured specific probes for the detection of the neurotoxin genes of various serotypes (A-G), toxin gene cluster components (ha70 and orfX1), and fldB as a marker for proteolytic clostridia (Group I). Eight pairs of strains selected from separate type A botulism outbreaks were included in the 27 subtype A1-A4 strains examined in this study. Each outbreak related strain pair consisted of strains isolated from different sources (stool and food). Of the eight outbreak related strain pairs, six groups of strains with indistinguishable hybridization patterns were identified. Outbreak related strains were shown to have identical hybridization patterns. Strain pairs from three separate outbreaks involving strains harboring both the type A neurotoxin gene (bont/A) and an unexpressed type B neurotoxin gene (bont/B) shared the same probe hybridization profile. The focused microarray format provides a rapid approach for neurotoxin gene detection and preliminary determination of the relatedness of strains isolated from different sources.

Phylogenetic analysis of Clostridium botulinum type A by multi-locus sequence typing.

The genus Clostridium comprises a heterogeneous group of organisms for which the phylogeny and evolutionary relationships are poorly understood. The elucidation of these evolutionary relationships necessitates the use of experimental methods that can distinguish Clostridium lineages that are time and cost effective, and can be accurately and reproducibly employed in different laboratories. Multi-locus sequence typing (MLST) has been successfully used as a reproducible and discriminating system in the study of eukaryotic and prokaryotic evolutionary biology, and for strain typing of various bacteria. In this study, MLST was applied to evaluate the evolutionary lineages in the serotype A group of Clostridium botulinum. C. botulinum type A has recently been shown to produce multiple subtypes, suggesting that it is not monophyletic as previously reported, but comprises distinct lineages. For MLST analysis, we initially evaluated 14 housekeeping genes (gapdh, tuf, sod, oppB, hsp60, dnaE, aroE, pta, 23S rDNA, aceK, rpoB, 16S rDNA, mdh and recA) for amplification and sequence analysis. In the first phase of the analysis, 30 C. botulinum type A strains producing botulinum neurotoxin subtypes A1-A4 were examined. Results of this pilot study suggested that seven of the genes (mdh, aceK, rpoB, aroE, hsp60, oppB and recA) could be used for elucidation of evolutionary lineages and strain typing. These seven housekeeping genes were successfully applied for the elucidation of lineages for 73 C. botulinum type A strains, which resulted in 24 distinct sequence types. This strategy should be applicable to phylogenetic studies and typing of other C. botulinum serotypes and Clostridium species.

Differentiation of Clostridium botulinum serotype A strains by multiple-locus variable-number tandem-repeat analysis.

Ten variable-number tandem-repeat (VNTR) regions identified within the complete genomic sequence of Clostridium botulinum strain ATCC 3502 were used to characterize 59 C. botulinum strains of the botulism neurotoxin A1 (BoNT/A1) to BoNT/A4 (BoNT/A1-A4) subtypes to determine their ability to discriminate among the serotype A strains. Two strains representing each of the C. botulinum serotypes B to G, including five bivalent strains, and two strains of the closely related species Clostridium sporogenes were also tested. Amplified fragment length polymorphism analyses revealed the genetic diversity among the serotypes and the high degree of similarity among many of the BoNT/A1 strains. The 10 VNTR markers amplified fragments within all of the serotype A strains but were less successful with strains of other serotypes. The composite multiple-locus VNTR analysis of the 59 BoNT/A1-A4 strains and 3 bivalent B strains identified 38 different genotypes. Thirty genotypes were identified among the 53 BoNT/A1 and BoNT/A1(B) strains, demonstrating discrimination below the subtype level. Contaminating DNA within crude toxin preparations of three BoNT/A subtypes (BoNT/A1 to BoNT/A3) also supported amplification of all of the VNTR regions. These markers provide clinical and forensics laboratories with a rapid, highly discriminatory tool to distinguish among C. botulinum BoNT/A1 strains for investigations of botulism outbreaks.

Plasmid encoded neurotoxin genes in Clostridium botulinum serotype A subtypes.

Clostridium botulinum, an important pathogen of humans and animals, produces botulinum neurotoxin (BoNT), the most poisonous toxin known. We have determined by pulsed-field gel electrophoresis (PFGE) and Southern hybridizations that the genes encoding BoNTs in strains Loch Maree (subtype A3) and 657Ba (type B and subtype A4) are located on large (approximately 280 kb) plasmids. This is the first demonstration of plasmid-borne neurotoxin genes in Clostridium botulinum serotypes A and B. The finding of BoNT type A and B genes on extrachromosomal elements has important implications for the evolution of neurotoxigenicity in clostridia including the origin, expression, and lateral transfer of botulinum neurotoxin genes.

A novel type A2 neurotoxin gene cluster in Clostridium botulinum strain Mascarpone.

The partial nucleotide sequence ( approximately 10 kb) of the cluster of genes encoding the botulinum neurotoxin complex in Clostridium botulinum type A strain Mascarpone was determined. The analysis revealed six ORFs (orfs), which were organized as in the type A2 and type A3 botulinum neurotoxin gene clusters of strains Kyoto-F and NCTC 2916, respectively. While the orfs at the proximal and distal ends of the sequence (orfX2 and bont/A genes) shared a high level of similarity with the corresponding sequences of strain Kyoto-F, the segment encompassing the orfX1 and botR/A genes within the sequence exhibited a higher degree of homology to the related region in strain NCTC 2916. The mosaic structure of the Mascarpone neurotoxin gene cluster suggests recombinational exchanges.

Genotyping Clostridium botulinum toxinotype A isolates from patients using amplified rDNA restriction analysis.

In this study, the application of amplified rDNA restriction analysis (ARDRA) for characterizing Clostridium botulinum toxinotype A strains isolated from individuals with botulism was evaluated. Ten restriction enzymes were tested for their suitability in ARDRA as a typing method and HhaI was selected for the best outcome. Analysis of HhaI restriction profiles of the amplified products divided C. botulinum isolates into three clusters. Non-toxigenic Clostridium sporogenes strains showed an ARDRA restriction pattern that was distinct from those observed for C. botulinum. The successful use of ARDRA for subdivision of C. botulinum in this study confirmed that this technique is a powerful method for typing of C. botulinum toxinotype A clonal diversity. In addition, it is rapid, sensitive and simple.

Differentiation of the gene clusters encoding botulinum neurotoxin type A complexes in Clostridium botulinum type A, Ab, and A(B) strains.

We describe a strategy to identify the clusters of genes encoding components of the botulinum toxin type A (boNT/A) complexes in 57 strains of Clostridium botulinum types A, Ab, and A(B) isolated in Italy and in the United States from different sources. Specifically, we combined the results of PCR for detecting the ha33 and/or p47 genes with those of boNT/A PCR-restriction fragment length polymorphism analysis. Three different type A toxin gene clusters were revealed; type A1 was predominant among the strains from the United States, whereas type A2 predominated among the Italian strains, suggesting a geographic distinction between strains. By contrast, no relationship between the toxin gene clusters and the clinical or food source of strains was evident. In two C. botulinum type A isolates from the United States, we recognized a third type A toxin gene cluster (designated type A3) which was similar to that previously described only for C. botulinum type A(B) and Ab strains. Total genomic DNA from the strains was subjected to pulsed-filed gel electrophoresis and randomly amplified polymorphic DNA analyses, and the results were consistent with the boNT/A gene clusters obtained.