Next Article in Journal
In Vitro and In Vivo Antifibrotic Effects of Fraxetin on Renal Interstitial Fibrosis via the ERK Signaling Pathway
Previous Article in Journal
Association of Sarcopenia and Gut Microbiota Composition in Older Patients with Advanced Chronic Kidney Disease, Investigation of the Interactions with Uremic Toxins, Inflammation and Oxidative Stress
Previous Article in Special Issue
Diversity of the Genomes and Neurotoxins of Strains of Clostridium botulinum Group I and Clostridium sporogenes Associated with Foodborne, Infant and Wound Botulism
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 13
Issue 7
10.3390/toxins13070473
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Theresa J. Smith
1,†,
Renmao Tian
2,†,
Behzad Imanian
2,3,
Charles H. D. Williamson
1,
Shannon L. Johnson
4,
Hajnalka E. Daligault
4 and
Kristin M. Schill
5,*,‡
1
Pathogen and Microbiome Institute, Northern Arizona University, Flagstaff, AZ 86011, USA
2
Institute for Food Safety and Health, Illinois Institute of Technology, Bedford Park, IL 60501, USA
3
Food Science and Nutrition, Illinois Institute of Technology, Chicago, IL 60616, USA
4
Los Alamos National Laboratory, Los Alamos, NM 87545, USA
5
Center for Food Safety and Applied Nutrition, Food and Drug Administration, Bedford Park, IL 60501, USA
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
‡
Current address: Food Research Institute, University of Wisconsin-Madison, Madison, WI 53706, USA.
Toxins 2021, 13(7), 473; https://doi.org/10.3390/toxins13070473
Submission received: 21 May 2021 / Revised: 30 June 2021 / Accepted: 5 July 2021 / Published: 8 July 2021
(This article belongs to the Special Issue New Challenges in Foodborne Botulism Outbreaks)
Browse Figure
Versions Notes

Abstract

:
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.
Keywords:
integration; plasmids; botulinum; clostridium; parabotulinum; sporogenes; argentinense
Key Contribution: Various clostridial species have acquired the ability to produce botulinum neurotoxins after acquisition of bont genes that are carried within large conjugative plasmids, often with subsequent integration of bont gene clusters into the chromosomes of these bacteria. This study describes the first report of rare integration events involving much larger complete plasmid sequences ranging from 140 to 250 kb.

1. Introduction

Botulinum neurotoxins (BoNTs) are extremely potent bacterial proteins produced in a variety of clostridial species. While the mechanism of activity of the toxins and the clinical presentation of botulism is consistent, the proteins themselves differ by as much as 70% in amino acid sequence. Currently, seven serotypes and over 40 toxin subtypes have been described [1].
Botulinum neurotoxin (bont) genes have been identified within seven distinct clostridial species [2]. The species designations in this article are based on phylogenetic relationships, not taxonomic entities. Accordingly, Clostridium parabotulinum, a previously used species designation for proteolytic clostridia that produce BoNT/A, BoNT/B, or BoNT/F, is used here to describe that particular genospecies. This species is sometimes listed as C. botulinum group I. The six additional genospecies are Clostridium botulinum (or C. botulinum group II), Clostridium novyi sensu lato (or C. botulinum group III), Clostridium argentinense, Clostridium baratii, Clostridium butyricum, and Clostridium sporogenes.
The bont genes are positioned in two types of toxin gene clusters ranging in size from approximately 12 to 17 kb that contain both structural genes and, typically, a transcriptional regulator. All gene clusters contain a gene that encodes a non-toxin non-hemagglutinin protein, which is believed to play a role in protecting the toxin molecule as it traverses the digestive system of the host. The two types of gene clusters differ in that one type contains three hemagglutinin genes and the other contains three orfX genes of unknown function [3].
The identification of prophage plasmids encoding bont/C and bont/D in C. botulinum group III [4,5] (now known as C. novyi sensu lato) initiated numerous research studies aimed at elucidating the genomic location of bont genes in other BoNT-producing clostridia. As a result, bont/G genes were reported to reside within plasmids [6]. However, bont/A, bont/B, bont/E, and bont/F genes were historically thought to be chromosomally located [7] due to lack of evidence of extrachromosomal elements prior to the development of genetic tools that are now used to identify them. Over the past two decades, plasmid-borne bont gene clusters have been discovered among BoNT/A and BoNT/B strains [8,9,10] BoNT/E strains [11,12], and BoNT/F strains [13,14,15]. For every toxin type, there are now examples of bont genes within mobile genetic elements and, with the exception of bont/C and bont/D, within the chromosome. Moreover, there are examples of strains producing toxin subtypes BoNT/A2, BoNT/A3, BoNT/B1, BoNT/B2, BoNT/E1, BoNT/E3, and BoNT/E10, where identical gene clusters reside within the chromosome and within plasmids.
These extrachromosomal elements tend to be species specific [16]. For example, the plasmids harbored by C. parabotulinum show no significant similarity with those of either C. botulinum [17] or C. argentinense [16], or with the bacteriophage sequences from C. novyi sensu lato that contain bont/C or bont/D genes [18]. However, plasmids found in numerous C. parabotulinum strains are homologous to plasmids found in the closely related C. sporogenes species.
These bont gene-encoding extrachromosomal elements provide a mechanism for the transfer of bont gene clusters from toxigenic strains to non-toxigenic bacteria, and it is hypothesized that integration of bont genes into the chromosome ensures retention of this feature, as the natural loss of plasmid-borne bont genes has been documented [19] and their artificial loss following repeated subculture or genetic disruption in the laboratory has also been demonstrated [20,21]. In contrast, there has been no demonstration of the natural loss of chromosomally located bont genes in the literature.
The bont gene clusters are integrated into the chromosome via homologous recombination between specific matching DNA sequences within mobile genetic elements and the chromosome. These homologous regions align, the DNA is split between specific nucleotides, and the extrachromosomal DNA sequence that contains the bont gene cluster is inserted at the excision site. Sometimes these insertions occur within genes, such as the insertion of bont/E genes via rarA gene excision [13], the insertion of bont/F6 gene clusters via topB gene excision [22], and the insertion of bont/F3 and bont/F4 genes via pulE gene excision [14]. In several of these cases, a duplicate intact copy of the split gene is also present. However, with other bont gene cluster integrations, the site lies within an intergenic sequence, as illustrated by the integration of the bont/A2 and bont/F1 gene clusters [13]. The split genes and/or neighboring genes often encode proteins that are involved in DNA manipulations, including insertion sequence elements, transposases, and resolvases.
The generation of complete, finished genome sequences allows for mapping of bont gene cluster locations. Comparisons of identical bont gene clusters within plasmids versus those that have integrated onto the chromosome affords us the opportunity to better understand the extent and composition of transferred material and identify surrounding genes or genetic elements that may be involved in the specificity of such transfers. With a few exceptions, the transferred material is limited to DNA sequences of ~17–33 kb containing the bont gene cluster and a few surrounding genes [3].
This study describes the rare integration of entire plasmid sequences of ~150–250 kb within the chromosomes of bacteria from three clostridial species (C. sporogenes, C. argentinense, and C. parabotulinum) using homologous recombination processes that are similar to those responsible for bont gene cluster integrations.

