Evolution of Chromosomal Clostridium botulinum Type E Neurotoxin Gene Clusters: Evidence Provided by Their Rare Plasmid-Borne Counterparts

The Light Chain Domain and Especially the C-Terminus of Receptor-Binding Domain of the Botulinum Neurotoxin (BoNT) Are the Hotspots for Amino Acid Variability and Toxin Type Diversity

Botulinum neurotoxins (BoNT) are the most potent toxins in the world. They are produced by a few dozens of strains within several clostridial species. The toxin that they produce can cause botulism, a flaccid paralysis in humans and other animals. With seven established serologically different types and over 40 subtypes, BoNTs are among the most diverse known toxins. The toxin, its structure, its function and its physiological effects on the neural cell and animal hosts along with its diversity have been the subjects of numerous studies. However, many gaps remain in our knowledge about the BoNT toxin and the species that produce them. One of these gaps involves the distribution and extent of variability along the full length of the gene and the protein as well as its domains and subdomains. In this study, we performed an extensive analysis of all of the available 143 unique BoNT-encoding genes and their products, and we investigated their diversity and evolution. Our results indicate that while the nucleotide variability is almost uniformly distributed along the entire length of the gene, the amino acid variability is not. We found that most of the differences were concentrated along the protein's light chain (LC) domain and especially, the C-terminus of the receptor-binding domain (H_CC). These two regions of the protein are thus identified as the main source of the toxin type differentiation, and consequently, this toxin's versatility to bind different receptors and their isoforms and act upon different substrates, thus infecting different hosts.

Endogenous CRISPR-Cas Systems in Group I Clostridium botulinum and Clostridium sporogenes Do Not Directly Target the Botulinum Neurotoxin Gene Cluster

Most strains of proteolytic group I Clostridium botulinum (G1 C. botulinum) and some strains of Clostridium sporogenes possess genes encoding botulinum neurotoxin (BoNT), a potent neuroparalytic agent. Within G1 C. botulinum, conserved bont gene clusters of three major toxin serotypes (bont/A/B/F) can be found on conjugative plasmids and/or within chromosomal pathogenicity islands. CRISPR-Cas systems enable site-specific targeting of previously encountered mobile genetic elements (MGE) such as plasmids and bacteriophage through the creation of a spacer library complementary to protospacers within the MGEs. To examine whether endogenous CRISPR-Cas systems restrict the transfer of bont gene clusters across strains we conducted a bioinformatic analysis profiling endogenous CRISPR-Cas systems from 241 G1 C. botulinum and C. sporogenes strains. Approximately 6,200 CRISPR spacers were identified across the strains and Type I-B, III-A/B/D cas genes and CRISPR array features were identified in 83% of the strains. Mapping the predicted spacers against the masked strain and RefSeq plasmid dataset identified 56,000 spacer-protospacer matches. While spacers mapped heavily to targets within bont(+) plasmids, no protospacers were identified within the bont gene clusters. These results indicate the toxin is not a direct target of CRISPR-Cas but the plasmids predominantly responsible for its mobilization are. Finally, while the presence of a CRISPR-Cas system did not reliably indicate the presence or absence of a bont gene cluster, comparative genomics across strains indicates they often occupy the same hypervariable loci common to both species, potentially suggesting similar mechanisms are involved in the acquisition and curation of both genomic features.

Novel Putative Transposable Element Associated with the Subtype E5 Botulinum Toxin Gene Cluster of Neurotoxigenic Clostridium butyricum Type E Strains from China