2. Results

2.1. Plasmid Integration in C. sporogenes CDC 1632

Currently, 24 C. sporogenes strains have been identified in the National Center for Biotechnology Information (NCBI) database that carry bont/B genes [20,23,24,25]. Seven of these strains contain bont/B1 gene clusters, 10 contain bont/B2 gene clusters, and five contain bont/B6 gene clusters. Table S1 provides details on these genomes. It is notable that for each of these genomes, including that of CDC 1632 as determined in this study, the putative location of the bont genes has been determined to be within plasmids. A BLASTn comparison of the chromosomal nucleotide sequence of CDC 1632 with that of CDC 67071 plasmid, pNPD7 showed 99.5% identity with 88% coverage, indicating that the bont/B1 gene cluster of CDC 1632 was within a plasmid. Alternatively, using pulsed-field gel electrophoresis (PFGE) followed by Southern blot hybridization with a DNA probe specific to bont/B [10], it was determined that the bont/B gene cluster of CDC 1632 was located within the chromosome. The seemingly contradictory evidence that the bont/B1 gene cluster was within a plasmid and also within the chromosome required further investigation. Whole genomic sequencing and assembly of C. sporogenes strain CDC 1632 verified that it contained no extrachromosomal elements and that the entire identified plasmid sequence was embedded within the chromosome. The depth of coverage for the reads at the chromosomal plasmid insertion sites was 44X and 69X, which is similar to the average read coverage of 52.6X for that genome (Table 1), providing evidence that this was not a mis-assembly. In addition, multiple programs were used to sequence, assemble, and verify the assembly of the closed and finished CDC 1632 genome. Details are included in the Materials and Methods section. This was the first discovery of a fully integrated plasmid within the chromosome of a neurotoxigenic clostridia.
This plasmid sequence was not inserted into the chromosome of CDC 1632 at the designated origin of reference plasmid pNPD7 but was instead inserted in the same orientation beginning at pNPD7 bp 142,310 (CDC 1632 chromosome position bp 3,983,188) and ending at pNPD7 bp 142,309 (CDC 1632 chromosome position bp 4,238,662) (Figure 1A). While pNPD7 is 235,650 bp in length, the inserted region within the chromosome of CDC 1632 is larger, at 255,474 bp. The inserted region within CDC 1632 showed significant similarity to the pNPD7 plasmid but various insertions and deletions have occurred between these DNA segments (Figure 1A), suggesting this plasmid integration event was not recent.
When bont gene clusters are integrated into the chromosome, the insertions are frequently found within a split gene sequence, but sometimes within intergenic sequences. The CDC 1632 plasmid sequence insertion was inserted within an intergenic region, with nearby ltrA genes present at the beginning and end of the inserted sequence. The ltrA genes encode self-splicing group II intron reverse transcriptase/maturase proteins (IEP) that are a component of Ll,ltrB introns. These introns insert foreign DNA from mobile genetic elements into the chromosome [27]. The ltrA protein is multifunctional, mediating RNA splicing, target site recognition, target DNA nicking, insertion of spliced RNA and synthesis of complementary DNA at the insertions site [28]. The Ll.ltrB intron is a component of the TargeTron and ClosTron technologies that are able to generate targeted mutants in E. coli and C. botulinum, respectively [28]. Their location adjacent to the plasmid insertion site within the chromosome of strain CDC 1632 indicates their involvement in facilitating the chromosomal integration of the pNPD7-like plasmid sequence. Five ltrA genes were found in the chromosome of CDC 1632 that were within or flanking the plasmid integration site, and two of these ltrA genes were also located within the CDC 67071 pNPD7 plasmid. The ltrA genes can be separated into two variants having 5% difference in nucleotide sequence (Table 2). One of these variants was exclusive to the CDC 1632 chromosome and included genes that flanked the plasmid insertion sequence (NPD5_3760; bp 3,981,276–3,983,102 and NPD5_4078; bp 4,239,250–4,241,076), while the other variant (NPD5_3758; bp 3,977,797–3,979,620) matched the two ltrA genes within the pNPD7 sequence (NPD7_3844 and NPD7_4079) (Figure 1A).
Two additional ltrA genes were located remotely from the plasmid insertion site within the CDC 1632 chromosome and a single ltrA gene was found within the chromosome of CDC 67071. A number of phage genes were positioned adjacent to one of these two ltrA genes in the CDC 1632 chromosome (NPD5_257; bp 264,722–266,545). The PHASTER program was used to examine the CDC 1632 and CDC 67071 chromosomal genomes for the presence of prophage sequences and complete, intact phage sequences matching that of PHAGE_Clostr_phiCD6356_NC_015626 were discovered adjacent to the ltrA gene located at CDC 1632 NPD5_257 and the lone CDC 67071 ltrA gene (NPD7_1244) (Table 3), indicating that these ltrA genes may be facilitators for the integration of both prophage and plasmid sequences.

2.2. Plasmid Integration in C. argentinense Strain CDC 2741 (GM 77/78)

BoNT/G-producing strains are quite rare. Since their discovery in 1969 [32], only two BoNT/G strains have been isolated in Argentina and six in Switzerland [33,34]. While originally considered to be C. botulinum strains, in 1988 it was proposed that BoNT/G-producing strains be considered a new species, C. argentinense, due to physiological and genetic differences with other C. botulinum strains and their similarity with several non-neurotoxigenic clostridia [34].
Evidence for a plasmid location for bont/G genes was initially provided by Eklund et al. [35] in 1988 and was confirmed by later studies [6]. Our analysis of the eight C. argentinense genomes available in the NCBI database revealed that seven of these genomes, including that of CDC 2741, contained bont/G genes, and comparisons with the pRSJ17 plasmid from strain 89G indicated that these bont/G genes are universally found within closely related plasmids (98% coverage with 100% identity) (Table S2).
Similar to the genome of CDC 1632, genomic analysis of CDC 2741 provided evidence that the plasmid in this strain was also integrated into the chromosome. The genome of strain CDC 2741 consists of 20 separate contigs (AYSO01000001-AYSO01000020) [16]. One of these contigs, AYSO01000020, was 1.1 Mb in length and contained the bont/G gene cluster, suggesting a chromosomal position within the genome. Comparison of CDC 2741 contig AYSO01000020 with reference plasmid pRSJ17_1 of strain 89G revealed the presence of a complete integrated ~140 kb plasmid sequence (100% coverage/100% identity). Read coverage of 156X and 128X were noted for the overlapping reads at the plasmid insertion site, which is similar to the average coverage of 172X for CDC 2741 listed in Table 1. In addition, multiple programs were used to sequence, assemble, and verify the assembly of the CDC 2741 genome.
Alignment of the 89G pRSJ17_1 plasmid sequence with CDC 2741 contig AYSO01000020 revealed that pRSJ17_1 was inserted in reverse orientation within this contig beginning at pRSJ17_1 bp 74,664 (AYSO01000020 bp 140,479) and ending at pRSJ17_1 bp 74,663 (AYSO01000020 bp 406) (Figure 1B). The inserted plasmid sequence began and ended with IS982 elements.
Insertion sequences (IS) are small bacterial transposons containing a transposase, the enzyme that catalyzes DNA movements, flanked by short terminal inverted repeats. IS982 insertion sequences are a class of IS elements that were first discovered in Lactococcus spp. and have subsequently been identified in multiple bacterial species [36]. The IS982 element at the beginning of the inserted plasmid sequence was truncated and the other part of this IS982 element was found at the 3′ region of contig AYSO01000003 (bp 6583–5949). In contrast, only one IS982 element was located within plasmid pRSJ17_1. Hence, the insertion of the plasmid sequence in CDC 2741 may have been facilitated through homologous recombination with the participation of these identical IS982 elements. The nucleotide sequence identity shared between plasmid pRSJ17_1 and the plasmid sequence in contig AYSO01000020, and the lack of DNA insertions or deletions, suggests that this was a recent integration event.

2.3. Plasmid Integration in C. parabotulinum DFPST0006 (NCTC 7273, Beans)

Genomic sequences representing more than 184 C. botulinum (genospecies C. parabotulinum) genomes are listed in the NCBI database as Clostridium botulinum. Unlike C. argentinense and C. sporogenes, the majority of bont gene clusters in C. parabotulinum are found within the chromosome. Plasmid-borne bont genes are found in 24 of the C. parabotulinum genomes listed in the NCBI database (Table S3). The majority of these plasmids contain more than one bont gene cluster (bont/A1 or bont/A2 plus bont/B; bont/B5 plus bont/A4; and bont/B5 plus bont/F2). The remaining genomes contain solitary bont/A3, bont/B, or bont/F5 genes.
The genome of BoNT/B1-producing C. parabotulinum strain DFPST0006, also known as NCTC 7273 or the Beans strain, was recently sequenced in 29 contigs. This strain was isolated from a case of foodborne botulism in the U. S. and its toxin was used to produce the WHO International BoNT/B Antitoxin Standard in 1962 in the UK [37]. The DFPST0006 genome is closely related to that of another U.S. isolate, C. parabotulinum subtype BoNT/B1 strain Okra. However, they differ in their bont gene locations. While the bont/B1 gene cluster in strain Okra is within a 149 kb plasmid (pCLD), in strain DFPST0006 this entire plasmid sequence, including the toxin gene cluster, was located within a 573 kb contig (JACBDK01000020) containing both chromosomal and plasmid sequence. The matching region between contig JACBDK01000020 and the pCLD plasmid sequence shows 100% identity with 99% coverage. Read coverage of 156X and 128X was noted for the overlapping reads at the plasmid insertion site, which is similar to the average read coverage of 100X for DFPST0006 listed in Table 1. Mapping of plasmid pCLD to contig JACBDK01000002 revealed that, as with the other two examples of complete plasmid integration, the insertion points were not aligned with the assembly-designated plasmid origin. The plasmid sequence was split at reference plasmid pCLD position bp 28,138–28,143 with a six-nucleotide overlap and was inserted in reverse orientation within the chromosome of strain DFPST0006 from JACBDK01000020 bp 421,015 to bp 570,638 (Figure 1C). The region in the chromosome where the split occurred (the plasmid insertion site) was within an intergenic region having no adjacent identified genes or repeat gene sequences, so that no facilitating genes or elements were located. The 100% identity of the Okra pCLD plasmid with the integrated plasmid in DFPST0006, coupled with a lack of insertions, deletions, or inversions, suggests that this was a recent integration event.

3. Discussion