Previously, a whole-genome comparison of three Clostridium butyricum type E strains from Italy and the United States with different C. botulinum type E strains indicated that the bont/e gene might be transferred between the two clostridia species through transposition. However, transposable elements (TEs) have never been identified close to the bont/e gene. Herein, we report the whole genome sequences for four neurotoxigenic C. butyricum type E strains that originated in China. An analysis of the obtained genome sequences revealed the presence of a novel putative TE upstream of the bont/e gene in the genome of all four strains. Two strains of environmental origin possessed an additional copy of the putative TE in their megaplasmid. Similar putative TEs were found in the megaplasmids and, less frequently, in the chromosomes of several C. butyricum strains, of which two were neurotoxigenic C. butyricum type E strains, and in the chromosome of a single C. botulinum type E strain. We speculate that the putative TE might potentially transpose the bont/e gene at the intracellular and inter-cellular levels. However, the occasional TE occurrence in the clostridia genomes might reflect rare transposition events.

Integration of Complete Plasmids Containing Bont Genes into Chromosomes of Clostridium parabotulinum, Clostridium sporogenes, and Clostridium argentinense

At least 40 toxin subtypes of botulinum neurotoxins (BoNTs), a heterogenous group of bacterial proteins, are produced by seven different clostridial species. A key factor that drives the diversity of neurotoxigenic clostridia is the association of bont gene clusters with various genomic locations including plasmids, phages and the chromosome. Analysis of Clostridium sporogenes BoNT/B1 strain CDC 1632, C. argentinense BoNT/G strain CDC 2741, and Clostridium parabotulinum BoNT/B1 strain DFPST0006 genomes revealed bont gene clusters within plasmid-like sequences within the chromosome or nested in large contigs, with no evidence of extrachromosomal elements. A nucleotide sequence (255,474 bp) identified in CDC 1632 shared 99.5% identity (88% coverage) with bont/B1-containing plasmid pNPD7 of C. sporogenes CDC 67071; CDC 2741 contig AYSO01000020 (1.1 MB) contained a ~140 kb region which shared 99.99% identity (100% coverage) with plasmid pRSJ17_1 of C. argentinense BoNT/G strain 89G; and DFPST0006 contig JACBDK0100002 (573 kb) contained a region that shared 100% identity (99%) coverage with the bont/B1-containing plasmid pCLD of C. parabotulinum Okra. This is the first report of full-length plasmid DNA-carrying complete neurotoxin gene clusters integrated in three distinct neurotoxigenic species: C. parabotulinum, C. sporogenes and C. argentinense.

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.

Pan-Genomic Analysis of Clostridium botulinum Group II (Non-Proteolytic C. botulinum) Associated with Foodborne Botulism and Isolated from the Environment

The neurotoxin formed by Clostridium botulinum Group II is a major cause of foodborne botulism, a deadly intoxication. This study aims to understand the genetic diversity and spread of C. botulinum Group II strains and their neurotoxin genes. A comparative genomic study has been conducted with 208 highly diverse C. botulinum Group II strains (180 newly sequenced strains isolated from 16 countries over 80 years, 28 sequences from Genbank). Strains possessed a single type B, E, or F neurotoxin gene or were closely related strains with no neurotoxin gene. Botulinum neurotoxin subtype variants (including novel variants) with a unique amino acid sequence were identified. Core genome single-nucleotide polymorphism (SNP) analysis identified two major lineages-one with type E strains, and the second dominated by subtype B4 strains with subtype F6 strains. This study revealed novel details of population structure/diversity and established relationships between whole-genome lineage, botulinum neurotoxin subtype variant, association with foodborne botulism, epidemiology, and geographical source. Additionally, the genome sequences represent a valuable resource for the research community (e.g., understanding evolution of C. botulinum and its neurotoxin genes, dissecting key aspects of C. botulinum Group II biology). This may contribute to improved risk assessments and the prevention of foodborne botulism.

Virulence Plasmids of the Pathogenic Clostridia