It is hypothesized that the acquisition of bont gene clusters by non-neurotoxigenic bacteria offers a distinct advantage to the bacteria, as it elevates these organisms from commensal soil inhabitants to pathogens that are capable of fatally intoxicating large numbers of fish, birds, and mammals, leading to better proliferation and spread. The bont gene clusters may be located within extrachromosomal plasmids or prophage sequences, or within the chromosome. While the presence of bont gene clusters within extrachromosomal elements may be transient due to plasmid or prophage DNA loss, bont gene clusters that reside within the chromosome are stably integrated and less subject to elimination. Thus, the integration of bont gene clusters into the bacterial chromosome provides a greater stability of this trait, ensuring its continuity in subsequent generations.
The integration of bont genes within chromosomes typically involves DNA sequences of 13–33 kb in size that contain the bont gene clusters and a few additional genes. Prior to this study, integration of larger sequences, such as 150–250 kb plasmid sequences, was unknown.
The processes that resulted in the insertion of entire plasmid sequences within the chromosome of C. sporogenes strain CDC 1632, C. argentinense strain CDC 2741, and C. parabotulinum strain DFPST0006 (Beans) are similar to those that facilitated the integration of individual bont genes or gene clusters into the chromosome. The relative rarity of integration of large extrachromosomal plasmids versus that of smaller bont gene clusters may simply be related to their size (140–240 kb versus 17–33 kb), but other as-yet-unknown influences may also be factors.
One unanswered question concerning the integration of these large plasmids is that of reversibility. Do these represent examples of mobile extrachromosomal elements that are forever integrated into chromosomes or do they represent a transitional state that is followed by excision from the chromosome and resumption as an intact mobile genetic element? While the answer to this question is not known, there is evidence that indicates the natural state of the plasmids is extrachromosomal. It is recognized that a major mechanism for the dissemination of virulence traits and/or antibiotic resistance in bacteria is through the transfer of these genes that are located within conjugative plasmids [38]. This is specifically true with both C. sporogenes and C. argentinense, where the overwhelming location for bont genes is within mobile plasmids. Of the 24 sequenced C. sporogenes strains that contain bont genes, there is evidence that the bont genes are within plasmids in 22 of these genomes (the remaining two genomes are too fractured to determine bont gene location). One of the genomes contains a plasmid that is integrated into the chromosome. There is evidence that all seven C. argentinense genomes contain their bont genes within plasmids; one of these contains a plasmid that is integrated into the chromosome. In addition, there is evidence of closely related plasmids containing identical bont/B2 genes in both C. parabotulinum and C. sporogenes, demonstrating that cross-species horizontal transfer of these genes likely occurred via conjugation of these extrachromosomal elements.
As cases where these plasmids have integrated into the chromosome are rarely observed, it is not possible to definitively state that their subsequent excision is not possible. However, with two of the chromosomally integrated plasmids in this study, facilitating genes have been identified that are associated with integration. The plasmid insertion in CDC 1632 involved ltrA genes, which are an essential component in the ClosTron system. The ClosTron system allows directed mutational analysis through targeted DNA insertions using specific introns. This system ensures stable integration of the mutated gene in the chromosome, suggesting the possibility that the chromosomal integration of the plasmid in CDC 1632 might be irreversible. With the chromosomally integrated plasmid in CDC 2741, IS982 insertion elements appear to be involved in facilitating integration. The IS982 elements that remain at the insertion sites are not identical, showing a sequence identity difference of over 2%. Mutations in insertion elements can disable their ability to facilitate horizontal transfer events, and the differences in sequence identity among these specific IS982 elements may be evidence that they are no longer able to facilitate the excision of plasmid DNA from the chromosome. While the above information suggests that chromosomal integration of these plasmids is irreversible, there is no definitive proof either for or against this theory.
The novel findings recorded here expand our knowledge about the movement of bont genes among various clostridial strains and species within extrachromosomal elements, and the processes involved in chromosomal integration of these genes that ensure persistence of neurotoxicity in these lineages.

4. Conclusions

Typically, chromosomal integration of bont genes involves acquisition of toxin gene cluster sequences from extrachromosomal plasmids via homologous recombination. However, the chromosomal integration of complete plasmids is rarely observed in bacteria, and it has been previously unknown among the BoNT-producing clostridia. Examination of the genomes of three different species of BoNT-producing clostridia revealed the first examples of chromosomal integration of large bont gene cluster-containing plasmids, illustrating a novel process for the potential generation of stable lineages of BoNT-producing clostridia. This process was not confined to a single bont serotype or bacterial species and utilized a variety of targeting and facilitating genes.

5. Materials and Methods

5.1. Whole-Genome Sequencing

Total genomic DNA (gDNA) extracted from C. sporogenes BoNT/B1-producing strains CDC 1632 and CDC 67071 was sequenced using Illumina [39] and PacBio [40] technologies. Genome assemblies were performed by the Los Alamos National Laboratory (LANL) Genome Science Group. The Illumina data were assembled using Velvet, version 1.2.08 [41], while the PacBio data were assembled using HGAP, version 3 [42]. Consensus sequences from both assemblers were computationally shredded and assembled using parallel Phrap, version SPS-4.24 [43,44]. The resulting assembly was brought to closed and finished status through both manual and computational finishing efforts using Consed [45] and in-house scripts. The assembled genomes were checked for misassemblies using BridgeMapper, a proprietary PacBio software package, and the assembled genome sequence was corrected by mapping Illumina reads (300x) back to the final consensus sequences using Burrows–Wheeler Alignment (BWA) [46], SAMtools [47] and in-house scripts. Annotations were completed at LANL using an automated system utilizing the Ergatis workflow manager [48] and in-house scripts. Detailed sequence information on these genomes is shown in Table 2.
The gDNA of C. parabotulinum BoNT/B1 strain DFPST0006 was extracted and prepared as described previously [49] and was sequenced using the Illumina MiSeq instrument. Read quality was assessed using FastQC v0.11.8 (Babraham Bioinformatics—FastQC A Quality Control tool for High Throughput Sequence Data). Paired-end reads were trimmed using trimmomatic v0.39 [50]. Reads were assembled into contigs using SPAdes v3.13.0 [51], a reference-free assembler, with k mer sizes of 21, 33, 55, 77, 99, and 127. The contigs were filtered to remove those with length <500 bp and k mer coverage <5. After quality trimming, paired-end reads were aligned to contigs using Bowtie2 v2.4.1 [52]. Base coverage was calculated using SAMtools.
Whole-genome sequencing of the following genomes has been previously described: C. argentinense BoNT/G strain CDC 2741 [16], C. argentinense BoNT/G strain 89G [26], and C. parabotulinum BoNT/B1 strain Okra [9]. Further information on these genomes is listed in Table 1.

5.2. Identification of Bont Genes and Their Genomic Locations

The presence and locations of bont gene-containing genetic elements (chromosomal versus plasmid) were identified using a combination of NCBI BLASTn searches against the bont/B2 gene and plasmid pNPD7 from C. sporogenes CDC 67071, the bont/G gene and C. argentinense 89G plasmid pRSJ17_1, or the bont/B1 gene and C. parabotulinum Okra plasmid pCLD [29]. Searches involved a total of 118 published C. sporogenes genomes, eight published C. argentinense genomes, and 185 published C. parabotulinum genomes, plus additional unpublished genomes in our collections.

5.3. Identification of Prophage Sequences

C. sporogenes chromosomes were also examined for the presence of prophages and prophage sequences using PHASTER online database searches [30,31].

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/toxins13070473/s1, Table S1: Characteristics of plasmids containing bont gene clusters within C. sporogenes strains; Table S2: Characteristics of plasmids containing bont gene clusters within C. argentinense strains; Table S3: Characteristics of plasmids containing bont gene clusters within C. parabotulinum strains.

Author Contributions

K.M.S. and T.J.S. conceptualized and designed experiments; K.M.S., S.L.J. and H.E.D. performed experiments; C.H.D.W., B.I. and R.T. provided comparative analyses of genomic data; T.J.S. contributed DNA and source information for strains and prepared the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the U.S. Food and Drug Administration collaborative grant number 5U19FD005322.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Publicly available datasets were analyzed in this study. The draft genome of newly sequenced strain DFPST0006 was deposited in the NCBI WGS database under Bioproject PRJNA643587 with an accession number of JACBDK000000000. Table 1 contains NCBI accession information for all genomes examined in this study.

Conflicts of Interest

The authors declare they have no conflict of interest.