The clostridia cause a spectrum of diseases in humans and animals ranging from life-threatening tetanus and botulism, uterine infections, histotoxic infections and enteric diseases, including antibiotic-associated diarrhea, and food poisoning. The symptoms of all these diseases are the result of potent protein toxins produced by these organisms. These toxins are diverse, ranging from a multitude of pore-forming toxins to phospholipases, metalloproteases, ADP-ribosyltransferases and large glycosyltransferases. The location of the toxin genes is the unifying theme of this review because with one or two exceptions they are all located on plasmids or on bacteriophage that replicate using a plasmid-like intermediate. Some of these plasmids are distantly related whilst others share little or no similarity. Many of these toxin plasmids have been shown to be conjugative. The mobile nature of these toxin genes gives a ready explanation of how clostridial toxin genes have been so widely disseminated both within the clostridial genera as well as in the wider bacterial community.

Historical Perspectives and Guidelines for Botulinum Neurotoxin Subtype Nomenclature

Botulinum neurotoxins are diverse proteins. They are currently represented by at least seven serotypes and more than 40 subtypes. New clostridial strains that produce novel neurotoxin variants are being identified with increasing frequency, which presents challenges when organizing the nomenclature surrounding these neurotoxins. Worldwide, researchers are faced with the possibility that toxins having identical sequences may be given different designations or novel toxins having unique sequences may be given the same designations on publication. In order to minimize these problems, an ad hoc committee consisting of over 20 researchers in the field of botulinum neurotoxin research was convened to discuss the clarification of the issues involved in botulinum neurotoxin nomenclature. This publication presents a historical overview of the issues and provides guidelines for botulinum neurotoxin subtype nomenclature in the future.

Diversity of the Germination Apparatus in Clostridium botulinum Groups I, II, III, and IV

Clostridium botulinum is a highly dangerous pathogen that forms very resistant endospores that are ubiquitous in the environment, and which, under favorable conditions germinate to produce vegetative cells that multiply and form the exceptionally potent botulinum neurotoxin. To improve the control of botulinum neurotoxin-forming clostridia, it is important to understand the mechanisms involved in spore germination. Here we present models for spore germination in C. botulinum based on comparative genomics analyses, with C. botulinum Groups I and III sharing similar pathways, which differ from those proposed for C. botulinum Groups II and IV. All spores germinate in response to amino acids interacting with a germinant receptor, with four types of germinant receptor identified [encoded by various combinations of gerA, gerB, and gerC genes (gerX)]. There are three gene clusters with an ABC-like configuration; ABC [gerX1], ABABCB [gerX2] and ACxBBB [gerX4], and a single CA-B [gerX3] gene cluster. Subtypes have been identified for most germinant receptor types, and the individual GerX subunits of each cluster show similar grouping in phylogenetic trees. C. botulinum Group I contained the largest variety of gerX subtypes, with three gerX1, three gerX2, and one gerX3 subtypes, while C. botulinum Group III contained two gerX1 types and one gerX4. C. botulinum Groups II and IV contained a single germinant receptor, gerX3 and gerX1, respectively. It is likely that all four C. botulinum Groups include a SpoVA channel involved in dipicolinic acid release. The cortex-lytic enzymes present in C. botulinum Groups I and III appear to be CwlJ and SleB, while in C. botulinum Groups II and IV, SleC appears to be important.

Impact of Clostridium botulinum genomic diversity on food safety

•

C. botulinum Groups I and II form botulinum neurotoxin and cause foodborne botulism.

•

Increased knowledge of C. botulinum Group I and II genomes and neurotoxin diversity.

•

Impact on food safety via improved surveillance and tracing/tracking during outbreaks.

•

New insights into C. botulinum biology, food chain transmission, evolution.

The deadly botulinum neurotoxin formed by Clostridium botulinum is the causative agent of foodborne botulism. The increasing availability of C. botulinum genome sequences is starting to allow the genomic diversity of C. botulinum Groups I and II and their neurotoxins to be characterised. This information will impact on microbiological food safety through improved surveillance and tracing/tracking during outbreaks, and a better characterisation of C. botulinum Groups I and II, including the risk presented, and new insights into their biology, food chain transmission, and evolution.