References

  1. Peck, M.W.; Smith, T.J.; Anniballi, F.; Austin, J.W.; Bano, L.; Bradshaw, M.; Cuervo, P.; Cheng, L.W.; Derman, Y.; Dorner, B.G.; et al. Historical Perspectives and Guidelines for Botulinum Neurotoxin Subtype Nomenclature. Toxins 2017, 9, 38. [Google Scholar] [CrossRef]
  2. Smith, T.J.; Williamson, C.H.; Hill, K.; Sahl, J.W.; Keim, P. Botulinum neurotoxin -producing bacteria—Isn’t it time we called a species a species? MBio 2018, 9. [Google Scholar] [CrossRef] [Green Version]
  3. Smith, T.J.; Williamson, C.H.; Hill, K.K.; Johnson, S.J.; Xie, G.; Anniballi, F.; Auricchio, B.; Fernandez, R.A.; Caballero, P.A.; Keim, P.; et al. 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. Front. Microbiol. 2021, 12, 1–19. [Google Scholar] [CrossRef]
  4. Eklund, M.W.; Poysky, F.T. Interconversion of type C and D strains of Clostridium botulinum by specific bacteriophages. Appl. Microbiol. 1974, 27, 251–258. [Google Scholar] [CrossRef] [PubMed]
  5. Oguma, K.; Iida, H.; Shiozaki, M.; Inoue, K. Antigenicity of converting phages obtained from Clostridium botulinum types C and D. Infect. Immun. 1976, 13, 855–860. [Google Scholar] [CrossRef] [Green Version]
  6. Zhou, Y.; Sugiyama, H.; Nakano, H.; Johnson, E.A. The genes for the Clostridium botulinum type G toxin complex are on a plasmid. Infect. Immun. 1995, 63, 2087–2091. [Google Scholar] [CrossRef] [Green Version]
  7. Popoff, M.R. Ecology of neurotoxigenic strains of clostridia. In Clostridial Neurotoxins, The Molecular Pathogenesis of Tetanus and Botulism; Montecucco, C., Ed.; Springer: Berlin, Germany, 1995. [Google Scholar]
  8. Marshall, K.M.; Bradshaw, M.; Pellett, S.; Johnson, E.A. Plasmid encoded neurotoxin genes in Clostridium botulinum serotype A subtypes. Biochem. Biophys. Res. Commun. 2007, 361, 49–54. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Smith, T.J.; Hill, K.K.; Foley, B.T.; Detter, J.C.; Munk, A.C.; Bruce, D.C.; Doggett, N.A.; Smith, L.A.; Marks, J.D.; Xie, G.; et al. 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Franciosa, G.; Maugliani, A.; Scalfaro, C.; Aureli, P. Evidence that plasmid-borne botulinum neurotoxin type B genes are widespread among Clostridium botulinum serotype B strains. PLoS ONE 2009, 4, e4829. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Zhang, Z.; Hintsa, H.; Chen, Y.; Korkeala, H.; Lindstrom, M. Plasmid-borne type E neurotoxin gene clusters in Clostridium botulinum strains. Appl. Environ. Microbiol. 2013, 79, 3856–3859. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Carter, A.T.; Austin, J.W.; Weedmark, K.A.; Peck, M.W. Evolution of Chromosomal Clostridium botulinum Type E Neurotoxin Gene Clusters: Evidence Provided by Their Rare Plasmid-Borne Counterparts. Genome Biol. Evol. 2016, 8, 540–555. [Google Scholar] [CrossRef]
  13. Hill, K.K.; Xie, G.; Foley, B.T.; Smith, T.J.; Munk, A.C.; Bruce, D.; Smith, L.A.; Brettin, T.S.; Detter, J.C. 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. [Google Scholar] [CrossRef] [PubMed]
  14. Dover, N.; Barash, J.R.; Hill, K.K.; Davenport, K.W.; Teshima, H.; Xie, G.; Arnon, S.S. Clostridium botulinum strain Af84 contains three neurotoxin gene clusters: Bont/A2, bont/F4 and bont/F5. PLoS ONE 2013, 8, e61205. [Google Scholar] [CrossRef]
  15. Raphael, B.H.; Bradshaw, M.; Kalb, S.R.; Joseph, L.A.; Luquez, C.; Barr, J.R.; Johnson, E.A.; Maslanka, S.E. Clostridium botulinum strains producing BoNT/F4 or BoNT/F5. Appl. Environ. Microbiol. 2014, 80, 3250–3257. [Google Scholar] [CrossRef] [Green Version]
  16. Smith, T.J.; Hill, K.K.; Xie, G.; Foley, B.T.; Williamson, C.H.D.; Foster, J.T.; Johnson, S.L.; Chertkov, O.; Teshima, H.; Gibbons, H.S.; et al. Genomic sequences of six botulinum neurotoxin-producing strains representing three clostridial species illustrate the mobility and diversity of botulinum neurotoxin genes. Infect. Genet. Evol. 2015, 30, 102–113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Carter, A.T.; Austin, J.W.; Weedmark, K.A.; Corbett, C.; Peck, M.W. Three classes of plasmid (47-63 kb) carry the type B neurotoxin gene cluster of group II Clostridium botulinum. Genome Biol. Evol. 2014, 6, 2076–2087. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Skarin, H.; Hafstrom, T.; Westerberg, J.; Segerman, B. Clostridium botulinum group III: A group with dual identity shaped by plasmids, phages and mobile elements. BMC Genom. 2011, 12, 185. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Yamakawa, K.; Karasawa, T.; Kakinuma, H.; Maruyama, H.; Takahashi, H.; Nakamura, S. Emergence of Clostridium botulinum type B-like nontoxigenic organisms in a patient with type B infant botulism. J. Clin. Microbiol. 1997, 35, 2163–2164. [Google Scholar] [CrossRef] [Green Version]
  20. Weigand, M.R.; Pena-Gonzalez, A.; Shirey, T.B.; Broeker, R.G.; Ishaq, M.K.; Konstantinidis, K.T.; Raphael, B.H. Implications of genome-based discrimination between Clostridium botulinum Group I and Clostridium sporogenes strains for bacterial taxonomy. Appl. Environ. Microbiol. 2015, 81, 5420–5429. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Pellett, S.; Tepp, W.H.; Stanker, L.H.; Band, P.A.; Johnson, E.A.; Ichtchenko, K. Neuronal targeting, internalization, and biological activity of a recombinant atoxic derivative of botulinum neurotoxin A. Biochem. Biophys. Res. Commun. 2011, 405, 673–677. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Carter, A.T.; Stringer, S.C.; Webb, M.D.; Peck, M.W. The type F6 neurotoxin gene cluster locus of group II Clostridium botulinum has evolved by successive disruption of two different ancestral precursors. Genome Biol. Evol. 2013, 5, 1032–1037. [Google Scholar] [CrossRef] [Green Version]
  23. Williamson, C.H.; Sahl, J.W.; Smith, T.J.; Xie, G.; Foley, B.T.; Smith, L.A.; Fernandez, R.A.; Lindstrom, M.; Korkeala, H.; Keim, P.; et al. Comparative genomic analyses reveal broad diversity in botulinum-toxin-producing Clostridia. BMC Genom. 2016, 17, 180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Smith, T.J.; Xie, G.; Williamson, C.H.D.; Hill, K.K.; Sahl, J.W.; Keim, P.; Johnson, S.L. Genomic characterization of newly completed genomes of botulinum neurotoxin-producing species from Argentina, Australia and Africa. Genome Biol. Evol. 2020. [Google Scholar] [CrossRef] [PubMed]
  25. Brunt, J.; van Vliet, A.H.M.; Carter, A.T.; Stringer, S.C.; Amar, C.; Grant, K.A.; Godbole, G.; Peck, M.W. Diversity of the Genomes and Neurotoxins of Strains of Clostridium botulinum Group I and Clostridium sporogenes Associated with Foodborne, Infant and Wound Botulism. Toxins 2020, 12, 586. [Google Scholar] [CrossRef] [PubMed]
  26. Halpin, J.L.; Hill, K.; Johnson, S.L.; Bruce, D.C.; Shirey, T.B.; Dykes, J.K.; Luquez, C. Finished Whole-Genome Sequence of Clostridium argentinense Producing Botulinum Neurotoxin Type G. Genome Announc. 2017, 5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Lambowitz, A.M.; Zimmerly, S. Group II introns: Mobile ribozymes that invade DNA. Cold Spring Harb. Perspect. Biol. 2011, 3, a003616. [Google Scholar] [CrossRef]
  28. Kuehne, S.A.; Minton, N.P. ClosTron-mediated engineering of Clostridium. Bioengineered 2012, 3, 247–254. [Google Scholar] [CrossRef] [Green Version]
  29. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  30. Arndt, D.; Grant, J.R.; Marcu, A.; Sajed, T.; Pon, A.; Liang, Y.; Wishart, D.S. PHASTER: A better, faster version of the PHAST phage search tool. Nucleic Acids Res. 2016, 44, W16–W21. [Google Scholar] [CrossRef] [Green Version]
  31. Zhou, Y.; Liang, Y.; Lynch, K.H.; Dennis, J.J.; Wishart, D.S. PHAST: A fast phage search tool. Nucleic Acids Res. 2011, 39, W347–W352. [Google Scholar] [CrossRef]
  32. Gimenez, D.F.; Ciccarelli, A.S. Another type of Clostridium botulinum. Zentralbl. Bakteriol. Parasitenk. Infektionskr. Hyg. Abt. I 1970, 215, 221–224. [Google Scholar]
  33. Sonnabend, O.; Sonnabend, W.; Heinzle, R.; Sigirst, T.; Dirnhofer, R.; Krech, U. Isolation of Clostridium botulinum type G and identification of type G botulinal toxin in humans: Report of five sudden unexpected deaths. J. Infect. Dis. 1981, 143, 22–27. [Google Scholar] [CrossRef]
  34. Suen, J.C.; Hatheway, C.L.; Steigerwalt, A.G.; Brenner, D.J. Clostridium argentinense sp. nov.: A genetically homogeneous group composed of all strains of Clostridium botulinum toxin type G and some nontoxigenic strains previously identified as Clostridium subterminale or Clostridium hastiforme. Int. J. Syst. Bacteriol. 1988, 38, 375–381. [Google Scholar] [CrossRef] [Green Version]
  35. Eklund, M.W.; Poysky, F.T.; Mseitif, L.M.; Strom, M.S. Evidence for plasmid-mediated toxin and bacteriocin production in Clostridium botulinum type G. Appl. Environ. Microbiol. 1988, 54, 1405–1408. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Filee, J.; Siguier, P.; Chandler, M. Insertion sequence diversity in archaea. Microbiol. Mol. Biol. Rev. 2007, 71, 121–157. [Google Scholar] [CrossRef] [Green Version]
  37. Bowmer, E.J. Antitoxins of Clostridium botulinum Types A, B, C, D, and E; Preparation and Assay of Proposed International Standards; University of Liverpool: Liverpool, UK, 1962. [Google Scholar]
  38. Guedon, G.; Libante, V.; Coluzzi, C.; Payot, S.; Leblond-Bourget, N. The Obscure World of Integrative and Mobilizable Elements, Highly Widespread Elements that Pirate Bacterial Conjugative Systems. Genes 2017, 8, 337. [Google Scholar] [CrossRef] [Green Version]
  39. Bennett, S. Solexa Ltd. Pharmacogenomics 2004, 5, 433–438. [Google Scholar] [CrossRef]
  40. Eid, J.; Fehr, A.; Gray, J.; Luong, K.; Lyle, J.; Otto, G.; Peluso, P.; Rank, D.; Baybayan, P.; Bettman, B.; et al. Real-time DNA sequencing from single polymerase molecules. Science 2009, 323, 133–138. [Google Scholar] [CrossRef]
  41. Zerbino, D.R.; Birney, E. Velvet: Algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18, 821–829. [Google Scholar] [CrossRef] [Green Version]
  42. Chin, C.S.; Alexander, D.H.; Marks, P.; Klammer, A.A.; Drake, J.; Heiner, C.; Clum, A.; Copeland, A.; Huddleston, J.; Eichler, E.E.; et al. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat. Methods 2013, 10, 563–569. [Google Scholar] [CrossRef]
  43. Ewing, B.; Green, P. Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8, 186–194. [Google Scholar] [CrossRef] [Green Version]
  44. Ewing, B.; Hillier, L.; Wendl, M.C.; Green, P. Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8, 175–185. [Google Scholar] [CrossRef] [Green Version]
  45. Gordon, D.; Green, P. Consed: A graphical editor for next-generation sequencing. Bioinformatics 2013, 29, 2936–2937. [Google Scholar] [CrossRef] [Green Version]
  46. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [Green Version]
  47. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; Genome Project Data Processing Subgroup. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [Green Version]
  48. Hemmerich, C.; Buechlein, A.; Podicheti, R.; Revanna, K.V.; Dong, Q. An Ergatis-based prokaryotic genome annotation web server. Bioinformatics 2010, 26, 1122–1124. [Google Scholar] [CrossRef] [PubMed]
  49. Schill, K.M.; Wang, Y.; Butler, R.R.; Pombert, J.-F.; Reddy, N.R.; Skinner, G.E.; Larkin, J.W. Genetic Diversity of Clostridium sporogenes PA 3679 Isolates Obtained from Different Sources as Resolved by Pulsed-Field Gel Electrophoresis and High-Throughput Sequencing. Appl. Environ. Microbiol. 2016, 82, 384–393. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  51. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [Green Version]
Toxins 13 00473 g001 550
Figure 1. Mapping of nucleotide sequences between bont gene-containing plasmids with matching chromosomal sequence within the related strains. Matching regions are shown with pink boxes and bont gene cluster positions are designated by gray boxes. The numbers indicate the base pair locations of matching sequences. (A) the bont/B1-containing plasmid, pNPD7, of C. sporogenes BoNT/B1 strain CDC 67071 was found within the chromosome of C. sporogenes BoNT/B1 strain CDC 1632, inserted in the same orientation as the chromosome. Two ltrA genes positioned prior to and within the inserted material in the CDC 1632 chromosome are identical to the two ltrA genes located within the pNPD7 plasmid sequence (black arrows). Two additional ltrA genes that flank the CDC 1632 chromosomal sequence and one internal ltrA gene are identical in sequence but differ by 5% from the other ltrA set (red lines). (B) Mapping of bont/G-containing plasmid, pRSJ17_1 of C. argentinense strain 89G to the chromosome of C. argentinense strain CDC 2741 contig AYSO01000020. The plasmid was inserted in the opposite orientation as the chromosome, splitting the plasmid at the IS982 site. An additional intact IS982 element was introduced at the beginning of the inserted sequence beginning at bp 406–1 of contig AYSO01000020 and ending in contig AYSO01000003 (bp 6583–5949). (C) Mapping of the bont/B1-containing plasmid of C. parabotulinum BoNT/B1 strain Okra (pCLD) to the chromosome of C. parabotulinum BoNT/B1 strain DFPST0006 (Beans) contig JACBDK01000020. The plasmid was inserted in the opposite orientation as the chromosome.
Figure 1. Mapping of nucleotide sequences between bont gene-containing plasmids with matching chromosomal sequence within the related strains. Matching regions are shown with pink boxes and bont gene cluster positions are designated by gray boxes. The numbers indicate the base pair locations of matching sequences. (A) the bont/B1-containing plasmid, pNPD7, of C. sporogenes BoNT/B1 strain CDC 67071 was found within the chromosome of C. sporogenes BoNT/B1 strain CDC 1632, inserted in the same orientation as the chromosome. Two ltrA genes positioned prior to and within the inserted material in the CDC 1632 chromosome are identical to the two ltrA genes located within the pNPD7 plasmid sequence (black arrows). Two additional ltrA genes that flank the CDC 1632 chromosomal sequence and one internal ltrA gene are identical in sequence but differ by 5% from the other ltrA set (red lines). (B) Mapping of bont/G-containing plasmid, pRSJ17_1 of C. argentinense strain 89G to the chromosome of C. argentinense strain CDC 2741 contig AYSO01000020. The plasmid was inserted in the opposite orientation as the chromosome, splitting the plasmid at the IS982 site. An additional intact IS982 element was introduced at the beginning of the inserted sequence beginning at bp 406–1 of contig AYSO01000020 and ending in contig AYSO01000003 (bp 6583–5949). (C) Mapping of the bont/B1-containing plasmid of C. parabotulinum BoNT/B1 strain Okra (pCLD) to the chromosome of C. parabotulinum BoNT/B1 strain DFPST0006 (Beans) contig JACBDK01000020. The plasmid was inserted in the opposite orientation as the chromosome.
Toxins 13 00473 g001
Table
Table 1. Genome sequence statistics for the strains included in this study.
Table 1. Genome sequence statistics for the strains included in this study.
Average Read Coverage
StrainGenomic ElementNCBI AccessionSize (Mb)ContigsIlluminaPacBio454
CDC 1632chromosomeCP0132434.393051319X52.6X
CDC 67071chromosomeCP0132424.116551649.91X4.34X
plasmid pNPD7CP0132410.235651664.75X5.05X
CDC 2741chromosomeAYSO0100000 [16]4.7425620171.8X-9.9x
89GchromosomeCP014176 [26]4.662991104.22X18.98X
plasmid pRSJ17_1CP014175 [26]0.140071128.79X36.39X
DFPST0006chromosomeJACBDK01000004.0678329100.5X
OkrachromosomeCP000939 [13]3.958231
plasmid pCLDCP000940 [13]0.148781
Table
Table 2. ltrA gene locations. The ltrA genes in bold font flank the insertion site of the CDC 67071-like plasmid sequence within the CDC 1632 chromosome. The ltrA genes in normal font are found within the integrated plasmid sequence in CDC 1632 and CDC 67071 plasmid pNPD7. The ltrA genes in red font are associated with the integration of a prophage sequence in both chromosomes. * Sequence identities are shown as percent identity compared to the CDC 67071 chromosomal ltrA sequence NPD7_1244 using BLASTn [29].
Table 2. ltrA gene locations. The ltrA genes in bold font flank the insertion site of the CDC 67071-like plasmid sequence within the CDC 1632 chromosome. The ltrA genes in normal font are found within the integrated plasmid sequence in CDC 1632 and CDC 67071 plasmid pNPD7. The ltrA genes in red font are associated with the integration of a prophage sequence in both chromosomes. * Sequence identities are shown as percent identity compared to the CDC 67071 chromosomal ltrA sequence NPD7_1244 using BLASTn [29].
CDC1632 Chromosome CDC 67071 Chromosome
Locus TagLocation (bp)Identity *Locus TagLocation (bp)Identity *
NPD5_257264,722–266,545100%NPD7_12441,338,963–1,340,786100%
NPD5_35983,820,672–3,822,49895%
NPD5_37583,977,797–3,979,620100%CDC 67071 Plasmid
NPD5_37603,981,276–3,983,10295%Locus tagLocation (bp)Identity *
NPD5_37884,000,956–4,002,78295%NPD7_384422,924–24,747100%
NPD5_39404,120,953–4,122,775100%NPD7_4079225,791–227,614100%
NPD5_40784,239,250–4,241,07695%
Table
Table 3. Identification and base pair location of intact prophage sequences in the chromosomes of strains CDC 1632 and CDC 67071 using PHASTER [30,31]. The prophage sequences in bold font are associated with CDC 1632 ltrA sequence NPD5_257 and CDC 67071 ltrA sequence NPD7_1244.
Table 3. Identification and base pair location of intact prophage sequences in the chromosomes of strains CDC 1632 and CDC 67071 using PHASTER [30,31]. The prophage sequences in bold font are associated with CDC 1632 ltrA sequence NPD5_257 and CDC 67071 ltrA sequence NPD7_1244.
DC 1632 Chromosome
SizeStateLocation (bp)Phage Name (Closest Match)
51.9 Kbintact261,289–313,268PHAGE_Clostr_phiCD6356_NC_015262(9)
50.2 Kbintact3,627,037–3,677,326PHAGE_Clostr_phiCTC2B_NC_030951(8)
CDC 67071 Chromosome
SizeStateLocation (bp)Phage Name (Closest Match)
46.5 Kbintact1,345,771–1,392,323PHAGE_Clostr_phiCD6356_NC_015262(9)
CDC 67071 Plasmid
No intact phage present; no incomplete phage associated with ltrA present.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Smith, T.J.; Tian, R.; Imanian, B.; Williamson, C.H.D.; Johnson, S.L.; Daligault, H.E.; Schill, K.M. Integration of Complete Plasmids Containing Bont Genes into Chromosomes of Clostridium parabotulinum, Clostridium sporogenes, and Clostridium argentinense. Toxins 2021, 13, 473. https://doi.org/10.3390/toxins13070473

AMA Style

Smith TJ, Tian R, Imanian B, Williamson CHD, Johnson SL, Daligault HE, Schill KM. Integration of Complete Plasmids Containing Bont Genes into Chromosomes of Clostridium parabotulinum, Clostridium sporogenes, and Clostridium argentinense. Toxins. 2021; 13(7):473. https://doi.org/10.3390/toxins13070473

Chicago/Turabian Style

Smith, Theresa J., Renmao Tian, Behzad Imanian, Charles H. D. Williamson, Shannon L. Johnson, Hajnalka E. Daligault, and Kristin M. Schill. 2021. "Integration of Complete Plasmids Containing Bont Genes into Chromosomes of Clostridium parabotulinum, Clostridium sporogenes, and Clostridium argentinense" Toxins 13, no. 7: 473. https://doi.org/10.3390/toxins13070473

APA Style

Smith, T. J., Tian, R., Imanian, B., Williamson, C. H. D., Johnson, S. L., Daligault, H. E., & Schill, K. M. (2021). Integration of Complete Plasmids Containing Bont Genes into Chromosomes of Clostridium parabotulinum, Clostridium sporogenes, and Clostridium argentinense. Toxins, 13(7), 473. https://doi.org/10.3390/toxins13070473

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop